Generalized linear inversion method for near‐surface modeling: A 3‐D processing case history

1992 ◽  
Author(s):  
Wendy R. Stresau ◽  
Lee Hansen ◽  
R. Wilson Rowe ◽  
Miguel Morales
Keyword(s):  
Case History ◽  
Surface Modeling ◽  
Inversion Method ◽  
Linear Inversion ◽  
Near Surface ◽  
Generalized Linear Inversion

A reality check on full-wave inversion applied to land seismic data for near-surface modeling

2022 ◽  
Vol 41 (1) ◽  
pp. 40-46
Author(s):  
Öz Yilmaz ◽  
Kai Gao ◽  
Milos Delic ◽  
Jianghai Xia ◽  
Lianjie Huang ◽  
...  
Keyword(s):  
Wave Velocity ◽  
Seismic Data ◽  
Surface Modeling ◽  
P Wave ◽  
Borehole Data ◽  
Near Surface ◽  
Traveltime Tomography ◽  
Wave Inversion ◽  
Full Wave ◽  
Elastic Inversion

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.

FULL-WAVEFORM INVERSION WITH SOURCE AND RECEIVER COUPLING EFFECTS CORRECTION

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

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

Near-surface modeling: A multiphysics approach

2016 ◽  
Vol 35 (11) ◽  
pp. 968-976 ◽  
Author(s):  
Marco Mantovani
Keyword(s):  
Surface Modeling ◽  
Near Surface

2-D and 3-D resistivity image reconstruction using crosshole data

Geophysics ◽  
1992 ◽  
Vol 57 (10) ◽  
pp. 1270-1281 ◽  
Author(s):  
Hiromasa Shima
Keyword(s):  
Field Test ◽  
Numerical Models ◽  
Electrical Potential ◽  
Nonlinear Least Squares ◽  
Three Dimensions ◽  
Inversion Method ◽  
Good Resolution ◽  
Inversion Algorithm ◽  
Linear Inversion ◽  
Resistivity Model

Theoretical changes in the distribution of electrical potential near subsurface resistivity anomalies have been studied using two resistivity models. The results suggest that the greatest response from such anomalies can be observed with buried electrodes, and that the resistivity model of a volume between boreholes can be accurately reconstructed by using crosshole data. The distributive properties of crosshole electrical potential data obtained by the pole‐pole array method have also been examined using the calculated partial derivative of the observed apparent resistivity with respect to a small cell within a given volume. The results show that for optimum two‐dimensional (2-D) and three‐dimensional (3-D) target imaging, in‐line data and crossline data should be combined, and an area outside the zone of exploration should be included in the analysis. In this paper, the 2-D and 3-D resistivity images presented are reconstructed from crosshole data by the combination of two inversion algorithms. The first algorithm uses the alpha center method for forward modeling and reconstructs a resistivity model by a nonlinear least‐squares inversion. Alpha centers express a continuously varying resistivity model, and the distribution of the electrical potential from the model can be calculated quickly. An initial general model is determined by the resistivity backprojection technique (RBPT) prior to the first inversion step. The second process uses finite elements and a linear inversion algorithm to improve the resolution of the resistivity model created by the first step. Simple 2-D and 3-D numerical models are discussed to illustrate the inversion method used in processing. Data from several field studies are also presented to demonstrate the capabilities of using crosshole resistivity exploration techniques. The numerical experiments show that by using the combined reconstruction algorithm, thin conductive layers can be imaged with good resolution for 2-D and 3-D cases. The integration of finite‐element computations is shown to improve the image obtained by the alpha center inversion process for 3-D applications. The first field test uses horizontal galleries to evaluate complex 2-D features of a zinc mine. The second field test illustrates the use of three boreholes at a dam site to investigate base rock features and define the distribution of an altered zone in three dimensions.

scholarly journals Estimation of the Rupture History of the 2008 Mw8.0 Wenchuan Earthquake by Joint Inversion of Teleseismic and Strong-Motion Records

