重磁异常三维物性反演随机子域法方法技术

本研究针对三维物性反演中存在的困难和问题,提出三维物性反演的随机子域方法技术,首先是将正反演中保持不变的几何格架分离计算并存储,避免重复计算,从而提高正反演计算速度;其次是利用对称性等实现等效计算,明显降低格架计算和存储要求;再通过随机子域方式,降低反演的维数问题;另外,通过概率方式控制子域生成的分布,实现约束新机制. 模型和实例计算表明了方法技术的效果,为大面积重磁数据的三维反演提供了有效的途径.     

频率域航空电磁数据的加权横向约束反演

传统的一维反演技术已经被广泛应用于航空电磁数据解释中.然而,利用单点水平层状介质模型模拟地下复杂地电结构有时会遇到困难.突出表现在反演参数的横向不连续性,即使相邻测点的反演结果也会出现突变.本文针对航空电磁直升机吊舱系统可进行密集采样,相邻测点地下电性结构应具有某种程度连续性的特点,研究航空电磁数据横向约束反演理论,并提出参数加权约束方法.首先阐述频率域航空电磁正演和加权横向约束反演理论,着重介绍这种拟二维反演方法的基本原理和实施步骤,以及将该方法成功应用于频率域航空电磁数据反演处理的方法技术.最后,通过对理论和实测数据反演处理,并与传统的一维反演结果进行对比,验证加权横向约束反演方法的有效性.     

晋东北及邻区重力三维正、反演的初步结果

本文根据地质和区域重力资料,利用重力三维正、反演方法研究得到了晋东北地区地壳深部构造。通过三维正演,得到了不同深度的剩余重力异常;用一个已知的莫霍界面的深度值,反演出了全区地壳厚度的平面分布。在此基础上,对一些活动断裂和构造的深度以及地震和第四纪火山活动的地壳深部构造背景进行了探讨。研究还发现,沿北纬40°线可能存在一条深部剪切构造带。     

基于剩余异常相关成像的重磁物性反演方法

将场源区剖分成长方体单元,通过采集的重磁数据反演出这些单元的密度或者磁化率变化,勾画出场源的分布图像,这种方式是重磁三维反演的重要方向.重磁相关成像通过计算测量的重磁异常与地下各点在测区上的重磁异常的归一化相关,显示出异常地质体的空间赋存状态和等效剩余重磁物性.该方法计算速度快,方法简单、稳定,但是反演的结果只是在-1到+1之间的等效物性,不能够直接反演剩余密度或者磁化率,并且无法引入已知的地质约束.本文通过对物性模型的正演和实测结果的残差进行相关成像,迭代更新物性模型实现对物性参数的反演过程.模型实验证明该方法相对相关成像不仅能提高分辨率,还能够得到真正的物性参数.     

基于方向tilt-Euler的三维磁数据快速反演

三维磁异常快速反演是磁法勘探进行地质解释的一项重要内容.本文在方向tilt梯度基础上进行导数处理,并进行合理改进,得到了6个改进的方向tilt梯度导数,丰富了不同方向上的磁异常信息.对常规欧拉反褶积方程进行三个方向求导,所构成的方程组与6个改进的方向tilt梯度导数进行结合,从而构建出了两个基于方向tilt梯度导数的快速反演方程组（方向tilt-Euler法）,用于反演场源的位置参数.另外,将三个方向导数下的欧拉反褶积进行平方和及均方根处理,得到了场源构造指数的求解公式.单一模型试验表明,6个改进的方向tilt梯度导数异常基本不受磁化方向影响,可从不同方向上反映磁异常信息,方向tilt-Euler能够准确地反演出场源的位置及构造指数.叠加模型试验证实了方向tilt-Euler法相对于tilt-Euler法,具有反演解聚集度更强、连续性更好、反演准确度更高、受噪声影响更小及无虚假反演解等优点.将本文方法应用于内蒙古塔木素航磁资料处理中,获得了丰富的地下磁源分布信息,为了解该地区隐伏断裂构造及岩体分布等提供了有效依据.

     

基于全卷积神经网络的磁异常及磁梯度异常反演

地球物理反演问题具有病态性、不适定性,传统的线性反演方法面临着次优逼近和初始模型选择等挑战,为了提高磁场数据反演的精度,受深度学习卓越的非线性映射能力的启发,本文提出了一种基于全卷积神经网络的磁异常及磁梯度异常反演方法.文中首先提出了一种基于网格点几何格架的磁异常及磁梯度异常的空间域快速正演算法,这为本文全卷积神经网络...     

基于地质体空间位置优化约束的航空重力梯度数据三维物性反演

