A hierarchical elastic full-waveform inversion scheme based on wavefield separation and the multistep-length approach

Geophysics ◽  
2016 ◽  
Vol 81 (3) ◽  
pp. R99-R123 ◽  
Author(s):  
Zhiming Ren ◽  
Yang Liu
Keyword(s):  
Wave Velocity ◽  
Waveform Inversion ◽  
P Wave ◽  
Model Parameters ◽  
Full Waveform Inversion ◽  
S Wave ◽  
Wave Velocities ◽  
Full Waveform ◽  
Wavefield Separation ◽  
P Wave Velocity

Elastic full-waveform inversion (FWI) updates model parameters by minimizing the residuals of the P- and S-wavefields, resulting in more local minima and serious nonlinearity. In addition, the coupling of different parameters degrades the inversion results. To address these problems, we have developed a hierarchical elastic FWI scheme based on wavefield separation and a multistep-length gradient approach. First, we have derived the gradients expressed by different wave modes; analyzed the crosstalk between various parameters; and evaluated the sensitivity of separated P-wave, separated S-wave, and P- and S-wave misfit functions. Then, a practical four-stage inversion workflow was developed. In the first stage, conventional FWI is used to achieve rough estimates of the P- and S-wave velocities. In the second stage, we only invert the P-wave velocity applying the separated P-wavefields when strong S-wave energy is involved, or we merely update the S-wave velocity by matching the separated S-wavefields for the weak S-wave case. The PP and PS gradient formulas are used in these two cases, respectively. Therefore, the nonlinearity of inversion and the crosstalk between parameters are greatly reduced. In the third stage, the multistep-length gradient scheme is adopted. The density structure can be improved owing to the use of individual step lengths for different parameters. In the fourth stage, we make minor adjustments to the recovered P- and S-wave velocities and density by implementing conventional FWI again. Synthetic examples have determined that our hierarchical FWI scheme with the aforementioned steps obtains more plausible models than the conventional method. Inversion results of each stage and any three stages reveal that wavefield decomposition and the multistep-length approach are helpful to improve the accuracy of velocities and density, respectively, and all the stages of our hierarchical FWI method are necessary to give a good recovery of P- and S-wave velocities and density.

Localizing CO2 at Sleipner — Seismic images versus P-wave velocities from waveform inversion

Geophysics ◽  
2013 ◽  
Vol 78 (3) ◽  
pp. B131-B146 ◽  
Author(s):  
Manuel Queißer ◽  
Satish C. Singh
Keyword(s):  
Wave Velocity ◽  
Waveform Inversion ◽  
P Wave ◽  
Full Waveform Inversion ◽  
Wave Velocities ◽  
Full Waveform ◽  
P Wave Velocity ◽  
Seismic Image ◽  
P Wave Velocities ◽  
Seismic Images

The presence of injected [Formula: see text] in the Utsira Sand at the Sleipner site, Norway, is associated with a high negative P-wave velocity anomaly; that is, a low postinjection velocity and a strong seismic response. Time-lapse seismic imaging of [Formula: see text] injection at Sleipner is thus a viable monitoring tool of the injected [Formula: see text]. The work flow usually involves conventional seismic processing, including stacking, and results in seismic images. Multiple reflections, interference effects such as tuning, and the velocity pushdown effect due to [Formula: see text] injection render these seismic images ambiguous in terms of the localization and the quantification of the [Formula: see text] in the Utsira Sand. Nonetheless, seismic images often form the basis for analyses that aim to quantify the injected [Formula: see text]. We employed elastic 2D full waveform inversion to invert prestack seismic Sleipner data from preinjection (1994) and postinjection (1999) and compared the resulting postinjection P-wave velocity model with the corresponding seismic image. We found that the high-amplitude reflections in the seismic image do not everywhere coincide with low postinjection P-wave velocities. Drawing extensive and integrated conclusions is out of our scope, because this would require full control over the seismic data processing and a more comprehensive forward modeling. For instance, modeling should be done in 3D and an adequate anelasticity formulation should be added. However, the waveform inversion scheme we used accounts for all the aforementioned elastic propagation effects. The results therefore suggested that the exclusive use of seismic images to quantify [Formula: see text] could be revised and full waveform inversion should be added to the analysis toolbox.

scholarly journals Interparameter trade-off quantification for isotropic-elastic full-waveform inversion with various model parameterizations

