scholarly journals spyro: a Firedrake-based wave propagation and full waveform inversion finite element solver

2021 ◽  
Author(s):  
Keith J. Roberts ◽  
Alexandre Olender ◽  
Lucas Franceschini ◽  
Robert C. Kirby ◽  
Rafael S. Gioria ◽  
...  
Keyword(s):  
Finite Element ◽  
Wave Propagation ◽  
Model Building ◽  
Seismic Velocity ◽  
Waveform Inversion ◽  
Higher Order ◽  
Full Waveform Inversion ◽  
State Equations ◽  
Full Waveform ◽  
Triangular Elements

Abstract. In this article, we introduce spyro, a software stack to solve acoustic wave propagation in heterogeneous domains and perform full waveform inversion (FWI) employing the finite element framework from Firedrake, a high-level Python package for the automated solution of partial differential equations using the finite element method. The capability of the software is demonstrated by using a continuous Galerkin approach to perform FWI for seismic velocity model building, considering realistic geophysics examples. A time-domain FWI approach is detailed that uses meshes composed of variably sized triangular elements to discretize the domain. To resolve both the forward and adjoint-state equations, and to calculate a mesh-independent gradient associated with the FWI process, a fully-explicit, variable higher-order (up to degree k = 5 in 2D and k = 3 in 3D) mass lumping method is used. We show that, by adapting the triangular elements to the expected peak source frequency and properties of the wavefield (e.g., local P-wavespeed) and by leveraging higher-order basis functions, the number of degrees-of-freedom necessary to discretize the domain can be reduced. Results from wave simulations and FWIs in both 2D and 3D highlight our developments and demonstrate the benefits and challenges with using triangular meshes adapted to the material properties.

scholarly journals Accelerating Numerical Wave Propagation by Wavefield Adapted Meshes, Part II: Full-Waveform Inversion

2019 ◽  
Author(s):  
Solvi Thrastarson ◽  
Martin van Driel ◽  
Lion Krischer ◽  
Dirk-Philip van Herwaarden ◽  
Christian Boehm ◽  
...  
Keyword(s):  
Wave Propagation ◽  
Waveform Inversion ◽  
Full Waveform Inversion ◽  
Full Waveform ◽  
Numerical Wave

Full-waveform inversion for full-wavefield imaging: Decades in the making

2021 ◽  
Vol 40 (5) ◽  
pp. 324-334
Author(s):  
Rongxin Huang ◽  
Zhigang Zhang ◽  
Zedong Wu ◽  
Zhiyuan Wei ◽  
Jiawei Mei ◽  
...  
Keyword(s):  
Model Building ◽  
Seismic Imaging ◽  
Structural Information ◽  
Waveform Inversion ◽  
Velocity Model ◽  
Data Sets ◽  
Full Waveform Inversion ◽  
Data Types ◽  
Velocity Models ◽  
Full Waveform

Seismic imaging using full-wavefield data that includes primary reflections, transmitted waves, and their multiples has been the holy grail for generations of geophysicists. To be able to use the full-wavefield data effectively requires a forward-modeling process to generate full-wavefield data, an inversion scheme to minimize the difference between modeled and recorded data, and, more importantly, an accurate velocity model to correctly propagate and collapse energy of different wave modes. All of these elements have been embedded in the framework of full-waveform inversion (FWI) since it was proposed three decades ago. However, for a long time, the application of FWI did not find its way into the domain of full-wavefield imaging, mostly owing to the lack of data sets with good constraints to ensure the convergence of inversion, the required compute power to handle large data sets and extend the inversion frequency to the bandwidth needed for imaging, and, most significantly, stable FWI algorithms that could work with different data types in different geologic settings. Recently, with the advancement of high-performance computing and progress in FWI algorithms at tackling issues such as cycle skipping and amplitude mismatch, FWI has found success using different data types in a variety of geologic settings, providing some of the most accurate velocity models for generating significantly improved migration images. Here, we take a step further to modify the FWI workflow to output the subsurface image or reflectivity directly, potentially eliminating the need to go through the time-consuming conventional seismic imaging process that involves preprocessing, velocity model building, and migration. Compared with a conventional migration image, the reflectivity image directly output from FWI often provides additional structural information with better illumination and higher signal-to-noise ratio naturally as a result of many iterations of least-squares fitting of the full-wavefield data.

Full waveform inversion: early start of model building process to resolve the complex near surface model in the Exmouth sub-basin, North West Shelf Australia

2018 ◽  
Vol 58 (2) ◽  
pp. 884
Author(s):  
Lianping Zhang ◽  
Haryo Trihutomo ◽  
Yuelian Gong ◽  
Bee Jik Lim ◽  
Alexander Karvelas
Keyword(s):  
Model Building ◽  
Waveform Inversion ◽  
Velocity Model ◽  
Surface Velocity ◽  
Image Point ◽  
Surface Model ◽  
Full Waveform Inversion ◽  
Near Surface ◽  
North West ◽  
Full Waveform

