Parametric convolutional neural network-domain full-waveform inversion

Geophysics ◽  
2019 ◽  
Vol 84 (6) ◽  
pp. R881-R896 ◽  
Author(s):  
Yulang Wu ◽  
George A. McMechan

Most current full-waveform inversion (FWI) algorithms minimize the data residuals to estimate a velocity model based on the assumption that the updated model is the sum of a background model and an estimated model perturbation. We have performed reparameterization of the initial velocity model, by the weights in a convolutional neural network (CNN), to automatically capture the salient features in the initial model, as a priori information. The prior information in CNN weights is iteratively updated as regularization to constrain the CNN-domain inversion to refine the features captured in CNN pretraining by reducing the data misfit. Synthetic examples using a 1D increasing velocity function v(z) and a 2D smoothed version of the correct Marmousi2 model as initial models indicate that the performance of the CNN-domain FWI depends on the existence and accuracy of the prior information in the initial velocity model (i.e., whether features whose positions, shapes, and values are present in the correct model are approximately included in the initial model). Forty different sets of randomly initialized CNN weights are used to parameterize and test CNN-domain FWI, using a 2D smoothed Sigsbee model as the initial velocity model. All 40 sets invert for the Sigsbee salt body more accurately (with a smaller standard deviation of the final rms model errors), by CNN-domain FWI, than FWI does. Features that are not represented within the CNN hidden layers in the initial velocity model, and so cannot be recovered by CNN-domain FWI, can be recovered using the final CNN-domain FWI velocity model as the starting model in a subsequent conventional FWI.

Author(s):  
Ehsan Jamali Hondori ◽  
Chen Guo ◽  
Hitoshi Mikada ◽  
Jin-Oh Park

AbstractFull-waveform inversion (FWI) of limited-offset marine seismic data is a challenging task due to the lack of refracted energy and diving waves from the shallow sediments, which are fundamentally required to update the long-wavelength background velocity model in a tomographic fashion. When these events are absent, a reliable initial velocity model is necessary to ensure that the observed and simulated waveforms kinematically fit within an error of less than half a wavelength to protect the FWI iterative local optimization scheme from cycle skipping. We use a migration-based velocity analysis (MVA) method, including a combination of the layer-stripping approach and iterations of Kirchhoff prestack depth migration (KPSDM), to build an accurate initial velocity model for the FWI application on 2D seismic data with a maximum offset of 5.8 km. The data are acquired in the Japan Trench subduction zone, and we focus on the area where the shallow sediments overlying a highly reflective basement on top of the Cretaceous erosional unconformity are severely faulted and deformed. Despite the limited offsets available in the seismic data, our carefully designed workflow for data preconditioning, initial model building, and waveform inversion provides a velocity model that could improve the depth images down to almost 3.5 km. We present several quality control measures to assess the reliability of the resulting FWI model, including ray path illuminations, sensitivity kernels, reverse time migration (RTM) images, and KPSDM common image gathers. A direct comparison between the FWI and MVA velocity profiles reveals a sharp boundary at the Cretaceous basement interface, a feature that could not be observed in the MVA velocity model. The normal faults caused by the basal erosion of the upper plate in the study area reach the seafloor with evident subsidence of the shallow strata, implying that the faults are active.


2018 ◽  
Vol 22 (4) ◽  
pp. 291-300
Author(s):  
Sagar Singh ◽  
Ali Ismet Kanli ◽  
Sagarika Mukhopadhyay

