Joint electromagnetic and seismic inversion using structural constraints

Geophysics ◽  
2009 ◽  
Vol 74 (6) ◽  
pp. R99-R109 ◽  
Author(s):  
Wenyi Hu ◽  
Aria Abubakar ◽  
Tarek M. Habashy
Keyword(s):  
Seismic Data ◽  
Joint Inversion ◽  
Monitoring And Evaluation ◽  
Structural Similarity ◽  
Compressional Wave ◽  
Seismic Inversion ◽  
P Wave ◽  
Optimization Approach ◽  
Structural Constraints ◽  
Inversion Algorithm

We have developed a frequency-domain joint electromagnetic (EM) and seismic inversion algorithm for reservoir evaluation and exploration applications. EM and seismic data are jointly inverted using a cross-gradient constraint that enforces structural similarity between the conductivity image and the compressional wave (P-wave) velocity image. The inversion algorithm is based on a Gauss-Newton optimization approach. Because of the ill-posed nature of the inverse problem, regularization is used to constrain the solution. The multiplicative regularization technique selects the regularization parameters automatically, improving the robustness of the algorithm. A multifrequency data-weighting scheme prevents the high-frequency data from dominating the inversion process. When the joint-inversion algorithm is applied in integrating marine controlled-source electromagnetic data with surface seismic data for subsea reservoir exploration applications and in integrating crosswell EM and sonic data for reservoir monitoring and evaluation applications, results improve significantly over those obtained from separate EM or seismic inversions.

Structural joint inversion on irregular meshes

2019 ◽  
Vol 220 (3) ◽  
pp. 1995-2008 ◽  
Author(s):  
C Jordi ◽  
J Doetsch ◽  
T Günther ◽  
C Schmelzbach ◽  
H Maurer ◽  
...  
Keyword(s):  
Electrical Resistivity ◽  
Joint Inversion ◽  
A Priori ◽  
Structural Similarity ◽  
P Wave ◽  
Parameter Distribution ◽  
Data Sets ◽  
Inversion Algorithm ◽  
Structural Joint ◽  
Irregular Meshes

SUMMARY Structural joint inversion of several data sets on an irregular mesh requires appropriate coupling operators. To date, joint inversion algorithms are primarily designed for the use on regular rectilinear grids and impose structural similarity in the direct neighbourhood of a cell only. We introduce a novel scheme for calculating cross-gradient operators based on a correlation model that allows to define the operator size by imposing physical length scales. We demonstrate that the proposed cross-gradient operators are largely decoupled from the discretization of the modelling domain, which is particularly important for irregular meshes where cell sizes vary. Our structural joint inversion algorithm is applied to a synthetic electrical resistivity tomography and ground penetrating radar 3-D cross-well experiment aiming at imaging two anomalous bodies and extracting the parameter distribution of the geostatistical background models. For both tasks, joint inversion produced superior results compared with individual inversions of the two data sets. Finally, we applied structural joint inversion to two field data sets recorded over a karstified limestone area. By including geological a priori information via the correlation-based operators into the joint inversion, we find P-wave velocity and electrical resistivity tomograms that are in accordance with the expected subsurface geology.

Joint Inversion of Electromagnetic and Seismic Data Based on Structural Constraints Using Variational Born Iteration Method

2018 ◽  
Vol 56 (1) ◽  
pp. 436-445 ◽  
Author(s):  
Tian Lan ◽  
Hai Liu ◽  
Na Liu ◽  
Jinghe Li ◽  
Feng Han ◽  
...  
Keyword(s):  
Seismic Data ◽  
Iteration Method ◽  
Joint Inversion ◽  
Structural Constraints

Simultaneous inversion of prestack seismic data for rock properties using simulated annealing

Geophysics ◽  
2002 ◽  
Vol 67 (6) ◽  
pp. 1877-1885 ◽  
Author(s):  
Xin‐Quan Ma
Keyword(s):  
Simulated Annealing ◽  
Seismic Data ◽  
Low Frequency ◽  
Optimization Procedure ◽  
P Wave ◽  
Rock Properties ◽  
Inversion Method ◽  
Inversion Algorithm ◽  
Simultaneous Inversion ◽  
Reflection Seismic

A new prestack inversion algorithm has been developed to simultaneously estimate acoustic and shear impedances from P‐wave reflection seismic data. The algorithm uses a global optimization procedure in the form of simulated annealing. The goal of optimization is to find a global minimum of the objective function, which includes the misfit between synthetic and observed prestack seismic data. During the iterative inversion process, the acoustic and shear impedance models are randomly perturbed, and the synthetic seismic data are calculated and compared with the observed seismic data. To increase stability, constraints have been built into the inversion algorithm, using the low‐frequency impedance and background Vs/Vp models. The inversion method has been successfully applied to synthetic and field data examples to produce acoustic and shear impedances comparable to log data of similar bandwidth. The estimated acoustic and shear impedances can be combined to derive other elastic parameters, which may be used for identifying of lithology and fluid content of reservoirs.

Accuracy Comparison of Elastic Impedance Formula

