Three-dimensional frequency-domain full-waveform inversion with an iterative solver

Geophysics ◽  
2009 ◽  
Vol 74 (6) ◽  
pp. WCC149-WCC157 ◽  
Author(s):  
René-Édouard Plessix

With the acquisition of wide-aperture seismic data sets, full-waveform inversion is an attractive method for deriving velocity models. Three-dimensional implementations require an efficient solver for the wave equation. Computing 3D time-harmonic responses with a frequency-domain solver is complicated because a large linear system with negative and positive eigenvalues must be solved. Time-domain schemes are an alternative. Nevertheless, existing frequency-domain iterative solvers with an efficient preconditioner are a viable option when full-waveform inversion is formulated in the frequency domain. An iterative solver with a multigrid preconditioner is competitive because of a high-order spatial discretization. Numerical examples illustrated the efficiency of the iterative solvers. Three dimensional full-waveform inversion was then studied in the context of deep-water ocean-bottom seismometer acquisition. Three dimensional synthetic data inversion results showed the behavior of full-waveform inversion with respect to the initial model and the minimum frequency available in the data set. Results on a 3D real ocean-bottom seismometer data set demonstrated the relevance of full-waveform inversion, especially to image the shallow part of the model.

Geophysics ◽  
2021 ◽  
pp. 1-52
Author(s):  
Yuzhu Liu ◽  
Xinquan Huang ◽  
Jizhong Yang ◽  
Xueyi Liu ◽  
Bin Li ◽  
...  

Thin sand-mud-coal interbedded layers and multiples caused by shallow water pose great challenges to conventional 3D multi-channel seismic techniques used to detect the deeply buried reservoirs in the Qiuyue field. In 2017, a dense ocean-bottom seismometer (OBS) acquisition program acquired a four-component dataset in East China Sea. To delineate the deep reservoir structures in the Qiuyue field, we applied a full-waveform inversion (FWI) workflow to this dense four-component OBS dataset. After preprocessing, including receiver geometry correction, moveout correction, component rotation, and energy transformation from 3D to 2D, a preconditioned first-arrival traveltime tomography based on an improved scattering integral algorithm is applied to construct an initial P-wave velocity model. To eliminate the influence of the wavelet estimation process, a convolutional-wavefield-based objective function for the preprocessed hydrophone component is used during acoustic FWI. By inverting the waveforms associated with early arrivals, a relatively high-resolution underground P-wave velocity model is obtained, with updates at 2.0 km and 4.7 km depth. Initial S-wave velocity and density models are then constructed based on their prior relationships to the P-wave velocity, accompanied by a reciprocal source-independent elastic full-waveform inversion to refine both velocity models. Compared to a traditional workflow, guided by stacking velocity analysis or migration velocity analysis, and using only the pressure component or other single-component, the workflow presented in this study represents a good approach for inverting the four-component OBS dataset to characterize sub-seafloor velocity structures.


2016 ◽  
Author(s):  
Mikhail Belonosov ◽  
Vladimir Tcheverda ◽  
Dmitry Neklyudov ◽  
Victor Kostin ◽  
Maxim Dmitriev

Solid Earth ◽  
2019 ◽  
Vol 10 (3) ◽  
pp. 765-784 ◽  
Author(s):  
Andrzej Górszczyk ◽  
Stéphane Operto ◽  
Laure Schenini ◽  
Yasuhiro Yamada

Abstract. Imaging via pre-stack depth migration (PSDM) of reflection towed-streamer multichannel seismic (MCS) data at the scale of the whole crust is inherently difficult. This is because the depth penetration of the seismic wavefield is controlled, firstly, by the acquisition design, such as streamer length and air-gun source configuration, and secondly by the complexity of the crustal structure. Indeed, the limited length of the streamer makes the estimation of velocities from deep targets challenging due to the velocity–depth ambiguity. This problem is even more pronounced when processing 2-D seismic data due to the lack of multi-azimuthal coverage. Therefore, in order to broaden our knowledge about the deep crust using seismic methods, we present the development of specific imaging workflows that integrate different seismic data. Here we propose the combination of velocity model building using (i) first-arrival tomography (FAT) and full-waveform inversion (FWI) of wide-angle, long-offset data collected by stationary ocean-bottom seismometers (OBSs) and (ii) PSDM of short-spread towed-streamer MCS data for reflectivity imaging, with the former velocity model as a background model. We present an application of such a workflow to seismic data collected by the Japan Agency for Marine-Earth Science and Technology (JAMSTEC) and the Institut Français de Recherche pour l'Exploitation de la Mer (IFREMER) in the eastern Nankai Trough (Tokai area) during the 2000–2001 Seize France Japan (SFJ) experiment. We show that the FWI model, although derived from OBS data, provides an acceptable background velocity field for the PSDM of the MCS data. From the initial PSDM, we refine the FWI background velocity model by minimizing the residual move-outs (RMOs) picked in the pre-stack-migrated volume through slope tomography (ST), from which we generate a better-focused migrated image. Such integration of different seismic datasets and leading-edge imaging techniques led to greatly improved imaging at different scales. That is, large to intermediate crustal units identified in the high-resolution FWI velocity model extensively complement the short-wavelength reflectivity inferred from the MCS data to better constrain the structural factors controlling the geodynamics of the Nankai Trough.


Sign in / Sign up

Export Citation Format

Share Document