韦红,白清云,张鹏志,甄宗玉

(中海石油(中国)有限公司 天津分公司渤海石油研究院，天津 300450)

0 引言

地震油水识别方法越来越广泛应用于勘探开发油田中，有关的油气检测技术快速发展。国内外研究人员利用叠前叠后地震资料研究流体检测方面取得了很多成果。Dilay A等讨论储层的内部及其顶、底的频率谱，分析含油、水对地震波频率的影响[1]。Stockwell R G等利用广义S变换时频分析法取得高精度的时频谱[2]。Goloshubin等[3-4]和Korneev V A等[5]建立考虑流体粘滞性和弥散性的声波方程，指导数值模拟，较好地描述了储层流体的频率依赖衰减特征。Korneev V A等在物理实验的基础上取得带衰减的地震记录，认为地层含流体后出现地震波的低频阴影、高频衰减异常和双程旅行时间延迟等等现象[5]。Goloshubin G等应用地震波的低频信息来识别孔隙介质中的流体性质[6]。国内研究人员利用对流体检测做了很多研究，并取得一些成果。卢明辉等基于双相介质理论进行数值模拟，分析油、气、水的地震波反射规律[7]。He Z H等寻找低频伴影及其能量强度随频率变化规律以识别流体[8]。孙万元等研究流体参数、频率对地震波衰减和频散的影响规律助力识别流体[9]。胡军辉等利用黏滞—弥散波动方程理论进行数值模拟，提出一种分析储层主频差、频带宽度差及频谱相似系数的流体检测方法[10]。张广智等在角道集资料基础上，基于Curvelet变换提取角度流体因子来区分流体信息[11]。黄捍东等和李红梅通过流体敏感参数分析，进行叠前反演，提高流体识别精度[12-13]。李相文等利用方位AVO分析技术削弱裂缝各向异性特征对油气响应的影响，提高油气预测精度[14]。李灿等利用纵横波速度比反演识别致密低渗储层中的流体[15]。但该油田是开发油田，钻井节奏很快，需要快速较准确地识别流体，帮助指导井位布署。由于该油田馆陶组的油属于较稠油，据文献研究其衰减吸收特征在叠后地震资料上会有较明显的特征，但是馆陶组的地震分辨率低，其储层特征是富砂特征，常规地震识别单砂体极困难，因此本文基于叠后地震资料的双相介质理论的地震波衰减理论，首先开展了反褶积广义S变换的时频分析，在高精度时频谱基础上利用油、水的频谱特征差异，提取单砂体储层顶、底的频率类属性，提高砂体识别精度，降低在时—空域中地震反演进行油、水层识别的多解性。

1 流体检测原理

1.1 反褶积广义S变换

由于渤海S油田馆陶组极其富砂的特征，要对单砂体进行准确识别顶底，当前原始地震资料的分辨率严重不足。因此，需要利用适合的时频分析方法进行改善。时频分析方法是进行地震信号的谱分解，检测频率信息的关键手段。在S变换的基础上，很多学者分别提出几种广义S变换方法，它们的区别在于窗函数不同，且不满足能量归一化条件，随频率增加，窗函数幅值会快速增大，时频分布的能量产生加权效应，从而得到错的能量分布特征。为此，本文利用陈学华等[16-17]代入高斯窗函数的调节参数λ>0、p>0，以改变高斯窗函数随着频率变化趋势，灵活分析和处理具体信号。

(1)

其中:

p(t,f)=x(T)e-i2πfT,

式中：T为高斯窗的控制时间，t为时间，f为频率。

广义S变换谱是高斯窗和原信号各自的Wigner-Vill分布的二维褶积，即：

WVDh(t-u,f-v)dudv

(2)

式中：WVDx、WVDh分别为原信号和高斯窗的Wigner-Vill分布。本文采用反褶积算法[18]计算广义S变换时频谱，即

(3)

1.2 双相介质理论

双相介质理论认为地下介质由固体和流体组成，且两种相态之间存在相对运动,并引起黏滞摩擦损耗[19]。据Biot理论，双相介质的波动方程表达式为[20-22]：

(4)

式(4)中考虑流体的影响。A、N为弹性参数，Q为常数，R为量度，θ为固体相应变，ε为流体相应变，t为时间，ρ11为单位体积中固体相的等效质量，u为介质的固体位移，ρ12为单位体积流体和固体间质量耦合系数，U为介质的液体位移，ρ22为单位体积流体相的等效质量，b为耗散系数,与流体黏滞系数、孔隙度和渗透率参数有关。