2012 ◽  
Vol 466-467 ◽  
pp. 400-404
Author(s):  
Jin Zhang ◽  
Huai Shan Liu ◽  
Si You Tong ◽  
Lin Fei Wang ◽  
Bing Xu
Keyword(s):  
Seismic Data ◽  
Poisson Ratio ◽  
Seismic Inversion ◽  
P Wave ◽  
S Wave ◽  
Elastic Parameters ◽  
Accuracy Comparison ◽  
Elastic Impedance ◽  
Prestack Seismic Inversion ◽  
Lamé Coefficients

Elastic impedance (EI) inversion is one of the prestack seismic inversion methods, which can obtain P-wave and S-wave velocity, density, Poisson ratio, Lame coefficients and other elastic parameters. But there have been many EI formulas nowadays, so which formula should be used in inversion is an urgent problem. This paper divides these formulas into two categories, and use several forward modeling to test the accuracy of these EI formulas. It shows that using the first kind of EI formulas in near offset seismic data can get high precision results.

Structurally constrained impedance inversion

2016 ◽  
Vol 4 (4) ◽  
pp. T577-T589 ◽  
Author(s):  
Haitham Hamid ◽  
Adam Pidlisecky
Keyword(s):  
Objective Function ◽  
Seismic Data ◽  
Field Data ◽  
A Priori ◽  
Structural Constraints ◽  
Inversion Algorithm ◽  
Data Set ◽  
Impedance Inversion ◽  
Single Trace ◽  
Complex Geology

In complex geology, the presence of highly dipping structures can complicate impedance inversion. We have developed a structurally constrained inversion in which a computationally well-behaved objective function is minimized subject to structural constraints. This approach allows the objective function to incorporate structural orientation in the form of dips into our inversion algorithm. Our method involves a multitrace impedance inversion and a rotation of an orthogonal system of derivative operators. Local dips used to constrain the derivative operators were estimated from migrated seismic data. In addition to imposing structural constraints on the inversion model, this algorithm allows for the inclusion of a priori knowledge from boreholes. We investigated this algorithm on a complex synthetic 2D model as well as a seismic field data set. We compared the result obtained with this approach with the results from single trace-based inversion and laterally constrained inversion. The inversion carried out using dip information produces a model that has higher resolution that is more geologically realistic compared with other methods.

A joint inversion algorithm to process geoelectric and sutface wave seismic data. Part I: basic ideas1

1995 ◽  
Vol 43 (2) ◽  
pp. 135-156 ◽  
Author(s):  
A. Hering ◽  
R. Misiek ◽  
A. Gyulai ◽  
T. Ormos ◽  
M. Dobroka ◽  
...  
Keyword(s):  
Seismic Data ◽  
Joint Inversion ◽  
Inversion Algorithm

Data-driven Prestack Waveform Inversion Using Genetic Algorithm- Methodology and Examples

2021 ◽  
pp. 1-97
Author(s):  
Lingxiao Jia ◽  
Subhashis Mallick ◽  
Cheng Wang
Keyword(s):  
Genetic Algorithm ◽  
Seismic Data ◽  
Waveform Inversion ◽  
Seismic Inversion ◽  
P Wave ◽  
Data Driven ◽  
Initial Model ◽  
Reservoir Properties ◽  
Well Control ◽  
Seismic Waveform

The choice of an initial model for seismic waveform inversion is important. In matured exploration areas with adequate well control, we can generate a suitable initial model using well information. However, in new areas where well control is sparse or unavailable, such an initial model is compromised and/or biased by the regions with more well controls. Even in matured exploration areas, if we use time-lapse seismic data to predict dynamic reservoir properties, an initial model, that we obtain from the existing preproduction wells could be incorrect. In this work, we outline a new methodology and workflow for a nonlinear prestack isotropic elastic waveform inversion. We call this method a data driven inversion, meaning that we derive the initial model entirely from the seismic data without using any well information. By assuming a locally horizonal stratification for every common midpoint and starting from the interval P-wave velocity, estimated entirely from seismic data, our method generates pseudo wells by running a two-pass one-dimensional isotropic elastic prestack waveform inversion that uses the reflectivity method for forward modeling and genetic algorithm for optimization. We then use the estimated pseudo wells to build the initial model for seismic inversion. By applying this methodology to real seismic data from two different geological settings, we demonstrate the usefulness of our method. We believe that our new method is potentially applicable for subsurface characterization in areas where well information is sparse or unavailable. Additional research is however necessary to improve the compute-efficiency of the methodology.

Simultaneous inversion of PP and PS seismic data

Geophysics ◽  
2006 ◽  
Vol 71 (3) ◽  
pp. R1-R10 ◽  
Author(s):  
Helene Hafslund Veire ◽  
Martin Landrø
Keyword(s):  
Seismic Data ◽  
Reservoir Characterization ◽  
Model Building ◽  
Joint Inversion ◽  
P Wave ◽  
Inversion Method ◽  
System Of Equations ◽  
S Wave ◽  
Data Set ◽  
Simultaneous Inversion

