Two-dimensional magnetotelluric smooth focusing inversion based on optimization strategy
BAI Ningbo1, ZHOU Junjun2, HU Xiangyun1,2
1. Key Laboratory of Geological Survey and Evaluation of Ministry of Education, China University of Geosciences(Wuhan), Wuhan, Hubei 430074, China; 2. Institute of Geophysics & Geomatics, China University of Geosciences(Wuhan), Wuhan, Hubei 430074, China
Abstract:On the basis of previous studies, this paper proposes a new inversion objective functional with the purposes of realizing rapid and stable inversion and obtaining clear geological interfaces. It adopts the smoothest model and the minimum support gradient model functional to constrain the data objective functional. Solved by the Gauss-Newton method, the new inversion objective functional enables the smooth focusing inversion of two-dimensional magnetotelluric data. The smooth focusing inversion can not only present clear geological interfaces but also avoid the distortion of structural morphology caused by focused inversion to a certain extent. In the process of inversion iteration, we adopt the optimization strategy of improving the Morozov discrepancy principle with the Nelder-Mead optimization algorithm to obtain the appropriate regularization factor, which greatly accelera-tes the inversion convergence. Finally, the proposed inversion method is verified with a typical model and real data and also compared with other inversion strategies. The inversion results show that for typical model inversion, the algorithm in this paper outperforms the others in agreeing with the model, with the convergence curve decreasing rapidly and the geological body interface being clear. The inversion results of real data further verify the reliability and effectiveness of this algorithm.
Zhang R,Li T,Deng X,et al.Two-dimensional data-space joint inversion of magnetotelluric,gravity,magnetic and seismic data with cross-gradient constraints[J].Geophysical Prospecting,2020,68(2):721-731.
Degroot-Hedlin C,Constable S.Occam's inversion to generate smooth,two-dimensional models from magnetotelluric data[J].Geophysics,1990,55(12):1613-1624.
[4]
Smith J T,Booker J R.Rapid inversion of two-and three-dimensional magnetotelluric data[J].Journal of Geophysical Research:Solid Earth,1991,96(B3):3905-3922.
Rodi W,Mackie R L.Nonlinear conjugate gradients algorithm for 2-D magnetotelluric inversion[J].Geophysics,2001,66(1):174-187.
[7]
Last B J,Kubik K.Compact gravity inversion[J].Geophysics,1983,48(6):713-721.
[8]
Portniaguine O,Zhdanov M S.Focusing geophysical inversion images[J].Geophysics,1999,64(3):874-887.
[9]
Zhang L,Koyama T,Utada H,et al.A regularized three-dimensional magnetotelluric inversion with a minimum gradient support constraint[J].Geophysical Journal International,2012,189(1):296-316.
[10]
Zhdanov M S.New advances in regularized inversion of gravity and electromagnetic data[J].Geophysical Prospecting,2009,57(4):463-478.
[11]
张罗磊,于鹏,王家林,等.光滑模型与尖锐边界结合的MT二维反演方法[J].地球物理学报,2009,52(6):1625-1632.ZHANG Luolei,YU Peng,WANG Jialin,et al.Smoothest model and sharp boundary based two-dimensional magnetotelluric inversion[J].Chinese Journal of Geophysics,2009,52(6):1625-1632.
[12]
Sun S,Yin C,Liu Y,et al.3D regularized focusing inversion of gravity data with a new stabilizing functional[C].International Workshop and Gravity & Magnetic Methods and Their Applications, Chengdu, China,2015,322-325.
[13]
Hu S,Zhao Y,Qin T,et al.Travel time tomography of cross hole ground-penetrating radar based on an arctangent functional with compactness constraints[J].Geophysics,2017,82(3):H1-H14.
[14]
Zhao C,Yu P,Zhang L.A new stabilizing functional to enhance the sharp boundary in potential field regularized inversion[J].Journal of Applied Geophysics,2016,135:356-366.
[15]
Hu M,Yu P,Rao C,et al.3D sharp-boundary inversion of potential-field data with an adjustable exponential stabilizing functional[J].Geophysics,2019,84(4):J1-J15.
[16]
Xiang Y,Yu P,Zhang L,et al.Regularized magnetotelluric inversion based on a minimum support gradient stabilizing functional[J].Earth,Planets and Space,2017,69(1):1-18.
[17]
Hansen P C, Leary D P O.The use of the L-curve in the regularization of discrete ill-posed problems[J].SIAM Journal on Scientific Computing,1993,14(6):1487-1503.
[18]
吴小平,徐果明.大地电磁数据的Occam反演改进[J].地球物理学报,1998,41(4):547-554.WU Xiaoping,XU Guoming.Improvement of Occam's inversion for MT data[J].Chinese Journal of Geophysics,1998,41(4):547-554.
[19]
陈小斌,赵国泽,汤吉,等.大地电磁自适应正则化反演算法[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.
[20]
朴英哲,李桐林,刘永亮.在大地电磁二维Occam反演中求取拉格朗日乘子方法改进[J].吉林大学学报(地球科学版),2014,44(2):660-667.PAK Yongchol,LI Tonglin,LIU Yongliang.Improvement of choosing Lagrange Multiplier on MT 2D Occam inversion[J].Journal of Jilin University(Earth Science Edition),2014,44(2):660-667.
[21]
向阳,于鹏,陈晓,等.大地电磁反演中改进的自适应正则化因子选取[J].同济大学学报(自然科学版),2013,41(9):1429-1434.XIANG Yang,YU Peng,CHEN Xiao,et al.An improve adaptive regularized parameter selection in magnetotelluric inversion[J].Journal of Tongji University (Natural Science),2013,41(9):1429-1434.
[22]
Zhdanov M S.Geophysical Inverse Theory and Regularization Problems[M].Elsevier,2002.
[23]
Zhdanov M S,Tolstaya E.Minimum support nonlinear parametrization in the solution of a 3D magnetotelluric inverse problem[J].Inverse Problems,2004,20(3):937-952.
[24]
Zhang Y M S,Li Q H,Wang L P,et al.Regularization particle size inversion of PCS based on fast algorithm[C].IEEE World Congress on Intelligent Control and Automation(WCIA),2014,5641-5645.
[25]
Yildiz A R,Yildiz B S,Sait S M,et al.A new hybrid Harris Hawks-Nelder-Mead optimization algorithm for solving design and manufacturing problems[J].Materials Testing, 2019, 61(8):735-743.
[26]
Sasaki Y.Two-dimensional joint inversion of magnetotelluric and dipole-dipole resistivity data[J].Geophysics, 1989,54(2):254-262.