基于OpenMP的重力张量并行正演

重力张量是重力位的二阶空间导数,对密度体的变化和细节部份反映更为灵敏.对于复杂的重力密度模型来说,张量的解析公式是很难推导的.为了模拟这种情况,将复杂模型进行有限元剖分,计算每个单元的重力异常对测点张量的影响,最后叠加得到整个复杂模型在测点处的重力张量.采用Delaunay四面体非结构化网格对密度体进行剖分,经分析表明,这种剖分方式具有较高的精度.有限元剖分的方式拟合复杂形体缺点之一是计算量比较大,而并行化能有效解决这个问题.这里基于OpenMP并行模型编写了重力张量并行正演程序,分析了不同情况下并行的执行性能,为大规模复杂模型的重力张量正演提供了一种并行策略和思路.     

井中重力矢量数据联合反演研究

井中重力矢量测量可以近距离感知地下目标体,获得地质体不同方向的重力异常特征,提高地下介质的纵向分辨率.这里首先基于点元法实现了井中重力矢量正演,通过模型正演计算,分析了井中三分量重力异常响应特征,然后基于相关性搜索的黄金分割算法实现了井中重力矢量联合反演,通过模型反演分析了井中重力不同分量数据联合反演的效果,验证了反演...     

有限元重力任意复杂地形校正方法研究

提出了一种有限元重力任意复杂地形的校正方法,以解决重力勘探中地形改正问题。利用改进的有限元地形体正演模拟算法,计算任意地形条件下地下密度异常体在地表所引起的重力响应。研究结果表明,在任意复杂地形影响下,经过有限元地形校正后,重力响应曲线可更加直观地反映出地下密度异常体赋存状态,与模型设置相符,表明了该方法的有效性。     

任意密度分布复杂地质体重力异常快速三维正演方法

解决任意密度分布复杂地质体重力异常三维正演快速、高精度计算问题,是实现重力三维反演、人机交互解释建模的关键。针对该问题,从积分方程出发,提出一种波数域重力异常三维正演方法,其关键环节包括三个方面:(1)将研究区域剖分成许多规则小棱柱体,每个小棱柱体密度值可以任意给定,以此刻画任意密度分布和起伏地形条件下的复杂地质体;(2)给出一种新的高精度均匀棱柱体重力异常二维波数域的计算公式,用于计算组合棱柱体模型的重力异常;(3)采用Gauss-FFT法将重力异常从波数域转换到空间域,保证计算效率的同时,有效克服了传统FFT法引起的边界效应问题。模型算例检验结果表明,该算法计算速度快、精度高,对于剖分为百万个棱柱体的模型,耗时只需几秒。     

基于OpenMP的重力张量并行正演CSCD

重力张量是重力位的二阶空间导数,对密度体的变化和细节部份反映更为灵敏。对于复杂的重力密度模型来说,张量的解析公式是很难推导的。为了模拟这种情况,将复杂模型进行有限元剖分,计算每个单元的重力异常对测点张量的影响,最后叠加得到整个复杂模型在测点处的重力张量。采用Delaunay四面体非结构化网格对密度体进行剖分,经分析表明,这种剖分方式具有较高的精度。有限元剖分的方式拟合复杂形体缺点之一是计算量比较大,而并行化能有效解决这个问题。这里基于OpenMP并行模型编写了重力张量并行正演程序,分析了不同情况下并行的执行性能,为大规模复杂模型的重力张量正演提供了一种并行策略和思路。     

基于深度神经网络的重力异常反演

为解决传统线性反演方法容易陷入局部极小,计算效率低等问题,本文提出了一种基于深度学习的重力异常反演方法.该方法首先构造不同形状的二维密度模型,正演得到重力异常,组成数据集;然后用该数据集训练深度神经网络;最后将重力异常数据输入到训练好的深度神经网络,直接得到反演结果.实验结果表明,该方法能快速、准确地反演出地下异常体的...     

中国境内沿北纬40°剖面的重力异常正演拟合及地质解释

