Progress on joint inversion of electromagnetic constraints combined with multiple methods and multiple information
ZHANG Shuoning1,2,3, YU Guo1,2,3, HE Zhanxiang1,2,3
1. Shenzhen Key Laboratory of Deep Offshore Oil and Gas Exploration Technology, Shenzhen, Guangdong 518055, China; 2. Guangdong Provincial Key Laboratory of Geophysical High-resolution Imaging Technology, Shenzhen, Guangdong 518055, China; 3. Department of Earth and Space Sciences, Southern University of Science and Technology, Shenzhen, Guangdong 518055, China
Abstract:Electromagnetic exploration plays an important role in oil and gas exploration,deep mineral resources survey,geothermal groundwater exploration,imaging of deep earth structures,and underwater military target reconnaissance.This paper systematically summarizes the research progress on joint inversion of electromagnetic constraints combined with multiple methods and information,introduces the research status from the perspective of method application,and focuses on the construction form,existing problems,and technical difficulties of the objective function for different inversion methods.In areas with more known information such as seismic data,different types of constraint functions are constructed by utilizing priori information to improve the tectonic imaging accuracy in saline and igneous rock development areas and reduce the inversion non-uniqueness.In areas where the magnetotelluric (MT) and controlled source electromagnetic (CSEM) data have been collected,the two methods can be employed to carry out joint inversion with their respective advantages.Various forms such as joint inversion with sequential constraints and joint inversion based on a unified objective function are developed.Additionally,CSEM and MT make up for each other in-depth inversion to improve the target inversion accuracy and reduce non-uniqueness.This paper discusses the selection of objective functions and weighting factors for this type of joint inversion.According to actual application effects,the most important and potential joint inversion is the inversion of electromagnetic and seismic data.Joint inversion has gone through from seismic constrained electromagnetic inversion to physical model-based and cross-gradient method-based joint inversion.This paper also analyzes the difficulty in constructing a reasonable connection between velocity and resistivity models and explains that the weighting factor setting of velocity and resistivity models in the objective functions is still the key and difficult points.Finally,this paper proposes to combine electromagnetic joint inversion,artificial intelligence algorithms,and three-dimensional visualization analysis.It also points out that large-scale parallel computation and other new methods and technologies can be adopted in the future,and joint acquisition systems can be developed to provide more effective data volumes for the multi-physical field joint inversion.
何展翔.电磁勘探技术的机遇与挑战及发展方向[J].物探化探计算技术,2019,41(4): 433-447.HE Zhanxiang.Opportunities,challenges and development directions of electromagnetic exploration today[J].Computing Techniques for Geophysical and Geochemical Exploration,2019,41(4): 433-447.
[2]
赵国泽,何展翔,魏文博.地球电磁学科的发展现状和展望[M].北京: 中国科学技术出版社,2009.
[3]
赵国泽,陈小斌,汤吉.中国地球电磁法新进展和发展趋势[J].地球物理学进展,2007,22(4): 1171-1180.ZHAO Guoze,CHEN Xiaobin,TANG Ji.Advanced geo-electromagnetic methods in China[J].Progress in Geophysics,2007,22(4): 1171-1180.
[4]
CAGNIARD L.Basic theory of the magnetotelluric method of geophysical prospecting[J].Geophysics,1953,18(3): 605-635.
[5]
SANDBERG S K,HOHMANN G W.Controlled-source audiomagnetotellurics in geothermal exploration[J].Geophysics,1982,47(1): 100-116.
[6]
薛国强,李貅,底青云.瞬变电磁法理论与应用研究进展[J].地球物理学进展,2007,22(4): 1195-1200.XUE Guoqiang,LI Xiu,DI Qingyun.The progress of TEM in theory and application[J].Progress in Geophysics,2007,22(4): 1195-1200.
何展翔,胡祖志,王志刚,等.时频电磁(TFEM)技术: 数据联合约束反演[J].石油地球物理勘探,2020,55(4): 898-905.HE Zhanxiang,HU Zuzhi,WANG Zhigang,et al.Time-frequency electromagnetic (TFEM) technique: step-by-step constraint inversion based on artificial fish swarm algorithm[J].Oil Geophysical Prospecting,2020,55(4): 898-905.
[9]
敬荣中.地球物理非线性联合反演方法研究[D].湖南长沙: 中南大学,2002.
[10]
杨文采.地球物理反演的理论与方法[M].北京: 地质出版社,1997.
[11]
王家映.地球物理反演理论[M].湖北武汉: 中国地质大学出版社,1998.
[12]
厉波.测井约束反演中利用声波测井进行精确时深转换[J].内蒙古石油化工,2005,31(11): 114-115.LI Bo.Accurate time-depth transformation using sonic logging in log-constraint inversion[J].Inner Mongolia Petrochemical Industry,2005,31(11): 114-115.
[13]
谢万学,王永刚,邢文军.井孔与地面地震数据联合反演方法综述[J].勘探地球物理进展,2006,29(1): 1-6.XIE Wanxue,WANG Yonggang,XING Wenjun.Review of joint inversion methods for borehole and surface seismic data[J].Progress in Exploration Geophysics,2006,29(1): 1-6.
[14]
聂凯轩,陆正元,冯斌.测井约束反演技术在LZ21-4构造带储盖层综合分析中的应用[J].石油地球物理勘探,2007,42(1): 98-101.NIE Kaixuan,LU Zhengyuan,FENG Bin.Application of log-constraint inversion technique in comprehensive analysis of reservoir and cap in LZ21-4 structural belt[J].Oil Geophysical Prospecting,2007,42(1): 98-101.
[15]
谭荣彪,李瑞,李珂泂,等.井—震约束反演在薄互层复杂构造地区的应用[J].物探化探计算技术,2008,30(1): 52-57.TAN Rongbiao,LI Rui,LI Kejiong,et al.The Application of logging-seismic constrained inversion in thin interbed complex structure area[J].Computing Techniques for Geophysical and Geochemical Exploration,2008,30(1): 52-57.
[16]
连太炜,樊洪海,于玲玲.测井约束地震反演在地层孔隙压力预测中的应用[J].西部探矿工程,2009,21(6): 64-65.LIAN Taiwei,FAN Honghai,YU Lingling.Application of the formation pore pressure prediction by log-constraint seismic inverison[J].West-China Exploration Engineering,2009,21(6): 64-65.
杨辉,王永涛,戴世坤,等.带地形的MT多参量二维快速模拟退火约束反演[J].石油地球物理勘探,2003,38(2): 213-217.YANG Hui,WANG Yongtao,DAI Shikun,et al.2-D inversion of magnetotelluric (MT) multi-parameters with topography using fast simulated annealing[J].Oil Geophysical Prospecting,2003,38(2): 213-217.
王有涛,何展翔,陈学国,等.基于物性分析建模的时频电磁反演及储层评价——以准北缘石北构造带为例[J].石油地球物理勘探,2022,57(6): 1489-1497.WANG Youtao,HE Zhanxiang,CHEN Xueguo,et al.Time-frequency electromagnetic data inversion and reservoir evaluation based on physical property analysis modeling: a case study of the Shibei structural belt on the northern margin of the Junggar Basin[J].Oil Geophysical Prospecting,2022,57(6): 1489-1497.
[24]
石艳玲,黄文辉,魏强,等.电磁井震约束反演识别川中深层裂谷[J].石油地球物理勘探,2016,51(6): 1233-1240.SHI Yanling,HUANG Wenhui,WEI Qiang,et al.Deep rift identification with MT inversion constrained by shallow logging and seismic data[J].Oil Geophysical Prospecting,2016,51(6): 1233-1240.
[25]
王传武,李术才,刘斌,等.基于参考模型的三维电阻率约束反演方法研究[J].岩土工程学报,2016,38(9): 1685-1694.WANG Chuanwu,LI Shucai,LIU Bin,et al.3D constrained electrical resistivity inversion method based on reference model[J].Chinese Journal of Geotechnical Engineering,2016,38(9): 1685-1694.
[26]
刘斌,聂利超,李术才,等.三维电阻率空间结构约束反演成像方法[J].岩石力学与工程学报,2012,31(11): 2258-2268.LIU Bin,NIE Lichao,LI Shucai,et al.3D electrical resistivity inversion tomography with spatial structural constraint[J].Chinese Journal of Rock Mechanics and Engineering,2012,31(11): 2258-2268.
[27]
DA SILVA PINHEIRO J L,SILVA M W C,DA PIEDADE A A,et al.3D joint inversion of MMT and MCSEM data[C].17th International Congress of the Brazilian Geophysical Society,Rio de Janeiro,Brazil,16-19 August 2021,1-7.
[28]
RODI W,MACKIE R L.Nonlinear conjugate gradients algorithm for 2-D magnetotelluric inversion[J].Geophysics,2001,66(1): 174-187.
[29]
COMMER M,NEWMAN G A.Three-dimensional controlled-source electromagnetic and magnetotelluric joint inversion[J].Geophysical Journal International,2009,178(3): 1305-1316.
[30]
MACKIE R,DONALD WATTS M,RODI W.Joint 3D inversion of marine CSEM and MT data[C].SEG Technical Program Expanded Abstracts,2007,26: 574-578.
[31]
GRIBENKO A V,ZHDANOV M S.Joint 3D inversion of marine CSEM and MT data[C].SEG Technical Program Expanded Abstracts,2011,30: 552-556.
[32]
KUMAR K,BORGOHAIN D,MORTEN J P,et al.Offshore sub-basalt exploration using CSEM and MT[C].International Petroleum Technology Conference,Kuala Lumpur,Malaysia,2014,IPTC-17917-MS.
[33]
SASAKI Y.3D inversion of marine CSEM and MT data: an approach to shallow-water problem[J].Geophysics,2013,78(1): E59-E65.
[34]
艾正敏,叶益信,汤文武,等.基于非结构三角网格的海洋CSEM和MT二维联合反演研究[J].物探与化探,2021,45(1): 149-158.AI Zhengmin,YE Yixin,TANG Wenwu,et al.A test study of 2D joint inversion of marine CSEM and MT based on unstructured triangular grid[J].Geophysical and Geochemical Exploration,2021,45(1): 149-158.
[35]
MEQBEL N,RITTER O.Joint 3D inversion of multiple electromagnetic datasets[J].Geophysical Prospecting,2015,63(6): 1450-1467.
[36]
LIAO W,PENG R,HU X,et al.3-D joint inversion of MT and CSEM data for imaging a high-temperature geothermal system in Yanggao region,Shanxi province,China[J].IEEE Transactions on Geoscience and Remote Sensing,2022,60: 1-13.
[37]
ABUBAKAR A,LI M,LIU J,et al.Simultaneous joint inversion of MT and CSEM data using a multiplicative cost function[C].SEG Technical Program Expanded Abstracts,2009,28: 719-723.
[38]
赵宁,王堃鹏,秦策.海洋CSEM和MT一维联合反演研究[J].物探化探计算技术,2014,36(2): 148-151.ZHAO Ning,WANG Kunpeng,QIN Ce.Research on 1-D joint inversion of marine CSEM and MT[J].Computing Techniques for Geophysical and Geochemical Exploration,2014,36(2): 148-151.
DE STEFANO M,GOLFRE ANDREASI F,RE S,et al.Multiple-domain,simultaneous joint inversion of geophysical data with application to subsalt imaging[J].Geophysics,2011,76(3): R69-R80.
[41]
GAO G,ABUBAKAR A,HABASHY T M.Joint petrophysical inversion of electromagnetic and full-waveform seismic data[J].Geophysics,2012,77(3): WA3-WA18.
HOVERSTEN G M,CASSASSUCE F,GASPERIKOVA E,et al.Direct reservoir parameter estimation using joint inversion of marine seismic AVA and CSEM data[J].Geophysics,2006,71(3): C1-C13.
[44]
ABUBAKAR A,GAO G,HABASHY T M,et al.Joint inversion approaches for geophysical electromagnetic and elastic full-waveform data[J].Inverse Problems,2012,28(5): 055016.
[45]
DU Z,MACGREGOR L M.Reservoir characterization from joint inversion of marine CSEM and seismic AVA data using genetic algorithms: a case study based on the Luva gas field[C].SEG Technical Program Expanded Abstracts,2010,29: 737-741.
[46]
徐凯军,杜润林,刘展.海洋可控源电磁与地震一维联合储层参数反演[J].石油地球物理勘探,2016,51(1): 197-203.XU Kaijun,DU Runlin,LIU Zhan.Joint reservoir parameter inversion of 1D marine controlled source electromagnetic and seismic data[J].Oil Geophysical Prospecting,2016,51(1): 197-203.
[47]
彭国民,徐凯军,杜润林,等.MCSEM和地震AVA数据联合反演储层物性[J].石油地球物理勘探,2018,53(5): 1110-1116.PENG Guomin,XU Kaijun,DU Runlin,et al.Reservoir petrophysical parameter estimation with joint inversion of MCSEM and seismic AVA data[J].Oil Geophysical Prospecting,2018,53(5): 1110-1116.
[48]
CHEN J,HOVERSTEN G M,VASCO D,et al.A Bayesian model for gas saturation estimation using marine seismic AVA and CSEM data[J].Geophysics,2007,72(2): WA85-WA95.
[49]
CHEN J,HOVERSTEN G M.Joint inversion of marine seismic AVA and CSEM data using statistical rock-physics models and Markov random fields[J].Geophysics,2012,77(1): R65-R80.
[50]
GALLARDO L A,MEJU M A.Joint two‐dimensional DC resistivity and seismic travel time inversion with cross-gradients constraints[J].Journal of Geophysical Research (Solid Earth),2004,109(B3): B03311.
[51]
彭淼,谭捍东,姜枚,等.基于交叉梯度耦合的大地电磁与地震走时资料三维联合反演[J].地球物理学报,2013,56(8): 2728-2738.PENG Miao,TAN Handong,JIANG Mei,et al.Three-dimensional joint inversion of magnetotelluric and seismic travel time data with cross-gradient constraints[J].Chinese Journal of Geophysics,2013,56(8): 2728-2738.
[52]
齐嘉慧.基于交叉梯度结构约束的大地电磁和地震全波形数据交替联合反演[D].吉林长春: 吉林大学,2020.QI Jiahui.Joint Iinversion of Alternating Magnetotelluric and Full-waveform Seismic Data Based on Cross-gradient Structure Constraints[D].Jilin University,Changchun,Jilin,2020.
GALLARDO L A,MEJU M A.Joint two-dimensional cross-gradient imaging of magnetotelluric and seismic traveltime data for structural and lithological classification[J].Geophysical Journal International,2007,169(3): 1261-1272.
[56]
GALLARDO L A.Multiple cross-gradient joint inversion for geospectral imaging[J].Geophysical Research Letters,2007,34(19): L19301.
[57]
GALLARDO L A,MEJU M A.Characterization of heterogeneous near‐surface materials by joint 2D inversion of DC resistivity and seismic data[J].Geophysical Research Letters,2003,30(13): 1658.
[58]
高级,张海江.基于交叉梯度交替结构约束的二维地震走时与全通道直流电阻率联合反演[J].地球物理学报,2016,59(11): 4310-4322.GAO Ji,ZHANG Haijiang.Two-dimensional joint inversion of seismic velocity and electrical resistivity using seismic travel times and full channel electrical measurements based on alternating cross-gradient structural constraint[J].Chinese Journal of Geophysics,2016,59(11): 4310-4322.
[59]
陈晓,于鹏,邓居智,等.地球物理联合反演新框架研究[J].石油地球物理勘探,2017,52(4): 851-858,883.CHEN Xiao,YU Peng,DENG Juzhi,et al.A new framework for geophysical joint inversion[J].Oil Geophysical Prospecting,2017,52(4): 851-858,883.
[60]
ISAEV I,OBORNEV E,OBORNEV I,et al.Neural network recognition of the type of parameterization scheme for magnetotelluric data[C].Advances in Neural Computation,Machine Learning,and Cognitive Research II,Springer,Cham,2018,176-183.
[61]
CONWAY D,ALEXANDER B,KING M,et al.Inverting magnetotelluric responses in a three-dimensional earth using fast forward approximations based on artificial neural networks[J].Computers & Geosciences,2019,127: 44-52.
[62]
LIU W,XI Z,WANG H,et al.Two-dimensional deep learning inversion of magnetotelluric sounding data[J].Journal of Geophysics and Engineering,2021,18(5): 627-641.
[63]
ALYOUSUF T,YAOGUO L.Inversion using adaptive physics-based neural network: application to magnetotelluric inversion[C].International Petroleum Technology Conference,Riyadh,Saudi Arabia,2022,IPTC-22504-EA.
[64]
SUN J,LI Y.Joint inversion of multiple geophysical data using guided fuzzy C-means clustering[J].Geophysics,2016,81(3): ID37-ID57.
郑元满,姚长利,张晨.重磁三维自动反演软件系统的分析与设计[J].现代地质,2012,26(6): 1225-1230.ZHENG Yuanman,YAO Changli,ZHANG Chen.Design a software for 3-D automatic gravity and magnetic inversion[J].Geoscience,2012,26(6): 1225-1230.
杨大方,翁爱华,杨悦,等.一个CSAMT非平面波效应的实例及数值分析[J].地质与勘探,2015,51(1): 151-156.YANG Dafang,WENG Aihua,YANG Yue,et al.An example of non-plane wave effect in CSAMT and numerical analysis[J].Geology and Exploration,2015,51(1): 151-156.