Numerical modeling of borehole-surface electromagnetic responses with 3-D finite difference method and comparison with physical simulations

详细信息    查看全文

• 作者：Hui Cao (1)
  Lifeng Mao (1)
  Xuben Wang (1)
  Zhanxiang He (2)
  Kunpeng Wang (3)
  
  1. Key Laboratory of Earth Exploration and Information Technology of MOE ; Chengdu University of Technology ; Chengdu ; Sichuan ; 610059 ; China
  2. BGP INC. ; China National Petroleum Corporation ; Zhuozhou ; Hebei ; 072750 ; China
  3. School of Geophysics and Information Technology ; China University of Geoscience ; Beijing ; 100083 ; China
  
• 关键词：borehole ; surface electromagnetic method ; finite difference method ; numerical simulation ; physical simulation
• 刊名：Studia Geophysica et Geodaetica
• 出版年：2015
• 出版时间：January 2015
• 年：2015
• 卷：59
• 期：1
• 页码：83-96
• 全文大小：699 KB
• 参考文献：1. An R., Li T. and Xu K., 2007. Well-surface 3-D resistivity inversion. / Progr. Geophys., 22 247鈥?49.
  2. Badea E.A., Everett M.E., Newman G.A. and Biro O., 2001. Finite-element analysis of controlledsource electromagnetic induction using Coulomb-gauged potentials. / Geophysics, 66, 786鈥?99. CrossRef
  3. Cao H., Wang X. and He Z., 2012. Borehole-to-surface electromagnetic responses of the horizontally stratified medium. / Oil Geophys. Prospect., 47, 338鈥?43.
  4. Champagne N.J., Berryman J.G. and Buettner H.M., 2001. FDFD: A 3D finite-difference frequency-domain code for electromagnetic induction tomography. / J. Comput. Phys., 170, 830鈥?48. CrossRef
  5. Freund R., 1993. A transpose-free quasi-minimal residual algorithm for non-Hermitian linear systems. / SIAM J. Sci. Comput, 14, 470鈥?82. CrossRef
  6. Frischknecht F.C., 1988. Electromagnetic physical scale modeling. In: Nabighian M.S. (Ed.), / Electromagnetic Methods in Applied Geophysics, Volume I, Theory. Society of Exploration Geophysicists, Tulsa, OK, 365鈥?41.
  7. Guptasarma D. and Singh B., 1997. New digital linear filters for Hankel J0 and J1 transforms. / Geophys. Prospect., 45, 745鈥?62. CrossRef
  8. Hohmann G.W., 1975. Three-dimensional induced polarization and EM modeling. / Geophysics, 40, 309鈥?24. CrossRef
  9. Ke G. and Huang Q., 2009. 3D forward and inversion problems of borehole-to-surface electrical method. / Acta Scientiarum Naturalium Universitatis Pekinensis, 45, 264鈥?72.
  10. Key K., 2009. 1D Inversion of multi-component, multi-frequency marine CSEM data, methodology and synthetic studies for resolving thin resistive layers. / Geophysics, 74, F9鈥揊20. CrossRef
  11. Liu C., Ren Z., Tang J. and Yan L., 2008. Three-dimensional magnetotellurics modeling using edgebased finite-element unstructured meshes. / Appl. Geophys., 5, 170鈥?80. CrossRef
  12. Liu X., Wang J., He Z. and Wang Z., 2006. Study of hydrocarbon accumulation by borehole-ground EM method. / Progr. Explor. Geophys., 29, 98鈥?01.
  13. Mitsuhata Y. and Uchida T., 2004. 3D magnetotelluric modeling using the T鈥?finite-element method. / Geophysics, 69, 108鈥?19. CrossRef
  14. Newman G.A. and Alumbaugh D.L., 1995. Frequency-domain modeling of airborne electromagnetic responses using staggered finite differences. / Geophys. Prospect., 43, 1021鈥?041. CrossRef
  15. Newman G.A., Hohmann G.W. and Anderson W.L., 1986. Transient electromagnetic response of the three-dimensional body in a layered Earth. / Geophysics, 51, 1608鈥?627. CrossRef
  16. Sleijpen G.L.G. and Fokkema D.R., 1993. BICGSTAB(L) for linear equations involving unsymmetric matrices with complex spectrum. / Electron. Trans. Numer. Anal., 1, 11鈥?2.
  17. Streich R., 2009. 3D Finite-difference frequency-domain modeling of controlled-source electromagnetic data. Direct solution and optimization for high accuracy. / Geophysics, 74, 95鈥?05. CrossRef
  18. Tripp A.C. and Hohmann G.W., 1984. Block diagonalization of the electromagnetic impedance matrix of a symmetric buried body using group theory. / IEEE Trans. Geosci. Remote Sens., 22, 62鈥?9. CrossRef
  19. Wang Z., He Z. and Liu Y., 2006. Research of 3D modeling and anomalous rule on boreholeground DC method. / Chinese J. Eng. Geophys., 3, 87鈥?2.
  20. Wang Z., He Z. and Wei W., 2007a. 3D modeling and Born approximation inversion for the borehole surface electromagnetic method. / Appl. Geophys., 4, 84鈥?8. CrossRef
  21. Wang Z., He Z. and Wei W., 2007b. Parallel algorithm research of integral equation method to 3D borehole-surface electromagnetic modeling. / Comput. Techn. Geophys. Geochem. Explor., 29, 425鈥?30.
  22. Wang Z., He Z., Wei W. and Deng M., 2005. 3D physical model experiments of well-to-ground electrical survey. / Oil Geophys. Prospect., 40, 594鈥?97.
  23. Xiong Z., 1992. Electromagnetic modeling of 3D structures by the method of system iteration using integral equations. / Geophysics, 57, 1556鈥?561. CrossRef
  24. Xu K. and Li T., 2006. The forward modeling of three-dimensional geoelectric field of vertical finite line source by finite-difference method. / J. Jilin Univ. (Earth Sci.), 36, 137鈥?41.
  25. Zhdanov M.S. and Fang S., 1997. Quasi-linear series in three-dimensional electromagnetic modeling. / Radio Sci., 31, 741鈥?54. CrossRef
  26. Zhang T. and Zhang B., 1995. Research on removing noises in DC resitivity trielectrode gradient measurement. / Oil Geophys. Prospect., 30, 100鈥?10.
  27. Yee K., 1966. Numerical solution of initial boundary value problems involving Maxwell鈥檚 equations in isotropic media. / IEEE Trans. Antennas Propag., 16, 302鈥?07.
  
