Viscoacoustic Wave Equation Research Articles

Modeling a pure qP wave in attenuative transverse isotropic media using a complex-valued fractional viscoacoustic wave equation

Seismic anisotropy, characterized by direction-dependent seismic velocity and attenuation, is a fundamental property of the earth, widely observed in field and laboratory measurements. This anisotropy significantly influences seismic wavefield attributes, such as amplitude, phase, and traveltime. Accurate quantification of anisotropy-related seismic wave propagation is essential for global- and regional-scale studies in seismology, such as imaging and inversion. Because the elastic assumption of the earth requires contending with S waves, which are often weak in our seismic recording, we derive a new general complex-valued spatial fractional Laplacian (SFL) viscoacoustic anisotropic pure quasi-P (qP) wave equation (WE) in attenuative transverse isotropic (A-TI) media under the acoustic assumption and use the fixed-order SFL operator to quantify seismic wave propagation in heterogeneous media. Our WE accurately captures anisotropic effects in seismic velocity and attenuation, even in strongly attenuative and anisotropic environments. It also describes the decoupled energy attenuation and velocity dispersion behaviors of the frequency-independent quality factor Q, beneficial for seismic imaging and inversion. Compared with existing viscoacoustic anisotropic pure qP WEs in A-TI media, our WE simplifies computational complexity and is easy to implement, as it decomposes into a complex-valued scalar operator and two differential operators. The evaluated dispersion curves further confirm the accuracy of our approach. Numerical examples in complex geologic settings demonstrate the new equation’s precision, stability, and efficiency in modeling seismic wave propagation.

A high-efficiency Q-compensated pure-viscoacoustic reverse time migration for tilted transversely isotropic media

The attenuation and anisotropy characteristics of real earth media give rise to amplitude loss and phase dispersion during seismic wave propagation. To address these effects on seismic imaging, viscoacoustic anisotropic wave equations expressed by the fractional Laplacian have been derived. However, the huge computational expense associated with multiple Fast Fourier transforms for solving these wave equations makes them unsuitable for industrial applications, especially in three dimensions. Therefore, we first derived a cost-effective pure-viscoacoustic wave equation expressed by the memory-variable in tilted transversely isotropic (TTI) media, based on the standard linear solid model. The newly derived wave equation featuring decoupled amplitude dissipation and phase dispersion terms, can be easily solved using the finite-difference method (FDM). Computational efficiency analyses demonstrate that wavefields simulated by our newly derived wave equation are more efficient compared to the previous pure-viscoacoustic TTI wave equations. The decoupling characteristics of the phase dispersion and amplitude dissipation of the proposed wave equation are illustrated in numerical tests. Additionally, we extend the newly derived wave equation to implement Q-compensated reverse time migration (RTM) in attenuating TTI media. Synthetic examples and field data test demonstrate that the proposed Q-compensated TTI RTM effectively migrate the effects of anisotropy and attenuation, providing high-quality imaging results.

A K-Space-Based Temporal Compensating Scheme for a First-Order Viscoacoustic Wave Equation with Fractional Laplace Operators

Inherent constant Q attenuation can be described using fractional Laplacian operators. Typically, the fractional Laplacian viscoacoustic or viscoelastic wave equations are addressed utilizing the staggered-grid pseudo-spectral (SGPS) method. However, this approach results in time numerical dispersion errors due to the low-order finite difference approximation. In order to address these time-stepping errors, a k-space-based temporal compensating scheme is established to solve the first-order viscoacoustic wave equation. This scheme offers the advantage of being nearly free from grid dispersion for homogeneous media and enhances simulation stability. Numerical examples indicate that the proposed k-space scheme aligns well with analytical solutions for homogeneous media. Additionally, this method demonstrates excellent applicability for complex models and is more efficient due to its ability to adopt a larger time step compared with conventional staggered-grid pseudo-spectral methods.

Source-independent Q-compensated viscoacoustic least-squares reverse time migration