重力数据的物性反演面临着严重的多解性问题,降低多解性的有效手段是加入约束条件.而边界识别、深度估计及成像方法可获取地质体的水平位置、深度范围等几何参数信息,本文将基于数据本身挖掘的地质体几何参数信息约束到物性反演中,以降低反演的多解性.通过引入基于深度信息的深度加权函数及基于水平位置的水平梯度加权函数建立优化约束条件,有效地提高了反演结果的横向及纵向分辨率.重力梯度数据包含更多的地质体空间特征信息,将优化约束反演方法应用到全张量数据的反演中,模型试验表明本文方法反演结果与理论模型更加吻合.最后对美国路易斯安那州文顿盐丘实测航空重力梯度数据的应用表明,本文方法在其他地球物理、地质资料不足的情况下获得更可靠的反演结果.

     

多层电磁波测井反演

针对轴对称二维非均匀介质地层,以两接收线圈的相位差和幅度比为附加条件,建立了反演地层介电常数和电导率的电磁波测井反演数学模型．针对反问题的不适定性,利用正则化方法建立了反演的选代法．为克服层位较多时整体反演在存贮量、计算量等方面的问题,考虑了一种每次反演一个目的层的多层反演方案,并针对6层和10层模型进行了数值模拟．模拟结果表明,以相位差和幅度比为测井响应来反演介质电性参数是可行的,并适合于薄层解释,但对多层模型不能采用整体反演．     

1D SNMR吉洪诺夫正则化反演方法研究

地面核磁共振(SNMR)技术是目前世界上直接用来寻找地下水的技术。本文在讨论反演基本问题的基础上,用共轭梯度法实现了核磁共振的一维正则化反演。通过反演理论数据、噪声数据和实测数据,说明该方法的可靠性。     

数据空间磁异常模量三维反演

强剩磁的存在通常导致了总磁化强度方向未知,进而影响了磁异常的反演和解释.磁异常模量是一种受磁化方向影响小的转换量,可以在强剩磁条件下通过反演三维磁化强度大小分布来推测场源分布状态.我们提出了一种数据空间磁异常模量反演算法来减少剩磁的影响.与标准的模型空间L2范数正则化反演方法相比,我们的方法有两个优点:一是无需搜索正则化参数(需要反复求解非线性反演问题),因而可以减少计算时间;二是反演结果更加聚焦,深度分辨率更高,我们对此进行了原因分析.通过模型和实测数据测试证明了该算法的有效性和更好的反演效果.     

重力和重力梯度数据联合聚焦反演方法

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

九瑞矿集区重磁三维约束反演及深部找矿意义

研究九瑞矿集区成矿地层和岩体的三维分布特征可为区域成矿背景和成矿规律研究提供新的信息,实现研究区深部及外围找矿突破.本文首先综合分析研究区区域地质及地层和岩石密度和磁化率特征,然后将区域地质和12条地质剖面信息利用自行开发的GIF2UBCmodel程序构建剩余密度和磁化率参考模型和上下边界约束模型,实现了带约束重磁数据三维反演.反演结果揭示九瑞矿集区剩余密度和磁化率三维分布特征,解译了矿集区主要地层和岩体分布特征.该区存在邓家山—东雷湾以及武山—丁家山两个大型岩体;盖层褶皱与结晶基底隆起方向一致,是深部基地隆起的继承.通过与矿集区已知矿床对比发现,高磁性岩体边缘和基底隆起(高密度体)两翼为成矿有利区.     

重力数据3D密度成像中EXTR方法的各参数变化对反演结果的影响(英文)

Tikhonov正则化(TR)方法在重磁数据处理中发挥了重要的作用,本文在研究如何利用Tikhonov正则化方法方法解决重力数据3D反演的同时,深入讨论了可进一步提高拟合误差的Extrapolation Tikhonov正则化方法(EXTR)的原理,并就其参数选择方法及各参数对拟合误差、迭代次数及反演结果的影响进行研究。常密度及变密度组合模型试算结果表明,与TR方法相比,EXTR方法不仅可以达到解释人员设定的先验拟合误差水平,在计算时间及迭代次数相应增加的前提下有更高的拟合精度;同时其反演结果也更加紧致,进一步改善了TR反演结果的发散性;并且其反演数据范围更贴近预设模型参数范围,模型特征与预设模型密度分布吻合较好。     

基于共轭梯度算法的重力梯度数据三维聚焦反演研究

