三维多层介质重力-地震同步联合反演
|
摘要
联合反演是地球物理勘探的重要解释手段,能够提高模型参数的反演精度.本文在归纳和分析重力与地震资料联合反演的研究和应用现状的基础上,利用三维多层介质模型的地震走时和重力正演公式,推导了地震走时和重力异常对界面深度的雅可比矩阵,实现了三维重力-地震同步联合反演界面成像.最后进行了数值理论模型模拟和实例计算,结果表明地震走时和重力同步联合反演很好的重建了三维多层介质界面.
Joint inversion is an important means of geophysical interpretation,which can improve the accuracy of the inversion for the model parameters.This paper summarizes and analyzes the current state of the research and application of joint inversion using gravity and seismic data.Then,we employ the forward formula of seismic travel time and gravity anomaly in a 3D Multi-layer media model to deduce the Jacobian matrix of travel time and gravity anomaly with interface.On the base of these results,we implement the joint inversion of gravity and seismic data in a 3D Multi-layer media model for imaging interface.We perform a numerical simulation and calculation on real data to verify this method.The results show that joint inversion using gravity and seismic data can reconstruct the model interface well.
Joint inversion is an important means of geophysical interpretation,which can improve the accuracy of the inversion for the model parameters.This paper summarizes and analyzes the current state of the research and application of joint inversion using gravity and seismic data.Then,we employ the forward formula of seismic travel time and gravity anomaly in a 3D Multi-layer media model to deduce the Jacobian matrix of travel time and gravity anomaly with interface.On the base of these results,we implement the joint inversion of gravity and seismic data in a 3D Multi-layer media model for imaging interface.We perform a numerical simulation and calculation on real data to verify this method.The results show that joint inversion using gravity and seismic data can reconstruct the model interface well.
引文
[1] 刘光鼎.论综合地球物理解释原则与实例[A].见:地球物理学报编辑委员会编.八十年代中国地球物理学进展[C].北京:学术书刊出版社,1989． 231～242． Liu G D.Principles and applications of integrated geophysical interpretation[A].In:Chinese J.Geophys.Eds.Progress of Chinese Geophysics in 1980s(in Chinese)[C].Beijing:Scientific Books and Periodicals Press,1989． 231～242．
[2] 杨辉,戴世坤,宋海斌,等.综合地球物理联合反演综述[J].地球物理学进展,2002,17(2) :262～271． Yang H,Dai S K,Song H B,et al.Overview of joint inversion of integrated geophysics[J].Progress in Geophysics(in Chinese),2002,17(2) :262～271．
[3] 敬荣中,鲍光淑,陈绍裘.地球物理联合反演研究综述[J].地球物理学进展,2003,18(3) :535～540． Jing R Z,Bao G S,Chen S Q.A review of the researches for geophysical combinative inversion[J].Progress in Geophysics(in Chinese),2003,18(3) :535～540．
[4] 于鹏,王家林,吴健生,等.地球物理联合反演的研究现状和分析[J].勘探地球物理进展,2006,29(2) :87～93． Yu P,Wang J L,Wu J S,et al.Review and discussions on geophysical joint inversion[J].Progress in Exploration Geophysics(in Chinese),2006,29(2) :87～ 93．
[S]Savino J M,Rodi W L,Masso J F.Simultaneous inversion of multiple geophysical data sets for earth structure[J].50~(th)Ann.Internat.Mtg,Soc.Expi.Geophys.Expanded Abstracts,1980,438～439．
[6] Golizdra G Y.Statement of the problem of comprehensive interpretation of gravity fields and seismic observations[J].Izvest.Earth Phys.,1980,16:535～539．
[7] 王一新,廉西京.重力与地震联合解释在石油勘探中的应用[J].石油物探,1983,22(3) :34～41． Wang Y X,Lian X J.Applications of gravity and seismic integrated interpretation in petroleum exploration[J].Geophysical Prospecting for Petroleum(in Chinese),1983,22(3) :34～41．
[8] Lines L,Schultz A K,Treitel S.Cooperative inversion of geophysicaldata[J].57~(th)Ann.Internat.Mtg.,Soc.Expi.Geophys.Expanded Abstracts,1987． 814～816．
[9] Lines L,Schultz A K,Treitel S.Cooperative inversion of geophysical data[J].Geophysics,1988,53(1) :8～20．
[10] 汪宏年,张文生.层状介质中重力、地震联合反演的迭代算法[J].石油地球物理勘探,1993,28(2) :153～165． Wang H N,Zhang W S.Iterative algorithm for joint inversion of gravimetric and seismic data in layered medium[[J].Oil Geophysical Prospecting(in Chinese),1993,28(2) :153～165．
[11] 冯锐,陶裕录.地震重力联合反演中的非块状一致性模型[J].地球物理学报,1993,36(4) :463～475． Feng R,Tao Y L.A non-block consistent model in seismo gravity inversion[J].Chinese J.Geophys.(in Chinese),1993,36(4) :463～475．
[12] 关小平,黄嘉正.谈重震联合反演问题[J].物探与化探,1994,18(6) :429～437． Guan X P,Huang J Z.On the combined gravity seismic inversion problem.Geophysical & Geochemical Exploration(in Chinese),1994,18(6) :429～437．
[13] 范兴才,王家林.松辽盆地北部杏山地区典型的地质地球物理模型及综合解释[A].见:王家林,吴健生主编.中国典型含油气盆地综合地球物理研究[C].上海:同济大学出版社,1995． 112～124． Fan X C,Wang J L.Typical geological and geophysical model and interpretation in Xing shan area of northern Song liao Basin[A].In:Wang J L,Wu J S eds.Integrated Geophysical Researches on Typical Petroleum Basins in China[C](in Chinese).Shanghai:Tongji University Press,19951112～124．
[14] 陈冰,王家林.剥离法与综合模型法的地球物理综合反演[A].见:王家林,吴健生主编.中国典型含油气盆地综合地球物理研究[C].上海:同济大学出版社,1995． 125～135． Chen B,Wang J L.Integrated geophysical inversion by stripping and integrated model methods[A].In:Wang J L,Wu J S eds.Integrated Geophysical Researches on Typical Petroleum Basins in China[C](in Chinese).Shanghai:Tongii University Press,1995． 125～135．
[15] 王西文,孟昭泰,张工会,等.重力、地震联合反演在研究含油气盆地构造中的应用[J].石油地球物理勘探,1991,26(2) :201～209． Wang X W,Meng Z T,Zhang G H,et al.Application of combined gravitational and seismic inversion to study geological structure inside a hydrocarbon bearing basin[J].Oil Geophysical Prospecting(in Chinese),1991,26(2) :201～209．
[16] 王西文.用重力、地震资料联合反演直接预测油气藏的方法[J].石油地球物理勘探,1997,32(2) :221～228． Wang X W.A method for direct hydrocarbon prediction using joint inversion of gravimetric and seismic data.Oil Geophysical Prospecting(in Chinese)[J],1997,32(2) :221～228．
[17] 杨辉.重力、地震联合反演基岩密度及综合解释[J].石油地球物理勘探,1998,33(4) :496～502． Yang H.Basement density inversion using gravimetric and seismic data and the integrative interpretation[J].Oil Geophysical Prospecting(in Chinese),1998,33(4) :496～ 502．
[18] AndersonB S,Weber M E,Bain J E.A report on the renaissance of gravity in the deep water Gulf of Mexico:A practical view of integration methods[J].68~(th)Ann.Internat Mtg.,Soc.Expi.Geophys.Expanded Abstracts,1998． 515～517．
[19] 贺日政,高锐,李秋生,等.新疆天山(独山子)-西昆仑(泉水沟)地学断面地震与重力联合反演地壳构造特征[J].地球学报,2001,22(6) :553～558． He R Z,Gao R,Li Q S,et al.Corridor gravity fields and crustal density structures in Tianshan(Dushanzi)～ West Kunlun(Quanshuigou)GGT[J].Acta Geoscience Tia Sinica(in Chinese),2001,22(6) :553～558．
[20] 刘天佑,杨宇山,李媛媛,等.重力、地震联合反演方法确定深层油气藏的分布 以锡林凹陷石炭系分布为例[J].天然气工业,2005,25(5) :34～36． Liu T Y,Yang Y S,Li Y Y,et al.Deciding on the distribution of deep reservoir by the joint inversion method of gravity and seismic data-taking the carboniferous in Xilin Sag as an example[J].Natural Gas Industry(in Chinese),2005,25(5) :34～36．
[21] 于鹏.王家林,吴健生,等.重力与地震资料的模拟退火约束联合反演[J].地球物理学报,2007,50(2) :529～538． Yu P,Wang J L,Wu J S,et al.Constrained joint inversion of gravity and seismic data using the simulated annealing algorithm[J].Chinese J.Geophys.(in Chinese),2007,50(3) :529～538．
[22] 于鹏,戴明刚,王家林,等.密度和速度随机分布共网格模型的 重力与地震联合反演[J].地球物理学报,2008,51(3) :845～852． Yu P,Dai M G,Wang J L,et al.Joint inversion of gravity and seismic data based on common gridded model with random density and velocity distributions[J].Chinese J.Geophys.(in Chinese),2008,51(3) :845～852．
[23] ColoMBO D,Stefano M D.Geophysical modeling via simultaneous joint inversion of seismic,gravity and electromagnetic data:Application to prestack depth imaging[J].The Leading Edge.2007,26(3) :326～331．
[24] 杨辉,戴世坤,牟永光.三维重力地震剥层联合反演[J].石油地球物理勘探,2004,39(4) :468～471． Yang H,Dai S K,Mu Y G.Joint layered-stripped inversion of 3 D gravity and seismic data[J].Oil Geophysical Prospecting(in Chinese),2004,39(4) :468 ～471．
[25] Vermeesch P M,Morgan J V,Christeson G Let al,.Three-dimensional joint inversion of traveltime and gravity data across the Chicxulub impact crater[J].Journal of Geophysical Research,2009． 114(doi:10． 1029/ 2008JB005776) :B02105．
[26] Bishop T N,Bube K P,Cutler R T,et al.Tomography determination of the velocity and depth in laterally varying media[J].Geophysics,1985,50:903 ～ 923．
[27] Zelt C A,Smith R B.Seismic traveltime inversion for 2-D crustal velocity structure[J].Geophys.J.Int.1992,108:16～34
[28] 华标龙,刘福田.界面和速度反射联合成像理论与方法[J].地球物理学报,1995,38(6) :750～756． Hua B L,Liu F T.Joint reflection image of velocity and interface-theory and method[J].Chinese J.Geophys.(in Chinese),1995,38(6) :750～756．
[29] 周龙泉,刘福田,陈晓非.三维介质中速度结构和界面的联合成像[J].地球物理学报,2006,49(4) :1062～1067． Zhou L Q,Liu F T,Chen X F.Simultaneous tomography of 3-D velocity structure and interface[J].Chinese J.Geophys.(in Chinese),2006,49(4) :1062～ 1067．
[30] 高尔根,席道瑛,曾昭发.三维结构下全路径迭代射线追踪方法[J].岩石力学与工程学报,2002,21(11) :1610～1614． Gao E G,Xi D Y,Zeng Z F.Whole-Path iterative ray tracing in 3D media[J].Chinese Journal of Rock Mechanics and Engineering.2002,21(11) :1610～1614．
[31] 王家映.地球物理反演理论(第二版)[M].北京.高等教育出版社.2002． Wang J Y.Inverse Theory in Geophysics(second edition)[M].Beijing.Higher education press,2002．
[2] 杨辉,戴世坤,宋海斌,等.综合地球物理联合反演综述[J].地球物理学进展,2002,17(2) :262～271． Yang H,Dai S K,Song H B,et al.Overview of joint inversion of integrated geophysics[J].Progress in Geophysics(in Chinese),2002,17(2) :262～271．
[3] 敬荣中,鲍光淑,陈绍裘.地球物理联合反演研究综述[J].地球物理学进展,2003,18(3) :535～540． Jing R Z,Bao G S,Chen S Q.A review of the researches for geophysical combinative inversion[J].Progress in Geophysics(in Chinese),2003,18(3) :535～540．
[4] 于鹏,王家林,吴健生,等.地球物理联合反演的研究现状和分析[J].勘探地球物理进展,2006,29(2) :87～93． Yu P,Wang J L,Wu J S,et al.Review and discussions on geophysical joint inversion[J].Progress in Exploration Geophysics(in Chinese),2006,29(2) :87～ 93．
[S]Savino J M,Rodi W L,Masso J F.Simultaneous inversion of multiple geophysical data sets for earth structure[J].50~(th)Ann.Internat.Mtg,Soc.Expi.Geophys.Expanded Abstracts,1980,438～439．
[6] Golizdra G Y.Statement of the problem of comprehensive interpretation of gravity fields and seismic observations[J].Izvest.Earth Phys.,1980,16:535～539．
[7] 王一新,廉西京.重力与地震联合解释在石油勘探中的应用[J].石油物探,1983,22(3) :34～41． Wang Y X,Lian X J.Applications of gravity and seismic integrated interpretation in petroleum exploration[J].Geophysical Prospecting for Petroleum(in Chinese),1983,22(3) :34～41．
[8] Lines L,Schultz A K,Treitel S.Cooperative inversion of geophysicaldata[J].57~(th)Ann.Internat.Mtg.,Soc.Expi.Geophys.Expanded Abstracts,1987． 814～816．
[9] Lines L,Schultz A K,Treitel S.Cooperative inversion of geophysical data[J].Geophysics,1988,53(1) :8～20．
[10] 汪宏年,张文生.层状介质中重力、地震联合反演的迭代算法[J].石油地球物理勘探,1993,28(2) :153～165． Wang H N,Zhang W S.Iterative algorithm for joint inversion of gravimetric and seismic data in layered medium[[J].Oil Geophysical Prospecting(in Chinese),1993,28(2) :153～165．
[11] 冯锐,陶裕录.地震重力联合反演中的非块状一致性模型[J].地球物理学报,1993,36(4) :463～475． Feng R,Tao Y L.A non-block consistent model in seismo gravity inversion[J].Chinese J.Geophys.(in Chinese),1993,36(4) :463～475．
[12] 关小平,黄嘉正.谈重震联合反演问题[J].物探与化探,1994,18(6) :429～437． Guan X P,Huang J Z.On the combined gravity seismic inversion problem.Geophysical & Geochemical Exploration(in Chinese),1994,18(6) :429～437．
[13] 范兴才,王家林.松辽盆地北部杏山地区典型的地质地球物理模型及综合解释[A].见:王家林,吴健生主编.中国典型含油气盆地综合地球物理研究[C].上海:同济大学出版社,1995． 112～124． Fan X C,Wang J L.Typical geological and geophysical model and interpretation in Xing shan area of northern Song liao Basin[A].In:Wang J L,Wu J S eds.Integrated Geophysical Researches on Typical Petroleum Basins in China[C](in Chinese).Shanghai:Tongji University Press,19951112～124．
[14] 陈冰,王家林.剥离法与综合模型法的地球物理综合反演[A].见:王家林,吴健生主编.中国典型含油气盆地综合地球物理研究[C].上海:同济大学出版社,1995． 125～135． Chen B,Wang J L.Integrated geophysical inversion by stripping and integrated model methods[A].In:Wang J L,Wu J S eds.Integrated Geophysical Researches on Typical Petroleum Basins in China[C](in Chinese).Shanghai:Tongii University Press,1995． 125～135．
[15] 王西文,孟昭泰,张工会,等.重力、地震联合反演在研究含油气盆地构造中的应用[J].石油地球物理勘探,1991,26(2) :201～209． Wang X W,Meng Z T,Zhang G H,et al.Application of combined gravitational and seismic inversion to study geological structure inside a hydrocarbon bearing basin[J].Oil Geophysical Prospecting(in Chinese),1991,26(2) :201～209．
[16] 王西文.用重力、地震资料联合反演直接预测油气藏的方法[J].石油地球物理勘探,1997,32(2) :221～228． Wang X W.A method for direct hydrocarbon prediction using joint inversion of gravimetric and seismic data.Oil Geophysical Prospecting(in Chinese)[J],1997,32(2) :221～228．
[17] 杨辉.重力、地震联合反演基岩密度及综合解释[J].石油地球物理勘探,1998,33(4) :496～502． Yang H.Basement density inversion using gravimetric and seismic data and the integrative interpretation[J].Oil Geophysical Prospecting(in Chinese),1998,33(4) :496～ 502．
[18] AndersonB S,Weber M E,Bain J E.A report on the renaissance of gravity in the deep water Gulf of Mexico:A practical view of integration methods[J].68~(th)Ann.Internat Mtg.,Soc.Expi.Geophys.Expanded Abstracts,1998． 515～517．
[19] 贺日政,高锐,李秋生,等.新疆天山(独山子)-西昆仑(泉水沟)地学断面地震与重力联合反演地壳构造特征[J].地球学报,2001,22(6) :553～558． He R Z,Gao R,Li Q S,et al.Corridor gravity fields and crustal density structures in Tianshan(Dushanzi)～ West Kunlun(Quanshuigou)GGT[J].Acta Geoscience Tia Sinica(in Chinese),2001,22(6) :553～558．
[20] 刘天佑,杨宇山,李媛媛,等.重力、地震联合反演方法确定深层油气藏的分布 以锡林凹陷石炭系分布为例[J].天然气工业,2005,25(5) :34～36． Liu T Y,Yang Y S,Li Y Y,et al.Deciding on the distribution of deep reservoir by the joint inversion method of gravity and seismic data-taking the carboniferous in Xilin Sag as an example[J].Natural Gas Industry(in Chinese),2005,25(5) :34～36．
[21] 于鹏.王家林,吴健生,等.重力与地震资料的模拟退火约束联合反演[J].地球物理学报,2007,50(2) :529～538． Yu P,Wang J L,Wu J S,et al.Constrained joint inversion of gravity and seismic data using the simulated annealing algorithm[J].Chinese J.Geophys.(in Chinese),2007,50(3) :529～538．
[22] 于鹏,戴明刚,王家林,等.密度和速度随机分布共网格模型的 重力与地震联合反演[J].地球物理学报,2008,51(3) :845～852． Yu P,Dai M G,Wang J L,et al.Joint inversion of gravity and seismic data based on common gridded model with random density and velocity distributions[J].Chinese J.Geophys.(in Chinese),2008,51(3) :845～852．
[23] ColoMBO D,Stefano M D.Geophysical modeling via simultaneous joint inversion of seismic,gravity and electromagnetic data:Application to prestack depth imaging[J].The Leading Edge.2007,26(3) :326～331．
[24] 杨辉,戴世坤,牟永光.三维重力地震剥层联合反演[J].石油地球物理勘探,2004,39(4) :468～471． Yang H,Dai S K,Mu Y G.Joint layered-stripped inversion of 3 D gravity and seismic data[J].Oil Geophysical Prospecting(in Chinese),2004,39(4) :468 ～471．
[25] Vermeesch P M,Morgan J V,Christeson G Let al,.Three-dimensional joint inversion of traveltime and gravity data across the Chicxulub impact crater[J].Journal of Geophysical Research,2009． 114(doi:10． 1029/ 2008JB005776) :B02105．
[26] Bishop T N,Bube K P,Cutler R T,et al.Tomography determination of the velocity and depth in laterally varying media[J].Geophysics,1985,50:903 ～ 923．
[27] Zelt C A,Smith R B.Seismic traveltime inversion for 2-D crustal velocity structure[J].Geophys.J.Int.1992,108:16～34
[28] 华标龙,刘福田.界面和速度反射联合成像理论与方法[J].地球物理学报,1995,38(6) :750～756． Hua B L,Liu F T.Joint reflection image of velocity and interface-theory and method[J].Chinese J.Geophys.(in Chinese),1995,38(6) :750～756．
[29] 周龙泉,刘福田,陈晓非.三维介质中速度结构和界面的联合成像[J].地球物理学报,2006,49(4) :1062～1067． Zhou L Q,Liu F T,Chen X F.Simultaneous tomography of 3-D velocity structure and interface[J].Chinese J.Geophys.(in Chinese),2006,49(4) :1062～ 1067．
[30] 高尔根,席道瑛,曾昭发.三维结构下全路径迭代射线追踪方法[J].岩石力学与工程学报,2002,21(11) :1610～1614． Gao E G,Xi D Y,Zeng Z F.Whole-Path iterative ray tracing in 3D media[J].Chinese Journal of Rock Mechanics and Engineering.2002,21(11) :1610～1614．
[31] 王家映.地球物理反演理论(第二版)[M].北京.高等教育出版社.2002． Wang J Y.Inverse Theory in Geophysics(second edition)[M].Beijing.Higher education press,2002．
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心