一种基于重力正演理论的海底地形反演迭代算法

1

作者 公维梁 屠泽杰 +3 位作者 孙月文 邢赛 赵福玺 阳凡林 《山东科技大学学报（自然科学版）》 CAS 北大核心 2024年第3期64-74,共11页

重力数据反演海底地形的方法是获取全球海底地形数据的重要途径。针对传统地形反演算法需要间接引入经验参数的问题,本研究以不采用经验公式为原则,从空间域方法入手,结合重力正演公式,建立重力异常与地形之间的解析观测方程,采用最小... 重力数据反演海底地形的方法是获取全球海底地形数据的重要途径。针对传统地形反演算法需要间接引入经验参数的问题,本研究以不采用经验公式为原则,从空间域方法入手,结合重力正演公式,建立重力异常与地形之间的解析观测方程,采用最小二乘方法进行求解。为了得到观测方程最优解,建立牛顿迭代关系,引入正则化参数,增强方程组的收敛性。顾及边缘效应对正演的影响,采用双线性插值算法对原始重力异常数据进行加密处理和无约束网格扩充以削弱误差。在太平洋海域(155°E~156°E,16°N~17°N)构建了1′×1′分辨率的海底地形模型,并通过船测水深数据验证。相比于传统空间域算法——重力地质法,本研究方法的均方根误差降低了12.7%,验证了其可行性与准确性。 展开更多

关键词 重力数据 重力正演 地形反演 正则化参数 边缘效应

下载PDF 职称材料

渤海重力三维正演研究的初步结果 被引量：5

2

作者 刘光夏 任文菊 +2 位作者 吴岫云 赵文俊 李志雄 《地震学报》 CSCD 北大核心 1997年第1期102-106,共5页

我们使用重力三维正演方法,把研究区内各时代地质体对布格异常的影响一一消除（正演剥皮）,最后得到深部重力异常.本区有6个深部重力高区,即渤中、渤海湾、辽东湾、渤南-莱州湾以及垦利和义和庄;一个深部重力低区。

关键词 渤海 重力正演 重力异常 地质构造

下载PDF 职称材料

起伏地表下的复杂三维地质模型建立与重力异常计算 被引量：5

3

作者 贾卓 刘四新 +3 位作者 赵雪然 鹿琪 李宏卿 王元新 《吉林大学学报（地球科学版）》 EI CAS CSCD 北大核心 2021年第1期277-285,共9页

由于地质体和矿体的形态非常复杂,使用长方体网格离散建立正演模型时可能和真实情况有很大差别,因此计算结果可靠性差。本文提出一种基于约束Delaunay网格剖分的方法对地质体进行离散并进行重力建模,在模型边界等复杂区域使用网格自适... 由于地质体和矿体的形态非常复杂,使用长方体网格离散建立正演模型时可能和真实情况有很大差别,因此计算结果可靠性差。本文提出一种基于约束Delaunay网格剖分的方法对地质体进行离散并进行重力建模,在模型边界等复杂区域使用网格自适应加密技术,将三维地质体离散为有限个四面体;并详细推导出针对四面体网格的重力正演公式,实现了基于约束Delaunay网格剖分技术的三维重力数值模拟;最后,针对一个合成数据模型,将计算解与解析解对比。结果表明,细化网格的模拟结果比粗糙网格更好,满足数值模拟的精度要求。将该方法应用到金川矿区实际地质体建模中,根据局部需要,建立各处网格密度不均匀的三维模型,并计算该模型的地表重力场,而后对比模拟数据与实测数据,结果表明Delaunay网格建模方法具有很强的适用性,能够模拟复杂的地质体重力异常。 展开更多

关键词 Delaunay网格 四面体 三维地质体模型 重力正演 金川矿区

下载PDF 职称材料

利用重力方法圈定巨野煤田成矿远景区 被引量：4

4

作者 李土雄 刘磊 +1 位作者 禚传源 姚敏 《山东国土资源》 2012年第7期28-31,共4页

