Transient-electromagnetic finite-difference time-domain earth modeling over steel infrastructure

Including highly conductive steel infrastructure into electromagnetic (EM) earth modeling is motivated by the fact that long metal-cased boreholes have the potential to be used as boosting antennas that enable larger source dipole moments and greater signal penetration depths. Unfortunately, geophysical algorithms designed to simulate EM responses over rather regional scales are complicated by material property contrasts and structure geometries that are more typical for EM engineering applications. Hence, the great majority of earth-modeling algorithms that consider EM responses from steel-cased boreholes use integral-equation methods. To be able to model complex casing scenarios, we revisited the finite-difference time-domain (FDTD) method to advance the modeling of transient-EM field responses from steel-cased boreholes. A time-dependent function that allows for larger FDTD time steps in the DuFort-Frankel method was developed, alleviating the generally large computational overhead. We compared our method against three different kinds of benchmark solutions to demonstrate the reliability of the FDTD field solutions. These test cases were carried out to check the feasibility of a final hydraulic fracturing study. Images of the electric current distribution in a sheetlike rock fracture were calculated for the cases with and without the presence of a connecting borehole casing, demonstrating the casing’s potential of illuminating deep target zones.

[1]  D. Rucker,et al.  Electrical-resistivity characterization of an industrial site using long electrodes , 2010 .

[2]  E. Haber,et al.  Hydro-frac Monitoring Using Ground Time-domain EM , 2014 .

[3]  G. W. Hohmann,et al.  A finite-difference, time-domain solution for three-dimensional electromagnetic modeling , 1993 .

[4]  Michael Commer,et al.  An accelerated time domain finite difference simulation scheme for three‐dimensional transient electromagnetic modeling using geometric multigrid concepts , 2006 .

[5]  Michael S. Zhdanov,et al.  Methods for modelling electromagnetic fields Results from COMMEMI—the international project on the comparison of modelling methods for electromagnetic induction , 1997 .

[6]  J. T. Smith,et al.  Resolution depths for some transmitter-receiver configurations , 2002, IEEE Transactions on Geoscience and Remote Sensing.

[7]  S. H. Ward,et al.  Electromagnetic scattering from cylinders of arbitrary cross-section in a conductive half-space , 1971 .

[9]  David L. Alumbaugh,et al.  An iterative finite element time-domain method for simulating three-dimensional electromagnetic diffusion in earth , 2012 .

[10]  James R. Wait,et al.  EM AND IP RESPONSE OF A STEEL WELL CASING FOR A FOUR-ELECTRODE SURFACE ARRAY. PART I: THEORY* , 1985 .

[11]  Gerald W. Hohmann,et al.  ELECTROMAGNETIC SCATTERING BY CONDUCTORS IN THE EARTH NEAR A LINE SOURCE OF CURRENT , 1971 .

[12]  Gregory A. Newman,et al.  Transient electromagnetic response of a three-dimensional body in a layered earth , 1986 .

[13]  A. Binley,et al.  Detection of leaks in underground storage tanks using electrical resistance methods. , 1996 .

[14]  W Hohmann Gerald,et al.  Numerical modeling for electromagnetic methods of geophysics , 1987 .

[15]  Tarek M. Habashy,et al.  Influence of steel casings on electromagnetic signals , 1994 .

[16]  David L. Alumbaugh,et al.  3D time-domain simulation of electromagnetic diffusion phenomena: A finite-element electric-field approach , 2010 .

[17]  C. Torres‐Verdín,et al.  1D subsurface electromagnetic fields excited by energized steel casing , 2009 .

[18]  Martin D. Müller,et al.  Understanding LOTEM data from mountainous terrain , 2000 .

[19]  H. F. Morrison,et al.  A theoretical study of surface-to-borehole electromagnetic logging in cased holes , 1989 .

[20]  Gerald W. Hohmann,et al.  Diffusion of electromagnetic fields into a two-dimensional earth; a finite-difference approach , 1984 .

[21]  J. Torquil Smith,et al.  Approximating spheroid inductive responses using spheres , 2006 .

[22]  Michael Commer,et al.  A parallel finite-difference approach for 3D transient electromagnetic modeling with galvanic sources , 2004 .

[23]  James R. Wait,et al.  EM AND IP RESPONSE OF A STEEL WELL CASING FOR A FOUR‐ELECTRODE SURFACE ARRAY. PART II: NUMERICAL RESULTS* , 1985 .

[24]  Gordon F. West,et al.  Effect of well casings on surface electrical surveys , 1984 .

[25]  Allen Q. Howard,et al.  THE ELECTROMAGNETIC FIELDS OF A SUBTERRANEAN CYLINDRICAL INHOMOGENEITY EXCITED BY A LINE SOURCE , 1972 .

[26]  James R. Wait The Effect of a Buried Conductor on the Subsurface Fields for Line Source Excitation , 1972 .

[27]  E. Haber,et al.  Hydro‐frac monitoring using ground time‐domain electromagnetics , 2015 .

[28]  H. F. Morrison,et al.  EFFECTS OF WELL CASING ON POTENTIAL FIELD MEASUREMENTS USING DOWNHOLE CURRENT SOURCES1 , 1990 .

[29]  James W. Haslett,et al.  Resistivity Response of a Homogeneous Earth with a Finite-Length Contained Vertical Conductor , 1987, IEEE Transactions on Geoscience and Remote Sensing.

[30]  James Macnae,et al.  6. Time Domain Electromagnetic Prospecting Methods , 1991 .

[31]  J. T. Smith,et al.  UXO detection and identification based on intrinsic target polarizabilities — A case history , 2009 .

[32]  Michael Commer,et al.  Finite element modelling of transient electromagnetic fields near steel-cased wells , 2015 .

[33]  F. Frischknecht,et al.  6. Electromagnetic Physical Scale Modeling , 1988 .

[34]  J. Rocroi,et al.  The use of vertical line sources in electrical prospecting for hydrocarbon , 1985 .

[35]  Douglas W. Oldenburg,et al.  Comparison of integral equation and physical scale modeling of the electromagnetic responses of models with large conductivity contrasts , 2006 .

[36]  Marion Jegen,et al.  The marine controlled source electromagnetic response of a steel borehole casing: applications for the NEPTUNE Canada gas hydrate observatory , 2013 .

[37]  B. Heincke,et al.  GPU parallelization of a three dimensional marine CSEM code , 2013, Comput. Geosci..

[38]  James R. Wait ELECTROMAGNETIC FIELDS OF CURRENT-CARRYING WIRES IN A CONDUCTING MEDIUM , 1952 .

[39]  Bingzhu Wang,et al.  2D and 3D potential‐field upward continuation using splines , 2006 .

[40]  Gerald W. Hohmann,et al.  Why time domain , 1980 .

[41]  H. Kim,et al.  Sensitivity analysis of marine controlled-source electromagnetic methods to a shallow gas-hydrate layer with 1D forward modeling , 2011 .

[42]  Robin L. Newmark,et al.  Low-cost reservoir tomographs of electrical resistivity , 2004 .

[43]  David L. Alumbaugh,et al.  Numerical modeling analysis of short-offset electric-field measurements with a vertical electric dipole source in complex offshore environments , 2012 .

[44]  B. M. A. Rahman,et al.  Elimination of Numerical Dispersion from Electromagnetic Time Domain Analysis by Using Resource Efficient Finite Element Technique , 2013 .

[45]  Yingkang Wei,et al.  Propagation of Electromagnetic Signal along a Metal Well in an Inhomogeneous Medium , 2013 .