An optimal transport distance for full-waveform inversion: Application to the 2014 Chevron benchmark data set

2016 ◽  
Author(s):  
Ludovic Metivier ◽  
Romain Brossier ◽  
Edouard Oudet ◽  
Quentin Mérigot ◽  
Jean Virieux
Keyword(s):  
Optimal Transport ◽  
Waveform Inversion ◽  
Full Waveform Inversion ◽  
Data Set ◽  
Benchmark Data ◽  
Transport Distance ◽  
Full Waveform

scholarly journals Increasing the robustness and applicability of full-waveform inversion: An optimal transport distance strategy

2016 ◽  
Vol 35 (12) ◽  
pp. 1060-1067 ◽  
Author(s):  
L. Métivier ◽  
R. Brossier ◽  
Q. Mérigot ◽  
E. Oudet ◽  
J. Virieux
Keyword(s):  
Optimal Transport ◽  
Waveform Inversion ◽  
Full Waveform Inversion ◽  
Transport Distance ◽  
Full Waveform

Optimal transport for mitigating cycle skipping in full-waveform inversion: A graph-space transform approach

Geophysics ◽  
2018 ◽  
Vol 83 (5) ◽  
pp. R515-R540 ◽  
Author(s):  
Ludovic Métivier ◽  
Aude Allain ◽  
Romain Brossier ◽  
Quentin Mérigot ◽  
Edouard Oudet ◽  
...  
Keyword(s):  
Seismic Data ◽  
Optimal Transport ◽  
Waveform Inversion ◽  
P Wave ◽  
Imaging Method ◽  
Convexity Property ◽  
Full Waveform Inversion ◽  
Low Frequencies ◽  
Transport Distance ◽  
Full Waveform

Optimal transport distance has been recently promoted as a tool to measure the discrepancy between observed and seismic data within the full-waveform-inversion strategy. This high-resolution seismic imaging method, based on a data-fitting procedure, suffers from the nonconvexity of the standard least-squares discrepancy measure, an issue commonly referred to as cycle skipping. The convexity of the optimal transport distance with respect to time shifts makes it a good candidate to provide a more convex misfit function. However, the optimal transport distance is defined only for the comparison of positive functions, while seismic data are oscillatory. A review of the different attempts proposed in the literature to overcome this difficulty is proposed. Their limitations are illustrated: Basically, the proposed strategies are either not applicable to real data, or they lose the convexity property of optimal transport. On this basis, we introduce a novel strategy based on the interpretation of the seismic data in the graph space. Each individual trace is considered, after discretization, as a set of Dirac points in a 2D space, where the amplitude becomes a geometric attribute of the data. This ensures the positivity of the data, while preserving the geometry of the signal. The differentiability of the misfit function is obtained by approximating the Dirac distributions through 2D Gaussian functions. The interest of this approach is illustrated numerically by computing misfit-function maps in schematic examples before moving to more realistic synthetic full-waveform exercises, including the Marmousi model. The better convexity of the graph-based optimal transport distance is shown. On the Marmousi model, starting from a 1D linearly increasing initial model, with data without low frequencies (no energy less than 3 Hz), a meaningful estimation of the P-wave velocity model is recovered, outperforming previously proposed optimal-transport-based misfit functions.

A new parameter set for anisotropic multiparameter full-waveform inversion and application to a North Sea data set

Geophysics ◽  
2016 ◽  
Vol 81 (4) ◽  
pp. U25-U38 ◽  
Author(s):  
Nuno V. da Silva ◽  
Andrew Ratcliffe ◽  
Vetle Vinje ◽  
Graham Conroy
Keyword(s):  
North Sea ◽  
Waveform Inversion ◽  
Real Data ◽  
Model Parameters ◽  
Full Waveform Inversion ◽  
Radiation Patterns ◽  
Data Set ◽  
Full Waveform ◽  
Seismic Image ◽  
Central North Sea

Parameterization lies at the center of anisotropic full-waveform inversion (FWI) with multiparameter updates. This is because FWI aims to update the long and short wavelengths of the perturbations. Thus, it is important that the parameterization accommodates this. Recently, there has been an intensive effort to determine the optimal parameterization, centering the fundamental discussion mainly on the analysis of radiation patterns for each one of these parameterizations, and aiming to determine which is best suited for multiparameter inversion. We have developed a new parameterization in the scope of FWI, based on the concept of kinematically equivalent media, as originally proposed in other areas of seismic data analysis. Our analysis is also based on radiation patterns, as well as the relation between the perturbation of this set of parameters and perturbation in traveltime. The radiation pattern reveals that this parameterization combines some of the characteristics of parameterizations with one velocity and two Thomsen’s parameters and parameterizations using two velocities and one Thomsen’s parameter. The study of perturbation of traveltime with perturbation of model parameters shows that the new parameterization is less ambiguous when relating these quantities in comparison with other more commonly used parameterizations. We have concluded that our new parameterization is well-suited for inverting diving waves, which are of paramount importance to carry out practical FWI successfully. We have demonstrated that the new parameterization produces good inversion results with synthetic and real data examples. In the latter case of the real data example from the Central North Sea, the inverted models show good agreement with the geologic structures, leading to an improvement of the seismic image and flatness of the common image gathers.