This paper investigates the capability of acoustic Full Waveform Inversion (FWI) in building Marmousi velocity model, in time and frequency domain. FWI is an iterative minimization of misfit between observed and calculated data which is generally solved in three segments: forward modeling, which numerically solves the wave equation with an initial model, gradient computation of the objective function, and updating the model parameters, with a valid optimization method. FWI codes developed in MATLAB herein FWISIMAT (Full Waveform Inversion in Seismic Imaging using MATLAB) are successfully implemented using the Marmousi velocity model as the true model. An initial model is obtained by smoothing the true model to initiate FWI procedure. Smoothing ensures an adequate starting model for FWI, as the FWI procedure is known to be sensitive on the starting model. The final model is compared with the true model to review the number of recovered velocities. FWI codes developed in MATLAB herein FWISIMAT (Full Waveform Inversion in Seismic Imaging using MATLAB) are successfully implemented usingMarmousi velocity model astrue model. An initial model is derived from smoothing the true model to initiate FWI procedure. Smoothing ensures an adequate starting model for FWI, as the FWI procedure is known to be sensitive onstarting model. The final model is compared with the true model to review theamount of recovered velocities. 


Geophysics ◽  
2014 ◽  
Vol 79 (2) ◽  
pp. R55-R61 ◽  
Author(s):  
Tariq Alkhalifah ◽  
Yunseok Choi

In full-waveform inversion (FWI), a gradient-based update of the velocity model requires an initial velocity that produces synthetic data that are within a half-cycle, everywhere, from the field data. Such initial velocity models are usually extracted from migration velocity analysis or traveltime tomography, among other means, and are not guaranteed to adhere to the FWI requirements for an initial velocity model. As such, we evaluated an objective function based on the misfit in the instantaneous traveltime between the observed and modeled data. This phase-based attribute of the wavefield, along with its phase unwrapping characteristics, provided a frequency-dependent traveltime function that was easy to use and quantify, especially compared to conventional phase representation. With a strong Laplace damping of the modeled, potentially low-frequency, data along the time axis, this attribute admitted a first-arrival traveltime that could be compared with picked ones from the observed data, such as in wave equation tomography (WET). As we relax the damping on the synthetic and observed data, the objective function measures the misfit in the phase, however unwrapped. It, thus, provided a single objective function for a natural transition from WET to FWI. A Marmousi example demonstrated the effectiveness of the approach.


2020 ◽  
Author(s):  
Andrzej Górszczyk ◽  
Ludovic Métivier ◽  
Romain Brossier

<p>Investigations of the deep lithosphere aiming at the reconstruction of the geological models remain one of the key sources of the knowledge about the processes shaping the outer shell of our planet. Among different methods, the active seismic Ocean-Bottom Seismometer (OBS) experiments conducted in wide-angle configuration are routinely employed to better understand these processes. Indeed, long-offset seismic data, combined with computationally efficient travetime tomographic methods, have a great potential to constrain the macro-scale subsurface velocity models at large depths. </p><p>On the other hand, decades of development of acquisition systems, more and more efficient algorithms and high-performance computing resources make it now feasible to move beyond the regional raytracing-based traveltime tomography. In particular, the waveform inversion methods, such as Full-Waveform Inversion (FWI), are able to exhaustively exploit the rich information collected along the long-offset diving and refraction wavepaths, additionally enriched with the wide-angle reflection arrivals. So far however, only a few attempts have been conducted in the academic community to combine wide-angle seismic data with FWI for high-resolution crustal-scale velocity model reconstruction. This is partially due to the non-convexity of FWI misfit function, which increases with the complexity of the geological setting reflected by the seismograms. </p><p>In its classical form FWI is a nonlinear least-squares problem, which is solved through the local optimization techniques. This imposes the strong constraint on the accuracy of the starting FWI model. To avoid cycle-skipping problem the initial model must predict synthetic data within the maximum error of half-period time-shift with respect to the observed data. The criterion is difficult to fulfil when facing the crustal-scale FWI, because the long-offset acquisition translates to the long time of wavefront propagation and therefore accumulation of the traveltime error along the wavepath simulated in the initial model. This in turns increases the possibility of the cycle-skipping taking into account large number of propagated wavelengths.</p><p>Searching to mitigate this difficulty, here we investigate FWI with a Graph-Space Optimal Transport (GSOT) misfit function. Comparing to the classical least-squares norm, GSOT is convex with respect to the patterns in the waveform which can be shifted in time for more than half-period. Therefore, with proper data selection strategy GSOT misfit-function has potential to reduce the risk of cycle-skipping. We demonstrate the robustness of this novel approach using 2D wide-angle OBS data-set generated in a GO_3D_OBS synthetic model of subduction zone (30 km x 175 km). We show that using GSOT cost-function combined with the multiscale FWI strategy, we reconstruct in details the highly complex geological structure starting from a simple 1D velocity model. We believe that further developments of OT-based misfit functions can significantly reduce the constraints on the starting model accuracy and reduce the overall risk of cycle-skipping during FWI of wide-angle OBS data.</p>


