scholarly journals Elastic reflection waveform inversion with variable density

Geophysics ◽  
2019 ◽  
Vol 84 (4) ◽  
pp. R553-R567 ◽  
Author(s):  
Yuanyuan Li ◽  
Qiang Guo ◽  
Zhenchun Li ◽  
Tariq Alkhalifah
Keyword(s):  
Waveform Inversion ◽  
Variable Density ◽  
Background Model ◽  
Reverse Time ◽  
S Wave ◽  
Low Frequencies ◽  
Wave Velocities ◽  
Reflection Data ◽  
Time Migration ◽  
Background Models

Elastic full-waveform inversion (FWI) provides a better description of the subsurface information than those given by the acoustic assumption. However, it suffers from a more serious cycle-skipping problem compared with the latter. Reflection waveform inversion (RWI) is able to build a good background model, which can serve as an initial model for elastic FWI. Because, in RWI, we use the model perturbation to explicitly fit reflections, such perturbations should include density, which mainly affects the dynamics. We applied Born modeling to generate synthetic reflection data using optimized perturbations of the P- and S-wave velocities and density. The inversion for the perturbations of the P- and S-wave velocities and density is similar to elastic least-squares reverse time migration. An incorrect background model will lead to misfits mainly at the far offsets, which can be used to update the background P- and S-wave velocities along the reflection wavepath. We optimize the perturbations and background models in an alternate way. We use two synthetic examples and a field-data case to demonstrate our proposed elastic RWI algorithm. The results indicate that our elastic RWI with variable density is able to build reasonably good background models for elastic FWI with the absence of low frequencies, and it can deal with the variable density, which is required in real cases.

Elastic reflection waveform inversion based on the decomposition of sensitivity kernels

Geophysics ◽  
2019 ◽  
Vol 84 (2) ◽  
pp. R235-R250 ◽  
Author(s):  
Zhiming Ren ◽  
Zhenchun Li ◽  
Bingluo Gu
Keyword(s):  
Wave Velocity ◽  
Waveform Inversion ◽  
Velocity Model ◽  
P Wave ◽  
Background Model ◽  
Reverse Time ◽  
S Wave ◽  
Velocity Models ◽  
Time Migration ◽  
Starting Model

Full-waveform inversion (FWI) has the potential to obtain an accurate velocity model. Nevertheless, it depends strongly on the low-frequency data and the initial model. When the starting model is far from the real model, FWI tends to converge to a local minimum. Based on a scale separation of the model (into the background model and reflectivity model), reflection waveform inversion (RWI) can separate out the tomography term in the conventional FWI kernel and invert for the long-wavelength components of the velocity model by smearing the reflected wave residuals along the transmission (or “rabbit-ear”) paths. We have developed a new elastic RWI method to build the P- and S-wave velocity macromodels. Our method exploits a traveltime-based misfit function to highlight the contribution of tomography terms in the sensitivity kernels and a sensitivity kernel decomposition scheme based on the P- and S-wave separation to suppress the high-wavenumber artifacts caused by the crosstalk of different wave modes. Numerical examples reveal that the gradients of the background models become sufficiently smooth owing to the decomposition of sensitivity kernels and the traveltime-based misfit function. We implement our elastic RWI in an alternating way. At each loop, the reflectivity model is generated by elastic least-squares reverse time migration, and then the background model is updated using the separated traveltime kernels. Our RWI method has been successfully applied in synthetic and real reflection seismic data. Inversion results demonstrate that the proposed method can retrieve preferable low-wavenumber components of the P- and S-wave velocity models, which are reliable to serve as a starting model for conventional elastic FWI. Also, our method with a two-stage inversion workflow, first updating the P-wave velocity using the PP kernels and then updating the S-wave velocity using the PS kernels, is feasible and robust even when P- and S-wave velocities have different structures.

True-amplitude linearized waveform inversion with the quasi-elastic wave equation

Geophysics ◽  
2019 ◽  
Vol 84 (6) ◽  
pp. R827-R844 ◽  
Author(s):  
Zongcai Feng ◽  
Gerard Schuster
Keyword(s):  
Wave Equation ◽  
Elastic Wave ◽  
Wave Equations ◽  
Waveform Inversion ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
S Wave ◽  
Elastic Wave Equation ◽  
P Waves ◽  
Time Migration

