航空重力梯度的模拟退火参数反演Inversion of airborne gravity gradient by using simulated annealing algorithm
刘金钊;王岩;高春春;李永昆;王勇;
摘要(Abstract):
针对目前地球物理问题反演时,通常利用线性或拟线性反演方法,从而导致反演过程不稳定或反演结果陷于局部最优的问题,该文采用了模拟退火算法的非线性全局反演技术,从航空重力梯度全张量数据中反演了地下密度异常体的几何参数信息。同时,利用无噪声和有噪声的两组航空重力梯度数据进行对比实验。结果表明,在重力梯度数据信噪比小于21.584 3的情况下,通过引入噪声加权的价值函数,可以提高地下密度异常体几何参数反演的准确性。
关键词(KeyWords): 重力梯度;模拟退火算法;参数反演;加权
基金项目(Foundation): 中国地震局星火计划项目(XH17058Y);中国地震局第一监测中心主任创新基金项目(FMC2016004);; 天津市自然科学基金项目(17JCYBJC21600)
作者(Authors): 刘金钊;王岩;高春春;李永昆;王勇;
DOI: 10.16251/j.cnki.1009-2307.2018.08.002
参考文献(References):
- [1]王浩然,陈超,杜劲松.重力梯度张量数据的三维反演方法与应用[J].石油地球物理勘探,2013,48(3):474-481.(WANG Haoran,CHEN Chao,DU Jinsong.3-D inversion of gravity gradient tensor data and its applications[J].Oil Geophysical Prospecting,2013,48(3):474-481.)
- [2]陈少华,朱自强,鲁光银,等.重力梯度张量的预条件共轭梯度法反演[J].中南大学学报(自然科学版),2013,44(2):619-625.(CHEN Shaohua,ZHU Ziqiang,LU Guangyin,et al.Inversion of gravity gradient tensor based on preconditioned conjugate gradient[J].Journal of Central South University(Science and Technology),2013,44(2):619-625.)
- [3]ZHDANOV M S,ELLIS R,MUKHERJEE S.Three-dimensional regularized focusing inversion of gravity gradient tensor component data[J].Geophysics,2004,69(4):925-937.
- [4]党亚民,陈俊勇,晁定波.模拟退火算法及其在大地测量反演中的应用[J].测绘科技动态,1999(1):8-14.(DANG Yamin,CHEN Junyong,CHAO Dingbo.Simulated annealing algorithm and its application in geodetic inversion[J].Surveying and Mapping Technology,1999(1):8-14.)
- [5]姚姚.地球物理反演基本理论与应用方法[M].武汉:中国地质大学出版社,2002.(YAO Yao.Basic theory and application method of geophysical inversion[M].Wuhan:China University of Geosciences Press,2002.)
- [6]NAGIHARA S,HALL S N.Three-dimensional gravity inversion using simulated annealing:constraints on the diapiric roots of allochthonous salt structures[J].Geophysics,2001,66(5):1438-1449.
- [7]陈华根,吴健生,王家林.改进的重力模拟退火反演研究[J].吉林大学学报(地球科学版),2002,32(3):294-298,303.(CHEN Huagen,WU Jiansheng,WANG Jialin.Research in modified simulated-annealing gravity inversion[J].Journal of Jilin University(Earth Science Edition),2002,32(3):294-298,303.)
- [8]于鹏,王家林,吴健生.二度半长方体组合模型的重力模拟退火反演[J].地球物理学报,2007,50(3):882-889.(YU Peng,WANG Jialin,WU Jiansheng.An inversion of gravity anomalies by using a 2.5dimensional rectangle gridded model and the simulated annealing algorithm[J].Chinese Journal of Geophysics,2007,50(3):882-889.)
- [9]于鹏,王家林,吴健生,等.重力与地震资料的模拟退火约束联合反演[J].地球物理学报,2007,50(2):529-538.(YU Peng,WANG Jialin,WU Jiansheng,et al.Constrained joint inversion of gravity and seismic data using the simulated annealing algorithm[J].Chinese Journal of Geophysics,2007,50(2):529-538.)
- [10]郭文斌,朱自强,鲁光银.重力梯度张量的拟BP神经网络反演[J].中南大学学报(自然科学版),2011,42(12):3797-3803.(GUO Wenbin,ZHU Ziqing,LU Guangyin.Quasi-BP neural network inversion of gravity gradient tensor[J].Journal of Central South University(Science and Technology),2011,42(12):3797-3803.)
- [11]GARCA-ABDESLEM J.Nonlinear 3-D inversion of gravity data over a sulfide ore body[J].Geofisica Internacional,2000,39(2):179-188.
- [12]BEIKI M,PEDERSEN L B,NAZI H.Interpretation of aeromagnetic data using eigenvector analysis of pseudogravity gradient tensor[J].Geophysics,2011,76(3):L1-L10.
- [13]OLOWOFELA J A,BADMUS B S,GANIYU S A,et al.Source location and depth estimation from digitized aeromagnetic data acquired from a basement complex formation[J].Earth Science India,2011,4(3):136-142.
- [14]高秀鹤,黄大年.基于共轭梯度算法的重力梯度数据三维聚焦反演研究[J].地球物理学报,2017,60(4):1571-1583.(GAO Xiuhe,HUANG Danian.Research on 3D focusing inversion of gravity gradient tensor data based on a conjugate gradient algorithm[J].Chinese J Geophys,2017,60(4):1571-1583.)
- [15]李红蕾,方剑,王新胜,等.重力及重力梯度联合反演青藏高原及邻区岩石圈三维密度结构[J].地球物理学报,2017,60(6):2469-2479.(LI Honglei,FANG Jian,WANG Xinsheng,et al.Lithoshperic 3Ddensity structure beneath the Tibetan plateau and adjacent areas derived from joint inversion of satellite gravity and gravity-gradient data[J].Chinese J Geophys,2017,60(6):2469-2479.)
- [16]UZUN S.Estimating parameters of subsurface structures from airborne gravity gradiometry data using a Monte-Carlo optimization method[D].Columbus,Ohio,USA:the Ohio State University,2013.
- [17]CORANA A,MARCHESI M,MARTINI C,et al.Minimizing multimodal functions of continuous variables with the“simulated annealing”algorithm[J].ACM Transactions on Mathematical Software,1987,13(3):262-280.
- [18]JEKELI C.Statistical analysis of moving-base gravimetry and gravity gradiometry[R/OL].(2003-09)[2017-02-08].https:∥kb.osu.edu/dspace/bitstream/handle/1811/78651/SES_GeodeticScience_Report_466.pdf?sequence=1.