scholarly journals Accelerated full-waveform inversion using dynamic mini-batches

2020 ◽  
Vol 221 (2) ◽  
pp. 1427-1438 ◽  
Author(s):  
Dirk Philip van Herwaarden ◽  
Christian Boehm ◽  
Michael Afanasiev ◽  
Solvi Thrastarson ◽  
Lion Krischer ◽  
...  
Keyword(s):  
Waveform Inversion ◽  
Real Data ◽  
Complete Data ◽  
Control Group ◽  
Full Waveform Inversion ◽  
African Plate ◽  
Data Set ◽  
Full Waveform ◽  
Spatial Coverage ◽  
Quasi Newton

SUMMARY We present an accelerated full-waveform inversion based on dynamic mini-batch optimization, which naturally exploits redundancies in observed data from different sources. The method rests on the selection of quasi-random subsets (mini-batches) of sources, used to approximate the misfit and the gradient of the complete data set. The size of the mini-batch is dynamically controlled by the desired quality of the gradient approximation. Within each mini-batch, redundancy is minimized by selecting sources with the largest angular differences between their respective gradients, and spatial coverage is maximized by selecting candidate events with Mitchell’s best-candidate algorithm. Information from sources not included in a specific mini-batch is incorporated into each gradient calculation through a quasi-Newton approximation of the Hessian, and a consistent misfit measure is achieved through the inclusion of a control group of sources. By design, the dynamic mini-batch approach has several main advantages: (1) The use of mini-batches with adaptive size ensures that an optimally small number of sources is used in each iteration, thus potentially leading to significant computational savings; (2) curvature information is accumulated and exploited during the inversion, using a randomized quasi-Newton method; (3) new data can be incorporated without the need to re-invert the complete data set, thereby enabling an evolutionary mode of full-waveform inversion. We illustrate our method using synthetic and real-data inversions for upper-mantle structure beneath the African Plate. In these specific examples, the dynamic mini-batch approach requires around 20 per cent of the computational resources in order to achieve data and model misfits that are comparable to those achieved by a standard full-waveform inversion where all sources are used in each iteration.

Evolutionary full-waveform inversion with dynamic mini-batches

2020 ◽  
Author(s):  
Dirk-Philip van Herwaarden ◽  
Christian Boehm ◽  
Michael Afanasiev ◽  
Solvi Thrastarson ◽  
Lion Krischer ◽  
...  
Keyword(s):  
Time Windows ◽  
Regional Scale ◽  
Waveform Inversion ◽  
Control Group ◽  
Full Waveform Inversion ◽  
African Plate ◽  
Full Waveform ◽  
Spatial Coverage ◽  
Gradient Calculation ◽  
Quasi Newton

<p>We present an evolutionary full-waveform inversion based on dynamic mini-batch optimization, which naturally exploits redundancies in observed data from different sources and allows the model to evolve along with the amount of available information in the data.</p><p>Quasi-random subsets (mini-batches) of sources are used to approximate the misfit and the gradient of the complete dataset. The size of the mini-batch is dynamically controlled by the desired quality of the approximation of the full gradient. Within each mini-batch, redundancy is minimized by selecting sources with the largest angular differences between their respective gradients, and spatial coverage is maximized by selecting candidate events with Mitchell’s best-candidate algorithm. Information from sources included in a previous mini-batch is incorporated into each gradient calculation through a quasi-Newton approximation of the Hessian, and a consistent misfit measure is achieved through the inclusion of a control group of sources.</p><p>By design, the dynamic mini-batch approach has several main advantages: (1) The use of mini-batches with adaptive sizes minimizes the number of redundant simulations per iteration, thus potentially leading to significant computational savings. (2) Curvature information is accumulated and used during the inversion, using a stochastic quasi-Newton method. (3) Data from new events or different time windows can seamlessly be incorporated during the iterations, thereby enabling an evolutionary mode of full-waveform inversion.</p><p>To illustrate our method, we start an inversion for upper mantle structure beneath the African plate. Starting from a smooth 1-D background model for a dataset recorded in the years 1990 to 1995, we then sequentially add more and more recent data into the inversion and show how the model can evolve as a function of data coverage. The mini-batch sampling approach allows us to incorporate data from several hundred earthquakes without increasing the computational burden, thereby going significantly beyond previous regional-scale full-waveform inversions.</p>

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.

FULL-WAVEFORM INVERSION WITH SOURCE AND RECEIVER COUPLING EFFECTS CORRECTION

Geophysics ◽  
2021 ◽  
pp. 1-37
Author(s):  
Xinhai Hu ◽  
Wei Guoqi ◽  
Jianyong Song ◽  
Zhifang Yang ◽  
Minghui Lu ◽  
...  
Keyword(s):  
Waveform Inversion ◽  
Nonlinear Least Squares ◽  
Real Data ◽  
Model Parameters ◽  
Full Waveform Inversion ◽  
Linear Inversion ◽  
Near Surface ◽  
Model Parameter ◽  
Full Waveform ◽  
Coupling Factors

Coupling factors of sources and receivers vary dramatically due to the strong heterogeneity of near surface, which are as important as the model parameters for the inversion success. We propose a full waveform inversion (FWI) scheme that corrects for variable coupling factors while updating the model parameter. A linear inversion is embedded into the scheme to estimate the source and receiver factors and compute the amplitude weights according to the acquisition geometry. After the weights are introduced in the objective function, the inversion falls into the category of separable nonlinear least-squares problems. Hence, we could use the variable projection technique widely used in source estimation problem to invert the model parameter without the knowledge of source and receiver factors. The efficacy of the inversion scheme is demonstrated with two synthetic examples and one real data test.

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.

