基于地震剖面瞬时峰值能量频率的衰减分析与气藏预测方法

由于流体本身的粘滞性和摩擦性,使得地震波在含油气地层中传播时衰减形成低频阴影.通过对叠后地震资料做时频分解,分析地层衰减特征,是当前除了AVO技术之外,比较流行的直接识别流体和气藏的技术.本文基于波动方程正演,模拟了地震波对含油气地层的响应.然后采用连续小波变换分析不同尺度剖面的瞬时能量特征,提取分频剖面上每一个点的峰值能量对应的频率,得到瞬时峰值能量频率剖面.通过分析瞬时峰值能量频率剖面,分析了地震波在地下介质中传播时的衰减特征.最后应用此方法对实际地震剖面进行了分析,验证了其有效性.     

背景为各向异性的含裂缝岩石频散和衰减计算方法研究

裂缝广泛分布于各类储层岩石中,并且会显著提高储层的渗流能力.因此,裂缝的评价和表征对于提高油气产能具有重要意义.由于裂缝与背景介质之间的波致流会显著影响地震波的频散和衰减特性,所以地震勘探是评价裂缝性储层的有效手段.裂缝地震定量表征的前提是要基于含裂缝岩石中波致流对频散和衰减的影响建立含裂缝岩石物理特性与地震性质的关系.然而,目前相关的理论研究大部分基于各向同性背景这一假设,难以有效应用于常见的各向异性储层.本文针对背景为各向异性的含裂缝岩石提出了频散和衰减的计算方法.该方法首先将含裂缝岩石中的各向异性背景介质等效为层状背景介质;然后,通过分析不同频率下层状含裂缝岩石中的流体压力分布,理论计算了两个特定的中间频率并求解得到两个中间频率下的弹性参数;进一步,以计算得到的两个特定中间频率以及高低频极限下的弹性参数为基础,应用数值方法求解得到弛豫函数中的未知参数,最终实现了背景为各向异性含裂缝岩石中频散和衰减的理论模拟.通过将理论预测结果与实验测量和数值模拟结果进行对比,验证了该方法在背景为各向异性含不同分布裂缝岩石中的有效性.本文提出的方法考虑了常见的各向异性背景对含裂缝岩石频散和衰减的影...     

气枪震源水中能量辐射特性研究

基于云南宾川气枪主动源激发实验数据,对气枪震源水中能量辐射特性进行了分析研究,得到了气枪震源水中激发脉动气泡周期和最大半径的经验计算公式。结果表明:(1)气枪震源的能量辐射效率在垂直向下方向要高于水平方向;(2)垂直向下的气泡能辐射效率随着供气压力的增加而增大;(3)供气压力为15MPa时,气枪震源垂直向下的能量辐射效率为9.4%;(4)供气压力越大,气枪震源压力波叠加增强的持续距离越远。     

近地表地震地质复杂性的一种定量分析方法

复杂近地表散射衰减所有的地面观测波场,形成半随机半相干的近地表强散射噪音背景,弥漫整个炮集,淹没深层反射信号,是导致地震资料极低信噪比的主要原因.如何研究和评价近地表散射强度一直是石油勘探未解决的问题,这与起伏地表的粗糙度、近地表速度横向变化和结构倾角分布密切相关.基于前期复杂近地表边界元法波动方程数值模拟研究,本文提出一种复杂近地表散射振幅矩阵方法来分析近地表散射强度.首先对复杂近地表结构进行边界元配置方法离散,根据边界积分方程生成矩阵方程.我们不求解该矩阵方程(涉及海量计算),只是利用矩阵分析技术来解析矩阵方程中的散射振幅系数矩阵,研究复杂近地表结构对不同频率波场的散射强度.该方法利用边界元对近地表结构几何特征的精确表征,研究起伏地表和非规则地质分界面对地震波传播的影响,由基本解及其在边界上的法向导数经过高斯数值积分计算得到的散射振幅系数矩阵,不仅描述了任意两点之间的相互影响,同时还刻画了边界形状特征的影响,为评价不同地质结构的散射强度提供了可能性.作为初步评价手段,我们采用矩阵元素总和与矩阵维数之比作为表征散射振幅系数矩阵散射特征的标量复杂系数,通过理论和实际模型测试,形成了一套行之有效、计算快速的近地表复杂性分析方法.     

基于修正伪Helmholtz分解算子的横向各向同性介质纵横波波场分离

纵横波波场分离是弹性波逆时偏移成像中的关键步骤,可以有效消除纵横波串扰和成像畸变,提高成像精度。在各向异性介质中,常用的波场分离法有非平稳滤波器或者低秩近似的方法,但由于使用多次傅里叶变换,导致波场分离计算成本较高。借鉴基于Helmhotz算子构建各向同性介质纵横波解耦方程的思路,提出消除振幅畸变的修正伪Helmholtz分解算子,通过待定系数方法求解横向各向同性介质中P波和S波2种模式的表达式,并将其变换为一阶解耦的伪弹性波方程,实现时空域纵横波波场解耦。通过简单模型的波场分离测试,得到分离后的P波与S波波场,验证了该方法的波场分离的有效性。此外,将解耦得到的矢量P波与矢量S波应用到弹性逆时偏移之中,利用矢量点积互相关成像条件得到清晰的弹性波逆时偏移成像结果,说明本方法在复杂介质中有着较好的适用性,同时可以有效地应用到VTI介质弹性逆时偏移之中。     

