Hierarchical approach of seismic full waveform inversion

被引：1

作者：

Asnaashari A. ^{[1
]}

Brossier R. ^{[1
]}

Castellanos C. ^{[2
]}

Dupuy B. ^{[1
]}

Etienne V. ^{[2
]}

Gholami Y. ^{[2
]}

Hu G. ^{[1
,2
]}

Métivier L. ^{[1
,2
]}

Operto S. ^{[2
]}

Pageot D. ^{[2
]}

Prieux V. ^{[2
]}

Ribodetti A. ^{[2
]}

Roques A. ^{[1
]}

Virieux J. ^{[1
]}

机构：

[1] ISTerre, Universit é de Grenoble I-CNRS, Universite Joseph Fourier-Grenoble I, Member of Institut Universitaire de France, IUF, Laboratory in Earth Sciences: ISTerre, Grenoble Cedex 9,38041

[2] Géoazur-Université Nice Sophia-Antipolis-CNRS, Nice

来源：

Asnaashari, A. (Amir.Asnaashari@obs.ujf-grenoble.fr) | 1600年 / Maik Nauka Publishing / Springer SBM卷 / 05期

关键词：

cycle skipping; optimization problem; quasi-Newton method; seismic traces;

D O I：

10.1134/S1995423912020012

中图分类号：

学科分类号：

摘要：

Full waveform inversion (FWI) of seismic traces recorded at the free surface allows the reconstruction of the physical parameters structure on the underlying medium. For such a reconstruction, an optimization problem is defined, where synthetic traces, obtained through numerical techniques as finite-difference or finite-element methods in a given model of the subsurface, should match the observed traces. The number of data samples is routinely around 1 billion for 2D problems and 1 trillion for 3D problems while the number of parameters ranges from 1 million to 10 million degrees of freedom. Moreover, if one defines the mismatch as the standard least-squares norm between values sampled in time/frequency and space, the misfit function has a significant number of secondary minima related to the ill-posedness and the nonlinearity of the inversion problem linked to the so-called cycle skipping. Taking into account the size of the problem, we consider a local linearized method where gradient is computed using the adjoint formulation of the seismic wave propagation problem. Starting for an initial model, we consider a quasi-Newtonian method, which allows us to formulate the reconstruction of various parameters such as P and S waves velocities or density or attenuation factors. A hierarchical strategy based on the incremental increase of the data complexity starting from low-frequency content to high-frequency content, from initial wavelets to later phases in the data space from narrow azimuths to wide azimuths and from simple observables to more complex ones. Different synthetic examples on realistic structures illustrate the efficiency of this strategy based on the data manipulation. This strategy related to the data space has to be inserted into a more global framework where we could improve significantly the probability to converge to the global minimum. When considering the model space, we may rely on the construction of the initial model or add constraints such as smoothness of the searched model and/or prior informations collected by other means. An alternative strategy concerns the building of the objective function and various possibilities must be considered, which may increase the linearity of the inversion procedure. © 2012 Pleiades Publishing, Ltd.

引用

页码：99 / 108

页数：9

共 30 条

[1] Asnaashari A., Brossier R., Garambois S., Audebert F., Thore P., Virieux J., Sensitivity Analysis of Time-Lapse Images Obtained by Differential Waveform Inversion with Respect to Reference Model, Expanded Abstr., Soc. Expl. Geophys., pp. 2482-2486, (2011)
[2] Beydoun W.B., Tarantola A., First Born and Rytov Approximation: Modeling and Inversion Conditions in a Canonical Example, J. Acoust. Soc. Am., 83, pp. 1045-1055, (1988)
[3] Billette F., Lambare G., Velocity Macro-Model Estimation from Seismic Reection Data by Stereotomography, Geophys. J. Int., 135, 2, pp. 671-680, (1998)
[4] Brenders A.J., Pratt R.G., Efficient Waveform Tomography for Lithospheric Imaging: Implications for Realistic 2D Acquisition Geometries and Low Frequency Data, Geophys. J. Int., 168, 1, pp. 152-170, (2007)
[5] Brossier R., Operto S., Virieux J., Seismic Imaging of Complex Onshore Structures by 2D Elastic Frequency-Domain Full-Waveform Inversion, Geophys., 74, 6, pp. 63-76, (2009)
[6] Brossier R., Operto S., Virieux J., Two-Dimensional Seismic Imaging of the Valhall Model from Synthetic OBC Data by Frequency-Domain Elastic Full-Waveform Inversion, SEG Tech. Program Expanded Abstr., 28, 1, pp. 2293-2297, (2009)
[7] Brossier R., Roux P., Seismic Imaging by Frequency-Domain Double-Beamforming Full-Waveform Inversion, Expanded Abstr., (2011)
[8] Bunks C., Salek F.M., Zaleski S., Chavent G., Multiscale SeismicWaveform Inversion, Geophys., 60, 5, pp. 1457-1473, (1995)
[9] Crase E., Pica A., Noble M., McDonald J., Tarantola A., Robust Elastic Non-Linear Waveform Inversion: Application to Real Data, Geophys., 55, pp. 527-538, (1990)
[10] Djikpesse H.A., Tarantola A., Multiparameter l <sub>1</sub> Norm Waveform Fitting: Interpretation of Gulf of Mexico Reflection Seismograms, Geophys., 64, 4, pp. 1023-1035, (1999)

← 1 2 3 →