[发明专利]一种基于条件递推的油藏描述的方法

说明书

技术领域

本发明涉及油藏描述领域，尤其涉及一种基于条件递推的油藏描 述的方法。

背景技术

油藏描述的理论及技术直接关系到是否能准确预测储量，对提高 采收率来说十分重要。现代油藏描述即是应用地质、物探、测井、测 试等多学科相关信息，以石油地质学、构造地质学、沉积学为理论基 础，以储层地质学、层序地层学、地震岩性学、油藏地球化学为方法， 以数据库为支柱，以计算机为手段，由复合型研究人员对油藏进行四 维定量化研究并进行可视化描述、表征及预测的技术。

由于油藏地质特征往往是千变万化的，要更加精确地反应真实的 地质特征，往往需要利用地质、物理、数学、概率统计等理论建立起 精确的油藏模型。

地震技术目前广泛应用于建立油藏描述模型中，它的基本原理是 利用记录地震波形，并反演出波阻抗值，根据反演波阻抗建立油藏描 述模型。目前通常利用测井数据、井间地震数据和三维地震数据来建 立油藏描述模型。但测井数据本身垂直分辨率高，但横向分辨率低； 而地震数据横向分辨率高，垂直分辨率低，为了匹配，通常将测井数 据进行粗化，这样建立的油藏描述模型分辨率相对不高，因此常常使 用插值法或模拟法进行处理，提高模型的分辨率。如图2所示的插值 法建立的油藏描述模型示意图和图3所示模拟法建立的油藏描述模型 示意图，与图1所示的真实油藏结构示意图相比，明显存在局部的较 大偏差。

发明内容

鉴于上述问题，本发明的目的是提供一种基于条件递推的油藏描 述的方法，使得油藏描述模型的分辨率更高，更加精确。

为了实现上述目的，本发明提供了一种油藏描述的方法，包括：

S11，获取初始油藏描述模型中各个网格的反演波阻抗；

S12，将所述各个网格细化成多个子网格并选定递推路径；

S13，根据当前待递推子网格所在的当前网格的反演波阻抗，以所 述当前网格的相邻网格的反演波阻抗和所述当前网格内邻近所述当前 待递推子网格的已递推子网格的波阻抗为约束条件，计算所述当前待 递推子网格的波阻抗；

S14，重复S13，直到递推得到所有子网格的波阻抗；

S15，根据所述所有子网格的波阻抗建立油藏描述模型；

所述S11进一步包括：

S111，根据地震数据得到地震反演波阻抗；

S112，根据测井数据得到测井波阻抗，并进行粗化；

S113，根据所述地震反演波阻抗和粗化后的测井波阻抗得到初始 油藏描述模型中各个网格的反演波阻抗；

S13中所述相邻网格的反演波阻抗包括所述当前网格的垂直上方 的相邻网格的反演波阻抗z₁和所述当前网格的垂直下方的相邻网格的 反演波阻抗z₂，所述当前网格内邻近所述当前待递推子网格的已递推 子网格的波阻抗包括与所述当前待递推子网格相距在预设距离内的已 递推子网格的波阻抗x_s；

所述S13进一步包括：

S131，对所述各个网格的反演波阻抗进行正态变换得到所述各个 网格的正态反演波阻抗；

S132，根据当前待递推子网格所在的当前网格的正态反演波阻抗 和所述波阻抗x_s，求解克里金系统得到高斯条件分布函数p(x_o|x_s)的均 值m_sk和方差其中x_o表示所述当前待递推子网格的正态波阻抗；

S133，根据所述当前网格的正态反演波阻抗和所述波阻抗x_s，求 解不同的克里金系统分别得到高斯条件分布函数f(z₁|x_s，x_o)的方差σ_f²和高斯条件分布函数g(z₂|x_s，x_o，z₁)的方差σ_g²；