基于谱元法的频率域三维海洋可控源电磁正演模拟

doi:10.6038/cjg2018L0308

摘要
图/表
参考文献
相关文章 (15)

全文: PDF (2220 KB) HTML (1 KB)
输出: BibTeX | EndNote (RIS)

摘要高精度、快速有效的正演模拟算法是三维电磁正反演的前提.为了提高海洋电磁三维数值模拟的精度和效率，本文提出利用基于Gauss-Lobatto-Chebyshev（GLC）基函数的谱元法进行海洋可控源三维电磁正演模拟.谱元法结合有限元法和谱方法的优点.我们通过应用伽辽金加权残差法离散二次电场矢量亥姆赫兹方程，在单元内选择混合阶GLC多项式的张量积作为高阶矢量插值基函数，在求解大型稀疏线性方程组时利用直接求解器进行快速求解，从而实现了三维海洋可控源电磁快速高精度正演模拟.一维和三维模型正演结果验证了本文算法的有效性和准确性.典型模型的数值结果表明谱元法是一种有效的三维海洋可控源电磁正演数值方法，能在稀疏网格剖分情况下获得精确的海洋电磁正演模拟响应.

	服务

	把本文推荐给朋友
	加入我的书架
	加入引用管理器
	E-mail Alert
	RSS
	作者相关文章
	刘玲
	殷长春
	刘云鹤
	邱长凯
	黄鑫
	张博

关键词 ：谱元法, 海洋可控源电磁, GLC多项式, 三维正演

Abstract：Marine controlled-source electromagnetic (MCSEM) method is an important predrill reservoir appraisal method to reduce exploration risk in detecting sub-seafloor hydrocarbon reservoirs. Most 3D forward modelings for MCSEM are based on conventional numerical methods like finite-difference and finite-element method. In this paper, we introduce spectral element method (SEM) based on Gauss-Lobatto-Chebyshev (GLC) polynomials for 3D frequency-domain MCSEM.
SEM is an accurate and efficient electromagnetic modeling methods due to its spectral accuracy and exponential convergence. The method combines the flexibility of finite-element method with the accuracy of spectral method by a simple application of spectral method to each element. Staring from the Maxwell's equations, we obtain a vector Helmholtz equation of electric field. Then we use Galerkin weighted residual method to discretize the vector Helmholtz equation, in which the curl-conforming Gauss-Lobatto-Chebyshev (GLC) polynomials are chosen as basis functions. As a kind of high-order complete orthogonal polynomials, the GLC polynomials have the characteristic of exponential convergence with the order and can derive the matrix elements analytically, which improves the modeling accuracy. Finally, we use the direct solver MUMPS to solve the final system of equations to speed up the modeling process. For numerical experiments, we first simulate responses of a 3-layer model and compare our SEM results with open-source software to check the accuracy of our SEM algorithm. After that, we analyze efficiency of SEM by comparing with finite-difference method. Finally, we apply SEM for typical 3D models and analyze the MCSEM responses. Numerical results show that SEM is an efficient and effective method for MSCEM modeling, and it can deliver very accurate results even with coarse meshes.

Key words： Spectral element method Marine controlled-source electromagnetic (MCSEM) GLC polynomials 3D modeling

收稿日期: 2017-05-13

PACS:

P631

基金资助:国家重点研发计划重点专项（2016YFC0303100，2017YFC0601903），国家自然科学基金重点项目（41530320），国家自然科学基金面上项目（41274121）和国家青年基金项目（41404093）联合资助.

通讯作者: 殷长春,男,1965年生,教授,国家"千人计划"特聘专家,主要从事电磁勘探理论,特别是航空和海洋电磁方面的研究.E-mail:yinchangchun@jlu.edu.cn E-mail: yinchangchun@jlu.edu.cn

作者简介: 刘玲,女,1993年生,硕士研究生,主要从事地球物理三维电磁正反演理论和方法技术研究.E-mail:liu_ling60@hotmail.com