The strong viscosity of the subsurface introduces amplitude absorption and phase-velocity dispersion. Incorrect compensation of the inherent attenuation (the strength of the seismic attenuation can be quantified by the inverse of the quality factor Q, which is defined as [Formula: see text] times the ratio of the stored energy to the lost energy in a single cycle of deformation) can significantly affect imaging quality. Although Q-least-squares reverse time migration ( Q-LSRTM) allows for the compensation of attenuation effects during the iterations, the traditional [Formula: see text]-norm minimization, which is highly sensitive to the source wavelet, poses a challenge in accurately estimating the source wavelet from the field data. Thus, we develop a source-independent Q-LSRTM, in which a convolutional objective function is introduced to replace the [Formula: see text]-norm constraint to mitigate the source wavelet effect. According to the Born approximation, we first linearize the constant-order decoupled fractional Laplacian viscoacoustic wave equation to derive the demigration operator and then construct the corresponding adjoint equation and gradient based on the convolutional objective function, iteratively estimating the reflectivity images. Our method relaxes the sensitivity to the wavelet compared with the conventional [Formula: see text]-norm scheme due to the convolutional objective function, which has the ability to construct the same new source for simulated and observed data. Numerical tests on a layered model, the Marmousi model, and field data demonstrate that our source-independent Q-LSRTM enables us to obtain high-quality reflectivity images even when using incorrect source wavelets.

Characterization of petrophysical and seismic properties for CO2 storage with sensitivity analysis

Saline aquifers are considered as highly favored reservoirs for CO2 sequestration due to their favorable properties. Understanding the impact of saline aquifer properties on the migration and distribution of CO2 plume is crucial. This study focuses on four key parameters—permeability, porosity, formation pressure, and temperature—to characterize the reservoir and analyse the petrophysical and elastic response of CO2. First, we performed reservoir simulations to simulate CO2 saturation, using multiple sets of these four parameters to examine their significance on CO2 saturation and the plume migration speed. Subsequently, the effect of these parameters on the elastic properties is tested using rock physics theory. We established a relationship of compressional wave velocity (Vp) and quality factor (Qp) with the four key parameters, and conducted a sensitivity analysis to test their sensitivity to Vp and Qp. Finally, we utilized visco-acoustic wave equation simulated time-lapse seismic data based on the computed Vp and Qp models, and analysed the impact of CO2 saturation changes on seismic data. As for the above numerical simulations and analysis, we conducted sensitivity analysis using both homogeneous and heterogeneous models. Consistent results are found between homogeneous and heterogeneous models. The permeability is the most sensitive parameter to the CO2 saturation, while porosity emerges as the primary factor affecting both Qp and Vp. Both Qp and Vp increase with the porosity, which contradicts the observations in gas reservoirs. The seismic simulations highlight significant variations in the seismic response to different parameters. We provided analysis for these observations, which serves as a valuable reference for comprehensive CO2 integrity analysis, time-lapse monitoring, injection planning and site selection.

Efficient modeling of fractional Laplacian viscoacoustic wave equation with fractional finite-difference method

The fractional viscoacoustic/viscoelastic wave equation, which accurately quantifies the frequency-independent anelastic effects, has been the focus of seismic industry in recent years. The pseudo-spectral (PS) method stands as one of the most widely used numerical methods for solving the fractional wave equation. However, the PS method often suffers from low accuracy and efficiency, particularly when modeling wave propagation in heterogeneous media. To address these issues, we propose a novel and efficient fractional finite-difference (FD) method for solving the wave equation with fractional Laplacian operators. This method develops an arbitrary high-order FD operator via the generating function of our fractional FD (F-FD) scheme, enhancing accuracy with L2-optimal FD coefficients. Similar to classic FD methods, our F-FD method is characterized by straightforward programming and excellent 3D extensibility. It surpasses the PS method by eliminating the need for Fast Fourier Transform (FFT) and inverse-FFT (IFFT) operations at each time step, offering significant benefits for 3D applications. Consequently, the F-FD method proves more adept for wave-equation-based seismic data processes like imaging and inversion. Compared with existing F-FD methods, our approach uniquely approximates the entire fractional Laplacian operator and stands as a local numerical algorithm, with an adjustable F-FD operator order based on model parameters for enhanced practicality. Accuracy analyses confirm that our method matches the precision of the PS method with a correctly ordered F-FD operator. Numerical examples show that the proposed method has good applicability for complex models. Finally, we have carried out reverse time migration on the Marmousi-2 model, and the imaging profiles indicate that the proposed method can be effectively applied to seismic imaging, demonstrating good practicability.

