6 results on '"Operto, Stéphane"'
Search Results
2. Building Initial Model for FWI From Ultra Long-Offset OBN Data by First-Arrival Traveltime + Slope Tomography
- Author
-
Sambolian, Serge, Górszczyk, A., Operto, Stéphane, Ribodetti, Alessandra, Tavakoli F., B., Virieux, Jean, Géoazur (GEOAZUR 7329), Institut national des sciences de l'Univers (INSU - CNRS)-Observatoire de la Côte d'Azur, COMUE Université Côte d'Azur (2015-2019) (COMUE UCA)-Université Côte d'Azur (UCA)-COMUE Université Côte d'Azur (2015-2019) (COMUE UCA)-Université Côte d'Azur (UCA)-Centre National de la Recherche Scientifique (CNRS)-Institut de Recherche pour le Développement (IRD [France-Sud]), Institut des Sciences de la Terre (ISTerre), Institut Français des Sciences et Technologies des Transports, de l'Aménagement et des Réseaux (IFSTTAR)-Institut national des sciences de l'Univers (INSU - CNRS)-Institut de recherche pour le développement [IRD] : UR219-Université Savoie Mont Blanc (USMB [Université de Savoie] [Université de Chambéry])-Centre National de la Recherche Scientifique (CNRS)-Université Grenoble Alpes [2016-2019] (UGA [2016-2019]), Institut national des sciences de l'Univers (INSU - CNRS)-Institut de Recherche pour le Développement (IRD [France-Sud])-Centre National de la Recherche Scientifique (CNRS)-Observatoire de la Côte d'Azur, Université Côte d'Azur (UCA)-Université Côte d'Azur (UCA), Université Grenoble Alpes (UGA)-Centre National de la Recherche Scientifique (CNRS)-Université Savoie Mont Blanc (USMB [Université de Savoie] [Université de Chambéry])-PRES Université de Grenoble-Institut de recherche pour le développement [IRD] : UR219-Institut national des sciences de l'Univers (INSU - CNRS)-Institut Français des Sciences et Technologies des Transports, de l'Aménagement et des Réseaux (IFSTTAR)-Université Joseph Fourier - Grenoble 1 (UJF), and Université Côte d'Azur (UCA)-Université Côte d'Azur (UCA)-Centre National de la Recherche Scientifique (CNRS)-Institut de Recherche pour le Développement (IRD [France-Sud])
- Subjects
Regional geology ,Offset (computer science) ,010504 meteorology & atmospheric sciences ,[SDU.STU.GP]Sciences of the Universe [physics]/Earth Sciences/Geophysics [physics.geo-ph] ,Engineering geology ,Gemology ,010502 geochemistry & geophysics ,01 natural sciences ,Glaciology ,14. Life underwater ,Economic geology ,Igneous petrology ,Geology ,Seismology ,ComputingMilieux_MISCELLANEOUS ,0105 earth and related environmental sciences ,Environmental geology - Abstract
Summary Ultra-long offset stationary-receiver acquisitions implemented with multi-component Ocean Bottom Nodes (OBN) are emerging as a suitable technology for deep offshore subsalt imaging by full waveform inversion (FWI). Ultra-long offsets allow to undershoot the deepest targeted structures with diving waves. In this framework, first-arrival traveltime + slope tomography (FASTT) provides a simple and efficient tool to build kinematically-accurate initial velocity model for full waveform inversion (FWI). In particular, the incorporation of slopes in the first-arrival tomography in addition to traveltimes removes artifacts in poorly illuminated areas and better reconstruct structural dips. Either the FASTT or the FWI velocity models can be used as reliable background models for prestack depth migration of co-incident towed-streamer data. However, the high resolution of the FWI velocity model facilitates the integrated interpretation of the reflectivity mapped by migration and the intermediate-scale velocity variations reconstructed by waveform inversion. We illustrate the workflow combining first-arrival tomography + FWI of ultra-long offset OBN data and migration of co-incident towed-streamer data with a challenging deep crustal case study in the eastern Nankai subduction zone.
- Published
- 2019
3. Velocity model building from seismic reflection data by full-waveform inversion.
- Author
-
Brossier, Romain, Operto, Stéphane, and Virieux, Jean
- Subjects
- *
SEISMOLOGY , *WAVE analysis , *QUANTITATIVE research , *IMAGING systems , *WAVENUMBER , *ALGORITHMS - Abstract
ABSTRACT Full-waveform inversion is re-emerging as a powerful data-fitting procedure for quantitative seismic imaging of the subsurface from wide-azimuth seismic data. This method is suitable to build high-resolution velocity models provided that the targeted area is sampled by both diving waves and reflected waves. However, the conventional formulation of full-waveform inversion prevents the reconstruction of the small wavenumber components of the velocity model when the subsurface is sampled by reflected waves only. This typically occurs as the depth becomes significant with respect to the length of the receiver array. This study first aims to highlight the limits of the conventional form of full-waveform inversion when applied to seismic reflection data, through a simple canonical example of seismic imaging and to propose a new inversion workflow that overcomes these limitations. The governing idea is to decompose the subsurface model as a background part, which we seek to update and a singular part that corresponds to some prior knowledge of the reflectivity. Forcing this scale uncoupling in the full-waveform inversion formalism brings out the transmitted wavepaths that connect the sources and receivers to the reflectors in the sensitivity kernel of the full-waveform inversion, which is otherwise dominated by the migration impulse responses formed by the correlation of the downgoing direct wavefields coming from the shot and receiver positions. This transmission regime makes full-waveform inversion amenable to the update of the long-to-intermediate wavelengths of the background model from the wide scattering-angle information. However, we show that this prior knowledge of the reflectivity does not prevent the use of a suitable misfit measurement based on cross-correlation, to avoid cycle-skipping issues as well as a suitable inversion domain as the pseudo-depth domain that allows us to preserve the invariant property of the zero-offset time. This latter feature is useful to avoid updating the reflectivity information at each non-linear iteration of the full-waveform inversion, hence considerably reducing the computational cost of the entire workflow. Prior information of the reflectivity in the full-waveform inversion formalism, a robust misfit function that prevents cycle-skipping issues and a suitable inversion domain that preserves the seismic invariant are the three key ingredients that should ensure well-posedness and computational efficiency of full-waveform inversion algorithms for seismic reflection data. [ABSTRACT FROM AUTHOR]
- Published
- 2015
- Full Text
- View/download PDF
4. A parametric analysis of two-dimensional elastic full waveform inversion of teleseismic data for lithospheric imaging.
- Author
-
Pageot, Damien, Operto, Stéphane, Vallée, Martin, Brossier, Romain, and Virieux, Jean
- Subjects
- *
WAVE analysis , *LITHOSPHERE , *SEISMOLOGY , *HIGH resolution imaging , *SHEAR waves , *QUASI-Newton methods , *SEISMIC tomography , *SCATTERING (Physics) , *THEORY of wave motion - Abstract
The development of dense networks of broad-band seismographs makes teleseismic data amenable to full-waveform inversion (FWI) methods for high-resolution lithospheric imaging. Compared to scattered-field migration, FWI seeks to involve the full seismic wavefield in the inversion. We present a parametric analysis of 2-D frequency-domain FWI in the framework of lithospheric imaging from teleseismic data to identify the main factors that impact on the quality of the reconstructed compressional (P)-wave and shear (S)-wave speed models. Compared to controlled-source seismology, the main adaptation of FWI to teleseismic configuration consists of the implementation with a scattered-filed formulation of plane-wave sources that impinge on the base of the lithospheric target located below the receiver network at an arbitrary incidence angle. Seismic modelling is performed with a hp-adaptive discontinuous Galerkin method on unstructured triangular mesh. A quasi-Newton inversion algorithm provides an approximate accounting for the Hessian operator, which contributes to reduce the footprint of the coarse acquisition geometry in the imaging. A versatile algorithm to compute the gradient of the misfit function with the adjoint-state method allows for abstraction between the forward-problem operators and the meshes that are during seismic modelling and inversion, respectively. An approximate correction for obliquity is derived for future application to real teleseismic data under the two-dimension approximation. Comparisons between the characteristic scales involved in exploration geophysics and in teleseismic seismology suggest that the resolution gain provided by full waveform technologies should be of the same order of magnitude for both applications. We first show the importance of the surface-reflected wavefield to dramatically improve the resolving power of FWI by combining tomography-like and migration-like imaging through the incorporation of the forward-scattered and the backscattered wavefields in the inversion. The resolution of FWI is assessed through checkerboard tests and confirms a resolution of the order of the wavelength for both the P and S speeds, when the full wavefield is incorporated in the inversion. Secondly, we show that computationally efficient strategies, which consist of decimating the number of frequency components involved in the inversion, do not apply to teleseismic acquisitions, because the scattering-angle bandwidth sampled by plane-wave sources can be narrow and coarsely sampled, compared to that provided by dense profiles of point sources in exploration seismology. The waveform inversion is less sensitive to the band of incidence angles spanned by the plane-wave sources and to the sampling of this band. However, the deficit of vertically propagating plane waves hampers the vertical resolution of planar layers. Aliasing artefacts created by coarse arrays of receivers are illustrated. We show how taking into account the Hessian in the inversion and the suitable management of frequencies in the inversion help to mitigate these artefacts. Acceptable reconstructions are shown for both the P- and S-wave speeds for a receiver spacing of up to 20 km in the 0.1–0.4 Hz frequency range. Building a reliable initial model for FWI is a highly non-linear problem in exploration seismology. We show how the low-frequency content of the teleseismic sources allow us to build accurate P- and S-wave speed models starting from simple vertical-gradient velocity models. All of these results are derived using a 2-D realistic synthetic experiment that was performed with noise-free data. Most of conclusions of this parametric study should apply in three dimensions. The impact of noise and the footprint of other experimental parameters ... [ABSTRACT FROM AUTHOR]
- Published
- 2013
- Full Text
- View/download PDF
5. On the footprint of anisotropy on isotropic full waveform inversion: the Valhall case study.
- Author
-
Prieux, Vincent, Brossier, Romain, Gholami, Yaser, Operto, Stéphane, Virieux, Jean, Barkved, O. I., and Kommedal, J. H.
- Subjects
ANISOTROPY ,ISOTROPY subgroups ,INVERSION (Geophysics) ,SEISMOLOGY ,TOMOGRAPHY ,ELASTIC wave propagation ,KINEMATICS ,MATHEMATICAL models - Abstract
SUMMARY The validity of isotropic approximation to perform acoustic full waveform inversion (FWI) of real wide-aperture anisotropic data can be questioned due to the intrinsic kinematic inconsistency between short- and large-aperture components of the data. This inconsistency is mainly related to the differences between the vertical and horizontal velocities in vertical-transverse isotropic (VTI) media. The footprint of VTI anisotropy on 2-D acoustic isotropic FWI is illustrated on a hydrophone data set of an ocean-bottom cable that was collected over the Valhall field in the North Sea. Multiscale FWI is implemented in the frequency domain by hierarchical inversions of increasing frequencies and decreasing aperture angles. The FWI models are appraised by local comparison with well information, seismic modelling, reverse-time migration (RTM) and source-wavelet estimation. A smooth initial VTI model parameterized by the vertical velocity V
0 and the Thomsen parameters δ and ε were previously developed by anisotropic reflection traveltime tomography. The normal moveout () and horizontal () velocity models were inferred from the anisotropic models to perform isotropic FWI. The VNMO models allows for an accurate match of short-spread reflection traveltimes, whereas the Vh model, after updating by first-arrival traveltime tomography (FATT), allows for an accurate match of first-arrival traveltimes. Ray tracing in the velocity models shows that the first 1.5 km of the medium are sampled by both diving waves and reflections, whereas the deeper structure at the reservoir level is mainly controlled by short-spread reflections. Starting from the initial anisotropic model and keeping fixed δ and ε models, anisotropic FWI allows us to build a vertical velocity model that matches reasonably well the well-log velocities. Isotropic FWI is performed using either the NMO model or the FATT model as initial model. In both cases, horizontal velocities are mainly reconstructed in the first 1.5 km of the medium. This suggests that the wide-aperture components of the data have a dominant control on the velocity estimation at these depths. These high velocities in the upper structure lead to low values of velocity in the underlying gas layers (either equal or lower than vertical velocities of the well log), and/or a vertical stretching of the structure at the reservoir level below the gas. This bias in the gas velocities and the mispositioning in depth of the deep reflectors, also shown in the RTM images, are required to match the deep reflections in the isotropic approximation and highlight the footprint of anisotropy in the isotropic FWI of long-offset data. Despite the significant differences between the anisotropic and isotropic FWI models, each of these models produce a nearly-equivalent match of the data, which highlights the ill-posedness of acoustic anisotropic FWI. Hence, we conclude with the importance of considering anisotropy in FWI of wide-aperture data to avoid bias in the velocity reconstructions and mispositioning in depth of reflectors. Designing a suitable parameterization of the VTI acoustic FWI is a central issue to manage the ill-posedness of the FWI. [ABSTRACT FROM AUTHOR]- Published
- 2011
- Full Text
- View/download PDF
6. FWT2D: A massively parallel program for frequency-domain full-waveform tomography of wide-aperture seismic data—Part 2: Numerical examples and scalability analysis
- Author
-
Sourbier, Florent, Operto, Stéphane, Virieux, Jean, Amestoy, Patrick, and L’Excellent, Jean-Yves
- Subjects
- *
SEISMOMETERS , *SYSTEMS design , *TOMOGRAPHY , *STOCHASTIC convergence , *SEISMOLOGY software , *MATHEMATICAL seismology , *SEISMOLOGY - Abstract
Abstract: This is the second paper in a two-part series that describes a massively parallel code that performs 2D full-waveform inversion of wide-aperture seismic data for imaging complex structures. We present several numerical validation of the full-waveform inversion code with both canonical and realistic synthetic examples. We illustrate how different multiscale strategies can be applied by either successive mono-frequency inversions or simultaneous multifrequency inversions and their impact on the convergence and the robustness of the inversion. We present a scalability analysis using a real marine data set recorded by a dense array of ocean bottom seismometers to image the crustal structure of a subduction zone. We obtained a speedup of 20 when using 50 processes on a PC cluster which allowed us to iteratively invert 13 frequencies of the full data set in less than 2 days. This computational performance will allow in the future more extensive analysis of full-waveform tomography methods when applied to representative case studies or when considering 3D geometries. [Copyright &y& Elsevier]
- Published
- 2009
- Full Text
- View/download PDF
Catalog
Discovery Service for Jio Institute Digital Library
For full access to our library's resources, please sign in.