基于非结构有限元的带地形海洋可控源电磁法二维Occam反演  

Two-dimensional Occam Inversion of Marine Controlled-source Electromagnetic Method in Terrain Based on Unstructured Finite Element Method

作　　者：肖俊 叶益信 薛海军[2,3] 艾正敏 杜家明 Xiao Jun;Ye Yixin;Xue Haijun;Ai Zhengmin;Du Jiaming(School of Geophysics and Measurement-Control Technology, East China University of Technology, Nanchang Jiangxi 330013, China;Key Laboratory of Coal Resources Exploration and Comprehensive Utilization, Ministry of Land and Resources,Xi′an Shanxi 710005, China;Shaanxi Coalfield Geophysical Prospecting and Surveying Co., Ltd, Xi′an Shanxi, 710005, China)

机构地区：[1]东华理工大学地球物理与测控技术学院,江西南昌330013 [2]国土资源部煤炭资源勘查与综合利用重点实验室,陕西西安710005 [3]陕西省煤田物探测绘有限公司,陕西西安710005

出　　处：《工程地球物理学报》2020年第2期217-224,共8页Chinese Journal of Engineering Geophysics

基　　金：国家自然科学基金(编号:41774078)。

摘　　要：基于非结构三角形网格的自适应有限元正演模拟,对带地形海洋可控源电磁法(MCSEM)二维Occam反演算法进行了研究。该算法正演采用基于非结构网格的自适应有限元法,能用较少的单元剖分精确地模拟复杂起伏地形和构造模型。反演运用基于高斯牛顿法改进的快速Occam算法,能快速地搜索拉格朗日乘子、减少模型的搜索量。对二维海洋油气藏模型进行了反演计算,结果表明,该算法能稳定收敛到真实模型附近,具有较高的计算精度和计算效率。通过对带地形模型的反演,讨论了起伏地形下海洋可控源电磁场二维反演效果,为准确解释MCSEM数据提供了参考。Based on the adaptive finite element forward simulation with unstructured triangular meshes,the two-dimensional Occam inversion algorithm of marine controlled-source electromagnetic method(MCSEM)in terrain is studied in this paper.The adaptive finite element method based on unstructured mesh is used in the forward modeling,which can accurately simulate rough topography and complex structures with fewer element subdivisions than with structured mesh.Fast Occam algorithm based on Gauss-Newton method is used in the inversion,which can search Lagrange multiplier quickly and reduce the search amount of models.The inversion of two-dimensional offshore oil and gas reservoir model shows that the inversion can converge steadily to the real model with greater precision and efficiency.The effectiveness of MCSEM data inversionwithrugged terrain is tested,which provides a reference for accurate interpretation of MCSEM data.

关 键 词：起伏地形 非结构网格 自适应有限元 MCSEM 反演 

分 类 号：P631.3[天文地球—地质矿产勘探]