三维重力、重力梯度和大地电磁数据的平滑聚焦结构约束联合反演

多种地球物理数据联合反演解释能够显著提高勘探分辨率,近年来,基于交叉梯度结构约束的联合反演得到了广泛研究.该方法有效克服了单一方法反演易受多解性影响和不同方法反演结果结构不一致等问题,已成为目前最流行的联合反演技术之一.传统的重力和大地电磁联合反演常因重力纵向分辨能力弱,导致密度模型可靠性较低,进而使联合反演结果有失准...     

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

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

电磁-重力联合建模约束反演圈定GS凹陷深层目标（英文）

GS凹陷深层上古生界煤系地层具有生烃能力,圈定出该套地层分布对凹陷油气评价意义重大。在缺少地震资料的情况下,依靠少量的钻井、地质资料,尝试利用电磁-重力联合建模约束反演圈定深层目标。首先根据三维大地电磁反演结果和已知钻井地质资料建立初始地质模型;再结合地层密度和重力异常特征,利用重力交互反演修正优化初始地质模型;然后对该模型进行大地电磁资料的电性约束反演,并进行拟合差评估,获得最优地质模型;最后将该地质解释方案与位于凹陷中部的一条地震剖面进行综合分析,确定地震剖面深层反射就是上古生界。该套烃源层主要分布于凹陷深部位,最大厚度达800m。实测资料处理结果表明,电磁-重力联合建模约束反演可以逐步降低单一地球物理方法的多解性,提高对深部地层的识别能力。非地震与地震结果相互印证,展示了电磁-重力联合建模约束反演技术在研究类似盆地深层目标中具有较好的应用潜力。     

大地电磁、重力、磁法和地震初至波走时的交叉梯度二维联合反演研究

为了准确的探测和描绘地下复杂的地质结构,同时克服地球物理单一方法反演的多解性和单一参数反演模型的不一致性等问题,近年来基于交叉梯度联合反演的综合地球物理解释已经得到了广泛的关注和应用.本文首先研究了两种地球物理方法的交叉梯度联合反演算法,在此基础上,推导并实现了多种地球物理方法(大地电磁,重力,磁法,地震初至波走时)的多交叉梯度约束的二维联合反演算法;其次,我们设计了结构不一致模型和复杂模型,针对多物性联合反演算法的准确性和有效性进行了模拟试算,并对复杂模型的单独反演结果和联合反演结果进行了交叉梯度值和物性交会图的对比;最后,本文将成熟的卫星资料多光谱综合分析技术应用到联合反演中,将多物性参数反演模型结果图通过RGB(红-绿-蓝)模式进行合成,得到融合的RGB合成图.结果表明:通过对结构不一致模型和复杂模型的联合反演结果和单独反演结果的对比分析,可以得出联合反演得到的结果更接近真实模型,并从得到的交叉梯度值进一步证明了联合反演模型相似度高,也从物性交会图中得到联合反演的物性相关性更好的结论,反向证明了算法的正确性.最终从得到的RGB合成图像,我们可以更直观的分析反演结果,更有利于准确划分地下模型结构.     

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

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

含噪声数据反演的概率描述

根据贝叶斯理论给出了对含噪声地球物理数据处理的具体流程和方法,主要包括似然函数估计和后验概率计算.我们将数据向量的概念扩展为数据向量的集合,通过引入数据空间内的信赖度,把数据噪声转移到模型空间的概率密度函数上,即获得了反映数据本身的不确定性的似然函数.该方法由于避免了处理阶段数据空间内的人工干预,因而可以保证模型空间中的概率密度单纯反映数据噪声,具有信息保真度高、保留可行解的优点.为了得到加入先验信息的后验分布,本文提出了使用加权矩阵的概率分析法,该方法在模型空间直接引入地质信息,对噪声引起的反演多解性有很强的约束效果.整个处理流程均以大地电磁反演为例进行了展示.     

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

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

大地电磁的人工鱼群最优化约束反演

大地电磁的反演问题是非线性,如果采用线性反演方法容易陷入局部极小,使得反演结果非唯一性严重.本文将人工鱼群算法引入到地球物理反演之中,提出了非线性的大地电磁人工鱼群最优化反演.该方法不需要进行偏导数的求取,可以对反演的范围进行约束,以减小反演结果的非唯一性.同时我们对搜索步长进行了改进,给出适用于大地电磁反演的人工鱼群参数.大量的理论数据试算表明,人工鱼群反演算法能够较好地寻找到全局最优解.实测数据的处理结果表明,该方法可以用来处理实际资料,并且能够取得很好的应用效果.     