重力方法是研究区域构造的重要方法,通过重力正演能准确地计算出具有一定密度差的地层分界面深度。通过收集巨野煤田早期的岩石密度资料以及钻孔资料,利用2.5D剖面重力正演计算方法,圈定了巨野煤田-1 000m以浅成矿远景区。该方法对其他... 重力方法是研究区域构造的重要方法,通过重力正演能准确地计算出具有一定密度差的地层分界面深度。通过收集巨野煤田早期的岩石密度资料以及钻孔资料,利用2.5D剖面重力正演计算方法,圈定了巨野煤田-1 000m以浅成矿远景区。该方法对其他地区盆地内沉积矿产找矿具有借鉴意义。 展开更多

关键词 重力方法 重力正演 巨野煤田

下载PDF 职称材料

基于起伏地层的曲面重力场快速高精度正演 被引量：2

5

作者 郑翾宇 柳建新 +3 位作者 陈波 赵广东 陈龙伟 郭荣文 《物探化探计算技术》 CAS 2019年第2期214-222,共9页

Parker公式在起伏地层重力场正演方面以其简洁、快速等优点被广泛应用于地形改正以及界面反演等,但传统的Parker公式正演方法存在一定的缺陷,在起伏地层模型重力场的正演计算中,由于正演式中e的指数项数值稳定性较差,导致正演结果精度... Parker公式在起伏地层重力场正演方面以其简洁、快速等优点被广泛应用于地形改正以及界面反演等,但传统的Parker公式正演方法存在一定的缺陷,在起伏地层模型重力场的正演计算中,由于正演式中e的指数项数值稳定性较差,导致正演结果精度不高。在前人研究的基础上,对Parker公式进行了改进,通过增加上、下界面平均值的方法,提高正演的数值稳定性,改进后的正演精度得到大幅提升。在此基础上,提出了一种曲面观测算法,将观测面由传统的水平面推广到任意起伏面,实现了上下界面起伏地层在起伏观测面上的快速高精度重力正演计算和带地形的密度界面起伏模型的快速高精度重力正演模拟。为了保障数值精度,采用了Gauss-FFT算法,模型实验均取得很好效果。 展开更多

关键词 重力正演 Parker公式 曲面观测 起伏地形

下载PDF 职称材料

重力—地震联合建模正演技术在SH地区潜山勘探中的应用 被引量：2

6

作者 李婧 崔永谦 +3 位作者 宋景明 赵荔 周赏 梁婉如 《石油地球物理勘探》 EI CSCD 北大核心 2016年第S1期152-160,11,共10页

渤海湾盆地油气藏勘探日趋成熟,寻找深层潜山油气目标成为今后勘探的方向,具有重要意义。本文阐述了在地震三维区应用重力勘探方法落实深部潜山目标的一种新技术,也是重力地震综合勘探的一种新应用:以钻井为约束,以地震为依托,联合构建... 渤海湾盆地油气藏勘探日趋成熟,寻找深层潜山油气目标成为今后勘探的方向,具有重要意义。本文阐述了在地震三维区应用重力勘探方法落实深部潜山目标的一种新技术,也是重力地震综合勘探的一种新应用:以钻井为约束,以地震为依托,联合构建地质模型,并以测井数据为纽带研究速度、密度、深度之间的相互转换关系。通过重力正演,计算各层位的异常响应,进而分析各地层与布格异常的关系,从而获得潜山构造的重力异常响应。通过研究模型规模与重力异常的对应关系,可以准确了解不同规模及埋深的潜山所引起的重力异常的幅值变化,从而推断地震解释正确与否,为地质解释提供合理可靠的辅助依据。 展开更多

关键词 重力正演 建模 潜山 综合勘探

下载PDF 职称材料

三维Tesseroid网格模型重力异常正演方法及并行算法 被引量：2

7

作者 王博 郭良辉 +1 位作者 崔亚彤 王祥 《物探与化探》 CAS 北大核心 2021年第6期1597-1605,共9页

