摘要 本文采用耦合势有限体积法建立了一套快速计算海洋移动导线源电磁响应的三维数值模拟技术。基于电场矢势标势分解将Maxwell方程转换为关于矢势和标势的亥姆霍兹方程；采用洛伦兹规范整理得到对称形式的离散系数矩阵。对导线源采用直接离散方法，基于ILUT （Incomplete LU Preconditioner with Threshold）预处理的BICGSTAB （Biconjugate Gradient Stabilized）迭代算法实现离散线性方程的求解。为了减少移动源总的计算时间，提出采用解的初值重置技术加快正演模拟的速度，即利用发射源移动前的电磁场响应作为发射源移动后的离散方程的初始解。数值计算结果表明：对导线源采用直接离散方法，当导线源的长度大于源周围网格长度的4倍以上时，源附近离散误差显著减小；相比常规的直接迭代算法，采用解的初值重置技术能够减少迭代次数、提高正演计算速度。
Abstract：We propose a marine controlled-source (moving-wire-source) electromagnetic (CSEM) 3D response simulation with a coupled potential finite volume method.First Maxwell's equations are transformed into coupled scalar-vector potential Helmholtz equations with the electric field scalar-vector potential decomposition.A discrete linear equation with symmetry coefficient matrix is obtained with Lorentz gauge.Then the direct discrete method is used for the source.The linear equation of marine CSEM is solved with BICGSTAB iteration after ILUT pre-conditioning.In order to speed up the moving-source iteration calculation and forward modeling,we develop a solution initial value reset technique,namely the response before source moving is used as a solution initial value for the discrete equation after source moving.Based on numerical model tests,the discrete error around the source is significantly reduced when the length of source is four times of the source-mesh length;and the forward modeling with initial value reset is faster than conventional forward modeling.All this proves the effectivity of the proposed approach.
Constable S.Ten years of marine CSEM for hydrocarbon exploration.Geophysics,2010,75(5):A67-A81.
[2]
何展翔,孙卫斌,孔繁恕等.海洋电磁法.石油地球物理勘探,2006,41(4):451-457.He Zhanxiang,Sun Weibin,Kong Fanshu et al.Marine electromagnetic approach.OGP,2006,41(4):451-457.
[3]
Constable S and Srnka L.An introduction to marine controlled-source electromagnetic methods for hydrocarbon exploration.Geophysics,2007,72(2):WA3-WA12.
[4]
Newman G A and Alumbaugh D L.Three-dimensio-nal induction logging problems,Part 2:A finite-diffe-rence solution.Geophysics,2002,67(2):484-491.
[5]
Puzyrev V.A parallel finite-element method for three-dimensional controlled-source electromagnetic forward modelling.Geophysical Journal International,2003,193(2):678-693.
[6]
Plessix R E,Darnet M and Mulder W A.An approach for 3D multisource,multi frequency CSEM modeling.Geophysics,2007,72(5):SM177-SM184.
[7]
殷长春,贲放,刘云鹤等.三维任意各向异性介质中海洋可控源电磁法正演研究.地球物理学报,2014,57(12):4110-4122.Yin Changchun,Ben Fang,Liu Yunhe et al.MCSEM 3D modeling for arbitrarily anisotropic media.Chinese Journal of Geophysics,2014,57(12):4110-4122.
严波,李予国,韩波等.任意方位电偶源的MCSEM电磁场三维正演.石油地球物理勘探,2017,52(4):859-868.Yan Bo,Li Yuguo,Han Bo et al.3D marine controlled-source electromagnetic forward modeling with arbitrarily orientated dipole source.OGP,2017,52(4):859-868.
[10]
朱成,李桐林,杨海斌等.带地形频率域可控源电磁法三维反演研究.石油地球物理勘探,2016,51(5):1031-1039.Zhu Cheng,Li Tonglin,Yang Haibin et al.3D controlled source electromagnetic inversion with topography in the frequency domain.OGP,2016,51(5):1031-1039.
[11]
杨波,徐义贤,何展翔等.考虑海底地形的三维频率域可控源电磁响应有限体积法模拟.地球物理学报,2012,55(4):1390-1399.Yang Bo,Xu Yixian,He Zhanxiang et al.3D frequency-domain modeling of marine controlled source electromagnetic responses with topography using finite volume method.Chinese Journal of Geophysics,2012,55(4):1390-1399.
[12]
杨军,刘颖,吴小平.海洋可控源电磁三维非结构矢量有限元数值模拟.地球物理学报,2015,58(8):2827-2838.Yang Jun,Liu Ying,Wu Xiaoping.3D simulation of marine CSEM using vector finite element method on unstructured grids.Chinese Journal of Geophysics,2015,58(8):2827-2838.
[13]
蔡红柱,熊彬,Michael Zhdanov.电导率各向异性的海洋电磁三维有限单元法正演.地球物理学报,2015,58(8):2839-2850.Cai Hongzhu,Xiong Bin,Zhdanov M.Three-dimensional marine controlled-source electromagnetic mode-lling in anisotropic medium using finite element method.Chinese Journal of Geophysics,2015,58(8):2839-2850.
[14]
韩波,胡祥云,黄一凡等.基于并行化直接解法的频率域可控源电磁三维正演.地球物理学报,2015,58(8):2812-2826.Han Bo,Hu Xiangyun,Huang Yifan et al.3-D frequency-domain CSEM modeling using a parallel direct solver.Chinese Journal of Geophysics,2015,58(8):2812-2826.
[15]
周建美,张烨,汪宏年等.耦合势有限体积法高效模拟各向异性地层中海洋可控源的三维电磁响应.物理学报,2014,63(15):159101-1.Zhou Jianmei,Zhang Ye,Wang Hongnian et al.Efficient simulation of three-dimensional marine controlled-source electromagnetic response in anisotropic formation by means of coupled potential finite volume method.Acta Physica Sinica,2014,63(15):159101-1.
[16]
Streich R.3D finite-difference frequency-domain mode-ling of controlled-source electromagnetic data:Direct solution and optimization for high accuracy.Geophy-sics,2009,74(5):F95.
[17]
Haber E,Ascher U,Aruliah D et al.Fast simulation of 3D electromagnetic problems using potentials.Journal of Computational Physics,2000,163(1):150-171.
[18]
Weiss C J. Project APhiD:a Lorenz-gauged A-decomposition for parallelized computation of ultra-broadband electromagnetic induction in a fully heterogeneous Earth.Computers & Geosciences,2013,58(1):40-52.
[19]
Key K.1D inversion of multi component,multi frequency marine CSEM data:Methodology and synthetic studies for resolving thin resistive layers.Geophysics,2009,74(2):F9-F20.
[20]
Streich R and Becken M.Electromagnetic fields ge-nerated by finite-length wire sources:comparison with point dipole solutions.Geophysical Prospecting,2011,59(2):361-374.
