曾昭发,,
孙思源,
于平
吉林大学地球探测科学与技术学院, 长春 130000
基金项目: 国家重点研发计划(2016YFC0601104)和国家自然科学基金项目(41574097)联合资助
详细信息
作者简介: 高秀鹤, 女, 1991年生, 博士, 主要从事重磁及张量梯度数据反演研究.E-mail:842091616@qq.com
通讯作者: 曾昭发, 男, 1966年生, 博士生导师, 主要从事地球物理方法理论与应用研究.E-mail:zengzf@jlu.edu.cn
中图分类号: P631收稿日期:2018-01-15
修回日期:2018-03-21
上线日期:2019-03-05
Joint inversion of gravity and gravity gradient data based on Cokriging method with the threshold constraint
GAO XiuHe,ZENG ZhaoFa,,
SUN SiYuan,
YU Ping
College of GeoExploration Science and Technology, Jilin University, Changchun 130000, China
More Information
Corresponding author: ZENG ZhaoFa,E-mail:zengzf@jlu.edu.cn
MSC: P631--> Received Date: 15 January 2018
Revised Date: 21 March 2018
Available Online: 05 March 2019
摘要
摘要:常规协克里金方法反演重力或重力梯度数据具有抗噪性好、加入先验信息容易等优点,其反演的地下密度分布能够识别异常体中心位置,还原异常体基本形态,但反演图像光滑,分辨率低,这是由于常规方法估计的密度协方差矩阵全局发散、平稳.为了通过协克里金方法获得聚焦的密度分布需要改善密度协方差矩阵的性质.首先,本文推导了理论密度协方差公式,其性质表明,当理论模型聚焦分布时,其密度协方差矩阵是非平稳且聚焦分布的.为了打破常规协方差矩阵全局平稳、发散的特征,本文设置密度阈值处理协方差矩阵,通过不断更新协方差矩阵来迭代实现协克里金反演,最终得到相对聚焦的反演结果.用本文方法处理重力与重力梯度数据恢复两种密度模型,均得到了与正演模型匹配的反演结果;再将方法运用于文顿盐丘的实际测量重力与重力梯度数据,反演结果与已知的地质情况匹配较好.
关键词: 协克里金/
联合反演/
重力/
重力梯度/
阈值
Abstract:The conventional Cokriging inversion of gravity and gravity gradient data has the advantages of good anti-noise performance and being easy to add prior information. And this method can identify the position of the anomaly body and to restore its basic shape. But the inversion image of this approach is smooth and its resolution is low which are due to the density covariance matrix estimated by the conventional method are of global divergence and stationary. So the performance of the density covariance matrix has to be improved for getting focusing density distribution by the Cokriging inversion. Firstly, this paper deduces the theoretical formula of density covariance matrix, which proves that the density covariance matrix is non-stationary and focused when the theoretical model is compact. In order to break the global stationary and divergence of the conventional covariance matrix, we set a density threshold to deal with the covariance matrix and update the covariance matrix to iteratively perform Cokriging inversion until the focused results are obtained. Then we set two kinds of forward models to testify our methods, showing the recovered models match with the forward models. Finally we use our method to process the gravity and gravity gradient data from the Vinton salt dome, and the inversion results are in good agreement with the known geological structures.
Key words:Cokriging/
Joint inversion/
Gravity/
Gravity gradient/
Threshold
PDF全文下载地址:
http://www.geophy.cn/data/article/export-pdf?id=dqwlxb_14908