在中国境内沿北纬40°截取了长剖面,对布格重力异常和卫星重力异常进行了正演拟合计算,并结合新的大地构造划分作了初步的地质解释。结果表明:古老的陆核、原地台、地台和地块等稳定地区内部的密度值横向变化小,年轻的褶皱带等活动地区内部的密度值横向变化大。各个构造域边界的地壳消减对接带,布格重力异常反映不明显,而卫星重力异常对应较好。布格重力异常与卫星重力异常的正演拟合结果不同的地区,其地壳和上地幔的密度分布有较大的差异。这些重要结论有助于寻求卫星重力异常与布格重力异常的在机制和应用上的共性与差异,以便达到对两种重力异常的本质有更进一步认识的目的。     

重力异常AlexNet深度神经网络反演

针对传统反演方法存在的初始模型依赖、计算时间较长等问题,提出了一种新的基于AlexNet深度神经网络的重力异常反演方法。该方法首先借鉴经典的深度神经网络AlexNet设计了一种用于重力异常反演的Alex反演网络(AlexInvNet),接着设计大量密度异常体模型并通过正演计算得到带标签的数据集,然后用该数据集训练AlexInvNet网络,最后将重力异常数据输入训练好的AlexInvNet网络直接得到反演结果。理论模型反演结果表明,该方法相较于全连接网络深度学习反演方法,能够更好地反演出异常体的位置和密度,具有较好的泛化能力和抗噪声能力。实测数据反演结果表明,该方法能有效解决重力异常反演问题。     

选择法在红旗岭矿区重力G8异常评价中的应用

为评价引起红旗岭矿区G8重力异常的原因,采用解重力勘探反演问题的选择法对其进行重新解释评价。根据剩余重力异常,结合异常地段及其外围地质、地球物理条件建立正演计算模型,进行正演计算并与实测重力异常曲线进行拟合。结果表明G8重力异常并非是由镁铁-超镁铁岩体引起,而是由地层密度局部增大及重力高背景场的局部重力低G11的衬托效应引起。     

用边界元法进行二维重力异常反演

重力异常正演计算的边界元方法是利用格林公式将通常的沿整个密度体的体积分转变为沿密度体边界的积分,这样使问题的维数降低了一维,从而加快了计算速度、提高了计算精度。本文利用边界积分公式解反问题,这样减少了求解参数,提高了稳定性。理论模型计算表明,该法适用于形状复杂均匀密度体的反演,特别是它能用于地下起伏密度界面的确定。     

突泉盆地磁性体重力场剥离的方法技术及岩体三维显示

各地层中赋存的岩体可能形成局部重力异常,对重力资料的密度界面反演计算可靠性产生影响,所以在进行反演计算之前去除岩体引起的异常影响,可增强密度界面反演的可靠性.笔者以突泉盆地重磁异常资料为基础,通过重磁异常特征分析,切取多条剖面进行2.5D正反演计算,进而对磁性岩体形成的重力场进行剥离,增加了重力场密度界面反演计算的可靠性;同时构建岩体的三维地质模型,清晰反映岩体的空间分布,为对目标地质体进行三维地质填图方法做了初步探索.     

复杂条件二维重力场及重力张量场空间波数域正演方法

随着传感器技术的发展,重力场与重力张量场测量技术发展迅速,为实现地下密度分布精细反演提供了数据保障。正演是反演的基础,解决任意密度分布复杂地质体重力场与重力张量正演高效、高精度计算问题,是实现重力高效、精细反演、人机交互反演解释的关键。针对起伏地形和任意密度分布这种复杂条件下二维重力场及重力张量场高效高精度正演问题,这里提出了一种空间波数混合域正演算法,其关键环节包括:①结合新的矩形二度体组合模型波数域表达式和一维Gauss-FFT算法,提出了一种任意密度分布和起伏地形下重力场及重力张量高效、高精度正演算法;②采用新的二维正演算法,计算观测最高点和最低点之间多个不同高度水平网格重力场及重力张量,结合三次样条插值方法,实现了起伏地形上重力场及重力张量场高效、高精度正演。模型算例结果表明,新方法具有高效、高精度的显著特点。     

二度重力人机结合选择法程序

一、功能本程序用于以人机结合选择法解二度重力勘探的反演问题。程序包括正演计算、绘图打印、物质带等效源反演计算三部分。正演计算采用水平多边形棱柱公式,可以有多个地质体,各地质体可以有不同的剩余密度。地形可以是平坦的或起伏的。在每次正演计算前,绘出本次计算所依据的地形地质剖面图。正演计算后,绘出异常剖面图。然后,由用户酌情选     

