Improvements to elastic full-waveform inversion using cross-gradient constraints

Geophysics ◽  
2018 ◽  
Vol 83 (2) ◽  
pp. R105-R115 ◽  
Author(s):  
Edgar Manukyan ◽  
Hansruedi Maurer ◽  
André Nuber

Seismic full-waveform inversion (FWI) is potentially a powerful method for obtaining high-resolution subsurface images, but the results are often distorted by nonlinear effects and parameter trade-offs. Such distortions can be particularly severe in the case of multiparameter FWI, such as elastic FWI, in which inversion is performed simultaneously for P- and S-wave velocities and density. The problem can be alleviated by adding constraints in the form of plausible a priori information. A usually well-justified constraint includes the structural similarity of different model parameters; i.e., an anomalous body likely exhibits variations in all elastic properties, although their magnitudes may be different. To consider such types of a priori information, we have developed a structurally constrained elastic FWI, which is based on minimization of the cross products of gradients of different model parameters. Our synthetic 2D experiments show that structurally constrained FWI can significantly improve model reconstruction. It is also demonstrated that our approach still leads to improved results, even when the structural similarity between the individual parameter types is not exactly met. Inversions of field data show that in comparison to conventional FWI, structurally constrained FWI is able to match the field data equally well while requiring less structural complexity of the subsurface.

Geophysics ◽  
2012 ◽  
Vol 77 (2) ◽  
pp. H19-H31 ◽  
Author(s):  
Knud Skou Cordua ◽  
Thomas Mejer Hansen ◽  
Klaus Mosegaard

We present a general Monte Carlo full-waveform inversion strategy that integrates a priori information described by geostatistical algorithms with Bayesian inverse problem theory. The extended Metropolis algorithm can be used to sample the a posteriori probability density of highly nonlinear inverse problems, such as full-waveform inversion. Sequential Gibbs sampling is a method that allows efficient sampling of a priori probability densities described by geostatistical algorithms based on either two-point (e.g., Gaussian) or multiple-point statistics. We outline the theoretical framework for a full-waveform inversion strategy that integrates the extended Metropolis algorithm with sequential Gibbs sampling such that arbitrary complex geostatistically defined a priori information can be included. At the same time we show how temporally and/or spatiallycorrelated data uncertainties can be taken into account during the inversion. The suggested inversion strategy is tested on synthetic tomographic crosshole ground-penetrating radar full-waveform data using multiple-point-based a priori information. This is, to our knowledge, the first example of obtaining a posteriori realizations of a full-waveform inverse problem. Benefits of the proposed methodology compared with deterministic inversion approaches include: (1) The a posteriori model variability reflects the states of information provided by the data uncertainties and a priori information, which provides a means of obtaining resolution analysis. (2) Based on a posteriori realizations, complicated statistical questions can be answered, such as the probability of connectivity across a layer. (3) Complex a priori information can be included through geostatistical algorithms. These benefits, however, require more computing resources than traditional methods do. Moreover, an adequate knowledge of data uncertainties and a priori information is required to obtain meaningful uncertainty estimates. The latter may be a key challenge when considering field experiments, which will not be addressed here.


Geophysics ◽  
2016 ◽  
Vol 81 (1) ◽  
pp. WA213-WA223 ◽  
Author(s):  
Lino Schmid ◽  
Jürg Schweizer ◽  
John Bradford ◽  
Hansruedi Maurer

Snow stratigraphy and liquid water content are key contributing factors to avalanche formation. Upward-looking ground-penetrating radar (upGPR) systems allow nondestructive monitoring of the snowpack, but deriving density and liquid water content profiles is not yet possible based on the direct analysis of the reflection response. We have investigated the feasibility of deducing these quantities using full-waveform inversion (FWI) techniques applied to upGPR data. For that purpose, we have developed a frequency-domain FWI algorithm in which we additionally took advantage of time-domain features such as the arrival times of reflected waves. Our results indicated that FWI applied to upGPR data is generally feasible. More specifically, we could show that in the case of a dry snowpack, it is possible to derive snow densities and layer thicknesses if sufficient a priori information is available. In case of a wet snowpack, in which it also needs to be inverted for the liquid water content, the algorithm might fail, even if sufficient a priori information is available, particularly in the presence of realistic noise. Finally, we have investigated the capability of FWI to resolve thin layers that play a key role in snow stability evaluation. Our simulations indicate that layers with thicknesses well below the GPR wavelengths can be identified, but in the presence of significant liquid water, the thin-layer properties may be prone to inaccuracies. These results are encouraging and motivate applications to field data, but significant issues remain to be resolved, such as the determination of the generally unknown upGPR source function and identifying the optimal number of layers in the inversion models. Furthermore, a relatively high level of prior knowledge is required to let the algorithm converge. However, we feel these are not insurmountable and the new technology has significant potential to improve field data analysis.


Geophysics ◽  
2014 ◽  
Vol 79 (6) ◽  
pp. R247-R261 ◽  
Author(s):  
Lisa Groos ◽  
Martin Schäfer ◽  
Thomas Forbriger ◽  
Thomas Bohlen

Full-waveform inversion (FWI) of Rayleigh waves is attractive for shallow geotechnical investigations due to the high sensitivity of Rayleigh waves to the S-wave velocity structure of the subsurface. In shallow-seismic field data, the effects of anelastic damping are significant. Dissipation results in a low-pass effect as well as frequency-dependent decay with offset. We found this by comparing recorded waveforms with elastic and viscoelastic wave simulation. The effects of anelastic damping must be considered in FWI of shallow-seismic Rayleigh waves. FWI using elastic simulation of wave propagation failed in synthetic inversion tests in which we tried to reconstruct the S-wave velocity in a viscoelastic model. To overcome this, [Formula: see text]-values can be estimated from the recordings to quantify viscoelasticity. Waveform simulation in the FWI then uses these a priori values when inferring seismic velocities and density. A source-wavelet correction, which is inevitable in FWI of field data, can compensate a significant fraction of the residuals between elastically and viscoelastically simulated data by narrowing the signals’ bandwidth. This way, elastic simulation becomes applicable in FWI of data from anelastic media. This approach, however, was not able to produce a frequency-dependent amplitude decay with offset. Reconstruction, therefore, was more accurate when using appropriate viscoelastic modeling in FWI of shallow-seismic Rayleigh waves. We found this by synthetic inversion tests using elastic forward simulation as well as viscoelastic simulation with different a priori values for [Formula: see text].


Geophysics ◽  
2016 ◽  
Vol 81 (4) ◽  
pp. U25-U38 ◽  
Author(s):  
Nuno V. da Silva ◽  
Andrew Ratcliffe ◽  
Vetle Vinje ◽  
Graham Conroy

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.


Geophysics ◽  
2021 ◽  
pp. 1-37
Author(s):  
Xinhai Hu ◽  
Wei Guoqi ◽  
Jianyong Song ◽  
Zhifang Yang ◽  
Minghui Lu ◽  
...  

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.


Geophysics ◽  
2019 ◽  
Vol 84 (1) ◽  
pp. R1-R10 ◽  
Author(s):  
Zhendong Zhang ◽  
Tariq Alkhalifah ◽  
Zedong Wu ◽  
Yike Liu ◽  
Bin He ◽  
...  

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.


Sign in / Sign up

Export Citation Format

Share Document