共查询到20条相似文献，搜索用时 831 毫秒
1.
2.
优化算法的选取在很大程度上影响着三维重力反演的计算效率,从而制约着三维重力反演的实用性.在复杂地质构造背景下,不同岩性单元之间可能会发生物性突变,产生尖锐边界.为此,本文提出了一种新的基于柯西分布约束和快速近端目标函数（Fast Proximal Objective Function,FPOF）优化的三维重力反演方法.FPOF优化方法的一个突出特点是在每一步迭代过程中逐一计算剖分网格内的未知密度参数,因此,有较低的计算复杂度和较高的计算效率.此外,目标函数中柯西范数（Cauchy norm）的引入会对反演结果施加稀疏性,有助于产生块状效果.理论模型测试表明,本文方法不仅能产生更加聚焦的反演效果,而且反演所需的时间也比传统的共轭梯度优化方法少.最后将本文方法应用于我国西部某地区实际重力数据,反演结果与已知的地质信息有较好的一致性. 相似文献
3.
考虑到地壳和上地幔密度异常对地表重力的影响,本文提出了修正单纯利用重力资料 反演地壳厚度的算法．该算法利用地震层析成像结果推算出地壳和上地幔密度异常分布,然 后计算由此分布引起的地表重力异常,再从观测布格重力异常中扣除上述密度异常引起的重 力效应,最后在地壳均衡假设的基础上顾及到对反演参考深度的修正,直接反演区域地壳厚 度．本文利用胥颐等地震层析成像的数据,修正了中国西北地区莫霍界面的深度．结果表明, 该区域地壳和上地幔密度异常引起的地表重力异常值为——60times;1O-5m／s2～30times;10-5m／s2． 与直接使用布格重力异常反演地壳厚度相比较,该密度异常对该区域莫霍界面将产生约6 km 的修正．本方法使得利用重力数据反演地壳厚度在理论上更为完备,实际应用中也比较可 行,为更进一步逼近莫霍界面的真实形态提供了一种途径． 相似文献
4.
三维重力反演是地质工作者了解地球深部构造,认知地下结构的重要手段.按照反演单元划分,三维重力反演有离散多面体(Discrete)反演和网格节点(Voxels)反演两种方式.离散多面体反演由于易于吸收先验地质信息得到的理论场能够很好地拟合观测场,因此,在实际重力反演中更受欢迎.目前离散多面体重力反演中初始模型的建立方法繁杂不一,实际应用受到很大的限制.本文本着充分挖掘利用先验信息和重力观测数据得到丰富可靠的反演结果这一原则,以离散多面体反演技术为基础,改进建模过程.在初始模型的建立中,吸收贝叶斯算法优势,采用隐马尔科夫链改善朴素贝叶斯方法的分类效果,通过最大似然函数算法求解,再采取模型降阶技术,固定所建模型中几何体的形态或密度,达到在几何体形态(x,y,z)、密度(σ)和重力值(g)五个参数中降低维数目的,从而减小高维不确定性和正演的计算量,由此反演计算的地质体密度和分布范围相对更准确,更利于重现重力模型结构.通过单位球体和任意形态几何体模拟实验,以及安徽省泥河矿区三维重力反演实践,得到非常接近实际的密度或重力值,大幅提高了三维重力反演的精度和效率,说明该方法是有效、实用的. 相似文献
5.
常规协克里金方法反演重力或重力梯度数据具有抗噪性好、加入先验信息容易等优点,其反演的地下密度分布能够识别异常体中心位置,还原异常体基本形态,但反演图像光滑,分辨率低,这是由于常规方法估计的密度协方差矩阵全局发散、平稳.为了通过协克里金方法获得聚焦的密度分布需要改善密度协方差矩阵的性质.首先,本文推导了理论密度协方差公式,其性质表明,当理论模型聚焦分布时,其密度协方差矩阵是非平稳且聚焦分布的.为了打破常规协方差矩阵全局平稳、发散的特征,本文设置密度阈值处理协方差矩阵,通过不断更新协方差矩阵来迭代实现协克里金反演,最终得到相对聚焦的反演结果.用本文方法处理重力与重力梯度数据恢复两种密度模型,均得到了与正演模型匹配的反演结果;再将方法运用于文顿盐丘的实际测量重力与重力梯度数据,反演结果与已知的地质情况匹配较好. 相似文献
6.
直接迭代法是利用重力异常反演密度界面的一种重要方法,对迭代公式分析的结果表明,可在反演公式中添加加速收敛因子以减少反演的迭代次数从而提高效率.对于平缓的密度界面采用较小的加速收敛因子、而对起伏较大的界面采用较大的加速收敛因子可有效减少反演迭代次数.为此,本文利用重力异常总水平导数判断待反演界面起伏平缓和起伏较大的区域,并以此为基础,将总水平导数计算结果归一化至1.4～2.2之间而形成全区变化的加速收敛因子.假设待反演密度界面上下的密度差沿横向变化,将加速收敛因子加入迭代公式中形成最终的反演迭代公式.模型试算结果证实了本文方法可有效减少直接迭代法的迭代次数,并且也适用于含噪重力数据反演.最后将该方法用于渭河盆地西安凹陷的实际资料处理之中,反演结果证实了该方法的实用性. 相似文献
7.
地震-重力联合反演中的非块状一致性模型 总被引：6，自引：3，他引：6
完善地震-重力联合反演的一个重要措施是建立一致性模型.本文按照地震测深的常用办法,采用二维四边形非块状模型,通过网格节点的密度值来刻划连续性或间断性的物性分布.四边形单元内部的密度分布为双线性函数,本文给出了重力异常的严谨解析式.计算中可采用模型分解技术,将四边形化为矩形与直角三角形的组合,简化程序设计.算法的实施有两种途径:独立单元法和整体模型法.在联合反演的数值模拟中,还推出了剪切模量μ模型以及纵波速度V_P,与密度ρ的分布关系. 相似文献
8.
本文针对蒸汽辅助重力泄油(SAGD)生产中开发监测问题,发展了综合应用地震及重力数据反演储层密度的联合反演算法.通过测井数据建立纵波阻抗与密度的直接关系,并推导出这种关系下重力与纵波阻抗数据联合反演的计算方法,从而计算出蒸汽腔体密度分布规律.文中应用密度反演后验约束正则化方法,采用Tikhonov正则化模型,通过波阻抗数据作为约束进行联合反演,在算法上提高了稳定性,同时得到较高的反演精度.文中对SAGD生产中的理论模型进行了方法试算,并分析了算法的误差,最终应用于SAGD生产的实际数据中,通过最终反演结果分析,该方法取得了很好的应用效果. 相似文献
9.
针对重力异常反演的现状,采用矩形单元体网格剖分的二度半体组合模型,模拟任意的多层密度界面或形体,实现了重力正演.该建模方法适应性强,适合于实际工作中常见的密度横向变化的复杂模型.使用改进的全局寻优的快速模拟退火算法,对重力异常进行反演,结合这种灵活的密度建模方法,反演过程中只需要反演各矩形单元体的密度参数,即可同时得到地质体的界面或形态以及密度值分布.通过复杂界面和形体的模型试验,证明了方法的效果.在江苏下扬子地质构造复杂区的实际应用说明,在地震等先验信息约束下,该重力反演方法提高了反演精度并减少了多解性,可有效解决古生界泥盆系-志留系目的层分布等地质问题. 相似文献
10.
由于目前计算机内存和速度的限制,在对大面积重力资料进行三维密度反演时,已有的反演方法很难奏效.文中提出了一种基于位场分离与延拓的视密度快速反演方法：首先应用场分离的切割法对平面上的重力场进行不同深度层源的切割分离;然后运用大深度向下延拓方法将各层的场延拓至相应的深度;最后反演出各深度层的密度.反演得到的密度是各深度层密度的近似分布,称为视密度反演.本反演方法克服了传统已有方法计算时间特别长、解稳定性差的缺点,在主频1.99 GHz的微机上,反演128×128×10个密度值的计算时间小于20 s,理论模型和实际资料的应用表明反演具有较好的效果. 相似文献
11.
地球物理反演是探索地下结构的最佳途径之一.地震波可以穿透到地球深部进行直接采样,是探测地球深部的主要方法.重力是结构体密度分布与地表观测点之间距离的体积积分效应,重力异常随着源深度的增加衰减很快,其对浅部结构的灵敏度明显优于地震数据.地震和重力联合反演能够相互补充和约束,提高空间分辨率,使反演结果更加稳定可靠.本文首先介绍了联系地震和重力数据的速度-密度经验关系,随后分别介绍了重震联合反演的3种常用方法—顺序反演、同步反演和交叉梯度反演,简要阐述了各种方法在国内外的应用情况.分析认为顺序反演将两类数据分开独立进行计算,原理简单,易于操作实现.但是该方法依赖于先验模型和速度-密度经验关系,可能存在分辨率较低区域(如模型边界)的误差放大效应.同步反演采用将地震和重力数据放在同一方程组中同时反演的策略,减弱了单一数据先验模型对结果的影响,但两种数据的同时运用势必引入数据权重分配问题.交叉梯度寻求不同物理参数模型在结构上的相似性,对潜在的岩石物性关系做了最少的假设,一定程度上降低了反演的非唯一性,但强制性地匹配模型的结构不一定完全符合地下介质的物性分布.因此使用交叉梯度方法反演时应注意模型的推导需要遵循客观标准,以控制模型的结构相似性和数据拟合度.最后指出重震联合反演中的速度-密度经验关系和数据的权重分配仍是值得探究的问题. 相似文献
12.
针对重力与地震联合反演存在的问题,结合已有的研究成果,本文研究实现了速度和密度随机分布共网格单元模型的建模技术,以适应密度和速度剧烈变化的复杂模型及联合反演的计算要求.重力正演利用了该网格的二度半体模型,并进一步改进了地震走时的二维射线追踪计算方法,以适用于速度随机分布的网格介质.结合改进的模拟退火算法,实现了这种共网格条件下的重力与地震资料的同步联合反演.模型试验证明了重力与地震联合反演可以准确确定复杂物性界面的密度和速度结构,适用于物性界面不完全一致和物性变化剧烈的复杂模型,并且联合反演结果要优于单独的重力反演.带先验信息约束下的实际资料的联合反演,进一步证明了该方法的适用性和效果,可提高反演精度并减少多解性. 相似文献
13.
将场源区剖分成长方体单元,通过采集的重磁数据反演出这些单元的密度或者磁化率变化,勾画出场源的分布图像,这种方式是重磁三维反演的重要方向.重磁相关成像通过计算测量的重磁异常与地下各点在测区上的重磁异常的归一化相关,显示出异常地质体的空间赋存状态和等效剩余重磁物性.该方法计算速度快,方法简单、稳定,但是反演的结果只是在-1到+1之间的等效物性,不能够直接反演剩余密度或者磁化率,并且无法引入已知的地质约束.本文通过对物性模型的正演和实测结果的残差进行相关成像,迭代更新物性模型实现对物性参数的反演过程.模型实验证明该方法相对相关成像不仅能提高分辨率,还能够得到真正的物性参数. 相似文献
14.
地质体的重力建模是正确解释和应用重力资料的关键问题之一.针对非规则形状变密度的三度体,本文提出了基于3D Delaunay剖分算法的重力建模方法.采用3D Delaunay剖分算法将三维目标地质体分解为若干变密度四面体体元,推导了基于四面体体元的重力正演公式,建立了剩余密度值与重力异常值的线性方程组;以变密度的长方体和倾斜台阶组合体为例,比较分析了常规块体算法和3D Delaunay剖分算法应用于重力正演的有效性,并采用共轭梯度法加密度约束条件对非规则形状变密度的倾斜台阶组合体进行了密度反演.计算结果验证了本文方法的正确性和有效性.基于3D Delaunay剖分算法的重力建模可应用于存在褶皱、断层、裂缝等复杂地质体的重力正反演计算. 相似文献
15.
16.
利用重力反演沉积盆地基底起伏是油气重力勘探研究的一个重要方面.随着地质目标勘探难度逐渐增大,亟需提高重力反演盆地基底的精细程度.本文研究了基于Lp范数约束的非线性重力反演方法,反演结果能更好地符合裂陷盆地基底起伏特征.反演时假设沉积盆地基底和沉积层均为单一岩性,即沉积层与基底的密度差为常数.将基底之上的沉积层剖分为垂直并置的棱柱体,用棱柱体底面深度近似沉积盆地基底深度,其为待反演的参数.利用Ekblom范数近似Lp范数的表达式,并建立Lp范数意义下的参数向量函数,将其与重力误差函数结合得到非线性反演的目标函数.利用二维裂陷盆地模型进行试算,结果表明当p取0.2时能较准确的反演裂陷盆地基底起伏形态,其应用效果优于基于L0范数和L1范数的约束反演方法.将本文方法应用于南海珠江口盆地实际重力资料反演,将反演的基底深度与地震剖面解释结果进行对比,结果表明该反演方法具有较好的实用性. 相似文献
17.
重力数据包含较多的低频信息,重力梯度数据包含较多的高频信息,将重力数据和重力梯度数据进行联合反演得到的结果更加可信.本文基于聚焦反演方法,实现了这一过程.因为联合反演中分量种类增加,所以计算灵敏度矩阵所需要的时间增加,为此,本文提出了一种快速计算灵敏度矩阵的方法.因为联合反演对内存的要求增大,本文选择有限内存BFGS拟牛顿法求解反演问题.本文通过再加权的方法实现深度加权.文中利用单一分量的反演结果来预测异常体的埋深信息,随后将埋深信息结合到深度加权函数中,将其用于多分量组合反演计算.给出了模型试验,发现预测得到的异常体的埋深信息与其实际埋深存在偏差,但是将这一信息应用到反演计算,能够得到与真实模型一致的结果.之后,本文通过模型试验来探究重力和重力梯度联合反演的优势,发现将重力和重力梯度数据联合,能够识别出额外的噪声,反演得到的模型更加合理.但是,对于不同分量组合得到的反演结果是相近的,反演模型的提高很小.最后,将联合反演方法应用到美国路易斯安那州Vinton岩丘的实际数据中,结果显示,将重力和重力梯度数据联合反演,反演模型得到了提高,反演得到的结果与地质资料吻合. 相似文献
18.
重力反演是恢复地下密度空间分布的有效工具,而选择合理的密度模型约束方法是提升重力反演分辨率和可靠性的关键.常规约束方法大多是从剖分网格空间中的密度模型出发,通过调整光滑或稀疏约束权重来匹配反演目标,但当地质体类型多样、异常分离不准确及网格剖分方案不合理时,模型约束的合理性与灵活性难以得到有效保证.为此,本文提出了一种基于密度模型稀疏表征的重力反演方法.首先假设待反演的密度模型表征为模型特征矩阵和稀疏分解系数的线性组合,之后重新推导了重力反演目标函数,并给出了分解系数的稀疏求解过程.相比现有重力反演方法,用于构建模型特征矩阵的特征模型可包含不同类型地质体的先验几何信息,分解系数的稀疏性保证了待反演目标来自于最典型的地质模式组合.最后,通过模型试验及实际资料验证了基于密度模型稀疏表征的重力反演方法的有效性. 相似文献
19.
重力数据的密度界面反演是位场数据解释中的一项主要工作,在区域构造演化、深部莫霍面确定等领域的研究中发挥重要作用.近年来,数据驱动的深度学习方法广泛地应用在地球物理数据处理与反演中,本文提出一种基于深度学习U-net网络的重力数据密度界面反演方法.首先,对半椭球体界面模型进行随机抽取和组合进而形成地下起伏界面数据集,并基于Parker正演理论对界面数据集进行重力异常正演计算,为深度学习网络模型的训练提供特征完备的数据源;其次,设计了基于U-net网络模型的深度学习界面反演算法,在传统的损失函数基础上增加光滑损失项和过拟合抑制项,提高重力界面反演结果的光滑性和收敛效率;最后通过测试样本集进行反演预测,验证建立深度学习网络模型的泛化性.本文通过理论模型和实际数据试验分析了本文方法在密度界面反演中的有效性和实用性,基于改进损失函数约束的深度学习界面反演方法有效地提高了密度界面反演的收敛效率和计算稳定性. 相似文献
20.
三维密度反演已经成为重力数据定量解释的常规方法,但由于重力数据本身并没有深度分辨率,为了减少由此引起的重力反演的非唯一性,常用的手段是引入额外的先验信息.本文提出了一种重力三维稀疏反演（以下简称稀疏反演）方法,该方法通过求解物性上下界约束时的Lp范数（0 ≤ p ≤ 1）稀疏优化问题,来获得具有尖锐边界的解.与传统的L2范数反演方法相比,稀疏反演方法可以更加有效地利用已知的物性信息,获得深度分辨率更高的反演结果.此外,我们也分析了稀疏反演方法与二值、三值反演算法的等价性以及在实际应用中需要注意的问题.最后,通过模型试验以及矿区实测数据反演验证了稀疏反演方法的有效性. 相似文献