三维黏声最小二乘逆时偏移方法模型研究

针对黏介质成像存在的问题及最小二乘逆时偏移方法的优势,通过广义标准线性固体(GSLS)的三维黏声波动方程,基于三维黏声波动方程的伴随算子及最小二乘逆时偏移框架,实现了三维黏声最小二乘逆时偏移方法.采用极性编码技术大幅度降低黏声最小二乘逆时偏移算法的计算量与内存,使三维黏声最小二乘逆时偏移算法的实用化成为可能.在模型试算部分,首先通过Marmousi模型验证了该方法对传统声波最小二乘逆时偏移方法的优势,最后对三维含Q平层模型及盐丘模型进行试算,证明了三维黏声最小二乘逆时偏移算法的正确性.     

绕射波与一次波联合黏声最小二乘逆时偏移

由于照明不足,小尺度断层和孔洞的成像是一个难题,地下衰减导致地震波振幅损失和相位畸变,在成像过程中忽略这种衰减会造成偏移振幅模糊。基于黏声最小二乘逆时偏移(QLSRTM)能够优化小尺度构造的成像,但这需要大量的迭代和计算成本。为了提高小尺度构造的成像效果,提出了一种充分使用绕射波的面向地质目标的的黏声LSRTM(J-QLSRTM)。在该方法中,构建了新的目标函数和梯度公式,并基于反演理论和伴随理论,推导了一次波和绕射波的Q补偿波场传播算子、Q补偿伴随算子和Q衰减反偏移算子。数值实例证明了J-QLSRTM比传统QLSRTM和声波J-LSRTM更有优势。     

南海天然气水合物储层多孔弹性力学研究

目前大多数对于水合物储层变形的研究都没有考虑毛管力的作用,将储层中的水合物压力和流体压力视为一致。但是在我国的南海海域,水合物储层多为黏土质粉砂,孔隙的孔径较小,因此毛细效应很明显。基于多孔介质力学以及水合物相平衡理论,建立了描述水合物分解过程中的储层变形模型。在该模型中,孔隙内的物质被分为两相,即水合物固体相和包裹不连续气泡的等效流体相。该模型重点研究了水合物固体和等效流体之间的压差,即毛管压力。最后针对我国南海神狐海域的数据具体分析了降压法与注热法两种不同的水合物开采方法。结果显示,对于南海黏土质粉砂储层,毛细效应会对水合物的相平衡条件以及储层变形产生较大的影响。对于给定的压力,水合物的熔点会随着孔径的减小而降低,这使得水合物的相平衡条件在压力-温度图像中表现为一个区域而不是一条曲线,并且如果忽略毛细效应,则会严重低估储层变形。     

基于互相关误差泛函的最小二乘逆时偏移方法

最小二乘偏移可以得到成像质量更高的剖面,有利于更好地进行岩性储层成像和储层参数反演。但是,最小二乘偏移在实用化过程中却很难实现其高精度岩性成像的目标。在分析传统的最小二乘偏移存在的问题后,发展了一种基于相位匹配的最小二乘偏移方法,该方法舍弃了传统的二次型泛函,采用互相关泛函,更多地强调相位信息在最小二乘偏移中的作用,在理论上可以更好地适用于实际资料。通过对互相关误差泛函构建、逆时反偏移合成预测数据、梯度预条件处理、共轭梯度法迭代求解等方面进行研究,探索实现了基于互相关误差泛函的最小二乘逆时偏移方法。简单凹陷模型和复杂salt模型测试表明,本文方法具有较好的适用性和有效性。     

基于稳定逆时传播算子的黏声介质最小二乘逆时偏移

基于GSLS模型黏声介质二阶拟微分方程,采用伪谱法进行数值模拟。针对黏声介质逆时传播过程中产生的高频不稳定问题,提出加入规则化算子对其进行消除的方法,构建了稳定的逆时传播算子。在最小二乘反演的基础上,将黏声介质逆时偏移与最小二乘思路相结合,发展了带有振幅补偿的黏声介质最小二乘逆时偏移( LSRTM)。 Marmousi模型结果表明:相对于常规最小二乘逆时偏移,黏声介质最小二乘逆时偏移校正了地层的黏滞性,得到了更加精确可靠的保幅成像剖面。     

弹性波叠前逆时偏移

利用弹性波一阶速度—应力方程,采用交错网格高阶有限差分法进行震源波场的模拟和接收波场的逆时外推,并结合了能流密度矢量的互相关成像条件,实现了弹性波叠前逆时偏移,且能很好的消除叠前逆时偏移中的假象干扰.笔者通过一个二层模型和一个金属矿模型进行了叠前逆时偏移成像,结果证明所述方法可以有效地消除假象干扰并能对复杂陡倾角界面进行清晰地成像.     

速度误差和地震噪声对最小二乘逆时偏移的影响分析