P-wave propagation in attenuative elliptical VTI media

Elliptical anisotropy is convenient to use as the reference medium in perturbation methods designed to study P-wave propagation for transverse isotropy (TI). We make the elliptically anisotropic TI model attenuative and discuss the corresponding P-wave dispersion relation and the wave equation. Our analysis leads to two conditions in terms of the Thomsen-type parameters, which guarantee that the P-wave slowness surface and the dispersion relation satisfy elliptical equations. We also obtain the viscoacoustic wave equation for such elliptically anisotropic media and solve it for point-source radiation using the correspondence principle. For the constant- Q TI model, we use the weighting function method to derive the viscoacoustic wave equation in differential form. Numerical examples validate the proposed elliptical conditions and illustrate the behavior of the P wavefield in attenuative elliptical TI models.

Power-law frequency-dependent Q simulations in viscoacoustic media using decoupled fractional Laplacians

Quantifying the seismic attenuation of wave propagation in the earth’s interior is essential for studying subsurface structures. Previous approaches for attenuation simulations (e.g., the standard linear solid and the fractional derivative model) are mainly based on the frequency-independent quality factor Q assumption. However, seismic attenuation in high-temperature and high-pressure regions usually exhibits power-law frequency-dependent Q characteristics. To simulate this Q effect in attenuative media, we derive a new viscoacoustic wave equation with decoupled fractional Laplacians in the time domain. Unlike the existing methods using relaxation functions to fit the power-law relationship in a specific frequency band, our equation is directly derived from the approximated complex modulus, which explicitly involves the reference quality factor and fractional exponent parameters. Furthermore, this equation contains two fractional Laplacians, which can easily simulate decoupled amplitude dissipation and phase distortion effects, making it amenable to Q-compensated reverse time migration. In the implementation, a Taylor-series expansion and a pseudospectral method are introduced to solve the fractional Laplacians with variable fractional exponents. Numerical experiments demonstrate the effectiveness of our method for power-law frequency-dependent Q simulations. As a forward-modeling engine, our derived viscoacoustic wave equation is a good supplement to the current Q simulation methods and it could be applied in many seismic applications, such as Q-compensated reverse time migration and full-waveform inversion.

Viscoacoustic VTI media reverse time migration method in the complex-domain based on the staining algorithm

Underground media exhibit characteristics of viscosity and anisotropy, with vertical transversely isotropic (VTI) media being a common example of anisotropic media. Consequently, it is crucial to address the viscoacoustic VTI wave equation during forward modeling and migration imaging processes to attain highly precise imaging outcomes. However, conventional subsalt imaging techniques yield low accuracy and suffer from a low signal-to-noise ratio. To overcome this challenge, this study proposes a staining algorithm-based viscoacoustic vertical transversely isotropic reverse time migration (VTI RTM) in the complex-domain method by deriving the complex-domain viscoacoustic VTI equation. This method requires both the traditional real velocity field and the imaginary velocity field as inputs. The imaginary velocity field contains stained geological bodies and significant areas without data. Seismic waves propagate normally in the real wavefield, while in the imaginary wavefield, reflections occur only when the wavefront encounters stained geological bodies. The proposed method is assessed through testing on both a horizontally layered model and a salt model. Through comparative analysis, it has been verified that the method fully incorporates reflected waves with distinct structural features into the imaginary wavefield. Numerical experiments demonstrate that this method can accurately image subsalt structures.

Q-compensated image-domain least-squares reverse time migration through preconditioned point-spread functions

