二维大地电磁尖锐边界反演研究

尖锐边界反演是大地电磁及其他地球物理反演中的一个较难解决的问题,目前有不少专家正对这方面进行研究.我们在分析了Zhdanov（2004）提出的基于最小支撑泛函聚焦反演方法的基础上,在目标函数中引入对角梯度支撑以改善倾斜电性分界面的反演效果.同时,为了提高计算效率,本文还使用互易定义计算灵敏度矩阵及使用正则化共轭梯度法优化目标函数.通过模型试验并与未引入对角梯度支撑的反演结果进行了比较,表明本文的方法反演效果更好而且算法稳定和实用.     

三维重力快速反演的优化算法（英文）

实用化的重力数据三维反演过程中需要大量的计算时间和存储空间,针对此问题本文提出一种组合优化算法:(1)利用空间域的精确性和波数域的快速性,提出一种高精度快速三维正演算法;(2)利用反演矩阵中的对称性,将重力共轭梯度反演(CG)中的主要计算分解成两次正演计算,从而优化了三维重力反演的计算效率。通过不同网格剖分数量模型的数值实验,验证了本文优化算法的计算精度和计算效率。     

基于模型降阶的贝叶斯方法在三维重力反演中的实践

三维重力反演是地质工作者了解地球深部构造,认知地下结构的重要手段.按照反演单元划分,三维重力反演有离散多面体(Discrete)反演和网格节点(Voxels)反演两种方式.离散多面体反演由于易于吸收先验地质信息得到的理论场能够很好地拟合观测场,因此,在实际重力反演中更受欢迎.目前离散多面体重力反演中初始模型的建立方法繁杂不一,实际应用受到很大的限制.本文本着充分挖掘利用先验信息和重力观测数据得到丰富可靠的反演结果这一原则,以离散多面体反演技术为基础,改进建模过程.在初始模型的建立中,吸收贝叶斯算法优势,采用隐马尔科夫链改善朴素贝叶斯方法的分类效果,通过最大似然函数算法求解,再采取模型降阶技术,固定所建模型中几何体的形态或密度,达到在几何体形态(x,y,z)、密度(σ)和重力值(g)五个参数中降低维数目的,从而减小高维不确定性和正演的计算量,由此反演计算的地质体密度和分布范围相对更准确,更利于重现重力模型结构.通过单位球体和任意形态几何体模拟实验,以及安徽省泥河矿区三维重力反演实践,得到非常接近实际的密度或重力值,大幅提高了三维重力反演的精度和效率,说明该方法是有效、实用的.     

基于数据空间和稀疏约束的三维重力和重力梯度数据联合反演

随着重力和重力梯度测量技术的日趋成熟,基于重力和重力梯度数据的反演技术得到了广泛关注.针对反演多解性严重、计算效率低和内存消耗大等难点问题,本文开展了三维重力和重力梯度数据的联合反演研究,该方法结合重力和重力梯度两种数据,将L0范数正则化项加入到目标函数中,并在数据空间下采用改进的共轭梯度算法求解反演最优化问题.同时,本文摒弃了依赖先验信息的深度加权函数,引入了自适应模型积分灵敏度矩阵,用来克服因重力和重力梯度数据核函数随深度增加而衰减引起的趋肤效应问题.为了提高反演计算效率,本文又推导出基于规则网格化的重力和重力梯度快速正演计算方法.模拟试算表明,改进的共轭梯度法可以降低反演的迭代次数,提高反演的收敛速度;自适应模型积分灵敏度矩阵,可以有效解决趋肤效应,提高反演纵向分辨能力;数据空间和改进的共轭梯度算法结合,可以更好地降低反演求解方程的维度,避免存储灵敏度矩阵,有效地降低反演计算时间和内存消耗量.野外实例表明,该算法可以在普通计算机下快速地获得地下密度分布模型,表现出较强的稳定性和适用性.     

基于优化有限差分和混合吸收边界条件的三维VTI介质声波和弹性波数值模拟

传统有限差分系数是通过泰勒级数展开求取的,这样导致所计算的频散曲线在大波数区域会产生较强的数值误差.针对二阶空间偏导数的显式有限差分离散,本文发展了一种新的优化差分系数方法：首先将泰勒级数展开与多点采样方法结合应用于空间频散关系,基于最大范数建立直观有效的优化目标函数,采用Remez算法求解该目标函数,从而获得最优化差分系数.利用优化有限差分方法求解三维垂直对称轴横向各向同性（VTI）介质中的声波和弹性波方程.另外,本文将二维混合吸收边界条件推广到三维VTI介质中,用于吸收人工截断边界反射;基于各向异性特征,合理调整了边界区域的速度值来提高吸收效果.考虑到三维情况下计算效率的问题,本文波场外推过程中采用图形处理器（GPU）取代传统的中央处理器（CPU）.数值精度分析表明,相比较于传统的泰勒级数展开方法,优化有限差分方法在大波数区域对频散误差的压制效果更明显.在三维均匀和修改的Hess VTI模型中的数值模拟实验证明了本文方法具有更高的精度与效率,混合吸收边界条件在三维VTI介质中具有良好的边界吸收效果.     

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

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

