大地电磁、重力、磁法和地震初至波走时的交叉梯度二维联合反演研究
|
摘要
为了准确的探测和描绘地下复杂的地质结构,同时克服地球物理单一方法反演的多解性和单一参数反演模型的不一致性等问题,近年来基于交叉梯度联合反演的综合地球物理解释已经得到了广泛的关注和应用.本文首先研究了两种地球物理方法的交叉梯度联合反演算法,在此基础上,推导并实现了多种地球物理方法(大地电磁,重力,磁法,地震初至波走时)的多交叉梯度约束的二维联合反演算法;其次,我们设计了结构不一致模型和复杂模型,针对多物性联合反演算法的准确性和有效性进行了模拟试算,并对复杂模型的单独反演结果和联合反演结果进行了交叉梯度值和物性交会图的对比;最后,本文将成熟的卫星资料多光谱综合分析技术应用到联合反演中,将多物性参数反演模型结果图通过RGB(红-绿-蓝)模式进行合成,得到融合的RGB合成图.结果表明:通过对结构不一致模型和复杂模型的联合反演结果和单独反演结果的对比分析,可以得出联合反演得到的结果更接近真实模型,并从得到的交叉梯度值进一步证明了联合反演模型相似度高,也从物性交会图中得到联合反演的物性相关性更好的结论,反向证明了算法的正确性.最终从得到的RGB合成图像,我们可以更直观的分析反演结果,更有利于准确划分地下模型结构.
To accurately detect complex geological structures in the subsurface and solve theproblems in geophysical inversion such as multiple solutions and inconsistency of the single parameter inversion model,the comprehensive geophysical interpretation based on the crossgradient joint inversion has been widely concerned and applied in recent years.In this paper,we first study the cross-gradient joint inversion algorithm of two geophysical methods.Then,we derive and implement the 2 Djoint inversion algorithm of multiple cross-gradients constraints of various geophysical methods(magnetotellurics(MT),gravity,magnetic,and seismic firstarrival traveltime).Next,we design models with inconsistent structure and a complex structure,respectively,and carry out simulation calculations focusing on the accuracy and effectiveness of the multiple joint inversion algorithm.Meanwhile,the results of separate inversion and joint inversion of complex model are compared with the cross-gradient value and cross-plots of pairs of different physical properties.Finally,this work applies the mature satellite data multi-spectral comprehensive analysis technology to the joint inversion,and synthesizes the multiple parameter inversion model results through RGB(red-green-blue)mode to obtain the fused RGB composite image.The results show that whether in the physical value or in the spatial geometry,the result of the joint inversion is closer to the real model.Furthermore,we also prove that this joint inversion model has a higher structure similarity and physical correlation by comparing the crossgradient value and crossplots of pairs of different physical properties.From the obtained RGB composite images,we can analyze the inversion results more intuitively,which permits to accurately classify models of underground models.
To accurately detect complex geological structures in the subsurface and solve theproblems in geophysical inversion such as multiple solutions and inconsistency of the single parameter inversion model,the comprehensive geophysical interpretation based on the crossgradient joint inversion has been widely concerned and applied in recent years.In this paper,we first study the cross-gradient joint inversion algorithm of two geophysical methods.Then,we derive and implement the 2 Djoint inversion algorithm of multiple cross-gradients constraints of various geophysical methods(magnetotellurics(MT),gravity,magnetic,and seismic firstarrival traveltime).Next,we design models with inconsistent structure and a complex structure,respectively,and carry out simulation calculations focusing on the accuracy and effectiveness of the multiple joint inversion algorithm.Meanwhile,the results of separate inversion and joint inversion of complex model are compared with the cross-gradient value and cross-plots of pairs of different physical properties.Finally,this work applies the mature satellite data multi-spectral comprehensive analysis technology to the joint inversion,and synthesizes the multiple parameter inversion model results through RGB(red-green-blue)mode to obtain the fused RGB composite image.The results show that whether in the physical value or in the spatial geometry,the result of the joint inversion is closer to the real model.Furthermore,we also prove that this joint inversion model has a higher structure similarity and physical correlation by comparing the crossgradient value and crossplots of pairs of different physical properties.From the obtained RGB composite images,we can analyze the inversion results more intuitively,which permits to accurately classify models of underground models.
引文
Bennington N L,Zhang H J,Thurber C H,et al.2015.Joint inversion of seismic and magnetotelluric data in the parkfield region of California using the normalized cross-gradient constraint.Pure and Applied Geophysics,172(5):1033-1052.
Colombo D,M De Stefano.2007.Geophysical modeling via simultaneous joint inversion of seismic,gravity,and electromagnetic data:Application to prestack depth imaging.The Leading Edge,26(3):326-331.
Constable S C,Parker R L,Constable C G.1987.Occam′s Inversion:Apractical algorithm for generating smooth models from electromagnetic sounding data.Geophysics,52(3):289-300.
De Groot-Hedlin C D,Constable S C.1990.Occam′s inversion to generate smooth two-dimensional models from magnetotelluric data.Geophysics,55(12):1613-1624.
Fregoso E,Gallardo L A.2009.Cross-gradients joint 3Dinversion with applications to gravity and magnetic data.Geophysics,74(4):L31-L42.
Gallardo L A,Meju M A.2003.Characterization of heterogeneous near-surface materials by joint 2Dinversion of DC resistivity and seismic data.Geophysical Research Letters,30(13):1658,doi:10.1029/2003GL017370.
Gallardo L A,Meju M A.2004.Joint two-dimensional DC resistivity and seismic travel time inversion with cross-gradients constraints.Journal of Geophysical Research:Solid Earth,109(B3):B03311,doi:10.1029/2003JB002716.
Gallardo L A.2007.Multiple cross-gradient joint inversion for geospectral imaging.Geophysical Research Letters,34(19):L19301,doi:10.1029/2007GL030409.
Gallardo L A,Fontes S L,Meju M A,et al.2012.Robust geophysical integration through structure-coupled joint inversion and multispectral fusion of seismic reflection,magnetotelluric,magnetic,and gravity images:Example from Santos Basin,offshore Brazil.Geophysics,77(5):B237-B251.
Heincke B,M Jegen,R Hobbs.2006.Joint inversion of MT,gravity and seismic data applied to sub-basalt imaging.76th Annual International Meeting.SEG,Expanded Abstracts,784-789.
Jones A G,Fishwick S,Evans R L,Spratt J.E.2009.Correlation of lithospheric velocity and electrical conductivity for Southern Africa.In:11th SAGA biennial technical meeting and exhibition,Swaziland,428-434.
Li T L,Zhang R Z,Piao Y C.2015.Joint Inversion of magnetotelluric and first-arrival wave seismic traveltime with cross-gradient constraints.Journal of Jilin University:Earth Science Edition,45(3),952-961,doi:10.13278/j.cnki.jjuese.201503304.
Li T L,Zhang R Z,Piao Y C,et al.2016.Multiple joint inversion of geophysical data with sub-region crossgradient constraints.Chinese Journal of Geophysics(in Chinese),59(8):2979-2988,doi:10.6038/cjg20160821.
Linde N,Tryggvason A,Peterson J E,et al.2008.Joint inversion of crosshole radar and seismic traveltimes acquired at the South Oyster Bacterial Transport Site.Geophysics,73(4):G29-G37.
Lochbühler T,Doetsch J,Brauchler R,et al.2013.Structurecoupled joint inversion of geophysical and hydrological data.Geophysics,78(3):ID1-ID14.
Moorkamp M,Jones A G,Eaton D W.2007.Joint inversion of teleseismic receiver functions and magnetotelluric data using a genetic algorithm:Are seismic velocities and electrical conductivities compatible?.Geophysical Research Letters,34(16):L16311.
Moorkamp M,Heincke B,Jegen M,et al.2011.A framework for3-D joint inversion of MT,gravity and seismic refraction data.Geophysical Journal International,184(1):477-493.
Moorkamp M,Roberts A W,Jegen M,et al.2013.Verification of velocity-resistivity relationships derived from structural joint inversion with borehole data.Geophysical Research Letters,40(14):3596-3601.
Peng M,Tan H D,Jiang M,et al.2013.Three-dimensional joint inversion of magnetotelluric and seismic travel time data with cross-gradient constraints.Chinese Journal of Geophysics(in Chinese),56(8):2728-2738,doi:10.6038/cjg20130821.
Singh B.2002.Simultaneous computation of gravity and magnetic anomalies resulting from a 2-D object.Geophysics,67(3):801-806.
Tarantola A.1987.Inverse problem theory:methods for data fitting and model parameter estimation.Physics of the Earth&Planetary Interiors,57(3):350-351.
Tillmann A,St9cker T.2000.A new approach for the joint inversion of seismic and geoelectric data.Paper Presented at 63rd EAGEConference and Technical Exhibition.Eur.Assoc.of Geosci.and Eng,Amsterdam.
Vidale J E.1990.Finite-difference calculation of travel-times in three dimensions.Geophysics,55(5):521-526.
Wang J,Meng X H,Chen Z X,et al.2013.The theory of crossgradient and its application in geophysical joint inversion.Progress in Geophysics(in Chinese),28(4):2094-2103,doi:10.6038/pg20130454.
Yang H,Dai S K,Song H B,et al.2002.Overview of joint inversion of integrated geophysics.Progress in Geophysics(in Chinese),17(2):262-271,doi:10.3969/j.issn.1004-2903.2002.02.011.
Yu P,Dai M G,Wang J L,et al.2009.Joint inversion of magnetotelluric and seismic data based on random resistivity and velocity distribution.Chinese Journal of Geophysics(in Chinese),52(4):1089-1097,doi:10.3969/j.issn.0001-5733.2009.04.026.
李桐林,张镕哲,朴英哲等.2015.大地电磁测深与地震初至波走时交叉梯度反演.吉林大学学报:地球科学版,45(3):952-961,doi:10.13278/j.cnki.jjuese.201503304.
李桐林,张镕哲,朴英哲等.2016.部分区域约束下的交叉梯度多重地球物理数据联合反演.地球物理学报,59(8):2979-2988,doi:10.6038/cjg20160821.
彭淼,谭捍东,姜枚等.2013.基于交叉梯度耦合的大地电磁与地震走时资料三维联合反演.地球物理学报,56(8):2728-2738,doi:10.6038/cjg20130821.
王俊,孟小红,陈召曦等.2013.交叉梯度理论及其在地球物理联合反演中的应用.地球物理学进展,28(4):2094-2103,doi:10.6038/pg20130454.
杨辉,戴世坤,宋海滨等.2002.综合地球物理联合反演综述.地球物理学进展,17(2):262-271,doi:10.3969/j.issn.1004-2903.2002.02.011.
于鹏,戴明刚,王家林等.2009.电阻率和速度随机分布的MT与地震联合反演.地球物理学报,52(4):1089-1097.doi:10.3969/j.issn.0001-5733.2009.04.026.
Colombo D,M De Stefano.2007.Geophysical modeling via simultaneous joint inversion of seismic,gravity,and electromagnetic data:Application to prestack depth imaging.The Leading Edge,26(3):326-331.
Constable S C,Parker R L,Constable C G.1987.Occam′s Inversion:Apractical algorithm for generating smooth models from electromagnetic sounding data.Geophysics,52(3):289-300.
De Groot-Hedlin C D,Constable S C.1990.Occam′s inversion to generate smooth two-dimensional models from magnetotelluric data.Geophysics,55(12):1613-1624.
Fregoso E,Gallardo L A.2009.Cross-gradients joint 3Dinversion with applications to gravity and magnetic data.Geophysics,74(4):L31-L42.
Gallardo L A,Meju M A.2003.Characterization of heterogeneous near-surface materials by joint 2Dinversion of DC resistivity and seismic data.Geophysical Research Letters,30(13):1658,doi:10.1029/2003GL017370.
Gallardo L A,Meju M A.2004.Joint two-dimensional DC resistivity and seismic travel time inversion with cross-gradients constraints.Journal of Geophysical Research:Solid Earth,109(B3):B03311,doi:10.1029/2003JB002716.
Gallardo L A.2007.Multiple cross-gradient joint inversion for geospectral imaging.Geophysical Research Letters,34(19):L19301,doi:10.1029/2007GL030409.
Gallardo L A,Fontes S L,Meju M A,et al.2012.Robust geophysical integration through structure-coupled joint inversion and multispectral fusion of seismic reflection,magnetotelluric,magnetic,and gravity images:Example from Santos Basin,offshore Brazil.Geophysics,77(5):B237-B251.
Heincke B,M Jegen,R Hobbs.2006.Joint inversion of MT,gravity and seismic data applied to sub-basalt imaging.76th Annual International Meeting.SEG,Expanded Abstracts,784-789.
Jones A G,Fishwick S,Evans R L,Spratt J.E.2009.Correlation of lithospheric velocity and electrical conductivity for Southern Africa.In:11th SAGA biennial technical meeting and exhibition,Swaziland,428-434.
Li T L,Zhang R Z,Piao Y C.2015.Joint Inversion of magnetotelluric and first-arrival wave seismic traveltime with cross-gradient constraints.Journal of Jilin University:Earth Science Edition,45(3),952-961,doi:10.13278/j.cnki.jjuese.201503304.
Li T L,Zhang R Z,Piao Y C,et al.2016.Multiple joint inversion of geophysical data with sub-region crossgradient constraints.Chinese Journal of Geophysics(in Chinese),59(8):2979-2988,doi:10.6038/cjg20160821.
Linde N,Tryggvason A,Peterson J E,et al.2008.Joint inversion of crosshole radar and seismic traveltimes acquired at the South Oyster Bacterial Transport Site.Geophysics,73(4):G29-G37.
Lochbühler T,Doetsch J,Brauchler R,et al.2013.Structurecoupled joint inversion of geophysical and hydrological data.Geophysics,78(3):ID1-ID14.
Moorkamp M,Jones A G,Eaton D W.2007.Joint inversion of teleseismic receiver functions and magnetotelluric data using a genetic algorithm:Are seismic velocities and electrical conductivities compatible?.Geophysical Research Letters,34(16):L16311.
Moorkamp M,Heincke B,Jegen M,et al.2011.A framework for3-D joint inversion of MT,gravity and seismic refraction data.Geophysical Journal International,184(1):477-493.
Moorkamp M,Roberts A W,Jegen M,et al.2013.Verification of velocity-resistivity relationships derived from structural joint inversion with borehole data.Geophysical Research Letters,40(14):3596-3601.
Peng M,Tan H D,Jiang M,et al.2013.Three-dimensional joint inversion of magnetotelluric and seismic travel time data with cross-gradient constraints.Chinese Journal of Geophysics(in Chinese),56(8):2728-2738,doi:10.6038/cjg20130821.
Singh B.2002.Simultaneous computation of gravity and magnetic anomalies resulting from a 2-D object.Geophysics,67(3):801-806.
Tarantola A.1987.Inverse problem theory:methods for data fitting and model parameter estimation.Physics of the Earth&Planetary Interiors,57(3):350-351.
Tillmann A,St9cker T.2000.A new approach for the joint inversion of seismic and geoelectric data.Paper Presented at 63rd EAGEConference and Technical Exhibition.Eur.Assoc.of Geosci.and Eng,Amsterdam.
Vidale J E.1990.Finite-difference calculation of travel-times in three dimensions.Geophysics,55(5):521-526.
Wang J,Meng X H,Chen Z X,et al.2013.The theory of crossgradient and its application in geophysical joint inversion.Progress in Geophysics(in Chinese),28(4):2094-2103,doi:10.6038/pg20130454.
Yang H,Dai S K,Song H B,et al.2002.Overview of joint inversion of integrated geophysics.Progress in Geophysics(in Chinese),17(2):262-271,doi:10.3969/j.issn.1004-2903.2002.02.011.
Yu P,Dai M G,Wang J L,et al.2009.Joint inversion of magnetotelluric and seismic data based on random resistivity and velocity distribution.Chinese Journal of Geophysics(in Chinese),52(4):1089-1097,doi:10.3969/j.issn.0001-5733.2009.04.026.
李桐林,张镕哲,朴英哲等.2015.大地电磁测深与地震初至波走时交叉梯度反演.吉林大学学报:地球科学版,45(3):952-961,doi:10.13278/j.cnki.jjuese.201503304.
李桐林,张镕哲,朴英哲等.2016.部分区域约束下的交叉梯度多重地球物理数据联合反演.地球物理学报,59(8):2979-2988,doi:10.6038/cjg20160821.
彭淼,谭捍东,姜枚等.2013.基于交叉梯度耦合的大地电磁与地震走时资料三维联合反演.地球物理学报,56(8):2728-2738,doi:10.6038/cjg20130821.
王俊,孟小红,陈召曦等.2013.交叉梯度理论及其在地球物理联合反演中的应用.地球物理学进展,28(4):2094-2103,doi:10.6038/pg20130454.
杨辉,戴世坤,宋海滨等.2002.综合地球物理联合反演综述.地球物理学进展,17(2):262-271,doi:10.3969/j.issn.1004-2903.2002.02.011.
于鹏,戴明刚,王家林等.2009.电阻率和速度随机分布的MT与地震联合反演.地球物理学报,52(4):1089-1097.doi:10.3969/j.issn.0001-5733.2009.04.026.