Smooth inversion of induction logs for conductivity models with mud filtrate invasion

Geophysics ◽  
2000 ◽  
Vol 65 (5) ◽  
pp. 1468-1475 ◽  
Author(s):  
Catherine D. de Groot—Hedlin

A robust, efficient method of inversion of induction logging data for smooth 2-D models, appropriate to an environment in which mud filtrate invades flat‐lying layers, is described. An infinite number of solutions exist to the problem of determining a conductivity structure from a finite number of imprecise induction data. Therefore, the inverse problem is regularized such that the smoothest model is sought subject to the condition that the resulting computed log agrees with the field log to a given preset level. At each iteration, the Jacobian sensitivities are approximated using the distorted Born approximation. In most cases, the algorithm converges in 3 to 4 iterations. The resulting maximally smooth models reflect the resolution power of the induction data and are unlikely to result in overinterpretation of the data. Inversion of both synthetic and field data indicates that layer boundaries are well resolved but radial boundaries are poorly resolved by conventional induction logging data.

2020 ◽  
Vol 222 (2) ◽  
pp. 1074-1089 ◽  
Author(s):  
Yajun Liu ◽  
Pritam Yogeshwar ◽  
Xiangyun Hu ◽  
Ronghua Peng ◽  
Bülent Tezkan ◽  
...  

SUMMARY Electrical anisotropy of formations has been long recognized by field and laboratory evidence. However, most interpretations of long-offset transient electromagnetic (LOTEM) data are based on the assumption of an electrical isotropic earth. Neglecting electrical anisotropy of formations may cause severe misleading interpretations in regions with strong electrical anisotropy. During a large scale LOTEM survey in a former mining area in Eastern Germany, data was acquired over black shale formations. These black shales are expected to produce a pronounced bulk anisotropy. Here, we investigate the effects of electrical anisotropy on LOTEM responses through numerical simulation using a finite-volume time-domain (FVTD) algorithm. On the basis of isotropic models obtained from LOTEM field data, various anisotropic models are developed and analysed. Numerical results demonstrate that the presence of electrical anisotropy has a significant influence on LOTEM responses. Based on the numerical modelling results, an isolated deep conductive anomaly presented in the 2-D isotropic LOTEM electric field data inversion result is identified as a possible artifact introduced by using an isotropic inversion scheme. Trial-and-error forward modelling of the LOTEM electric field data using an anisotropic conductivity model can explain the data and results in a reasonable quantitative data fit. The derived anisotropic 2-D model is consistent with the prior geological information.


Geophysics ◽  
2005 ◽  
Vol 70 (3) ◽  
pp. G51-G58 ◽  
Author(s):  
Tae Jong Lee ◽  
Toshihiro Uchida

Electromagnetic (EM) traveltime tomography has been applied for reservoir characterization at the Lost Hills oil field, California. Four data sets at frequencies of 24, 90, 370, and 1000 Hz were obtained along a pair of monitoring boreholes located 80.5 m apart. Traveltime information was first extracted from these EM data sets using a wavefield transform with a ray series approximation. The conductivity contrast of each layer is no greater than two in the region of interest, so the first arrivals can be estimated within 5% error by the approximate scheme. A nonlinear traveltime tomography algorithm adopting a Fresnel zone concept was then applied to obtain the conductivity model between the boreholes. The resultant conductivity image represents the conductivity structure between the boreholes. This image is consistent with the results of both a finite-difference inversion and the induction log obtained prior to waterflooding. Comparing the two conductivity images with the induction logs, we observe major differences in the fracture-dominant resistive reservoir layer, which may have been caused by changes in reservoir condition before and after waterflooding.


Geophysics ◽  
2012 ◽  
Vol 77 (4) ◽  
pp. WB219-WB231 ◽  
Author(s):  
P. Kaikkonen ◽  
S. P. Sharma ◽  
S. Mittal

Three-dimensional linearized nonlinear electromagnetic inversion is developed for revealing the subsurface conductivity structure using isolated very low frequency (VLF) and VLF-resistivity anomalies due to conductors that may be arbitrarily directed towards the measuring profiles and the VLF transmitter. We described the 3D model using a set of variables in terms of geometric and physical parameters. These model parameters were then optimized (parametric inversion) to obtain their best estimates to fit the observations. Two VLF transmitters, i.e., the [Formula: see text], [Formula: see text] (“E”) and the [Formula: see text], [Formula: see text] (“H”) polarizations, respectively, can be considered jointly in inversion. After inverting several noise-free and noisy synthetic data, the results revealed that the estimated model parameters and the functionality of the approach were very good and reliable. The inversion procedure also worked well for the field data. The reliability and validity of the results after the field data inversion have been checked using data from a shear zone associated with uranium mineralization.


Geophysics ◽  
2000 ◽  
Vol 65 (4) ◽  
pp. 1113-1123 ◽  
Author(s):  
Andreas Hördt ◽  
Martin Müller

Long‐offset transient electromagnetic (LOTEM) data from the Vesuvius volcano, in Italy, show that the EM response of the topography is a potential cause of data distortions. A modeling study was carried out to simulate the effect of mountainous terrain on vertical magnetic‐field time derivatives using a 3-D finite‐difference code. The objectives were to assess the importance of topographic effects and to help identify them in existing field data. The total effect of topography on the LOTEM response can be considered as a combination of four distortions of the corresponding responses for a flat terrain. First, the receiver is at some height above the flat surface. Second, the mountain acts as a conductive body displacing air. Third, large loop receivers are nonhorizontal and sense a combination of horizontal and vertical magnetic fields. Finally, the electromagnetic coupling between the mountain and deeper‐lying structure modifies the structure response. Each of the effects can be identified in field data recorded at Mount Vesuvius. The topographic induced distortions for the model used in this study are moderate in the sense that 1-D inversions of the theoretical data still recover the gross conductivity structure, albeit with small deviations from the true parameters. Although this result might imply that topography might be ignored during the first stage of an interpretation, no simple correction method is evident, so topography will have to be included in any 2-D or 3-D inversion attempt.


2018 ◽  
Vol 2018 (1) ◽  
pp. 1-4
Author(s):  
Zhaohai Meng ◽  
Fengting Li ◽  
Hao Yu ◽  
Lin Ma ◽  
Zhongli Li

Geophysics ◽  
1997 ◽  
Vol 62 (1) ◽  
pp. 253-265 ◽  
Author(s):  
Enru Liu ◽  
Stuart Crampin ◽  
John A. Hudson

We describe a method of modeling seismic waves interacting with single liquid‐filled large cracks based on the Kirchhoff approximation and then apply it to field data in an attempt to estimate the size of a hydraulic fracture. We first present the theory of diffraction of seismic waves by fractures using a Green's function representation and then compute the scattered radiation patterns and synthetic seismograms for fractures with elliptical and rectangular shapes of various dimensions. It is shown that the characteristics of the diffracted wavefield from single cracks are sensitive to both crack size and crack shape. Finally, we compare synthetic waveforms to observed waveforms recorded during a hydraulic fracturing experiment and are able to predict successfully the size of a hydraulically induced fracture (length and height). In contrast to previously published work based on the Born approximation, we model both phases and amplitudes of observed diffracted waves. Our modeling has resulted in an estimation of a crack length 1.1 to 1.5 times larger than previously predicted, whereas the height remains essentially the same as that derived using other techniques. This example demonstrates that it is possible to estimate fracture dimensions by analyzing diffracted waves.


Sign in / Sign up

Export Citation Format

Share Document