Image-domain least-squares reverse time migration (IDLSRTM) through point-spread functions (PSFs) has been proven to be a feasible approach to improve the spatial resolution and amplitude fidelity of reflection images recovered by reverse time migration (RTM). However, it usually ignores the earth’s Q-effects, which may lead to an unfocused reflection image with an undesired spatial resolution. In this paper, we develop a Q-compensated IDLSRTM approach (denoted as Q-IDLSRTM) through PSFs, in which we use the viscoacoustic wave equation based on the generalized standard linear solid model to simulate inherent subsurface attenuation and the linear inversion to compensate for the amplitude attenuation. The PSFs are estimated by a round of modeling-migration computation and spatial interpolation on the fly. There are two key points in the developed Q-IDLSRTM approach. The first is that we must apply the deblurring filter as a preconditioner to compensate for the attenuation of image amplitude of PSFs and RTM in a viscoacoustic medium, before the iterative solution. The preconditioned PSFs and RTM images can help us to construct a less ill-posed image-domain inverse problem that can produce an improved image quality and a faster convergence rate, compared with the conventional Q-IDLSRTM approach without the deblurring filter. The second key point is that we can impose the L1-norm constraint and total variation regularization on the reflection image to stabilize the solution of the ill-posed inverse problem. Several 2D and 3D experiments verify that the developed approach can achieve better imaging quality in terms of amplitude fidelity and spatial resolution relative to the conventional Q-IDLSRTM and acoustic IDLSRTM approaches.

Multidimensional Q‐compensated reverse time migration using a high‐efficient decoupled viscoacoustic wave equation

AbstractSeismic waves propagating through attenuating media induce amplitude loss and phase dispersion. Neglecting the attenuation effects during seismic processing results in the imaging profiles with weakened energy, mispositioned interfaces and reduced resolution. To obtain high‐quality imaging results, Q‐compensated reverse time migration is developed. The kernel of the Q‐compensated reverse time migration algorithm is a viscoacoustic wave equation with decoupled amplitude loss and phase dispersion terms. However, the majority of current decoupled viscoacoustic wave equations are solved using the computationally expensive pseudo‐spectral method. To enhance computational efficiency, we initiate our approach from the dispersion relation of a single standard linear solid model. Subsequently, we derive a novel decoupled viscoacoustic wave equation by separating the amplitude loss and phase dispersion terms, previously coupled in the memory variable. The newly derived decoupled viscoacoustic wave equation can be efficiently solved using the finite‐difference method. Then, we reverse the sign of the amplitude loss term of the newly derived viscoacoustic wave equation to implement high‐efficient Q‐compensated reverse time migration based on the finite‐difference method. In addition, we design a regularization term to suppress the high‐frequency noise for stabilizing the wavefield extrapolation. Forward modelling tests validate the decoupled amplitude loss and phase dispersion characteristics of the newly derived viscoacoustic wave equation. Numerical examples in both two‐dimensional and three‐dimensional confirm the effectiveness of the Q‐compensated reverse time migration based on the finite‐difference algorithm in mitigating the attenuation effects in subsurface media and providing high‐quality imaging results.

Viscoacoustic generalized screen propagator in constant-Q model

When seismic waves propagate through the geological formation, there is a significant loss of energy and a decrease in imaging resolution, because of the viscoacoustic properties of subsurface medium. This profoundly impacts seismic wavefield propagation, imaging and interpretation. To accurately image the true structure of subsurface medium, the consensus among geophysicists is to no longer treat subsurface medium as ideal homogeneous medium, but rather to incorporate the viscoacoustic properties of subsurface medium. Based on the generalized screen propagator using conventional acoustic wave equation (acoustic GSP), our developed method introduces viscoacoustic compensation strategy, and derives a one-way wave generalized screen propagator based on time-fractional viscoacoustic wave equation (viscoacoustic GSP). In numerical experiments, we conducted tests on two-dimensional multi-layer model and the Marmousi model. When comparing with the acoustic GSP using the acoustic data, we found that the imaging results of the viscoacoustic GSP using the viscoacoustic data showed a significant attenuation compensation effect, and achieved imaging results for both algorithms were essentially consistent. However, the imaging results of acoustic GSP using viscoacoustic data showed significant attenuation effects, especially for deep subsurface imaging. This indicates that we have proposed an effective method to compensate the attenuated seismic wavefield. Our application on a set of real seismic data demonstrated that the imaging performance of our proposed method in local areas surpassed that of the conventional acoustic GSP. This suggests that our proposed method holds practical value and can more accurately image real subsurface structures while enhancing imaging resolution compared with the conventional acoustic GSP. Finally, with respect to computational efficiency, we gathered statistics on running time to compare our proposed method with conventional Q-RTM, and it is evident that our method exhibits higher computational efficiency. In summary, our proposed viscoacoustic GSP method takes into account the true properties of the medium, still achieves migration results comparable to conventional acoustic GSP.

