[发明专利]利用多方向变差函数分析进行储层属性模拟的方法

权利要求书

1.利用多方向变差函数分析进行储层属性模拟的方法，其特征在于，该利用多方向变差函数分析进行储层属性模拟的方法包括：

步骤1，根据区域沉积特征进行沉积微相划分；

步骤2，在沉积微相划分的基础上，结合储层发育方向进行分区；

步骤3，针对每个分区分别分析不同相带的变差函数；以及

步骤4，在分区分相带的精细变差函数的基础上，进行储层孔隙度、渗透率的模拟，得到储层的属性参数模型；

在步骤1中，运用贝叶斯判别分析的方法定量划分沉积微相，分别建立每种沉积微相与特征参数之间的判别函数，对未知样品代入判别函数，可求得样品属于每个微相的条件概率，根据概率的相对大小判定样品属于哪个微相；

在进行贝叶斯判别分析时：

如果从G(G>2)种沉积微相a₁、a₂…a_G中分别取出n₁、n₂…n_G个样品，并且每个样品有m个特征参数，那么就构成了一个g×k的矩阵：

Xgk=xgk(1)xgk(2)...xgk(m),(g=1,2,...G;k=1,2,...,ng)]]>

式中为总体a_g(g＝1,2,…,G)中第k(k＝1,2,…,n_g)个样品第i(i＝1,2,…,m)个特征参数的观测值，

根据样品数据可求得各类沉积微相的判别函数：

Fg(X)=lnQg+Σk=1mCg(k)xg(k)+Cog,(g=1,2,...,G)]]>

式中：G—沉积微相的类别数；g—沉积微相类别号；Q_g—样本属第g类微相的先验概率，设各组先验概率相等，即Q_g＝1/G；C_og—常数项；—第g类微相中特征参数x^(k)前的系数；—第g类微相中第k个特征参数的平均值，对于新的待判样本X＝(x(1)x(2)…x(m))，将其代入Fg(X)，得F1(x)，F2(x)，…，FG(x)，若则判定样品X为判别函数值最大的第k类微相，它属于k类微相的条件概率为：式中(k＝1,2,…,G)；

选出5个特征参数建立了四种优势沉积微相的判别函数：

主河道：y＝2.707h+9.422φ-0.003k+279.007Md+6.048Vsh-211.458

河道边缘：y＝1.430h+9.648φ-0.003k+216.838Md+6.471Vsh-206.613

天然堤：y＝0.760h+9.093φ-0.003k+171.541Md+6.723Vsh-186.935

废弃河道：y＝0.991h+9.173φ-0.004k+284.500Md+6.191Vsh-196.771

式中：h—砂厚，m；φ—孔隙度，％；K—渗透率，10^-3μm²；Md—粒度中值，mm；Vsh—泥质含量，％；

在步骤2中，同一相带内砂体根据展布方向进行分区，同一分布方向作为一个分区；

在步骤3中，在分析变差函数时，设Z(x)是一个随机函数，如果差函数Z(x+h)-Z(x)的一阶矩和二阶矩仅依赖于点x+h和点x之差h，那么定义这一差函数的方差之半为变差函数γ(h)，

当Z(x)是一阶平稳时，变差函数写成下式：

γ(h)＝1/2*E[Z(x+h)-Z(x)]²

式中x——空间中的一个点；

h——其中的一个向量。

2.根据权利要求1所述的利用多方向变差函数分析进行储层属性模拟的方法，其特征在于，在步骤1中，在验证判别函数的性能时，将取心井样品代入判别函数，根据判别函数值确定每个样品分属哪一类，通过以上判别函数对研究区馆上段多个主力层进行微相定量划分，根据井点判别结果勾绘了沉积微相平面图，对与地质认识不一致的井点进行了人工校正。