受里奇理论的启示，如果我们把量子力学的态迭加、坍缩现象引入摄影理论观测，那么可以比较容易发现和理解当代数码摄影的量子—数码世界观转向。令人惊奇的是，后现代理论家们如罗兰·巴特、雅克·德里达、让·鲍德里亚等卓越学者，早就站在新世界观（量子—数码世界观）的高点前瞻摄影的发展，牛顿式—传统摄影理论已经悄然向量子—数码世界观转换。这也是数码时代摄影实践转向的理论凝结。

(5)

式中：u0和U0为初始值，k′为复波数，k为波数，α为衰减系数，ω为圆频率，将式(5)代入式(4)，并简化之后得到：

(6)

其中，P=A+2N。

求解该式得到：

(7)

其中：

H=P+R+2Q，

ρ=ρ11+ρ22+2ρ12，

式(7)可见，衰减系数与圆频率近似成正比。式(4)为非黏弹性方程，说明地震波的衰减是双相介质固有特性。而地下流体的介质具有黏弹性，它会引起能量的衰减，使高频能量衰减更大。

2 模型试算

通过合成信号来进行时频分析，并证明广义S变换与反褶积广义S变换的不同的效果。设合成信号为x(t)，即

本次利用广义S变换和反褶积广义S变换得到相应的时频谱(图1)。从图中可见，反褶积广义S变换得到的时频谱能更好地区分频率发生突变的位置，其边界也更清晰，也即更高的分辨率和时频能量聚集。

图1 广义S变换时频谱(a)和反褶积广义S变换时频谱(b)

利用双相介质理论进行黏滞摩擦损耗的数值模拟(图2)，基于反褶积广义S时频分析后的地震剖面能较好地体现含流体储层在频率域的地震波衰减特性。图2b油层地震剖面显示，其顶底的频率差异较大；图2c水层的地震剖面显示，其顶底的频率差异弱于油层的。

图2 透镜体储层模型(a)及含油透镜体(b)、含水透镜体(c)

3 应用实例

3.1 油田概况

为了证实文中方法的有效性，本次对渤海S油田馆陶组的储层进行了反褶积广义S变换时频分析后，再进行频率域的储层顶底的频率差属性进行油水识别。

S油田的构造是一个在古潜山背景上发育起来的被断层复杂化的大型低幅度披覆构造，形成于渐新世，定形于上新世末期。馆陶组构造高点、主体部位构造特征与潜山构造面貌相近。馆陶组地层厚度137～459 m，极富砂，为辨状河沉积，其中馆陶组二段的主力油组为一套巨厚砂砾岩沉积，单井平均钻遇砂岩8、9层，平均砂岩含量89%。纵向上呈砂岩夹泥岩特征，储层平面分布稳定。储层岩性以中—细砂岩及粉砂岩为主；泥岩夹层的厚度约10 m，储层厚度约10～26 m不等。从井震标定结果来看，该主力油组的地震反射特征为较强振幅，地震同相轴较连续(见图3)。现有的地震资料品质较好，主频为40 Hz，地层平均速度为3 000 m/s，但是对于馆陶组的复合叠置砂体来说，其地震资料垂向分辨率不能准确识别储层或泥岩夹层的顶底界面，且振幅受干涉影响，并不能真实反映储层及流体的响应。

图3 常规的连井地震剖面

3.2 流体识别的现状

渤海S油田的馆陶组的复合叠置砂体油水关系复杂，严重制约井位部署。如馆陶组二段的主力油组1640砂体属于典型的复合叠置砂体，地质储量较大，目前综合调整部井10余口，但是第1口最有把握的I38P1遭遇失利，它离探井的距离仅仅300 m远，它的地震剖面与平面特征也与探井10井的非常相似(见图4、5)，10井为油层，而I38P1井确钻遇水层，两种不同的流体类型有相似的地震反射特征，这给后续井位部署带来极大不确定性。

图4 过I38P1-10井的地震剖面