On the numerical resolution of the second-order viscoacoustic and viscoelastic anisotropic wave equations using the recursive convolution method

Accurate and efficient seismic wave modeling is fundamental for the high-resolution seismic full-waveform inversion and correct interpretation of seismic data. Based on the viscoelastic mechanism of the Generalized Standard Linear Solid and the constitutive relationship between stress and strain, we demonstrate a generalized recursive formula using the Taylor series expansion to directly compute the temporal convolution in the constitutive equation, rather than solving auxiliary differential equations of the memory variables in the conventional method. New sets of the second order viscoacoustic and viscoelastic anisotropic wave equations are established with the new recursive temporal convolutions, which completely replaces the traditional auxiliary differential equations of the memory variables. The error analysis shows the proposed method can offer better accuracy than the existing recursive methods, and is very competitive with the widely used memory variable method. The numerical results comparing the different methods show the accuracy and feasibility of the generalized recursive convolution method to simulate the viscoacoustic and viscoelastic anisotropic waves.

Viscoacoustic least squares reverse-time migration using L1-2 norm sparsity constraint

Abstract Least-squares reverse-time migration (LSRTM) has become an advanced technique for complex structures imaging of the subsurface, as it can provide a higher resolution and more balanced amplitude migrated image than conventional reverse-time migration (RTM). However, the intrinsic attenuation of subsurface introduces amplitude attenuation and phase dispersion of seismic wavefield, which leads to the inverted image kinematically and dynamically inexactitude. Moreover, the imperfect geometry, limited bandwidth of seismic data, and inappropriate modeling kernel etc., would inevitably introduce two side-effects in migrated image, resulting in degradation of LSRTM imaging potential. To alleviate above issues, we present a data-domain sparsity constraint viscoacoustic least-squares reverse-time migration algorithm in this paper. In particular, we utilize the decoupled constant Q fractional Laplacians (DFLs) viscoacoustic wave equation as the modeling kernel to describe the attenuation effects of the subsurface, while a model constraint constructed in the misfit function via L1-2 norm is carried out to clear the migrated artefacts and boost the imaging resolution. Thanks to the excellent performance in sparsity, the drawbacks of unconstraint LSRTM can be effectively mitigated by the L1-2 norm-based regularization. In this paper, we adopt the alternating direction of multipliers method (ADMM) to iteratively address the constrained L1-2 minimization problem by implementing a proximal operator, and three synthetic examples are hired to evaluate the effectiveness and practicability of the proposed strategy. Migration results prove that the proposed scheme can effectively compensate the attenuation effects, improve the resolution, and suppress the migration artifacts of inverted images even in the complex imaging situations.

Seismic reverse time migration based on biaxial wavefield decomposition

