A multiple level-set method for 3D boundary inversion of magnetic data
XIAO Xiao1,2,3,4, DUAN Ya-ting3, HU Shuanggui3, TANG Jingtian1,2,3,4, XIE Yong5, LIU Changsheng5
1. Key Laboratory of Metallogenic Prediction of Nonferrous Metals and Geological Environment Monitoring (Central South University), Ministry of Education, Changsha, Hunan 410083, China; 2. Hunan Key Laboratory of Nonferrous Resources and Geological Hazards Exploration, Changsha, Hunan 410083, China; 3. School of Geosciences and Info-physics, Central South University, Changsha, Hunan 410083, China; 4. Technical Innovation Center of Coverage Area Deep Resources Exploration, Ministry of Natural Resources, Changsha, Hunan 410083, China; 5. Changsha Aeronautical Vocational and Technical College, Changsha, Hunan 410124, China
Abstract:Present algorithms for inverting the boundaries of magnetic targets use only two level sets. There are usually multiple magnetic geological bodies with different susceptibility in actual exploration. This paper proposes a new multiple level-set inversion algorithm for 3D inverting the position and geometry of magnetic targets with known susceptibility. First, based on the principle of multiple level sets, an objective function based on a multiple level sets function is established. Then, during the inversion process, an arbitrary tetrahedron element magnetic analytical solution algorithm is used for high-precision forward calculation, and the physical property mapping between forward and inverse grids is used to realize the independent operation of the forward grids and the inverse grids. The weighted essentially non-oscillatory scheme (WENO) is introduced to update and reinitialize the level set function, thereby improving the reliability and efficiency of the inversion. Finally, the effectiveness of the algorithm is verified by theoretical models with various numbers of level sets. The results show that the inversion based on the multiple level-set method has strong flexibility, and can automatically merge and separate regions in the model, thereby changing the model's topology with no need to manually reparameterize it. The accuracy of the inversion results is improved, and the boundary of the anomalous body obtained from the inversion accords well with the true boundary.
冯杰,欧洋,赵勇,等.三维井地磁测联合约束反演[J]. 地球物理学报, 2019, 62(10):3686-3698.FENG Jie, OU Yang, ZHAO Yong, et al. 3D joint constrained inversion of borehole and ground magne-tic data[J]. Chinese Journal of Geophysics, 2019, 62(10):3686-3698.
李志鹏,高嵩,王绪本.特殊区域旋翼无人机航磁测量研究[J].地球物理学报, 2018, 61(9):317-326.LI Zhipeng, GAO Song, WANG Xuben. New me-thod of aeromagnetic surveys with rotorcraft UAV in particular areas[J]. Chinese Journal of Geophysics, 2018, 61(9):317-326.
[5]
Cribb J. Application of the generalized linear inverse to the inversion of static potential data[J]. Geophy-sics, 1976, 41(6):1365-1369.
[6]
Li Y and Oldenburg D W. 3-D inversion of magnetic data[J]. Geophysics, 1996, 61(2):394-408.
[7]
Pilkington M. 3-D magnetic imaging using conjugate gradients[J]. Geophysics, 1997, 62(4):1132-1142.
[8]
Miguel B J. Joint inversion of gravity and magnetic data under lithologic constraints[J]. The Leading Edge, 2001, 20(8):877-881.
[9]
Fregoso E and Gallardo L A. Cross-gradients joint 3D inversion with applications to gravity and magnetic data[J]. Geophysics, 2009, 74(4):L31-L42.
[10]
管志宁,安玉林.频率域磁性单界面反演方法[J].地球物理学报, 1984, 27(2):190-204.GUAN Zhining, An Yulin. Inversion method of a single magnetic boundary in the frequency domain[J]. Chinese Journal of Geophysics, 1984, 27(2):190-204.
[11]
Pilkington M, Crossley D J. Determination of crustal interface topography from potential fields[J]. Geophysics, 1986, 51(6):1277-1284.
[12]
Wang X, Hansen R O. Inversion for magnetic ano-malies of arbitrary three-dimensional bodies[J]. Geophysics, 1990, 55(10):1321-1326.
[13]
刘沈衡,吴燕岗,李永朴.利用欧拉算法反演磁性界面[J].石油物探, 2000, 39(2):101-106.LIU Shenheng, WU Yan'gang, LI Yongpu. Magne-tic interface inversion by Euler algorithm[J]. Geophysical Prospecting for Petroleum, 2000, 39(2):101-106.
[14]
Fullagar P K, Pears G A, Mcmonnies B. Constrained inversion of geologic surfaces:Pushing the boundaries[J]. The Leading Edge, 2008, 27(1):98-105.
[15]
赵文举,刘云祥,陶德强,等. BP神经网络磁性体顶面埋深预测方法[J].石油地球物理勘探, 2020, 55(5):1139-1148.ZHAO Wenju, LIU Yunxiang, TAO Deqiang, et al. Prediction of magnetic body top based on BP neural network[J]. Oil Geophysical Prospecting, 2020, 55(5):1139-1148.
[16]
郑强,郭华,张贵宾,等.基于磁力梯度全张量特征值的均衡边界识别方法[J].石油地球物理勘探, 2020, 55(2):454-464.ZHENG Qiang, GUO Hua, ZHANG Guibin, et al. Balanced edge detection based on eigenvalues of full magnetic gradient tensor[J]. Oil Geophysical Prospecting, 2020, 55(2):454-464.
[17]
Chang C J, Hsieh, Jun W C, et al. Tracking multiple moving objects using a level-set method[J]. International Journal of Pattern Recognition and Artificial Intelligence, 2004, 18(2):101-125.
[18]
Mansouri A R, Konrad J. Multiple motion segmentation with level sets[J]. IEEE Transactions on Image Processing, 2003, 12(2):201-220.
[19]
Sussman M, Smereka P, Osher S. A level set ap-proach for computing solutions to incompressible two-phase flow[J]. Journal of Computational Phy-sics, 1994, 114:146-159.
[20]
Tsai Y H R, Osher S. Total variation and level set methods in image science[J]. Acta Numerica, 2005, 14(2):509-573.
[21]
Zhao H K, Chan T, Merriman B, et al. A variational level set approach to multiphase motion[J]. Journal of Computational Physics, 1996, 127:179-195.
[22]
Isakov V, Leung S, Qian J. A fast local level set method for inverse gravimetry[J]. Communications in Computational Physics, 2011, 10(4):1044-1070.
[23]
Zheglova P, Farquharson C G and Hurich C. 2-D re-construction of boundaries with level set inversion of traveltimes[J]. Geophysical Journal International, 2013, 192(2):688-698.
[24]
Li W B, Lu W T, Qian J L, et al. A multiple level set method for three dimensional inversion of magne-tic data[J]. Geophysics, 2017, 82(5):J61-J81.
[25]
Li W B, Qian J L, Li Y G. Joint inversion of surface and borehole magnetic data:A level-set approach[J]. Geophysics, 2020, 85(1):J15-J32.
[26]
Ren Z, Chen C, Tang J, et al. Closed-form formula of magnetic gradient tensor for a homogeneous polyhedral magnetic target:A tetrahedral grid example[J]. Geophysics, 2017, 82(6):V1-V34.
[27]
Osher S and Fedkiw R. Level set methods and dynamic implicit surfaces[J]. Springer Science & Business Media, 2006, 153:230-240.
[28]
Lu W T, Leung S Y, Qian J L. An improved fast local level set method for three-dimensional inverse gravimetry[J]. Inverse Problems and Imaging, 2015, 9(1):479-509.
[29]
Jiang G S and Peng D. Weighted ENO schemes for Hamilton-Jacobi equations[J]. SIAM Journal on Scientific Computing, 2000, 21(6):2126-2143.
[30]
Guang S J, Chi W S. Efficient implementation of weighted ENO schemes[J]. Journal of Computational Physics, 1996, 126:202-228.
[31]
Liu X D, Osher S, Chan T. Weighted essentially non-oscillatory schemes[J]. Journal of Computational Phy-sics, 1996, 115:200-212.
[32]
Shu C W, Osher S. Efficient implementation of essentially non-oscillatory shock-capturing schemes[J]. Journal of Computational Physics, 1988, 77:439-471.
[33]
Frederic G R F, Osher S. A review of level-set me-thods and some recent applications[J]. Journal of Computational Physics, 2018, 353:82-109.
[34]
Courant R, Friedrichs K, Lewy H. On the partial difference equations of mathematical physics[J]. IBM Journal of Research & Development, 1967, 11(2):215-234.
[35]
Li W, Leung S. A fast local level set adjoint state method for first arrival transmission traveltime tomography with discontinuous slowness[J]. Geophysical Journal International, 2013, 195:582-596.
[36]
Li W, Leung S and Qian J L. A level-set adjoint-state method for crosswell transmission reflection traveltime tomography[J]. Geophysical Journal International, 2014, 199:348-367.