谭捍东,,
彭淼,
孔文新,
吴萍萍
中国地质大学(北京)地球物理与信息技术学院, 北京 100083
基金项目: 国家重点研发计划课题"综合航空地球物理数据处理、解释和管理软件平台研发"（2017YFC0602204）和国家自然科学基金项目（41830429，41804130，41704055）资助
详细信息
作者简介: 闫政文, 男, 1993年生, 硕士研究生, 主要从事地球物理联合反演和电磁学算法理论研究.E-mail:11930895@mail.sustech.edu.cn
通讯作者: 谭捍东, 男, 1966年生, 教授, 博导, 主要从事电法勘探理论及应用研究.E-mail:thd@cugb.edu.cn
中图分类号: P631收稿日期:2018-06-11
修回日期:2019-06-12
上线日期:2020-02-05
Three-dimensional joint inversion of gravity, magnetic and magnetotelluric data based on cross-gradient theory
YAN ZhengWen,TAN HanDong,,
PENG Miao,
KONG WenXin,
WU PingPing
School of Geophysics and Information Technology, China University of Geosciences, Beijing 10083, China
More Information
Corresponding author: TAN HanDong,E-mail:thd@cugb.edu.cn
MSC: P631--> Received Date: 11 June 2018
Revised Date: 12 June 2019
Available Online: 05 February 2020
摘要
摘要:为了降低单一地球物理方法反演的多解性及受噪声的影响程度，本文围绕重力、磁法和大地电磁法开展了三维联合反演的研究.重、磁采用基于对数障碍法的正则化反演算法，大地电磁使用limited-memory Broyden-Fletcher-Goldfarb-Shanno（L-BFGS）反演算法，引入交叉梯度函数实现了三种物性结构的相互耦合，最终开发出一套重磁电三维联合反演算法，并实现MPI并行加速计算.通过理论模型算例验证了算法的准确性，结果表明：不论是单棱柱体模型还是组合棱柱体模型，联合反演结果相较单独反演对于异常体的空间形态刻画以及物性数值恢复具有较好的提升；单棱柱体模型算例使得异常体的物性参数（密度、磁化率和电阻率）更加接近于真实的物性参数；组合棱柱体模型的联合反演结果不仅仅消除了围岩物性参数的假异常，而且还增强了异常体边界结构的恢复程度.
关键词: 联合反演/
三维交叉梯度函数/
L-BFGS反演/
大地电磁测深
Abstract:In order to reduce the uncertainty of inversion from a single geophysical method and the influence of noise, this paper conducts three-dimensional joint inversion of the gravity, magnetic and magnetotelluric data. We use regularization algorithm for Gravity and magnetic inversions based on logarithm barrier method. We use limited-memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) algorithm for Magnetotellurics inversion. In the joint inversion, we add cross-gradient function to achieve the structural coupling of these three kinds of physical models (density, magnetic susceptibility and resistivity). Finally, we developed three-dimensional joint inversion of gravity, magnetic and magnetoelectric data. Furthermore, we used the MPI parallel algorithm to speed up the calculation. The accuracy of the algorithm is verified by the inversions of several theoretical models. By inverting synthetic data from a single-prism model and a multiple-prism model, we conclude that joint inversion has better performance in recovering anomalous prism body's spatial morphology. The resistivity value of the single-prism is much closer to the true model in joint inversion. Besides, joint inversion not only removes electrical false anomaly of backgrounds to a certain degree, but also largely improves the resolution of density and magnetic susceptibility at deep depth in the multiple prism models. The results show that the joint inversion is better than the single inversion.
Key words:Joint inversion/
Three-dimensional cross gradient function/
L-BFGS inversion/
Magnetotelluric
PDF全文下载地址:
http://www.geophy.cn/data/article/export-pdf?id=dqwlxb_15352