We present a quasi-elastic wave equation as a function of the pressure variable, which can accurately model PP reflections with elastic amplitude variation with offset effects under the first-order Born approximation. The kinematic part of the quasi-elastic wave equation accurately models the propagation of P waves, whereas the virtual-source part, which models the amplitudes of reflections, is a function of the perturbations of density and Lamé parameters [Formula: see text] and [Formula: see text]. The quasi-elastic wave equation generates a scattering radiation pattern that is exactly the same as that for the elastic wave equation, and only requires the solution of two acoustic wave equations for each shot gather. This means that the quasi-elastic wave equation can be used for true-amplitude linearized waveform inversion (also known as least-squares reverse time migration) of elastic PP reflections, in which the corresponding misfit gradients are with respect to the perturbations of density and the P- and S-wave impedances. The perturbations of elastic parameters are iteratively updated by minimizing the [Formula: see text]-norm of the difference between the recorded PP reflections and the predicted pressure data modeled from the quasi-elastic wave equation. Numerical tests on synthetic and field data indicate that true-amplitude linearized waveform inversion using the quasi-elastic wave equation can account for the elastic PP amplitudes and provide a robust estimate of the perturbations of P- and S-wave impedances and, in some cases, the density. In addition, true-amplitude linearized waveform inversion provides images with a wider bandwidth and fewer artifacts because the PP amplitudes are accurately explained. We also determine the 2D scalar quasi-elastic wave equation for P-SV reflections and the 3D vector equation for PS reflections.

Elastic model low- to intermediate-wavenumber inversion using reflection traveltime and waveform of multicomponent seismic data

Geophysics ◽  
2019 ◽  
Vol 84 (1) ◽  
pp. R109-R123 ◽  
Author(s):  
Wencai Xu ◽  
Tengfei Wang ◽  
Jiubing Cheng
Keyword(s):  
Waveform Inversion ◽  
Elastic Model ◽  
S Wave ◽  
Reflected Waves ◽  
Wave Velocities ◽  
Reflection Data ◽  
Model Background ◽  
Mode Decomposition ◽  
Model Components ◽  
Two Stages

Low-, intermediate-, and high-wavenumber components of P- and S-wave velocities jointly influence the elastic wave propagation and scattering in an isotropic medium. By taking advantage of all information in the data, elastic full-waveform inversion (E-FWI) has the potential to recover these model components. However, if the transmitted wave data are insufficient to illuminate the deeper part of the subsurface, we should rely on the solutions using reflection data. To reduce the nonlinearity of waveform inversion, we choose to decouple the effects of the model background and perturbation on the reflected waves within a linearized inversion framework. This resorts to three stages aiming to gradually fit the traveltimes and waveforms of the reflected PP and PS waves based on data or gradient preconditioning through P/S mode decomposition. For the first two stages, once the multicomponent seismograms have been separated into PP and PS reflection recordings, reflection traveltime inversion using an acoustic wave propagator (A-RTI) can successively recover the low-wavenumber components of P- and S-wave velocities. In the last stage, starting from the models having reliable low-wavenumber components, elastic reflection waveform inversion (E-RWI) can easily get out of the local minima and continue to retrieve the increasing wavenumber features sensitive to the waveform and amplitude variations. This is supported by gradient preconditioning through P/S mode decomposition of the extrapolated normal and adjoint wavefields, and alternately updating model background and high-wavenumber components in terms of linearized least-squares inversion. Numerical examples have demonstrated the performance of our E-RWI approach and the validity of the three-stage inversion workflow.

Breakthrough acquisition and technologies for subsalt imaging

Geophysics ◽  
2011 ◽  
Vol 76 (5) ◽  
pp. WB41-WB51 ◽  
Author(s):  
Denes Vigh ◽  
Jerry Kapoor ◽  
Nick Moldoveanu ◽  
Hongyan Li
Keyword(s):  
Velocity Field ◽  
Inversion Layer ◽  
Waveform Inversion ◽  
Full Waveform Inversion ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Data Set ◽  
Low Frequencies ◽  
Full Waveform ◽  
Time Migration

The recently introduced method of wide-azimuth data acquisition offers better illumination, noise attenuation, and lower frequencies to more accurately determine a velocity field for imaging. For the field data experiment to demonstrate the technologies, we used a Gulf of Mexico (GOM) wide-azimuth data set that allows us to take advantage of possible low frequencies, relatively large crossline offsets, and increased illumination. The input data was processed through true 3D azimuthal surface-related multiple elimination (SRME) with zero-phasing and debubble. Eliminating the surface-related multiples aids the velocity determination and helps uncover the subsalt sediments at the final imaging stage. After the initial velocity derivation, which was constrained to wells and geology, full-waveform inversion (FWI) was used to further update the velocity field to achieve an enhanced image. The methodology used follows the top-down approach where suprasalt sediment model is developed followed by the top of salt, salt flanks, base of salt, and finished with a limited subsalt update. To approximate the observed data by using an acoustic inversion procedure, the propagator includes effects of attenuation, anisotropy, acquisition source, and receiver depth. The geological environment is salt related, which implies that the observed data is highly elastic, even though it is input to an acoustic full waveform inversion. To use the proper constraints for the inversion, layer-stripping method is used to develop the high-resolution velocity field. The inversion stages were carefully quality controlled through gather displays to ensure the kinematics were honored. We then demonstrated the benefit of the FWI velocity field by comparing the images derived with the traditional ray-based tomographic velocity field versus the velocity field derived by FWI performing reverse time migration to produce these images. Finally, the images were compared at key well locations to evaluate the robustness of the workflow.