凤凰─茶陵地学断面重力异常机制分析

对凤凰─茶陵地学断面重力异常,根据爆炸地震测深资料,确定密度分层,用变密度水平梯形体构制剖面密度模型,计算理论重力异常和实测布格重力异常吻合较好。根据剖面密度模型,可以合理地解释布格重力异常。并将地壳划分为沉积盖层、变质岩层、麻粒岩层三层结构,分别计算三层的重力效应曲线及上地幔低密度体重力效应曲线。还计算了莫霍面起伏引起的重力异常曲线,分析了引起地壳重力异常的各种因素及特点。     

凤凰─茶陵地学断面重力异常机制分析

对凤凰─茶陵地学断面重力异常，根据爆炸地震测深资料，确定密度分层，用变密度水平梯形体构制剖面密度模型，计算理论重力异常和实测布格重力异常吻合较好。根据剖面密度模型，可以合理地解释布格重力异常。并将地壳划分为沉积盖层、变质岩层、麻粒岩层三层结构，分别计算三层的重力效应曲线及上地幔低密度体重力效应曲线。还计算了莫霍面起伏引起的重力异常曲线，分析了引起地壳重力异常的各种因素及特点。     

重力异常半定量解释中的若干问题

重力异常的半定量解释,即根据重力异常的特征及工区具体的地质地球物理条件,运用正演概念确定与之有关的地质体的总体形态,产状及分布范围等.目前尚无可供参考的实测数据资料.为此我们想通过对一些简单规则形体的重力异常进行计算,为金属矿区的重力异常半定量解释提供一些参考数据.     

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

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

人工质量体产生引力变化的正演计算

将人工质量体产生的引力以及引力梯度变化的正演计算方法和结论作为精确给定重力和重力梯度变化量的数学基础,可用于测试和标定重力梯度测量装置。与正方形人工质量体相比,球形质量体产生的引力变化量近似于球体全部质量位于球心位置上的点源作用,便于计算和应用。为此,采用球状人工质量体作为引力激励的分析对象,形成正演计算方法。     

重力高阶方向导数的可视化拟合及其应用

针对目前人们所使用的交互式重力剖面拟合软件多以正演布格重力异常为主,在应用效果方面存有欠缺的情况,提出了在作重力剖面布格异常拟合的同时,增加高阶方向导数的拟合。其目的是发挥布格重力异常的高阶方向导数在突出浅部地质体的特征,压制区域性深部地质因素,区分相邻地质体引起的叠加异常等方面的优势。这里采用最小二乘法演算了布格重力异常高阶方向导数的计算过程,对计算程序编制和使用中涉及的相关问题进行了阐述,并例举了在地球物理勘探中,利用高阶方向导数提高重力资料解释效果的实例。并在最后对重力剖面进行可视化拟合解释时应注意的要点做了分析,强调了兼顾布格重力异常及其高阶方向导数二者拟合曲线均基本相似,而不应追求其中某一方面精确拟合的解释原则。     

基于卫星重力异常的南海南部中生界反演

已有钻井以及地震资料揭示南海南部海域广泛发育中生代地层。为了进一步了解南海南部中生代地层的分布特征及油气地质意义,本文利用最新的卫星重力数据反演中生界的深度和厚度。首先,采用重力场正演技术消除海水的影响,获得南海南部布格重力异常。其次,为了消除新生代沉积层的影响,将新生代沉积层划分为0～3、3～6、6～10 km 3层,并根据前人在南海获得的密度与深度的关系,采用变密度Parker正演方法计算新生代沉积层产生的重力影响,并将其从布格重力异常中减去,从而获得前新生代重力异常。在此基础上,采用小波多尺度分解技术,消除深部莫霍面以及局部岩体重力的影响,从而计算得到反映中生界的重力异常。最后采用三维Parker变密度界面快速反演技术获得南海南部中生界深度和厚度。反演结果与已知中生界钻井具有较好的对应。南海南部地区中生界主要分布在礼乐滩、巴拉望岛北部和万安地区,厚度分布小于9 km,其他大部分地区中生界厚度小于1 km或者不存在中生界。其中中生界在礼乐滩地区最为发育,其次在巴拉望岛北部也广泛发育。结合前人在该区域的油气地质条件研究成果,认为南海南部海域礼乐滩地区中生界具有较好的油气勘探前景。