Modeling and inversion of the CSEM field in 2.5D anisotropic media with an adaptive finite element method
Su Qi1,2, Dai Shikun1,2, Zhao Dongdong1,2
1. School of Geosciences and Info-Physics, Central South University, Changsha, Hu'nan 410083, China; 2. Key Laboratory of Metallogenic Prediction of Nonferrous Metals and Geological Environment Monitoring, Ministry of Education, Central South University, Changsha, Hu'nan 410083, China
Abstract:We conduct a forward calculation of anisotropic media with the adaptive finite element method based on unstructured triangular meshes.In the simulation of tough topography and complex geological structures, the adaptive mesh strategy can ensure not only computational efficiency but also computational accuracy.The improved Occam algorithm based on Gauss Newton method can find quickly and appropriately the Lagrange multiplier under the condition of little dependence on the initial conditions.The forward modelings of a 1D VTI model and a classical 2D anisotropic ocean reservoir model prove the computational accuracy and efficiency.For a complex 2D anisotropic seafloor model, the inversion algorithm shows resistivity and buried depth of surrounding rocks, bedrocks and anisotropic anomalous bodies.Finally, tests on measured data shows the feasibility of the proposed algorithm.
粟琪, 戴世坤, 赵东东. 可控源电磁法2.5维自适应有限元各向异性正反演[J]. 石油地球物理勘探, 2018, 53(2): 418-429.
Su Qi, Dai Shikun, Zhao Dongdong. Modeling and inversion of the CSEM field in 2.5D anisotropic media with an adaptive finite element method. Oil Geophysical Prospecting, 2018, 53(2): 418-429.
Coggon J H.Electromagnetic and electrical modeling by the finite element method.Geophysics,1971,36(2):132-155.
[2]
Stoyer C H and Greenfield R J.Numerical solutions of the response of a two-dimensional earth to an oscillating magnetic dipole source.Geophysics,1976,41(3):519-530.
[3]
Doherty J.EM modeling using surface integral equations.Geophysics,1988,53(6):644-668.
[4]
Unsworth J M,Bryan J T,Alan D C.Electromagnetic induction by a finite electric dipole source over a 2-D earth.Geophysics,1993,58(2):198-214.
[5]
Torres V C,Habashy T M.Rapid 2.5D forward mode-ling and inversion via a new nonlinear scattering approximation.Radio Science,1994,29(3):1051-1079.
[6]
Mistsuhata Y.2-D electromagnetic modeling by finite-element method with a dipole source and topography.Geophysics,2000,65(2):465-475.
[7]
Li Y.A finite-element algorithm for electromagnetic induction in two dimensional anisotropic conductivity structures.Geophysical Journal International,2002,148(3):389-401.
Li Y,Constable S.2D marine controlled-source electro-magnetic modeling (Part 2):The effect of bathymetry.Geophysics,2007,72(2):WA63-WA71.
[10]
柳建新,汤文武,童孝忠.基于双二次插值的2.5维FCSEM有限元正演模拟.中南大学学报(自然科学版),2014,45(2):474-482.Liu Jianxin,Tang Wenwu,Tong Xiaozhong.Finite element forward modeling of 2.5-D FCSEM based on biquadratic interpolation.Journal of Central South University (Science and Technology),2014(2):474-482.
[11]
Unsworth M,Oldenburg D.Subspace inversion of electromagnetic data:Application to mid-ocean-ridge exploration.Geophysical Journal International,1995,123(1):161-168.
[12]
Abubakar A,Habashy T M,Druskin V L et al.2.5D forward and inverse modeling for interpreting low-frequency electromagneticmeasurements.Geophysics,2008,73(4):F165-F177.
[13]
Key K,Ovall J.A parallel goal-oriented adaptive finite element method for 2.5-D electromagnetic modelling.Geophysical Journal International,2011,186(1):137-154.
[14]
Key K.Marine EM inversion using unstructured grids:a 2D parallel adaptive finite element algorithm.SEG Technical Program Expanded Abstracts,2012,31:1-5.
[15]
Myer D,Key K,Constable S.Marine CSEM of the Scarborough gas field,Part 2:2D inversion.Geophy-sics,2015,80(3):E187-E196.
[16]
陈光源,杜立彬,景建恩等.2.5维海洋可控源电磁反演算法及影响参数研究.地球物理学进展,2016,31(4):1796-1802.Chen Guangyuan,Du Libin,Jing Jian'en et al.2.5D inversion algorithm and parameters of marine controlled-source electromagnetic method.Progress in Geophysics,2016,31(4):1796-1802.
[17]
Li X and Pedersen L B.The electromagnetic response of an azimuthally anisotropic half-space.Geophysics,1991,56(9):1462-1473.
[18]
Li X and Pedersen L B.Controlled-source tensor magnetotelluric responses of a layered earth with azimu-thal anisotropy.Geophysical Journal International,1992,111(1):91-103.
[19]
Loseth L O,Ursin B.Electromagnetic fields in planarly layered anisotropic media.Geophysical Journal International,2007,170(1):44-80.
[20]
罗鸣,李予国.一维电阻率各向异性对海洋可控源电磁响应的影响研究.地球物理学报,2015,58(8):2851-2861.Luo Ming,Li Yuguo.Effects of the electric anisotropy on marine controlled-source electromagnetic responses.Chinese Journal of Geophysics,2015,58(8):2851-2861.
[21]
罗鸣,李予国,李刚.一维垂直各向异性介质频率域海洋可控源电磁资料反演方法.地球物理学报,2016,59(11):4349-4359.Luo Ming,Li Yuguo,Li Gang.Frequency domain inversion of marine CSEM data in one dimensional vertically anisotropic structures.Chinese Journal of Geophysics,2016,59(11):4349-4359.
[22]
Kong F N,Johnstad S E,Rosten T et al.A 2.5D finite element-modeling difference method for marine CSEM modeling in stratified anisotropic media.Geophysics,2008,73(1):F9-F19.
[23]
Li Y G,Dai S K.Finite element modelling of marine controlled-source electromagnetic responses in two-dimensional dipping anisotropic conductivity structures.Geophysical Journal International,2011,185(2):622-636.
[24]
刘颖,李予国.层状各向异性介质中任意取向电偶源的海洋电磁响应.石油地球物理勘探,2015,50(4):755-765.Liu Ying,Li Yuguo.Marine controlled-source electromagnetic fields of an arbitrary electric dipole over a layered anisotropic medium.OGP,2015,50(4):755-765.
[25]
赵东东.各向异性介质可控源电磁法2.5维数值模拟及观测系统研究[学位论文].湖南长沙:中南大学,2016.Zhao Dongdong.Research on 2.5D Modeling and Observation System for Controlled-Source Electromagnetic Method in Anisotropic Medium[D].Central South University,Changsha,Hu'nan,2016.
[26]
Kittur M G,Huston R L,Oswald F B.Mesh refinement in finite element analysis by minimization of the stiffness matrix trace.Computers & Structures,1989,31(6):891-896.
[27]
Shewchuk J R.Reprint of delaunay refinement algorithms for triangular mesh generation.Computational Geometry,2014,47(7):741-778.
[28]
刘颖.海洋可控源电磁法二维有限元正演及反演[学位论文].山东青岛:中国海洋大学,2014.Liu Ying.2D Finite Element Modeling and Inversion for Marine Controlled-Source Electromagnetic Fields[D].Chinese Marine University,Qingdao,Shandong,2014.
[29]
McGillivray P R,Oldenburg D W,Ellis R G et al.Calculation of sensitivities for the frequency-domain electromagnetic problem.Geophysical Journal International,1994,116(1):1-4.
[30]
Farquharson C G,Oldenburg D W.Approximate sensi-tivities for the electromagnetic inverse problem.Geophysical Journal International,1996,126(1):235-252.
[31]
Key K.MARE2DEM:a 2-D inversion code for con-trolled-source electromagnetic and magnetotelluric data.Geophysical Journal International,2016,207(1):571-588.
[32]
Constable S C,Parker R L,Constable C G.Occam'sinversion-A practical algorithm for generating smooth models from electromagnetic sounding data.Geophysics,1987,52(1):289-300.
[33]
Weitemeyer K,Constable S,Key K et al.First results from a marine controlled-source electromagnetic survey to detect gas hydrates offshore Oregon.Geophysical Research Letters,2006,33(3):L03304.
[34]
Weitemeyer K,Gao G,Constable S et al.The practical application of 2D inversion to marine controlled-source electromagnetic data.Geophysics,2010,75(6):F199-F211.