Seismic images generated by reverse time migration (RTM) are often contaminated by low-wavenumber artifacts. Decomposition of up- and downgoing waves during propagation can effectively reduce these artifacts. However, decomposition along the vertical axis is not sufficient to eliminate the migration artifacts. We introduce biaxial wavefield decomposition to decompose seismic wavefields in the lateral and vertical directions. After decomposition, there are 16 crosscorrelation terms between different components, and there are only four crosscorrelation terms that can be used separately to image flat layers and tilted interfaces. Within these four terms, the remaining artifacts can be easily identified and reduced by other methods, and the resolution can be further improved by adjusting the weights of the four terms. The decomposition method based on the analytic wavefield via the Hilbert transform has been shown to agree with the viscoacoustic wave equation with fractional spatial derivatives. Therefore, this biaxial wavefield decomposition method is suitable for decomposing the wavefields for the attenuation-compensated RTM. The decomposed and recombined migration results can better image a complex subsurface structure with high resolution.

A Multi-Task Learning Framework of Stable Q-Compensated Reverse Time Migration Based on Fractional Viscoacoustic Wave Equation

Q-compensated reverse time migration (Q-RTM) is a crucial technique in seismic imaging. However, stability is a prominent concern due to the exponential increase in high-frequency ambient noise during seismic wavefield propagation. The two primary strategies for mitigating instability in Q-RTM are regularization and low-pass filtering. Q-RTM instability can be addressed through regularization. However, determining the appropriate regularization parameters is often an experimental process, leading to challenges in accurately recovering the wavefield. Another approach to control instability is low-pass filtering. Nevertheless, selecting the cutoff frequency for different Q values is a complex task. In situations with low signal-to-noise ratios (SNRs) in seismic data, using low-pass filtering can make Q-RTM highly unstable. The need for a small cutoff frequency for stability can result in a significant loss of high-frequency signals. In this study, we propose a multi-task learning (MTL) framework that leverages data-driven concepts to address the issue of amplitude attenuation in seismic records, particularly when dealing with instability during the Q-RTM (reverse time migration with Q-attenuation) process. Our innovative framework is executed using a convolutional neural network. This network has the capability to both predict and compensate for the missing high-frequency components caused by Q-effects while simultaneously reconstructing the low-frequency information present in seismograms. This approach helps mitigate overwhelming instability phenomena and enhances the overall generalization capacity of the model. Numerical examples demonstrate that our Q-RTM results closely align with the reference images, indicating the effectiveness of our proposed MTL frequency-extension method. This method effectively compensates for the attenuation of high-frequency signals and mitigates the instability issues associated with the traditional Q-RTM process.

Frequency-dependent Q simulation and viscoacoustic reverse time migration based on the fractional Zener model

Seismic attenuation is a basic physical property of the earth, which significantly affects the characteristics of seismic wavefields. Accurately simulating wave propagation in the earth is essential to image subsurface structures. Some prevailing methods (e.g., the standard linear solid and fractional Laplacian equation) to describe seismic wave propagation in attenuating media are mainly based on the constant- Q model (CQM), which is valid at room temperature and pressure. However, laboratory measurements suggest that the quality factor Q is a function of frequencies in some regions. To simulate the frequency-dependent Q effect, we derive a viscoacoustic wave equation from the stress-strain relationship of the fractional Zener model (FZM) with variable fractional orders. During the implementation, we separate the real and imaginary parts of the modulus and introduce a low-rank decomposition method to solve the FZM equation. Because the amplitude dissipation and phase dispersion are decoupled, we establish a compensated reverse time migration ( Q-RTM) algorithm to mitigate adverse effects caused by seismic attenuation and improve the quality of seismic migration in frequency-dependent attenuating media. A two-layer model and the BP gas chimney model are used to perform Q-RTM tests. A low-pass filter with a Tukey window function is applied to suppress numerical instability during the compensation. Numerical results demonstrate that our FZM Q-RTM approach can produce high-resolution images with corrected reflector positions and amplitudes. Because the CQM equation ignores the frequency dependence of Q, it may lead to overcompensation in Q-RTM.

A stable Q-compensated reverse-time migration method using a modified fractional viscoacoustic wave equation

