【数值模拟】各种方法汇总.docVIP

7 -（共 NUMPAGES 10页） 油藏数值模拟基本过程 一、数值模拟发展概况 30年代人们开始研究地下流体渗流规律并将理论用于石油开发； 50年代在模似计算的方法方面，取得较大进展； 60年代起步，人们开始用计算机解决油田开发上的一些较为简单间题，由于当时计算机的速度只有每秒几万到几十万次，实际上只能做些简单的科学运算； 70 年后主要体现于计算机的快速升级带动了油藏数模的迅猛发展，大型标量机计算速度达到100--500万次，内存也高增主约16兆字节。在理论上黑油模型计算方法更趋成熟，D. W. Peaceman的油藏数值模似基础以及K. Aziz和A. Settari的油藏模似等主要著作都是在这个阶段出版的，但仍受到计算机速度和内存的限制，使用的方法一般仅限于IMPES及半隐式等，只能解决中小型油藏的模拟应用问题； 80年代则是油藏数值模似技术飞跃发展的年代，解决不同类型油藏的数模计算方法及软件相应问世，同时超级向量机的诞生，使计算机速度达到亿次，甚至几十亿次，内存高达10—20亿字节。 90年代特别是后期，油藏模似软件各模块功能也有了惊人的发展，主要体现为向一体化方面发展；即集地震、测井、油藏工程（数模）、工艺及地面集输、经济评价等为一体的大型软件方面发展。 目前油藏数值模似软件基本上形成了一套能处理各种类型油气藏和各种不同开采方式的软件系列。 ?黑油模型已被广泛用于各种常规油气藏的模拟； ?裂缝模型可用来解决除砂岩以外的灰岩、花岗岩、凝灰岩和变质岩的裂缝性油气藏开发问题； ?组分模型用于凝析气藏、轻质油、挥发油藏的开发设计和混相驱的研究； ?热采模型用于稠(重)油油藏蒸气吞吐、蒸汽驱和就地燃烧的设计； ?化学驱模型用于在注入水中添加聚合物、表面活性剂、碱等各种化学剂进行三次采油提高采收率的计算和设计。 油藏数值模拟方法的新突破随着计算机运算速度的提高,向量算法的出现和应用是软件设计上一个划时代的发展。预处理共轭梯度法更快速、有效地解各种更为复杂和困难的大型稀疏线性方程组。网格化方面不局限于静态和动态的局部网格加密技术,不规则网格、PEBI网格的出现更好的解决了在边界、断层插值计算以及面与面垂直正交的新型数模计算方法，更快速收敛。此外多重网格法、混合有限元法、流线法等都在逐步完善和发展。同时，并行处理技术给大中型油田数模工作带来了生机。 二、数值模拟基本原理 以渗流力学为基础建立数值模型，即通过一组方程组，在一定的假设条件下，考虑油藏构造形态、断层位置、砂体分布、储层孔渗饱等参数的变化；流体高压物性变化；不同岩石类型；不同渗流驱替特征曲线（相渗）；井筒垂直管流等描述油藏真实的物理过程。主要包括：运动方程、状态方程和连续方程。 1、运动方程 引入单相（或多相）流的达西定律方程： 单相流达西定律：  多相流达西定律： 其中：K:空气渗透率；达西 A:截面积；厘米2 L:长度；厘米 △P:压差；大气压 Kro、Krw:油、水相对渗透率 2、连续性方程 研究流入流出单元体中质量的变化方程。根据物质平衡原理，流入单元体中的流体流量减去流出单元体流体流量等于单元体流体质量变化。 单元体中流体质量发生变化的速率 油 水 将达西定律代入连续性方程得出油、水的流动方程： 油 水 3、状态方程 为求解上述两个方程中的两个未知数P、Sw，引入状态方程：（1）（2）h、 、 为常数（3）（4） 、 为 的函数 So=1-Sw 由连续性方程推导出油、水的流动方程，从而引入压力方程： 饱和度方程引入含水百分数： 最后把两维两相流动方程写成如下形式： 解上述方程通过把微分方程离散化变代数形式来求解。 三、数值模拟的目的 （一）、为什么开展油藏数值模拟工作 研究和开发一个油田是一个复杂的综合性的科技问题，高精度的地震资料的处理解释提供研究区域的构造、断层、边界及其走向，但地震纵向分辨率受到限制，不能很好的反映一个同相轴（地震道）中沉积砂体的物性变化特征；测井可较好的反映到小于1米以下沉积砂体的物性特征，提供可靠的地层对比结果。但作为新老油田开发方案的研究及剩余油分布的研究，是地震、地质、测井理论方法都无法做到的。地质上仅定性或半定量分析，测井用于生产监测不能以点带面。惟独油藏数值模拟工作可再现生产历史，定量分析剩余油潜力；并做到室内研究投入少、时间短，还可进行开发方案优选及经济评价工作。所以总公司强调开发方案的部署一定要开展数值模拟工作。值得强调