链接本文:

http://manu39.magtech.com.cn/Geophy/CN/10.6038/cjg2018L0308 或 http://manu39.magtech.com.cn/Geophy/CN/Y2018/V61/I2/756

Amestoy P R, Duff I S, L'Excellent J Y, et al. 2001. A fully asynchronous multifrontal solver using distributed dynamic scheduling. SIAM Journal on Matrix Analysis and Applications, 23(1):15-41, doi:10.1137/S0895479899358194.
Avdeev D B. 2005. Three-dimensional electromagnetic modelling and inversion from theory to application. Surveys in Geophysics, 26(6):767-799, doi:10.1007/s10712-005-1836-x.
Börner R U. 2010. Numerical modelling in geo-electromagnetics:advances and challenges. Surveys in Geophysics, 31(2):225-245, doi:10.1007/s10712-009-9087-x.
Chen G B, Wang H N, Yao J J, et al. 2009. Modeling of electromagnetic responses of a 3-D electrical anomalous body in a layered anisotropic earth using integral equations. Chinese J. Geophys. (in Chinese), 52(8):2174-2181, doi:10.3969/j.issn.0001-5733.2009.08.028.
Constable S. 2010. Ten years of marine CSEM for hydrocarbon exploration. Geophysics, 75(5):75A67-75A81, doi:10.1190/1.3483451.
Ellingsrud S, Eidesmo T, Johansen S, et al. 2002. Remote sensing of hydrocarbon layers by seabed logging (SBL):Results from a cruise offshore Angola. The Leading Edge,21(10):972-982, doi:10.1190/1.1518433.
Han B, Hu X Y, Huang Y F, et al. 2015. 3-D frequency-domain CSEM modeling using a parallel direct solver. Chinese J. Geophys. (in Chinese), 58(8):2812-2826, doi:10.6038/cjg20150816.
Hesthaven J S, Warburton T. 2002. Nodal high-order methods on unstructured grids:I. Time-domain solution of Maxwell's equations. Journal of Computational Physics,181(1):186-221, doi:10.1006/jcph.2002.7118.
Key K. 2009. 1D inversion of multicomponent, multifrequency marine CSEM data:Methodology and synthetic studies for resolving thin resistive layers. Geophysics, 74(2):F9-F20, doi:10.1190/1.3058434.
Key K. 2012. Marine electromagnetic studies of seafloor resources and tectonics. Surveys in Geophysics, 33(1):135-167, doi:10.1007/s10712-011-9139-x.
Komatitsch D, Tromp J. 1999. Introduction to the spectral-element method for 3-D seismic wave propagation. Geophysical journal international, 139:806-822, doi:10.1046/j.1365-246x.1999.00967.x.
Lee J H, Xiao T, Liu Q H. 2006. A 3-D spectral-element method using mixed-order curl conforming vector basis functions for electromagnetic fields. IEEE Transactions on Microwave Theory and Techniques, 54(1):437-444, doi:10.1109/TMTT.2005.860502.
Li L, Liu T, Hu T Y. 2014. Spectral element method with triangular mesh and its application in seismic modeling. Chinese J. Geophys. (in Chinese), 57(4):1224-1234, doi:10.6038/cjg20140419.
Li Y, Key K. 2007. 2D marine controlled-source electromagnetic modeling:Part 1-An adaptive finite-element algorithm. Geophysics, 72(2):WA51-WA62, doi:10.1190/1.2432262.
Liu Y S, Teng J W, Xu T, et al. 2014. Numerical modeling of seismic wavefield with the SEM based on Triangles. Progress in Geophysics (in Chinese), 29(4):1715-1726, doi:10.6038/pg20140430.
Maday Y, Patera A T, Ronquist E M. 1987. A well-posed optimal spectral element approximation for the Stokes problem.
Patera A T. 1984. A spectral element method for fluid dynamics:laminar flow in a channel expansion. Journal of Computational Physics, 54(3):468-488, doi:10.1016/0021-9991(84)90128-1.
Ren Q, Tobón L E, Sun Q, et al. 2015. A new 3-D nonspurious discontinuous Galerkin spectral element time-domain (DG-SETD) method for Maxwell's equations. IEEE Transactions on Antennas and Propagation, 63(6):2585-2594, doi:10.1109/TAP.2015.2417891.
Runge C. 1901. Vber empirische Funktionen und die Interpolation zwischen äquidistanten Ordinaten. Zeitschrift für Mathematik und Physik, 46:224-243.
Seriani G, Priolo E. 1991. High-order spectral element method for acoustic wave modeling. SEG Technical Program Expanded Abstracts 1991. Society of Exploration Geophysicists,1561-1564, doi:10.1190/1.1888989.
Sun X Y, Nie Z P, Zhao Y W, et al. 2008. The electromagnetic modeling of logging-while-drilling tool in tilted anisotropic formations using vector finite element method. Chinese J. Geophys. (in Chinese), 51(5):1600-1607.
Tape C, Liu Q, Maggi A, et al. 2010. Seismic tomography of the southern California crust based on spectral-element and adjoint methods. Geophysical Journal International, 180(1):433-462, doi:10.1111/j.1365-246X.2009.04429.x.
Wang X M, Seriani G, Lin W J. 2007. Some theoretical aspects of elastic wave modeling with a recently developed spectral element method. Science in China Series G:Physics Mechanics and Astronomy, 50(2):185-207, doi:10.1007/s11433-007-0022-1.
Yang J, Liu Y, Wu X P. 2015. 3D simulation of marine CSEM using vector finite element method on unstructured grid. Chinese J. Geophys. (in Chinese), 58(8):2827-2838, doi:10.6038/cjg20150817.
Yin C C, Ben F, Liu Y H, et al. 2014. MCSEM 3D modeling for arbitrarily anisotropic media. Chinese J. Geophys. (in Chinese), 57(12):4110-4122, doi:10.6038/cjg20141222.
Yin C, Huang X, Liu Y, et al. 2017. 3-D Modeling for Airborne EM using the Spectral-element Method. Journal of Environmental and Engineering Geophysics,22(1):13-23, doi:10.2113/JEEG22.1.13.
Zhdanov M S, Lee S K, Yoshioka K. 2006. Integral equation method for 3D modeling of electromagnetic fields in complex structures with inhomogeneous background conductivity. Geophysics, 71(6):G333-G345, doi:10.1190/1.2358403.
陈桂波, 汪宏年, 姚敬金等. 2009. 用积分方程法模拟各向异性地层中三维电性异常体的电磁响应. 地球物理学报, 52(8):2174-2181, doi:10.3969/j.issn.0001-5733.2009.08.028.
韩波, 胡祥云, 黄一凡等. 2015. 基于并行化直接解法的频率域可控源电磁三维正演. 地球物理学报, 58(8):2812-2826, doi:10.6038/cjg20150816.
李琳, 刘韬, 胡天跃. 2014. 三角谱元法及其在地震正演模拟中的应用. 地球物理学报, 57(4):1224-1234, doi:10.6038/cjg20140419.
刘有山, 滕吉文, 徐涛等. 2014. 三角网格谱元法地震波场数值模拟. 地球物理学进展, 29(4):1715-1726, doi:10.6038/pg20140430.
王秀明, 林伟军. 2007. 利用谱元法计算弹性波场的若干理论问题. 中国科学:G 辑, 37(1):41-59.
杨军, 刘颖, 吴小平. 2015. 海洋可控源电磁三维非结构矢量有限元数值模拟. 地球物理学报, 58(8):2827-2838, doi:10.6038/cjg20150817.
殷长春, 贲放, 刘云鹤等. 2014. 三维任意各向异性介质中海洋可控源电磁法正演研究. 地球物理学报, 57(12):4110-4122, doi:10.6038/cjg20141222.