[发明专利]流体活动因子属性体的确定方法和装置

权利要求书

1.一种流体活动因子属性体的确定方法，其特征在于，包括：

获取目标区域的叠后地震资料；

根据所述叠后地震资料，求取瞬时振幅数据体、瞬时频率数据体；

通过差分法，分别计算瞬时振幅数据体对于时间的偏导数和瞬时频率数据体对于时间的偏导数；

根据改进的流体活动因子计算公式、所述瞬时振幅数据体对于时间的偏导数、所述瞬时频率数据体对于时间的偏导数、所述瞬时振幅数据体、所述瞬时频率数据体，确定所述流体活动因子属性体；

其中，根据改进的流体活动因子计算公式、所述瞬时振幅数据体对于时间的偏导数、所述瞬时频率数据体对于时间的偏导数、所述瞬时振幅数据体、所述瞬时频率数据体，确定所述流体活动因子属性体，包括：按照以下改进的流体活动因子计算公式，计算所述流体活动因子属性体：

上式中，M为流体活动因子属性体，为瞬时振幅数据体对于时间t的偏导数，为瞬时频率数据体对于时间t的偏导数，f为瞬时频率数据体中时间为t的数据。

2.根据权利要求1所述的方法，其特征在于，根据所述叠后地震资料，求取瞬时振幅数据体、瞬时频率数据体，包括：

对所述叠后地震资料进行希尔伯特变换，以求取所述瞬时振幅数据体、所述瞬时频率数据体。

3.根据权利要求1所述的方法，其特征在于，通过差分法，计算瞬时振幅数据体对于时间的偏导数，包括：

按照以下公式，计算所述瞬时振幅数据体对于时间的偏导数：

上式中，为瞬时振幅数据体对于时间t的偏导数，A(t+Δt)为瞬时振幅数据体中时间为t+Δt的数据，A(t-Δt)为瞬时振幅数据体中时间为t-Δt的数据，t为时间，Δt为相邻两个时间点之间的时间间隔。

4.根据权利要求1所述的方法，其特征在于，通过差分法，计算瞬时频率数据体对于时间的偏导数，包括：

按照以下公式，计算所述瞬时频率数据体对于时间的偏导数：

上式中，为瞬时频率数据体对于时间t的偏导数，f(t+Δt)为瞬时频率数据体中时间为t+Δt的数据，f(t-Δt)为瞬时频率数据体中时间为t-Δt的数据，t为时间，Δt为相邻两个时间点之间的时间间隔，f为瞬时频率数据体中时间为t的数据。

5.根据权利要求1所述的方法，其特征在于，在确定所述流体活动因子属性体后，所述方法还包括：

对所述目标区域进行井震标定，以确定目的层段的纵向位置；

根据所述目的层段的纵向位置，确定提取时窗；

利用所述提取时窗，从所述流体活动因子属性体中提取目的层段的流体活动性属性数据；

根据所述目的层段的流体活动性属性数据，确定目标区域中的油气展布。

6.根据权利要求1所述的方法，其特征在于，在根据所述叠后地震资料，求取瞬时振幅数据体、瞬时频率数据体后，所述方法还包括：

确定处理时窗；

利用所述处理时窗对所述瞬时振幅数据体、所述瞬时频率数据体中信噪比低于预设阈值的数据分别进行平滑处理。

7.一种流体活动因子属性体的确定装置，其特征在于，包括：

获取模块，用于获取目标区域的叠后地震资料；

求取模块，用于根据所述叠后地震资料，求取瞬时振幅数据体、瞬时频率数据体；

计算模块，用于通过差分法，分别计算瞬时振幅数据体对于时间的偏导数和瞬时频率数据体对于时间的偏导数；

确定模块，用于根据改进的流体活动因子计算公式、所述瞬时振幅数据体对于时间的偏导数、所述瞬时频率数据体对于时间的偏导数、所述瞬时振幅数据体、所述瞬时频率数据体，确定所述流体活动因子属性体；

其中，所述确定模块用于按照以下改进的流体活动因子计算公式，计算所述流体活动因子属性体中的数据：

上式中，M为流体活动因子属性体，为瞬时振幅数据体对于时间t的偏导数，为瞬时频率数据体对于时间t的偏导数，f为瞬时频率数据体中时间为t的数据。

8.根据权利要求7所述的装置，其特征在于，所述计算模块包括第一计算单元和第二计算单元，其中：

所述第一计算单元用于按照以下公式计算所述瞬时振幅数据体对于时间的偏导数：

上式中，为瞬时振幅数据体对于时间t的偏导数，A(t+Δt)为瞬时振幅数据体中时间为t+Δt的数据，A(t-Δt)为瞬时振幅数据体中时间为t-Δt的数据，t为时间，Δt为相邻两个时间点之间的时间间隔；

所述第二计算单元用于按照以下公式计算所述瞬时频率数据体对于时间的偏导数：

上式中，为瞬时频率数据体对于时间t的偏导数，f(t+Δt)为瞬时频率数据体中时间为t+Δt的数据，f(t-Δt)为瞬时频率数据体中时间为t-Δt的数据，t为时间，Δt为相邻两个时间点之间的时间间隔，f为瞬时频率数据体中时间为t的数据。