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

权利要求书

1.一种计算Gassmann流体项的方法，其特征在于，包括以下步骤，

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

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

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

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

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

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

2.如权利要求1所述的一种计算Gassmann流体项的方法，其特征在于，

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

S11、对获取的覆盖研究区的叠前地震角度道集进行解释性预处理；

S12、对步骤S11的叠前地震角度道集进行分角度叠加处理，得到第一至第三角度的叠加地震数据；

S13、基于井旁目的层段频谱分析结果，对步骤S12中第一至第三角度的叠加地震数据进行分频处理，获得第一至第三尺度的第一至第三角度的叠加地震数据体。

3.如权利要求2所述的一种计算Gassmann流体项的方法，其特征在于，步骤S12中，分角度叠加处理过程，包括以下步骤：

S121、获取研究区的地层速度模型数据；

S122、基于地层速度模型将叠前地震角度道集由时间域转化为深度域，利用目的层深度与最大偏移距长度进行偏移距与入射角度的转换，转化依据为最大偏移距长度的一半与目的深度的比值为入射角度的正切数值；

S123、综合分析理论最大入射角度、覆盖次数以及近道的多次波残留的因素，明确角度地震道集的有效入射角度范围区间，将有效入射角度区间进行三等分，形成第一入射角度等分范围、第二入射角度等分范围和第三入射角度等分范围，并将三个等分范围的中值确定为第一角度、第二角度、和第三角度；

S124、将位于第一入射角度等分范围内的叠前地震角度道集对应叠加、将位于第二入射角度等分范围内的叠前地震角度道集对应叠加、位于第三入射角度等分范围的叠前地震角度道集对应叠加，构成第一至第三角度的叠加地震数据。

4.如权利要求2所述的一种计算Gassmann流体项的方法，其特征在于，在步骤S13中，分频处理过程如下：

S131、基于井旁目的层段的地震频谱分析结果，明确有效频带宽度以及主频，以有效频带宽度的低频与主频的中值频率为第一尺度，主频为第二尺度，有效频宽的高频与主频的中值频率为第三尺度；

S132、利用匹配追踪时频分解方法对第一至第三角度的叠加地震数据进行分频处理，获得第一至第三角度的叠加地震数据的第一至第三尺度地震数据，分别为第一尺度第一角度地震数据、第二尺度第一角度地震数据、第三尺度第一角度地震数据、第一尺度第二角度地震数据、第二尺度第二角度地震数据、第三尺度第二角度地震数据、第一尺度第三角度地震数据、第二尺度第三角度地震数据、第三尺度第三角度地震数据。