[发明专利]一种计算Gassmann流体项的方法及系统

说明书

本发明公开了一种计算Gassmann流体项的方法及系统，包括以下步骤：S1获取不同尺度的不同角度的叠加地震数据体；S2建立第一至第三角度的流体弹性阻抗低频模型；S3按由第一至第三尺度顺序逐渐进行第一至第三角度的流体弹性阻抗反演，并以第一尺度反演结果为先验模型约束第二尺度反演的寻优过程，以第二尺度的反演结果为先验模型约束第三尺度的寻优过程，最终获得第一至第三角度的流体弹性阻抗反演结果；S4对步骤S3中获得的第一至第三角度流体弹性阻抗反演结果进行对数变换，结合测井统计信息构建Gassmann流体项的求解方程，求解方程组并获得地震尺度的Gassmann流体项。

技术领域

本发明涉及一种基于逐尺度渐进反演策略计算Gassmann流体项的方法及系统，属于地球物理勘探技术领域。

背景技术

随着岩石物理理论、地震处理技术以及叠前地震反演算法的发展进步，利用弹性参数进行储层烃类检测已经成为油气工业的常用方法。在实际生产应用中，为了提高烃类检测的精度，如何构建对地下孔隙介质流体类型敏感的弹性参数以及对其稳定可靠的计算一直是学者们研究的方向。

Gassmann流体项是Russell等人在2006基于Biot-Gassmann方程推导出的一种弹性参数，该弹性参数可以较好的降低地下孔隙介质的孔隙度参数造成的烃类检测假象，从而较为直接指示地下介质孔隙流体类型，具有敏感的储层烃类判识敏感性。目前，相关业内多以叠前地震为主体，以已知地质规律、测井资料为先验约束，采用流体弹性阻抗反演技术计算Gassmann流体项。然而，随着地下目标越来越复杂，常规的叠前反演方法已经难以满足实际生产的需求，尤其是针对隐蔽岩性、深层以及河流-三角洲相薄储层等类型的油气藏，这些类型的储层厚度一般小于目的层段的地震分辨率，且空间非均质性比较强，加之常规反演方法的不适定性进一步影响了全局最优解的计算，使得利用常规反演方法计算的Gassmann流体项参数难以满足该类复杂地质目标的烃类检测需求，进而影响储层评价质量。

发明内容

针对上述问题，本发明的目的是提供一种基于逐尺度渐进反演策略计算Gassmann流体项的方法及系统。

为实现上述目的，本发明采用以下技术方案，一种计算Gassmann流体项的方法，包括以下步骤，

S1、获取不同尺度的不同角度的叠加地震数据体，其中不同尺度包括依次递减的第一、第二、第三尺度，不同角度包括依次递增的第一、第二、第三角度；

S2、建立第一至第三角度的流体弹性阻抗低频模型；

S3、按由第一至第三尺度顺序逐渐进行第一至第三角度的流体弹性阻抗反演，并以第一尺度反演结果为先验模型约束第二尺度反演的寻优过程，以第二尺度反演结果为先验模型约束第三尺度的寻优过程，最终获得第一至第三角度的流体弹性阻抗反演结果；

S4、对步骤S3中获得的第一至第三角度流体弹性阻抗反演结果进行对数变换，结合测井统计信息构建Gassmann流体项的求解方程，求解方程组并获得地震尺度的Gassmann流体项

构建Gassmann流体项的求解方程如下：

式中，EI表示流体弹性阻抗；t表示时间；θ₁表示第一角度；θ₂表示第二角度；θ₃表示第三角度；θ表示入射角度；表示饱和流体岩石纵横波速度比的平方；表示干岩石纵横波速度比的平方；f、μ、ρ分别表示Gassmann流体项、剪切模量和密度参数；f₀、μ₀和ρ₀分别定义为Gassmann流体项，剪切模量和密度参数的平均值；EI₀是标准化因子。

进一步地，在步骤S1中，获取不同尺度的不同角度的叠加地震数据体，具体过程包括：