[发明专利]一种基于地震波形结构属性的河流相砂体叠置区识别方法

权利要求书

1.一种基于地震波形结构属性的河流相砂体叠置区识别方法，其特征在于：其方法步骤为：

步骤一：进行资料准备，准备地质资料、地震资料和测井资料；

步骤二：筛选一周期时窗的目标层位，提取积分类、统计类和差分类的波形结构属性；

步骤三：通过积分类属性识别砂体分布较厚区域，找到砂体较厚的分布区域；

步骤四：凭借统计类属性识别砂体叠置区，如果厚砂体区域存在叠置砂体，那么在峰度、变异系数和偏度属性切片上，厚砂体区域中就会出现一些线条或纽带，而这些线条、纽带多是砂体叠置区反映出来的；

步骤五：依据差分类属性再次验证识别砂体叠置区，通过复合包络差、半时弯曲度差和峰谷峰度差属性切片再次验证识别河道砂体叠置形成的线条或纽带，同时和统计类属性切片作联立分析以确定出河道叠置区。

2.根据权利要求1所述的一种基于地震波形结构属性的河流相砂体叠置区识别方法，其特征在于：其中积分类包括波形面积、波形长度和波形平均弯曲度。

3.根据权利要求1所述的一种基于地震波形结构属性的河流相砂体叠置区识别方法，其特征在于：其中统计类包括波形结构峰度、波形变异系数和波形结构偏度。

4.根据权利要求1所述的一种基于地震波形结构属性的河流相砂体叠置区识别方法，其特征在于：其中差分类包括波形复合包络差、波形半时弯曲度差和峰谷峰度差。

5.根据权利要求2所述的一种基于地震波形结构属性的河流相砂体叠置区识别方法，其特征在于：其中波形面积为时窗内波形曲线的近似面积，计算公式为：

其中波形长度，时窗内波形曲线的近似长度，计算公式为：

其中，波形平均弯曲度指时窗内波形的长度与时窗长度之比，计算公式

为：

6.根据权利要求3所述的一种基于地震波形结构属性的河流相砂体叠置区识别方法，其特征在于：其中波形结构峰度是四阶中心矩除以标准差的四次方，用来度量数据在中心聚集程度，计算公式为：

其中波形变异系数为标准差与平均数的比值，计算公式为：

其中波形结构偏度是三阶中心矩除以标准差的三次方，是描述分布偏离对称性程度，计算公式为：

7.根据权利要求4所述的一种基于地震波形结构属性的河流相砂体叠置区识别方法，其特征在于：其中复合包络差是指上半时窗的平均反射强度与下半时窗的平均反射强度之差，计算公式为：

S_udr＝S_u-S_d。

8.根据权利要求4所述的一种基于地震波形结构属性的河流相砂体叠置区识别方法，其特征在于：其中半时弯曲度差是指上半时窗的波形平均弯曲度与下半时窗的波形平均弯曲度之差，计算公式为：

L_u＝(N/2)Δt；B_u＝W_u/L_u

L_d＝(N/2)Δt；B_d＝W_d/L_d；B_r＝B_u-B_d。

9.根据权利要求4所述的一种基于地震波形结构属性的河流相砂体叠置区识别方法，其特征在于：其中峰谷峰度差是波峰的峰度与波谷的峰度之差，计算公式为：

Ku_r＝W_P-Kurtosis-W_T-Kurtosis。