基于二维图像的数字岩心电导率计算方法研究

电导率是表征岩石电学性质的重要物理参数,在地质资源勘查和测井解释等领域发挥着巨大作用.快速、准确地确定岩石电导率具有重要的理论和实践意义.作为近年来发展的一种岩石物理数值模拟工具,数字岩心技术在定量计算电导率等物性参数方面应用广泛.三维微观结构的准确获取是数字岩心技术计算岩石电导率的关键,但传统获取岩石三维微观结构的方法较为复杂费时.为了方便快速地通过数字岩心技术计算岩石的电导率,本文研究了岩石二维与三维数字岩心的电导率联系.我们基于微米级X射线CT扫描得到的三个砂岩样品的微观结构信息建立了三维数字岩心,并通过有限元法计算的三维数字岩心电导率与实验数据的对比验证数值计算方法的有效性.随后我们数字地扩展了岩石的孔隙,产生了较大孔隙度的三维数字岩心样本,在此基础上,计算了三维数字岩心和相应二维数字岩心的电导率,并通过Archie公式分别拟合了电导率与孔隙度之间的关系,得到了相应的胶结系数.结果表明,三维数字岩心的胶结系数小于二维数字岩心的胶结系数,且二者的比值与岩石实测孔隙度呈线性负相关关系.以该联系为纽带,通过二维图像快速计算得到的电导率与孔隙度关系,确定了三维数字岩心的电导率与孔隙度关系,并进一步通过三维数字岩心的孔隙度计算其电导率.该方法计算得到的人工砂岩样品的电导率与其三维数字岩心电导率相关系数高于96%,验证了基于二维图像的数字岩心电导率计算方法的有效性.本文的研究结果为快速、准确地计算岩石电导率提供了新的思路,在油气勘探开发中有广阔的应用前景.     

热弹性介质中波传播特征

深层-超深层油气地震勘探涉及高温介质地震波传播问题,热弹介质参数对地震波传播有重要影响.含弛豫时间修正项的Lord-Shulman双曲型耦合热弹波动方程从理论上预测了热弹性介质中存在快纵波、慢纵波(一种准静态慢纵波,简称热波)和横波的传播,两个纵波为热耗散衰减波而横波不受介质热特性的影响.本文结合平面波频散分析和格林函数法数值模拟,详细研究两个热耗散衰减波的频散和衰减特征,着重分析热导率、热膨胀系数及比热的变化对波速和衰减的影响.研究表明热导率作为主要参数决定了波速与衰减的临界变化,热膨胀系数对波速和衰减的幅度有明显影响,比热则兼顾了前两个热弹系数的影响特征.最后,利用热弹性动力学频率域的二阶格林函数进行波场快照数值模拟,展示热弹性介质中纵波、横波和热波的传播行为.     

基于耦合反射/透射系数单程波传播算子的地震波模拟研究

单程波近似实际上是一种多次前向散射和单次后向散射近似.利用单程波近似来描述波传播可以极大地节省地震数值模拟的计算时间和内存,实现地震波长距离传播模拟和三维地震模拟快速计算.本文基于单程波近似和波动积分方程的分离变量逼近,从广义Lippmann-Schwinger波动积分方程推导出耦合反射/透射系数的单程波传播算子.该算子由两部分构成：分离变量Fourier单程波传播算子和薄板间的反射/透射系数表达.前者将常规的Fourier分裂步单程波传播算子(SSF)推广适应横向强速度变化介质和大角度传播波场.后者是利用垂直波数来表示反射/透射系数,自然耦合到波场传播的计算过程中,其为地质界面倾角的隐式表达,精确描述振幅随入射角的变化,能适应任意复杂的模型.通过两个数值算例和一个实际地质模型的计算,本文将该方法和边界元法进行了比较,结果表明：在算例给出的介质横向速度变化情况下,本文提出的方法在相位和振幅方面与全波数值方法基本吻合.     

复杂地表低信噪比地震数据处理研究

焉耆盆地地质条件复杂,侏罗纪末燕山运动和第三纪喜山运动的挤压推覆作用在焉耆盆地南缘霍拉山前形成一系列推覆构造.由于地形起伏大,近地表结构复杂,采集的地震资料信噪比较低;以推覆构造为主体的地下地质复杂性给地震成像造成极大的困难.针对该复杂地区低信噪比地震资料,在偏移前资料预处理中,通过采取一系列旨在提高信噪比、重建反射信号的技术,为后续的复杂构造成像处理奠定基础.针对断裂交错的复杂推覆构造成像,采用增强高波数波场成像的半解析Fourier地震偏移方法,能很好地成像强速度对比下的陡倾角大断裂地质构造,有效突出推覆体构造的空间形态,改善推覆体下的成像效果.     

海上时移地震中多道匹配的观测系统重复性研究

观测系统的重复性是时移地震的关键问题.而对于全面评价,多道重复性度量又是很有意义的.基于归一化均方根差异计算方法,地震多道数据的重复性一般为单道数据重复性的加权均方根.针对海上拖缆监测数据与基线数据多道匹配不确定的问题,采取了定向最佳匹配原则和失配虚拟技术方法.基于地震数据重复性与观测系统重复性的线性相关模型,研究了时移地震观测系统的多道加权均方根重复性特征.结合实际数据研究表明,基于面元的观测系统最佳定向匹配的多道加权均方根重复性基本反映全部观测系统的重复性和覆盖程度.