随着勘探开发目标的日益精细化和复杂化,基于反演理论的最小二乘逆时偏移逐渐成为目前的发展趋势。为了使最小二乘逆时偏移更好地为油气藏的勘探开发服务,分析其不同条件下的成像效果,有助于在实际资料处理过程中指导处理人员做出适当的调整。为此,针对实际勘探过程中出现的偏移速度误差以及地震资料中存在噪声的情况,对比了最小二乘逆时偏移在不同条件下的成像效果。结果表明:在偏移速度存在系统误差的情况下,界面成像深度会出现误差,成像剖面难以反映实际构造形态;平滑后的偏移速度误差对最小二乘逆时偏移影响远小于偏移速度的系统误差,适当平滑可改善成像质量,但是过分平滑会影响成像质量;最小二乘逆时偏移能够压制部分噪声,但随着地震记录中噪声增强,弱反射逐渐被淹没,难以反映地下构造形态。     

横向各向同性介质中的弹性波方程逆时偏移及其成像条件

作者在本文中导出了在横向各向同性介质情况下，二维弹性波逆时传播的高阶差分格式，实现了弹性波在数值空间中的逆时延拓，从snell定量出发，采用插值法求取网格空间中各点的直达波旅行时，以此作为弹性波逆时偏移的成像条件，实现多波多分量资料的逆时偏移，数值试验表明，逆时偏移能够使地表接收到的波场准确归位，提高资料处理的精度。     

基于偏最小二乘与神经网络耦合的储层参数预测

将偏最小二乘回归（PLS）与神经网络（NN）耦合，建立了储层参数预报模型。利用偏最小二乘对影响储层参数的诸多因素进行分析，提取对因变量影响强的成分，从而克服了变量间的多重相关性问题，降低了神经网络的输入维数；同时，利用神经网络建模可以较好地解决非线性的储层参数预测问题。计算实例表明，本耦合模型的拟合和预报精度优于独立使用神经网络模型的精度。     

起伏地表条件下的弹性波叠前逆时偏移

目前,地震勘探条件复杂、地下地质构造多样的“双复杂”地区的地震勘探面临着诸多挑战,特别是地表起伏条件给地下复杂构造精确成像带来了困难。通过研究起伏地表条件下的弹性波叠前逆时偏移方法,采用简洁有效的处理方法解决起伏地表的自由边界条件;利用有效边界的边界存储策略解决弹性波叠前逆时偏移中存储量大的问题。数值实验结果证明该方法能够准确地重建震源波场,SEG加拿大起伏地表逆掩断层模型进一步验证了提出方法的有效性。     

基于解耦传播的波场分解方法在VTI介质弹性波逆时偏移中的应用

由于弹性波逆时偏移更符合实际情况,而且转换波的成像结果有更高的分辨率,因此弹性波逆时偏移越来越受到重视。然而,弹性波逆时偏移需要多波多分量数据,为了减少成像结果中的串扰假象,在逆时偏移过程中进行P波和S波分解就变得非常必要。结合基于向量的激发振幅成像条件,我们把基于解耦传播的波场分解方法应用到弹性波逆时偏移中,并对比了其在各向同性介质和横向各向同性（VTI）介质中的应用效果。结果说明,该方法可以在各向同性介质中完全分解P波和S波,并保留波场的向量信息。虽然在VTI介质中有较小的分解残余,但是该分解残余不会在逆时偏移结果中产生明显的串扰;因此,这种波场分解方法可以应用于各向同性介质和VTI介质弹性波逆时偏移。该方法是在时间空间域实现的,可以在波场传播过程中直接对P波和S波进行分离,应用方便,计算效率高。与利用Helmholtz分解的弹性波逆时偏移相比,该方法避免了在PS波成像结果中的极性反转问题。在复杂Hess VTI模型的逆时偏移结果中,高速岩体和断层的成像清晰,甚至是两个低速薄夹层也能较好成像;这说明该方法对复杂介质具有较好的适应性。PS波成像结果中的各向异性体成像清晰,说明各向异性介质弹性波逆时偏移可以对传统逆时偏移不能很好成像的构造进行成像。     

基于P＋S波震源的弹性波叠前逆时偏移方法

基于弹性波动方程,推导了二维非均匀介质条件下纵横波解耦且无转换波的一阶速度一应力弹性波方程,使得弹性波在延拓过程中实现矢量偏移保证偏移质量,且不产生其他类型的转换干扰波。从激发时间成像条件以及反射系数理论出发,并与纵横波解耦且无转换波的一阶速度一应力弹性波方程结合,给出了适用于P＋s波震源模拟数据的成像条件以及偏移方法。文中给出了两个典型模型的理论记录的偏移结果,试验表明,笔者提出的方法能够对P＋s波震源的产生的波场准确归位,提高资料处理精度。     

基于双相介质弹性波波场分离的逆时偏移方法