三维网格模型的正演计算是重力资料反演的基础。高精度、高效率的正演有利于提高反演解释质量。针对大尺度、地表观测面研究区的高精度、高效率重力正演问题,本文给出球坐标系三维Tesseroid网格模型重力异常正演方法及并行算法。其中,... 三维网格模型的正演计算是重力资料反演的基础。高精度、高效率的正演有利于提高反演解释质量。针对大尺度、地表观测面研究区的高精度、高效率重力正演问题,本文给出球坐标系三维Tesseroid网格模型重力异常正演方法及并行算法。其中,正演算法采用改进的高斯—勒让德积分法实现大尺度、地表观测面的重力异常高精度计算,并行算法采用基于OpenMP的MATLAB任务并行算法实现高效率计算。理论模型和华东岩石圈三维模型数据试验验证了本文方法的有效性。本文方法为高效的大尺度重力场模拟和三维反演提供技术支撑。 展开更多

关键词 球坐标系 Tesseroid 重力正演 并行计算 地表观测面

下载PDF 职称材料

三维各向异性重力正演模型的高阶展开项评估

8

作者 肖翼忻 霍学深 申文斌 《地球科学前沿（汉斯）》 2018年第3期500-512,共13页

全球范围地表浅层质量分布的重力正演,在计算评估布格异常、地壳均衡补偿、超高精度地球重力场构建及地球内部结构反演时具有重要意义。本文在考虑质量分布的三维各向异性的重力正演方法的基础上,考察GEMMA地壳结构模型下各个大地构造... 全球范围地表浅层质量分布的重力正演,在计算评估布格异常、地壳均衡补偿、超高精度地球重力场构建及地球内部结构反演时具有重要意义。本文在考虑质量分布的三维各向异性的重力正演方法的基础上,考察GEMMA地壳结构模型下各个大地构造分区的重力正演各阶结果对总体重力正演结果的贡献率。以仅考虑重力正演算法的误差因素的影响量级为0.1 mGal精度为判断准则,通过考察在空域下各个大地构造分区按照级数展开阶数的结果与总体重力正演的结果的差值,结果表明,其中地壳厚度较浅的大地构造分区(洋壳、大洋中脊)所需截断阶次为6阶,而地壳厚度较厚的大地构造分区(造山带、盆地、火成岩省、地盾、地台)所需截断阶次为12阶。 展开更多

关键词 重力正演 球谐分析 地壳结构模型 级数展开

下载PDF 职称材料

东营凹陷南部重力异常地质成因 被引量：4

9

作者 胡加山 隋志强 刘成斋 《油气地质与采收率》 CAS CSCD 北大核心 2009年第2期39-42,共4页

引起重力异常的地质因素非常复杂,要对其做出合理地质解释,不仅需要精细分离重力异常,而且需要与其他资料结合。从重力正演模拟出发,分析了可能产生重力异常的地质因素,以及不同地质体的重力响应形态。在此基础上,利用垂向二阶导数及其... 引起重力异常的地质因素非常复杂,要对其做出合理地质解释,不仅需要精细分离重力异常,而且需要与其他资料结合。从重力正演模拟出发,分析了可能产生重力异常的地质因素,以及不同地质体的重力响应形态。在此基础上,利用垂向二阶导数及其他信息增强技术,有效提取东营凹陷南部地区局部重力异常,结合地震、钻井资料,进行了地质成因推断分析,认为影响该区重力局部正异常的地质因素主要是古潜山(包括断块山和残丘山)、盖层构造以及其他如陡坡带砂砾岩体等,古潜山的分布多受大断裂控制,平面上成带连片分布;而区域重力异常则反映了该区的区域构造特征。应用效果表明,在地震成熟勘探区,高精度重力不仅可以发现一些新的地质目标,而且与地震相结合,在诸如深层构造、潜山、有规模的隐蔽性断裂等目标的落实方面,也可以发挥积极作用。 展开更多

