A synthetic 2-D-elastic full-waveform inversion study for the Ecuadorian margin—resolution capability of a dense onshore/offshore network

2020 ◽  
Vol 222 (2) ◽  
pp. 1236-1244
Author(s):  
L Gassner ◽  
N Thiel ◽  
A Rietbrock
Keyword(s):  
Subduction Zone ◽  
Waveform Inversion ◽  
Wave Model ◽  
P Wave ◽  
Full Waveform Inversion ◽  
S Wave ◽  
Traveltime Tomography ◽  
Full Waveform ◽  
Local Earthquake ◽  
Source Orientation

SUMMARY Subduction zones are the places on the Earth where the greatest earthquakes occur. It is now widely accepted that seismic asperities at the interface of subducting plates play a major role in whether a region of a subduction zone behaves seismically, creating strong earthquakes or exhibits aseismic slip. In the last decades, huge advances have been made to decipher the underlying processes; however, the physical parameters along the subduction zone interfaces are still not very well known due to a sparsity of high-resolution experiments and significant costs associated with amphibious seismic experiments. Therefore, synthetic tests are needed to investigate the potential of currently possible high density seismic deployments and to aid future experiment design. As standard local earthquake traveltime tomography in a subduction zone setting cannot resolve structures on a kilometre scale at depth, we explore the suitability of full-waveform inversion (FWI) to increase resolution by using amplitude and phase information in the recorded earthquake seismograms. We apply 2-D-elastic FWI to synthetic earthquake data, using vertical and horizontal receivers, and utilize a realistic model of the seismic velocities at the Ecuadorian margin. We add perturbations within the subducting plates of 4×4 km and 2×2 km in P- and S-wave velocities, respectively, such that potential crosstalk between the two models can be identified. Our results show that the location and amplitude of the perturbations can be reconstructed in high quality down to approximately 70 km depth. We find that the inversion of the S-wave velocity prior to the inversion of the P-wave velocity is necessary to guarantee a good reconstruction of both models; however, the spatial resolution of the S-wave model is superior to the P-wave model. We also show that frequencies up to 1 Hz are sufficient to achieve high resolution. Further tests demonstrate how results depend on the accuracy of the estimated source orientation. Resulting models do not suffer in quality as artefacts near the source positions compensate for the inaccuracy of source orientation. If sources are located within the subducted plate instead of beneath, resulting models are comparable and the convergence of the inversion scheme is sped up. The accuracy of the source position within the model compared to the true earthquake location is critical and implies that earthquake relocation during the inversion process is necessary, in a similar way as in local earthquake traveltime tomography.

Amplitude-based misfit functions in viscoelastic full-waveform inversion applied to walk-away vertical seismic profile data

Geophysics ◽  
2019 ◽  
Vol 84 (5) ◽  
pp. B335-B351 ◽  
Author(s):  
Wenyong Pan ◽  
Kristopher A. Innanen
Keyword(s):  
Amplitude Ratio ◽  
Waveform Inversion ◽  
Seismic Profile ◽  
P Wave ◽  
Western Canada ◽  
Full Waveform Inversion ◽  
Walk Away ◽  
Profile Data ◽  
S Wave ◽  
Full Waveform

Viscoelastic full-waveform inversion is applied to walk-away vertical seismic profile data acquired at a producing heavy-oil field in Western Canada for the determination of subsurface velocity models (P-wave velocity [Formula: see text] and S-wave velocity [Formula: see text]) and attenuation models (P-wave quality factor [Formula: see text] and S-wave quality factor [Formula: see text]). To mitigate strong velocity-attenuation trade-offs, a two-stage approach is adopted. In Stage I, [Formula: see text] and [Formula: see text] models are first inverted using a standard waveform-difference (WD) misfit function. Following this, in Stage II, different amplitude-based misfit functions are used to estimate the [Formula: see text] and [Formula: see text] models. Compared to the traditional WD misfit function, the amplitude-based misfit functions exhibit stronger sensitivity to attenuation anomalies and appear to be able to invert [Formula: see text] and [Formula: see text] models more reliably in the presence of velocity errors. Overall, the root-mean-square amplitude-ratio and spectral amplitude-ratio misfit functions outperform other misfit function choices. In the final outputs of our inversion, significant drops in the [Formula: see text] to [Formula: see text] ratio (~1.6) and Poisson’s ratio (~0.23) are apparent within the Clearwater Formation (depth ~0.45–0.50 km) of the Mannville Group in the Western Canada Sedimentary Basin. Strong [Formula: see text] (~20) and [Formula: see text] (~15) anomalies are also evident in this zone. These observations provide information to help identify the target attenuative reservoir saturated with heavy-oil resources.