基于GPU并行的重力、重力梯度三维正演快速计算及反演策略

利用NVIDIA CUDA编程平台,实现了基于GPU并行的重力、重力梯度三维快速正演计算方法.采用当前在重力数据约束反演或联合反演中流行的物性模型(密度大小不同、规则排列的长方体单元)作为地下剖分单元,对任意三维复杂模型体均可用很多物性模型进行组合近似,利用解析方法计算出所有物性模型在计算点的异常值并累加求和,得到整个模型体在某一计算点引起的重力(或重力梯度)值.针对精细的复杂模型体产生的问题,采用GPU并行计算技术,主要包括线程有效索引与优化的并行归约技术进行高效计算.在显卡型号为NVIDIA Quadro 2000相对于单线程CPU程序,重力和重力梯度U_xx、U_xy正演计算可以分别达到60与50倍的加速.本文还讨论了GPU并行计算在两种反演方法中的策略,为快速三维反演技术提供了借鉴.     

基于物性参数梯度累加约束的电阻率法和重力勘探数据二维联合反演

对不同地球物理方法的数据进行联合反演,通过模型参数的相互约束可以减小反演的多解性.本文对电阻率法和重力勘探数据开展了基于电阻率和剩余密度梯度累加约束的二维联合反演研究.电阻率法采用数据空间Occam反演,重力勘探采用基于对数障碍法的正则化反演.通过在电阻率法和重力勘探反演的目标函数中引入电阻率和剩余密度梯度累加约束项,利用交替迭代方式实现了电阻率法和重力勘探数据的二维联合反演.对理论模型合成数据进行了单方法反演和联合反演,对比了反演效果.结果表明：联合反演结果优于单方法反演结果,联合反演结果对异常体的物性值恢复和形态刻画效果更好.     

基于深度学习U-net网络的重力数据界面反演方法

重力数据的密度界面反演是位场数据解释中的一项主要工作,在区域构造演化、深部莫霍面确定等领域的研究中发挥重要作用.近年来,数据驱动的深度学习方法广泛地应用在地球物理数据处理与反演中,本文提出一种基于深度学习U-net网络的重力数据密度界面反演方法.首先,对半椭球体界面模型进行随机抽取和组合进而形成地下起伏界面数据集,并基于Parker正演理论对界面数据集进行重力异常正演计算,为深度学习网络模型的训练提供特征完备的数据源;其次,设计了基于U-net网络模型的深度学习界面反演算法,在传统的损失函数基础上增加光滑损失项和过拟合抑制项,提高重力界面反演结果的光滑性和收敛效率;最后通过测试样本集进行反演预测,验证建立深度学习网络模型的泛化性.本文通过理论模型和实际数据试验分析了本文方法在密度界面反演中的有效性和实用性,基于改进损失函数约束的深度学习界面反演方法有效地提高了密度界面反演的收敛效率和计算稳定性.     

基于Lp范数(0＜p＜0.5)约束的裂陷盆地基底非线性重力反演

利用重力反演沉积盆地基底起伏是油气重力勘探研究的一个重要方面.随着地质目标勘探难度逐渐增大,亟需提高重力反演盆地基底的精细程度.本文研究了基于Lp范数约束的非线性重力反演方法,反演结果能更好地符合裂陷盆地基底起伏特征.反演时假设沉积盆地基底和沉积层均为单一岩性,即沉积层与基底的密度差为常数.将基底之上的沉积层剖分为垂直并置的棱柱体,用棱柱体底面深度近似沉积盆地基底深度,其为待反演的参数.利用Ekblom范数近似Lp范数的表达式,并建立Lp范数意义下的参数向量函数,将其与重力误差函数结合得到非线性反演的目标函数.利用二维裂陷盆地模型进行试算,结果表明当p取0.2时能较准确的反演裂陷盆地基底起伏形态,其应用效果优于基于L0范数和L1范数的约束反演方法.将本文方法应用于南海珠江口盆地实际重力资料反演,将反演的基底深度与地震剖面解释结果进行对比,结果表明该反演方法具有较好的实用性.     

重磁数据三维物性反演方法进展

综述了重磁数据三维物性反演方法中的几个关键问题.主要包括正演快速算法、反演框架、约束因子讨论、反演算法实现等方面.正演快速算法主要讨论了等效存储几何格架技术、基于GPU加速的并行计算技术以及小波压缩技术.三维物性反演则是在最小二乘意义下使目标函数达到极小的线性或非线性反演.指出,对于特定地质问题需要谨慎选择不同且合适的约束方法乃至反演算法,才能达到好的效果.最后讨论了重磁数据三维物性反演较好的应用前景及发展方向.     

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

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

基于3DDelaunay剖分算法的重力建模与分析

