Multi-parameter constrained three-dimensional shape inversion of magnetic bodies
LI Jinpeng1, FAN Hongbo2, LIU Li1, ZHANG Yingtang2, LI Zhining2, WU Bingyang3
1. Unit 93114 of PLA, Beijing 100195, China; 2. Army Engineering University(Shijiazhuang), Shijiazhuang, Hebei 050003, China; 3. Beijing Institute of Remote Sensoring Information, Beijing 100192, China
Abstract：Traditionally, the magnetic and geometric parameters for shape inversion of magnetic bodies are defined manually. We propose a three-dimensional shape inversion method of magnetic bodies based on multi-parameter constraints. First, an objective function is established through L1 norm and regularization constraint function, and the optimal solution is selected in the module to be updated to update the inversion model; then the method for estimating the multiple parameters of the magnetic target is proposed, which can directly estimate the position, magnetization direction and magnetization intensity of the magnetic target without prior information; and finally, the magnetic component (Bz) along z and magnetic gradient tensor are used for joint inversion, and the inversion result is compared with normalized source strength and total magnitude. The simulation and experimental results show that, even without prior information, the method we proposed can effectively estimate the magnetic parameters of the magnetic target and accurately obtain the inversion results of the magne-tic target. At the same time, under the condition of remanent magnetization, compared with the traditional normalized source strength and total magnitude, this method can effectively improve the accuracy of three-dimensional inversion of magnetic bodies.
郑建拥,范红波,张琪,等.利用磁梯度张量识别地下小型目标体[J].石油地球物理勘探,2019,54(3):692-699.ZHENG Jianyong,FAN Hongbo,ZHANG Qi,et al. Underground small target recognition using magnetic gradient tensor[J].Oil Geophysical Prospecting,2019,54(3):692-699.
[2]
李青竹,李志宁,张英堂,等.基于二阶磁张量欧拉反褶积的磁源单点定位方法[J].石油地球物理勘探,2019,54(4):915-924.LI Qingzhu,LI Zhining,ZHANG Yingtang,et al. Magnetic source single-point positioning based on se-cond-order magnetic tensor Euler deconvolution[J].Oil Geophysical Prospecting,2019,54(4):915-924.
[3]
李华东,于鹏,刘振友.基于随机分布共网格模型的重磁电震联合反演技术及应用[J].石油地球物理勘探,2015,50(4):742-748.LI Huadong,YU Peng,LIU Zhenyou. Joint inversion of gravity,magnetotelluric and seismic data based on common gridded model with random distributions[J].Oil Geophysical Prospecting,2015,50(4):742-748.
[4]
习宇飞,刘天佑,刘双.井中磁测三分量联合反演[J].石油地球物理勘探,2012,47(2):344-352.XI Yufei,LIU Tianyou,LIU Shuang. A joint inversion method for borehole magnetic three-component data[J].Oil Geophysical Prospecting,2012,47(2):344-352.
[5]
Schneider M,Stolz R,Linzen S,et al. Inversion of geo-magnetic full-tensor gradiometer data[J].Journal of Applied Geophysics,2013,92:57-67.
[6]
李金朋,张英堂,范红波,等.基于χ2准则的磁梯度张量3D聚焦反演方法[J].武汉大学学报(信息科学版),2018,43(2):255-261.LI Jinpeng,ZHANG Yingtang,FAN Hongbo,et al. Three-dimensional focusing inversion of magnetic gradient tensor data based on the χ2 principle[J].Geomatics and Information Science of Wuhan University,2018,43(2):255-261.
[7]
索奎,张贵宾,梅岩辉,等.重磁三维反演伊犁盆地中部密度和磁性结构[J].地球物理学报,2018,61(8):3410-3419.SUO Kui,ZHANG Guibin,MEI Yanhui,et al. Density and magnetic susceptibility distribution of central Yili Basin by three-dimensional of gravity and magnetic data[J].Chinese Journal of Geophysics,2018,61(8):3410-3419.
[8]
Li S L,Li Y. Inversion of magnetic anomaly on rugged observation surface in the presence of strong re-manent magnetization[J].Geophysics,2014,79(2):J11-J19.
[9]
Antonio G,Camacho,Fuensanta G,et al. A 3-D gra-vity inversion tool based on exploration of model possibilities[J].Computer & Geosciences,2002,28(1):191-204.
[10]
Uieda L,Barbosa V C F. Use of the "shape-of-ano-maly" data misfit in 3D inversion by planting anomalous densities[C].SEG Technical Program Expanded Abstracts,2012,31:3160-3164.
[11]
曹书锦,朱自强,鲁光银.重力梯度张量种植反演[J].中国有色金属学报,2017,27(5):997-1005.CAO Shujin,ZHU Ziqiang,LU Guangyin. Planting inversion of tensor gravity data[J].The Chinese Journal of Nonferrous Metals,2017,27(5):997-1005.
[12]
Uieda L,Barbosa V C F. Robust 3D gravity gradient inversion by planting anomalous densities[J].Geophysics,2012,77(4):G55-G66.
[13]
Cardarelli E, Fischanger F. 2D data modelling by electrical resistivity tomography for complex subsurface geology[J].Geophysical Prospecting,2006,54(2):121-133.
[14]
Zhdanov M S,Liu X,Wilson G A,et al. 3D migration for rapid imaging of total-magnetic-intensity data[J].Geophysics,2012,77(2):J1-J5.
[15]
Zhdanov M S,Liu X,Wilson G A,et al. Potential field migration for rapid imaging of gravity gradiometry data[J].Geophysical Prospecting,2011,59:1052-1071.
[16]
尹刚,张林.基于形态反演的小尺度磁性目标三维重建方法[J].上海交通大学学报,2019,53(9):1074-1083.YIN Gang,ZHANG Lin. Three-dimensional reconstruction method for small scale magnetic target based on shape inversion[J].Journal of Shanghai Jiao Tong University,2019,53(9):1074-1083.
[17]
Beiki M,Clark D A,Austin J R,et al. Estimating source location using normalized magnetic source strength calculated from magnetic gradient tensor data[J].Geophysics,2012,77(6):J23-J37.
[18]
Zhou J,Meng X. Three-dimensional cross-gradient joint inversion of gravity and normalized magnetic source strength data in the presence of remanent magnetization[J].Journal of Applied Geophysics,2015,119:51-60.
[19]
Guo L,Meng X,Zhang G. Three-dimensional correlation imaging for total amplitude magnetic anomaly and normalized source strength in the presence of strong remanent magnetization[J].Journal of Applied Geophysics,2014,111:121-128.
[20]
李金朋,张英堂,范红波,等.基于磁传感器阵列的多磁源参数反演方法[J].仪器仪表学报,2019,40(10):28-37.LI Jinpeng,ZHANG Yingtang,FAN Hongbo,et al. Multi-magnetic source parameter inversion method based on magnetic sensor array[J].Chinese Journal of Scientific Instrument,2019,40(10):28-37.
[21]
谢汝宽,王平,刘浩军.基于最小反演拟合差的重磁场源深度计算方法[J].地球物理学报,2016,59(2):711-720.XIE Rukuan,WANG Ping,LIU Haojun. Depth estimation of potential field by minimum inversion fitting error[J].Chinese Journal of Geophysics,2016,59(2):711-720.
[22]
Blakely R J. Potential Theory in Gravity and Magne-tic Applications[M].Cambridge University Press,Cambridge,1996.
[23]
Yatong C,Lianghui G. A wavenumber-domain iterative approach for 3D imaging of magnetic anomalies and gradients with depth constraints[J].Journal of Geophysics and Engineering,2019,16(6):1-16.