scholarly journals 3D elastic full-waveform inversion using P-wave excitation amplitude: Application to ocean bottom cable field data

Geophysics ◽  
2018 ◽  
Vol 83 (2) ◽  
pp. R129-R140 ◽  
Author(s):  
Ju-Won Oh ◽  
Mahesh Kalita ◽  
Tariq Alkhalifah
Keyword(s):  
Waveform Inversion ◽  
Excitation Amplitude ◽  
Maximum Energy ◽  
P Wave ◽  
Full Waveform Inversion ◽  
Ocean Bottom ◽  
Wave Excitation ◽  
S Wave ◽  
Full Waveform ◽  
Gradient Direction

We have developed an efficient elastic full-waveform inversion (FWI) based on the P-wave excitation amplitude (maximum energy arrival) approximation in the source wavefields. Because, based on the P-wave excitation approximation (ExA), the gradient direction is approximated by the crosscorrelation of source and receiver wavefields at only excitation time, it estimates the gradient direction faster than its conventional counterpart. In addition to this computational speedup, the P-wave ExA automatically ignores SP and SS correlations in the approximated gradient direction. In elastic FWI for ocean bottom cable (OBC) data, the descent direction for the S-wave velocity is often degraded by undesired long-wavelength features from the SS correlation. For this reason, the P-wave excitation approach increases the convergence rate of multiparameter FWI compared with the conventional approach. The modified 2D Marmousi model with OBC acquisition is used to verify the differences between the conventional method and ExA. Finally, the feasibility of the proposed method is demonstrated on a real OBC data from the North Sea.

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.

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.

Parameterization analysis and field validation of VTI-elastic full-waveform inversion in a walk-away vertical seismic profile configuration

Geophysics ◽  
2020 ◽  
Vol 85 (3) ◽  
pp. B87-B107 ◽  
Author(s):  
Wenyong Pan ◽  
Kristopher A. Innanen ◽  
Yanfei Wang
Keyword(s):  
Wave Velocity ◽  
Waveform Inversion ◽  
Seismic Profile ◽  
P Wave ◽  
Full Waveform Inversion ◽  
Walk Away ◽  
S Wave ◽  
Full Waveform ◽  
Trade Offs ◽  
Vsp Data

Elastic full-waveform inversion (FWI) in transversely isotropic media with a vertical symmetry axis (VTI) is applied to field walk-away vertical seismic profile (W-VSP) data acquired in Western Canada. The performance of VTI-elastic FWI is significantly influenced by the model parameterization choice. Synthetic analysis based on specific field survey configuration is carried out to evaluate three different VTI-elastic model parameterizations. Interparameter trade-offs are quantified using the recently introduced interparameter contamination sensitivity kernel approach. Synthetic results suggest that neglecting anisotropy leads to inaccurate velocity estimations. For the conventional vertical velocity-Thomsen’s parameter parameterization (i.e., vertical P-wave velocity, vertical S-wave velocity, Thomsen’s parameters [Formula: see text] and [Formula: see text]), a sequential inversion strategy is designed to reduce strong natural interparameter trade-offs. The model parameterizations of elastic-constant ([Formula: see text], [Formula: see text], [Formula: see text], and [Formula: see text]) and velocity-based (vertical, horizontal, and normal move-out P-wave velocities and vertical S-wave velocity) models appear to suffer from fewer interparameter trade-offs, providing more reliable velocity and anisotropy models. Results derived from application of VTI-elastic FWI to the field W-VSP data set tend to support the synthetic conclusions. Multiparameter point spread functions are calculated to quantify the local interparameter trade-offs of the inverted models. The output inversion results are interpreted to provide valuable references regarding the target hydrocarbon reservoir.

