Forward modeling formulae without analytic singularities for third-order magnetic gradient tensor of uniformly magnetized vertical cuboid
WANG Zhen1, DU Jinsong1,2,3, QIAN Bohao1, YUAN Changqing1, HU Zhengwang1,2
1. School of Geophysics and Geomatics, China University of Geosciences, Wuhan, Hubei 430074, China; 2. Hubei Subsurface Multiscale Imaging Key Laboratory, China University of Geosciences(Wuhan), Wuhan, Hubei 430074, China; 3. State Key Laboratory of Geological Processes and Mineral Resources, China University of Geosciences(Wuhan), Wuhan, Hubei 430074, China
Abstract:As one of the most basic and most commonly used source models, the uniformly magnetized vertical cuboid is usually adopted for forward modeling and three-dimensional inversion of magnetic anomalies. Thus, its forward modeling formulae are of great importance both in theory and practice. The third-order magnetic gradient tensor (MGT) is a newly proposed concept and may be an important development direction of satellite magnetic and aeromagnetic surveys. Furthermore, many magnetic transforms for anomaly interpretation involve the forward modeling of several elements of the third-order MGT. Therefore, given the previous studies, the forward modeling formulae without analytic singularities for third-order MGT of the uniformly magnetized vertical cuboid are derived. Then, the correctness of the derived formulae is verified. Finally, the corresponding formulae are applied to the forward modeling of the analytic signal (AS) of the first-order vertical derivatives of total magnetic intensity (TMI) anomalies, first-order derivatives of the invariants of second-order MGT, and the first-order derivatives of normalized source strength (NSS). With high accuracy and high efficiency, the analytic formulae for forward modeling lay the foundation for future magnetic exploration utilizing the third-order MGT and the magnetic anomaly interpretation by magnetic transforms.
钟炀,管彦武,石甲强,等.长方体全张量磁梯度的两种正演公式对比[J].世界地质,2019,38(4):1073-1081,1090.ZHONG Yang,GUAN Yanwu,SHI Jiaqiang,et al.Comparison of two forward modeling formulas for full tensor magnetic gradient of cuboid[J].Global Geology,2019,38(4):1073-1081,1090.
[2]
PEDERSEN L B,RASMUSSEN T M.The gradient tensor of potential field anomalies:some implications on data collection and data processing of maps[J].Geophysics,1990,55(12):1558-1566.
[3]
SCHMIDT P W,CLARK D A.The magnetic gradient tensor:its properties and uses in source characterization[J].The Leading Edge,2006,25(1):75-78.
[4]
FOSS C.Improvements in source resolution that can be expected from inversion of magnetic field tensor data[J].The Leading Edge,2006,25(1):81-84.
[5]
张昌达.航空磁力梯度张量测量——航空磁测技术的最新进展[J].工程地球物理学报,2006,3(5):354-361.ZHANG Changda.Airborne tensor magnetic gradio-metry:the latest progress of airborne magnetometric technology[J].Chinese Journal of Engineering Geophysics,2006,3(5):354-361.
[6]
吴招才,刘天佑.磁力梯度张量测量及应用[J].地质科技情报,2008,27(3):107-110.WU Zhaocai,LIU Tianyou.Magnetic gradient tensor:its properties and uses in geophysics[J].Geological Science and Technology Information,2008,27(3):107-110.
[7]
DENG X,SHEN W.Evaluation of gravitational curvatures of a tesseroid in spherical integral kernels[J].Journal of Geodesy,2018,92(4):415-429.
[8]
杜劲松,邱峰.重力位三阶梯度张量及其勘探能力初步分析[J].大地测量与地球动力学,2019,39(4):331-338.DU Jinsong,QIU Feng.Third-order gradient tensor of gravitational potential and preliminary analysis of its exploration capacity[J].Journal of Geodesy and Geodynamics,2019,39(4):331-338.
[9]
YIN G,ZHANG Y,FAN H,et al.Magnetic dipole localization based on magnetic gradient tensor data at a single point[J].Journal of Applied Remote Sensing,2014,8(1):083596.
[10]
DENG X,SHEN W,KUHN M,et al.Magnetic curvatures of a uniformly magnetized tesseroid using the Cartesian kernels[J].Surveys in Geophysics,2020,41(5):1075-1099.
[11]
SUI Y,LESLIE K,CLARK D.Multiple-order magnetic gradient tensors for localization of a magnetic dipole[J].IEEE Magnetics Letters,2017,8:1-5.
[12]
ZHANG H,RAVAT D,MARANGONI Y R,et al.Improved total magnetization direction determination by correlation of the normalized source strength derivative and the reduced-to-pole fields[J].Geophy-sics,2018,83(6):J75-J85.
[13]
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.
[14]
CLARK D A.New methods for interpretation of magnetic vector and gradient tensor data Ⅱ:application to the Mount Leyshon anomaly,Queensland, Australia [J].Exploration Geophysics,2013,43(4):114-127.
[15]
CLARK D A.New methods for interpretation of magnetic vector and gradient tensor data I: eigenvector analysis and the normalised source strength[J].Exploration Geophysics,2012,43(4):267-282.
[16]
管志宁.地磁场与磁力勘探[M].北京:地质出版社,2005.GUAN Zhining.Geomagnetic Field and Magnetic Exploration[M].Geological Publishing House,Beijing,2005.
李金朋,范红波,刘利,等.多参数约束磁性体三维形态反演[J].石油地球物理勘探,2021,56(2):407-418.LI Jinpeng,FAN Hongbo,LIU Li,et al.Multi-para-meter constrained three-dimensional shape inversion of magnetic bodies[J].Oil Geophysical Prospecting,2021,56(2):407-418.
[19]
LI Y,OLDENBURG D W.3-D inversion of magnetic data[J].Geophysics,1996,61(2):394-408.
[20]
BHATTACHARYYA B K.Magnetic anomalies due to prism-shaped bodies with arbitrary polarization[J].Geophysics,1964,29(4):517-531.
[21]
GOODACRE A K.Some comments on the calculation of the gravitational and magnetic attraction of a homo-geneous rectangular prism[J].Geophysical Prospecting,1973,21(1):66-69.
[22]
PLOUFF D.Gravity and magnetic fields of polygonal prisms and application to magnetic terrain corrections[J].Geophysics,1976,41(4):727-741.
[23]
KUNARATNAM K.Simplified expressions for the magnetic anomalies due to vertical rectangular prisms[J].Geophysical Prospecting,1981,29(6):883-890.
[24]
RAO D B,BABU N R.A rapid method for three-dimensional modeling of magnetic anomalies[J].Geophysics,1991,56(11):1729-1737.
[25]
郭志宏,管志宁,熊盛青.长方体ΔT场及其梯度场无解析奇点理论表达式[J].地球物理学报,2004,47(6):1131-1138.GUO Zhihong,GUAN Zhining,XIONG Shengqing.Cuboid ΔT and its gradient forward theoretical expressions without analytical odd points[J].Chinese Journal of Geophysics,2004,47(6):1131-1138.
[26]
骆遥,姚长利.长方体磁场及其梯度无解析奇点表达式理论研究[J].石油地球物理勘探,2007,42(6):714-719.LUO Yao,YAO Changli.Theoretical study on cuboid magnetic field and its gradient expression without ana-lytic singular point[J].Oil Geophysical Prospecting,2007,42(6):714-719.
[27]
KUANG X,YANG H,ZHU X,et al.Singularity-free expression of magnetic field of cuboid under undula-ting terrain[J].Applied Geophysics,2016,13(2):238-248.
[28]
邰振华,柴琳,黄德智,等.均匀磁化长方体的磁张量无解析奇点正演方法[J].黑龙江科技大学学报,2022,32(3):286-292.TAI Zhenhua,CHAI Lin,HUANG Dezhi,et al.Singularity-free forward calculation method of uniformly magnetized cuboid magnetic tensors[J].Journal of Heilongjiang University of Science and Technology,2022,32(3):286-292.
[29]
黄远生,王彦国,罗潇.基于归一化磁源强度垂向差分的磁源参数快速估计方法[J].物探与化探,2021,45(6):1588-1596.HUANG Yuansheng,WANG Yanguo,LUO Xiao.A fast estimation method of magnetic-source parameters based on the vertical difference of normalized source strength[J].Geophysical and Geochemical Exploration,2021,45(6):1588-1596.
[30]
BLAKELY R J.Potential Theory in Gravity and Magnetic Applications[M].Cambridge University Press,Cambridge,1995.
[31]
NABIGHIAN M N.The analytic signal of two-dimensional magnetic bodies with polygonal cross-section:its properties and use for automated anomaly interpretation[J].Geophysics,1972,37(3):507-517.
[32]
NABIGHIAN M N.Toward a three-dimensional automatic interpretation of potential field data via genera-lized Hilbert transforms:fundamental relations[J].Geophysics,1984,49(6):780-786.
[33]
AGARWAL B N P,SHAW R K.Comment on 'An analytic signal approach to the interpretation of total field magnetic anomalies' by Shuang Qin[J].Geophy-sical Prospecting,1996,44(5):911-914.
[34]
LI X.Understanding 3D analytic signal amplitude[J].Geophysics,2006,71(2):L13-L16.
[35]
HSU S K,SIBUET J C,SHYU C T.High-resolution detection of geologic boundaries from potential-field anomalies:an enhanced analytic signal technique[J].Geophysics,1996,61(2):373-386.
[36]
BAUR O,SNEEUW N,GRAFAREND E W.Methodo-logy and use of tensor invariants for satellite gravity gradiometry[J].Journal of Geodesy,2008,82(4):279-293.
[37]
DENG X,SHEN W,YANG M,et al.First-order derivatives of principal and main invariants of magnetic gradient tensor of a uniformly magnetized tesseroid and spherical shell[J].Surveys in Geophysics,2022,43(4):1233-1262.
[38]
GUO L,GAO R,ZHANG G.Estimating the magnetization direction of sources through correlation between reduced-to-pole anomaly and normalized source strength[J].Applied Mechanics and Materials,2014,644-650:3793-3796.
[39]
PILKINGTON M,BEIKI M.Mitigating remained magnetization effects in magnetic data using the normalized source strength[J].Geophysics,2013,78(3):J25-J32.