2.5-D modeling of cross-hole electromagnetic measurement by finite element method

Elsevier BV - Tập 5 - Trang 126-134 - 2008
Jinsong Shen1,2, Wenbo Sun1
1Department of Geophysics, School of Resources and Information Technology, China University of Petroleum, Beijing, China
2CNPC Key Lab of Geophysical Exploration, China University of Petroleum, Beijing, China

Tóm tắt

A finite element method is developed for simulating frequency domain electromagnetic responses due to a dipole source in the 2-D conductive structures. Computing costs are considerably minimized by reducing the full three-dimensional problem to a series of two-dimensional problems. This is accomplished by transforming the problem into y-wave number (K y ) domain using Fourier transform and the y-axis is parallel to the structural strike. In the K y domain, two coupled partial differential equations for magnetic field Hy and electric field Ey are derived. For a specific value of K y , the coupled equations are solved by the finite element method with isoparametric elements in the x-z plane. Application of the inverse Fourier transform to the K y domain provides the electric and magnetic fields in real space. The equations derived can be applied to general complex two-dimensional structures containing either electric or magnetic dipole source in any direction. In the modeling of the electromagnetic measurement, we adopted a pseudo-delta function to distribute the dipole source current and circumvent the problem of singularity at the source point. Moreover, the suggested method used isoparametric finite elements to accommodate the complex subsurface formation. For the large scale linear system derived from the discretization of the Maxwell’s equations, several iterative solvers were used and compared to select the optimal one. A quantitative test of accuracy was presented which compared the finite element results with analytic solutions for a dipole source in homogeneous space for different ranges and different wave numbers K y . to validate the code and check its effectiveness. In addition, we addressed the effects of the distribution range τ of the pseudo-delta function on the numerical results in homogeneous medium.

Tài liệu tham khảo

