专利名称：一种高时间分辨率的地层含油气性预测方法
技术领域：
本发明涉及地震勘探数据处理技术，是一种利用地震数据每个样点的频率属性，具有高时间分辨率的地层含油气性预测方法。
背景技术：
衰减是地震波在地下介质传播中总能量的损失，是介质内在的属性。在Domenic (1982)通过实验发现油气饱和度增大导致频率衰减这一现象后，Dilay和Eastwood(1995)分析了含油气层对地震信号频谱的影响。当地层含有油气后，会导致地震波的各种频率成分的能量都发生明显的衰减，高频成分衰减更快，且随含油气饱和度的增大，这种衰减程度更加明显。地震波各种频率的能量衰减程度反映地层的含油气情况，从而可以利用其进行含油气性预测。当前利用时频谱进行含油气性预测的主要步骤是对原始数据在频率域进行广义·S变换得到其每个样点的时频谱，将时频谱中总能量的65%和85%对应的频率间的时频谱进行E指数拟合，指数系数函数中的一次项系数即为该样点的吸收衰减梯度；依次求取每个样点的吸收衰减梯度，并最终确定吸收衰减梯度大的样点对应为含油气预测的有利区域。由于地震记录中高频信号的能量较弱，当有很强的高频干扰时，高频信号的信噪比较低，而当前吸收衰减梯度的计算又仅利用高频信号，仅靠高频吸收衰减梯度并不能很好的指示油气；同时，由于地震记录主要表现为中低频，目前的时频分析方法在低频段不能取得很高的时间分辨率。为突显低频信号的吸收衰减梯度，并获得很高的时间分辨率，需要一种对低频检测灵敏的时频分析方法。
发明内容
本发明目的是提供一种高计算效率的高时间分辨率的地层含油气性预测方法。本发明具体实现步骤包括I)采集记录地震原始数据；2)对原始数据在频率域进行广义S变换得到其每时间点的时频谱，并利用中值滤波对时频谱进行平滑，确定时频谱中瞬时主频和计算高频段吸收衰减梯度的起始频率和终止频率；步骤2)所述的时频谱中瞬时主频是每个时间点的时频谱中瞬时主频中时频谱最大值对应的频率，并将其作为计算高频段吸收衰减梯度的起始频率。步骤2)所述的计算高频段吸收衰减梯度的终止频率是时频谱中总能量的85%对应的频率。步骤2)所述的中值滤波是以误差的绝对值之和达到最小来确定滤波器输出响应的方法。3)对原始数据在频率域进行改进的广义S变换得到其每个点的时频谱，并利用中值滤波对时频谱进行平滑，确定时频谱中总能量的10-20%对应的频率为计算低频段吸收衰减梯度的起始频率，瞬时主频为计算低频段吸收衰减梯度的终止频率；步骤3)所述的改进的广义S变换的时窗函数G(t，f)为
_4]咖= *-(_☆)式中s和r均为调节因子,且s > 0, r > 0 ;f为频率，t为时间。4)对起始频率和终止频率间的时频谱p (Tk，f)取自然对数后利用最小二乘法进行 拟合，得到该时间点的低频段吸收衰减梯度和高频段吸收衰减梯度；步骤4)所述的吸收衰减梯度的计算公式为
I
_9] p{TkJ) = p(0J) exp[-\
\Jr J HeqX^k)
I^
\Jr J ^ieq k )式中p (Tk, f)为Tk时间，对应频率f的地震信号的时频谱；p (0, f)为地震子波的波峰处,对应频率f的地震信号的时频谱；
TkQeq(Tk) ~ ^ ^ dz ’ Qeq (Tk)是 Tk 处的等效 Q 值；
J。Q(T)Q(T)为时间T处对应的层Qfr为参考频率；Q1 为参考频率fj 的品质因子；s和r均为比例因子,取正数；Uin为地震子波的主频。5)重复步骤2)-4)直到所有道中每个时间点的低频段吸收衰减梯度和高频段吸收衰减梯度都计算完毕；6)确定低频段吸收衰减梯度和高频段吸收衰减梯度值都大的位置为含油气有利区域，得到油气预测结果。本发明在低频段具有更好的时间分辨率，并避开了高频干扰的影响，同时，其在高频段具有更好的时间分辨率，并避开低频干扰的影响，油气预测时间分辨率大大提高，且算法简单，计算效率高。
图I是传统的广义S变换时窗；图2是改进的广义S变换时窗；图3是计算吸收衰减梯度计算示意图；图4是动校后的VSP上行波；
图5是瞬时主频；图6是层速度；图7是高频吸收衰减梯度；图8是低频吸收衰减梯度。
具体实施例方式本发明高分辨率的油气检测主要利用一种改进的广义S变换，采用时窗宽度随频率f呈正比变化高斯窗函数，在低频段时窗较窄，获得很高的时间分辨率，求取地震波低频段吸收衰减梯度；同时利用广义S变换，采用时窗宽度随频率f呈反比变化高斯窗函数，在高频段时窗较窄，获得很高的时间分辨率求得高频吸收衰减梯度，并综合利用低频吸收衰减梯度和闻频吸收裳减梯度，从而提闻时间分辨率，以更精细的展不油气的空间展布。以下结合附图详细说明本发明。利用广义S变换和改进的广义S变换进行地震波吸收衰减梯度计算，并用于含油气性预测的具体实施步骤为I)采集记录地震原始数据；图4是井的VSP资料经过处理后的上行波拉平剖面，据钻井显示，含气层出现在
I.78-1. 79S 处；2)对原始数据在频率域进行广义S变换得到其每个时间点的时频谱，如图I所示，并利用中值滤波对时频谱进行平滑，确定时频谱中瞬时主频和计算高频段吸收衰减梯度的起始频率和终止频率；步骤2)所述的时频谱中瞬时主频是每个时间点的时频谱中瞬时主频中时频谱最大值对应的频率，并将其作为计算高频段吸收衰减梯度的起始频率。步骤2)所述的计算高频段吸收衰减梯度的终止频率是时频谱中总能量的85%对应的频率。步骤2)所述的中值滤波是以误差的绝对值之和达到最小来确定滤波器输出响应的方法。3)对原始数据在频率域进行改进的广义S变换得到其每个时间点的时频谱，如图2所示，并利用中值滤波对时频谱进行平滑，确定时频谱中总能量的15%或20%对应的频率为计算低频段吸收衰减梯度的起始频率，瞬时主频为计算低频段吸收衰减梯度的终止频率；如图3所示。步骤3)所述的改进的广义S变换的时窗函数G(t，f)为 °{tJ) = V2^yfexp(_2 77)式中s和r均为调节因子,且s > 0, r > 0 ;f为频率，t为时间。4)对起始频率和终止频率间的时频谱p (Tk, f)取自然对数后利用最小二乘法进行拟合，得到该时间点的低频段吸收衰减梯度和高频段吸收衰减梯度；步骤4)所述的吸收衰减梯度的计算公式为
权利要求
1.一种高时间分辨率的地层含油气性预测方法，特点是实现步骤包括 1)采集记录地震原始数据； 2)对原始数据在频率域进行广义S变换得到其每时间点的时频谱，并利用中值滤波对时频谱进行平滑，确定时频谱中瞬时主频和计算高频段吸收衰减梯度的起始频率和终止频率； 3)对原始数据在频率域进行改进的广义S变换得到其每个点的时频谱，并利用中值滤波对时频谱进行平滑，用确定的频率为计算低频段吸收衰减梯度的起始频率，瞬时主频为计算低频段吸收衰减梯度的终止频率； 4)对起始频率和终止频率间的时频谱P(Tk，f)取自然対数后利用最小二乗法进行拟合，得到该时间点的低频段吸收衰减梯度和高频段吸收衰减梯度； 5)重复步骤2)-4)直到所有道中每个时间点的低频段吸收衰减梯度和高频段吸收衰减梯度都计算完毕； 6)确定低频段吸收衰减梯度和高频段吸收衰减梯度值都大的位置为含油气有利区域，得到油气预测結果。
2.根据权利要求I所述的方法，特点是步骤2)所述的时频谱中瞬时主频是每个时间点的时频谱中瞬时主频中时频谱最大值对应的频率，并将其作为计算高频段吸收衰减梯度的起始频率。
3.根据权利要求I所述的方法，特点是步骤2)所述的计算高频段吸收衰减梯度的終止频率是时频谱中总能量的85%对应的频率。
4.根据权利要求I所述的方法，特点是步骤2)所述的利用中值滤波对时频谱进行平滑，中值滤波是以误差的绝对值之和达到最小来确定滤波器输出响应的方法。
5.根据权利要求I所述的方法，特点是步骤3)所述的改进的广义S变换的时窗函数G(t，f)为
6.根据权利要求I所述的方法，特点是步骤3)所述的用确定的频率是采用时频谱中总能量10-25%对应的频率。
7.根据权利要求I所述的方法，特点是步骤4)所述的吸收衰减梯度的计算公式为
8.根据权利要求I所述的方法，特点是步骤6)综合利用高频吸收衰减梯度和低频吸收衰减梯度进行地层含油气性预測。
全文摘要
本发明是一种高时间分辨率的地层含油气性预测方法，利用改进的广义S变换，采用时窗宽度随频率f呈正比变化高斯窗函数，在低频段时窗较窄，获得很高的时间分辨率，求取地震波低频段吸收衰减梯度；同时利用广义S变换，采用时窗宽度随频率f呈反比变化高斯窗函数，在高频段时窗较窄，获得很高的时间分辨率求得高频吸收衰减梯度，并综合利用低频吸收衰减梯度和高频吸收衰减梯度。本发明在低频段具有更好的时间分辨率，并避开了高频干扰的影响，同时，其在高频段具有更好的时间分辨率，并避开低频干扰的影响，油气预测时间分辨率大大提高，且算法简单，计算效率高。
文档编号G01V1/32GK102692647SQ20111007110
公开日2012年9月26日 申请日期2011年3月23日 优先权日2011年3月23日
发明者容娇君, 张固澜, 张彦斌, 李彦鹏 申请人:中国石油天然气集团公司, 中国石油集团东方地球物理勘探有限责任公司