Geophysics ◽  
2019 ◽  
Vol 84 (2) ◽  
pp. R185-R206 ◽  
Author(s):  
Wenyong Pan ◽  
Kristopher A. Innanen ◽  
Yu Geng ◽  
Junxiao Li
Keyword(s):  
Wave Velocity ◽  
Waveform Inversion ◽  
P Wave ◽  
Physical Parameters ◽  
Full Waveform Inversion ◽  
S Wave ◽  
Trade Off ◽  
Full Waveform ◽  
Trade Offs ◽  
P Wave Velocity

Simultaneous determination of multiple physical parameters using full-waveform inversion (FWI) suffers from interparameter trade-off difficulties. Analyzing the interparameter trade-offs in different model parameterizations of isotropic-elastic FWI, and thus determining the appropriate model parameterization, are critical for efficient inversion and obtaining reliable inverted models. Five different model parameterizations are considered and compared including velocity-density, modulus-density, impedance-density, and two velocity-impedance parameterizations. The scattering radiation patterns are first used for interparameter trade-off analysis. Furthermore, a new framework is developed to evaluate the interparameter trade-off based upon multiparameter Hessian-vector products: Multiparameter point spread functions (MPSFs) and interparameter contamination sensitivity kernels (ICSKs), which provide quantitative, second-order measurements of the interparameter contaminations. In the numerical experiments, the interparameter trade-offs in various model parameterizations are evaluated using the MPSFs and ICSKs. Inversion experiments are carried out with simple Gaussian-anomaly models and a complex Marmousi model. Overall, the parameterization of the P-wave velocity, S-wave velocity, and density, and the parameterization of the P-wave velocity, S-wave velocity, and S-wave impedance perform best for reconstructing all of the physical parameters. Isotropic-elastic FWI of the Hussar low-frequency data set with various model parameterizations verifies our conclusions.

Multi-parameter model building for the Qiuyue structure using four-component ocean-bottom seismometer data

Geophysics ◽  
2021 ◽  
pp. 1-52
Author(s):  
Yuzhu Liu ◽  
Xinquan Huang ◽  
Jizhong Yang ◽  
Xueyi Liu ◽  
Bin Li ◽  
...  
Keyword(s):  
Wave Velocity ◽  
Waveform Inversion ◽  
Velocity Model ◽  
P Wave ◽  
Velocity Analysis ◽  
Full Waveform Inversion ◽  
Ocean Bottom ◽  
Ocean Bottom Seismometer ◽  
Full Waveform ◽  
P Wave Velocity

Thin sand-mud-coal interbedded layers and multiples caused by shallow water pose great challenges to conventional 3D multi-channel seismic techniques used to detect the deeply buried reservoirs in the Qiuyue field. In 2017, a dense ocean-bottom seismometer (OBS) acquisition program acquired a four-component dataset in East China Sea. To delineate the deep reservoir structures in the Qiuyue field, we applied a full-waveform inversion (FWI) workflow to this dense four-component OBS dataset. After preprocessing, including receiver geometry correction, moveout correction, component rotation, and energy transformation from 3D to 2D, a preconditioned first-arrival traveltime tomography based on an improved scattering integral algorithm is applied to construct an initial P-wave velocity model. To eliminate the influence of the wavelet estimation process, a convolutional-wavefield-based objective function for the preprocessed hydrophone component is used during acoustic FWI. By inverting the waveforms associated with early arrivals, a relatively high-resolution underground P-wave velocity model is obtained, with updates at 2.0 km and 4.7 km depth. Initial S-wave velocity and density models are then constructed based on their prior relationships to the P-wave velocity, accompanied by a reciprocal source-independent elastic full-waveform inversion to refine both velocity models. Compared to a traditional workflow, guided by stacking velocity analysis or migration velocity analysis, and using only the pressure component or other single-component, the workflow presented in this study represents a good approach for inverting the four-component OBS dataset to characterize sub-seafloor velocity structures.

scholarly journals Time–frequency windowing in multiparameter elastic FWI of shallow seismic wavefield

2020 ◽  
Vol 222 (2) ◽  
pp. 1164-1177
Author(s):  
Nikolaos Athanasopoulos ◽  
Edgar Manukyan ◽  
Thomas Bohlen ◽  
Hansruedi Maurer
Keyword(s):  
Wave Velocity ◽  
Time Window ◽  
Mass Density ◽  
Waveform Inversion ◽  
P Wave ◽  
S Wave ◽  
Time Frequency ◽  
Full Waveform ◽  
Shallow Seismic ◽  
P Wave Velocity

