[发明专利]一种基于正态分布的薄互层砂厚预测方法

权利要求书

1.一种基于正态分布的薄互层砂厚预测方法，其特征在于，包括步骤：

利用油田已有钻井资料建立井点砂厚-地震敏感属性交会图，其中横轴为井上统计的砂厚，纵轴为提取的井点地震敏感属性；

根据交会图分别统计砂厚均值、地震敏感属性均值、砂厚与地震敏感属性的协方差矩阵信息；

利用所述砂厚均值、地震敏感属性均值以及砂厚与地震敏感属性的协方差信息绘制正态分布函数；

输入新的地震敏感属性值至所述正态分布函数，根据所述正态分布函数逐样点计算每个样点的地震敏感属性值所对应的具有最大正态分布值的砂厚以及最大正态分布概率；

按照样点位置依次输出具有最大正态分布值的砂厚以及最大正态分布概率，形成最大概率砂厚属性与砂厚可信度属性。

2.根据权利要求1所述的方法，其特征在于，所述利用统计得到的砂厚均值、地震敏感属性均值、砂厚与地震敏感属性的协方差矩阵绘制正态分布函数的具体方法为：

以砂厚-地震敏感属性交会图为底图，将砂厚与地震敏感属性均匀剖分成M×N的网格，其中M为砂厚剖分的段数，N为地震敏感属性剖分的段数；

利用所述正态分布函数计算所述网格中每个样点的正态分布值并以散点形式绘制在底图上，根据散点的正态分布值勾绘出K条等值线，其中K为预先定义的等值线条数，即得到高斯概率分布平面图。

3.根据权利要求1所述的方法，其特征在于，所述砂厚与地震敏感属性的协方差矩阵的计算公式为：

式(1)中Σ为协方差矩阵，x_thick为砂体砂厚组成的列向量，x_attr为地震地震敏感属性组成的列向量，cov(·)为协方差计算函数：

cov(x_i,x_j)＝E[(x_i-E|x_i|)(x_j-E|x_j|)] (2)

式(2)中x_i,x_j代表x_thick或x_attr变量，E(·)为求取期望值的函数。

4.根据权利要求3所述的方法，其特征在于，所述正态分布函数

其中f(x)为正态分布函数，x＝[x_thick,x_attr]^T为由砂体砂厚与地震敏感属性组成的列向量，μ＝[μ_thick,μ_attr]^T为由砂厚均值μ_thick与地震敏感属性均值μ_attr组成的列向量，Σ为协方差矩阵。

5.根据权利要求1所述的方法，其特征在于，所述根据正态分布函数逐样点计算每个样点的地震敏感属性值所对应的具有最大正态分布值的砂厚，以及最大正态分布概率包括；

砂厚-地震敏感属性交会图，预先确定地震敏感属性值上限a与下限b，以小步长s将a至b的属性值剖分成L段，则对应L+1个节点；

通过遍历的方式计算每个节点地震敏感属性值所对应的最大概率的砂厚与最大正态分布概率；

逐样点计算新输入地震属性中任一样点c的最大概率的砂厚与最大正态分布概率时，取第(c-a)/L+1个所述节点对应的最大概率的砂厚与最大正态分布概率作为样点c的预测结果，从而实现最大概率的砂厚与最大正态分布概率的快速预测。