油藏Voronoi网格化的研究.pdfVIP

第32卷 第 1期 西南石油大学学报 (自然科学版) Vo1．32 No．1 2010年 2月 JournalofSouthwestPetroleumUniversity(ScienceTechnologyEdition) Feb． 201O 文章编号：1674—5086(2010)01—0121—04 油藏 Voronoi网格化的研究 陈晓军h，陈伟 ，段永刚k，刘柏峰 ，张 静岫 (1．西南石油大学：a．现代教育技术中心，b．计算机科学学院，c．石油工程学院，四川 成都 610500 2．西南石油大学高等职业技术学院，四川 南充637001) 摘 要：Voronoi网格具有很强的灵活性和局部正交等特点，可以把局部细分网格与基本粗化网格有效地结合在一起， 满足了油藏数值模拟技术对于大范围粗格化模拟、局部精细化描述的需要。对Voronoi网格的生成技术进行了分析和研 究，提出了一种新的生成方法：把油藏区域按断层、井层、基层分层，针对不同层，采用不同的几何模型布点，最后采用间 接法构造vomnoi网格。从计算几何的角度把层干扰归纳为点与点、点与线和线与线 3类，并给出了相关的解决措施，重 点对断层与垂直井、水平井的干扰进行处理。通过Vomnoi网格化模型，验证 了该算法的正确性和有效性。 关键词 ：油藏数值模拟 ；Voronoi网格 ；几何模型；间接法；层干扰 中图分类号：TE319 文献标识码 ：A DOI：10．3863／j．issn．1674—5086．2010．01．023 引 言 油藏数值模拟是应用计算机研究油、气藏中多 相流体渗流规律的新型数值计算方法，在油、气 田开 发方面有着非常广泛的应用，是提高油、气 田开发水 平的一项重要技术 j。油藏数值模拟利用的是偏 微分方程数值解法求解模型方程 ，涉及数值离散化 的问题，即把 求解 区域 分解 成有 限 的求解 网 格 。因此分解网格技术即网格剖分技术在油 图 1 Voronoi网格 藏数值模拟中占有非常重要的地位。随着计算机技 Fig．1 Voronoigrid 术的飞速发展以及油藏数值模拟研究的进展，网格 技术也进一步得到发展。由于结构化网格技术不能 1 Voronoi网格的常规生成方法 很好地模拟断层，同时不能很好地体现流体流向，非 结构化网格技术越来越得到人们的重视。 Voronoi网格本质上就是计算几何中的Voronoi Heinemann… 1988年提出了Voronoi网格概念并 图。构造Voronoi图的算法很多文献 [12—14]都有 对其可行性进行了研究，采用Voronoi网格模拟油藏 介绍，大致可分为直接法和间接法。直接法包括半 流动。如图1所示，Voronoi网格定义了这样一个空间 平面的交、增量法、分治法、减量算法和平面扫描算 区域：在此空间区域内的任意一点到网格中心点(节 法等；间接法的实质是根据Voronoi图和Delaunay三 点)的距离均小于此点到其他节点的距离。在此定义 角形的对偶性，先对离散点进行Delaunay三角剖分， 下，任意两节点之间的连线必然被两节点所在网格边 然后在求出的三角形集 中求出外心，并连接相邻三 界垂直平分。结果表明，Voronoi网格对于复杂形态 角形 (两三角形有共 同的边 )的外心从而构成 的油藏描述较灵活，模拟结果较为准确，可减少网格 Voronoi图。在实际应用中由于直接法没有固定的 取向效应；其正交特点，符合径向流的方向，Voronoi网 数据结构，效率比较低，而构造Delaunay三角形的算 格有望在数值试井中得以较多的应用。 法相当成熟，因此从复杂性 、计算效率和实现难度方 收稿 日期 ：2009—05—02 作者简介：陈晓军(198l一)，男(汉族)，江苏大丰人，硕士，主要从事计算机应用方面的工作。 122 西南石油大学学报 (自然科学版) 面