随着地震勘探的发展,以单相介质模型为基础的声波逆时偏移和弹性波逆时偏移很难比较精确地描述地下介质和满足复杂的工程要求,尤其是在处理地下含流体多孔介质时,所以需要引入双相介质模型。首先基于Biot理论双相介质一阶速度-应力弹性波方程,利用散度场和旋度场的波场分离理论将固相和流相中的纵横波场解耦分离,并由波场快照图得出固相与流相中的纵波波场与横波波场均完成分离,而快纵波与慢纵波相互耦合均存在于纵波波场中。采用高阶交错网格有限差分方法,结合PML吸收边界条件,推导了具有波场分离的双相介质弹性波正传与反传的偏移延拓算子,结合互相关成像条件,可以得到双相介质弹性波固相与流相中的各分量偏移成像剖面。最后对模型使用该双相介质弹性波波场分离的逆时偏移方法进行偏移成像试算并取得成功,效果较常规双相介质弹性波逆时偏移成像效果有了明显提升。     

晶体点阵参数最小二乘法计算程序

作者以最小二乘法理法导出了晶体点阵参数计算的公式,并顾及了X射线粉晶衍射图测定所引起的系统误差修正,从而得到点阵参数的最佳估值。以导出的公式为数学模型编制了PC-1500微机的扩展BASIC语言程序,可适用于各晶系点阵参数的计算,使用方便。     

木里地区天然气水合物储层参数计算

测井作为一种重要的地球物理勘探手段,在天然气水合物资源定量评价中发挥着重要的作用,随着木里地区天然气水合物勘探井的日益增多,需建立一套适合本地区的泥质含量、含水饱和度和孔隙度等储层参数的计算方法,这对储层优选、资源量的评估等具有重要意义.笔者基于国内外文献常用的天然气水合物测井评价模型,选用经过环境校正和标准化处理后的测井数据,定量计算天然气水合物的孔隙度与饱和度等参数,结合岩心测试和岩电资料,修正并完善储层参数的计算模型,分析各种测井评价方法在水合物储层计算当中的优劣.研究结果表明:泥质含量选用自然伽马值求取为宜;孔隙度选用受水合物分解所产生气体影响较小的密度值求取;饱和度的求取需根据水合物含量的不同,分别选用印度尼西亚公式、修正的wood方程或修正的阿尔奇公式.     

最小二乘曲线拟合在原子荧光光谱分析中应用

曲线拟合是目前数据处理的基本方法之一,通过在微机上用Visual Basic语言编制程序,将实验数据使用最小二乘曲线拟合,构造一个最佳通用多项式拟合模型。再利用得到的函数关系实现对原子荧光光谱仪的实验数据进行多次曲线拟合与仿真,以扩大仪器的测量范围。     

二维视电阻率断面的快速最小二乘反演

快速最小二乘反演是以平滑限定的最小二乘方法为基础,是对二维视电阻率断面进行反演的一种方法。反演过程不需要提供初始模型,在首次迭代时使用一均匀介质地下模型作为初始模型,该模型的视电阻率偏导数值可以用解析法得到。在后面的迭代中,使用了拟牛顿法去修改每一次迭代的偏导数矩阵,避免了偏导数矩阵的直接计算,从而减少了计算时间和存储空间。同时运用牛顿矩阵校正技术解最小二乘方程组也减少了大量的计算时间。总之,该方     

分散型天然气水合物对储层的物性影响特征分析

天然气水合物以不同的赋存形态存在于温压条件适合的海底沉积物及永久冻土带中。目前对于天然气水合物的勘探以地震方法为主,而随着天然气水合物钻探工作的不断开展,地球物理测井成为必不可少的工作,测井通过不同的方法采集储集层的地球物理参数,由于天然气水合物的存在对储层物性的特殊影响,在利用测井方法评价储集层时需要综合考虑这些影响因素,才能获得准确的储层参数。本文详细阐述了由于天然气水合物的存在对地层的声学、电学及密度特性的影响及其在利用这些参数评价储层时需注意的问题。     

频率域最小二乘反演

本研究目的在于利用频率域的全波场有限元最小二乘反演求取地下更多的岩性信息如各向异性弹性参数、粘性参数等物性参数,对石油勘探中油气藏的直接检测和油藏描述等岩性解释很有帮助。由于直接求取所有的各向异性及粘性参数是高度非线性的、非常困难的,所以先选择了比较简单一些的弹性波速度信息：纵波速度与横波速度进行求取,这比常规声波方程资料处理能够多得到横波的信息,对岩性解释很有帮助。在频率域弹性波有限元正演方程的基础上,依据使目标函数(也就是观测数据和正演数据残差的二次范数)最小的准则,导出了频率域弹性波有限元最小二乘反演算法。其中用矩阵压缩存储与LU分解技术存储与求解正演方程中的大型稀疏复系数矩阵,用Green函数更快速的求取Jacobian矩阵,用共轭梯度算法求解大型稀疏矩阵的最小二乘解。