[发明专利]一种碳酸盐岩油气储层预测的方法

权利要求书

1.一种碳酸盐岩油气储层预测的方法，其特征在于,包括：

(1)确定流体储层的低频响应特征：从弥散粘滞性波动方程确定低频强能量，低频强能量是主频向低频端移动的结果，这是由流体的流动性和弥散系数之间的关系决定的；

(2)提取地震信号低频能量异常，自适应提取高分辨低频异常强能量：将地震信号时频域能量的单频能量时间序列重新看成一组时间信号，在此时间信号上作用一个复合滤波器；

(3)避开低频分量的低时间分辨率，选择低频段：在瞬时谱上，油气储层与非储层高频段差异最大的频段是累加谱能量65％—85％对应的频段，将谱累加能量的35％到65％对应的频率段选择为提取低频强能量异常的默认低频段，

所述的步骤(1)进一步包括：

弥散粘滞性波动方程的表达式为：

$\frac{{\∂}^{2}u}{\∂t^2}+f\frac{\∂u}{\∂t}-\mathrm{\η}\frac{{\∂}^{3}u}{\∂z^2\∂t}-v^2\frac{{\∂}^{2}u}{\∂z^2}=0---\left(1\right)$

其中，u为地震波场，f为弥散系数，表征流体的内摩擦系数，η为粘滞系数，t，z分别为时间和空间坐标；这是个变系数的偏微分方程，结合实际正演的需要，在很小的Δz内波场延拓，这样f，η，v就是常量；下面将此时间，空间域的波动方程转到频率、空间域，转换的算子为如下形式：

$\frac{\∂}{\∂t}\→-iw---\left(2\right)$

(1)可以写为：

$\frac{{\∂}^{2}u}{\∂z^2}+\left(\frac{w^2+iwf}{v^2-iw\mathrm{\η}}\right)u=0---\left(3\right)$

在Δz内

其中k_z＝k+iα (5)

$k=sqrt\left(\frac{w^2{v}^2-{\mathrm{wv}}^{2}f}{v^4+w^2{\mathrm{\η}}^2}\right)---\left(6\right)$

$\mathrm{\α}=\left(\frac{w^3\mathrm{\η}+{\mathrm{wv}}^{2}f}{v^4+w^2{\mathrm{\η}}^2}\right)/2k---\left(7\right)$

频散公式为：

$\frac{v\left(w_r\right)}{v\left(w\right)}=1-\frac{1}{\mathrm{\π}Q}\ln \left(\frac{w}{w_r}\right)---\left(8\right)$

其中，v表示地震波在地层中传播的速度；w表示地震波的角频率；k_z表示求解公式(3)过程中的中间数学变量，是具有实部和虚部的复数，如公式(5)所示；α表示k_z的虚部的系数，为中间数学变量；w_r表示地震波的参考频率；Q表示地层对地震波的吸收衰减作用，也即，品质因子；

所述的步骤(2)所述的滤波器由局部平均能量滤波器与局部高频滤波器的乘积组成，具体表达式为：

$f\[x\left(n\right)\]=\frac{1}{3}\[x(n-1)+x\left(n\right)+x(n+1)\]\[x\left(n\right)-x(n-1)+x\left(n\right)-x(n+1)\]---\left(9\right)$

对时频域低频段能量每一单频时间序列都作用(9)式，这样得到一个新的地震信号时频域能量，取每一个时刻取最大能量即可。

2.根据权利要求1所述的一种碳酸盐岩油气储层预测的方法,其特征在于，具体的技术实现包括：根据井上信息恰当选择小波基对每一道地震信号s(t)做时频分析,包括：小波变换，广义S变换，得到时频域能量分布，E(t,f)；

(4.1)选择提取低频强能量异常的瞬时低频段，得到E[t,(f₁,f₂)]；

(4.2)对E[t,(f₁-f₂)]每一个单频能量时间序列e(t)带入(9)式中的滤波器计算新的时频域能量E′[t,(f₁-f₂)]；

(4.3)取新能量E′[t,(f₁-f₂)]的瞬时最大值；

(4.4)取下一道返回第(4.1)步骤，

其中，f₁、f₂分别表示低频段的起始频率和终了频率。