Abubakar A., Habashy T., Druskin V., et al. Two-and-half dimensional Forward and Inverse Modeling for Marine CSEM Problems. 76th Ann. Internat. Mtg. Soc. Expl. Geophys. Expanded Abstracts. 2006: 750–754 Avdeeva A., Commer M., and Newman, G. A. Hydrocarbon reservoir detectability study for marine CSEM methods: time domain versus frequency domain, 77th Ann. Internat. Mtg. Soc. Expl. Geophys., Expanded Abstracts. 2007: 628–632 Buell J C. A hybrid numerical method for three-dimensional spatially developing free shear flows. Journal of Computational Physics. 1991. 95(2): 313–338 Constable K S. and Weiss C. Mapping 3D salt using 2D marine MT: Case study from Gemini Prospect, Gulf of Mexico: 74th Annual International Mtg. Soc. Expl. Geophys., Expanded Abstracts, 2004: 596–599 Constable, S C. Marine electromagnetic methods—Anew tool for offshore exploration: The Leading Edge. 2006. 25: 438–444 Constable S C and Srnka L J. An introduction to marine controlledsource electromagnetic methods for hydrocarbon exploration. Geophysics. 2007.72(2): WA3–WA12 David L A and Morrison H F. Theoretical and practical considerations for cross well electromagnetic tomography assuming a cylindrical geometry. Geophysics. 1995. 60(3): 846–870 Davydycheva S., Druskin V. and Habashy T. An efficient finitedifference scheme for electromagnetic logging in 3D anisotropic inhomogeneous media. Geophysics. 2003. 68(5): 1525–1536 Druskin V and Knizhnerman L. Spectral approach to solving threedimensional Maxwell’s diffusion equations in the time and frequency domains. Radio Science. 1994. 29(4): 937–953 Edwards N. Marine controlled source electromagnetics: Principles and methodologies-future commercial applications: Surveys in Geophysics. 2005. 26(2): 675–700 Ellingsrud S., Eidesmo T., Johansen S., et al. Remote sensing of hydrocarbon layers by seabed logging (SBL): Results from a cruise offshore Angola. The Leading Edge. 2002. 7: 972–982 Everett M E and Edwards R N. Transient marine electromagnetics: the 2.5-D forward problem. Geophysical Journal International. 1993. 113(3): 545–561 Farrelly B., Ringstad C., Johnstad S E., et al. Remote characterization of hydrocarbon filled reservoirs at the Troll Field by Sea Bed Logging. EAGE Fall Research Workshop. Rhodes. Greece 19th–23rd September. 2004: 662–668 Furumura T and Takenaka H. 2.5-D modeling of elastic waves using the pseudo-spectral method. Geophysical Journal International. 1996. 124: 820–832 Gallardo L A. and Meju M A. Characterization of heterogeneous nearsurface materials by joint 2D inversion of DC resistivity and seismic data. Geophysical Research Letters. 2003. 30: 1658–1670 Goldman Y, Hubans C, Nicoletis S, et al. A finite element solution for the transient electromagnetic response of an arbitrary two-dimensional resistivity distribution. Geophysics. 1986. 51: 1450–1461 Herrmann R B. SH-wave generation by dislocation source-A numerical study. Bulletin of the Seismological Society of America. 1979. 69(1): 1–15 Hou J. and Torres-Verdín C. Finite-difference modeling of EM fields using coupled potentials in 3D anisotropic media: Application to borehole logging: 73rd Annual International Mtg., Soc. Expl. Geophys., Expanded Abstracts. 2003: 522–525 Hoversten G M., Gritto R., Washbourne J. et al. Pressure and fluid saturation prediction in a multi-component reservoir, using combined seismic and electromagnetic imaging. Geophysics. 2003. 68(5): 1580–1591 Hoversten G M. and Newman G A. 3D numerical simulation of a deepwater EM exploration survey. 74th Annual International Mtg., Soc. Expl. Geophys., Expanded Abstracts. 2004: 672–675 Hoversten G M., Newman G. A., Geier N., et al. 3D modeling of a deepwater EM exploration survey. Geophysics. 2006. 71(5): G239–G248 Johansen S E., Amundsen H E., Rosten T., et al. Subsurface hydrocarbons detected by electromagnetic sounding. First Break. 2005. 23: 31–36 Kong F N., Johnstad S E., Røsten T. et al. A 2.5D finite-element modeling difference method for marine CSEM modeling in stratified anisotropic media. 2008. Geophysics. 73(1): F9–F19 Leppin M. Electromagnetic modeling of 3-D sources over 2-D inhomogeneities in the time domain. Geophysics. 1992. 57(8): 994–1003 Liu S X and Sato M. Application of borehole radar for subsurface physical measurement. Journal of Geophysics and Engineering. 2004. 1(3): 221–227 MacGregor L M., Andreis D., Tomlinson J., et al. Controlled-source electromagnetic imaging on the Nuggets-1 reservoir. The Leading Edge. 2006. 25: 984–992 Smith J T. Conservative modeling of 3-D electromagnetic fields, Part2: Biconjugate gradient solution and an accelerator. Geophysics. 1996. 61(5): 1319–1324 Spies B R and Frischknecht F C. Electromagnetic sounding. Electromagnetic methods in applied geophysics: Application. in Nabighian, M. N., Ed., Electromagnetic methods in applied geophysics-theory. 1991. Vol.2: Soc. Explor. Geophys: 431–443 Spies B R and Habashy T M. Sensitivity analysis of cross-well electromagnetics. Geophysics. 1995. 60(3): 834–845 Srnka L J., Carrazone J J., Eriksen E A. et al. Remote reservoir resistivity mapping — An overview. The Leading Edge. 2006. 25: 972–975 Stoyer C H. and Greenfield R J. Numerical solutions of the response of a two-dimensional earth to an oscillating magnetic dipole source, Geophysics. 1976. 41(2): 519–530 Sugeng F, Raiche A, and Xiong Z. Modeling topographic effects using an isoparametric finite element method. The 13th Workshop on Electromagnetic Induction in the Earth. Abstracts. 1996: 109–110 Tompkins M J. Marine controlled source electromagnetic imaging for hydrocarbon exploration: Interpreting subsurface electrical properties. First Break. 2004. 22: 45–51 Ueda T and Zhdanov M S. Fast numerical modeling of marine controlled source electromagnetic data using quasi-liner method. 75th Annual International Mtg., Soc. Expl. Geophys., Expanded Abstracts. 2005: 231–235 Unsworth M J, Travis B J, and Chave A D. Electromagnetic induction by a finite electric dipole source over a 2-D earth. Geophysics. 1993. 58(2): 198–214 Ward S H. and Hohmann G W. Electromagnetic theory for geophysical applications. in Nabighian M. N., Ed., Electromagnetic methods in applied geophysics-theory. Soc. Explor. Geophys. 1988. 1: 131–311 Weiss C J. and Newman G A. Electromagnetic induction in a generalized 3D anisotropic earth, part 2: The LIN preconditioner. Geophysics, 2003. 68(3): 922–930 Zaslavsky M, Davdycheva S., Druskin V. et al. Finite difference solution of the three-dimensional electromagnetic problem using divergence free preconditioners. 76th Ann. Internet. Mtg. Soc. Expl. Geophys., Expanded Abstracts. 2006: 775–778 Zhdanov M.S., Gribenco A. and Cuma M. Regularized focusing inversion of marine CSEM data using minimum vertical support stabilizer, 77th Ann. Internet. Mtg. Soc. Expl. Geophys., Expanded Abstracts. 2007: 579–583 Zienkiewicz O C. The finite element method: third edition. McGraw-Hill. 1977