2020 ◽  
Vol 2020 ◽  
pp. 1-13
Author(s):  
Deyu Yin ◽  
Yun Dong ◽  
Qifang Liu ◽  
Yuexin She ◽  
Jingke Wu ◽  
...  
Keyword(s):  
Wenchuan Earthquake ◽  
Strong Motion ◽  
Inversion Method ◽  
Near Surface ◽  
Strong Motion Records ◽  
The Wenchuan Earthquake ◽  
Finite Fault ◽  
History Of ◽  
Beichuan Fault ◽  
Dip Angle

In order to reproduce the rupture history of the 2008 Mw8.0 Wenchuan earthquake, the teleseismic and strong-motion records are adopted. Based on a multiple-segment, variable-slip model, the finite fault inversion method is utilized to recover the rupture process. The results are as follows: (1) the rupture duration of the Wenchuan earthquake is about 100 s, and the released seismic moment is 1.24 × 1021 N·m, equal to the moment magnitude Mw8.0. There are 5 asperities on the fault plane, indicating that the earthquake is composed of at least 5 subevents. (2) The slip is mainly distributed on the Beichuan fault, indicating that the Beichuan fault is the main rupture fault. On the southern part of the Beichuan fault, the dislocation underside the Longmenshan area and Hongkou-Yingxiu near-surface area is dominated by thrust, and the maximum slip is 11.8 m. Slip between the Yuejiashan and Qingping area is dominated by thrust. On the northern part of the Beichuan fault, the area under Beichuan is dominated by thrust, the slip under Nanba is thrust and strike, near Qingchuan, the slip turns into the strike slip, and the maximum slip is 13.1 m. The dislocation under Bailu is also dominated by thrust, with maximum slip 8.9 m. (3) The rupture of the Wenchuan earthquake is mainly a unilateral rupture to the northeast. The rupture started at the low dip angle part of the Beichuan fault, and after 3 s, it propagated to the Pengguan fault. After 10 s, the largest asperity under Longmenshan in the south section of the Beichuan fault began to break, lasting for about 24 s. Then, the Xiaoyudong fault was triggered by the Pengguan fault, and the bilateral rupture of the high dip angle part of the Beichuan fault started at about 6 s. South section of the Beichuan fault began to break at about 35 s, and at 43 s, 63 s, and 80 s, the rupture extended to Beichuan, Nanba, and Qingchuan areas.

Turning‐ray tomography

Geophysics ◽  
1995 ◽  
Vol 60 (6) ◽  
pp. 1917-1929 ◽  
Author(s):  
Joseph P. Stefani
Keyword(s):  
Velocity Structure ◽  
Initial Point ◽  
Real Data ◽  
Point Sources ◽  
Surface Velocity ◽  
Inversion Algorithm ◽  
Low Velocity ◽  
Linear Inversion ◽  
Near Surface ◽  
Velocity Inversion

Turning‐ray tomography is useful for estimating near‐surface velocity structure in areas where conventional refraction statics techniques fail because of poor data or lack of smooth refractor/velocity structure. This paper explores the accuracy and inherent smoothing of turning‐ray tomography in its capacity to estimate absolute near‐surface velocity and the statics times derived from these velocities, and the fidelity with which wavefields collapse to point diffractors when migrated through these estimated velocities. The method comprises nonlinear iterations of forward ray tracing through triangular cells linear in slowness squared, coupled with the LSQR linear inversion algorithm. It is applied to two synthetic finite‐ difference data sets of types that usually foil conventional refraction statics techniques. These models represent a complex hard‐rock overthrust structure with a low‐velocity zone and pinchouts, and a contemporaneous near‐shore marine trench filled with low‐ velocity unconsolidated deposits exhibiting no seismically apparent internal structure. In both cases velocities are estimated accurately to a depth of one‐ fifth the maximum offset, as are the associated statics times. Of equal importance, the velocities are sufficiently accurate to correctly focus synthetic wavefields back to their initial point sources, so migration/datuming applications can also use these velocities. The method is applied to a real data example from the Timbalier Trench in the Gulf of Mexico, which exhibits the same essential features as the marine trench synthetic model. The Timbalier velocity inversion is geologically reasonable and yields long and short wavelength statics that improve the CMP gathers and stack and that correctly align reflections to known well markers. Turning‐ray tomography estimates near‐surface velocities accurately enough for the three purposes of lithology interpretation, statics calculations, and wavefield focusing for shallow migration and datuming.

Sign in / Sign up

Export Citation Format

Share Document