基于重启多项式Krylov子空间模型降阶的全波形瞬变电磁三维正演  

Three-dimensional full-time transient electromagnetic modeling using restarting polynomial Krylov subspace model order reduction method

作　　者：周建美[1,2,3] 景旭 鲁凯亮 刘文韬 李貅[1,2,3] ZHOU JianMei;JING Xu;LU KaiLiang;LIU WenTao;LI Xiu(College of Geology Engineering and Geomatics,Chang′an University,Xi′an 710054,China;Integrated Geophysical Simulation Lab of Chang′an University(Key Laboratory of Chinese Geophysical Society),Xi′an 710054,China;National Engineering Research Center of Offshore Oil and Gas Exploration,Beijing100028,China)

机构地区：[1]长安大学地质工程与测绘学院,西安710054 [2]长安大学地球物理场多参数综合模拟实验室(中国地球物理学会重点实验室),西安710054 [3]海洋油气勘探国家工程研究中心,北京100028

出　　处：《地球物理学报》2023年第6期2646-2657,共12页Chinese Journal of Geophysics

基　　金：国家自然科学基金项目(42274092,41830101);陕西省自然科学基金项目(2022JQ-285);长安大学中央高校基本科研业务费专项资金(300102261201);甘肃煤田地质局科研项目(2017119);甘肃省自然资源厅科技创新项目(202238)联合资助。

摘　　要：瞬变电磁三维正演在内存占用和求解时间上都是大的挑战.瞬变电磁响应受到发射波形的影响,正演计算中还需考虑发射波形.基于重启多项式Krylov子空间模型降阶方法,本文实现了一种新的全波形瞬变电磁三维正演算法.采用规则六面体网格的拟态有限体积方法对全波形控制方程进行空间离散,引入一阶φ函数推导得到全时段响应的统一表达式,利用重启多项式Krylov子空间算法求解φ函数表达式,对于给定的重启子空间维度,基于残差公式得到满足给定精度的任意时刻的正演响应.该算法内存占用小,不需要求解大型线性方程组,内存占用主要由空间离散网格数量和重启子空间维度决定.数值模型正演结果验证了本文算法的有效性.Transient electromagnetic(TEM)modeling is a challenge in terms of memory and computing time.TEM responses are affected by the full transmitting-current waveforms and needs to be considered in the forward calculation.In this paper,we present a restarting polynomial Krylov method for modeling three-dimensional full-time TEM responses.The mimetic finite volume method is carried out for spatial discretization.The full-time TEM response then can be expressed as a matrixφfunction.The restarting polynomial Krylov method is used to solve the matrixφfunction.For a given restart subspace dimension,the residual is used to obtain the forward response at any time that meet the given accuracy.This method does not need to solve large-scale linear equations.Furthermore,the memory usage is mainly determined by the number of spatial discrete grids and the restart subspace order.The numerical experiments demonstrated that the method is accurate and memory efficient.

关 键 词：瞬变电磁 三维正演 发射波形 模型降阶 

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