一种高时间分辨率的地层含油气性预测方法

摘要

本发明是一种高时间分辨率的地层含油气性预测方法，利用改进的广义S变换，采用时窗宽度随频率f呈正比变化高斯窗函数，在低频段时窗较窄，获得很高的时间分辨率，求取地震波低频段吸收衰减梯度；同时利用广义S变换，采用时窗宽度随频率f呈反比变化高斯窗函数，在高频段时窗较窄，获得很高的时间分辨率求得高频吸收衰减梯度，并综合利用低频吸收衰减梯度和高频吸收衰减梯度。本发明在低频段具有更好的时间分辨率，并避开了高频干扰的影响，同时，其在高频段具有更好的时间分辨率，并避开低频干扰的影响，油气预测时间分辨率大大提高，且算法简单，计算效率高。

主权项

一种高时间分辨率的地层含油气性预测方法，特点是实现步骤包括：1)采集记录地震原始数据；2)对原始数据在频率域进行广义S变换得到其每时间点的时频谱，并利用中值滤波对时频谱进行平滑，确定时频谱中瞬时主频和计算高频段吸收衰减梯度的起始频率和终止频率；所述的时频谱中瞬时主频是每个时间点的时频谱中最大值对应的频率，并将其作为计算高频段吸收衰减梯度的起始频率；所述的计算高频段吸收衰减梯度的终止频率是时频谱中总能量的85％对应的频率；3)对原始数据在频率域进行改进的广义S变换得到其每个点的时频谱，并利用中值滤波对时频谱进行平滑，用确定的频率为计算低频段吸收衰减梯度的起始频率，瞬时主频为计算低频段吸收衰减梯度的终止频率；4)对起始频率和终止频率间的时频谱p(T_k,f)取自然对数后利用最小二乘法进行拟合，得到该时间点的低频段吸收衰减梯度和高频段吸收衰减梯度；所述的吸收衰减梯度的计算公式为：$\begin{array}{c}p(T_k,f)=p(0,f)\exp [-{\left(\frac{f}{f_r}\right)}^{-\frac{1}{\pi Q_r}}\frac{\mathrm{\pi f}{T}_k}{Q_{\mathrm{eq}}\left(T_k\right)}]\\ =\exp (-2{\pi }^2{s}^2{f}^r{(f_{\mathrm{main}}-f)}^2)*\exp [-{\left(\frac{f}{f_r}\right)}^{-\frac{1}{\pi Q_r}}\frac{\mathrm{\pi f}{T}_k}{Q_{\mathrm{eq}}\left(T_k\right)}]\end{array}$式中：p(T_k,f)为T_k时间，对应频率f的地震信号的时频谱；p(0,f)为地震子波的波峰处，对应频率f的地震信号的时频谱；Q_eq(T_k)是T_k处的等效Q值；Q(t)为时间t处对应的层Qf_r为参考频率；Q_r为参考频率f_r的品质因子；s和r均为比例因子，取正数；f_main为地震子波的主频；5)重复步骤2)‑4)直到所有道中每个时间点的低频段吸收衰减梯度和高频段吸收衰减梯度都计算完毕；6)确定低频段吸收衰减梯度和高频段吸收衰减梯度值都大的位置为含油气有利区域，得到油气预测结果。