地质体的重力建模是正确解释和应用重力资料的关键问题之一.针对非规则形状变密度的三度体,本文提出了基于3D Delaunay剖分算法的重力建模方法.采用3D Delaunay剖分算法将三维目标地质体分解为若干变密度四面体体元,推导了基于四面体体元的重力正演公式,建立了剩余密度值与重力异常值的线性方程组;以变密度的长方体和倾斜台阶组合体为例,比较分析了常规块体算法和3D Delaunay剖分算法应用于重力正演的有效性,并采用共轭梯度法加密度约束条件对非规则形状变密度的倾斜台阶组合体进行了密度反演.计算结果验证了本文方法的正确性和有效性.基于3D Delaunay剖分算法的重力建模可应用于存在褶皱、断层、裂缝等复杂地质体的重力正反演计算.     

基于密度模型稀疏表征的重力反演方法

重力反演是恢复地下密度空间分布的有效工具,而选择合理的密度模型约束方法是提升重力反演分辨率和可靠性的关键.常规约束方法大多是从剖分网格空间中的密度模型出发,通过调整光滑或稀疏约束权重来匹配反演目标,但当地质体类型多样、异常分离不准确及网格剖分方案不合理时,模型约束的合理性与灵活性难以得到有效保证.为此,本文提出了一种基于密度模型稀疏表征的重力反演方法.首先假设待反演的密度模型表征为模型特征矩阵和稀疏分解系数的线性组合,之后重新推导了重力反演目标函数,并给出了分解系数的稀疏求解过程.相比现有重力反演方法,用于构建模型特征矩阵的特征模型可包含不同类型地质体的先验几何信息,分解系数的稀疏性保证了待反演目标来自于最典型的地质模式组合.最后,通过模型试验及实际资料验证了基于密度模型稀疏表征的重力反演方法的有效性.     

基于Extrapolation Tikhonov正则化算法的重力数据三维约束反演

通过研究重力数据三维反演解的病态性,利用基于拉格朗日插值方法的Extrapolation Tikhonov正则化方法来解决反演中解的不唯一性和不稳定性问题,该方法最大限度的减少了因正则化参数的引入而在反演结果中介入的误差,同时详细讨论了基于三种选择原则的正则化双参数的具体选择方法,模型试算结果表明,与原Tikhonov方法相比,该方法提高了反演的拟合精度.其次,为了消除核函数随深度增加而快速衰减对反演结果的影响,本文改进了前人的重力数据三维反演深度加权函数,改进后的加权函数与原函数相比能更好的识别异常体底部密度分布特征,对于埋深较深的异常体具有较好的识别效果,更好的解决了由近地面趋肤效应作用引起的密度分布不均的问题.同时,利用上下限约束函数限制每一个立方体的密度差范围,并应用于多组人工合成模型.结果表明:该反演方法能准确地获得正演模型的预设参数范围和位置.     

二度半长方体组合模型的重力模拟退火反演

针对重力异常反演的现状,采用矩形单元体网格剖分的二度半体组合模型,模拟任意的多层密度界面或形体,实现了重力正演.该建模方法适应性强,适合于实际工作中常见的密度横向变化的复杂模型.使用改进的全局寻优的快速模拟退火算法,对重力异常进行反演,结合这种灵活的密度建模方法,反演过程中只需要反演各矩形单元体的密度参数,即可同时得到地质体的界面或形态以及密度值分布.通过复杂界面和形体的模型试验,证明了方法的效果.在江苏下扬子地质构造复杂区的实际应用说明,在地震等先验信息约束下,该重力反演方法提高了反演精度并减少了多解性,可有效解决古生界泥盆系-志留系目的层分布等地质问题.     

基于抛物线密度模型的频率域三维界面反演及其在川滇地区的应用

密度界面反演作为了解地球内部结构的一种重要方法,长期以来都是重力学研究的主要内容.本文结合抛物线密度模型及频率域算法的优点,将抛物线密度函数应用于Parker-Oldenburg算法,经过理论推导得到了抛物线密度模型的频率域公式,从而建立了基于抛物线密度模型的三维密度界面重力异常正反演的算法和流程.理论模型数据试验表明本方法快速、有效,适用于大多数浅部比深部增加更快的实际地壳密度.研究中还利用该方法对川滇地区重力异常进行了反演,获得了该区的莫霍面深度分布,并与接收函数研究结果进行对比分析,进一步验证了本文方法的正确性和有效性.     

New advances in regularized inversion of gravity and electromagnetic data

The interpretation of potential and electromagnetic fields observed over 3D geological structures remains one of the most challenging problems of exploration geophysics. In this paper I present an overview of novel methods of inversion and imaging of gravity and electromagnetic data, which are based on new advances in the regularization theory related to the application of special stabilizing functionals, which allow the reconstruction of both smooth images of the underground geological structures and models with sharp geological boundaries. I demonstrate that sharp-boundary geophysical inversion can improve the efficiency and resolution of the inverse problem solution. The methods are illustrated with synthetic and practical examples of the 3D inversion of potential and electromagnetic field data.