Constraint inversion of three-dimensional borehole-to-surface resistivity based on unstructured finite element
WANG Zhi1(), WANG Cheng2, FANG Si-Nan3
1. Electronics and Information School, Yangtze University, Jingzhou 434023,China 2. Xi’an Research Institute Co. Ltd., China Coal Technology and Engineering Group Corp., Xi’an 710077, China 3. College of Geophysics and Petroleum Resources,Yangtze University,Wuhan 430100,China
The inversion of electromagnetic detection data is a typical ill-posed problem and is prone to cause a multiplicity of solutions of the inversion results. The ill-posedness is an inherent characteristic of inversion and is difficult to overcome without additional information. An effective way to solve this problem is constrained inversion. In this study, the Gauss-Newton - conjugate gradient (GN-CG) method was used to directly impose constraints on the inversion objective function. Specifically, the dielectric resistivity range was introduced into the inversion objective function as the prior information and constraints using the exterior penalty function method. Compared with the conventional three-dimensional resistivity inversion objective function, the objective function with inequality constraints can suppress the multiplicity of solutions in theory. As revealed by the testing results of various theoretical models, the three-dimensional borehole-to-surface resistivity inversion algorithm based on inequality constraints effectively improves the precision of inversion results, and the way of imposing inequality constraints using the penalty function method is feasible and effective.
Liu J X, Zhao R, Guo Z W. Research progress of electromagnetic methods in the exploration of metal deposits[J]. Progress in Geophysics, 2019, 34(1):151-160.
Lyu Q T, Zhang X P, Tang J T, et al. Review on advancement in technology and equipment of geophysical exploration for metallic deposits in China[J]. Chinese Journal of Geophysics, 2019, 62(10):3629-3664.
Wu X P, Liu Y, Wang W. 3D resistivity inversion incorporating topography based on unstructured meshes[J]. Chinese Journal of Geophysics, 2015, 58(8):2706-2717.
Wang Z, Wu A P, Li G. Forward modeling of borehole-ground induced polarization method under undulating topography[J]. Geophyscial Prospecting for Petroleum, 2018, 57(6):927-935,951.
Pan H P. The Role of Borehole induced polarization/resistivity method in the exploration of mineral resoueces[J]. Geophysics and Geochemical Exploration, 2013, 37(4):620-626.
Tang J T, Zhang J F, Feng B, et al. Detemination of borders for resistive oil and gas reservoirs by deviation rate using the hole-to-surface resistivity method[J]. Chinese Journal of Geophysics, 2007, 50(3):926-931.
Huang J G, Wang J L, Ruan B Y. A study on FEM modeling of anomalies of 3-D high-density E-SCAN resistivity survey[J]. Chinese Journal of Geophysics, 2006, 49(4):1206-1214.
[8]
Li Y, Spitzer K. Finite element resistivity modelling for three-dimensional structures with arbitrary anisotropy[J]. Physics of the Earth & Planetary Interiors, 2005, 150(1-3):15-27.
[9]
Wu X. A 3-D finite-element algorithm for DC resistivity modelling using the shifted incomplete Cholesky conjugate gradient method[J]. Geophysical Journal International, 2003, 154(3):947-956.
Wu X P, Wang T T. A 3-D finite-element resistivity forward modeling using conjugate gradient algorithm[J]. Chinese Journal of Geophysics, 2003, 46(3):428-432.
[11]
Li Y, Klaus S. Three-dimensional DC resistivity forward modelling using finite elements in comparison with finite-difference solutions[J]. Geophysical Journal International, 2002, 51(3):924-934.
Huang J G, Ruan B Y, Bao G S. Fem under quantic-boundary condition for modeling resistivity on 3-D geoelectric section[J]. Journal of Guilin Institute of Technology, 2002, 22(1):11-14.
[13]
Ren Z, Qiu L, Tang J, et al. 3-D direct current resistivity anisotropic modelling by goal-oriented adaptive finite element methods[J]. Geophysical Journal International, 2018, 212(1):76-87.
[14]
Ren Z, Tang J. A goal-oriented adaptive finite-element approach for multi-electrode resistivity system[J]. Geophysical Journal International, 2014, 199(1):136-145.
[15]
Wei W, Xiaoping W, Spitzer K. Three-dimensional DC anisotropic resistivity modelling using finite elements on unstructured grids[J]. Geophysical Journal International, 2013, 193(2):734-746.
[16]
Ren Z, Jingtian T. 3D direct current resistivity modeling with unstructured mesh by adaptive finite-element method[J]. Geophysics, 2010, 75(1):H7-H17.
[17]
Blome M, Maurer H R, Schmidt K. Advances in three-dimensional geoelectric forward solver techniques[J]. Geophysical Journal International, 2009, 176(3):740-752.
[18]
Rücker C, Günther T, Spitzer K. Three-dimensional modelling and inversion of dc resistivity data incorporating topography-I. Modelling[J]. Geophysical Journal International, 2006, 166(2):495-505.
[19]
Günther T, Rücker C, Spitzer K. Three-dimensional modelling and inversion of DC resistivity data incorporating topography-II. Inversion[J]. Geophysical Journal of the Royal Astronomical Society, 2006, 166(2):506-517.
[20]
Zhou B, Greenhalgh S A. Finite element three-dimensional direct current resistivity modelling; accuracy and efficiency considerations[J]. Geophysical Journal International, 2001, 145(3):679-688.
[21]
Sasaki Y. 3-D resistivity inversion using finite-element method[J]. Geophysics, 1994, 59(12):1839.
[22]
Wu X, Xiao Y, Qi C, et al. Computations of secondary potential for 3-D DC resistivity modelling using an incomplete Choleski conjugate-gradient method[J]. Geophysical Prospecting, 2003, 51(6):567-577.
Wu X P, Xu G M, Li S C. The calculation of three-dimensional geoelectric field of point source by incomplete cholesky conjugate gradient method[J]. Chinese Journal of Geophysics, 1998, 41(6):848-855.
[25]
Zhao S, Yedlin M J. Some refinements on the finite-difference method for 3-D dc resistivity modeling[J]. Geophysics, 1996, 61(5):1301-1307.
[26]
Zhang J, Mackie R L, Madden T R. 3-D resistivity forward modeling and inversion using conjugate gradients[J]. Geophysics, 1995, 60(5):1313-1325.
[27]
Spitzer K. A 3-D finite-difference algorithm for DC resistivity modelling using conjugate gradient methods[J]. Geophysical Journal International, 1995, 123(3):903-914.
Tang J T, Wang F Y, Ren Z Y. 2.5-D DC resistivity modeling by adaptive finite-element method with unstructured triangulation[J]. Chinese Journal of Geophysics, 2010, 53(3):708-716.
Ren Z Y, Tang J T. Finite element modeling of 3-D DC resistivity using locally refined unstructured meshes[J]. Chinese Journal of Geophysics, 2009, 52(10):2627-2634.
Peng R H, Hu X Y, Li J H, et al. 3D inversion of frequency-domain marine CSEM data in VTI media[J]. Chinese Journal of Geophysics, 2019, 62(6):2165-2175.
[31]
Gundogdu N Y, Candansayar M E. Three-dimensional regularized inversion of DC resistivity data with different stabilizing functionals[J]. Geophysics, 2018, 83(6):E399-E407.
Guo L G, Dai G L, Yang B C. et al. 3D resistivity inversion with multiple priori-information constraint[J]. Oil Geophysical Prospecting, 2018, 53(6):1333-1340.
Peng R H, Hu X Y, Han B. 3D inversion of frequency-domain CSEM data based on Gauss-Newton optimization[J]. Chinese Journal of Geophysics, 2016, 59(9):3470-3481.
[34]
Oldenburg D W, Haber E, Shekhtman R. Three dimensional inversion of multisource time domain electromagnetic data[J]. Geophysics, 2013, 78(1):E47-E57.
[35]
Pidlisecky A, Haber E, Knight R J. RESINVM3D: A 3D resistivity inversion package[J]. Geophysics, 2007, 72(2):H1-H10.
Xu K J, Li T L, Zhang H, et al. 3D resistivity inversion of vertical finite line source using conjugate gradients[J]. Coal Geology and Exploration, 2006, 34(3):68-71.
[38]
Cao X Y, Huang X, Yin C C, et al. 3D MT anisotropic inversion based on unstructured finite-element method[J]. Journal of Environmental & Engineering Geophysics, 2021, 26(1):49-60.
Hui Z J, Yin C C, Liu Y H, et al. 3D inversion of time-domain marine EM data based on unstructured finite-element method[J]. Chinese Journal of Geophysics, 2020, 63(8):3167-3179.
Yu H, Deng J Z, Chen H, et al. Three-dimensional magnetotelluric inversion under topographic relief based on the limited-memory quasi-Newton algorithm(L_BFGS)[J]. Chinese Journal of Geophysics, 2019, 62(8):3175-3188.
Deng Y, Tang J, Ruan S. Adaptive regularized three-dimensional magnetotelluric inversion based on the LBFGS quasi-Newton method[J]. Chinese Journal of Geophysics, 2019, 62(9):3601-3614.
Yin C C, Zhu J, Qiu C K, et al. Spatially constrained inversion for airborne EM data using quasi-3D models[J]. Chinese Journal of Geophysics, 2018, 61(6):2537-2547.
Qin C, Wang X B, Zhao N. Parallel three-dimensional forward modeling and inversion of magnetotelluric based on a secondary field approach[J]. Chinese Journal of Geophysics, 2017, 60(6):2456-2468.
Liu Y H, Yin C C. 3D inversion for frequency-domain HEM data[J]. Chinese Journal of Geophysics, 2013, 56(12):4278-4287.
[46]
Avdeev D, Avdeeva A. 3D Magnetotelluric inversion using a limited-memory quasi-Newton optimization[J]. Geophysics, 2009, 74(3):F45.
[47]
Xiao Y, Wei Z, Wang Z. A limited memory BFGS-type method for large-scale unconstrained optimization[J]. Computers & Mathematics with Applications, 2008, 56(4):1001-1009.
[48]
Eldad, Haber, Douglas, et al. Inversion of time domain three-dimensional electromagnetic data[J]. Geophysical Journal International, 2007, 171(2):B23-B34.
[49]
Avdeeva A, Avdeev D B. A limited memory quasi-Newton inversion for 1D magnetotellurics[J]. Geophysics, 2006, 71(5):G191-G196.
[50]
Haber E. Quasi-Newton methods for large-scale electromagnetic inverse problems[J]. Inverse Problems, 2005, 21(1):305-323.
[51]
Newman G A, Boggs P T. Solution accelerators for large-scale 3D electromagnetic inverse problems[J]. Inverse Problems, 2004, 20(6):S151-S170.
[52]
Nash S G, Nocedal J. A numerical study of the limited memory BFGS method and the truncated-Newton method for large scale optimization[J]. SIAM Journal on Optimization, 1991, 1(3):358-372.
[53]
Nocedal J. Updating quasi-Newton matrices with limited storage[J]. Mathematics of Computation, 1980, 35(151):773-782.
Tang C Z, Cheng J Z, Yan L J, et al. LBFGS CSAMT 1D inversion of limited memory based on boundary constraint and its application[J]. Coal Geology and Exploration, 2019, 47(5):193-200.
Ma H, Guo Y, Wu P P, et al. 3-D joint inversion of multi-array data set in the resistivity method based on MPI parallel algorithm[J]. Chinese Journal of Geophysics, 2018, 61(12):5052-5065.
Dong H, Wei W B, Ye G F, et al. Study of Three-dimensional magnetotelluric inversion including surface topography based on Finite-difference method[J]. Chinese Journal of Geophysics, 2014, 57(3):939-952.
Zhang K, Dong H, Yan J Y, et al. A NLCG inversion method of magnetotellurics with parallel structure[J]. Chinese Journal of Geophysics, 2013, 56(11):3922-3931.
Lin C H, Tan H D, Shu Q, et al. Three-dimensional conjugate gradient inversion of CSAMT data[J]. Chinese Journal of Geophysics, 2012, 55(11):3829-3838.
Hu Z Z, Hu X Y, He Z X. Pseudo-three-dimensional magnetotelluric using nonlinear conjugate gradients[J]. Chinese Journal of Geophysics, 2006, 49(4):1226-1234.
[60]
Rodi W L, Mackie R L. Nonlinear conjugate gradients algorithm for 2-D magnetotelluric inversion[J]. Geophysics, 2001, 66(1):174-187.
[61]
Newman G A, Alumbaugh D L. Three-dimensional magnetotelluric inversion using non-linear conjugate gradients[J]. Geophysical Journal International, 2000, 140(2):410-424.
Wang Z, Pan H P, Luo Y H, et al. 3-D hole-to-surface resistivity inversion with nonlinear conjugate gradients method under the constraint of inequality[J]. Progress in Geophysics, 2016, 31(1):360-370.
Wang Z, Pan H P, Wu A P, et al. 3D inversion of borehole induced polarization under the inequality constraint[J]. Geophysical Prospecting for Petroleum, 2016, 55(3):455-466.
[64]
Kim H J, Song Y, Lee K H. Inequality constraint in least-squares inversion of geophysical data[J]. Earth, Planets and Space, 1999, 51(4):255-259.
[65]
Li Y, Oldenburg D W. 3-D Inversion of induced polarization data[J]. Geophysics, 2000, 65(6):1931-1945.
Huang J G, Ruan B Y, Bao G S. Resistivity inversion on 3-D section based on FEM[J]. Journal of Central South University:Science and Technology, 2004, 35(2):295-299.
Wan X L, Xi D Y, Gao E G, et al. 3-D resistivity inversion by the least-squares QR factorization method under improved smoothness constraint condition[J]. Chinese Journal of Geophysics, 2005, 48(2):439-444.
Liu B, Li S C, Nie L C, et al. Inversion imaging of 3D resistivity detection using adaptive-weighted smooth constraint and PCG algorithm[J]. Chinese Journal of Geotechnical Engineering, 2012, 34(9):1646-1653.
Liu B, Li S C, Li S C, et al. 3D electrical resistivity inversion with least-squares method based on inequality constraint and its computation effciency optimization[J]. Chinese Journal of Geophysics, 2012, 55(1):260-268.
Liu B, Nie L C, Li S C, et al. 3D electrical resistivity inversion tomography with spatial structural constraint[J]. Chinese Journal of Rock Mechanics and Engineering, 2012, 31(11):2258-2268.
Di Q Y, Xue G Q, Yin C C, et al. New methods of controlled-source electromagnetic detection in China[J]. Science China Earth Sciences, 2020, 50(9):1219-1227.
Han B, Hu X Y, He Z X, et al. Mathematical classification of magnetotelluric inversion methods[J]. Oil Geophysical Prospecting, 2012, 47(1):177-188.
[76]
Siripunvaraporn W. Three-dimensional magnetotelluric inversion: An introductory guide for developers and users[J]. Surveys in Geophysics, 2012, 33(1):5-27.
[77]
Geuzaine C, Remacle J F. Gmsh: A three-dimensional finite element mesh generator with built-in pre-and post-processing facilities[J]. International Journal for Numerical Methods in Engineering, 2009, 79(1):1309-1331.
[78]
Ahrens J, Geveci B, Law C. ParaView: An end-user tool for large data visualization[M]. New York: Academic Press, 2005.
[79]
Ellis R G, Oldenburg D W. The pole-pole 3-D Dc-resistivity inverse problem:A conjugategradient approach[J]. Geophysical Journal of the Royal Astronomical Society, 1994, 119(1):187-194.
[80]
Li Y, Oldenburg D W. Inversion of 3-D DC resisitivity data using an approximate inverse mapping[J]. Geophysical Journal International, 1994, 116(4):527-537.