基于交叉梯度结构约束的大地电磁主轴各向异性并行三维反演

各向异性介质对大地电磁观测数据的影响往往不可忽略,因此需要提高大地电磁各向异性三维反演的可靠性和有效性.为了满足大地电磁各向异性三维反演的需求,本文研究了一种基于交叉梯度结构约束的大地电磁主轴各向异性并行三维反演算法.根据大地电磁平面波理论假设,正演方程采用背景场与二次场分离的计算方式,二次场利用交错网格有限差分法求解.由于各向异性反演的多解性,本文将各向异性介质简化为主轴各向异性,并在此基础上进一步采用有限内存拟牛顿LBFGS法实现三维各向异性反演.为了提高各向异性反演的分辨率,反演目标函数中引入交叉梯度项,利用先验的结构信息,对三个方向的电阻率参数进行结构约束,最终的反演进一步利用MPI(Message Passing Interface,消息传递接口)技术实现分频并行计算,测试结果显示并行接近线性加速比.     

基于改进差分进化算法的大地电磁测深和重力宽范围物性约束联合反演

宽范围物性约束技术容易实现、具有一定容错性,目前已在大地电磁测深（MT）和地震、MT和重力联合反演中实现,但该技术是结合模拟退火算法实现的.差分进化算法（DE）是一种全局优化算法,但该算法在地球物理联合反演领域应用较少.基于此,本文以双种群设置方案为框架改进了DE算法,并提出了基于改进DE算法的宽范围物性约束技术.MT和重力联合反演的模型试验表明：与传统的DE算法相比,改进的DE算法收敛速度更快,寻优能力更强;基于改进DE算法的宽范围物性约束技术可以促进不同岩石物性参数在一定"范围"内实现耦合,既可以利用岩石物性关联的导向作用,又可以发挥优化算法的寻优能力,进而降低地球物理联合反演对先验信息的要求;此外,该技术的实现也验证了宽范围物性约束思想在联合反演领域中的适用性,具有进一步推广至其他优化算法中的潜质.     

部分区域约束下的交叉梯度多重地球物理数据联合反演

交叉梯度联合反演方法通过对多种地球物理模型实现结构耦合,在岩石物性关系不确定的情况下,既能提高反演结果的可靠性,又能减少反演的多解性,还能减少不同方法解释结果之间的矛盾.当不同的模型观测数据覆盖范围不一致时,交叉梯度联合反演通常需要取出重叠区域数据进行联合反演,并且建模时还要扩展一些模型范围.本文首先提出并实现了部分区域约束下的交叉梯度多重地球物理数据联合反演算法;接着进行了算法的模型试算;最后,我们将该反演算法用于本溪—集安深部地质调查重磁电综合地质地球物理解释中.结果表明:该算法不但能在重叠区域内很好地恢复结构相似的模型,而且在非重叠区域与重叠区域的边界处仍然可以得到平滑变化的模型;在本溪—集安10号剖面所获得的结构上相似的电阻率、密度及磁化率模型较好地反映了该区的深部地质结构,对于确定深部地质体的性质提供了有力的证据.     

一种高效的基于交叉梯度结构约束的三维地震走时与直流电阻率联合反演策略

基于交叉梯度结构约束的不同数据的联合反演可以提高地球物理成像的可靠度,但是由于不同观测数据对地下介质的灵敏度不同以及不同反演系统灵敏度矩阵元素的数值可能存在较大的差异,导致构建和求解联合反演系统存在很大的挑战.针对以上问题,本文提出一种新的基于单独反演模型更新量与交叉梯度结构约束相结合的联合反演策略.该策略利用单独反演系统分别确定出各个模型的更新量,然后利用它们约束交叉梯度系统的反演,得出新的模型更新量.通过这样的方式,有效实现了数据拟合与结构约束的平衡,实现了基于交叉梯度联合反演的目标.新的联合反演策略不需要对原来单独反演程序修改即可实现联合反演,减小了联合反演实现的难度,极大地提高了联合反演的易实现性,而且避免了联合反演矩阵存储及结构过于复杂难以求解的问题.基于新的策略,本文首次实现了基于交叉梯度结构约束的三维地震走时与直流电阻率联合反演.合成模型测试表明,与单一成像相比,联合成像减少了地震走时反演中出现的干扰异常并提高了电阻率反演的分辨率.     

