二度半长方体组合模型的重力模拟退火反演
|
摘要
针对重力异常反演的现状,采用矩形单元体网格剖分的二度半体组合模型,模拟任意的多层密度界面或形体,实现了重力正演.该建模方法适应性强,适合于实际工作中常见的密度横向变化的复杂模型.使用改进的全局寻优的快速模拟退火算法,对重力异常进行反演,结合这种灵活的密度建模方法,反演过程中只需要反演各矩形单元体的密度参数,即可同时得到地质体的界面或形态以及密度值分布.通过复杂界面和形体的模型试验,证明了方法的效果.在江苏下扬子地质构造复杂区的实际应用说明,在地震等先验信息约束下,该重力反演方法提高了反演精度并减少了多解性,可有效解决古生界泥盆系-志留系目的层分布等地质问题.
In accordance with the present state of gravity anomaly inversion, we use the 2.5 dimensional model construction method to divide the underground density structure into rectangle gridded units to simulate density models with arbitrary interfaces or bodies, then perform gravity forward modeling. This model construction is flexible and adaptable especially for the complicated geological models with horizontal varying densities usually appearing in practical situations. We employ the improved very fast simulated annealing algorithm (VFSA), which has obvious advantage of finding the overall optimal solution in inverse problems to perform gravity anomaly inversion. Only the value of each rectangle unit needs to be inversed, the interfaces or shapes and the density distributions can be simultaneously determined during gravity anomaly inversion.Through the complicated model test, this approach is proved effective and practicable. The application in the lower Yangtze area of Jiangsu shows that this inversion approach could increase the accuracy of inversion under the constrained conditions of known seismic data, the objective strata of Devonian-Silurian and related geological problems could be solved efficiently.
In accordance with the present state of gravity anomaly inversion, we use the 2.5 dimensional model construction method to divide the underground density structure into rectangle gridded units to simulate density models with arbitrary interfaces or bodies, then perform gravity forward modeling. This model construction is flexible and adaptable especially for the complicated geological models with horizontal varying densities usually appearing in practical situations. We employ the improved very fast simulated annealing algorithm (VFSA), which has obvious advantage of finding the overall optimal solution in inverse problems to perform gravity anomaly inversion. Only the value of each rectangle unit needs to be inversed, the interfaces or shapes and the density distributions can be simultaneously determined during gravity anomaly inversion.Through the complicated model test, this approach is proved effective and practicable. The application in the lower Yangtze area of Jiangsu shows that this inversion approach could increase the accuracy of inversion under the constrained conditions of known seismic data, the objective strata of Devonian-Silurian and related geological problems could be solved efficiently.
引文
[1]Chakravarthi V,Singh S B,Babu G A.INVER2DBASE-Aprogram to compute basement depths of density interface above which the density contrast varies with depth.Computers&Geosciences,2001,27:1127~1133
[2]王家林,王一新,林桂康.利用场变换的分离场法反演多层密度界面.石油物探,1986,25(2):69~80Wang J L,Wang YX,Lin GK.Inversionof multi-densityinterfaces by fields-separation of field-transformation.Geophysical Prospecting for Petroleum(in Chinese),1986,25(2):69~80
[3]王西文,孟昭泰,张工会等.重力、地震联合反演在研究含油气盆地构造中的应用.石油地球物理勘探,1991,26(2):201~209Wang X W,Meng Z T,Zhang G H,et al.Application of combinedgravitational and seismic inversionto study geological structure inside a hydrocarbon-bearing basin.Oil Geophysical Prospecting(in Chinese),1991,26(2):201~209
[4]杨辉.重力、地震联合反演基岩密度及综合解释.石油地球物理勘探,1998,33(4):496~502Yang H.Basement density inversion using gravimetric and seismic data and the integrative interpretation.Oil Geophysical Prospecting(in Chinese),1998,33(4):496~502
[5]杨辉,戴世坤,牟永光.三维重力地震剥层联合反演.石油地球物理勘探,2004,39(4):468~471Yang H,Dai S K,Mu Y G.Joint layered-stripped inversion of3-D gravity and seismic data.Oil Geophysical Prospecting(in Chinese),2004,39(4):468~471
[6]王家林,王一新,万明浩等.用重力归一化总梯度法确定密度界面.石油地球物理勘探,1987,22(6):684~690Wang J L,Wang Y X,Wan M H,et al.The determination of densityinterface using the normalized total gravity gradient method.Oil Geophysical Prospecting(in Chinese),1987,22(6):684~690
[7]王一新,王家林,张曙明等.研究多层密度界面的正则化非线性反演方法.石油物探,1987,26(1):78~90Wang Y X,Wang J L,Zhang S M,et al.Normalized nonlinear inversionfor the multi-density-interface.Geophysical Prospectingfor Petroleum(in Chinese),1987,26(1):78~90
[8]杨辉.重力异常多层界面反演的应用.石油地球物理勘探,1992,27(1):140~146Yang H.Application of gravity anomaly inversion for multi-interfaces.Oil Geophysical Prospecting(in Chinese),1992,27(1):140~146
[9]刘云峰,沈晓华.二维密度界面的遗传算法反演.物探化探计算技术,1997,19(2):138~142Liu YF,Shen X H.Nonlinear inversion of gravity anomalies caused by2-Dsurface of geologic structures genetic algorithms.Computing Techniques for Geophysical and Geochemical Exploration(in Chinese),1997,19(2):138~142
[10]陈军,王家林,吴健生等.应用改进的遗传算法反演多层密度界面.地球科学—中国地质大学学报,2000,25(6):651~655Chen J,Wang J L,Wu J S,et al.Application of improved genetic algorithm to inversion of multi-layered density interface.Earth Science-Journal of China University of Geosciences(in Chinese),2000,25(6):651~655
[11]陈军,王家林,吴健生等.用于重力反演的“自适应多段交换”遗传算法.石油物探,2001,40(1):102~106ChenJ,Wang J L,Wu J S,et al.Self-adapting and multi-parts crossing genetic algorithm for gravitational inversion.Geophysical Prospectingfor Petroleum(in Chinese),2001,40(1):102~106
[12]陈东敬,张新兵.带模拟退火的拟BP神经网络在伊朗某地区重力资料反演中的应用.勘探地球物理进展,2005,28(3):215~218Chen D J,Zhang X B.Application of SA-BPNNin inversion of gravity data from Iran.Progress in Exploration Geophysics(in Chinese),2005,28(3):215~218
[13]张新兵,王家林,陈冰等.混合最优化方法在潮汕坳陷重力资料反演中的应用.热带海洋学报,2005,24(2):86~93Zhang X B,Wang J L,Chen B,et al.Application of hybid optimizationin gravityinversion of Chaoshan col.Journal ofTropical Oceanography(in Chinese),2005,24(2):86~93
[14]王家林,王一新,万明浩.石油重磁解释.北京:石油工业出版社,1991Wang J L,Wang Y X,Wan M H.Petroleum Interpretation by Gravity and Magnetic data(in Chinese).Beijing:PetroleumIndustry Press,1991
[15]杨长福.用脊回法反演重力异常的多层密度及其界面.西北地震学报,2004,26(4):293~297Yang C F.Inversion of gravity anomaly for multi-layered densities and their interfaces by ridge regression procedure.Northwestern Seismological Journal(in Chinese),2004,26(4):293~297
[16]杨文采.非线性地球物理反演方法回顾与展望.地球物理学进展,2002,17(2):255~261Yang WC.Non-linear geophysical inversion methods:review and perspective.Progress in Geophysics(in Chinese),2002,17(2):255~261
[17]陈华根,吴健生,王家林.改进的重力模拟退火反演研究.吉林大学学报(地球科学版),2002,32(3):294~298Chen H G,Wu J S,Wang J L.Research in modified simulated annealing gravityinversion.Journal ofJilin University(Earth Science Edition)(in Chinese),2002,32(3):294~298
[18]张霖斌,姚振兴,纪晨等.快速模拟退火算法及应用.石油地球物理勘探,1997,32(5):654~660Zhang LB,Yao Z X,Ji C,et al.Fast simulated annealingalgorithm and its application.Oil Geophysical Prospecting(in Chinese),1997,32(5):654~660
[19]杨辉,王家林,吴健生等.大地电磁与地震资料仿真退火约束联合反演.地球物理学报,2002,45(5):723~734Yang H,Wang J L,Wu J S,et al.Constrained joint inversion of magnetotelluric and seismic data usingsimulated annealing algorithm.Chinese J.Geophys.(in Chinese),2002,45(5):723~734
[20]陈华根.基于GIS技术的大地电磁-重力模拟退火联合反演研究[博士论文].上海:同济大学地球物理系,2001Chen HG.GIS based joint inversion of magnetotelluric and gravity data using simulated annealing algorithm[Ph.D.thesis](in Chinese).Shanghai:Geophysics Department of Tongji University,2001
[21]陈东敬,于鹏.大地电磁与地震联合反演技术在西亚某区的综合应用.石油物探,2005,44(3):251~256Chen DJ,Yu P.Integrated applicationinthe area of west Asia using joint inversion technique of magnetotelluric and seismic data.Geophysical Prospectingfor Petroleum(in Chinese),2005,44(3):251~256
[22]陈永杰,于鹏.大地电磁、地震与重力资料综合反演在徐闻某区的应用.勘探地球物理进展,2005,28(4):271~275Chen Y J,Yu P.Application of integrated joint inversion of magnetotelluric,seismic and gravity data in Xuwen area.Progress in Exploration Geophysics(in Chinese),2005,28(4):271~275
[23]于鹏,王家林,吴健生等.重力与地震资料模拟退火约束联合反演.地球物理学报,2007,50(2):529~538Yu P,WangJ L,WuJ S,et al.Constrainedjoint iversion of gravity and seismic data usingthe simulated annealing algorithm.Chinese J.Geophys.(in Chinese),2007,50(2):529~538
[24]于鹏,王家林,吴健生.大地电磁全频段偏移成像结合反演技术.同济大学学报,2005,33(4):540~544Yu P,WangJ L,WuJ S.Magnetotelluric soubding whole frequency bands migration imaging versus inversion technique.Journal of Tongji University(Natural Science)(in Chinese),2005,33(4):540~544
[2]王家林,王一新,林桂康.利用场变换的分离场法反演多层密度界面.石油物探,1986,25(2):69~80Wang J L,Wang YX,Lin GK.Inversionof multi-densityinterfaces by fields-separation of field-transformation.Geophysical Prospecting for Petroleum(in Chinese),1986,25(2):69~80
[3]王西文,孟昭泰,张工会等.重力、地震联合反演在研究含油气盆地构造中的应用.石油地球物理勘探,1991,26(2):201~209Wang X W,Meng Z T,Zhang G H,et al.Application of combinedgravitational and seismic inversionto study geological structure inside a hydrocarbon-bearing basin.Oil Geophysical Prospecting(in Chinese),1991,26(2):201~209
[4]杨辉.重力、地震联合反演基岩密度及综合解释.石油地球物理勘探,1998,33(4):496~502Yang H.Basement density inversion using gravimetric and seismic data and the integrative interpretation.Oil Geophysical Prospecting(in Chinese),1998,33(4):496~502
[5]杨辉,戴世坤,牟永光.三维重力地震剥层联合反演.石油地球物理勘探,2004,39(4):468~471Yang H,Dai S K,Mu Y G.Joint layered-stripped inversion of3-D gravity and seismic data.Oil Geophysical Prospecting(in Chinese),2004,39(4):468~471
[6]王家林,王一新,万明浩等.用重力归一化总梯度法确定密度界面.石油地球物理勘探,1987,22(6):684~690Wang J L,Wang Y X,Wan M H,et al.The determination of densityinterface using the normalized total gravity gradient method.Oil Geophysical Prospecting(in Chinese),1987,22(6):684~690
[7]王一新,王家林,张曙明等.研究多层密度界面的正则化非线性反演方法.石油物探,1987,26(1):78~90Wang Y X,Wang J L,Zhang S M,et al.Normalized nonlinear inversionfor the multi-density-interface.Geophysical Prospectingfor Petroleum(in Chinese),1987,26(1):78~90
[8]杨辉.重力异常多层界面反演的应用.石油地球物理勘探,1992,27(1):140~146Yang H.Application of gravity anomaly inversion for multi-interfaces.Oil Geophysical Prospecting(in Chinese),1992,27(1):140~146
[9]刘云峰,沈晓华.二维密度界面的遗传算法反演.物探化探计算技术,1997,19(2):138~142Liu YF,Shen X H.Nonlinear inversion of gravity anomalies caused by2-Dsurface of geologic structures genetic algorithms.Computing Techniques for Geophysical and Geochemical Exploration(in Chinese),1997,19(2):138~142
[10]陈军,王家林,吴健生等.应用改进的遗传算法反演多层密度界面.地球科学—中国地质大学学报,2000,25(6):651~655Chen J,Wang J L,Wu J S,et al.Application of improved genetic algorithm to inversion of multi-layered density interface.Earth Science-Journal of China University of Geosciences(in Chinese),2000,25(6):651~655
[11]陈军,王家林,吴健生等.用于重力反演的“自适应多段交换”遗传算法.石油物探,2001,40(1):102~106ChenJ,Wang J L,Wu J S,et al.Self-adapting and multi-parts crossing genetic algorithm for gravitational inversion.Geophysical Prospectingfor Petroleum(in Chinese),2001,40(1):102~106
[12]陈东敬,张新兵.带模拟退火的拟BP神经网络在伊朗某地区重力资料反演中的应用.勘探地球物理进展,2005,28(3):215~218Chen D J,Zhang X B.Application of SA-BPNNin inversion of gravity data from Iran.Progress in Exploration Geophysics(in Chinese),2005,28(3):215~218
[13]张新兵,王家林,陈冰等.混合最优化方法在潮汕坳陷重力资料反演中的应用.热带海洋学报,2005,24(2):86~93Zhang X B,Wang J L,Chen B,et al.Application of hybid optimizationin gravityinversion of Chaoshan col.Journal ofTropical Oceanography(in Chinese),2005,24(2):86~93
[14]王家林,王一新,万明浩.石油重磁解释.北京:石油工业出版社,1991Wang J L,Wang Y X,Wan M H.Petroleum Interpretation by Gravity and Magnetic data(in Chinese).Beijing:PetroleumIndustry Press,1991
[15]杨长福.用脊回法反演重力异常的多层密度及其界面.西北地震学报,2004,26(4):293~297Yang C F.Inversion of gravity anomaly for multi-layered densities and their interfaces by ridge regression procedure.Northwestern Seismological Journal(in Chinese),2004,26(4):293~297
[16]杨文采.非线性地球物理反演方法回顾与展望.地球物理学进展,2002,17(2):255~261Yang WC.Non-linear geophysical inversion methods:review and perspective.Progress in Geophysics(in Chinese),2002,17(2):255~261
[17]陈华根,吴健生,王家林.改进的重力模拟退火反演研究.吉林大学学报(地球科学版),2002,32(3):294~298Chen H G,Wu J S,Wang J L.Research in modified simulated annealing gravityinversion.Journal ofJilin University(Earth Science Edition)(in Chinese),2002,32(3):294~298
[18]张霖斌,姚振兴,纪晨等.快速模拟退火算法及应用.石油地球物理勘探,1997,32(5):654~660Zhang LB,Yao Z X,Ji C,et al.Fast simulated annealingalgorithm and its application.Oil Geophysical Prospecting(in Chinese),1997,32(5):654~660
[19]杨辉,王家林,吴健生等.大地电磁与地震资料仿真退火约束联合反演.地球物理学报,2002,45(5):723~734Yang H,Wang J L,Wu J S,et al.Constrained joint inversion of magnetotelluric and seismic data usingsimulated annealing algorithm.Chinese J.Geophys.(in Chinese),2002,45(5):723~734
[20]陈华根.基于GIS技术的大地电磁-重力模拟退火联合反演研究[博士论文].上海:同济大学地球物理系,2001Chen HG.GIS based joint inversion of magnetotelluric and gravity data using simulated annealing algorithm[Ph.D.thesis](in Chinese).Shanghai:Geophysics Department of Tongji University,2001
[21]陈东敬,于鹏.大地电磁与地震联合反演技术在西亚某区的综合应用.石油物探,2005,44(3):251~256Chen DJ,Yu P.Integrated applicationinthe area of west Asia using joint inversion technique of magnetotelluric and seismic data.Geophysical Prospectingfor Petroleum(in Chinese),2005,44(3):251~256
[22]陈永杰,于鹏.大地电磁、地震与重力资料综合反演在徐闻某区的应用.勘探地球物理进展,2005,28(4):271~275Chen Y J,Yu P.Application of integrated joint inversion of magnetotelluric,seismic and gravity data in Xuwen area.Progress in Exploration Geophysics(in Chinese),2005,28(4):271~275
[23]于鹏,王家林,吴健生等.重力与地震资料模拟退火约束联合反演.地球物理学报,2007,50(2):529~538Yu P,WangJ L,WuJ S,et al.Constrainedjoint iversion of gravity and seismic data usingthe simulated annealing algorithm.Chinese J.Geophys.(in Chinese),2007,50(2):529~538
[24]于鹏,王家林,吴健生.大地电磁全频段偏移成像结合反演技术.同济大学学报,2005,33(4):540~544Yu P,WangJ L,WuJ S.Magnetotelluric soubding whole frequency bands migration imaging versus inversion technique.Journal of Tongji University(Natural Science)(in Chinese),2005,33(4):540~544
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心