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

摘要

本发明提供一种利用多方向变差函数分析进行储层属性模拟的方法，该利用多方向变差函数分析进行储层属性模拟的方法包括步骤1，根据区域沉积特征进行沉积微相划分；步骤2，在沉积微相划分的基础上，结合储层发育方向进行分区；步骤3，针对每个分区分别分析不同相带的变差函数；以及步骤4，在分区分相带的精细变差函数的基础上，进行储层孔隙度、渗透率等属性参数的模拟，得到储层的属性参数模型。该利用多方向变差函数分析进行储层属性模拟的方法解决了曲流河沉积砂体方向多变而变差函数单一的储层属性模拟问题，为河流相储层建模提供了切实可行的方法。

主权项

利用多方向变差函数分析进行储层属性模拟的方法，其特征在于，该利用多方向变差函数分析进行储层属性模拟的方法包括：步骤1，根据区域沉积特征进行沉积微相划分；步骤2，在沉积微相划分的基础上，结合储层发育方向进行分区；步骤3，针对每个分区分别分析不同相带的变差函数；以及步骤4，在分区分相带的精细变差函数的基础上，进行储层孔隙度、渗透率的模拟，得到储层的属性参数模型；在步骤1中，运用贝叶斯判别分析的方法定量划分沉积微相，分别建立每种沉积微相与特征参数之间的判别函数，对未知样品代入判别函数，可求得样品属于每个微相的条件概率，根据概率的相对大小判定样品属于哪个微相；在进行贝叶斯判别分析时：如果从G(G>2)种沉积微相a1、a2…aG中分别取出n1、n2…nG个样品，并且每个样品有m个特征参数，那么就构成了一个g×k的矩阵：Xgk=xgk(1)xgk(2)...xgk(m),(g=1,2,...G;k=1,2,...,ng)]]>式中为总体ag(g＝1,2,…,G)中第k(k＝1,2,…,ng)个样品第i(i＝1,2,…,m)个特征参数的观测值，根据样品数据可求得各类沉积微相的判别函数：Fg(X)=lnQg+Σk=1mCg(k)xg(k)+Cog,(g=1,2,...,G)]]>式中：G—沉积微相的类别数；g—沉积微相类别号；Qg—样本属第g类微相的先验概率，设各组先验概率相等，即Qg＝1/G；Cog—常数项；—第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μm2；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)]2式中x——空间中的一个点；h——其中的一个向量。