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

摘要

本发明涉及一种利用低频信息进行碳酸盐岩油气储层预测的方法，首先利用弥散—粘滞性波动方程对流体层正演，得到流体层的低频地震反射特征即低频强能量异常；其次将地震信号的单频时频域能量时间序列作为一个新的时间信号对其作用一个混合滤波器，再取瞬时最大时频域能量，这样得到的低频能量既异常突出又具有很高的时间分辨率。以上都是在低频段进行的，低频段的确定是以瞬时谱累加能量的15%—35%对应的频段。本发明利用弥散粘滞性波动方程做正演得出流体储层的低频地震反射特征即低频强能量异常；相比现有的低频技术本发明提取的低频强能量异常不但非常突出，而且具有很高的时间分辨，提高了碳酸盐岩油气储层预测的精度。

主权项

一种碳酸盐岩油气储层预测的方法，其特征在于,包括：(1)确定流体储层的低频响应特征：从弥散粘滞性波动方程确定低频强能量，低频强能量是主频向低频端移动的结果，这是由流体的流动性和弥散系数之间的关系决定的；(2)提取地震信号低频能量异常，自适应提取高分辨低频异常强能量：将地震信号时频域能量的单频能量时间序列重新看成一组时间信号，在此时间信号上作用一个复合滤波器；(3)避开低频分量的低时间分辨率，选择低频段：在瞬时谱上，油气储层与非储层高频段差异最大的频段是累加谱能量65％—85％对应的频段，将谱累加能量的35％到65％对应的频率段选择为提取低频强能量异常的默认低频段，所述的步骤(1)进一步包括：弥散粘滞性波动方程的表达式为：∂2u∂t2+f∂u∂t-η∂3u∂z2∂t-v2∂2u∂z2=0---(1)]]>其中，u为地震波场，f为弥散系数，表征流体的内摩擦系数，η为粘滞系数，t，z分别为时间和空间坐标；这是个变系数的偏微分方程，结合实际正演的需要，在很小的Δz内波场延拓，这样f，η，v就是常量；下面将此时间，空间域的波动方程转到频率、空间域，转换的算子为如下形式：∂∂t→-iw---(2)]]>(1)可以写为：∂2u∂z2+(w2+iwfv2-iwη)u=0---(3)]]>在Δz内其中kz＝k+iα   (5)k=sqrt(w2v2-wv2fv4+w2η2)---(6)]]>α=(w3η+wv2fv4+w2η2)/2k---(7)]]>频散公式为：v(wr)v(w)=1-1πQln(wwr)---(8)]]>其中，v表示地震波在地层中传播的速度；w表示地震波的角频率；kz表示求解公式(3)过程中的中间数学变量，是具有实部和虚部的复数，如公式(5)所示；α表示kz的虚部的系数，为中间数学变量；wr表示地震波的参考频率；Q表示地层对地震波的吸收衰减作用，也即，品质因子；所述的步骤(2)所述的滤波器由局部平均能量滤波器与局部高频滤波器的乘积组成，具体表达式为：f[x(n)]=13[x(n-1)+x(n)+x(n+1)][x(n)-x(n-1)+x(n)-x(n+1)]---(9)]]>对时频域低频段能量每一单频时间序列都作用(9)式，这样得到一个新的地震信号时频域能量，取每一个时刻取最大能量即可。