基于交叉梯度交替结构约束的二维地震走时与全通道直流电阻率联合反演

在利用不同的地球物理勘探方法对地下复杂介质成像时,因观测系统的非完备性及数据本身对某些岩石物性的不敏感性,单独成像的结果存在较大的不确定性和不一致性.对于地震体波走时成像与直流电阻率成像,均面临着成像阴影区问题.对于地震走时成像,地震射线对低速区域覆盖较差形成阴影区,造成低速区域分辨率降低.对于电阻率成像,电场线在高阻区域分布较少,造成高阻区域分辨率较低.为了提高地下介质成像的精度,Gallado和Meju(2003)提出了基于交叉梯度结构约束的联合地球物理成像方法.在要求不同的物性模型拟合各自对应的数据同时,模型之间的结构要求一致,即交叉梯度趋于零.为了更有效地实现基于交叉梯度的结构约束,我们提出了一种新的交替结构约束的联合反演流程,即交替反演不同的数据而且在反演一种数据时要求对应的模型与另一个模型结构一致.新的算法能够更容易地把单独的反演系统耦合在一起,而且也更容易建立结构约束和数据拟合之间的平衡.基于新的联合反演流程,我们测试了基于交叉梯度结构约束的二维跨孔地震走时和直流电阻率联合成像.合成数据测试表明,我们提出的交替结构约束流程能够很好地实现基于交叉梯度结构约束的联合成像.与单独成像结果相比,地震走时和全通道电阻率联合成像更可靠地确定了速度和电阻率异常.     

航空电磁拟三维模型空间约束反演

为了克服时间域航空电磁数据单点反演结果中常见的电阻率或层厚度横向突变造成数据难以解释的问题,通过引入双向约束实现航空电磁拟三维空间约束反演.除考虑沿测线方向相邻测点之间的横向约束外,同时还考虑了垂直测线方向测点在空间上的相互约束.为此,首先设计拟三维模型中固定层厚和可变层厚两种空间约束反演方案,然后通过在目标函数中引入沿测线和垂直测线方向上的模型参数约束矩阵,并使用L-BFGS算法使目标函数最小化,获得最优拟三维模型空间反演解.基于理论模型和实测数据反演,对单点反演与两种空间约束反演方案的有效性进行比较,证明本文空间约束反演算法对于噪声的压制效果好,反演的界面连续光滑,同时内存需求和反演时间少,是一种快速有效的反演策略.     

基于宽范围岩石物性约束的大地电磁和地震联合反演

岩石物性约束是地球物理联合反演重要的实现方式.在以往联合反演研究中,这种约束通常以点对点的映射方式出现,要求经验性的岩石物性关联是相对精确的,无形中提高了岩石物性约束应用于联合反演的条件.基于此,本文结合全局寻优的模拟退火算法,提出了宽范围岩石物性约束技术.该技术将岩石物性约束与模拟退火模型扰动、模型接受准则相结合,实现了空间对空间的映射,将电阻率和速度耦合在一起.文章以MT和地震正则化同步联合反演为例,检验了不同的先验信息条件下的宽范围岩石物性约束技术的适用性.模型试验表明:新技术容易实现,容错性高,可降低"不精确"岩石物性关联信息引入带来的风险,提高岩石物性约束应用的灵活性.     

混合范数下的最优化反演方法

在求解地球物理反问题时，通常根据最小二乘准则构造目标函数进行反演，并在实践中得到了广泛的应用.为进一步增强反演的稳健性及减少多解性，不损失反演结果的分辨率，本文提出了混合范数下的最优化反演方法，它根据数据和模型可能服从不同的概率分布，对数据空间和模型空间采用不同的范数来构造目标函数.在给出目标函数的基础上，导出了混合范数下的线性反演方程.由于该线性反演方程的复杂性，我们采用混合范数下迭代再加权共轭梯度法进行求解.最后，通过对模拟的电阻率数据进行反演，验证了本文计算方法是可行的.