SUMMARY Full-waveform inversion of shallow seismic wavefields is a promising method to infer multiparameter models of elastic material properties (S-wave velocity, P-wave velocity and mass density) of the shallow subsurface with high resolution. Previous studies used either the refracted Pwaves to reconstructed models of P-wave velocity or the high-amplitude Rayleigh waves to infer the S-wave velocity structure. In this work, we propose a combination of both wavefields using continuous time–frequency windowing. We start with the contribution of refracted P waves and gradually increase the time window to account for scattered body waves, higher mode Rayleigh waves and finally the fundamental Rayleigh wave mode. The opening of the time window is combined with opening the frequency bandwidth of input signals to avoid cycle skipping. Synthetic reconstruction tests revealed that the reconstruction of P-wave velocity model and mass density can be improved. The S-wave velocity reconstruction is still accurate and robust and is slightly benefitted by time–frequency windowing. In a field data application, we observed that time–frequency windowing improves the consistency of multiparameter models. The inferred models are in good agreement with independent geophysical information obtained from ground-penetrating radar and full-waveform inversion of SH waves.

scholarly journals Semiglobal viscoacoustic full-waveform inversion

Geophysics ◽  
2019 ◽  
Vol 84 (2) ◽  
pp. R271-R293 ◽  
Author(s):  
Nuno V. da Silva ◽  
Gang Yao ◽  
Michael Warner
Keyword(s):  
Physical Properties ◽  
Wave Velocity ◽  
Seismic Data ◽  
Waveform Inversion ◽  
P Wave ◽  
Full Waveform Inversion ◽  
Data Set ◽  
Intrinsic Attenuation ◽  
Full Waveform ◽  
P Wave Velocity

Full-waveform inversion deals with estimating physical properties of the earth’s subsurface by matching simulated to recorded seismic data. Intrinsic attenuation in the medium leads to the dispersion of propagating waves and the absorption of energy — media with this type of rheology are not perfectly elastic. Accounting for that effect is necessary to simulate wave propagation in realistic geologic media, leading to the need to estimate intrinsic attenuation from the seismic data. That increases the complexity of the constitutive laws leading to additional issues related to the ill-posed nature of the inverse problem. In particular, the joint estimation of several physical properties increases the null space of the parameter space, leading to a larger domain of ambiguity and increasing the number of different models that can equally well explain the data. We have evaluated a method for the joint inversion of velocity and intrinsic attenuation using semiglobal inversion; this combines quantum particle-swarm optimization for the estimation of the intrinsic attenuation with nested gradient-descent iterations for the estimation of the P-wave velocity. This approach takes advantage of the fact that some physical properties, and in particular the intrinsic attenuation, can be represented using a reduced basis, substantially decreasing the dimension of the search space. We determine the feasibility of the method and its robustness to ambiguity with 2D synthetic examples. The 3D inversion of a field data set for a geologic medium with transversely isotropic anisotropy in velocity indicates the feasibility of the method for inverting large-scale real seismic data and improving the data fitting. The principal benefits of the semiglobal multiparameter inversion are the recovery of the intrinsic attenuation from the data and the recovery of the true undispersed infinite-frequency P-wave velocity, while mitigating ambiguity between the estimated parameters.

scholarly journals 2-D multiparameter viscoelastic shallow-seismic full-waveform inversion: reconstruction tests and first field-data application

2020 ◽  
Vol 222 (1) ◽  
pp. 560-571
Author(s):  
Lingli Gao ◽  
Yudi Pan ◽  
Thomas Bohlen
Keyword(s):  
Wave Velocity ◽  
Wave Attenuation ◽  
Waveform Inversion ◽  
P Wave ◽  
S Wave ◽  
Velocity Models ◽  
Full Waveform ◽  
Shallow Seismic ◽  
P Wave Velocity ◽  
Data Application

SUMMARY 2-D full-waveform inversion (FWI) of shallow-seismic wavefields has recently become a novel way to reconstruct S-wave velocity models of the shallow subsurface with high vertical and lateral resolution. In most applications, seismic wave attenuation is ignored or considered as a passive modelling parameter only. In this study, we explore the feasibility and performance of multiparameter viscoelastic 2-D FWI in which seismic velocities and attenuation of P and S waves, respectively, and mass density are inverted simultaneously. Synthetic reconstruction experiments reveal that multiple crosstalks between all viscoelastic material parameters may occur. The reconstruction of S-wave velocity is always robust and of high quality. The parameters P-wave velocity and density exhibit weaker sensitivity and can be reconstructed more reliably by multiparameter viscoelastic FWI. Anomalies in S-wave attenuation can be recovered but with limited resolution. In a field-data application, a small-scale refilled trench is nicely delineated as a low P- and S-wave velocity anomaly. The reconstruction of P-wave velocity is improved by the simultaneous inversion of attenuation. The reconstructed S-wave attenuation reveals higher attenuation in the shallow weathering zone and weaker attenuation below. The variations in the reconstructed P- and S-wave velocity models are consistent with the reflectivity observed in a ground penetrating radar (GPR) profile.

