快速模拟退火法重力及其梯度异常与地震数据的联合反演
|
摘要
本文提出采用快速模拟退火方法进行重力及其梯度数据与地震数据的联合反演,可以改善以往该类方法对界面的突变区域无法识别的缺陷。首先推导出密度界面起伏与重力导数之间的对应关系,为更好描述实际地层的水平密度变化,将地层之间的密度差设为随位置变化的函数;然后通过理论模型试验证明该方法能很好地进行界面反演,且对于界面突变区也有很好的效果;最后应用该方法对M区实测数据进行处理。结果表明,利用改进的快速模拟退火法进行重力及其梯度数据与地震数据的联合反演,能很好地揭示物性界面的分布特征,且计算稳定、精度高。
Conventional inversion methods cannot recognize sudden change surfaces.To solve this problem,we present in this paper a joint inversion method of gravity,gravity gradient and seismic data using the very fast simulated annealing(VFSA).First,we derive the corresponding relation between density interface and gravity derivative.The density contrasts are set to be functions variable with locations to delineate the horizontal variation of density.Then,we use the VFSA method on model data,and accomplish successfully the inversion of sudden change surfaces.At last,this VFSA method is applied to real data acquired in M area.Application results show that our VFSA method can accomplish successfully the inversion,and effectively display the surface distribution with the advantages of computation stable and high precision.
Conventional inversion methods cannot recognize sudden change surfaces.To solve this problem,we present in this paper a joint inversion method of gravity,gravity gradient and seismic data using the very fast simulated annealing(VFSA).First,we derive the corresponding relation between density interface and gravity derivative.The density contrasts are set to be functions variable with locations to delineate the horizontal variation of density.Then,we use the VFSA method on model data,and accomplish successfully the inversion of sudden change surfaces.At last,this VFSA method is applied to real data acquired in M area.Application results show that our VFSA method can accomplish successfully the inversion,and effectively display the surface distribution with the advantages of computation stable and high precision.
引文
[1]蒋龙聪,刘江平.模拟退火算法及其改进.工程地球物理学报,2007,4(2):135～140Jiang Longcong,Liu Jiangping.Revised simulatedannealing algorithm.Chinese Journal of EngineeringGeophysics,2007,4(2):135～140
[2]敬荣中,鲍光淑,陈绍裘.地球物理联合反演研究综述.地球物理学进展,2003,18(3):535～540Jing Rongzhong,Bao Guangshu,Chen Shaoqiu.Areview of the researches for geophysical combinativeinversion.Progress in Geophysics,2003,18(3):535～540
[3]刘天佑.位场勘探数据处理新方法.北京:科学出版社,2007
[4]于鹏,王家林,吴健生.地球物理联合反演的研究现状和分析.勘探地球物理进展,2006,29(2):87～93Yu Peng,Wang Jialin,Wu Jiansheng.Review and dis-cussions on geophysical joint inversion.PEG,2006,29(2):87～93
[5]杨辉.重力、地震联合反演基岩密度及综合解释.石油地球物理勘探,1998,33(4):496～502Yang Hui.Basement density inversion using gravi-metric and seismic data and the integrative interpreta-tion.OGP,1998,33(4):496～502
[6]Vozoff K,Jupp D L B.Effective search for a buried-layer:An approach to experimental design in geo-physics.Expl Geophys,1977,8(1):6～15
[7]Savino J M,Rodi W L,Masso J F.Simultaneous in-version of multiple geophysical data sets for earthstructure.The 50th Ann Internat Mtg of SEG,1980,438～439
[8]Menichetti V,Guillen A.Simultaneous interactivemagnetic and gravity inversion.Geophysical Prospec-ting,1983,31(6):929～944
[9]Gomez-Trevino E,Edwards R N.Electromagneticsoundings in the sedimentary basin of southern On-tario-a case history.Geophysics,1983,48(3):311～330
[10]Mottl J,Mottlova L.The simultaneous solution ofthe inverse problem of gravimetry and magnetics bymeans of nonliner programming.Geophysical Journalof the Royal Astronomical Society,1984,76(3):563～579
[11]Chavez R E et al.Joint interpretation of gravity andmagnetic data over axial symmetric bodies with appli-cation to the Darnley bay anomaly,Nwt Canada.Geo-physical Prospecting,1987,35(4):374～392
[12]Lines L,Schultz A K,Treitel S.Cooperative inver-sion of geophysical data.SEG Technical ProgramExpanded Abstracts,1987,6:814～816
[13]Lines L,Schultz A K,Treitel S.Cooperative inver-sion of geophysical data.Geophysics,1988,53(1):8～20
[14]Sasaki Y.Two-dimensional joint inversion of magne-totel-luric and dipoledipole resistivity data.Geophys-ics,1989,54(2):254～262
[15]Dobroka M,Gyulai A,Ormos T et al.Joint inver-sion of seismic and geoelectric data recorded in an un-derground coal mine.Geophysics Prospecting,1991,39(5):643～665
[16]Rasmussen H J,Hohmann G W.Cooperative inver-sion of transient electromagnetic and gravity data fordepth determinations in a basin.SEG Technical Pro-gram Expanded Abstracts,1991,10:614～617
[17]Zeyen H,Pous J.3-D joint inversion of magnetic andgravimetric data with a priori information.Geophysi-cal Journal International,1993,112(2):244～256
[18]Zhao S.Joint inversion of observed gravity and GPSbaseline changes for the detection of the active faultsegment at the Red River fault zone.GeophysicalJournalInternational,1995,122(1):70～88
[19]Hering A,Misiek R,Gyulai et al.A joint inversionalgorithm to process geoelectric and surfacewave seis-mic data,partⅠ:Basic ideas.Geophysics Prospec-ting,1995,43(2):135～156
[20]Vasco D W,Peterson J E,Majer E L.A simultane-ous inversion of seismic traveltimes and amplitudesfor velocity and attenuation.Geophysics,1996,61(6):1738～1757
[21]Tezkan B et al.A joint application of radio magneto-tellurics and transient electromagnetics to the investi-gation of a waste deposit in Cologne.Journal of Ap-plied Geophysics,1996,34(3):199～212
[22]Grechka V,Theophanis S,Tsvankin I.Joint inver-sion of P-and PS-waves in orthorhombic media:The-ory and a Physical-modeling study.SEG TechnicalProgram Expanded Abstracts,1997,16:1250～1253
[23]Roth M and Holliger K.Joint inversion of Rayleighand guided waves in high-resolution seismic data u-sing agenetic algorithm.SEG Technical ProgramExpanded Abstracts,1998,17:1570～1573
[24]Chunduru R,Zhang Z,Busch R et al.Joint inversionof MWD and wireline measurements.SEG TechnicalProgram Expanded Abstracts,1999,18:132～135
[25]Nielsen L,Jacobsen B H.Integrated gravity andwide-angle seismic inversion for two-dimensionalcrustal modeling.Geophysical Journal Internation-al,2000,140(1):222～232
[26]Bosch M,McGaughey J.Joint inversion of gravityand magnetic data under lithologic constraints.TheLeading Edge,2001,20(8):877～881
[27]Kuhn E,Lavielle M and Zeyen H.Joint inversion ofteleseismic delay times and gravity anomaly data,Anew approach.Geophysical Journal International,2002,(3):1～8
[28]杨辉,戴世坤,牟永光.三维重力地震剥层联合反演.石油地球物理勘探,2004,39(4):468～471Yang Hui,Dai Shikun and Mu Yongguang.Joint lay-er-stripped inversion of 3-D gravity and seismic data.OGP,2004,39(4):468～471
[29]于鹏,王家林,吴健生等.重力与地震资料的模拟退火约束联合反演.地球物理学报,2007,50(2):529～538Yu Peng,Wang Jialin,Wu Jiansheng et al.Con-strained joint inversion of gravity and seismic data u-sing the simulated annealing algorithm.Chinese JGeophys,2007,50(2):529～538
[30]于鹏,戴明刚,王家林等.密度和速度随机分布共网格模型的重力与地震联合反演.地球物理学报,2008,51(3):845～852Yu Peng,Dai Minggang,Wang Jialin et al.Joint in-version of gravity and seismic data based on commongridded model with random density and velocity dis-tributions.Chinese J Geophys,2008,51(3):845～852
[31]Parker R L.The rapid calculation of Potential anoma-lies.Geophys Jour Roy Astro Soc,1972,31(2):447～455
[32]Oldenburg D W.The inversion and interpretation ofgravity anomalies.Geophysics,1974,9(4):526～536
[33]张霖斌,姚振兴,纪晨等.快速模拟退火算法及应用.石油地球物理勘探,1997,32(5):654～660Zhang Linbin,Yao Zhenxing,Ji Chen et al.Fastsimulated annealing algorithm and its application.OGP,1997,32(5):654～660
[34]任义庆,徐仲达,马在田.应用模拟退火法反演横波速度.石油地球物理勘探,1996,31(5):677～684Ren Yiqing,Xu Zhongda and Ma Zaitian.S-wave ve-locity inversion uing simulated annealing method.OGP,1996,31(5):677～684
[35]Szu H and Hartley R.Nonconvex optimization byfast simulated annealing.Proc IEEE,1987,75(11):1538～1540
[36]lngber L.Very fast simulated re-annealing.MathlComput Modeling,1989,12(8):967～993
[2]敬荣中,鲍光淑,陈绍裘.地球物理联合反演研究综述.地球物理学进展,2003,18(3):535～540Jing Rongzhong,Bao Guangshu,Chen Shaoqiu.Areview of the researches for geophysical combinativeinversion.Progress in Geophysics,2003,18(3):535～540
[3]刘天佑.位场勘探数据处理新方法.北京:科学出版社,2007
[4]于鹏,王家林,吴健生.地球物理联合反演的研究现状和分析.勘探地球物理进展,2006,29(2):87～93Yu Peng,Wang Jialin,Wu Jiansheng.Review and dis-cussions on geophysical joint inversion.PEG,2006,29(2):87～93
[5]杨辉.重力、地震联合反演基岩密度及综合解释.石油地球物理勘探,1998,33(4):496～502Yang Hui.Basement density inversion using gravi-metric and seismic data and the integrative interpreta-tion.OGP,1998,33(4):496～502
[6]Vozoff K,Jupp D L B.Effective search for a buried-layer:An approach to experimental design in geo-physics.Expl Geophys,1977,8(1):6～15
[7]Savino J M,Rodi W L,Masso J F.Simultaneous in-version of multiple geophysical data sets for earthstructure.The 50th Ann Internat Mtg of SEG,1980,438～439
[8]Menichetti V,Guillen A.Simultaneous interactivemagnetic and gravity inversion.Geophysical Prospec-ting,1983,31(6):929～944
[9]Gomez-Trevino E,Edwards R N.Electromagneticsoundings in the sedimentary basin of southern On-tario-a case history.Geophysics,1983,48(3):311～330
[10]Mottl J,Mottlova L.The simultaneous solution ofthe inverse problem of gravimetry and magnetics bymeans of nonliner programming.Geophysical Journalof the Royal Astronomical Society,1984,76(3):563～579
[11]Chavez R E et al.Joint interpretation of gravity andmagnetic data over axial symmetric bodies with appli-cation to the Darnley bay anomaly,Nwt Canada.Geo-physical Prospecting,1987,35(4):374～392
[12]Lines L,Schultz A K,Treitel S.Cooperative inver-sion of geophysical data.SEG Technical ProgramExpanded Abstracts,1987,6:814～816
[13]Lines L,Schultz A K,Treitel S.Cooperative inver-sion of geophysical data.Geophysics,1988,53(1):8～20
[14]Sasaki Y.Two-dimensional joint inversion of magne-totel-luric and dipoledipole resistivity data.Geophys-ics,1989,54(2):254～262
[15]Dobroka M,Gyulai A,Ormos T et al.Joint inver-sion of seismic and geoelectric data recorded in an un-derground coal mine.Geophysics Prospecting,1991,39(5):643～665
[16]Rasmussen H J,Hohmann G W.Cooperative inver-sion of transient electromagnetic and gravity data fordepth determinations in a basin.SEG Technical Pro-gram Expanded Abstracts,1991,10:614～617
[17]Zeyen H,Pous J.3-D joint inversion of magnetic andgravimetric data with a priori information.Geophysi-cal Journal International,1993,112(2):244～256
[18]Zhao S.Joint inversion of observed gravity and GPSbaseline changes for the detection of the active faultsegment at the Red River fault zone.GeophysicalJournalInternational,1995,122(1):70～88
[19]Hering A,Misiek R,Gyulai et al.A joint inversionalgorithm to process geoelectric and surfacewave seis-mic data,partⅠ:Basic ideas.Geophysics Prospec-ting,1995,43(2):135～156
[20]Vasco D W,Peterson J E,Majer E L.A simultane-ous inversion of seismic traveltimes and amplitudesfor velocity and attenuation.Geophysics,1996,61(6):1738～1757
[21]Tezkan B et al.A joint application of radio magneto-tellurics and transient electromagnetics to the investi-gation of a waste deposit in Cologne.Journal of Ap-plied Geophysics,1996,34(3):199～212
[22]Grechka V,Theophanis S,Tsvankin I.Joint inver-sion of P-and PS-waves in orthorhombic media:The-ory and a Physical-modeling study.SEG TechnicalProgram Expanded Abstracts,1997,16:1250～1253
[23]Roth M and Holliger K.Joint inversion of Rayleighand guided waves in high-resolution seismic data u-sing agenetic algorithm.SEG Technical ProgramExpanded Abstracts,1998,17:1570～1573
[24]Chunduru R,Zhang Z,Busch R et al.Joint inversionof MWD and wireline measurements.SEG TechnicalProgram Expanded Abstracts,1999,18:132～135
[25]Nielsen L,Jacobsen B H.Integrated gravity andwide-angle seismic inversion for two-dimensionalcrustal modeling.Geophysical Journal Internation-al,2000,140(1):222～232
[26]Bosch M,McGaughey J.Joint inversion of gravityand magnetic data under lithologic constraints.TheLeading Edge,2001,20(8):877～881
[27]Kuhn E,Lavielle M and Zeyen H.Joint inversion ofteleseismic delay times and gravity anomaly data,Anew approach.Geophysical Journal International,2002,(3):1～8
[28]杨辉,戴世坤,牟永光.三维重力地震剥层联合反演.石油地球物理勘探,2004,39(4):468～471Yang Hui,Dai Shikun and Mu Yongguang.Joint lay-er-stripped inversion of 3-D gravity and seismic data.OGP,2004,39(4):468～471
[29]于鹏,王家林,吴健生等.重力与地震资料的模拟退火约束联合反演.地球物理学报,2007,50(2):529～538Yu Peng,Wang Jialin,Wu Jiansheng et al.Con-strained joint inversion of gravity and seismic data u-sing the simulated annealing algorithm.Chinese JGeophys,2007,50(2):529～538
[30]于鹏,戴明刚,王家林等.密度和速度随机分布共网格模型的重力与地震联合反演.地球物理学报,2008,51(3):845～852Yu Peng,Dai Minggang,Wang Jialin et al.Joint in-version of gravity and seismic data based on commongridded model with random density and velocity dis-tributions.Chinese J Geophys,2008,51(3):845～852
[31]Parker R L.The rapid calculation of Potential anoma-lies.Geophys Jour Roy Astro Soc,1972,31(2):447～455
[32]Oldenburg D W.The inversion and interpretation ofgravity anomalies.Geophysics,1974,9(4):526～536
[33]张霖斌,姚振兴,纪晨等.快速模拟退火算法及应用.石油地球物理勘探,1997,32(5):654～660Zhang Linbin,Yao Zhenxing,Ji Chen et al.Fastsimulated annealing algorithm and its application.OGP,1997,32(5):654～660
[34]任义庆,徐仲达,马在田.应用模拟退火法反演横波速度.石油地球物理勘探,1996,31(5):677～684Ren Yiqing,Xu Zhongda and Ma Zaitian.S-wave ve-locity inversion uing simulated annealing method.OGP,1996,31(5):677～684
[35]Szu H and Hartley R.Nonconvex optimization byfast simulated annealing.Proc IEEE,1987,75(11):1538～1540
[36]lngber L.Very fast simulated re-annealing.MathlComput Modeling,1989,12(8):967～993
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心