Full-waveform inversion (FWI) faces the persistent challenge of cycle skipping, which can result in stagnation of the iterative methods at uninformative models with poor data fit. Extended reformulations of FWI avoid cycle skipping through adding auxiliary parameters to the model so that a good data fit can be maintained throughout the inversion process. The volume-based matched source waveform inversion algorithm introduces source parameters by relaxing the location constraint of source energy: It is permitted to spread in space, while being strictly localized at time [Formula: see text]. The extent of source energy spread is penalized by weighting the source energy with distance from the survey source location. For transmission data geometry (crosswell, diving wave, etc.) and transparent (nonreflecting) acoustic models, this penalty function is stable with respect to the data-frequency content, unlike the standard FWI objective. We conjecture that the penalty function is actually convex over much larger region in model space than is the FWI objective. Several synthetic examples support this conjecture and suggest that the theoretical limitation to pure transmission is not necessary: The inversion method can converge to a solution of the inverse problem in the absence of low-frequency data from an inaccurate initial velocity model even when reflections and refractions are present in the data along with transmitted energy.
Full-waveform inversion produces highly resolved images of the subsurface and quantitative estimation of seismic wave velocity, provided that its initial model is kinematically accurate at the longest data wavelengths. If this initialization constraint is not satisfied, iterative model updating tends to stagnate at kinematically incorrect velocity models producing suboptimal images. The source-receiver extension overcomes this “cycle-skip” pathology by modeling each trace with its own proper source wavelet, permitting a good data fit throughout the inversion process. Because source wavelets should be constant (or vary systematically) across a shot gather, a measure of source trace dependence, for example, the mean square of the signature-deconvolved wavelet scaled by time lag, can be minimized to update the velocity model. For kinematically simple data, such measures of wavelet variance are mathematically equivalent to traveltime misfit. Thus, the model obtained by source-receiver extended inversion is close to that produced by traveltime tomography, even though the process uses no picked times. For more complex data, in which energy travels from source to receiver by multiple raypaths, Green’s function spectral notches may lead to slowly decaying trace-dependent wavelets with energy at time lags unrelated to traveltime error. Tikhonov regularization of the data-fitting problem suppresses these large-lag signals. Numerical examples suggest that this regularized formulation of source-receiver extended inversion is capable of recovering reasonably good velocity models from synthetic transmission and reflection data without stagnation at suboptimal models encountered by standard full-waveform inversion, but with essentially the same computational cost.
Source signature estimation from seismic data is a crucial ingredient for successful application of seismic migration and full-waveform inversion (FWI). If the starting velocity deviates from the target velocity, FWI method with on-the-fly source estimation may fail due to the cycle-skipping problem. We have developed a source-based extended waveform inversion method, by introducing additional parameters in the source function, to solve the FWI problem without the source signature as a priori. Specifically, we allow the point source function to be dependent on spatial and time variables. In this way, we can easily construct an extended source function to fit the recorded data by solving a source matching subproblem; hence, it is less prone to cycle skipping. A novel source focusing annihilator, defined as the distance function from the real source position, is used for penalizing the defocused energy in the extended source function. A close data fit avoiding the cycle-skipping problem effectively makes the new method less likely to suffer from local minima, which does not require extreme low-frequency signals in the data. Numerical experiments confirm that our method can mitigate cycle skipping in FWI and is robust against random noise.
No abstract
SUMMARYThis abstract presents a computationally efficient method to approximate the inverse of the Hessian or normal operator arising in a linearized inverse problem for constant density acoustics model of reflection seismology. Solution of the linearized inverse problem problem involves construction of an image via prestack depth migration, then correction of the image amplitudes via application of the inverse of the normal operator. The normal operator acts by dip-dependent scaling of the amplitudes of its input vector. This property permits us to efficiently approximate the normal operator, and its inverse, from the result of its application to a single input vector, for example the image, and thereby approximately solve the linearized inverse scattering problem. We validate the method on a 2D section of the Marmousi model to correct the amplitudes of the migrated image.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.