Articoli di riviste sul tema "Inversion full wave"
Segui questo link per vedere altri tipi di pubblicazioni sul tema: Inversion full wave.
Cita una fonte nei formati APA, MLA, Chicago, Harvard e in molti altri stili
Vedi i top-50 articoli di riviste per l'attività di ricerca sul tema "Inversion full wave".
Accanto a ogni fonte nell'elenco di riferimenti c'è un pulsante "Aggiungi alla bibliografia". Premilo e genereremo automaticamente la citazione bibliografica dell'opera scelta nello stile citazionale di cui hai bisogno: APA, MLA, Harvard, Chicago, Vancouver ecc.
Puoi anche scaricare il testo completo della pubblicazione scientifica nel formato .pdf e leggere online l'abstract (il sommario) dell'opera se è presente nei metadati.
Vedi gli articoli di riviste di molte aree scientifiche e compila una bibliografia corretta.
1
Luo, Yi, Yue Ma, Yan Wu, Hongwei Liu e Lei Cao. "Full-traveltime inversion". GEOPHYSICS 81, n. 5 (settembre 2016): R261—R274. http://dx.doi.org/10.1190/geo2015-0353.1.
Abstract (sommario):
Many previously published wave-equation-based methods, which attempt to automatically invert traveltime or kinematic information in seismic data or migrated gathers for smooth velocities, suffer a common and severe problem — the inversions are involuntarily and unconsciously hijacked by amplitude information. To overcome this problem, we have developed a new wave-equation-based traveltime inversion methodology, referred to as full-traveltime (i.e., fully dependent on traveltime) inversion (FTI), to automatically estimate a kinematically accurate velocity model from seismic data. The key idea of FTI is to make the inversion fully dependent on traveltime information, and thus prevent amplitude interference during inversion. Under the assumption that velocity perturbations cause only traveltime changes, we have derived the FTI method in the data and image domains, which are applicable to transmitted arrivals and reflected waves, respectively. FTI does not require an accurate initial velocity model or low-frequency seismic data. Synthetic and field data tests demonstrate that FTI produces satisfactory inversion results, even when using constant velocity models as initials.
2
Brossier, Romain, Stéphane Operto e Jean Virieux. "Seismic imaging of complex onshore structures by 2D elastic frequency-domain full-waveform inversion". GEOPHYSICS 74, n. 6 (novembre 2009): WCC105—WCC118. http://dx.doi.org/10.1190/1.3215771.
Abstract (sommario):
Quantitative imaging of the elastic properties of the subsurface at depth is essential for civil engineering applications and oil- and gas-reservoir characterization. A realistic synthetic example provides for an assessment of the potential and limits of 2D elastic full-waveform inversion (FWI) of wide-aperture seismic data for recovering high-resolution P- and S-wave velocity models of complex onshore structures. FWI of land data is challenging because of the increased nonlinearity introduced by free-surface effects such as the propagation of surface waves in the heterogeneous near-surface. Moreover, the short wavelengths of the shear wavefield require an accurate S-wave velocity starting model if low frequencies are unavailable in the data. We evaluated different multiscale strategies with the aim of mitigating the nonlinearities. Massively parallel full-waveform inversion was implemented in the frequency domain. The numerical optimization relies on a limited-memory quasi-Newton algorithm thatoutperforms the more classic preconditioned conjugate-gradient algorithm. The forward problem is based upon a discontinuous Galerkin (DG) method on triangular mesh, which allows accurate modeling of free-surface effects. Sequential inversions of increasing frequencies define the most natural level of hierarchy in multiscale imaging. In the case of land data involving surface waves, the regularization introduced by hierarchical frequency inversions is not enough for adequate convergence of the inversion. A second level of hierarchy implemented with complex-valued frequencies is necessary and provides convergence of the inversion toward acceptable P- and S-wave velocity models. Among the possible strategies for sampling frequencies in the inversion, successive inversions of slightly overlapping frequency groups is the most reliable when compared to the more standard sequential inversion of single frequencies. This suggests that simultaneous inversion of multiple frequencies is critical when considering complex wave phenomena.
3
Zhang, Chao, Ting Lei e Yi Wang. "Two-Dimensional Full-Waveform Joint Inversion of Surface Waves Using Phases and Z/H Ratios". Applied Sciences 11, n. 15 (22 luglio 2021): 6712. http://dx.doi.org/10.3390/app11156712.
Abstract (sommario):
Surface-wave dispersion and the Z/H ratio are important parameters used to resolve the Earth’s structure, especially for S-wave velocity. Several previous studies have explored using joint inversion of these two datasets. However, all of these studies used a 1-D depth-sensitivity kernel, which lacks precision when the structure is laterally heterogeneous. Adjoint tomography (i.e., full-waveform inversion) is a state-of-the-art imaging method with a high resolution. It can obtain better-resolved lithospheric structures beyond the resolving ability of traditional ray-based travel-time tomography. In this study, we present a systematic investigation of the 2D sensitivities of the surface wave phase and Z/H ratio using the adjoint-state method. The forward-modeling experiments indicated that the 2D phase and Z/H ratio had different sensitivities to the S-wave velocity. Thus, a full-waveform joint-inversion scheme of surface waves with phases and a Z/H ratio was proposed to take advantage of their complementary sensitivities to the Earth’s structure. Both applications to synthetic data sets in large- and small-scale inversions demonstrated the advantage of the joint inversion over the individual inversions, allowing for the creation of a more unified S-wave velocity model. The proposed joint-inversion scheme offers a computationally efficient and inexpensive alternative to imaging fine-scale shallow structures beneath a 2D seismic array.
4
Yilmaz, Öz, Kai Gao, Milos Delic, Jianghai Xia, Lianjie Huang, Hossein Jodeiri e Andre Pugin. "A reality check on full-wave inversion applied to land seismic data for near-surface modeling". Leading Edge 41, n. 1 (gennaio 2022): 40–46. http://dx.doi.org/10.1190/tle41010040.1.
Abstract (sommario):
We evaluate the performance of traveltime tomography and full-wave inversion (FWI) for near-surface modeling using the data from a shallow seismic field experiment. Eight boreholes up to 20-m depth have been drilled along the seismic line traverse to verify the accuracy of the P-wave velocity-depth model estimated by seismic inversion. The velocity-depth model of the soil column estimated by traveltime tomography is in good agreement with the borehole data. We used the traveltime tomography model as an initial model and performed FWI. Full-wave acoustic and elastic inversions, however, have failed to converge to a velocity-depth model that desirably should be a high-resolution version of the model estimated by traveltime tomography. Moreover, there are significant discrepancies between the estimated models and the borehole data. It is understandable why full-wave acoustic inversion would fail — land seismic data inherently are elastic wavefields. The question is: Why does full-wave elastic inversion also fail? The strategy to prevent full-wave elastic inversion of vertical-component geophone data trapped in a local minimum that results in a physically implausible near-surface model may be cascaded inversion. Specifically, we perform traveltime tomography to estimate a P-wave velocity-depth model for the near-surface and Rayleigh-wave inversion to estimate an S-wave velocity-depth model for the near-surface, then use the resulting pairs of models as the initial models for the subsequent full-wave elastic inversion. Nonetheless, as demonstrated by the field data example here, the elastic-wave inversion yields a near-surface solution that still is not in agreement with the borehole data. Here, we investigate the limitations of FWI applied to land seismic data for near-surface modeling.
5
Tran, Khiem T., Michael McVay, Michael Faraone e David Horhota. "Sinkhole detection using 2D full seismic waveform tomography". GEOPHYSICS 78, n. 5 (1 settembre 2013): R175—R183. http://dx.doi.org/10.1190/geo2013-0063.1.
Abstract (sommario):
We have developed an application of 2D time-domain waveform tomography for detection of embedded sinkholes and anomalies. The measured seismic surface wavefields were inverted using a full-waveform inversion (FWI) technique, based on a finite-difference solution of 2D elastic wave equations and the Gauss-Newton inversion method. The key advantage of this approach is the ability to generate all possible wave propagation modes of seismic wavefields (body waves and Rayleigh waves) that are then compared with measured data to infer complex subsurface properties.The pressure-wave (P-wave) and shear-wave (S-wave) velocities are inverted independently and simultaneously. The FWI was applied to one synthetic and two real experimental data sets. The inversion results of synthetic data showed the useful capability of the waveform analysis in identifying an embedded void. The inversion results of real data sets showed that the waveform analysis was able to delineate (1) an embedded concrete culvert and (2) a complex profile with an embedded void and highly variable bedrock laterally and vertically. An independent invasive test (standard penetration test) was also conducted to verify the seismic test results.
6
Barnes, Christophe, e Marwan Charara. "The domain of applicability of acoustic full-waveform inversion for marine seismic data". GEOPHYSICS 74, n. 6 (novembre 2009): WCC91—WCC103. http://dx.doi.org/10.1190/1.3250269.
Abstract (sommario):
Marine reflection seismic data inversion is a compute-intensive process, especially in three dimensions. Approximations often are made to limit the number of physical parameters we invert for, or to speed up the forward modeling. Because the data often are dominated by unconverted P-waves, one popular approximation is to consider the earth as purely acoustic, i.e., no shear modulus. The material density sometimes is taken as a constant. Nonlinear waveform seismic inversion consists of iteratively minimizing the misfit between the amplitudes of the measured and the modeled data. Approximations, such as assuming an acoustic medium, lead to incorrect modeling of the amplitudes of the seismic waves, especially with respect to amplitude variation with offset (AVO), and therefore have a direct impact on the inversion results. For evaluation purposes, we have performed a series of inversions with different approximations and different constraints whereby the synthetic data set to recover is computed for a 1D elastic medium. A series of numerical experiments, although simple, help to define the applicability domain of the acoustic assumption. Acoustic full-wave inversion is applicable only when the S-wave velocity and the density fields are smooth enough to reduce the AVO effect, or when the near-offset seismograms are inverted with a good starting model. However, in many realistic cases, acoustic approximation penalizes the full-wave inversion of marine reflection seismic data in retrieving the acoustic parameters.
7
Biondi, Biondo, e Ali Almomin. "Simultaneous inversion of full data bandwidth by tomographic full-waveform inversion". GEOPHYSICS 79, n. 3 (1 maggio 2014): WA129—WA140. http://dx.doi.org/10.1190/geo2013-0340.1.
Abstract (sommario):
The convergence of full-waveform inversion can be improved by extending the velocity model along either the subsurface-offset axis or the time-lag axis. The extension of the velocity model along the time-lag axis enables us to linearly model large time shifts caused by velocity perturbations. This linear modeling was based on a new linearization of the scalar wave equation in which perturbation of the extended slowness squared was convolved in time with the second time derivative of the background wavefield. The linearization was accurate for reflected events and transmitted events. We determined that it can effectively model conventional reflection data as well as modern long-offset data containing diving waves. It also enabled the simultaneous inversion of reflections and diving waves, even when the starting velocity model was far from being accurate. We solved the optimization problem related to the inversion with a nested algorithm. The inner iterations were based on the proposed linearization and on a mixing of scales between the short- and long-wavelength components of the velocity model. Numerical tests performed on synthetic data modeled on the Marmousi model and on the “Caspian Sea” portion of the well-known BP model demonstrated the global-convergence properties as well as the high-resolution potential of the proposed method.
8
da Silva, Nuno V., Gang Yao e Michael Warner. "Semiglobal viscoacoustic full-waveform inversion". GEOPHYSICS 84, n. 2 (1 marzo 2019): R271—R293. http://dx.doi.org/10.1190/geo2017-0773.1.
Abstract (sommario):
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.
9
Luo, Y., e G. T. Schuster. "Wave‐equation traveltime inversion". GEOPHYSICS 56, n. 5 (maggio 1991): 645–53. http://dx.doi.org/10.1190/1.1443081.
Abstract (sommario):
This paper presents a new traveltime inversion method based on the wave equation. In this new method, designated as wave‐equation traveltime inversion (WT), seismograms are computed by any full‐wave forward modeling method (we use a finite‐difference method). The velocity model is perturbed until the traveltimes from the synthetic seismograms are best fitted to the observed traveltimes in a least squares sense. A gradient optimization method is used and the formula for the Frechét derivative (perturbation of traveltimes with respect to velocity) is derived directly from the wave equation. No traveltime picking or ray tracing is necessary, and there are no high frequency assumptions about the data. Body wave, diffraction, reflection and head wave traveltimes can be incorporated into the inversion. In the high‐frequency limit, WT inversion reduces to ray‐based traveltime tomography. It can also be shown that WT inversion is approximately equivalent to full‐wave inversion when the starting velocity model is “close” to the actual model. Numerical simulations show that WT inversion succeeds for models with up to 80 percent velocity contrasts compared to the failure of full‐wave inversion for some models with no more than 10 percent velocity contrast. We also show that the WT method succeeds in inverting a layered velocity model where a shooting ray‐tracing method fails to compute the correct first arrival times. The disadvantage of the WT method is that it appears to provide less model resolution compared to full‐wave inversion, but this problem can be remedied by a hybrid traveltime + full‐wave inversion method (Luo and Schuster, 1989).
10
Dettmer, Jan, Stan E. Dosso e Charles W. Holland. "Full wave-field reflection coefficient inversion". Journal of the Acoustical Society of America 122, n. 6 (dicembre 2007): 3327–37. http://dx.doi.org/10.1121/1.2793609.
11
Xue, Zhiguang, Junzhe Sun, Sergey Fomel e Tieyuan Zhu. "Accelerating full-waveform inversion with attenuation compensation". GEOPHYSICS 83, n. 1 (1 gennaio 2018): A13—A20. http://dx.doi.org/10.1190/geo2017-0469.1.
Abstract (sommario):
The calculation of the gradient in full-waveform inversion (FWI) usually involves crosscorrelating the forward-propagated source wavefield and the back-propagated data residual wavefield at each time step. In the real earth, propagating waves are typically attenuated due to the viscoelasticity, which results in an attenuated gradient for FWI. Replacing the attenuated true gradient with a [Formula: see text]-compensated gradient can accelerate the convergence rate of the inversion process. We have used a phase-dispersion and an amplitude-loss decoupled constant-[Formula: see text] wave equation to formulate a viscoacoustic FWI. We used this wave equation to generate a [Formula: see text]-compensated gradient, which recovers amplitudes while preserving the correct kinematics. We construct an exact adjoint operator in a discretized form using the low-rank wave extrapolation technique, and we implement the gradient compensation by reversing the sign of the amplitude-loss term in the forward and adjoint operators. This leads to a [Formula: see text]-dependent gradient preconditioning method. Using numerical tests with synthetic data, we demonstrate that the proposed viscoacoustic FWI using a constant-[Formula: see text] wave equation is capable of producing high-quality velocity models, and our [Formula: see text]-compensated gradient accelerates its convergence rate.
12
Bleibinhaus, Florian, e Stéphane Rondenay. "Effects of surface scattering in full-waveform inversion". GEOPHYSICS 74, n. 6 (novembre 2009): WCC69—WCC77. http://dx.doi.org/10.1190/1.3223315.
Abstract (sommario):
In full-waveform inversion of seismic body waves, often the free surface is ignored on grounds of computational efficiency. A synthetic study was performed to investigate the effects of this simplification. In terms of size and frequency, the test model and data conform to a real long-offset survey of the upper crust across the San Andreas fault. Random fractal variations are superimposed on a background model with strong lateral and vertical velocity variations ranging from 1200 to 6800 m/s. Synthetic data were computed and inverted for this model and different topographies. A fully viscoelastic time-domain code was used to synthesize the seismograms, and a viscoacoustic frequency-domain code was utilized to invert them. The inversion was focused on early arrivals, which are dominated by P-waves but also contain strong P-Rayleigh wave conversions from the near-field of the receiver. Resulting waveform models show artifacts and a loss of resolution from neglecting the free surface in the inversion, but the inversions are stable, and they still improve the resolution of kinematic models. The extent of deterioration depends more on the subsurface than on the surface structure. Inversion results were improved at no additional expense by introducing a weak contrast along a staircase function above shots and receivers.
13
Guddati, Murthy, Tuhin Roy, Abdelrahman M. Elmeliegy e Matthew W. Urban. "Shear wave elastography: From dispersion matching to full waveform inversion". Journal of the Acoustical Society of America 153, n. 3_supplement (1 marzo 2023): A265. http://dx.doi.org/10.1121/10.0018796.
Abstract (sommario):
Shear Wave Elastography (SWE) involves estimating mechanical properties through inversion, i.e., matching measured and simulated propagation characteristics of shear waves in the tissue. The accuracy of the estimated properties depends significantly on the specific characteristics/responses that are being matched. These could range from simple group velocity to dispersion curves and to full-wave response (particle velocity measurements). Using specific applications of arterial, liver, and tumor elstography, we illustrate that effective SWE is performed by resorting to an inversion approach, or combination of inversion approaches, guided by the underlying physics. To this end, we present inversion approaches ranging from matching dispersion characteristics to matching full waveform responses and provide rationale for choosing the appropriate technique(s) depending on the problem at hand.
14
Vigh, Denes, e E. William Starr. "3D prestack plane-wave, full-waveform inversion". GEOPHYSICS 73, n. 5 (settembre 2008): VE135—VE144. http://dx.doi.org/10.1190/1.2952623.
Abstract (sommario):
Prestack depth migration has been used for decades to derive velocity distributions in depth. Numerous tools and methodologies have been developed to reach this goal. Exploration in geologically more complex areas exceeds the abilities of existing methods. New data-acquisition and data-processing methods are required to answer these new challenges effectively. The recently introduced wide-azimuth data acquisition method offers better illumination and noise attenuation as well as an opportunity to more accurately determine velocities for imaging. One of the most advanced tools for depth imaging is full-waveform inversion. Prestack seismic full-waveform inversion is very challenging because of the nonlinearity and nonuniqueness of the solution. Combined with multiple iterations of forward modeling and residual wavefield back propagation, the method is computer intensive, especially for 3D projects. We studied a time-domain, plane-wave implementation of 3D waveform inversion. We found that plane-wave gathers are an attractive input to waveform inversion with dramatically reduced computer run times compared to traditional shot-gather approaches. The study was conducted on two synthetic data sets — Marmousi2 and SMAART Pluto 1.5 — and a field data set. The results showed that a velocity field can be reconstructed well using a multiscale time-domain implementation of waveform inversion. Although the time-domain solution does not take advantage of wavenumber redundancy, the method is feasible on current computer architectures for 3D surveys. The inverted velocity volume produces a quality image for exploration geologists by using numerous iterations of waveform inversion.
15
Sears, Timothy J., Penny J. Barton e Satish C. Singh. "Elastic full waveform inversion of multicomponent ocean-bottom cable seismic data: Application to Alba Field, U. K. North Sea". GEOPHYSICS 75, n. 6 (novembre 2010): R109—R119. http://dx.doi.org/10.1190/1.3484097.
Abstract (sommario):
Elastic full waveform inversion of multichannel seismic data represents a data-driven form of analysis leading to direct quantification of the subsurface elastic parameters in the depth domain. Previous studies have focused on marine streamer data using acoustic or elastic inversion schemes for the inversion of P-wave data. In this paper, P- and S-wave velocities are inverted for using wide-angle multicomponent ocean-bottom cable (OBC) seismic data. Inversion is undertaken using a two-dimensional elastic algorithm operating in the time domain, which allows accurate modeling and inversion of the full elastic wavefield, including P- and mode-converted PS-waves and their respective amplitude variation with offset (AVO) responses. Results are presented from the application of this technique to an OBC seismic data set from the Alba Field, North Sea. After building an initial velocity model and extracting a seismic wavelet, the data are inverted instages. In the first stage, the intermediate wavelength P-wave velocity structure is recovered from the wide-angle data and then the short-scale detail from near-offset data using P-wave data on the [Formula: see text] (vertical geophone) component. In the second stage, intermediate wavelengths of S-wave velocity are inverted for, which exploits the information captured in the P-wave’s elastic AVO response. In the third stage, the earlier models are built on to invert mode-converted PS-wave events on the [Formula: see text] (horizontal geophone) component for S-wave velocity, targeting first shallow and then deeper structure. Inversion of [Formula: see text] alone has been able to delineate the Alba Field in P- and S-wave velocity, with the main field and outlier sands visible on the 2D results. Inversion of PS-wave data has demonstrated the potential of using converted waves to resolve shorter wavelength detail. Even at the low frequencies [Formula: see text] inverted here, improved spatial resolution was obtained by inverting S-wave data compared with P-wave data inversion results.
16
Witte, Philipp, Mathias Louboutin, Keegan Lensink, Michael Lange, Navjot Kukreja, Fabio Luporini, Gerard Gorman e Felix J. Herrmann. "Full-waveform inversion, Part 3: Optimization". Leading Edge 37, n. 2 (febbraio 2018): 142–45. http://dx.doi.org/10.1190/tle37020142.1.
Abstract (sommario):
This tutorial is the third part of a full-waveform inversion (FWI) tutorial series with a step-by-step walkthrough of setting up forward and adjoint wave equations and building a basic FWI inversion framework. For discretizing and solving wave equations, we use Devito ( http://www.opesci.org/devito-public ), a Python-based domain-specific language for automated generation of finite-difference code ( Lange et al., 2016 ). The first two parts of this tutorial ( Louboutin et al., 2017 , 2018 ) demonstrated how to solve the acoustic wave equation for modeling seismic shot records and how to compute the gradient of the FWI objective function using the adjoint-state method. With these two key ingredients, we will now build an inversion framework that can be used to minimize the FWI least-squares objective function.
17
Ma, Yong, e Dave Hale. "Wave-equation reflection traveltime inversion with dynamic warping and full-waveform inversion". GEOPHYSICS 78, n. 6 (1 novembre 2013): R223—R233. http://dx.doi.org/10.1190/geo2013-0004.1.
Abstract (sommario):
In reflection seismology, full-waveform inversion (FWI) can generate high-wavenumber subsurface velocity models but often suffers from an objective function with local minima caused mainly by the absence of low frequencies in seismograms. These local minima cause cycle skipping when the low-wavenumber component in the initial velocity model for FWI is far from the true model. To avoid cycle skipping, we discovered a new wave-equation reflection traveltime inversion (WERTI) to update the low-wavenumber component of the velocity model, while using FWI to only update high-wavenumber details of the model. We implemented the low- and high-wavenumber inversions in an alternating way. In WERTI, we used dynamic image warping (DIW) to estimate the time shifts between recorded data and synthetic data. When compared with correlation-based techniques often used in traveltime estimation, DIW can avoid cycle skipping and estimate the time shifts accurately, even when shifts vary rapidly. Hence, by minimizing traveltime shifts estimated by dynamic warping, WERTI reduces errors in reflection traveltime inversion. Then, conventional FWI uses the low-wavenumber component estimated by WERTI as a new initial model and thereby refines the model with high-wavenumber details. The alternating combination of WERTI and FWI mitigates the velocity-depth ambiguity and can recover subsurface velocities using only high-frequency reflection data.
18
Agudo, Òscar Calderón, Nuno Vieira da Silva, George Stronge e Michael Warner. "Mitigating elastic effects in marine 3-D full-waveform inversion". Geophysical Journal International 220, n. 3 (18 dicembre 2019): 2089–104. http://dx.doi.org/10.1093/gji/ggz569.
Abstract (sommario):
SUMMARY The potential of full-waveform inversion (FWI) to recover high-resolution velocity models of the subsurface has been demonstrated in the last decades with its application to field data. But in certain geological scenarios, conventional FWI using the acoustic wave equation fails in recovering accurate models due to the presence of strong elastic effects, as the acoustic wave equation only accounts for compressional waves. This becomes more critical when dealing with land data sets, in which elastic effects are generated at the source and recorded directly by the receivers. In marine settings, in which sources and receivers are typically within the water layer, elastic effects are weaker but can be observed most easily as double mode conversions and through their effect on P-wave amplitudes. Ignoring these elastic effects can have a detrimental impact on the accuracy of the recovered velocity models, even in marine data sets. Ideally, the elastic wave equation should be used to model wave propagation, and FWI should aim to recover anisotropic models of velocity for P waves (vp) and S waves (vs). However, routine three-dimensional elastic FWI is still commercially impractical due to the elevated computational cost of modelling elastic wave propagation in regions with low S-wave velocity near the seabed. Moreover, elastic FWI using local optimization methods suffers from cross-talk between different inverted parameters. This generally leads to incorrect estimation of subsurface models, requiring an estimate of vp/vs that is rarely known beforehand. Here we illustrate how neglecting elasticity during FWI for a marine field data set that contains especially strong elastic heterogeneities can lead to an incorrect estimation of the P-wave velocity model. We then demonstrate a practical approach to mitigate elastic effects in 3-D yielding improved estimates, consisting of using a global inversion algorithm to estimate a model of vp/vs, employing matching filters to remove elastic effects from the field data, and performing acoustic FWI of the resulting data set. The quality of the recovered models is assessed by exploring the continuity of the events in the migrated sections and the fit of the latter with the recovered velocity model.
19
Wen, Jiawei, Can Jiang e Hao Chen. "High-Precision Corrosion Detection via SH1 Guided Wave Based on Full Waveform Inversion". Sensors 23, n. 24 (18 dicembre 2023): 9902. http://dx.doi.org/10.3390/s23249902.
Abstract (sommario):
Corrosion detection for industrial settings is crucial for safe and efficient operations. Due to its high imaging resolution, the guided–wave full–waveform inversion tomography technique has significant potential for corrosion detection of plate metals. Limited by the long wavelengths of A0 and S0 mode waves, this method exhibits inadequate detection resolution for the earlier shallow and small corrosion defects. Based on the relatively short wavelength characteristics of the SH1 mode wave, we propose a high–precision corrosion detection method via SH1 guided wave using the full waveform inversion algorithms. By conducting finite element simulations of ultrasonic–guided waves on aluminum plates with varying corrosion defects, a comparison was made to assess the detection precision across A0, S0, and SH1 modes. The comparison results showed that, whether for regular or irregular defects, the SH1 mode wave always exhibited higher imaging accuracy than the A0 and S0 mode waves for shallow and small–sized defects. The corresponding experiments were conducted on an aluminum plate with simple or complex defects. The results of the experiments reconfirmed that the full waveform inversion method using SH1 guided wave can effectively reconstruct the shape and size of small and shallow corrosion defects within aluminum plates.
20
Kwon, Taekhyun, Joongmoo Byun, Byoung Yeop Kim e Sik Huh. "Elastic full waveform inversion using plane-wave". Journal of Applied Geophysics 170 (novembre 2019): 103826. http://dx.doi.org/10.1016/j.jappgeo.2019.103826.
21
Sjögreen, Björn, e N. Anders Petersson. "Source Estimation by Full Wave Form Inversion". Journal of Scientific Computing 59, n. 1 (7 agosto 2013): 247–76. http://dx.doi.org/10.1007/s10915-013-9760-6.
22
Oh, Ju-Won, Youngjae Shin, Tariq Alkhalifah e Dong-Joo Min. "Multistage elastic full-waveform inversion for tilted transverse isotropic media". Geophysical Journal International 223, n. 1 (20 giugno 2020): 57–76. http://dx.doi.org/10.1093/gji/ggaa295.
Abstract (sommario):
SUMMARY Seismic anisotropy is an important physical phenomenon that significantly affects wave propagation in complex sedimentary basins. When geological structures exhibit steep dips or severe folding, the symmetry axis of the transversely isotropic (TI) representation of the region can be rotated, leading to tilted transversely isotropic (TTI) media. We seek to find the optimal full-waveform inversion (FWI) strategy to estimate both the seismic velocities and the anisotropic parameters, including the tilt angle, in the presence of elastic TTI media. We first formulate the forward and inverse problems for elastic TTI media and analyse the radiation patterns of the model parameters. Based on the analyses of the radiation patterns, we propose two similar multistage FWI strategies that add inversion parameters over three stages, beginning with the isotropic parameters (horizontal P- and vertical S-wave velocity) and moving to the anisotropic parameters; the tilt angle is directly inverted in the last stage. Since diving waves, which are useful for providing long-wavelength updates, are mainly controlled by horizontal motion in anisotropic media, it is reasonable to choose the horizontal P-wave velocity rather than the vertical P-wave velocity. Then, the anisotropic parameters are inverted mainly using the reflected waves based on the isotropic background model built in the first stage. The main difference between the two multistage FWI strategies is whether the anisotropic parameter η is inverted. Comparing the two multistage FWI strategies with the simultaneous inversion strategy for a downsized version of the synthetic BP TTI model, we confirm that the multistage FWI strategies yield better inversion results than the simultaneous inversion strategy. When we compare the two multistage FWI strategies with each other for surface seismic data, ignoring η during the FWI process (focused multistage FWI) yields better inversion results for the tilt angle than those obtained with the inversion of η because η has less influence on the FWI than the other parameters and is not recovered well, which plays a role in degrading the tilt angle. Numerical examples support our conclusions that the focused multistage FWI strategy (neglecting η) is the optimal FWI strategy for TTI media and achieves computational efficiency for surface seismic data.
23
Borisov, Dmitry, Ryan Modrak, Fuchun Gao e Jeroen Tromp. "3D elastic full-waveform inversion of surface waves in the presence of irregular topography using an envelope-based misfit function". GEOPHYSICS 83, n. 1 (1 gennaio 2018): R1—R11. http://dx.doi.org/10.1190/geo2017-0081.1.
Abstract (sommario):
Full-waveform inversion (FWI) is a powerful method for estimating the earth’s material properties. We demonstrate that surface-wave-driven FWI is well-suited to recovering near-surface structures and effective at providing S-wave speed starting models for use in conventional body-wave FWI. Using a synthetic example based on the SEG Advanced Modeling phase II foothills model, we started with an envelope-based objective function to invert for shallow large-scale heterogeneities. Then we used a waveform-difference objective function to obtain a higher-resolution model. To accurately model surface waves in the presence of complex tomography, we used a spectral-element wave-propagation solver. Envelope misfit functions are found to be effective at minimizing cycle-skipping issues in surface-wave inversions, and surface waves themselves are found to be useful for constraining complex near-surface features.
24
Louboutin, Mathias, Philipp Witte, Michael Lange, Navjot Kukreja, Fabio Luporini, Gerard Gorman e Felix J. Herrmann. "Full-waveform inversion, Part 2: Adjoint modeling". Leading Edge 37, n. 1 (gennaio 2018): 69–72. http://dx.doi.org/10.1190/tle37010069.1.
Abstract (sommario):
This is the second part of a three-part tutorial series on full-waveform inversion (FWI) in which we provide a step-by-step walk through of setting up forward and adjoint wave equation solvers and an optimization framework for inversion. In Part 1 ( Louboutin et al., 2017 ), we showed how to use Devito ( http://www.opesci.org/devito-public ) to set up and solve acoustic wave equations with (impulsive) seismic sources and sample wavefields at the receiver locations to forward model shot records. Here in Part 2, we will discuss how to set up and solve adjoint wave equations with Devito and, from that, how we can calculate gradients and function values of the FWI objective function.
25
Zhang, Zhendong, Tariq Alkhalifah, Zedong Wu, Yike Liu, Bin He e Juwon Oh. "Normalized nonzero-lag crosscorrelation elastic full-waveform inversion". GEOPHYSICS 84, n. 1 (1 gennaio 2019): R1—R10. http://dx.doi.org/10.1190/geo2018-0082.1.
Abstract (sommario):
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.
26
Qi, Qiyuan, Wensha Huang, Donghao Zhang e Liguo Han. "Robust Elastic Full-Waveform Inversion Based on Normalized Cross-Correlation Source Wavelet Inversion". Applied Sciences 13, n. 24 (6 dicembre 2023): 13014. http://dx.doi.org/10.3390/app132413014.
Abstract (sommario):
The elastic full-waveform inversion (EFWI) method efficiently utilizes the amplitude, phase, and travel time information present in multi-component seismic recordings to create detailed parameter models of subsurface structures. Within full-waveform inversion (FWI), accurate source wavelet estimation significantly impacts both the convergence and final result quality. The source wavelet, serving as the initial condition for the wave equation’s forward modeling algorithm, directly influences the matching degree between observed and synthetic data. This study introduces a novel method for estimating the source wavelet utilizing cross-correlation norm elastic waveform inversion (CNEWI) and outlines the EFWI algorithm flow based on this CNEWI source wavelet inversion. The CNEWI method estimates the source wavelet by employing normalized cross-correlation processing on near-offset direct waves, thereby reducing the susceptibility to strong amplitude interference such as bad traces and surface wave residuals. The proposed CNEWI method exhibits a superior computational efficiency compared to conventional L2-norm waveform inversion for source wavelet estimation. Numerical experiments, including in ideal scenarios, with seismic data with bad traces, and with multi-component data, validate the advantages of the proposed method in both source wavelet estimation and EFWI compared to the traditional inversion method.
27
Li, Jing, Sherif Hanafy, Zhaolun Liu e Gerard T. Schuster. "Wave-equation dispersion inversion of Love waves". GEOPHYSICS 84, n. 5 (1 settembre 2019): R693—R705. http://dx.doi.org/10.1190/geo2018-0039.1.
Abstract (sommario):
We present a theory for wave-equation inversion of Love-wave dispersion curves, in which the misfit function is the sum of the squared differences between the wavenumbers along the predicted and observed dispersion curves. Similar to inversion of Rayleigh-wave dispersion curves, the complicated Love-wave arrivals in traces are skeletonized as simpler data, namely, the picked dispersion curves in the [Formula: see text] domain. Numerical solutions to the SH-wave equation and an iterative optimization method are then used to invert these dispersion curves for the S-wave velocity model. This procedure, denoted as wave-equation dispersion inversion of Love waves (LWD), does not require the assumption of a layered model or smooth velocity variations, and it is less prone to the cycle-skipping problems of full-waveform inversion. We demonstrate with synthetic and field data examples that LWD can accurately reconstruct the S-wave velocity distribution in a laterally heterogeneous medium. Compared with Rayleigh waves, inversion of the Love-wave dispersion curves empirically exhibits better convergence properties because they are completely insensitive to the P-velocity variations. In addition, Love-wave dispersion curves for our examples are simpler than those for Rayleigh waves, and they are easier to pick in our field data with a low signal-to-noise ratio.
28
Wang, Haiyang, Satish C. Singh e Henri Calandra. "Integrated inversion using combined wave-equation tomography and full waveform inversion". Geophysical Journal International 198, n. 1 (21 maggio 2014): 430–46. http://dx.doi.org/10.1093/gji/ggu138.
29
Xu, Tong, George A. McMechan e Robert Sun. "3-D prestack full‐wavefield inversion". GEOPHYSICS 60, n. 6 (novembre 1995): 1805–18. http://dx.doi.org/10.1190/1.1443913.
Abstract (sommario):
A full‐wavefield inversion algorithm for direct imaging of a 3-D compressional wave velocity distribution is based on the full 3-D scalar wave equation and operates on common‐source data recorded by areal arrays. For each source, the method involves reverse‐time extrapolation of the residual wavefield. Application of the image condition by crosscorrelation with the source wavefield at each time step produces a 3-D image whose amplitude at each point is proportional to the required velocity update at that point. Convergence to local minima is mitigated against by gradually increasing the wavenumber bandwidth in the estimated 3-D velocity distribution as iterations proceed, starting from the smallest wavenumber. The algorithm is illustrated by successful application to synthetic data for a multilayered monocline, and for a multilayered structure with the geometry of the standard French model. The latter demonstrates good performance with noisy, unequally spaced data with significant elevation statics.
30
Cheng, Guangsen, Xingyao Yin, Zhaoyun Zong, Tongxing Xia, Jianli Wang e Haojie Liu. "Seismic inversion using complex spherical-wave reflection coefficient at different offsets and frequencies". GEOPHYSICS 87, n. 2 (10 gennaio 2022): R183—R192. http://dx.doi.org/10.1190/geo2020-0787.1.
Abstract (sommario):
Compared with the plane-wave reflection coefficient, the spherical-wave reflection coefficient (SRC) can more accurately describe the reflected wavefield excited by a point source, especially in the case of low seismic frequency and short travel distance. However, unlike the widely used plane-wave amplitude-variation-with-offset/frequency (AVO/AVF) inversion, the practical application of spherical-wave AVO/AVF inversion in multilayer elastic media is still in the exploratory stage. One of the difficulties is how to fully use the amplitude and phase information of the complex-valued SRC and the spherical-wave response property of each frequency component to obtain the spherical-wave synthetic seismogram in multilayer elastic media. In view of this, we have developed a complex convolution model considering the amplitude and phase information of an SRC to obtain the complex synthetic seismogram of a certain frequency component. A simple harmonic superposition method is further developed. By superposing the complex synthetic seismograms of different frequency components, the synthetic seismogram of the full-frequency band can be obtained. In addition, a novel three-parameter SRC in terms of P- and S-wave moduli and density is derived. Based on the SRC and complex seismic traces with different offsets (or incidence angles) and frequency components, an inversion approach of complex spherical-wave amplitude and phase variation with offset and frequency is proposed. A noisy synthetic data example verifies the robustness of our complex spherical-wave inversion approach. Field data examples indicate that the P- and S-wave moduli estimated by the complex spherical-wave inversion approach can reasonably match the filtered well-logging data. Considering spherical waves rather than plane waves can improve the accuracy of seismic inversion results.
31
Borisov, Dmitry, Fuchun Gao, Paul Williamson e Jeroen Tromp. "Application of 2D full-waveform inversion on exploration land data". GEOPHYSICS 85, n. 2 (9 gennaio 2020): R75—R86. http://dx.doi.org/10.1190/geo2019-0082.1.
Abstract (sommario):
Estimating subsurface seismic properties is an important topic in civil engineering, oil and gas exploration, and global seismology. We have developed an application of 2D elastic waveform inversion with an active-source on-shore data set, as is typically acquired in exploration seismology on land. The maximum offset is limited to 12 km, and the lowest available frequency is 5 Hz. In such a context, surface waves are generally treated as noise and are removed as a part of data processing. In contrast to the conventional approach, our workflow starts by inverting surface waves to constrain shallow parts of the shear wavespeed model. To mitigate cycle skipping, frequency- and offset-continuation approaches are used. To accurately take into account free-surface effects (and irregular topography), a spectral-element-based wave propagation solver is used for forward modeling. To reduce amplitude influences, a normalized crosscorrelation (NC) objective function is used in conjunction with systematic updates of the source wavelet during the inversion process. As the inversion proceeds, body waves are gradually incorporated in the process. At the final stage, surface and body waves are inverted together using the entire offset range over the band between 5 and 15 Hz. The inverted models include high-resolution features in the first 500 m of compressional and shear wavespeeds, with some model updates down to 4.0 km in the first parameter. The inversion results confirmed by well-log information, indicate a better fit of compressional to shear wavespeeds ratios compared with the initial model. The final data fit is also noticeably improved compared to the initial one. Although our results confirm previous studies demonstrating that an NC norm combined with a source time function correction can partly stabilize purely elastic inversions of viscoelastic data, we believe that including an attenuation depth model in the forward simulation gives better results.
32
Groos, Lisa, Martin Schäfer, Thomas Forbriger e Thomas Bohlen. "Application of a complete workflow for 2D elastic full-waveform inversion to recorded shallow-seismic Rayleigh waves". GEOPHYSICS 82, n. 2 (1 marzo 2017): R109—R117. http://dx.doi.org/10.1190/geo2016-0284.1.
Abstract (sommario):
The S-wave velocity of the shallow subsurface can be inferred from shallow-seismic Rayleigh waves. Traditionally, the dispersion curves of the Rayleigh waves are inverted to obtain the (local) S-wave velocity as a function of depth. Two-dimensional elastic full-waveform inversion (FWI) has the potential to also infer lateral variations. We have developed a novel workflow for the application of 2D elastic FWI to recorded surface waves. During the preprocessing, we apply a line-source simulation (spreading correction) and perform an a priori estimation of the attenuation of waves. The iterative multiscale 2D elastic FWI workflow consists of the preconditioning of the gradients in the vicinity of the sources and a source-wavelet correction. The misfit is defined by the least-squares norm of normalized wavefields. We apply our workflow to a field data set that has been acquired on a predominantly depth-dependent velocity structure, and we compare the reconstructed S-wave velocity model with the result obtained by a 1D inversion based on wavefield spectra (Fourier-Bessel expansion coefficients). The 2D S-wave velocity model obtained by FWI shows an overall depth dependency that agrees well with the 1D inversion result. Both models can explain the main characteristics of the recorded seismograms. The small lateral variations in S-wave velocity introduced by FWI additionally explain the lateral changes of the recorded Rayleigh waves. The comparison thus verifies the applicability of our 2D FWI workflow and confirms the potential of FWI to reconstruct shallow small-scale lateral changes of S-wave velocity.
33
Rao, Jing, Madis Ratassepp e Zheng Fan. "Guided Wave Tomography Based on Full Waveform Inversion". IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control 63, n. 5 (maggio 2016): 737–45. http://dx.doi.org/10.1109/tuffc.2016.2536144.
34
Minkoff, Susan E., e William W. Symes. "Full waveform inversion of marine reflection data in the plane‐wave domain". GEOPHYSICS 62, n. 2 (marzo 1997): 540–53. http://dx.doi.org/10.1190/1.1444164.
Abstract (sommario):
Full waveform inversion of a p‐τ marine data set from the Gulf of Mexico provides estimates of the long‐wavelength P‐wave background velocity, anisotropic seismic source, and three high‐frequency elastic parameter reflectivities that explain 70% of the total seismic data and 90% of the data in an interval around the gas sand target. The forward simulator is based on a plane‐wave viscoelastic model for P‐wave propagation and primary reflections in a layered earth. Differential semblance optimization, a variant of output least‐squares inversion, successfully estimates the nonlinear P‐wave background velocity and linear reflectivities. Once an accurate velocity is estimated, output least‐squares inversion reestimates the reflectivities and an anisotropic seismic source simultaneously. The viscoelastic model predicts the amplitude‐versus‐angle trend in the data more accurately than does an elastic model. Simultaneous inversion for reflectivities and source explains substantially more of the actual data than does inversion for reflectivities with fixed source from an air‐gun modeler. The best reflectivity estimates conform to widely accepted lithologic relationships and closely match the filtered well logs.
35
Agudo, Òscar Calderón, Nuno Vieira da Silva, Michael Warner e Joanna Morgan. "Acoustic full-waveform inversion in an elastic world". GEOPHYSICS 83, n. 3 (1 maggio 2018): R257—R271. http://dx.doi.org/10.1190/geo2017-0063.1.
Abstract (sommario):
Full-waveform inversion (FWI) is a technique used to obtain high-quality velocity models of the subsurface. Despite the elastic nature of the earth, the anisotropic acoustic wave equation is typically used to model wave propagation in FWI. In part, this simplification is essential for being efficient when inverting large 3D data sets, but it has the adverse effect of reducing the accuracy and resolution of the recovered P-wave velocity models, as well as a loss in potential to constrain other physical properties, such as the S-wave velocity given that amplitude information in the observed data set is not fully used. Here, we first apply conventional acoustic FWI to acoustic and elastic data generated using the same velocity model to investigate the effect of neglecting the elastic component in field data and we find that it leads to a loss in resolution and accuracy in the recovered velocity model. Then, we develop a method to mitigate elastic effects in acoustic FWI using matching filters that transform elastic data into acoustic data and find that it is applicable to marine and land data sets. Tests show that our approach is successful: The imprint of elastic effects on the recovered P-wave models is mitigated, leading to better-resolved models than those obtained after conventional acoustic FWI. Our method requires a guess of [Formula: see text] and is marginally more computationally demanding than acoustic FWI, but much less so than elastic FWI.
36
Groos, Lisa, Martin Schäfer, Thomas Forbriger e Thomas Bohlen. "The role of attenuation in 2D full-waveform inversion of shallow-seismic body and Rayleigh waves". GEOPHYSICS 79, n. 6 (1 novembre 2014): R247—R261. http://dx.doi.org/10.1190/geo2013-0462.1.
Abstract (sommario):
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].
37
Nguyen, Trung Dung, e Khiem T. Tran. "Site characterization with 3D elastic full-waveform tomography". GEOPHYSICS 83, n. 5 (1 settembre 2018): R389—R400. http://dx.doi.org/10.1190/geo2017-0571.1.
Abstract (sommario):
We have developed a 3D elastic full-waveform inversion (FWI) method for geotechnical site characterization. The method is based on a solution of 3D elastic-wave equations for forward modeling to simulate wave propagation and a local optimization approach based on the adjoint-state method to update the model parameters. The staggered-grid finite-difference technique is used to solve the wave equations together with implementation of the perfectly matched layer condition for boundary truncation. Seismic wavefields are acquired from geophysical testing using sensors and sources located in uniform 2D grids on the ground surface, and they are then inverted for the extraction of 3D subsurface wave velocity structures. The capability of the presented FWI method is tested on synthetic and field data sets. The inversion results from synthetic data indicate the ability of characterizing laterally variable low- and high-velocity layers. Field experimental data were collected using 96 receivers and a propelled energy generator to induce seismic wave energy. The field data result indicates that the waveform analysis was able to delineate variable subsurface soil layers. The seismic inversion results are generally consistent with invasive standard penetration test [Formula: see text]-values, including identification of a low-velocity zone.
38
Um, Evan Schankee, Michael Commer e Gregory A. Newman. "A strategy for coupled 3D imaging of large-scale seismic and electromagnetic data sets: Application to subsalt imaging". GEOPHYSICS 79, n. 3 (1 maggio 2014): ID1—ID13. http://dx.doi.org/10.1190/geo2013-0053.1.
Abstract (sommario):
Offshore seismic and electromagnetic (EM) imaging for hydrocarbons can require up to tens of millions of parameters to describe the 3D distribution of complex seabed geology and relevant geophysical attributes. The imaging and data volumes for such problems are enormous. Descent-based methods are the only viable imaging approach, where it is often challenging to manage the convergence of stand-alone seismic and EM inversion experiments. When a joint seismic-EM inversion is implemented, convergence problems with descent-based methods are further aggravated. Moreover, resolution mismatches between seismic and EM pose another challenge for joint inversion. To overcome these problems, we evaluated a coupled seismic-EM inversion workflow and applied it to a set of full-wave-seismic, magnetotelluric (MT) and controlled-source electromagnetic (CSEM) data for subsalt imaging. In our workflow, we address disparate resolution properties between seismic and EM data by implementing the seismic inversion in the Laplace domain, where the wave equation is transformed into a diffusion equation. The resolution of seismic data thus becomes comparable to that of EM data. To mitigate the convergence problems, the full joint seismic-EM inverse problem is split into manageable components: separate seismic and EM inversions and an intermediate step that enforces structural coupling through a cross-gradient-only inversion and resistivity-velocity crossplots. In this workflow, stand-alone seismic and MT inversion are performed first. The cross-gradient-only inversion and the crossplots are used to precondition the resistivity and velocity models for subsequent stand-alone inversions. By repeating the sequence of the stand-alone seismic, MT, and cross-gradient-only inversions along with the crossplots, we introduce the seismic structural information into the resistivity model, and vice versa, significantly improving the salt geometry in both resistivity and velocity images. We conclude that the improved salt geometry can then be used to precondition a starting model for CSEM inversions, yielding significant improvement in the resistivity images of hydrocarbon reservoirs adjacent to the salt.
39
Athanasopoulos, Nikolaos, Edgar Manukyan, Thomas Bohlen e Hansruedi Maurer. "Time–frequency windowing in multiparameter elastic FWI of shallow seismic wavefield". Geophysical Journal International 222, n. 2 (15 maggio 2020): 1164–77. http://dx.doi.org/10.1093/gji/ggaa242.
Abstract (sommario):
SUMMARY Full-waveform inversion of shallow seismic wavefields is a promising method to infer multiparameter models of elastic material properties (S-wave velocity, P-wave velocity and mass density) of the shallow subsurface with high resolution. Previous studies used either the refracted Pwaves to reconstructed models of P-wave velocity or the high-amplitude Rayleigh waves to infer the S-wave velocity structure. In this work, we propose a combination of both wavefields using continuous time–frequency windowing. We start with the contribution of refracted P waves and gradually increase the time window to account for scattered body waves, higher mode Rayleigh waves and finally the fundamental Rayleigh wave mode. The opening of the time window is combined with opening the frequency bandwidth of input signals to avoid cycle skipping. Synthetic reconstruction tests revealed that the reconstruction of P-wave velocity model and mass density can be improved. The S-wave velocity reconstruction is still accurate and robust and is slightly benefitted by time–frequency windowing. In a field data application, we observed that time–frequency windowing improves the consistency of multiparameter models. The inferred models are in good agreement with independent geophysical information obtained from ground-penetrating radar and full-waveform inversion of SH waves.
40
Wang, Kai, Meiyan Guo, Qingxia Xiao, Chuanyi Ma, Lingli Zhang, Xinji Xu, Ming Li e Ningbo Li. "Frequency Domain Full Waveform Inversion Method of Acquiring Rock Wave Velocity in Front of Tunnels". Applied Sciences 11, n. 14 (8 luglio 2021): 6330. http://dx.doi.org/10.3390/app11146330.
Abstract (sommario):
Ahead geological prospecting, which can estimate adverse geology ahead of the tunnel face, is necessary in the process of tunnel construction. Due to its long detection range and good recognition effect on the interface, the seismic method is widely used in tunnel ahead prospecting. However, the observation space in tunnels is quite narrow compared to ground seismic prospecting, which leads to some problems in the acquisition of wave velocity, including: the velocity of the direct wave is used to replace the wave velocity of the forward rock approximately; the arrival time information of seismic waves is the main factor in time-travel inversion or the tomography method, which is sufficient to provide a simple model rather than deal with complex geological conditions. In view of the above problems, the frequency domain full waveform inversion method in ground prospecting is introduced to tunnel seismic prospecting. In addition, the optimized difference format is given according to the particularity of the tunnel environment. In this method, the kinematics and dynamics of the seismic wavefield are fully used to obtain more accurate wave velocity results. Simultaneously, forward modeling and inversion simulations on tunnel samples with typical adverse geological bodies are given here, which verified the validity and reliability of the proposed method.
41
Zhang, Wensheng. "Frequency-domain elastic full-waveform inversion based on Green functions". Journal of Physics: Conference Series 2444, n. 1 (1 febbraio 2023): 012003. http://dx.doi.org/10.1088/1742-6596/2444/1/012003.
Abstract (sommario):
Abstract In this paper, the elastic full-waveform inversion in the frequency domain based on Green functions is investigated. The method allows to image media velocities of compressional wave and shear wave by using vertical and horizontal components. The forward problem is solved by the finite difference scheme with the perfectly matched layer in the frequency domain. The inversion is an optimization iterative process to minimize the residual between the synthetic data and the observed data on the surface. A gradient method is used to find the optimization direction, preconditioned with the diagonal part of the inverse Hessian matrix. The gradient is constructed based on Green functions, which serves to compute the gradient of the objective function with respect to model parameters effectively. The inversion is completed from low frequency to high frequency so that smaller wavelengths are progressively introduced in inversion. Numerical computations for the benchmark Overthrust model demonstrate the full-waveform inversion method in this paper can inverse media velocities of complex structures accurately.
42
Feng, Zongcai, e Lianjie Huang. "Shear reflectivity compensation in full-waveform inversion using least-squares reverse-time migration". Geophysical Journal International 227, n. 1 (19 maggio 2021): 1–15. http://dx.doi.org/10.1093/gji/ggab193.
Abstract (sommario):
SUMMARY The computational cost of elastic-waveform inversion is too high for inverting PP reflections, while using acoustic full-waveform inversion (FWI) is inaccurate because it does not depend on the shear modulus/velocity/impedance that affects elastic PP wavefield amplitudes. To solve this problem, we develop a waveform inversion method that uses acoustic least-squares reverse-time migration (LSRTM) to compensate the shear reflectivity for acoustic FWI. Our method is based on the quasi-elastic-wave equation developed by Chapman et al. (2014). The quasi-elastic-wave equation uses a linearized acoustic-wave equation with shear modulus μ as a virtual source to correct the acoustic PP wavefield amplitudes toward elastic ones. Our waveform inversion method inverts for elastic parameters by minimizing the L2 norm of the difference between recorded and predicted PP reflections modelled using the quasi-elastic-wave equation. Numerical tests on synthetic and field data show that our method can properly handle the amplitudes of elastic PP reflections and provides an accurate estimate of the P- and S-wave velocities/impedances and, in some cases, the density. The method does not need the computationally expensive numerical solution to the elastic-wave equation. It also gives a better estimate of elastic parameters than a pure LSRTM method for elastic PP reflections.
43
Xing, Zhen, e Alfredo Mazzotti. "Two-grid full-waveform Rayleigh-wave inversion via a genetic algorithm — Part 2: Application to two actual data sets". GEOPHYSICS 84, n. 5 (1 settembre 2019): R815—R825. http://dx.doi.org/10.1190/geo2018-0800.1.
Abstract (sommario):
We have applied our two-grid genetic-algorithm Rayleigh-wave full-waveform inversion (FWI) to two actual data sets acquired in Luni (Italy) and Grenoble (France), respectively. Because our technique used 2D elastic finite-difference modeling for solving the forward problem, the observed data were 3D to 2D corrected prior to the inversion. To limit the computing time, both inversions focused on predicting low-resolution, smooth models by using quite coarse inversion grids. The wavelets for FWI were estimated directly from the observed data by using the Wiener method. In the Luni case, due to the strong dispersion effects on the data, to strengthen the inversion, envelopes and waveforms were considered in the objective function and an offset-marching strategy was applied. Though no a priori information was exploited, the outcomes of the Luni and Grenoble data inversion were fair. The predicted Luni [Formula: see text] model indicates a strong velocity increase from approximately 3 to 6 m, and velocity inversions have been detected at approximately 2 and 9 m depths. Analyzing the dispersion spectra, it results that the predicted Luni data reasonably reproduced the waveforms related to the fundamental mode and, likely, a small part of those related to the first higher mode. Concerning the Grenoble example, the predicted [Formula: see text] model coincides reasonably well with the long-wavelength structures presented in the [Formula: see text] profiles obtained from nearby boreholes. The data reconstruction is generally satisfactory, and when mismatches occur between the predicted and observed traces, the phase differences are always within half-periods. The fair inversion outcomes suggest that the predicted Luni and Grenoble models would likely be adequate initial models for local FWI, which could further increase the resolution and the details of the estimated [Formula: see text] models.
44
Zhang, Zhen-Dong, e Tariq Alkhalifah. "Wave-equation Rayleigh-wave dispersion inversion using fundamental and higher modes". GEOPHYSICS 84, n. 4 (1 luglio 2019): EN57—EN65. http://dx.doi.org/10.1190/geo2018-0506.1.
Abstract (sommario):
Recorded surface waves often provide reasonable estimates of the S-wave velocity in the near surface. However, existing algorithms are mainly based on the 1D layered-model assumption and require picking the dispersion curves either automatically or manually. We have developed a wave-equation-based inversion algorithm that inverts for S-wave velocities using fundamental and higher mode Rayleigh waves without picking an explicit dispersion curve. Our method aims to maximize the similarity of the phase velocity spectrum ([Formula: see text]) of the observed and predicted surface waves with all Rayleigh-wave modes (if they exist) included in the inversion. The [Formula: see text] spectrum is calculated using the linear Radon transform applied to a local similarity-based objective function; thus, we do not need to pick velocities in spectrum plots. As a result, the best match between the predicted and observed [Formula: see text] spectrum provides the optimal estimation of the S-wave velocity. We derive S-wave velocity updates using the adjoint-state method and solve the optimization problem using a limited-memory Broyden-Fletcher-Goldfarb-Shanno algorithm. Our method excels in cases in which the S-wave velocity has vertical reversals and lateral variations because we used all-modes dispersion, and it can suppress the local minimum problem often associated with full-waveform inversion applications. Synthetic and field examples are used to verify the effectiveness of our method.
45
Liu, Zhaolun, Jing Li, Sherif M. Hanafy e Gerard Schuster. "3D wave-equation dispersion inversion of Rayleigh waves". GEOPHYSICS 84, n. 5 (1 settembre 2019): R673—R691. http://dx.doi.org/10.1190/geo2018-0543.1.
Abstract (sommario):
The 2D wave-equation dispersion (WD) inversion method is extended to 3D wave-equation dispersion inversion of surface waves for the shear-velocity distribution. The objective function of 3D WD is the frequency summation of the squared wavenumber [Formula: see text] differences along each azimuth angle of the fundamental or higher modes of Rayleigh waves in each shot gather. The S-wave velocity model is updated by the weighted zero-lag crosscorrelation between the weighted source-side wavefield and the back-projected receiver-side wavefield for each azimuth angle. A multiscale 3D WD strategy is provided, which starts from the pseudo-1D S-velocity model, which is then used to get the 2D WD tomogram, which in turn is used as the starting model for 3D WD. The synthetic and field data examples demonstrate that 3D WD can accurately reconstruct the 3D S-wave velocity model of a laterally heterogeneous medium and has much less of a tendency to getting stuck in a local minimum compared with full-waveform inversion.
46
Tao, Yi, e Mrinal K. Sen. "Frequency-domain full waveform inversion with plane-wave data". GEOPHYSICS 78, n. 1 (1 gennaio 2013): R13—R23. http://dx.doi.org/10.1190/geo2012-0267.1.
Abstract (sommario):
We derived an efficient frequency-domain full waveform inversion (FWI) method using plane-wave encoded shot records. The forward modeling involved application of position dependent linear time shifts at all source locations. This was followed by propagation of wavefields into the medium from all shotpoints simultaneously. The gradient of the cost function needed in the FWI was calculated first by transforming the densely sampled seismic data into the frequency-ray parameter domain and then backpropagating the residual wavefield using an adjoint-state approach. We used a Gauss-Newton framework for model updating. The approximate Hessian matrix was formed with a plane-wave encoding strategy, which required a summation over source and receiver ray parameters of the Green’s functions. Plane-wave encoding considerably reduces the computational burden and crosstalk artifacts are effectively suppressed by stacking over different ray parameters. It also has the advantage of directional illumination of the selected targets. Numerical examples show the accuracy and efficiency of our method.
47
Tang, Zhiyuan. "Frequency Domain Wave Equation Inversion and Its Application on the Heterogeneous Reservoir Model Data". Earth Science Research 6, n. 1 (12 dicembre 2016): 55. http://dx.doi.org/10.5539/esr.v6n1p55.
Abstract (sommario):
Seismic full waveform inversion seeks to make use of the full information based on full wave field modeling to extract quantitative information from seismograms. Its serious nonlinearity and high dependence on initial velocity model often results in unsatisfactory inversion results in paleo-karsts carbonate reservoir characterized by strong heterogeneity. The paper presents an improved strategy of multi-scale inversion to establish velocity field model of waveform tomography. the forward wave equation algorithm was derived in frequency domain, and then the Matrix formalism for the iterative inverse methods is derived by gradient methods to speed up calculation and to avoid convergence to local minimum value. After massive amount of frequencies tests, the appropriate bandwidth are extracted, and the velocity field calculated at low frequency is used as the input of the high frequency. After the iteration, the accurate velocity field is inverted. Finally, frequency domain wave equation full waveform inversion in mathematical and physical models is conducted in order to verify the inverse program. The method of selecting the inverse frequencies is proved to be effective.
48
Jaimes-Osorio, Ligia Elena, Alison Malcolm e Ali Gholami. "Amplitude variation with offset inversion using acoustic-elastic local solver". GEOPHYSICS 85, n. 3 (1 aprile 2020): R251—R262. http://dx.doi.org/10.1190/geo2019-0108.1.
Abstract (sommario):
Conventional amplitude variation with offset (AVO) inversion analysis uses the Zoeppritz equations, which are based on a plane-wave approximation. However, because real seismic data are created by point sources, wave reflections are better modeled by spherical waves than by plane waves. Indeed, spherical reflection coefficients deviate from planar reflection coefficients near the critical and postcritical angles, which implies that the Zoeppritz equations are not applicable for angles close to critical reflection in AVO analysis. Elastic finite-difference simulations provide a solution to the limitations of the Zoeppritz approximation because they can handle near- and postcritical reflections. We have used a coupled acoustic-elastic local solver that approximates the wavefield with high accuracy within a locally perturbed elastic subdomain of the acoustic full domain. Using this acoustic-elastic local solver, the local wavefield generation and inversion are much faster than performing a full-domain elastic inversion. We use this technique to model wavefields and to demonstrate that the amplitude from within the local domain can be used as a constraint in the inversion to recover elastic material properties. Then, we focus on understanding how much the amplitude and phase contribute to the reconstruction accuracy of the elastic material parameters ([Formula: see text], [Formula: see text], and [Formula: see text]). Our results suggest that the combination of amplitude and phase in the inversion helps with the convergence. Finally, we analyze elastic parameter trade-offs in AVO inversion, from which we find that to recover accurate P-wave velocities we should invert for [Formula: see text] and [Formula: see text] simultaneously with fixed density.
49
Geary, Andrew. "Seismic Soundoff: A reality check on full-wave inversion". Leading Edge 42, n. 8 (agosto 2023): 580. http://dx.doi.org/10.1190/tle42080580.1.
50
Zhang, Tuo, e Christoph Sens-Schönfelder. "Adjoint envelope tomography for scattering and absorption using radiative transfer theory". Geophysical Journal International 229, n. 1 (11 novembre 2021): 566–88. http://dx.doi.org/10.1093/gji/ggab457.
Abstract (sommario):
SUMMARY To investigate the small-scale elastic structure of the subsurface at length scales below the resolution limits of waveform tomography, envelopes of high-frequency scattered seismic waveforms have been used with a variety of approaches. However, a rigorous framework for the iterative inversion of seismogram envelopes to image heterogeneity and high-frequency attenuation comparable to full waveform inversion (FWI) is missing. We present the mathematical framework for an iterative full envelope inversion using forward and adjoint simulations of the radiative transfer equations, in full analogy to FWI that is based on the wave equation. The forward and adjoint problems are solved by modelling 2-D multiple non-isotropic scattering in a random elastic medium with spatially variable heterogeneity and attenuation using the Monte Carlo method. Sensitivity kernels are derived for the squared difference between the full observed and modelled envelopes which is iteratively minimized with the L-BFGS method. We apply this algorithm in numerical tests in the acoustic approximation and show that it is possible to image the spatial distribution of small-scale heterogeneity and attenuation in iterative inversions. Our analysis shows that the relative importance of scattering and attenuation anomalies needs to be considered when the model resolution is assessed. The inversions confirm that the early coda is important for imaging the distribution of heterogeneity while later coda waves are more sensitive to intrinsic attenuation and we show that this dependency can be used to cope with the trade-off that exists between both material properties.