The Schlumberger Multiclient Exmouth 3D survey was acquired over the Exmouth sub-basin, North West Shelf Australia and covers 12 600 km2. One of the primary objectives of this survey was to produce a wide coverage of high quality imaging with advanced processing technology within an agreed turnaround time. The complexity of the overburden was one of the imaging challenges that impacted the structuration and image quality at the reservoir level. Unlike traditional full-waveform inversion (FWI) workflow, here, FWI was introduced early in the workflow in parallel with acquisition and preprocessing to produce a reliable near surface velocity model from a smooth starting model. FWI derived an accurate and detailed near surface model, which subsequently benefitted the common image point (CIP) tomography model updates through to the deeper intervals. The objective was to complete the FWI model update for the overburden concurrently with the demultiple stages hence reflection time CIP tomography could start with a reasonably good velocity model upon completion of the demultiple process.

Assisting salt model building with reflection full-waveform inversion

2019 ◽  
Vol 7 (2) ◽  
pp. SB43-SB52 ◽  
Author(s):  
Adriano Gomes ◽  
Joe Peterson ◽  
Serife Bitlis ◽  
Chengliang Fan ◽  
Robert Buehring
Keyword(s):  
Model Building ◽  
Waveform Inversion ◽  
Low Frequency ◽  
Velocity Model ◽  
Full Waveform Inversion ◽  
Reflection Data ◽  
Full Waveform ◽  
Velocity Model Building ◽  
Wave Limit ◽  
Rabbit Ears

Inverting for salt geometry using full-waveform inversion (FWI) is a challenging task, mostly due to the lack of extremely low-frequency signal in the seismic data, the limited penetration depth of diving waves using typical acquisition offsets, and the difficulty in correctly modeling the amplitude (and kinematics) of reflection events associated with the salt boundary. However, recent advances in reflection FWI (RFWI) have allowed it to use deep reflection data, beyond the diving-wave limit, by extracting the tomographic term of the FWI reflection update, the so-called rabbit ears. Though lacking the resolution to fully resolve salt geometry, we can use RFWI updates as a guide for refinements in the salt interpretation, adding a partially data-driven element to salt velocity model building. In addition, we can use RFWI to update sediment velocities in complex regions surrounding salt, where ray-based approaches typically struggle. In reality, separating the effects of sediment velocity errors from salt geometry errors is not straightforward in many locations. Therefore, iterations of RFWI plus salt scenario tests may be necessary. Although it is still not the fully automatic method that has been envisioned for FWI, this combined approach can bring significant improvement to the subsalt image, as we examine on field data examples from the Gulf of Mexico.

New tools for 2D full-waveform inversion: applications on Brazilian Pre-Salt velocity model from Santos Basin

2020 ◽  
Author(s):  
Henrique Santos ◽  
Claus Eikmeier ◽  
Ernani Volpe
Keyword(s):  
Oil And Gas ◽  
Model Building ◽  
Waveform Inversion ◽  
Numerical Algorithms ◽  
Velocity Model ◽  
Southeastern Brazil ◽  
Full Waveform Inversion ◽  
Case Scenario ◽  
Velocity Inversion ◽  
Full Waveform

<p>In this work, we present full-waveform inversion (FWI) results of a typical Brazilian Pre-Salt model (Santos Basin) using new open-source tools. The large accumulations of oil with excellent quality and high commercial value discovered in the pre-salt carbonates of southeastern Brazil, especially in the Santos Basin, have made this province one of the most prospective in the world. Velocity model building in areas of highly complex geology (like the Santos Basin) remains a challenging step in seismic processing. FWI proved to be an efficient tool for the determination of high-resolution details in multiparameter models of complex subsurface structures, and it has been applied in different geophysical problem scales. However, since FWI is a computationally and mathematically challenging problem, many issues remain open, such as more efficient ways to deal with multiparameter inversion problems such crosstalk and different orders of magnitude in the seismic signal for different classes of parameters. Inversions for more than one class of parameters are of particular importance in the estimation of the physical properties of rocks (poroacoustic or poroelastic applications), for example, to monitoring oil and gas reservoirs and for monitoring the injection of carbon dioxide into geological structures. Also, programming complex numerical algorithms for each application is time-consuming and often evades the expertise of researchers from the geoscientific community. In this sense, a high-level computational tool for simulations and inversions would greatly improve the working time for researchers. Existing finite difference based FWI tools such as Devito, and finite elements based partial differential equations (PDE) solvers tools such as FEniCS and Firedrake are being explored and used for these purposes. In this work, we initially present an FWI acoustic isotropic inversion test (velocity inversion only), performed with the Devito software while a particular code is being developed in FEniCS and Firedrake computer programs. Devito is also a new and under development software and therefore must be tested under different conditions. Our first numerical results indicate the potential of using freely available computational programs in a real case scenario.</p>

Sign in / Sign up

Export Citation Format

Share Document