scholarly journals Optimal full-waveform inversion strategy for marine data in azimuthally rotated elastic orthorhombic media

Geophysics ◽  
2018 ◽  
Vol 83 (4) ◽  
pp. R307-R320 ◽  
Author(s):  
Ju-Won Oh ◽  
Tariq Alkhalifah

The orthorhombic (ORT) anisotropic description of earth layers can allow the capture of much of the earth’s anisotropic complexity. The inversion for high-resolution azimuthal variation of anisotropy is important for reservoir characterization, among other applications. A high-resolution description of the azimuth of fractures can help us to predict flow preferences. To verify the feasibility of multiparameter full-waveform inversion (FWI) for marine data assuming azimuthally rotated elastic ORT media, we have analyzed the radiation patterns and gradient directions of ORT parameters to the reflection data. First, we express the gradient direction of the ORT parameters considering the azimuthal rotation of the symmetric planes. Then, to support our observations in the gradient direction, the radiation patterns of the partial derivative wavefields from each parameter perturbation are also derived under the rotated elastic ORT assumption. To find an optimal parameterization, we compare three different parameterizations: monoclinic, velocity-based, and hierarchical parameterizations. Then, we suggest an optimal multistage update strategy by analyzing the behavior of the rotation angle as a FWI target. To analyze the trade-off among parameters in different parameterizations, we calculate the gradient direction from a hockey-puck model, in which each parameter is perturbed at the different location on a horizontal layer. The trade-off analysis supports that the hierarchical parameterization provides us with more opportunities to build up subsurface models with less trade-off between parameters and less influence of the azimuthal rotation of ORT anisotropy. The feasibility of the proposed FWI strategy is examined using synthetic marine streamer data from a simple 3D reservoir model with a fractured layer.

Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. WA137-WA146
Author(s):  
Zhen-dong Zhang ◽  
Tariq Alkhalifah

Reservoir characterization is an essential component of oil and gas production, as well as exploration. Classic reservoir characterization algorithms, deterministic and stochastic, are typically based on stacked images and rely on simplifications and approximations to the subsurface (e.g., assuming linearized reflection coefficients). Elastic full-waveform inversion (FWI), which aims to match the waveforms of prestack seismic data, potentially provides more accurate high-resolution reservoir characterization from seismic data. However, FWI can easily fail to characterize deep-buried reservoirs due to illumination limitations. We have developed a deep learning-aided elastic FWI strategy using observed seismic data and available well logs in the target area. Five facies are extracted from the well and then connected to the inverted P- and S-wave velocities using trained neural networks, which correspond to the subsurface facies distribution. Such a distribution is further converted to the desired reservoir-related parameters such as velocities and anisotropy parameters using a weighted summation. Finally, we update these estimated parameters by matching the resulting simulated wavefields to the observed seismic data, which corresponds to another round of elastic FWI aided by the a priori knowledge gained from the predictions of machine learning. A North Sea field data example, the Volve Oil Field data set, indicates that the use of facies as prior knowledge helps resolve the deep-buried reservoir target better than the use of only seismic data.


Geophysics ◽  
2016 ◽  
Vol 81 (4) ◽  
pp. U25-U38 ◽  
Author(s):  
Nuno V. da Silva ◽  
Andrew Ratcliffe ◽  
Vetle Vinje ◽  
Graham Conroy

Parameterization lies at the center of anisotropic full-waveform inversion (FWI) with multiparameter updates. This is because FWI aims to update the long and short wavelengths of the perturbations. Thus, it is important that the parameterization accommodates this. Recently, there has been an intensive effort to determine the optimal parameterization, centering the fundamental discussion mainly on the analysis of radiation patterns for each one of these parameterizations, and aiming to determine which is best suited for multiparameter inversion. We have developed a new parameterization in the scope of FWI, based on the concept of kinematically equivalent media, as originally proposed in other areas of seismic data analysis. Our analysis is also based on radiation patterns, as well as the relation between the perturbation of this set of parameters and perturbation in traveltime. The radiation pattern reveals that this parameterization combines some of the characteristics of parameterizations with one velocity and two Thomsen’s parameters and parameterizations using two velocities and one Thomsen’s parameter. The study of perturbation of traveltime with perturbation of model parameters shows that the new parameterization is less ambiguous when relating these quantities in comparison with other more commonly used parameterizations. We have concluded that our new parameterization is well-suited for inverting diving waves, which are of paramount importance to carry out practical FWI successfully. We have demonstrated that the new parameterization produces good inversion results with synthetic and real data examples. In the latter case of the real data example from the Central North Sea, the inverted models show good agreement with the geologic structures, leading to an improvement of the seismic image and flatness of the common image gathers.


Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE101-VE117 ◽  
Author(s):  
Hafedh Ben-Hadj-Ali ◽  
Stéphane Operto ◽  
Jean Virieux

We assessed 3D frequency-domain (FD) acoustic full-waveform inversion (FWI) data as a tool to develop high-resolution velocity models from low-frequency global-offset data. The inverse problem was posed as a classic least-squares optimization problem solved with a steepest-descent method. Inversion was applied to a few discrete frequencies, allowing management of a limited subset of the 3D data volume. The forward problem was solved with a finite-difference frequency-domain method based on a massively parallel direct solver, allowing efficient multiple-shot simulations. The inversion code was fully parallelized for distributed-memory platforms, taking advantage of a domain decomposition of the modeled wavefields performed by the direct solver. After validation on simple synthetic tests, FWI was applied to two targets (channel and thrust system) of the 3D SEG/EAGE overthrust model, corresponding to 3D domains of [Formula: see text] and [Formula: see text], respectively. The maximum inverted frequencies are 15 and [Formula: see text] for the two applications. A maximum of 30 dual-core biprocessor nodes with [Formula: see text] of shared memory per node were used for the second target. The main structures were imaged successfully at a resolution scale consistent with the inverted frequencies. Our study confirms the feasibility of 3D frequency-domain FWI of global-offset data on large distributed-memory platforms to develop high-resolution velocity models. These high-velocity models may provide accurate macromodels for wave-equation prestack depth migration.


Geophysics ◽  
2020 ◽  
Vol 85 (6) ◽  
pp. R553-R563
Author(s):  
Sagar Singh ◽  
Ilya Tsvankin ◽  
Ehsan Zabihi Naeini

The nonlinearity of full-waveform inversion (FWI) and parameter trade-offs can prevent convergence toward the actual model, especially for elastic anisotropic media. The problems with parameter updating become particularly severe if ultra-low-frequency seismic data are unavailable, and the initial model is not sufficiently accurate. We introduce a robust way to constrain the inversion workflow using borehole information obtained from well logs. These constraints are included in the form of rock-physics relationships for different geologic facies (e.g., shale, sand, salt, and limestone). We develop a multiscale FWI algorithm for transversely isotropic media with a vertical symmetry axis (VTI media) that incorporates facies information through a regularization term in the objective function. That term is updated during the inversion by using the models obtained at the previous inversion stage. To account for lateral heterogeneity between sparse borehole locations, we use an image-guided smoothing algorithm. Numerical testing for structurally complex anisotropic media demonstrates that the facies-based constraints may ensure the convergence of the objective function towards the global minimum in the absence of ultra-low-frequency data and for simple (even 1D) initial models. We test the algorithm on clean data and on surface records contaminated by Gaussian noise. The algorithm also produces a high-resolution facies model, which should be instrumental in reservoir characterization.


Sign in / Sign up

Export Citation Format

Share Document