关键词 重力正演模拟 重力局部正异常 构造 地质成因 东营凹陷

下载PDF 职称材料

二维重力正演在六盘山盆地西南缘构造样式分析中的应用 被引量：3

10

作者 夏鲁 刘震 +4 位作者 崔红庄 赵海华 相鹏 李廷安 孙迪 《科学技术与工程》 北大核心 2016年第22期183-186,212,共5页

六盘山盆地南部处于旋扭构造的收敛端,地表和地下构造极其复杂,构造样式认识不清成为制约勘探的关键问题。盆地西南缘构造样式问题的焦点在于沙沟断阶带是否存在大规模的逆冲推覆体。针对这一问题,采用二维地震剖面重力正演的方法来分... 六盘山盆地南部处于旋扭构造的收敛端,地表和地下构造极其复杂,构造样式认识不清成为制约勘探的关键问题。盆地西南缘构造样式问题的焦点在于沙沟断阶带是否存在大规模的逆冲推覆体。针对这一问题,采用二维地震剖面重力正演的方法来分析构造样式。首先确定这一地区的密度分层方案;其次,根据两种可能的构造样式以X测线为基础分别建立地质模型,然后进行重力正演模拟;最后将得到的结果与实测数据进行比较。结果表明,与逆冲推覆样式相比,走滑斜冲样式的模拟结果与实测重力剖面更吻合。同时,综合最近的野外构造调查成果和最新地震剖面解释方案,认为六盘山盆地西南缘不存在大规模的逆冲推覆体,主要发育走滑构造。 展开更多

关键词 构造样式 逆冲推覆 走滑斜冲 密度分层方案 地质模型 重力正演模拟

下载PDF 职称材料

井巷环境下微重力异常研究 被引量：1

11

作者 强建科 李俊营 +2 位作者 鲁凯 朱裕振 毛先成 《地球物理学进展》 CSCD 北大核心 2018年第6期2575-2580,共6页

井巷环境相比于地面来说,更容易使物探传感器接近目标异常体,在井巷中进行地球物理勘探,可以提高发现异常的能力.本文采用高精度快速三维重力场模拟算法,研究了钻孔、巷道周围密度异常体的剩余重力异常特征,并分析了钻孔侧向目标体的产... 井巷环境相比于地面来说,更容易使物探传感器接近目标异常体,在井巷中进行地球物理勘探,可以提高发现异常的能力.本文采用高精度快速三维重力场模拟算法,研究了钻孔、巷道周围密度异常体的剩余重力异常特征,并分析了钻孔侧向目标体的产状以及观测点上方或下方异常体的判别方法.模型算例表明:(1)井中微重力方法可以发现井旁一定范围的隐伏密度异常,且重力值明显;(2)井中重力异常极大值与极小值的绝对值大小可以帮助确定异常体倾向,异常绝对值越大,越靠近井轴,否则远离井轴;极值点之间的距离可以确定异常体在深度方向的延伸长度;(3)巷道中下方盈余质量或上方亏损质量会产生相同类型的重力高异常曲线,本文提出通过垂直梯度观测方式来区分异常体的属性,物理依据是:当观测点位于盈余质量之上时,随着观测点靠近目标体,重力异常增大,否则相反.以上研究结果将为未来在复杂井巷环境下开展微重力探测工作提供理论基础. 展开更多

关键词 三维重力正演模拟 井中微重力 巷道微重力 积分方程法

原文传递

基于随机介质模型的二维重力正演

12

作者 刘长扬 罗润林 +5 位作者 李亚南 邓海鹏 彭诚 钟伟 谢淑萍 连参军 《矿产与地质》 2018年第3期528-533,共6页

