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

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

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

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

基于对称多项式约束的重力反演

由于地球物理场的等效性、观测数据集的有限性和数据观测存在误差等问题,地球物理反演通常为不适定问题.Tikhonov正则化是应用最广泛的减少反演不适定性的方法之一,稳定函数的引入一方面提高了反演稳定性;另一方面将先验信息融入反演,提高了反演的效果.稳定函数的构建是正则化反演的关键之一,本文提出一种可以融入任意数量先验物性信息的对称多项式模型约束稳定函数.通过引入该函数,达到充分利用岩石先验物性信息、提高反演效果的目的.重力模型及实测数据反演表明,基于对称多项式的模型约束反演,可以实现任意数量岩石统计物性的利用,并获得理想的异常体边界和密度范围约束效果.

     

基于Tikhonov正则化方法的背景噪声源反演

在平面波假设下,基于Tikhonov正则化方法对背景噪声源能量随方位角的分布进行反演;从数据分辨率矩阵、模型分辨率矩阵和单位协方差矩阵的角度对解进行了评价;通过数值实验,分别研究了速度模型、噪声源能量分布模型、周期对解的估计的影响;从不同的角度,说明了Tikhonov正则化方法在背景噪声能量随方位角分布的反演中是一种有...     

基于Lp范数正则化的V字型密度界面重力反演

V字型密度界面是一类常见的密度界面,如海沟、半地堑以及俯冲带之下的莫霍面,利用重力数据刻画此类密度界面形态对于区域构造研究、油气勘探以及物理海洋学等都具有重要意义.本文首先建立了L_p-范数形式的模型约束函数,并利用正则化原理将其与重力数据误差函数和已知深度约束函数结合形成V字型密度界面反演的目标函数,推导了目标函数的梯度表达式,并以非线性共轭梯度法为核心给出了反演流程.二维简单模型试算结果表明p=5时该方法能准确地刻画V字型密度界面起伏特征,且亦能准确地应用于二维复杂密度界面和三维界面的反演.最后将反演方法应用于挑战者深渊及邻区的实际资料处理之中,利用研究区海底地形数据和沉积层厚度数据对自由空间重力异常逐层剥离而得到莫霍面引起的重力异常,用本文方法对此重力异常进行反演,结果呈现了板块俯冲作用引起的V字型莫霍面起伏特征.     

地震波形反演的稀疏约束正则化方法

本文考虑地震波形反演问题.为了克服传统的Tikhonov正则化方法过度光滑的弊端,引入了非线性稀疏约束正则化方法,并采用对偶方法求解稀疏约束泛函的极小点.基于二维声波方程波形反演问题进行了数值模拟,针对不同模型对稀疏约束正则化方法进行了测试.结果表明,稀疏约束正则化方法对不连续介质模型的介质边缘具有良好的识别能力.     

基于q-高斯分布和零阶最小熵正则化的三维重力聚焦反演

观测的重力数据通常含有误差,传统的三维重力反演方法假设数据残差服从高斯分布,数据误差的大小和可靠估计对于三维重力反演的精度和可靠性是非常重要的.在模型正则化方面,传统的三维重力反演主要是基于Tikhonov正则化理论来获得稳定解,但是当地下密度异常体存在尖锐边界时,基于L2范数正则化的三维重力反演得到的密度异常体边界比...     

一维声波方程联合反演的分步正则化法

本文用一个纵波信息，对一维声波方程的速度和源函数进行联合反演，并考虑到声波方程的反问题是一个不适应问题，对源函数和波速分别和正则化法分步迭代求解，减少反问题的计算工作量，改善该问题的计算稳定性，为计算实际工程和岩性学问题供了一种方法。文中给出只用一个反问题补充条件同时进行多参数反演的公式，并对相应的数值算例进行分析和比较。     

正则化等效层重力向下延拓方法

近年来,利用时移微重力技术进行储层开发监测受到国内外学者广泛关注.时移微重力观测数据存在信噪比低,信号弱的问题,难以实现储层内物质运移的定量解释.为压制数据噪声,增强有效弱信号,本文研究了利用Tikhonov正则化方法反演等效层（源）,并由等效源实现重力场向下延拓的方法;在此基础上,本文推导了波数域正则化等效源向下延拓算子.针对向下延拓场幅值衰减问题,提出了正则化等效源迭代补偿算法.通过模拟数据实验研究了不同深度正则化等效源滤波算子及向下延拓算子的波数响应;与波数域Tikhonov正则化向下延拓方法相比,正则化等效源向下延拓方法的延拓精度更高、更稳定.最后,将基于迭代补偿的正则化等效源向下延拓技术应用于实测时移微重力数据证实了该方法能够有效增强局部异常,实现时移微重力数据大深度稳定向下延拓.

     

地震偏移反演成像的迭代正则化方法研究

利用伴随算子L*,直接的偏移方法通常导致一个低分辨率或模糊的地震成像.线性化偏移反演方法需求解一个最小二乘问题.但直接的最小二乘方法的数值不稳定,为目视解译带来困难.本文建立约束正则化数学模型,研究了地震偏移反演成像问题的迭代正则化求解方法.首先对最小二乘问题施加正则化约束,接着利用梯度迭代法求解反演成像问题,特别是提出了共轭梯度方法的混合实现技巧.为了表征该方法的可实际利用性,分别对一维,二维和三维地震模型进行了数值模拟.结果表明该正则偏移反演成像方法是有效的,对于实际的地震成像问题有着良好的应用前景.     

水平层状介质中双侧向测井资料的迭代Tikhonov正则化反演