scholarly journals Normalized nonzero-lag crosscorrelation elastic full-waveform inversion

Geophysics ◽  
2019 ◽  
Vol 84 (1) ◽  
pp. R1-R10 ◽  
Author(s):  
Zhendong Zhang ◽  
Tariq Alkhalifah ◽  
Zedong Wu ◽  
Yike Liu ◽  
Bin He ◽  
...  
Keyword(s):  
Field Data ◽  
Extreme Values ◽  
Time Lag ◽  
Waveform Inversion ◽  
Optimal Time ◽  
Full Waveform Inversion ◽  
Data Set ◽  
Velocity Models ◽  
Full Waveform ◽  
The Earth

Full-waveform inversion (FWI) is an attractive technique due to its ability to build high-resolution velocity models. Conventional amplitude-matching FWI approaches remain challenging because the simplified computational physics used does not fully represent all wave phenomena in the earth. Because the earth is attenuating, a sample-by-sample fitting of the amplitude may not be feasible in practice. We have developed a normalized nonzero-lag crosscorrelataion-based elastic FWI algorithm to maximize the similarity of the calculated and observed data. We use the first-order elastic-wave equation to simulate the propagation of seismic waves in the earth. Our proposed objective function emphasizes the matching of the phases of the events in the calculated and observed data, and thus, it is more immune to inaccuracies in the initial model and the difference between the true and modeled physics. The normalization term can compensate the energy loss in the far offsets because of geometric spreading and avoid a bias in estimation toward extreme values in the observed data. We develop a polynomial-type weighting function and evaluate an approach to determine the optimal time lag. We use a synthetic elastic Marmousi model and the BigSky field data set to verify the effectiveness of the proposed method. To suppress the short-wavelength artifacts in the estimated S-wave velocity and noise in the field data, we apply a Laplacian regularization and a total variation constraint on the synthetic and field data examples, respectively.

A Graph-Space Optimal Transport Approach for Full Waveform Inversion

2018 ◽  
Author(s):  
L. Metivier ◽  
A. Allain ◽  
R. Brossier ◽  
Q. Merigot ◽  
E. Oudet ◽  
...  
Keyword(s):  
Optimal Transport ◽  
Waveform Inversion ◽  
Full Waveform Inversion ◽  
Full Waveform ◽  
Transport Approach

Tackling cycle skipping in full-waveform inversion with intermediate data

Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. R411-R427 ◽  
Author(s):  
Gang Yao ◽  
Nuno V. da Silva ◽  
Michael Warner ◽  
Di Wu ◽  
Chenhao Yang
Keyword(s):  
Objective Function ◽  
Waveform Inversion ◽  
Synthetic Data ◽  
Velocity Model ◽  
Full Waveform Inversion ◽  
Data Set ◽  
Intermediate Data ◽  
Full Waveform ◽  
L2 Norm ◽  
Recorded Data

Full-waveform inversion (FWI) is a promising technique for recovering the earth models for exploration geophysics and global seismology. FWI is generally formulated as the minimization of an objective function, defined as the L2-norm of the data residuals. The nonconvex nature of this objective function is one of the main obstacles for the successful application of FWI. A key manifestation of this nonconvexity is cycle skipping, which happens if the predicted data are more than half a cycle away from the recorded data. We have developed the concept of intermediate data for tackling cycle skipping. This intermediate data set is created to sit between predicted and recorded data, and it is less than half a cycle away from the predicted data. Inverting the intermediate data rather than the cycle-skipped recorded data can then circumvent cycle skipping. We applied this concept to invert cycle-skipped first arrivals. First, we picked up the first breaks of the predicted data and the recorded data. Second, we linearly scaled down the time difference between the two first breaks of each shot into a series of time shifts, the maximum of which was less than half a cycle, for each trace in this shot. Third, we moved the predicted data with the corresponding time shifts to create the intermediate data. Finally, we inverted the intermediate data rather than the recorded data. Because the intermediate data are not cycle-skipped and contain the traveltime information of the recorded data, FWI with intermediate data updates the background velocity model in the correct direction. Thus, it produces a background velocity model accurate enough for carrying out conventional FWI to rebuild the intermediate- and short-wavelength components of the velocity model. Our numerical examples using synthetic data validate the intermediate-data concept for tackling cycle skipping and demonstrate its effectiveness for the application to first arrivals.

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.

A graph-space approach to optimal transport for full-waveform inversion

2018 ◽  
Author(s):  
L. Métivier ◽  
A. Allain ◽  
R. Brossier ◽  
Q. Mérigot ◽  
E. Oudet ◽  
...  
Keyword(s):  
Optimal Transport ◽  
Waveform Inversion ◽  
Full Waveform Inversion ◽  
Full Waveform ◽  
Space Approach
Sign in / Sign up

Export Citation Format

Share Document