Joint petrophysical full-waveform inversion of the shallow-seismic and multi-offset GPR data

2021 ◽  
Author(s):  
Tan Qin ◽  
Thomas Bohlen ◽  
Yudi Pan
Keyword(s):  
Dielectric Permittivity ◽  
Wave Velocity ◽  
Waveform Inversion ◽  
P Wave ◽  
Relative Dielectric Permittivity ◽  
Full Waveform Inversion ◽  
S Wave ◽  
Petrophysical Parameters ◽  
Full Waveform ◽  
Shallow Seismic

<p>Shallow-seismic surface wave and ground penetrating radar (GPR) are employed in a wide range of engineering and geosciences applications. Full-waveform inversion (FWI) of either seismic or multi-offset GPR data are able to provide high-resolution subsurface characterization and have received particular attention in the past decade. Those two geophysical methods are involved in the increasing requirements of comprehensive site and material investigations. However, it is still challenging to provide an effective integration between seismic data and electromagnetic data. In this paper, we investigated the joint petrophysical inversion (JPI) of shallow-seismic and multi-offset GPR data for more consistent imaging of near surface. As a bridge between the seismic parameters (P-wave velocity, S-wave velocity, and density) and GPR parameters (relative dielectric permittivity and electric conductivity), the petrophysical relationships with the parameters namely porosity and saturation are employed to link two data sets. We first did a sensitivity analysis of the petrophysical parameters to the seismic and GPR parameters and then determined an efficient integration of using shallow-seismic FWI to update porosity and GPR FWI to update saturation, respectively. A comparison of several parameterisation combinations shows that the seismic velocity parameterisation in shallow-seismic FWI and a modified logarithm parameterisation in GPR FWI works well in reconstructing reliable S-wave velocity and relative dielectric permittivity models, respectively. With the help from the petrophysical links, we realized JPI by transforming those well recovered parameters to the petrophysical parameters and then to other seismic and GPR parameters. A synthetic test indicates that, compared with the individual petrophysical inversion and individual FWI, JPI outperforms in simultaneously reconstructing all seismic, GPR, and petrophysical parameters with higher resolution and improved details. It is proved that JPI would be a potential data integration approach for the shallow subsurface investigation.</p>

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.

scholarly journals Reflection Full Waveform Inversion with Decoupled Elastic-wave Equations in Inhomogeneous Medium

2021 ◽  
Vol 64 (1) ◽  
Author(s):  
Zhanyuan Liang ◽  
Guochen Wu ◽  
Xiaoyu Zhang ◽  
Qingyang Li
Keyword(s):  
Elastic Wave ◽  
Inhomogeneous Medium ◽  
Wave Velocity ◽  
Wave Equations ◽  
Waveform Inversion ◽  
P Wave ◽  
Full Waveform Inversion ◽  
S Wave ◽  
Full Waveform ◽  
Elastic Wave Equations

Reflection full-waveform inversion (RFWI) can reduce the nonlinearity of inversion providing an accurate initial velocity model for full-waveform inversion (FWI) through the tomographic components (low-wavenumber). However, elastic-wave reflection full-waveform inversion (ERFWI) is more vulnerable to the problem of local minimum due to the complicated multi-component wavefield. Our algorithm first divides kernels of ERFWI into four subkernels based on the theory of decoupled elastic-wave equations. Then we try to construct the tomographic components of ERFWI with only single-component wavefields, similarly to acoustic inversions. However, the S-wave velocity is still vulnerable to the coupling effects because of P-wave components contained in the S-wavefield in an inhomogeneous medium. Therefore we develop a method for decoupling elastic- wave equations in an inhomogeneous medium, which is applied to the decomposition of kernels in ERFWI. The new decoupled system can improve the accuracy of S-wavefield and hence further reduces the high-wavenumber crosstalk in the subkernel of S-wave velocity after kernels are decomposed. The numerical examples of Sigsbee2A model demonstrate that our ERFWI method with decoupled elastic-wave equations can efficiently recover the low-wavenumber components of the model and improve the precision of the S-wave velocity.