Targeted stacking of weak seismic phases for improving mantle discontinuity imaging using full-waveform inversion

2021 ◽  
Author(s):  
Maria Koroni ◽  
Andreas Fichtner
Keyword(s):  
Seismic Waves ◽  
Waveform Inversion ◽  
Real Data ◽  
Full Waveform Inversion ◽  
Finite Frequency ◽  
Adjoint Methods ◽  
Frequency Sensitivity ◽  
Full Waveform ◽  
Discontinuity Structure ◽  
Mantle Discontinuity

<p>This study is a continuation of our efforts to connect adjoint methods and full-waveform inversion to common beamforming techniques, widely used and developed for signal enhancement. Our approach is focusing on seismic waves traveling in the Earth's mantle, which are phases commonly used to image internal boundaries, being however quite difficult to observe in real data. The main goal is to accentuate precursor waves arriving in well-known times before some major phase. These waves generate from interactions with global discontinuities in the mantle, thus being the most sensitive seismic phases and therefore most suitable for better understanding of discontinuity seismic structure. </p><p>Our work is based on spectral-element wave propagation which allows us to compute exact synthetic waveforms and adjoint methods for the calculation of sensitivity kernels. These tools are the core of full-waveform inversion and by our efforts we aim to incorporate more parts of the waveform in such inversion schemes. We have shown that targeted stacking of good quality waveforms arriving from various directions highlights the weak precursor waves. It additionally makes their traveltime finite frequency sensitivity prominent. This shows that we can benefit from using these techniques and exploit rather difficult parts of the seismogram.  It was also shown that wave interference is not easily avoided, but coherent phases arriving before the main phase also stack well and show on the sensitivity kernels. This does not hamper the evaluation of waveforms, as in a misfit measurement process one can exploit more phases on the body wave parts of seismograms.</p><p>In this study, we go a step forward and present recent developments of the approach relating to the effects of noise and a real data experiment. Realistic noise is added to synthetic waveforms in order to assess the methodology in a more pragmatic scenario. The addition of noise shows that stacking of coherent seismic phases is still possible and the sensitivity kernels of their traveltimes are not largely distorted, the precursor waves contribute sufficiently to their traveltime finite-frequency sensitivity kernels.<br>Using a well-located seismic array, we apply the method to real data and try to examine the possibilities of using non-ideal waveforms to perform imaging of the mantle discontinuity structure on the specific areas. In order to make the most out of the dense array configuration, we try subgroups of receivers for the targeted stacking and by moving along the array we aim at creating a cluster of stacks. The main idea is to use the subgroups as single receivers and create an evaluation of seismic discontinuity structure using information from each stack belonging to a subgroup. <br>Ideally, we aim at improving the tomographic images of discontinuities of selected regions by exploiting weaker seismic waves, which are nonetheless very informative.</p>

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.

Which parameterization is suitable for acoustic vertical transverse isotropic full waveform inversion? Part 2: Synthetic and real data case studies from Valhall

Geophysics ◽  
2013 ◽  
Vol 78 (2) ◽  
pp. R107-R124 ◽  
Author(s):  
Yaser Gholami ◽  
Romain Brossier ◽  
Stéphane Operto ◽  
Vincent Prieux ◽  
Alessandra Ribodetti ◽  
...  
Keyword(s):  
Waveform Inversion ◽  
Real Data ◽  
Velocity Model ◽  
Oil Field ◽  
Full Waveform Inversion ◽  
Trade Off ◽  
Wave Speeds ◽  
Wide Aperture ◽  
Full Waveform ◽  
Transverse Isotropic

It is necessary to account for anisotropy in full waveform inversion (FWI) of wide-azimuth and wide-aperture seismic data in most geologic environments, for correct depth positioning of reflectors, and for reliable estimations of wave speeds as a function of the direction of propagation. In this framework, choosing a suitable anisotropic subsurface parameterization is a central issue in monoparameter and multiparameter FWI. This is because this parameterization defines the influence of each physical parameter class on the data as a function of the scattering angle, and hence the resolution of the parameter reconstruction, and on the potential trade-off between different parameter classes. We apply monoparameter and multiparameter frequency-domain acoustic vertical transverse isotropic FWI to synthetic and real wide-aperture data, representative of the Valhall oil field. We first show that reliable monoparameter FWI can be performed to build a high-resolution velocity model (for the vertical, the horizontal, or normal move-out velocity), provided that the background models of two Thomsen parameters describe the long wavelengths of the subsurface sufficiently accurately. Alternatively, we show the feasibility of the joint reconstruction of two wave speeds (e.g., the vertical and horizontal wave speeds) with limited trade-off effects, while Thomsen parameter [Formula: see text] is kept fixed during the inversion. The influence of the wave speeds on the data for a limited range of scattering angles when combined each other can, however, significantly hamper the resolution with which the two wave speeds are imaged. These conclusions inferred from the application to the real data are fully consistent with those inferred from the theoretical parameterization analysis of acoustic vertical transverse isotropic FWI performed in the companion report.

Sign in / Sign up

Export Citation Format

Share Document