重力梯度数据相对于传统重力数据,能够更细致、准确地描述地球浅部构造和研究矿产资源分布等信息.本文采用共轭梯度算法,在加权密度域求解重力梯度数据三维聚焦反演最优化问题,以恢复地下三维密度分布,目标函数包括数据不拟合函数和最小支撑稳定函数.首先,在推导目标函数对加权密度的一阶导数时,为了得到更合理的计算公式,我们考虑变加权函数中含有密度变量;此外,本文通过密度上下限约束,改善了传统聚焦反演中聚焦因子选取困难的问题.新算法获得的反演结果,对聚焦因子的选择约束较少,相比传统聚焦算法,能够更容易的获得理想结果.将方法应用于理论模型验证其有效性和正确性,并应用本文方法处理文顿盐丘地区的航空全张量重力梯度数据,得到了与已知地质信息匹配的密度分布,表明本文方法具有处理实际数据的能力.

     

大地电磁三维快速松弛反演并行算法研究(英文)

为了进一步提高大地电磁三维快速松弛反演的计算效率,在深入分析大地电磁三维快速松弛反演算法的基础上,结合MPI自身的优越性,确定了并行计算的思路,实现了三维快速松弛反演的并行计算。通过理论模型合成数据和实测数据对实现的三维快速松弛反演并行程序进行了试算,分析对比了在多种情况下程序的执行效率。测试结果表明,所实现的三维快速松弛反演并行程序运行结果正确,效率提高明显。此成果在普通微机上实现,推动了大地电磁三维反演技术的实用化,可为其它地球物理三维正反演研究所借鉴。     

基于局部相关性约束的三维大地电磁数据和重力数据的联合反演

为解决地球物理反演中多解性的问题,综合多种地球物理信息的联合反演受到了广泛的关注.本文依据不同地球物理响应可能由相同异常体引起,而不同地球物理分布参数之间存在相关性等特点,提出了一种基于局部Pearson相关系数约束的联合反演方法.该方法假设每个局部区域模型参数的分布具有线性相关特性,在拟合不同类型观测数据时,对局部模型参数施加相关性约束,进行联合反演以减少多解性.本文采用交替迭代联合反演流程,改善了同一目标函数下联合反演收敛性和速度问题.基于新的联合反演方法和流程,我们测试了三维大地电磁和重力仿真数据的联合反演.结果表明,本文提出的基于局部相关性约束的联合反演方法,能充分利用大地电磁和重力观测数据信息,有效改善单一地球物理反演收敛性和多解性的问题,反演效果得到明显提升.

     

基于交叉梯度约束的重力、磁法和大地电磁三维联合反演

为了降低单一地球物理方法反演的多解性及受噪声的影响程度,本文围绕重力、磁法和大地电磁法开展了三维联合反演的研究.重、磁采用基于对数障碍法的正则化反演算法,大地电磁使用limited-memory Broyden-Fletcher-Goldfarb-Shanno（L-BFGS）反演算法,引入交叉梯度函数实现了三种物性结构的相互耦合,最终开发出一套重磁电三维联合反演算法,并实现MPI并行加速计算.通过理论模型算例验证了算法的准确性,结果表明：不论是单棱柱体模型还是组合棱柱体模型,联合反演结果相较单独反演对于异常体的空间形态刻画以及物性数值恢复具有较好的提升;单棱柱体模型算例使得异常体的物性参数（密度、磁化率和电阻率）更加接近于真实的物性参数;组合棱柱体模型的联合反演结果不仅仅消除了围岩物性参数的假异常,而且还增强了异常体边界结构的恢复程度.

     

基于Lp范数稀疏优化算法的重力三维反演

三维密度反演已经成为重力数据定量解释的常规方法,但由于重力数据本身并没有深度分辨率,为了减少由此引起的重力反演的非唯一性,常用的手段是引入额外的先验信息.本文提出了一种重力三维稀疏反演（以下简称稀疏反演）方法,该方法通过求解物性上下界约束时的L_p范数（0 ≤ p ≤ 1）稀疏优化问题,来获得具有尖锐边界的解.与传统的L₂范数反演方法相比,稀疏反演方法可以更加有效地利用已知的物性信息,获得深度分辨率更高的反演结果.此外,我们也分析了稀疏反演方法与二值、三值反演算法的等价性以及在实际应用中需要注意的问题.最后,通过模型试验以及矿区实测数据反演验证了稀疏反演方法的有效性.

     

重力与地震资料的模拟退火约束联合反演

联合反演是综合地球物理研究的重要定量解释手段.本文在总结和分析重力与地震资料联合反演的研究现状基础上,利用改进的全局寻优的快速模拟退火算法,实现了重力和地震资料的约束同步联合反演.针对性地设计了密度和速度界面不完全一致的模型,理论模型的试验说明了方法的效果和适用性.结合最近完成的广东徐闻地区实际资料的处理和解释,表明该方法可准确确定复杂构造物性界面的密度和速度结构,在该地区的油气勘探中发挥了作用.在先验信息约束下,该联合反演方法要明显优于单独的重力反演.