据胡军辉等[10]基于黏滞—弥散波动方程理论，分析不同流体的黏滞系数对储层底的地震波形及频谱的影响，从数值模拟结果来看，在时间域里，黏滞系数越大，储层底的地震波的波形变化越明显；在频率域里，黏滞系数越大，储层底地震波频谱变化曲线越明显，且曲线差异大于时间域的差异[10]。该油田的原油具有密度高、黏度高、胶质沥青含量高、含蜡量中等、含硫量低及凝固点低的特征，属重质稠油。因此，在该油田的流体识别方法优选中，选取基于黏滞—弥散波动方程理论的频谱分析，提取储层顶底的频率信息差异，会对流体识别有较好的效果。然而由于地震分辨率低，较难准确识别储层顶底。

3.3 流体识别方法的应用

由于该油田的油层的黏度是变化的，为了减小流体识别的多解性，本文在黏滞—弥散波动方程理论的基础上，再结合双相介质理论同时研究流体变化。首先进行反褶积广义S变换时频分析，在频率域，计算高频属性体和低频属性体，据双相介质理论，低频吸收，高频衰减的动力学特征，利用已钻井进行反复尝试，计算得到地震波的一个低频体，其地震波衰减最小，振幅最大；同时计算一个高频体，地震波衰减最大，振幅最小。再将高频属性体与低频属性体作比值，得到一个新的频率属性体，若含油气，则该属性值会变小。同时基于高频属性体，较清楚地识别储层顶底，利用已钻井反复计算，最终提取得到储层顶底的主频差异属性，主频差异越大，含油气可能性越大。最后将储层顶底的主频差属性值比上新的频率属性体得到一个最终的频率属性体，以此放大流体差异，本文称这个最终的频率属性体为构建的该主力油组的流体识别因子，如式(11):

I=(fmb-fmt)/(Fh/Fl)

(11)

式中：fmb为储层底的主频，fmt为储层顶的主频，Fh为高频属性体，Fl为低频属性体。

本文利用构建的流体识别因子，首先检验失利井I38P1，从图6中可以看到，虽然I38P1井与10井距离很近，且从时间域里其地震波形相似，但是在频率域里，通过构建的流体识别因子—频率属性体里可以清楚地发现I38P1井的地震波的频谱变化特征截然不同，I38P1井的流体因子显示为弱属性值(0～35之间)，10井的为强属性值(35～55之间)，地震剖面和平面属性均有明显的区别(见图7)，这与实钻结果非常吻合。并且从流体识别因子的平面属性表明，I38P1井位于一条明显的分支河道中，与主体河道有不同的油水系统。这一发现为后续的井位布署提供了很多帮助，并在I38P1井北边500 m的位置布署一口水平井I38H井，并取得了成功。显然，I38H井位在常规地震振幅属性中(值35～55之间)也是属于有利区域(见图5)，不足以证明流体识别因子的可靠性。为此，我们在另一口探井附近布署一口评价性质的井J56H井，若评价成功，则转为采油水平井。钻前分析，从常规地震剖面和平面属性来看，该井离探井11井距离约60 m，构造位于较高部位，着陆点位置储层特征在常规地震剖面和平面属性中表现为较弱振幅属性(见图8a，图9a)，若着陆失败，则进行侧钻备选井位。而基于流体识别因子的地震剖面与平面属性中，其储层特征则表现为较好的地震属性(见图8b，图9b)。最终，J56H井成功着陆并钻遇油层，并且水平段也钻遇80%的优质储层。这口井的成功充分证明了，本次构建的流体识别因子在该油田的准确性和实用性，为后续调整井的实施提供了巨大帮助。

图5 主力油层的平面振幅属性

图6 基于流体识别因子过I38P1-10井地震剖面

图7 基于流体识别因子平面属性

图8 常规地震剖面(a)与基于流体识别因子的过11-J56H井的地震剖面(b)

图9 常规振幅平面属性(a)与基于流体因子的平面属性(b)

4 结论与认识

1)反褶积广义S变换在本文的时频分析中起到了提高分辨率和频率能量聚集的效果，适用于该油田的地震资料解释性处理。

2)在反褶积时频分析的基础上进行双相介质理论的应用，寻找到含油储层和含水储层的频率衰减特征，能较好地区分流体。并且利用新钻井的开发井再次验证了本次研究的实用性，助力该油田后期10余口调整井的布署和实施。

致谢:感谢本次一起工作的流体检测小组成员的辛勤付出，感谢专家和编辑部提出的宝贵意见。