scholarly journals Elastic least-squares reverse time migration via linearized elastic full-waveform inversion with pseudo-Hessian preconditioning

Geophysics ◽  
2017 ◽  
Vol 82 (5) ◽  
pp. S341-S358 ◽  
Author(s):  
Ke Chen ◽  
Mauricio D. Sacchi
Keyword(s):  
Least Squares ◽  
Waveform Inversion ◽  
Wave Impedance ◽  
Full Waveform Inversion ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
S Wave ◽  
Full Waveform ◽  
Time Migration ◽  
Hessian Operator

Time-domain elastic least-squares reverse time migration (LSRTM) is formulated as a linearized elastic full-waveform inversion problem. The elastic Born approximation and elastic reverse time migration (RTM) operators are derived from the time-domain continuous adjoint-state method. Our approach defines P- and S-wave impedance perturbations as unknown elastic images. Our algorithm is obtained using continuous functional analysis in which the problem is discretized at the final stage (optimize-before-discretize approach). The discretized numerical versions of the elastic Born operator and its adjoint (elastic RTM operator) pass the dot-product test. The conjugate gradient least-squares method is used to solve the least-squares migration quadratic optimization problem. In other words, the Hessian operator for elastic LSRTM is implicitly inverted via a matrix-free algorithm that only requires the action of forward and adjoint operators on vectors. The diagonal of the pseudo-Hessian operator is used to design a preconditioning operator to accelerate the convergence of the elastic LSRTM. The elastic LSRTM provides higher resolution images with fewer artifacts and a superior balance of amplitudes when compared with elastic RTM. More important, elastic LSRTM can mitigate crosstalk between the P- and S-wave impedance perturbations given that the off-diagonal elements of the Hessian are attenuated via the inversion.

Vector-based elastic reverse time migration based on scalar imaging condition

Geophysics ◽  
2017 ◽  
Vol 82 (2) ◽  
pp. S111-S127 ◽  
Author(s):  
Qizhen Du ◽  
ChengFeng Guo ◽  
Qiang Zhao ◽  
Xufei Gong ◽  
Chengxiang Wang ◽  
...  
Keyword(s):  
Alternative Methods ◽  
Polarity Reversal ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
S Wave ◽  
Converted Wave ◽  
Imaging Condition ◽  
Time Migration ◽  
Wavefield Separation ◽  
Wave Modes

The scalar images (PP, PS, SP, and SS) of elastic reverse time migration (ERTM) can be generated by applying an imaging condition as crosscorrelation of pure wave modes. In conventional ERTM, Helmholtz decomposition is commonly applied in wavefield separation, which leads to a polarity reversal problem in converted-wave images because of the opposite polarity distributions of the S-wavefields. Polarity reversal of the converted-wave image will cause destructive interference when stacking over multiple shots. Besides, in the 3D case, the curl calculation generates a vector S-wave, which makes it impossible to produce scalar PS, SP, and SS images with the crosscorrelation imaging condition. We evaluate a vector-based ERTM (VB-ERTM) method to address these problems. In VB-ERTM, an amplitude-preserved wavefield separation method based on decoupled elastic wave equation is exploited to obtain the pure wave modes. The output separated wavefields are both vectorial. To obtain the scalar images, the scalar imaging condition in which the scalar product of two vector wavefields with source-normalized illumination is exploited to produce scalar images instead of correlating Cartesian components or magnitude of the vector P- and S-wave modes. Compared with alternative methods for correcting the polarity reversal of PS and SP images, our ERTM solution is more stable and simple. Besides these four scalar images, the VB-ERTM method generates another PP-mode image by using the auxiliary stress wavefields. Several 2D and 3D numerical examples are evaluated to demonstrate the potential of our ERTM method.

From acoustic to elastic inverse extended Born modeling: a first insight in the marine environment

Geophysics ◽  
2021 ◽  
pp. 1-73
Author(s):  
Milad Farshad ◽  
Hervé Chauris
Keyword(s):  
Least Squares ◽  
Marine Environment ◽  
Computational Cost ◽  
Physical Parameters ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
S Wave ◽  
Time Migration ◽  
Physical Density ◽  
Acoustic Approximation