Elastic parameters derived from seismic data are valuable input for reservoir characterization because they can be related to lithology and fluid content of the reservoir through empirical relationships. The relationship between physical properties of rocks and fluids and P-wave seismic data is nonunique. This leads to large uncertainties in reservoir models derived from P-wave seismic data. Because S- waves do not propagate through fluids, the combined use of P-and S-wave seismic data might increase our ability to derive fluid and lithology effects from seismic data, reducing the uncertainty in reservoir characterization and thereby improving 3D reservoir model-building. We present a joint inversion method for PP and PS seismic data by solving approximated linear expressions of PP and PS reflection coefficients simultaneously using a least-squares estimation algorithm. The resulting system of equations is solved by singular-value decomposition (SVD). By combining the two independent measurements (PP and PS seismic data), we stabilize the system of equations for PP and PS seismic data separately, leading to more robust parameter estimation. The method does not require any knowledge of PP and PS wavelets. We tested the stability of this joint inversion method on a 1D synthetic data set. We also applied the methodology to North Sea multicomponent field data to identify sand layers in a shallow formation. The identified sand layers from our inverted sections are consistent with observations from nearby well logs.

Analysis of prior models for a blocky inversion of seismic AVA data

Geophysics ◽  
2010 ◽  
Vol 75 (3) ◽  
pp. C25-C35 ◽  
Author(s):  
Ulrich Theune ◽  
Ingrid Østgård Jensås ◽  
Jo Eidsvik
Keyword(s):  
Synthetic Data ◽  
Iterative Solution ◽  
Seismic Inversion ◽  
Inversion Method ◽  
Optimization Approach ◽  
Additional Parameter ◽  
Inversion Algorithm ◽  
Borehole Data ◽  
Solution Approach ◽  
Exploration And Production

Resolving thinner layers and focusing layer boundaries better in inverted seismic sections are important challenges in exploration and production seismology to better identify a potential drilling target. Many seismic inversion methods are based on a least-squares optimization approach that can intrinsically lead to unfocused transitions between adjacent layers. A Bayesian seismic amplitude variation with angle (AVA) inversion algorithm forms sharper boundaries between layers when enforcing sparseness in the vertical gradients of the inversion results. The underlying principle is similar to high-resolution processing algorithms and has been adapted from digital-image-sharpening algorithms. We have investigated the Cauchy and Laplace statistical distributions for their potential to improve contrasts betweenlayers. An inversion algorithm is derived statistically from Bayes’ theorem and results in a nonlinear problem that requires an iterative solution approach. Bayesian inversions require knowledge of certain statistical properties of the model we want to invert for. The blocky inversion method requires an additional parameter besides the usual properties for a multivariate covariance matrix, which we can estimate from borehole data. Tests on synthetic and field data show that the blocky inversion algorithm can detect and enhance layer boundaries in seismic inversions by effectively suppressing side lobes. The analysis of the synthetic data suggests that the Laplace constraint performs more reliably, whereas the Cauchy constraint may not find the optimum solution by converging to a local minimum of the cost function and thereby introducing some numerical artifacts.

Joint acoustic full-waveform inversion of crosshole seismic and ground-penetrating radar data in the frequency domain

Geophysics ◽  
2017 ◽  
Vol 82 (6) ◽  
pp. H41-H56 ◽  
Author(s):  
Xuan Feng ◽  
Qianci Ren ◽  
Cai Liu ◽  
Xuebing Zhang
Keyword(s):  
Frequency Domain ◽  
Ground Penetrating Radar ◽  
Joint Inversion ◽  
Waveform Inversion ◽  
Structural Similarity ◽  
P Wave ◽  
Full Waveform Inversion ◽  
Full Waveform ◽  
Frequency Domain Methods ◽  
Ground Penetrating

Integrating crosshole ground-penetrating radar (GPR) with seismic methods is an efficient way to reduce the uncertainty and ambiguity of data interpretation in shallow geophysical investigations. We have developed a new approach for joint full-waveform inversion (FWI) of crosshole seismic and GPR data in the frequency domain to improve the inversion results of both FWI methods. In a joint objective function, three geophysical parameters (P-wave velocity, permittivity, and conductivity) are effectively connected by three weighted cross-gradient terms that enforce the structural similarity between parameter models. Simulation of acoustic seismic and scalar electromagnetic problems is implemented using 2D finite-difference frequency-domain methods, and the inverse problems of seismic FWI and GPR FWI are solved using a matrix-free truncated Newton algorithm. The joint inversion procedure is performed in several hierarchical frequencies, and the three parameter models are sequentially inverted at each frequency. The joint FWI approach is illustrated using three numerical examples. The results indicate that the joint FWI approach can effectively enhance the structural similarity among the models, modify the structure of each model, and improve the accuracy of inversion results compared with those of individual FWI approaches. Moreover, joint inversion can reduce the trade-off between permittivity and conductivity in GPR FWI, leading to an improved conductivity model in which artifacts are significantly decreased.

Sign in / Sign up

Export Citation Format

Share Document