Geophysics ◽  
2016 ◽  
Vol 81 (4) ◽  
pp. R173-R184 ◽  
Author(s):  
Angelo Sajeva ◽  
Mattia Aleardi ◽  
Eusebio Stucchi ◽  
Nicola Bienati ◽  
Alfredo Mazzotti

We have developed a stochastic full-waveform inversion that uses genetic algorithms (GA FWI) to estimate acoustic macro models of the P-wave velocity field. Stochastic methods such as GA severely suffer the curse of dimensionality, meaning that they require unaffordable computer resources for inverse problems with many unknowns and expensive forward modeling. To mitigate this issue, we have proposed a two-grid technique with a coarse grid to represent the subsurface for the GA inversion and a finer grid for the forward modeling. We have applied this procedure to invert synthetic acoustic data of the Marmousi model, and we have developed three different tests. The first two tests use a velocity model derived from standard stacking velocity analysis as prior information and differ only for the parameterization of the coarse grid. Their comparison indicates that a smart parameterization of the coarse grid may significantly improve the final result. The third test uses a linearly increasing 1D velocity model as prior information, a layer-stripping procedure, and a large number of model evaluations. All three tests return velocity models that fairly reproduce the long-wavelength structures of the Marmousi. First-break cycle skipping related to the seismograms of the final GA-FWI models is significantly reduced compared with that computed on the models used as prior information. Descent-based FWIs starting from final GA-FWI models yield velocity models with low and comparable model misfits and with an improved reconstruction of the structural details. The quality of the models obtained by GA FWI plus descent-based FWI is benchmarked against the models obtained by descent-based FWI started from a smoothed version of the Marmousi and started directly from the prior information models. Our results are promising and demonstrate the ability of the two-grid GA FWI to yield velocity models suitable as input to descent-based FWI.


2021 ◽  
Author(s):  
Ehsan Jamali Hondori ◽  
Chen Guo ◽  
Hitoshi Mikada ◽  
Jin-Oh Park

Abstract Full waveform inversion (FWI) of limited-offset marine seismic data is a challenging task due to the lack of refracted energy and diving waves from the shallow sediments, which are fundamentally required to update the long-wavelength background velocity model through a tomographic fashion. When these events are absent, a reliable initial velocity model is necessary to assure that the observed and simulated waveforms kinematically fit within an error less than half a wavelength to protect the FWI iterative local optimization scheme from cycle skipping. We use a migration-based velocity analysis (MVA) method, including a combination of the layer stripping approach and iterations of Kirchhoff prestack depth migration (KPSDM) to build an accurate initial velocity model for the FWI application on 2D seismic data with a maximum offset of 5.8 km. The data is acquired in the Japan Trench subduction zone, and we focus on the area where the shallow sediments overlying a highly reflective basement on top of the Cretaceous erosional unconformity are severely faulted and deformed. Despite the limited offsets available in the seismic data, our carefully designed workflow for data preconditioning, initial model building, and waveform inversion provides a velocity model which could improve the depth images down to a depth of almost 3.5 km. We present several quality control measures to assess the reliability of the resulting FWI model, including ray path illuminations, sensitivity kernels, reverse time migration (RTM) images, and KPSDM common image gathers. A direct comparison between the FWI and MVA velocity profiles reveals a sharp boundary at the Cretaceous basement interface, a feature which could not be observed in the MVA velocity model. The normal faults caused by the basal erosion of the upper plate in the study area reach the seafloor with evident subsidence of the shallow strata, implying that the faults are active.


Sign in / Sign up

Export Citation Format

Share Document