《油藏数值模拟基础资料》由会员分享,可在线阅读,更多相关《油藏数值模拟基础资料(85页珍藏版)》请在金锄头文库上搜索。
1、油藏数值模拟基础 刘 月 田主编 石油大学(北京)石油天然气工程学院 目目录录 第一章第一章 绪论绪论 1 1 1 1 1.1油藏数值模拟的概念和作用 1 1.2油藏数值模拟研究的主要内容 2 1.3课程的特点、目的和要求 2 第二章第二章 油藏渗流的基本数学模型油藏渗流的基本数学模型 3 3 3 3 2.1数学模型的分类 3 2.2建立数学模型的步骤 3 2.3单相流的数学模型 4 2.4两相流的数学模型 6 2.5渗流数学方程的通式 8 2.6多组分模型 9 2.7黑油模型 11 2.8定解条件 13 第三章第三章 有限差分与差分方程有限差分与差分方程 15151515 3.1离散化概念
2、15 3.2有限差分方法 15 3.3差分方程的建立 17 3.4差分格式的稳定性分析 21 3.5边界条件的处理方法 22 3.6差分方程组的线性化 23 第四章第四章 线性代数方程组的数值解法线性代数方程组的数值解法25252525 4.1三对角方程组的直接解法追赶法 25 4.2GAUSS 消元法和直接分解法 26 4.3交替方向隐式方法(ADIP) 32 4.4简单迭代法和 Seidel 迭代法 33 4.5松驰迭代法(Successive Over Relaxation) 34 4.6交替方向隐式迭代方法(ADIIP) 36 4.7*强隐式方法 38 4.8迭代法和直接法的总结比较
3、44 第五章第五章 单相渗流数值模拟方法单相渗流数值模拟方法 45454545 5.1径向渗流问题的数值模拟方法 45 5.2平面渗流问题的数值模拟方法 51 第六章第六章 两相渗流数值模拟方法两相渗流数值模拟方法 56565656 6.1两相流完整的数学模型 56 6.2二维两相渗流差分方程的建立 56 6.3两相渗流差分方程的 IMPES 分离解法 58 6.4半隐式(联立)解法 61 6.5IMPIMS(分离)解法 63 6.6边界条件的统一化处理 64 6.7两相渗流数模程序编制 65 第七章第七章 现代油藏数值模拟技术及发展趋势现代油藏数值模拟技术及发展趋势68686868 7.1油
4、藏数值模拟技术发展简介 68 7.2现代油藏数模技术与大型数值模拟软件 68 7.3油藏数值模拟技术发展趋势 68 第八章第八章 油藏数值模拟应用油藏数值模拟应用 69696969 8.1问题分析及模型与软件的选择 69 8.2油藏几何参数处理与网格系统的建立 69 8.3地质参数的处理 69 8.4流体资料的处理 71 8.5油藏初始资料的处理 72 8.6生产动态资料 72 8.7数值模拟输出数据 73 8.8历史拟合与动态预测 73 第九章第九章 实用模拟软件及油藏模拟实例实用模拟软件及油藏模拟实例74747474 第一章第一章绪绪论论 1.11.11.11.1油藏数值模拟的概念和作用油
5、藏数值模拟的概念和作用 一、油藏数值模拟在油气田开发中的作用 油气藏存在和运动形式及其开发过程都十分复杂, 包括油藏构造、 岩石物性、 孔隙结构 、 非渗透夹层等静态地层描述的复杂性,包括流体形态、性质及其与岩石的相互作用,还包括 油田生产过程中油气产量、压力、温度和饱和度变化的复杂性,尤其是热力驱油、化学驱油 等“三次采油”过程的复杂性。 油气田开发的目的,就是排除各种复杂因素影响,克服各种困难,将地下岩石中的油气 流体开采到地面,并不断降低成本,提高最终采收率。经历过自喷、机械采油、注水开发、 化学及热力采油等阶段,油气田采收率不断提高。 油气田开发进程离不开科学技术。 利用现代科技手段不
6、断深入研究油气田、 认识油气田 , 是提高油气田开发效益的必经之路,研究油气田开发的方法有以下几类: 1直接法 如钻井、取芯、试采、试井、井组开发试验、先导开发区试验、工业化开发区试验等。 其优点是:直接、直观看得见,摸得着;准确本身即客观存在,不存在人为误 差。其缺点是:成本高,周期长;对油藏变化的不可逆性;只能在井点或油藏局部进 行。 2模拟法 数值模拟 解析方法 数学 物理 (1)物理模拟。如渗透率测试、相渗测试、毛管力测试、流体性质试验以及三采驱油 机理试验等。其作用是:发现微观运动规律,验证理论和数值研究结果;其缺点是: 很难进行复杂或整体油藏研究,成本较高,周期较长。 (2)解析方
7、法。通过数学方程组的解析解来研究其所描述的物理过程变化规律的方法。 优点是:简便,物理变化规律明显。缺点是:只能研究理想的、较简单的问题,很难直 接解决实际油气藏开发问题。 (3)数值模拟。根据油气藏地质及开发实际情况,建立物理数学模型并利用计算机求 得数值解来研究其运动变化规律。 其实质就是利用数学、地质、物理、计算机等理论方法技术对实际油藏的复制。 油藏数值模拟技术出现于 50年代, 到90 年代就已发展成熟。 现在处于另一个发展周期 。 油藏数值模拟功能包括两大部分: 复杂渗流力学研究, 实际油气藏开发过程整体模 拟研究,且可重复、周期短、费用低。 油藏数值模拟已经普遍应用于各种油气藏开
8、发过程, 成为油气田开发不可或缺的方法和 工具,被称作“现代油藏工程”。跟其它油藏研究方法相比,有着不可替代的优势。 油藏数值模拟的局限性主要在于:模拟误差,结果不唯一。误差主要来自两方面, 一是模型本身有误差, 二是油气藏资料不全或不准。 结果不唯一问题跟研究者的水平直接相 关。 二、从学科方向和课程设置看 科学研究主要靠三大手段:实验、理论、数值(模拟、实验、仿真) ,随着计算机技术 和计算科学的发展, 数值研究的地位越来越重要, 油藏数值模拟的功能和作用也随之不断增 强,被誉为“现代油藏工程”。 本课程是本专业同学在大学期间学习用计算机手段解决油气田开发问题的唯一的一门 课程,也是所有课
9、程中综合性最强的课程之一。 本课程虽没有被列入必修课, 但每个同学都应该修此课程。 一个石油工程专业的大学生 不懂得油藏数值模拟, 在今后的工作和学习中就缺少三大支柱研究手段之一, 且是非常重要 的一个,这将是一个很大的缺憾。 1.21.2 油藏数值模拟研究的主要油藏数值模拟研究的主要内容内容 计算软件编制 线性方程组求解 线性化 离散化 数值模型建立 数学模型建立 理论方法 1.31.3 课程的特点、目的和要求课程的特点、目的和要求 一、特点 (1)综合性强,需要宽厚的基础理论和专业知识,同时以方法为主线、应用为目的, 理论、方法和应用统一形成完整的有机体系。 (2)需较多的先修课程,故一般
10、在最后学习本课程。数学类:微分方程、线性代数, 专业类:石油地质、油藏物理、渗流力学、油藏工程和采油工程等;计算机类:计算机原理 、 程序语言等。 二、目的 (1)为了进一步利用数值手段研究复杂油气田开发问题打下基础理论方法。 (2)为从事实际工作做准备数模应用 三、课程要求 (1)了解掌握油藏数模的基本原理和方法 (2)能够编制简单的渗流计算程序 (3)掌握实际油藏数模及软件应用的方法步骤 (4)平时作业约 10 题,两个大作业:调试两相流程序和用现有软件进行实例模拟计算 第二章第二章油藏渗流的基本数学模型油藏渗流的基本数学模型 油藏内各种力的作用 质、各向异性)性质油藏介质(形状、非均 换
11、流体相态组合及物质交 油藏数学模型考虑因素 定解条件 微分方程(组) 数学模型的组成及要求 2.12.1数学模型的分类数学模型的分类 一、按流体相态数目分类 单相、两相、三相流模型 二、按油藏介质空间分类 1.零维模型: 油藏压力、 饱和度等量不随任何空间变量变化, 如油藏工程中物质平衡法 ; 2.一维模型:平行流(x) 、径向对称流(r) ; 3.二维模型:水平面内流动 x-y,垂直面内流动 x-z,柱状对称流动 r-z; 4.三维模型:实际油藏流动(x-y-z 或 r-z)。 三、按模型使用功能分类 1气藏模型:描述天然气藏及其运动(开发)过程的数学模型; 2黑油模型:油气水三相三组分,气
12、油间可发生溶解和分离,气水和油水间无交换; 3多组分模型:多组分三相模型,各相间有交换,如凝析气藏模型、聚合物驱模型; 4其它:热采模型(考虑温度变化、热能传递) , 聚合物驱(流体粘度发生变换) , 多重介质模型(裂缝和孔隙) , 各向异性模型(渗流介质的复杂性、张量渗透率) ; 总之,并没有严格统一的分类标准,约定俗成而已。 2.22.2建立数学模型的步骤建立数学模型的步骤 一、确定所要求解的问题(未知量 P、S) 二、问题满足的(隐含)条件:往往可以简化问题 三、基本方程和物理原理 四、定解条件 运动方程(达西定律) (非线性二项式) 状态方程(流体) (孔隙) 辅助方程(流动、参量、)
13、 化学、物理) 质量守恒方程 能量守恒方程 渗流偏微 分方程组 五、量纲: 如达西定律P K V= 六、适用范围 七、适定性 2 2.3 .3 .3 .3单相流的数学模型单相流的数学模型 推导单相流体三维渗流方程(数学模型) ,可压和微可压,恒等温。未知量(只有一个) P,运动方程和连续性方程、状态方程、无能量方程。 一、单相可压流体三维渗流 连续性方程 (1)()()()( t V z V y V x zyx = 运动方程(达西定律) (2) = = = )( )( )( z D g z pk V y D g y pk V x D g x pk V z y x 两方程联立,得 (3) )()( )()( tz D g z pk z y D g y pk yx D g x pk x = + + 渗流控制方程,D 为深度,垂直向下增加;g 重力加速度。 当 Z 垂直向下,Z 垂直向上,;同时有1= z D 1= z D 0= = y D x D 用 Hamilton 算子形式并考虑注入项写成: (4))()(