[发明专利]一种基于混沌属性的砂砾岩体空间展布预测方法及装置

权利要求书

1.一种基于混沌属性的砂砾岩体空间展布预测方法，其特征在于，该预测方法包括以下步骤：

1)将三维地震数据体中的地震反射振幅作为三维空间变量，计算任意一点P(x、y、t)地震反射振幅A(x、y、t)的梯度；任意一点P(x、y、t)地震反射振幅A(x、y、t)的梯度为：

▽A＝(▽x,▽y,▽t)^T

其中A表示三维地震数据、▽表示求梯度，上标T表示矩阵转置；

2)在给定范围内，计算出每个点的梯度向量，并建立相应的协方差矩阵；在定给范围内各点的梯度向量为：

其中N为时窗长度的样点数；

3)计算协方差矩阵的特征值，根据协方差矩阵的特征值确定混沌属性值J，根据混沌属性值J确定砂砾岩体的大小，以实现对砂砾岩体空间展布的预测；根据混沌属性值J确定砂砾岩体的大小的步骤为：

A)在相应的地质层位上、统计研究区内钻遇砂砾岩与未钻遇砂砾岩的井段，利用合成记录对目的层进行标定，确定砂砾岩体与非砂砾岩体在混沌属性体上对应的属性值；B)通过多口井的统计，将统计井的混沌属性值利用直方图进行分析，对于未能完全区分的混沌属性值范围，即重叠区，根据重叠区间求取二者的均值作为门槛值；C)根据求取的门槛值大小及混沌属性含义，大于门槛值的混沌属性部分为预测的砂砾岩体。

2.根据权利要求1所述的基于混沌属性的砂砾岩体空间展布预测方法，其特征在于，该方法还包括利用合成记录对目的层进行标定，对属性体资料进行解释，通过岩层属性提取对砂砾岩体进行描述。

3.根据权利要求1或2所述的基于混沌属性的砂砾岩体空间展布预测方法，其特征在于，混沌属性值J的计算公式为：

其中λ_max为协方差矩阵的最大特征值，λ_min为协方差矩阵的最小特征值，λ_mid为协方差矩阵的中间特征值。

4.根据权利要求1所述的基于混沌属性的砂砾岩体空间展布预测方法，其特征在于，横向上J值较大的区域反映砂体的平面展布，纵向上J值较小的区间反映砂砾岩体形成的间歇期。

5.根据权利要求1所述的基于混沌属性的砂砾岩体空间展布预测方法，其特征在于，在利用三维地震数据计算地震反射振幅梯度时需对三维地震数据进行叠后去噪处理。

6.一种基于混沌属性的砂砾岩体空间展布预测装置，其特征在于，该装置包括梯度计算模块、协方差矩阵构建模块和混沌属性值确定模块，

所述的梯度计算模块用于将三维地震数据体中的地震反射振幅作为三维空间变量，计算任意一点P(x、y、t)地震反射振幅A(x、y、t)的梯度；任意一点P(x、y、t)地震反射振幅A(x、y、t)的梯度为：

▽A＝(▽x,▽y,▽t)^T

其中A表示三维地震数据、▽表示求梯度，上标T表示矩阵转置；

所述的协方差矩阵构建模块用于在给定范围内，计算出每个点的梯度向量，并建立相应的协方差矩阵；在定给范围内各点的梯度向量为：

其中N为时窗长度的样点数；

所述的混沌属性值确定模块用于计算协方差矩阵的特征值，根据协方差矩阵的特征值确定混沌属性值J，并根据混沌属性值J确定砂砾岩体的大小；根据混沌属性值J确定砂砾岩体的大小的步骤为：

A)在相应的地质层位上、统计研究区内钻遇砂砾岩与未钻遇砂砾岩的井段，利用合成记录对目的层进行标定，确定砂砾岩体与非砂砾岩体在混沌属性体上对应的属性值；B)通过多口井的统计，将统计井的混沌属性值利用直方图进行分析，对于未能完全区分的混沌属性值范围，即重叠区，根据重叠区间求取二者的均值作为门槛值；C)根据求取的门槛值大小及混沌属性含义，大于门槛值的混沌属性部分为预测的砂砾岩体。

7.根据权利要求6所述的基于混沌属性的砂砾岩体空间展布预测装置，其特征在于，该装置还包括用于利用合成记录对目的层进行标定，对属性体资料进行解释，通过岩层属性提取对砂砾岩体进行描述的模块。

8.根据权利要求6或7所述的基于混沌属性的砂砾岩体空间展布预测装置，其特征在于，所述混沌属性值确定模块采用的混沌属性值计算公式为：

其中λ_max为协方差矩阵的最大特征值，λ_min为协方差矩阵的最小特征值，λ_mid为协方差矩阵的中间特征值。