Seismic wave propagation inside the earth is usually accompanied by amplitude dissipation and velocity dispersion. It is essential to accurately characterize these effects and make the corresponding compensations in seismic reverse-time migration (RTM) and the following interpretation. Due to the approximation of [Formula: see text] in the derivation, the conventional fractional viscoacoustic wave equation in strongly attenuative media is unsatisfactory. This paper develops a modified viscoacoustic wave equation by combining the relationship between angular frequency and complex wavenumber. The numerical simulation demonstrates its advantage of high accuracy in describing constant- Q effects, especially in strongly attenuative media with small Q values. A truncated Taylor-series expansion algorithm and a pseudospectral method are developed to solve this equation during wave propagation. To address the issue of numerical instability in Q-compensated RTM ( Q-RTM), a stabilized technology with regularization is further developed. Unlike the conventional filtering method implemented in the wavenumber domain, our approach only adds an explicit regularization term, which can be directly computed in the space domain. Furthermore, this regularization depends on velocities and quality factors, significantly improving its applicability in complicated structures. Noise-free data tests demonstrate that our Q-RTM method can effectively correct phase distortion, compensate for amplitude attenuation, and produce high-resolution migrated images without any cutoff-frequency artifacts caused by filtering. Noisy data experiments further verify the antinoise performance and robustness of our method.

Q-compensated wavefield depth extrapolation-based migration using a viscoacoustic wave equation

In a viscoacoustic medium, intrinsic attenuation causes seismic wavefields to attenuate in amplitude and become dispersed in phase, leading to distorted structural imaging and inaccurate migrated amplitudes. To address this problem, viscoacoustic reverse time migration corrects for dispersion and amplitude attenuation effects in wavefield propagation in forward and backward directions. To provide a parallel alternative approach, the time fractional derivative viscoacoustic wave equation in the depth domain is solved and a Q-compensated wavefield depth extrapolation scheme is established to compensate for viscoacoustic effects during recursive wavefield depth extrapolation. This approach decouples the amplitude attenuation and phase dispersion effects from the viscoacoustic vertical wavenumber. To suppress high wavenumber components and address the problem of exponential amplitude growth during wavefield depth extrapolation, we limit the imaginary part of the vertical wavenumber in the frequency wavenumber domain and use an adaptive stabilization solution. Numerical experiments of impulse responses in a viscoacoustic isotropic medium demonstrate that our proposed scheme can observe calculated wavefields exhibiting amplitude attenuation and phase dispersion effects in comparison to wavefields of the acoustic medium, indicating good agreement with theoretical expectations. We also demonstrate the capability of our proposed scheme to recover imaging amplitudes through a variety of numerical experiments, such as imaging of three-layer, Marmousi, and BP gas reservoir models. The results indicate that our proposed scheme can recover amplitude attenuation and phase dispersion effects more accurately than acoustic migration. We also use marine seismic data featuring natural gas hydrates to indicate that our proposed Q-compensated scheme can generate an enhanced imaging result, especially for the bottom simulating reflectors, compared with the conventional imaging algorithm without Q-compensation.

A simple and high-efficiency viscoacoustic reverse time migration calculated by finite difference

The viscoacoustic wave equation with decoupled amplitude attenuation and phase dispersion terms is convenient for implementing reverse time migration (RTM) in attenuating media. However, the traditional viscoacoustic wave equations are expressed by the fractional Laplacian, which requires computationally expensive fast Fourier transforms to perform numerical solutions. This is not conducive for commercial applications with seismic data in terabyte, particularly in three dimensions. Based on the memory-variable-represented viscoacoustic wave equation, the new dispersion- and dissipation-dominated wavefield extrapolation operators that can be solved using the economical finite-difference method (FDM) are developed. In numerical simulation, the staggered-grid FDM is used to solve the temporal and spatial derivatives. The computational efficiency analysis indicates that, compared to the traditional viscoacoustic wave equation, the new viscoacoustic wave equation has an efficiency improvement of approximately more than two times. The accuracy of simulations of amplitude attenuation and phase dispersion is evaluated in viscoacoustic media with varying Q values. The proposed viscoacoustic wave equation is also extended to Q-compensated RTM. Two synthetic tests and one field data test indicate that the Q-compensated RTM can produce high-resolution images by correcting phase dispersion and compensating amplitude loss.