Tunnel detection at Yuma Proving Ground, Arizona, USA — Part 1: 2D full-waveform inversion experiment

Geophysics ◽  
2019 ◽  
Vol 84 (1) ◽  
pp. B95-B105 ◽  
Author(s):  
Yao Wang ◽  
Richard D. Miller ◽  
Shelby L. Peterie ◽  
Steven D. Sloan ◽  
Mark L. Moran ◽  
...  
Keyword(s):  
Surface Wave ◽  
Waveform Inversion ◽  
Velocity Profiles ◽  
Infinite Length ◽  
Full Waveform Inversion ◽  
S Wave ◽  
Near Surface ◽  
Tunnel Detection ◽  
Full Waveform ◽  
Wave Methods

We have applied time domain 2D full-waveform inversion (FWI) to detect a known 10 m deep wood-framed tunnel at Yuma Proving Ground, Arizona. The acquired seismic data consist of a series of 2D survey lines that are perpendicular to the long axis of the tunnel. With the use of an initial model estimated from surface wave methods, a void-detection-oriented FWI workflow was applied. A straightforward [Formula: see text] quotient masking method was used to reduce the inversion artifacts and improve confidence in identifying anomalies that possess a high [Formula: see text] ratio. Using near-surface FWI, [Formula: see text] and [Formula: see text] velocity profiles were obtained with void anomalies that are easily interpreted. The inverted velocity profiles depict the tunnel as a low-velocity anomaly at the correct location and depth. A comparison of the observed and simulated waveforms demonstrates the reliability of inverted models. Because the known tunnel has a uniform shape and for our purposes an infinite length, we apply 1D interpolation to the inverted [Formula: see text] profiles to generate a pseudo 3D (2.5D) volume. Based on this research, we conclude the following: (1) FWI is effective in near-surface tunnel detection when high resolution is necessary. (2) Surface-wave methods can provide accurate initial S-wave velocity [Formula: see text] models for near-surface 2D FWI.

3D elastic full-waveform inversion of surface waves in the presence of irregular topography using an envelope-based misfit function

Geophysics ◽  
2018 ◽  
Vol 83 (1) ◽  
pp. R1-R11 ◽  
Author(s):  
Dmitry Borisov ◽  
Ryan Modrak ◽  
Fuchun Gao ◽  
Jeroen Tromp
Keyword(s):  
Surface Wave ◽  
Objective Function ◽  
Surface Waves ◽  
Large Scale ◽  
Body Wave ◽  
Waveform Inversion ◽  
Full Waveform Inversion ◽  
S Wave ◽  
Near Surface ◽  
Full Waveform

Full-waveform inversion (FWI) is a powerful method for estimating the earth’s material properties. We demonstrate that surface-wave-driven FWI is well-suited to recovering near-surface structures and effective at providing S-wave speed starting models for use in conventional body-wave FWI. Using a synthetic example based on the SEG Advanced Modeling phase II foothills model, we started with an envelope-based objective function to invert for shallow large-scale heterogeneities. Then we used a waveform-difference objective function to obtain a higher-resolution model. To accurately model surface waves in the presence of complex tomography, we used a spectral-element wave-propagation solver. Envelope misfit functions are found to be effective at minimizing cycle-skipping issues in surface-wave inversions, and surface waves themselves are found to be useful for constraining complex near-surface features.

Sign in / Sign up

Export Citation Format

Share Document