岩(矿)石密度的空间分布具有随机性和相关性,采用随机过程的谱分解理论和指数型椭圆自相关函数理论,结合江西德兴某重力勘探剖面资料,建立了随机介质模型,并进行重力异常计算。算例表明,采用随机建模的方法可以很好地模拟实际勘探断面... 岩(矿)石密度的空间分布具有随机性和相关性,采用随机过程的谱分解理论和指数型椭圆自相关函数理论,结合江西德兴某重力勘探剖面资料,建立了随机介质模型,并进行重力异常计算。算例表明,采用随机建模的方法可以很好地模拟实际勘探断面的密度结构,据此计算得到的重力异常特征也与实测的基本一致。因此,在重力资料的解释中,根据实际的地质和密度测试资料建立随机介质模型进行正演计算,有助于对实测异常的识别、分析和解释,提高解释的可靠性。 展开更多

关键词 随机介质模型 二维重力正演 岩(矿)石密度

下载PDF 职称材料

基于三维地质建模的点元法重力正演——以某大型铁矿为例

13

作者 周旋 金瑜 许春梅 《物探与化探》 CAS CSCD 2015年第3期490-495,共6页

三维地质建模是目前西方地质学和矿业界通行的一种地质研究和储量计算方法,它运用地质统计学原理,引入区域化变量概念,对矿体三维模型合理划分块段,并进行品位估值和储量计算。以吉尔吉斯斯坦某大型铁矿为研究对象,利用矿山矿石品位与... 三维地质建模是目前西方地质学和矿业界通行的一种地质研究和储量计算方法,它运用地质统计学原理,引入区域化变量概念,对矿体三维模型合理划分块段,并进行品位估值和储量计算。以吉尔吉斯斯坦某大型铁矿为研究对象,利用矿山矿石品位与密度之间的相关函数,将品位估值的矿体三维模型转换为密度三维模型(密度点元),最后通过点元法重力正演计算,并将正演结果与矿区实测重力数据进行对比分析,提出铁矿矿山重力勘探的新思路,为矿山增储及勘探工作部署提出依据。 展开更多

关键词 点元法重力正演 三维地质模型 品位估值 密度品位相关性

下载PDF 职称材料

基于Matlab平台实现二维复杂地质模型多边形网格建模及重力异常正演计算 被引量：3

14

作者 王芃 张忠杰 +4 位作者 张晰 韩颜颜 王敏玲 梁晓峰 滕吉文 《桂林理工大学学报》 CAS 北大核心 2014年第2期254-259,共6页

对已知地质模型进行重力异常正演计算是利用重力场数据对地下介质的密度异常信息进行解释的重要手段。针对三维复杂地质体的建模困难且不便,而且多数三维线性地质体可以近似看作二维地质体以实现快速方便处理的情况,对二维地质体进行多... 对已知地质模型进行重力异常正演计算是利用重力场数据对地下介质的密度异常信息进行解释的重要手段。针对三维复杂地质体的建模困难且不便,而且多数三维线性地质体可以近似看作二维地质体以实现快速方便处理的情况,对二维地质体进行多边形网格建模,并利用Matlab平台编写了适用于复杂二维地质模型的重力异常计算的软件系统。该软件系统对模型设定简单、运算速度快。选择华北地区的剖面使用目前业内广为使用的IGMAS正演软件和文中的软件系统进行了对比试验,结果显示两种方法计算的重力异常结果相当。 展开更多

关键词 重力异常正演 二维密度模型 正演程序 MATLAB平台

下载PDF 职称材料

利用残差地形模型空域法精化局部地球重力场 被引量：1

15

作者 刘聪 王正涛 +1 位作者 张华伟 许智铭 《武汉大学学报（信息科学版）》 EI CAS CSCD 北大核心 2022年第3期369-376,共8页

