[发明专利]储层地震波衰减特性分析方法及储层反射系数的反演方法

权利要求书

1.一种储层地震波衰减特性分析方法，其特征在于，包括以下步骤：

基于目标储层的非稳态地震反射系数反演模型，根据目标储层的地震记录和与各预设品质因子分别对应的时变子波库，确定与各预设品质因子分别对应的反射系数序列，并确定各反射系数序列的Lp范数；

从各反射系数序列的Lp范数中确定最小Lp范数，将最小Lp范数的反射系数序列所对应的预设品质因子作为目标品质因子；

利用所述目标品质因子分析目标储层的地震波衰减特性。

2.根据权利要求1所述的方法，其特征在于，在基于目标储层的非稳态地震反射系数反演模型，根据目标储层的地震记录和与各预设品质因子分别对应的时变子波库，确定与各预设品质因子分别对应的反射系数序列，并确定各反射系数序列的Lp范数之前，还包括以下步骤：

构建目标储层的非稳态地震反射系数反演模型。

3.根据权利要求2所述的方法，其特征在于，构建目标储层的非稳态地震反射系数反演模型，具体包括以下步骤：

通过非稳态地震数据的正演模拟，将地层Q滤波效应加入多道稳态地震反射波褶积模型中，以获取目标储层的非稳态地震反射系数反演模型。

4.根据权利要求3所述的方法，其特征在于，所述目标储层的非稳态地震反射系数反演模型如下：

s＝G(Q)r

其中，s为检波器检测的地震记录，G(Q)为与各预设品质因子分别对应的时变子波库，r为非稳态地震反射系数。

5.根据权利要求4所述的方法，其特征在于，通过以下表达式获取与各预设品质因子分别对应的时变子波库：

G(Q)＝F^HWA

其中，G(Q)为与各预设品质因子分别对应的时变子波库，F^H为反傅立叶算子，W为频率初始地子波，A为与各预设品质因子有关的衰减因子。

6.根据权利要求1所述的方法，其特征在于，通过以下表达式确定各反射系数序列的Lp范数：

其中，L_p为反射系数序列的Lp范数，r_j,i为反射系数序列，L为反射系数序列的横向样点总个数，N为反射系数序列的纵向样点总个数，p为范数。

7.根据权利要求1所述的方法，其特征在于，从各反射系数序列的Lp范数中确定最小Lp范数，将最小Lp范数的反射系数序列所对应的预设品质因子作为目标品质因子，具体包括以下步骤：

第一步骤，判断当前反射系数序列的Lp范数是否小于上一个反射系数序列的Lp范数：

若是，则存储当前反射系数序列的Lp范数；

若否，则存储上一个反射系数序列的Lp范数；

第二步骤，判断是否已遍历所有反射系数序列的Lp范数；

若是，则将存储的Lp范数的反射系数序列所对应的预设品质因子作为目标品质因子；

若否，则将下一个反射系数序列作为当前反射系数序列，并返回第一步骤。

8.根据权利要求1所述的方法，其特征在于，预设品质因子是通过以下步骤获取的：

在给定的预设品质因子的取值范围内，以给定的间隔步长选取多个间隔均匀的品质因子作为预设品质因子。

9.一种储层反射系数的反演方法，其特征在于，包括以下步骤：

按照权利要求1至8中任一项所述的储层地震波衰减特性分析方法获得目标品质因子；

根据目标品质因子对目标储层的反射系数进行反演。

10.一种存储介质，其特征在于，其上存储有用于实现如权利要求1至8中任一项所述储层地震波衰减特性分析方法的步骤的计算机程序，或者用于实现如权利要求9所述储层反射系数的反演方法的步骤。