利用共轭梯度方法的激发极化三维快速反演

利用共轭梯度方法实现了激发极化(IP)三维快速反演。首先,利用共轭梯度方法反演电位数据,得到地下的三维电阻率模型,(由于避免了直接求偏导数矩阵,因此反演计算速度很快。)然后,以此电阻率模型为地下电导率分布,再反演IP数据得到三维极化率分布理论模型。试算结果表明其效果较好。     

激发极化数据的最小二乘二维反演方法

给出了一种电阻雍激发极化法数据的最小二乘二维反演方法,该方法与以往这类方法的主要差别在于网格单元中电导率和极化率参数呈线性变化而不是均匀值,从而使电经和极化率的反演结果更准确,更容易用等值线成图,方法在二维有限元正演计算中,采用三角单元,使实测数据反演前不需要进行地形改正;在目标函数中加入最简单模型以及背景场等先验信息,既民反演问题的多解性,又使反结果更接近实际情况,还采用了前人用电位函数与模型参     

井—地三维激发极化数据快速反演

为实现井—地激发极化数据的快速反演,分析三维有限元正演四面体网格剖分形成的系数矩阵的元素规律,应用MSR非零元素压缩存储和SSOR-PCG方程求解技术,提高了正演执行效率。在正演基础上,利用Jacobian矩阵与正演方程的简单关系,采用共轭梯度法求解最小二乘目标函数的模型修改量,避开Jacobian矩阵的直接计算和存储,实现了三维井-地激电数据的分钟级快速反演。模型正演合成数据的反演结果表明,能基本再现已知的模型。     

电阻率/激发极化率数据的二维反演程序

本程序在目前日本流行的电阻率／激发极化法二维反演程序［４］基础上,将电导率和极化率分块均匀变为分块连续变化,并在目标函数中加入了更多的先验信息,使反演结果的唯一性和真实性大大增加。在有限元正演模拟中,又采用三角单元剖分,使实测数据在反演以前不需要进行地形改正。本文介绍了该方法的基本原理和程序的使用方法,并给出了几个理论算例和野外实例     

三维重力资料反演

我们介绍两种三维空间密度分布重力资料反演方法。在第一种方法中，我们根据重、磁间的泊松关系将重力资料转换成伪磁力资料并用三维磁法反演算法进行反演；在第二种方法中，我们研制了直接反演方法，由重力资料直接反演出密度分布的极小结构模型。在以上两种方法，三维地质体由大量长方体元密度分布构成，每一体元中密度值为常数。     

三维地电断面激发极化法有限元数值模拟

用有限单元法求解三维地电断面激发极化法正演模拟算法.首先给出了三维构造中点源电场异常电位的边值问题与变分问题, 简化无穷远边界上的边界条件以提高计算速度及计算精度.以此为基础计算视电阻率对模型电阻率的偏导数矩阵, 并进行三维地电断面激发极化法正演模拟, 与等效电阻率法相比节约了计算时间.对几例较典型的地电模型进行计算, 结果表明本方法是正确可行的.      

三维井地电阻率有限元数值模拟及反演

油气资源和矿产资源勘探中,常需要较大的勘探深度和高的深部分辨能力,充分利用现有钻孔的三维井地电阻率成像技术能够适应于这种要求。通过设计一个高阻水平板状体模型,利用Ansys有限元软件对模型进行三维井地电阻率法有限元正演计算,并将模拟的响应视电阻率作为反演的输入,结合平滑约束最小二乘反演,分析了异常特征和分布规律,其结论和认识具有一定实用价值。     

岩石激发极化衰减谱的多指数反演

激发极化衰减谱比极化率具有更为丰富的储层信息,与地层的孔隙结构、渗透率、地层水电阻率以及泥质含量密切相关。衰减谱是一系列单指数衰减的线性叠加,因此,可以对衰减谱进行多指数反演,得到弛豫时间谱,从而可以更好地利用激发极化衰减谱来进行储层评价。一般的非负最小二乘反演的结果振荡剧烈,与实际情况不相符。这里采用阻尼最小二乘方法,结合奇异值分解技术对泥质砂岩的激发极化衰减谱进行了多指数反演,得到了光滑连续的弛豫时间谱,弛豫时间分布点数为32点~64点较合适,算法和软件对实际的岩心具有较好的效果。     

基于模拟测井资料的波阻抗反演

波阻抗反演技术把具有高纵向分辨率的已知测井资料与连续观测的地震资料联系起来,实行优势互补,是煤田岩性地震勘探的重要手段。今天,模拟测井资料占主导地位是不容忽视的现实,因此,关键问题是如何利用模拟测井资料做好波阻抗反演。这里针对上述问题进行了有益的尝试,提出了解决这一问题的方法。从对实际资料进行处理后显示,该方法提高了地震资料的分辨率,能够推断地层岩性在平面上的变化及顶、底板岩性情况。     

交叉梯度约束的三维极化率/电阻率联合反演

为减少反演的多解性,多数情况下期望电阻率和极化率异常具有相同结构。笔者利用交叉梯度算法实现了联合反演。在计算交叉梯度项过程中,提出一种新的离散化网格边界处理方法,提高了交叉梯度项计算的速度和精度。试算结果表明三维电阻率/极化率联合反演能有效地提高同结构异常体的反演效果。     

使用共轭梯度的三维电阻率正演模拟与反演

研制了利用共轭梯度松驰法的三维直流电阻率正演模拟与反演算法。在正演网络的模拟计算中，将一种用于预处理的不完全Ｃｈｏｌｅｓｋｙ分解法与稀疏矩阵程序结合起来，便可得到一种快速有效的算法。侧面和底部的边界条件用阻抗条件来标度，而阻抗条件考虑了由于电源流的一些配置而在边界上产生的电流。     

三维直流电阻率的非线性反演成像

研制了直流电经数据解释的全三维反演成像方法。由于反演算子是严重不适定的，所以我们采用正则化算法来解释这些反演问题，为了得到高分辨率成像和一种精确而的反演方法，利用非线性正则化反问题的准凸面特性，并选取最佳正则化参数，反演算法能单独解释地面数据，单孔数据，跨孔数据，或者这些数据的任意组合，得到了综合数据组的高分辨率的成像。     

泥质砂岩激发极化弛豫时间谱的正则化反演

激发极化衰减谱是一系列单指数衰减的线性叠加,因此可以对衰减谱进行多指数反演,得到的弛豫时间谱可用于储层孔隙结构和渗透性评价,是一种非常有前途的技术。由于实际测井的信噪比相对较低,常规的奇异值分解方法难以获得合理的弛豫时间谱,为此,采用正则化约束,结合奇异值分解,对合成数据和实际测量的衰减谱进行反演,得到连续的弛豫时间谱;同时,讨论了弛豫时间谱的时间常数分布、正则化因子以及信噪比对反演结果的影响。研究结果表明:该方法适用于信噪比较低的实际测量结果的反演,弛豫时间常数的分布点数为32~64较为合适;随着正则化因子的增加,反演结果逐渐变得平滑,存在一个最优的正则化因子,用于获得最为合理的弛豫时间谱;在双对数图中,最优正则化因子随着信噪比的增加而线性降低,因此通过测量结果的信噪比,可以求取最优正则化因子,进而将该技术应用于实际。     

用于三维层析成像反演的适应几何准则

确定深度波速是油气勘探的一个战略域，为此目的，我们发展三维地震资料层析成像反演程序，它用了ＧＯＣＡＤ三模拟器作反射布置，以及样条张量积作为层速度。反资料为叠加时间及叠加速度。此方法应用于实际数据，证明是高效的