[发明专利]一种基于沉积地层单元的地震相分析方法

权利要求书

1.一种基于沉积地层单元的地震相分析方法，其特征在于，本发明利用三维地震属性数据体进行储层沉积特征的分析，运用高分辨率层序地层学原理和方法，以及地震属性分析方法和模式识别方法；

方法中采用K均值聚类算法，在沉积等时界面的控制下，进行沉积地层单元的沉积特征和沉积演化的分析过程。

2.根据权利要求1所述的一种基于沉积地层单元的地震相分析方法，其特征在于，

所述方法包括，

①.三维地震数据体的变换步骤，

在进行地震相分析之前，根据地震反射同相轴的参考地质年代，进行三维的等时沉积面的变换，首先根据井中地层界面的地质年代与深度的关系，在地震数据体中定义等时沉积面的追踪种子点；其次，根据种子点的地震属性特征，进行地震同相轴的追踪；然后，在地震同相轴追踪的基础上，构建等时沉积面的空间曲面；最后，根据等时沉积面的空间曲面，按照一定的地质年代步长，对三维地震数据体进行重采样和插值，得到能够反映同一沉积地层单元内的地震特征的数据体；

②.建立地震相的分类模式，

根据三维地震测网内不同时期沉积地层单元的特点，确定三维地震测网中地震相分析的地质年代分析时间步长，即时窗，结合井孔沉积相和成因地层单元的结果，建立不同沉积相带的地质年代域地震道的分类模式；根据研究区的地质特点，定义不同类型的沉积相的地震相模式，通过分类模式的建立，将地震属性特征与沉积相建立关联，并利用地震属性数据体根据分类模式进行同沉积时期的沉积相的自动分析；

③.采用K均值聚类算法实现沉积地层单元的地震相分析步骤，

K均值聚类算法用于确定时窗在三维属性体自动分析每一道的地震属性特征，从而得到地震相分类数据体；

其包括以下五个步骤：

i.根据井孔解释的沉积地层单元，为每个聚类确定一个初始的聚类中心，k个聚类存在k个聚类中心；

ii.将样本集中的每一个样本按照最小距离原则

D_i＝min{‖x-c_i‖}，x∈DataSet，i＝1，2，...，k    (1)

分配到k个聚类中的某一个类别；

iii.使用每个聚类中所有样本的均值作为新的聚类中心；

iv.若聚类中心在地质年代域的属性空间中产生位置的变化，则重复ii、iii步直到聚类中心不再变化为止；

v.最后得到k个聚类中心就是聚类的结果。

3.根据权利要求2所述的一种基于沉积地层单元的地震相分析方法，其特征在于，

所述K均值聚类算法的聚类准则为

J=Σj=1kΣi=1n||xi(j)-cj||2,xi(j)∈Sj,]]>c_j是聚类S_j的聚类中心                    (2)