1. School of Geophysics and Measurement-control Technology, East China University of Technology, Nanchang, Jiangxi 330013, China;
2. School of Geosciences and Info-Physics, Central South University, Changsha, Hunan 410083, China
Abstract:We use the unstructured adaptive grid algorithm based on model sensitivity information to generate high-quality inversion-oriented meshes,and adopt them to 3D DC resistivity inversion.The model sensitivity is a measure of the overall response of underground model change to real datasets.So,the optimized mesh generated based on the model sensitivity information has higher quality than the ordinary mesh.This new mesh can reduce the dependence of inversion on model regularization constraints and improve the quality of inversion results.We use the minimum-structure inversion object function based on a tetrahedral element mesh.The Gauss-Newton method optimizes the inversion object function.The Gauss-Newton equation is solved by the stable double conjugate gradient method in order to increase the stability of 3D DC resistivity inversion.The calculating results of both synthetic and field data prove the validity of the proposed approach.
吴小平,刘洋,王威.基于非结构网格的电阻率三维带地形反演[J].地球物理学报,2015,58(8):2706-2717.WU Xiaoping,LIU Yang,WANG Wei.3D resistivity inversion incorporating topography based on unstructured meshes[J].Chinese Journal of Geophysics,2015,58(8):2706-2717.
[3]
Zhang J,Mackie R L,Madden T R.Three-dimensional resistivity forward modeling and inversion using conjugate gradients[J].Geophysics,1995,60(5):1313-1325.
[4]
Loke M H,and Barker R D.Rapid least-squares inversion of apparent conductivity pseudosection using a quasi-Newton method[J].Geophysical Prospecting,1996,44(1):131-152.
[5]
刘洋.基于非结构网格的电阻率三维正反演及其应用硏究[D].安徽合肥:中国科学技术大学,2016.LIU Yang.Study of 3D Resistivity Modeling and Inversion Using Unstructured Grids and Their Applications[D].University of Science and Technology of China,Hefei,Anhui,2016.
[6]
Coggon J H.Electromagnetic and electrical modeling by the finite element method[J].Geophysics,1971,36(1):132-155.
[7]
周熙襄,钟本善.电法勘探数值模拟技术[M].四川成都:四川科学技术出版社,1986.
[8]
徐世浙.地球物理中的有限单元法[M].北京:科学出版社,1994.
[9]
Sasaki Y.3D resistivity inversion using the finite-element method[J].Geophysics,1997,59(11):1839-1848.
[10]
Rücker C,Günther T,Spitzer K.Three dimensional modeling and inversion of DC resistivity data incorporating topography I:Modeling[J].Geophysical Journal International,2006,166:495-505.
[11]
任政勇,汤井田.基于局部加密非结构化网格的三维电阻率法有限元数值模拟[J].地球物理学报,2009,52(10):2627-2634.REN Zhengyong,TANG Jingtian.Finite element mo-deling of 3-D resistivity using locally refined unstructured meshes[J].Chinese Journal of Geophysics,2009,52(10):2627-2634.
[12]
Lowry T,Allen M B,and Shive P N.Singularity removal:a refinement of resistivity modeling techniques[J].Geophysics,1989,54(6):766-774.
[13]
Williams N C.Geologically-Constrained Ubc-Gif Gravity and MagneticInversions with Examples from the Agnew-Wiluna Greenstone Belt[D].University of British Columbia,2008.
[14]
Blome M,Maurer H R,Schmidt K.Advances in three-dimensional geoelectric forward solver techniques[J].Geophysical Journal International,2009,176:740-752.
[15]
Menke W.Geophysical Data Analysis:Discrete Inverse:Theory[M].Academic Press,San Diego,1989.
[16]
Bohm G and Vesnaver A.Inquest of the grid[J].Geophysics,1999,64(4):1116-1125.
[17]
Kim M,Kim S,Lee K,et al.Improvement of the electrical impedance tomographic image for the two-phase system with adaptive element grouping technique[J].Measurement Science & Technology,2004,15(7):1391-1401.
[18]
Pessel M and Gibert D.Multiscale electrical impedancetomography[J].Journal of Geophysical Research,2003,108:2054.
[19]
Molinar M,Blott B,Cox S,et al.Optimal imaging with adaptive mesh refinement in electrical impedance tomography[J].Physiological Measurement,2002,23(1):121-128.
[20]
Günther.Three-dimensional modelling and inversion of DC resistivity data incorporating toporating Ⅱ:Inversion[J].Geophysical Journal International,2006,166:506-517.
[21]
Plattner.3-D electrical resistivity tomography using adaptive waveletparameter grids[J].Geophysical Journal International,2012,189:317-330.
[22]
Wang W,Wu X P,Spitzer K.Three-dimensional DC anisotropic resistivity modelling using finite elements on unstructured grids[J].Geophysical Journal International,2013,193:734-746.
[23]
吉洪诺夫·阿尔先宁著,王秉忱译.不适定问题的解法[M].北京:地质出版社,1979.
[24]
Mallows C L.Some comments on Cp[J].Technome-trics,1973,15(4):661-675.
[25]
Golub G H,Michael H,Grace W.Generalized cross-validation as a method for choosing a good Ridge parameter[J].Technometrocs,1979,21(2):215-223.
[26]
Colin G,Farquharson and Douglas W,et al.Automatic estimation of the trade-off parameter in nonlinear inverse problems using the GCV and L-curve criteria[C].SEG Technical Program Expanded Abstracts,2000,19:265-268.
[27]
Hansen P C,and O'Leary D P.The use of the L-curve in the regularization of discrete ill-posed problems[J].SIAM Journal on Numerical Analysis,1993,14(6):1487-1530.
[28]
陈小斌.赵国泽,汤吉,等.大地电磁自适应正则化反演算法[J].地球物理学报,2005,48(4):937-946.CHEN Xiaobin,ZHAO Guoze,TANG Ji,et al.An adaptive regularized inversion algorithm for magnetotelluric data[J].Chinese Journal of Geophysics,2005,48(4):937-946.
[29]
李曼,林文东.基于最小支持的2.5维直流电阻率正则化反演研究[J].东华理工大学学报(自然科学版),2014,37(3):292-298.LI Man,LIN Wendong.The study of regularization inversion for 2.5 dimension DC resistivity based on minimum support stabilizing factor[J].Journal of East China Institute of Technology (Natural Science),2014,37(3):292-298.
[30]
Constable S C,Parker R C and Constable G G.Occam's inversion:a practical algorithm for generating smooth models from EM sounding data[J].Geophysics,1987,52(3):289-300.
[31]
Smith J T,and Booker J R.Rapid inversion of two and three-dimensional magnetotelluric data[J].Journal of Geophysical Research,1991,96:3905-3922.
[32]
Oldenburg D,Li Y.Inversion of induced polarization data[J].Geophysics,1994,59(9):1327-1341.
[33]
Last B J and Kubik K.Compact gravity inversion[J].Geophysics,1983,48(6):713-721.
[34]
Guillen A and Menichetti V.Gravity and magnetic inversion with minimization of a specific functional[J].Geophysics,1984,49(8):1354-1360.
[35]
Portniaguine O,and Zhdanov M S.Focusing geophysical inversion images[J].Geophysics,1999,64(3):874-887.
Zhdanov M S.Geophysical Inverse Theory and Regularization Problems[M].Elsevier Science,Netherlands,2002.
[38]
Mackie R L,Madden T R.Three-dimensional magnetotelluric inversion using conjugate gradients[J].Geophysical Journal International,1993,115:215-229.
[39]
Rodi W L,and Mackie R L.Nonlinear conjugate gradients algorithm for 2-D magnetotelluric inversion[J].Geophysics,2001,66(1):174-187.
[40]
McGillivray P R,and Oldenburg D W.Methods for calculating Fréchet derivatives and sensitivities for the non-linear inverse problems:a comparative study[J].Geophysical Prospecting,1990,38(5):499-524.
[41]
Colin G.Approximate Sensitivities for the Multi-Dimensional Electromagnetic Inverse Problem[D].University of Edinburgh,1995.
[42]
柳建新.不完全LU分解预处理的BICGSTAB算法在大地电磁二维正演模拟中的应用[J].中南大学学报(自然科学版),2009,40(2):484-491.LIU Jianxin.Application of BICGSTAB algorithm with incomplete LU decomposition preconditioning to two-dimensional magnetotelluric forward modeling[J].Journal of Central South University (Science and Technology),2009,40(2):484-491.
[43]
Saad Y.Iterative Methods for Sparse Linear Systems[M].PWS Pub,2009.
[44]
Sethian J A,and Vladimirsky A.Fast methods for the Eikonal and related Hamilton-Jacobi equations on unstructured meshes[J].Proceedings of the National Acdemy of Sciences of the United States of America.2000,97(11):5699-5703.
[45]
Peter G,Leliéver,and Colin G.Farquharson gradient and smoothness regularization operators for geophysical inversion on unstructured meshes[J].Geophysical Journal International,2013,195:330-341.