Elastic full waveform inversion of multicomponent ocean-bottom cable seismic data: Application to Alba Field, U. K. North Sea

Geophysics ◽  
2010 ◽  
Vol 75 (6) ◽  
pp. R109-R119 ◽  
Author(s):  
Timothy J. Sears ◽  
Penny J. Barton ◽  
Satish C. Singh
Keyword(s):  
Wave Velocity ◽  
Seismic Data ◽  
Waveform Inversion ◽  
P Wave ◽  
Full Waveform Inversion ◽  
Ocean Bottom ◽  
S Wave ◽  
Wave Data ◽  
Data Inversion ◽  
Full Waveform

Elastic full waveform inversion of multichannel seismic data represents a data-driven form of analysis leading to direct quantification of the subsurface elastic parameters in the depth domain. Previous studies have focused on marine streamer data using acoustic or elastic inversion schemes for the inversion of P-wave data. In this paper, P- and S-wave velocities are inverted for using wide-angle multicomponent ocean-bottom cable (OBC) seismic data. Inversion is undertaken using a two-dimensional elastic algorithm operating in the time domain, which allows accurate modeling and inversion of the full elastic wavefield, including P- and mode-converted PS-waves and their respective amplitude variation with offset (AVO) responses. Results are presented from the application of this technique to an OBC seismic data set from the Alba Field, North Sea. After building an initial velocity model and extracting a seismic wavelet, the data are inverted instages. In the first stage, the intermediate wavelength P-wave velocity structure is recovered from the wide-angle data and then the short-scale detail from near-offset data using P-wave data on the [Formula: see text] (vertical geophone) component. In the second stage, intermediate wavelengths of S-wave velocity are inverted for, which exploits the information captured in the P-wave’s elastic AVO response. In the third stage, the earlier models are built on to invert mode-converted PS-wave events on the [Formula: see text] (horizontal geophone) component for S-wave velocity, targeting first shallow and then deeper structure. Inversion of [Formula: see text] alone has been able to delineate the Alba Field in P- and S-wave velocity, with the main field and outlier sands visible on the 2D results. Inversion of PS-wave data has demonstrated the potential of using converted waves to resolve shorter wavelength detail. Even at the low frequencies [Formula: see text] inverted here, improved spatial resolution was obtained by inverting S-wave data compared with P-wave data inversion results.

scholarly journals Elastic reflection-based waveform inversion with a nonlinear approach

Geophysics ◽  
2017 ◽  
Vol 82 (6) ◽  
pp. R309-R321 ◽  
Author(s):  
Qiang Guo ◽  
Tariq Alkhalifah
Keyword(s):  
Wave Velocity ◽  
Equivalent Stress ◽  
Waveform Inversion ◽  
P Wave ◽  
S Wave ◽  
Velocity Models ◽  
Wave Velocities ◽  
Full Waveform ◽  
Highly Nonlinear ◽  
Elastic Reflection

Full-waveform inversion (FWI) is a highly nonlinear problem due to the complex reflectivity of the earth, and this nonlinearity only increases under the more expensive elastic assumption. In elastic media, we need a good initial P-wave velocity and even better initial S-wave velocity models with accurate representation of the low model wavenumbers for FWI to converge. However, inverting for the low-wavenumber components of P- and S-wave velocities using reflection waveform inversion (RWI) with an objective to fit the reflection shape, rather than produce reflections, may mitigate the limitations of FWI. Because FWI, performing as a migration operator, is preferred of the high-wavenumber updates along reflectors. We have developed an elastic RWI that inverts for the low-wavenumber and perturbation components of the P- and S-wave velocities. To generate the full elastic reflection wavefields, we derive an equivalent stress source made up by the inverted model perturbations and incident wavefields. We update the perturbation and propagation parts of the velocity models in a nested fashion. Applications on the synthetic isotropic models and field data indicate that our method can efficiently update the low- and high-wavenumber parts of the models.

Sign in / Sign up

Export Citation Format

Share Document