一种基于MRF的油气储层岩相随机模拟及实现方法

摘要

本发明公开了一种基于MRF的油气储层岩相随机模拟及实现方法，利用MRF的条件概率公式，选取势函数为对数化的转移概率函数得到目标位置的岩相发生概率，利用邻域内的已知点和已模拟点，采用序贯模拟和条件模拟相结合的思想，依据转移概率函数图计算待模拟区域各类岩相的发生概率，利用蒙特卡洛随机模拟，得到该网格的岩相模拟结果。利用应用程序接口将所述方法以算法插件的形式整合到SGeMS平台中，可以实现在三维空间目标储层任意角度旋转并进行任意位置剖面切割的效果。由于转移概率方向性的存在，对于刻画岩相转移的单向分布趋势和各向异性特征，基于MRF的石油储层岩相预测的准确率得到显著提高。

主权项

一种基于MRF的油气储层岩相随机模拟及实现方法，其特征在于，包括以下几个步骤：步骤1：获取待模拟地区所选井段的油气储层的探井岩相属性数据；步骤2：对待模拟地区进行网格化处理，并进行岩相数字标记；根据待模拟地区的空间大小及岩相变化规律，对待模拟地区进行网格化处理，使得每个网格具有一种岩相属性，即每个网格具有唯一的岩相数字标记；步骤3：对待模拟地区网格，计算笛卡尔方向的各位置上岩相之间的转移概率，获取待模拟地区网格的笛卡尔方向的转移概率函数图；步骤4：依据步骤3得到的转移概率函数图，按照以下的目标条件概率公式计算待模拟区域网格中各目标位置属于各类岩相的发生概率：$P\{x_s=k|x_{s_1}=l_1,...,x_{s_N}=l_N\}=\frac{p_{l_{1}k}\left(h_1\right){\Pi }_{i=2}^N{p}_{{\mathrm{kl}}_i}\left(h_i\right)}{{\Sigma }_{f=1}^K\left[p_{l_{1}f}\left(h_1\right){\Pi }_{i=2}^N{p}_{{\mathrm{fl}}_i}\left(h_i\right)\right]}---\left(1\right)$所述目标条件概率公式是利用MRF的条件概率公式，选取势函数为对数化的转移概率函数得到；其中，s₁,…,s_N是目标位置s的邻居，且是按照与目标位置之间的距离从近到远排列的序列；l₁,…l_N,k,f分别为岩相类别数字标记，取值范围为{1,2,3,4,…,K}，对应于步骤1中获取的研究区域的岩相类型；h_i表示目标位置与目标位置的邻居点之间的距离，称为滞后距，按照距离从近到远依次取值；表示由目标位置的岩相类型k转化成间隔为h_i的邻居所对应的岩相类型l_i的转移概率；表示与目标位置距离最近的钻井位置s₁对应的岩相l₁转化为目标位置s所对应岩相k的转移概率，所述钻井位置s₁与目标位置s之间的距离为h₁；K表示所有的岩相种类总数；N表示目标位置的邻居个数；x_s是目标位置s处的岩相类别标志，x_r是位置r处的岩相类别标志。步骤5：根据步骤4计算得到的待模拟区域网格中各位置属于各类岩相的发生概率，利用蒙特卡洛随机模拟，得到岩相分布模拟图。