陆地高分辨率重力数据是超高阶重力场模型及其应用研究的基础,但现有的观测技术和手段限制了陆地重力测量的覆盖区域,全球仍有大量的重力测量空白地区。采用残差地形模型空域法,利用高通滤波技术提取航天飞机雷达地形测绘任务(shuttle r... 陆地高分辨率重力数据是超高阶重力场模型及其应用研究的基础,但现有的观测技术和手段限制了陆地重力测量的覆盖区域,全球仍有大量的重力测量空白地区。采用残差地形模型空域法,利用高通滤波技术提取航天飞机雷达地形测绘任务(shuttle radar topography mission,SRTM)分辨率3″×3″的V4.1数据短波信号,结合EGM2008重力场模型,反演精化1’×1’分辨率的局部地球重力场参量。利用美国地面实测重力观测数据NGS99共57387个点检验基于该方法恢复的1’×1’分辨率的美国局部重力场,相比于EGM2008模型,其误差均方根(root mean square,RMS)由15.470 mGal降低到12.999 mGal,表明该方法可以在保证精度的同时有效地提高陆地重力场的分辨率。将该方法应用于中国华北平原(114°E~117°E,32°N~35°N),计算了该地区1’×1’分辨率的重力扰动和垂线偏差,重力扰动的RMS为25.311 mGal,垂线偏差南北分量和东西分量的RMS分别为2.215″和8.452″,填补了该地区1’×1’分辨率重力测量的空白。 展开更多

关键词 残差地形模型 数字高程模型 重力场正演 陆地重力场

原文传递

基于MPICH环境下的复杂重力异常体并行正演模拟 被引量：2

16

作者 陈石 张健 《地震》 CSCD 北大核心 2013年第4期182-189,共8页

任意多面体重力异常正演公式常用于解决复杂几何形态地质体的正演问题。本文以均匀物性多面体重力异常正演公式为基础,应用有限元技术中的网格离散化思想,以任意四面体为基本单元,通过并行计算技术在MPICH环境下实现了任意连续空间物性... 任意多面体重力异常正演公式常用于解决复杂几何形态地质体的正演问题。本文以均匀物性多面体重力异常正演公式为基础,应用有限元技术中的网格离散化思想,以任意四面体为基本单元,通过并行计算技术在MPICH环境下实现了任意连续空间物性分布复杂异常体网格模型的重力异常正演模拟,通过并行处理可以有效加速正演计算速度。本文研究结果对于联合重力异常场正演建模和开展复杂模型网格的重力场计算有一定参考意义。 展开更多

关键词 重力异常正演 复杂异常体 并行计算 MPICH

下载PDF 职称材料

基于MatlabGUI的多球体重力正演 被引量：1

17

作者 张博汉 张怀榜 《工程地球物理学报》 2016年第5期580-585,共6页

为了更好地解决重力勘探中的正演、反演问题,提高重力勘探资料的解释精度,在MatlabGUI程序开发环境下对地球物理多个球体重力异常正演问题进行了研究与程序编写,重点研究了一个均匀球体、两个及两个以上均匀球体模型的正演模拟,使一些... 为了更好地解决重力勘探中的正演、反演问题,提高重力勘探资料的解释精度,在MatlabGUI程序开发环境下对地球物理多个球体重力异常正演问题进行了研究与程序编写,重点研究了一个均匀球体、两个及两个以上均匀球体模型的正演模拟,使一些较复杂的重力问题能够通过多球体模型近似模拟得以解决。该正演模拟程序界面简洁,操作方便,可交互输入球体模型深度、密度、半径等参数,模拟结果可进行三维、二维等多方式显示,三维曲面可任意调整观测视角,二维曲线可任意选择测线位置进行显示。本程序的模拟结果可以为复杂球体重力勘探资料解释与反演提供参考,也是球体重力问题正演研究的有效手段。 展开更多

关键词 多球体模型 重力异常正演 MatlabGUI.程序模拟

下载PDF 职称材料

典型地质模型2D与3D重力正演模拟对比研究

18

作者 张龙杰 《地球科学前沿（汉斯）》 2015年第4期291-301,共11页

本文典型地质模型三维重力正演模拟研究是在二维重力正演方法总结的基础上,对地质模型作三维重力正演探究。其中所建的地质模型以球体为主。分别导出了二维与三维重力布格异常,重力异常各阶导数公式。在MATLAB环境下实现了模型的二维与... 本文典型地质模型三维重力正演模拟研究是在二维重力正演方法总结的基础上,对地质模型作三维重力正演探究。其中所建的地质模型以球体为主。分别导出了二维与三维重力布格异常,重力异常各阶导数公式。在MATLAB环境下实现了模型的二维与三维重力异常曲线图。结合二维重力正演的分析方法,本文对模型进行了三维重力正演分析,对比讨论了三维重力正演与二维重力正演的优缺点,得到了二维重力正演只能反应地质体垂直切片的重力信息,三维重力正演可以就反映出地质体的整个重力信息,小的局部地质体可以通过高阶导数分析得出其规模形状的结论。 展开更多

关键词 地质模型 重力异常正演 重力异常各阶导数 二维 三维

下载PDF 职称材料

重力和重力梯度数据联合聚焦反演方法 被引量：39

19

作者 秦朋波 黄大年 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2016年第6期2203-2224,共22页

重力数据包含较多的低频信息,重力梯度数据包含较多的高频信息,将重力数据和重力梯度数据进行联合反演得到的结果更加可信.本文基于聚焦反演方法,实现了这一过程.因为联合反演中分量种类增加,所以计算灵敏度矩阵所需要的时间增加,为此,... 重力数据包含较多的低频信息,重力梯度数据包含较多的高频信息,将重力数据和重力梯度数据进行联合反演得到的结果更加可信.本文基于聚焦反演方法,实现了这一过程.因为联合反演中分量种类增加,所以计算灵敏度矩阵所需要的时间增加,为此,本文提出了一种快速计算灵敏度矩阵的方法.因为联合反演对内存的要求增大,本文选择有限内存BFGS拟牛顿法求解反演问题.本文通过再加权的方法实现深度加权.文中利用单一分量的反演结果来预测异常体的埋深信息,随后将埋深信息结合到深度加权函数中,将其用于多分量组合反演计算.给出了模型试验,发现预测得到的异常体的埋深信息与其实际埋深存在偏差,但是将这一信息应用到反演计算,能够得到与真实模型一致的结果.之后,本文通过模型试验来探究重力和重力梯度联合反演的优势,发现将重力和重力梯度数据联合,能够识别出额外的噪声,反演得到的模型更加合理.但是,对于不同分量组合得到的反演结果是相近的,反演模型的提高很小.最后,将联合反演方法应用到美国路易斯安那州Vinton岩丘的实际数据中,结果显示,将重力和重力梯度数据联合反演,反演模型得到了提高,反演得到的结果与地质资料吻合. 展开更多

关键词 重力和重力梯度数据正演 重力和重力梯度联合反演 有限内存BFGS拟牛顿法 深度加权函数 最小梯度支撑函数

下载PDF 职称材料

基于单元细分H-自适应有限元全张量重力梯度正演 被引量：7

20

作者 曹书锦 朱自强 +2 位作者 鲁光银 曾思红 郭文斌 《地球物理学进展》 CSCD 北大核心 2010年第3期1015-1023,共9页

在重力梯度各阶张量正演中,难于求取复杂地质体的解析式,转而利用已知解析解的简单规则几何形体近似剖分地质体,其会引起较大的几何拟合误差.本文首先从全张量重力梯度解析公式出发结合等参变换基于有限元分析给出了满足地下复杂地质体... 在重力梯度各阶张量正演中,难于求取复杂地质体的解析式,转而利用已知解析解的简单规则几何形体近似剖分地质体,其会引起较大的几何拟合误差.本文首先从全张量重力梯度解析公式出发结合等参变换基于有限元分析给出了满足地下复杂地质体正演的型函数,引入体积分的误差指示器用以估计全张量重力梯度误差,提出了适用于全张量重力梯度正演的自适应计算的策略及迭代算法,数值算例证明了本文方法的正确性和有效性. 展开更多

关键词 全张量重力梯度正演 非结构化网格 自适应有限元 复杂形体

下载PDF 职称材料