王绪本1,2,,,
曹辉1,2,
蓝星1,4,
段长生5,
罗威1,4,
原健龙1
1. 成都理工大学地球物理学院, 成都 610059
2. 油气藏地质及开发工程国家重点实验室(成都理工大学), 成都 610059
3. 中国地质调查局成都地质调查中心, 成都 610081
4. 四川省冶勘设计集团有限公司, 成都 610051
5. 赣中南地质矿产勘查研究院, 南昌 330029
基金项目: 国家自然科学基金重点研究计划(91755215),国家自然科学基金面上项目(41974090),国家重点研发计划项目(2017YFC0601504),油气藏地质及开发工程国家重点实验室(成都理工大学)开放基金(PLC2020003)联合资助
详细信息
作者简介: 王堃鹏, 男, 1989年生, 讲师, 主要从事电磁勘探方法正反演研究工作.E-mail: xfnwkp@163.com
通讯作者: 王绪本, 男, 1956年生, 教授, 博士生导师, 主要从事地球物理电磁法理论及应用研究.E-mail: wxb@cdut.edu.cn
中图分类号: P631;P319收稿日期:2020-08-27
修回日期:2020-11-08
上线日期:2021-04-10
Magnetotelluric axial anisotropic parallelized 3D inversion based on cross-gradient structural constraint
WANG KunPeng1,2,3,,WANG XuBen1,2,,,
CAO Hui1,2,
LAN Xing1,4,
DUAN ChangSheng5,
LUO Wei1,4,
YUAN JianLong1
1. College of Geophysics, Chengdu University of Technology, Chengdu 610059, China
2. State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation(Chengdu University of Technology), Chengdu 610059, China
3. Chengdu Center, China Geological Survey, Chengdu 610081, China
4. Sichuan Metallurgical Geological Survey and Design Group Co., Ltd., Chengdu 610051, China
5. Ganzhongnan Institute of Geology and Mineral Exploration, Nanchang 330029, China
More Information
Corresponding author: WANG XuBen,E-mail:wxb@cdut.edu.cn
MSC: P631;P319--> Received Date: 27 August 2020
Revised Date: 08 November 2020
Available Online: 10 April 2021
摘要
摘要:各向异性介质对大地电磁观测数据的影响往往不可忽略,因此需要提高大地电磁各向异性三维反演的可靠性和有效性.为了满足大地电磁各向异性三维反演的需求,本文研究了一种基于交叉梯度结构约束的大地电磁主轴各向异性并行三维反演算法.根据大地电磁平面波理论假设,正演方程采用背景场与二次场分离的计算方式,二次场利用交错网格有限差分法求解.由于各向异性反演的多解性,本文将各向异性介质简化为主轴各向异性,并在此基础上进一步采用有限内存拟牛顿LBFGS法实现三维各向异性反演.为了提高各向异性反演的分辨率,反演目标函数中引入交叉梯度项,利用先验的结构信息,对三个方向的电阻率参数进行结构约束,最终的反演进一步利用MPI(Message Passing Interface,消息传递接口)技术实现分频并行计算,测试结果显示并行接近线性加速比.
关键词: 大地电磁法/
二次场/
主轴各向异性/
交叉梯度/
并行
Abstract:The influence of anisotropic media on the magnetotelluric observation data is often not negligible, so it is of great significance to improve the reliability and effectiveness of the 3D inversion of the magnetotelluric anisotropy. In order to meet the needs of 3D inversion of magnetotelluric anisotropy, this paper studies a parallel 3D inversion algorithm of magnetotelluric anisotropy based on cross-gradient structural constraint. Based on the assumption of the magnetotelluric plane wave theory, the forward equation uses the separation of the primary field and the secondary field, and the secondary field uses the staggered-grid finite difference method. Due to the multi-solution of anisotropic inversion, this paper simplifies the anisotropic medium to the axial anisotropy, and on this basis, further adopts the LBFGS (Limited-memory Broyden-Fletcher-Goldfarb-Shanno) method to realize the 3D anisotropic inversion. In order to improve the resolution of the anisotropic inversion, the inversion objective function introduces the cross-gradient term, and uses a priori structural information to constrain structure of the resistivity parameters in the three directions. The final inversion further uses MPI (Message Passing Interface) technology to realize the frequency division parallel. The test results show that the parallel is close to the linear acceleration ratio.
Key words:Magnetotelluric/
Secondary field/
Axial anisotropy/
Cross-gradient/
Parallel
PDF全文下载地址:
http://www.geophy.cn/data/article/export-pdf?id=6076842dc23e6710c26d5b41