Elastic least-squares reverse time migration is the state-of-the-art linear imaging technique to retrieve high-resolution quantitative subsurface images. A successful application requires many migration/modeling cycles. To accelerate the convergence rate, various pseudoinverse Born operators have been proposed, providing quantitative results within a single iteration, while having roughly the same computational cost as reverse time migration. However, these are based on the acoustic approximation, leading to possible inaccurate amplitude predictions as well as the ignorance of S-wave effects. To solve this problem, we extend the pseudoinverse Born operator from acoustic to elastic media to account for the elastic amplitudes of PP reflections and provide an estimate of physical density, P- and S-wave impedance models. We restrict the extension to marine environment, with the recording of pressure waves at the receiver positions. Firstly, we replace the acoustic Green's functions by their elastic version, without modifying the structure of the original pseudoinverse Born operator. We then apply a Radon transform to the results of the first step to calculate the angle-dependent response. Finally, we simultaneously invert for the physical parameters using a weighted least-squares method. Through numerical experiments, we first illustrate the consequences of acoustic approximation on elastic data, leading to inaccurate parameter inversion as well as to artificial reflector inclusion. Then we demonstrate that our method can simultaneously invert for elastic parameters in the presence of complex uncorrelated structures, inaccurate background models, and Gaussian noisy data.

Reducing artifacts of elastic reverse time migration by the deprimary technique

Geophysics ◽  
2018 ◽  
Vol 83 (6) ◽  
pp. S569-S577 ◽  
Author(s):  
Yang Zhao ◽  
Houzhu Zhang ◽  
Jidong Yang ◽  
Tong Fei
Keyword(s):  
Complex Function ◽  
Noise Removal ◽  
Pure Mode ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
S Wave ◽  
Imaging Condition ◽  
Velocity Models ◽  
Time Migration ◽  
Wavefield Decomposition

Using the two-way elastic-wave equation, elastic reverse time migration (ERTM) is superior to acoustic RTM because ERTM can handle mode conversions and S-wave propagations in complex realistic subsurface. However, ERTM results may not only contain classical backscattering noises, but they may also suffer from false images associated with primary P- and S-wave reflections along their nonphysical paths. These false images are produced by specific wave paths in migration velocity models in the presence of sharp interfaces or strong velocity contrasts. We have addressed these issues explicitly by introducing a primary noise removal strategy into ERTM, in which the up- and downgoing waves are efficiently separated from the pure-mode vector P- and S-wavefields during source- and receiver-side wavefield extrapolation. Specifically, we investigate a new method of vector wavefield decomposition, which allows us to produce the same phases and amplitudes for the separated P- and S-wavefields as those of the input elastic wavefields. A complex function involved with the Hilbert transform is used in up- and downgoing wavefield decomposition. Our approach is cost effective and avoids the large storage of wavefield snapshots that is required by the conventional wavefield separation technique. A modified dot-product imaging condition is proposed to produce multicomponent PP-, PS-, SP-, and SS-images. We apply our imaging condition to two synthetic models, and we demonstrate the improvement on the image quality of ERTM.

Imaging conditions for elastic reverse time migration

Geophysics ◽  
2019 ◽  
Vol 84 (2) ◽  
pp. S95-S111 ◽  
Author(s):  
Wei Zhang ◽  
Ying Shi
Keyword(s):  
Decomposition Method ◽  
Inner Product ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
S Wave ◽  
Imaging Condition ◽  
Subsurface Structures ◽  
Time Migration ◽  
Imaging Conditions ◽  
Wave Modes

Elastic reverse time migration (RTM) has the ability to retrieve accurately migrated images of complex subsurface structures by imaging the multicomponent seismic data. However, the imaging condition applied in elastic RTM significantly influences the quality of the migrated images. We evaluated three kinds of imaging conditions in elastic RTM. The first kind of imaging condition involves the crosscorrelation between the Cartesian components of the particle-velocity wavefields to yield migrated images of subsurface structures. An alternative crosscorrelation imaging condition between the separated pure wave modes obtained by a Helmholtz-like decomposition method could produce reflectivity images with explicit physical meaning and fewer crosstalk artifacts. A drawback of this approach, though, was that the polarity reversal of the separated S-wave could cause destructive interference in the converted-wave image after stacking over multiple shots. Unlike the conventional decomposition method, the elastic wavefields can also be decomposed in the vector domain using the decoupled elastic wave equation, which preserves the amplitude and phase information of the original elastic wavefields. We have developed an inner-product imaging condition to match the vector-separated P- and S-wave modes to obtain scalar reflectivity images of the subsurface. Moreover, an auxiliary P-wave stress image can supplement the elastic imaging. Using synthetic examples with a layered model, the Marmousi 2 model, and a fault model, we determined that the inner-product imaging condition has prominent advantages over the other two imaging conditions and generates images with preserved amplitude and phase attributes.

Sign in / Sign up

Export Citation Format

Share Document