基于交叉梯度耦合的大地电磁与地震走时资料三维联合反演
|
摘要
为了有效解决目前大地电磁和地震走时资料单方法反演结果一致性不好的问题,同时克服基于岩石不同物性参数间关系耦合约束联合反演的局限性,本文研究了基于交叉梯度耦合约束的大地电磁与地震走时资料的三维联合反演算法.以较为成熟的天然地震走时资料三维正反演和大地电磁三维正反演算法为基础,实现了具有共同的反演网格,以交叉梯度结构耦合约束,并能同时获得电阻率和速度模型的三维联合反演算法.分别利用单棱柱体模型和双棱柱体模型合成数据进行了联合反演试算.结果表明:无论是单棱柱体模型还是双棱柱体模型,联合反演结果比单独反演对异常体的空间形态都有更好的恢复,其中单棱柱体模型反演的异常体电阻率更接近于真实电阻率,双棱柱体模型的联合反演结果不仅消除了围岩的部分电阻率假异常,而且增强了对异常体深部速度结构特征的恢复程度.联合反演还能同时改善电阻率和速度反向变化异常体的单独反演结果,进一步证明交叉梯度耦合不依赖于岩石物性关系,而强调地下结构的相似性,具有更普遍的适用性.
To effectively overcome the weak consistency of single inversion of magnetotelluric(or teleseismic travel time data),and also to solve the limitations of petrophysical approach,we study a three-dimensional joint inversion algorithm of magnetotelluric and seismic travel time data with cross-gradient constraints.Based on well-developed three-dimensional teleseismic and magnetotelluric forward and inversion algorithms,a three-dimensional joint inversion algorithm having common inversion grids,with cross-gradient as structural constraint,and with both resistivity and velocity models able to be simultaneously obtained has been achieved.With synthetic data from a single-prism model and two double-prism models,trial computation of joint inversion has been carried out,and the results indicate that,compared with single inversion,joint inversion has better effect in recovering both single prism's and double-prism's spatial morphology.The resistivity value of the single-prism is much closer to the true model in joint inversion,while by double-prism models,joint inversion not only removes false resistivity anomaly of wall rock to a certain degree,but also has large improvement on increasing the resolution of P-wave deep structure of anomalous body.Moreover,joint inversion improve the separate inversion results when resistivity and velocity is incompatible,indicating the cross-gradient is not be restricted on the relationship of different petrophysical parameters,but rather focuses on structural similarity and thus has wide applicability.
To effectively overcome the weak consistency of single inversion of magnetotelluric(or teleseismic travel time data),and also to solve the limitations of petrophysical approach,we study a three-dimensional joint inversion algorithm of magnetotelluric and seismic travel time data with cross-gradient constraints.Based on well-developed three-dimensional teleseismic and magnetotelluric forward and inversion algorithms,a three-dimensional joint inversion algorithm having common inversion grids,with cross-gradient as structural constraint,and with both resistivity and velocity models able to be simultaneously obtained has been achieved.With synthetic data from a single-prism model and two double-prism models,trial computation of joint inversion has been carried out,and the results indicate that,compared with single inversion,joint inversion has better effect in recovering both single prism's and double-prism's spatial morphology.The resistivity value of the single-prism is much closer to the true model in joint inversion,while by double-prism models,joint inversion not only removes false resistivity anomaly of wall rock to a certain degree,but also has large improvement on increasing the resolution of P-wave deep structure of anomalous body.Moreover,joint inversion improve the separate inversion results when resistivity and velocity is incompatible,indicating the cross-gradient is not be restricted on the relationship of different petrophysical parameters,but rather focuses on structural similarity and thus has wide applicability.
引文
[1]Aric K,Adam A,Smythe D K.Combined seismic andmagnetotelluric imaging of upper crystalline crust in theSouthern Bohemian Massif.First Break,1997,15(8):265-271.
[2]Favetto A,Pomposiello C,Booker J,et al.Magnetotelluricinversion constrained by seismic data in the Tucumán Basin(Andean Foothills,27°S,NW Argentina).J.Geophys.Res.,2007,112(B9),B09104.
[3]Muoz,Bauer K,Moeck I,et al.Exploring the GroSchnebeck(Germany)geothermal site using a statisticaljoint interpretation of magnetotelluric and seismic tomographymodels.Geothermics,2010,39(1):35-45.
[4]高德章,赵金海,薄玉玲等.东海重磁地震综合探测剖面研究.地球物理学报,2004,47(5):853-861.Gao D Z,Zhao J H,Bo Y L,et al.A profile study ofgravitative-magnetic and seismic comprehensive survey in theEast China Sea.Chinese J.Geophys.(in Chinese),2004,47(5):853-861.
[5]胡平,刘保金,白立新等.奥林匹克公园地区隐伏断裂综合探测.地球物理学报,2010,53(6):1486-1494.Hu P,Liu B J,Bai L X,et al.Synthetic exploration of theburied faults in Olympic park area.Chinese J.Geophys.(inChinese),2010,53(6):1486-1494.
[6]李德春,杨书江,胡祖志等.三维重磁电震资料的联合解释———以库车大北地区山前砾石层为例.石油地球物理勘探,2012,47(2):353-359.Li D C,Yang S J,Hu Z Z,et al.Integrated interpretation of3Dgravity,magnetic,electromagnetic and seismic data:acase study of conglomerate mass investigation in piedmontarea of Kuche Depression.Oil Geophysical Prospecting(inChinese),2012,47(2):353-359.
[7]姜枚,谭捍东,钱辉等.金川铜镍矿床的地球物理深部结构与成因模式.矿床地质,2012,31(2):207-215.Jiang M,Tan H D,Qian H,et al.Geophysical deepstructure and genetic model of Jinchuan copper-nickeldeposit.Mineral Deposits(in Chinese),2012,31(2):207-215.
[8]姜枚,彭淼,王有学等.喜马拉雅东构造结岩石圈板片深俯冲的地球物理证据.岩石学报,2012,28(6):1755-1764.Jiang M,Peng M,Wang Y X,et al.Geophysical evidencefor deep subduction of Indian lithospheric plate beneathEastern Himalayan Syntaxis.Acta Petrologica Sinica(inChinese),2012,28(6):1755-1764.
[9]Moorkamp M,Jones A G,Eaton D W.Joint inversion ofteleseismic receiver functions and magnetotelluric data using agenetic algorithm:Are seismic velocities and electricalconductivities compatible?Geophys.Res.Lett.,2007,34(16):L16311.
[10]Moorkamp M,Jones A G,Fishwick S.Joint inversion ofreceiver functions,surface wave dispersion,and magnetotelluricdata.J.Geophys.Res.,2010,115(B4):B04318.
[11]Roux E,Moorkamp M,Jones A G,et al.Joint inversion oflong-period magnetotelluric data and surface-wave dispersioncurves for anisotropic structure:Application to data fromCentral Germany.Geophys.Res.Lett.,2011,38(5):L05304.
[12]彭淼,谭捍东,姜枚等.利用接收函数和大地电磁数据联合反演南迦巴瓦构造结中部地区壳幔结构.地球物理学报,2012,55(7):2281-2291.Peng M,Tan H D,Jiang M,et al.Joint inversion of receiverfunctions and magnetotelluric data:Application to crustal andmantle structure beneath central Namche Barwa,easternHimalayan syntaxis.Chinese J.Geophys.(in Chinese),2012,55(7):2281-2291.
[13]Kozlovskaya E,Hjelt S E.Modeling of elastic and electricalproperties of solid-liquid rock system with fractal microstructure.Physics and Chemistry of the Earth,Part A:Solid Earth andGeodesy,2000,25(2):195-200.
[14]Heincke B,Jegen M,Hobbs R.Joint inversion of MT,gravity and seismic data applied to sub-basalt imaging.//SEG Technical Program Expanded Abstracts,2006,25(1):784-789.
[15]Colombo D,Stefano M D.Geophysical modeling viasimultaneous joint inversion of seismic,gravity,andelectromagnetic data:application to prestack depth imaging.The Leading Edge,2007,26(3):326-331.
[16]Vermeesch P M,Morgan J V,Christeson G L,et al.Three-dimensional joint inversion of traveltime and gravity dataacross the Chicxulub impact crater.J.Geophys.Res.,2009,114(6):B02105.
[17]杨辉,王家林,吴健生等.大地电磁与地震资料仿真退火约束联合反演.地球物理学报,2002,45(5):723-734.Yang H,Wang J L,Wu J S,et al.Constrained jointinversion of magnetotelluric and seismic data using simulatedannealing algorithm.Chinese J.Geophys.(in Chinese),2002,45(5):723-734.
[18]于鹏,戴明刚,王家林等.电阻率和速度随机分布的MT与地震联合反演.地球物理学报,2009,52(4):1089-1097.Yu P,Dai M G,Wang J L,et al.Joint inversion ofmagnetotelluric and seismic data based on random resistivityand velocity distributions.Chinese J.Geophys.(in Chinese),2009,52(4):1089-1097.
[19]陈高,于鹏,陈晓等.改进的大地电磁与地震资料联合反演方法在黔中隆起区的适用性研究.石油物探,2010,49(2):158-163,197.Chen G,Yu P,Chen X,et al.Adaptability research onimproved joint inversion method of magnetotelluric andseismic data in middle Guizhou Uplift.Geophys.Pros.forPetroleum.(in Chinese),2010,49(2):158-163,197.
[20]陈晓,于鹏,张罗磊等.地震与大地电磁测深数据的自适应正则化同步联合反演.地球物理学报,2011,54(10):2673-2681.Chen X,Yu P,Zhang L L,et al.Adaptive regularizedsynchronous joint inversion of MT and seismic data.ChineseJ.Geophys.(in Chinese),2011,54(10):2673-2681.
[21]陈华根,李嘉虓,吴健生等.MT-重力模拟退火联合反演研究.地球物理学报,2012,55(2):663-670.Chen H G,Li J X,Wu J S,et al.Study on simulated-annealing MT-gravity joint inversion.Chinese J.Geophys.(in Chinese),2012,55(2):663-670.
[22]何委徽,王家林,于鹏.地球物理联合反演研究的现状与趋势分析.地球物理学进展,2009,24(2):530-540.He W H,Wang J L,Yu P.Overview of the status andprospect of geophysical joint inversion.Progress in Geophysics.(in Chinese),2009,24(2):530-540.
[23]刘彦,吕庆田,孟贵祥等.大地电磁与地震联合反演研究现状与展望.地球物理学进展,2012,27(6):2444-2451.Liu Y,Lu Q T,Meng G X,et al.Joint electromagnetic andseismic inversion survey:status and prospect.Progress inGeophysics(in Chinese),2012,27(6):2444-2451.
[24]Meju M A,Gallardo L A,Mohamed A K.Evidence forcorrelation of electrical resistivity and seismic velocity inheterogeneous near-surface materials.Geophys.Res.Lett.,2003,30(7):1373-1376.
[25]Haber E,Oldenburg D.Joint inversion:a structuralapproach.Inverse Problems,1997,13(1):63-77.
[26]Gallardo L A,Meju M A.Characterization of heterogeneousnear-surface materials by joint 2Dinversion of DC resistivityand seismic data.Geophys.Res.Lett.,2003,30(13):1658.
[27]Gallardo L A,Meju M A.Joint two-dimensional cross-gradient imaging of magnetotelluric and seismic traveltimedata for structural and lithological classification.Geophys.J.Int.,2007,169(3):1261-1272.
[28]Linde N,Binley A,Tryggvason A,et al.Improvedhydrogeophysical characterization using joint inversion ofcross-hole electrical resistance and ground-penetrating radartraveltime data.Water Resources Research,2006,42(12):W04410.
[29]Hu W Y,Abubakar A,Habasshy T M.Jointelectromagnetic and seismic inversion using structuralconstraints.Geophysics,2009,74(6):R99-R109.
[30]Moorkamp M,Heincke B,Jegen M,et al.A framework for3-D joint inversion of MT,gravity and seismic refractiondata.Geophys.J.Int.,2011,184(1):477-493.
[31]谭捍东,余钦范,Booker J等.大地电磁法三维交错采样有限差分数值模拟.地球物理学报,2003,46(5):705-711.Tan H D,Yu Q F,Booker J,et al.Magnetotelluric three-dimensional modeling using the staggered-grid finite differencemethod.Chinese J.Geophys.(in Chinese),2003,46(5):705-711.
[32]Siripunvaraporn W,Egbert G.An efficient data-subspaceinversion method for 2-D magnetotelluric data.Geophysics,2000,65(3):791-803.
[33]Siripunvaraporn W,Egbert G,Lenbury Y,et al.Three-dimensional magnetotelluric inversion:data-space method.Phys.Earth Planet.Interiors,2005,150(1-3):3-14.
[34]Sethian J A,Popovici A M.3-D traveltime computation usingthe fast marching method.Geophysics,1999,64(2):516-523.
[35]De Kool M,Rawlinson N,Sambridge M.A practical grid-based method for tracking multiple refraction and reflectionphases in three-dimensional heterogeneous media.Geophys.J.Int.,2006,167(1):253-270.
[36]Kennett B L N,Sambridge M S,Williamson P R.Subspacemethods for large inverse problems with multiple parameterclasses.Geophys.J.Int.,1988,94(2):237-247.
[37]Sambridge M S.Non-linear arrival time inversion:constraining velocity anomalies by seeking smooth models in3-D.Geophys.J.Int.,1990,102(3):653-677.
[38]Rawlinson N,Reading A M,Kennett B L N.Lithosphericstructure of Tasmania from a novel form of teleseismictomography.J.Geophys.Res.,2006,111(B2):B02301.
[2]Favetto A,Pomposiello C,Booker J,et al.Magnetotelluricinversion constrained by seismic data in the Tucumán Basin(Andean Foothills,27°S,NW Argentina).J.Geophys.Res.,2007,112(B9),B09104.
[3]Muoz,Bauer K,Moeck I,et al.Exploring the GroSchnebeck(Germany)geothermal site using a statisticaljoint interpretation of magnetotelluric and seismic tomographymodels.Geothermics,2010,39(1):35-45.
[4]高德章,赵金海,薄玉玲等.东海重磁地震综合探测剖面研究.地球物理学报,2004,47(5):853-861.Gao D Z,Zhao J H,Bo Y L,et al.A profile study ofgravitative-magnetic and seismic comprehensive survey in theEast China Sea.Chinese J.Geophys.(in Chinese),2004,47(5):853-861.
[5]胡平,刘保金,白立新等.奥林匹克公园地区隐伏断裂综合探测.地球物理学报,2010,53(6):1486-1494.Hu P,Liu B J,Bai L X,et al.Synthetic exploration of theburied faults in Olympic park area.Chinese J.Geophys.(inChinese),2010,53(6):1486-1494.
[6]李德春,杨书江,胡祖志等.三维重磁电震资料的联合解释———以库车大北地区山前砾石层为例.石油地球物理勘探,2012,47(2):353-359.Li D C,Yang S J,Hu Z Z,et al.Integrated interpretation of3Dgravity,magnetic,electromagnetic and seismic data:acase study of conglomerate mass investigation in piedmontarea of Kuche Depression.Oil Geophysical Prospecting(inChinese),2012,47(2):353-359.
[7]姜枚,谭捍东,钱辉等.金川铜镍矿床的地球物理深部结构与成因模式.矿床地质,2012,31(2):207-215.Jiang M,Tan H D,Qian H,et al.Geophysical deepstructure and genetic model of Jinchuan copper-nickeldeposit.Mineral Deposits(in Chinese),2012,31(2):207-215.
[8]姜枚,彭淼,王有学等.喜马拉雅东构造结岩石圈板片深俯冲的地球物理证据.岩石学报,2012,28(6):1755-1764.Jiang M,Peng M,Wang Y X,et al.Geophysical evidencefor deep subduction of Indian lithospheric plate beneathEastern Himalayan Syntaxis.Acta Petrologica Sinica(inChinese),2012,28(6):1755-1764.
[9]Moorkamp M,Jones A G,Eaton D W.Joint inversion ofteleseismic receiver functions and magnetotelluric data using agenetic algorithm:Are seismic velocities and electricalconductivities compatible?Geophys.Res.Lett.,2007,34(16):L16311.
[10]Moorkamp M,Jones A G,Fishwick S.Joint inversion ofreceiver functions,surface wave dispersion,and magnetotelluricdata.J.Geophys.Res.,2010,115(B4):B04318.
[11]Roux E,Moorkamp M,Jones A G,et al.Joint inversion oflong-period magnetotelluric data and surface-wave dispersioncurves for anisotropic structure:Application to data fromCentral Germany.Geophys.Res.Lett.,2011,38(5):L05304.
[12]彭淼,谭捍东,姜枚等.利用接收函数和大地电磁数据联合反演南迦巴瓦构造结中部地区壳幔结构.地球物理学报,2012,55(7):2281-2291.Peng M,Tan H D,Jiang M,et al.Joint inversion of receiverfunctions and magnetotelluric data:Application to crustal andmantle structure beneath central Namche Barwa,easternHimalayan syntaxis.Chinese J.Geophys.(in Chinese),2012,55(7):2281-2291.
[13]Kozlovskaya E,Hjelt S E.Modeling of elastic and electricalproperties of solid-liquid rock system with fractal microstructure.Physics and Chemistry of the Earth,Part A:Solid Earth andGeodesy,2000,25(2):195-200.
[14]Heincke B,Jegen M,Hobbs R.Joint inversion of MT,gravity and seismic data applied to sub-basalt imaging.//SEG Technical Program Expanded Abstracts,2006,25(1):784-789.
[15]Colombo D,Stefano M D.Geophysical modeling viasimultaneous joint inversion of seismic,gravity,andelectromagnetic data:application to prestack depth imaging.The Leading Edge,2007,26(3):326-331.
[16]Vermeesch P M,Morgan J V,Christeson G L,et al.Three-dimensional joint inversion of traveltime and gravity dataacross the Chicxulub impact crater.J.Geophys.Res.,2009,114(6):B02105.
[17]杨辉,王家林,吴健生等.大地电磁与地震资料仿真退火约束联合反演.地球物理学报,2002,45(5):723-734.Yang H,Wang J L,Wu J S,et al.Constrained jointinversion of magnetotelluric and seismic data using simulatedannealing algorithm.Chinese J.Geophys.(in Chinese),2002,45(5):723-734.
[18]于鹏,戴明刚,王家林等.电阻率和速度随机分布的MT与地震联合反演.地球物理学报,2009,52(4):1089-1097.Yu P,Dai M G,Wang J L,et al.Joint inversion ofmagnetotelluric and seismic data based on random resistivityand velocity distributions.Chinese J.Geophys.(in Chinese),2009,52(4):1089-1097.
[19]陈高,于鹏,陈晓等.改进的大地电磁与地震资料联合反演方法在黔中隆起区的适用性研究.石油物探,2010,49(2):158-163,197.Chen G,Yu P,Chen X,et al.Adaptability research onimproved joint inversion method of magnetotelluric andseismic data in middle Guizhou Uplift.Geophys.Pros.forPetroleum.(in Chinese),2010,49(2):158-163,197.
[20]陈晓,于鹏,张罗磊等.地震与大地电磁测深数据的自适应正则化同步联合反演.地球物理学报,2011,54(10):2673-2681.Chen X,Yu P,Zhang L L,et al.Adaptive regularizedsynchronous joint inversion of MT and seismic data.ChineseJ.Geophys.(in Chinese),2011,54(10):2673-2681.
[21]陈华根,李嘉虓,吴健生等.MT-重力模拟退火联合反演研究.地球物理学报,2012,55(2):663-670.Chen H G,Li J X,Wu J S,et al.Study on simulated-annealing MT-gravity joint inversion.Chinese J.Geophys.(in Chinese),2012,55(2):663-670.
[22]何委徽,王家林,于鹏.地球物理联合反演研究的现状与趋势分析.地球物理学进展,2009,24(2):530-540.He W H,Wang J L,Yu P.Overview of the status andprospect of geophysical joint inversion.Progress in Geophysics.(in Chinese),2009,24(2):530-540.
[23]刘彦,吕庆田,孟贵祥等.大地电磁与地震联合反演研究现状与展望.地球物理学进展,2012,27(6):2444-2451.Liu Y,Lu Q T,Meng G X,et al.Joint electromagnetic andseismic inversion survey:status and prospect.Progress inGeophysics(in Chinese),2012,27(6):2444-2451.
[24]Meju M A,Gallardo L A,Mohamed A K.Evidence forcorrelation of electrical resistivity and seismic velocity inheterogeneous near-surface materials.Geophys.Res.Lett.,2003,30(7):1373-1376.
[25]Haber E,Oldenburg D.Joint inversion:a structuralapproach.Inverse Problems,1997,13(1):63-77.
[26]Gallardo L A,Meju M A.Characterization of heterogeneousnear-surface materials by joint 2Dinversion of DC resistivityand seismic data.Geophys.Res.Lett.,2003,30(13):1658.
[27]Gallardo L A,Meju M A.Joint two-dimensional cross-gradient imaging of magnetotelluric and seismic traveltimedata for structural and lithological classification.Geophys.J.Int.,2007,169(3):1261-1272.
[28]Linde N,Binley A,Tryggvason A,et al.Improvedhydrogeophysical characterization using joint inversion ofcross-hole electrical resistance and ground-penetrating radartraveltime data.Water Resources Research,2006,42(12):W04410.
[29]Hu W Y,Abubakar A,Habasshy T M.Jointelectromagnetic and seismic inversion using structuralconstraints.Geophysics,2009,74(6):R99-R109.
[30]Moorkamp M,Heincke B,Jegen M,et al.A framework for3-D joint inversion of MT,gravity and seismic refractiondata.Geophys.J.Int.,2011,184(1):477-493.
[31]谭捍东,余钦范,Booker J等.大地电磁法三维交错采样有限差分数值模拟.地球物理学报,2003,46(5):705-711.Tan H D,Yu Q F,Booker J,et al.Magnetotelluric three-dimensional modeling using the staggered-grid finite differencemethod.Chinese J.Geophys.(in Chinese),2003,46(5):705-711.
[32]Siripunvaraporn W,Egbert G.An efficient data-subspaceinversion method for 2-D magnetotelluric data.Geophysics,2000,65(3):791-803.
[33]Siripunvaraporn W,Egbert G,Lenbury Y,et al.Three-dimensional magnetotelluric inversion:data-space method.Phys.Earth Planet.Interiors,2005,150(1-3):3-14.
[34]Sethian J A,Popovici A M.3-D traveltime computation usingthe fast marching method.Geophysics,1999,64(2):516-523.
[35]De Kool M,Rawlinson N,Sambridge M.A practical grid-based method for tracking multiple refraction and reflectionphases in three-dimensional heterogeneous media.Geophys.J.Int.,2006,167(1):253-270.
[36]Kennett B L N,Sambridge M S,Williamson P R.Subspacemethods for large inverse problems with multiple parameterclasses.Geophys.J.Int.,1988,94(2):237-247.
[37]Sambridge M S.Non-linear arrival time inversion:constraining velocity anomalies by seeking smooth models in3-D.Geophys.J.Int.,1990,102(3):653-677.
[38]Rawlinson N,Reading A M,Kennett B L N.Lithosphericstructure of Tasmania from a novel form of teleseismictomography.J.Geophys.Res.,2006,111(B2):B02301.
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心