根据非线性反演理论与Morozov偏差原理研究建立从双侧向测井(DLL)资料中同时重构地层原状电阻率、侵入带电阻率、侵入半径、层界面位置以及井眼泥浆电阻率的迭代正则化算法.首先利用Tikhonov正则化反演理论将双侧向测井资料的反演问题转化为含有稳定泛函的非线性目标函数的极小化问题,并利用Gauss-Newton算法确定极小化解.为得到稳定的反演结果并有效实现测井资料的最佳拟合,在迭代过程中将Morozov偏差原理和Cholesky分解技术相结合,建立了一套后验选择正则化因子的方法.最后通过理论模型和大庆油田实际测井资料的处理结果,验证了该算法能够取得更为满意的反演效果.     

基于软约束模式的波阻抗反演

通过在目标函数中增加先验知识约束,我们可以解决波阻抗反演的病态问题及带限问题.基于马尔柯夫随机域（MRF）的邻域系统及势函数构建一种软性约束模式,将井眼约束信息外推到井间,避免过往强制约束的不足.马尔柯夫随机域的一阶邻域偏好于常数值区域,而模型中常梯度值区域和常曲率值区域分别需要用到二阶和三阶邻域,综合使用三种类型的邻域可以改善反演效果.通过理论模型试算和实际资料处理,证明了本文提出的约束模式具有非强制性约束特点,实用性强.     

页岩储层叠前密度稳定反演方法

密度参数是页岩储层评价的重要敏感参数,但通常认为准确获得密度信息难度较大.本文探索提出了一种稳定的基于弹性阻抗的叠前密度反演方法,以提高优质页岩的识别精度.首先,基于不同的参数化弹性阻抗方程开展密度敏感性分析,优选了敏感性最高的Aki-Richards近似弹性阻抗方程作为反演方程.然后,在弹性阻抗方程线性化的基础上,假设反演方程系数矩阵的逆存在,将密度（取对数）表示为多个角度弹性阻抗（取对数）的加权和,通过井旁道弹性阻抗反演结果与测井数据的回归求取加权系数,避免了常规反演方法对大型系数矩阵的求逆计算,从而提高了密度反演的稳定性,同时由于利用测井资料作为匹配目标,密度反演的精度也得到了提高.模型试算和涪陵页岩气田的实际资料应用实践均表明了该方法能够有效地提高密度反演的精度与稳定性.

     

二维波动方程速度的正则化－同伦－测井约束反演

针对二维波动方程反问题,将大范围收敛的同伦方法引入速度参数的反演过程中,并将其与求解不适定问题的Tikhonov正则化有机结合,提出了一种新的、特别适用于非线性的、不适定的、多极值的地震勘探反演问题的反演策略：正则化－同伦方法. 为了充分利用测井资料和地震资料的互补特征,进一步提高反演分辨率并压制噪声,设计了正则化-同伦－测井约束联合反演方法. 大量数值试验结果表明了这两种方法的有效性.     

基于深度学习的2D-3D地下密度异常体反演方法

地下密度异常体反演需要由二维数据反演三维结果,为了高效、高精度地反演地下异常体的位置及密度信息,本文提出2D-3D InvNet深度学习反演方法.方法的编-解码器结构,在编码阶段利用二维卷积网络结构提取地表重力异常及重力梯度异常数据信息,解码阶段利用三维卷积网络结构恢复异常体地下形态,实现了二维和三维网络结构的结合.为了精确地反演异常体密度,提出利用加权均方误差（WtdMSE）作为损失函数,同时为了更好地评估反演结果,引入核密度函数作为评价手段.与均方误差（MSE）相比,利用WtdMSE作为损失函数对异常体密度的反演结果更为准确,异常体所在区域的密度误差减少50%以上.对理论样本的反演结果表明,2D-3D InvNet在准确反演异常体地下位置的同时,也能给出准确的密度信息.应用此方法对西澳大利亚Kauring地区实测数据进行反演,我们成功获得了此区域地下异常体的密度分布.理论与实际应用结果表明,2D-3D InvNet深度学习方法稳定且具有较强的泛化性能,能在无需规定限制条件的情况下快速获得准确的反演结果.

     

约束等效源稳定向下延拓方法研究

在等效源向下延拓方法的理论基础上,研究提出加入观测平面以下空间实测数据作为约束的等效源约束向下延拓法,用于二维位场数据向下延拓.与传统的向下延拓方法和等效源法延拓相比,等效源约束向下延拓法抗噪声性强,能延拓出由于延拓深度过大而衰减至难以探测的高频信号,使延拓结果更接近真实值.理论重力模型数据的延拓结果显示该方法有效,延拓结果优于传统波数域迭代法和等效源法.利用等效源约束延拓法对我国某海域水面实测重力值进行向下延拓并与水下实测值进行对比,结果表明该方法延拓结果最接近实测值.本文方法为我国深水油气近海底重磁高精度探测研究提供一定的技术支撑.

     

基于阈值约束的协克里金法联合反演重力与重力梯度数据

常规协克里金方法反演重力或重力梯度数据具有抗噪性好、加入先验信息容易等优点,其反演的地下密度分布能够识别异常体中心位置,还原异常体基本形态,但反演图像光滑,分辨率低,这是由于常规方法估计的密度协方差矩阵全局发散、平稳.为了通过协克里金方法获得聚焦的密度分布需要改善密度协方差矩阵的性质.首先,本文推导了理论密度协方差公式,其性质表明,当理论模型聚焦分布时,其密度协方差矩阵是非平稳且聚焦分布的.为了打破常规协方差矩阵全局平稳、发散的特征,本文设置密度阈值处理协方差矩阵,通过不断更新协方差矩阵来迭代实现协克里金反演,最终得到相对聚焦的反演结果.用本文方法处理重力与重力梯度数据恢复两种密度模型,均得到了与正演模型匹配的反演结果;再将方法运用于文顿盐丘的实际测量重力与重力梯度数据,反演结果与已知的地质情况匹配较好.