• 刊物类别：Earth and Environmental Science
• 刊物主题：Earth sciences
  Geophysics and Geodesy
  Structural Geology
  Meteorology and Climatology
  
• 出版者：Springer Netherlands
• ISSN：1573-1626

文摘

Borehole-surface electromagnetic method (BSEM) is a high-accuracy electromagnetic prospecting method that uses AC-powered vertical finite line source as the excitation source. The observed surface electromagnetic field is inverted for the subsurface conductivity structure. In this paper, the total field is separated into a primary part, due to a horizontally layered host medium, and a secondary part due to 3-D heterogeneities. After solving for the primary field with an analytical method, the frequency-domain second-order differential equation for the secondary field is discretized by finite differences. A sparse matrix storage scheme is employed and a BiCGSSTAB(m) method with a diagonal matrix preconditioner is used to obtain the secondary field as well as the 3-D BSEM response of the model. We compare the result of anomalous responses of three-layered medium derived by 3-D forward modeling with the result of semi-analytical solution. We also perform physical simulation and 3-D numerical forward modeling based on similarity criterion. As a result, the shape of both anomalous response curves are the same, which validates the 3-D numerical simulation method. The anomalous fields of 3-D numerical forward and physical simulation share similar anomalous feature of symmetric bimodal structure that is consistent with its harmonic response curve. It demonstrates that borehole-surface electromagnetic method can be used not only for prospection by employing a multi-frequency response, but also can provide multi-angle information about subsurface anomaly by varying relative depths of vertical finite line source in the borehole.