油藏数值模拟方法流程

油藏数值模拟方法流程（新手上路）

油藏数值模拟

实用方法文字

前言：

油藏数值模拟是随着计算机的发展，而在石油行业中逐步成为一门成熟的技术。追溯油藏数值模拟的发展史，从30年代开始研究渗流力学到50年代在石油工业方面得以应用，到70年代进入商品化阶段，而80年代油藏数值模拟又向完善、配套、大型多功能一体化综合性软件飞跃发展。近十年油藏数值模拟已成为油田开发研究，解决油田开发决策问题的有力工具。在衡量油田开发好坏、预测投资、对比油田开发方案、评价提高采收率方法等方面应用都极为广泛。 油藏数值模拟就是应用数学模型再现实际油田生产动态。具体通过渗流力学方程借用大型计算机，结合地震、地质、测井、油藏工程学等方法在建立的三维地层属性参数场中，对数学方程进行求解，实现再现油田生产历史，解决油田实际问题。

油藏数值模拟是一门综合性很强的科学技术，涉及油田地质、油层物理、油藏工程、采油工程、测井、数学、计算机及系统等学科。而油藏数值模拟工作又以其繁重的前期准备和上机历史拟合运算工作让人望而生畏。

那么如何做好前期资料准备工作和尽快掌握模拟技巧？使得今后的油藏数值模拟工作在作业区顺利开展，便是出此书的目的所在。

本书结合以往工作中的实际经验教训，成功与失败，参考诸多资料从前期数据准备工作开始到模拟技巧做了较为的详细介绍，以舐读者。有不妥之处，请予指证。同时，今后不定期的将更新的模拟技术及方法推荐给大家。

目 录

一、数值模拟发展概况 二、数值模拟的基本原理

二、选择适当的数值模型及相类

三、数据录取准备工作

（一） 建立油藏地质模型 （二） 网格选择 （三） 数据录入准备

四、历史拟合方法及技巧

1

（一）确定模型参数的可调范围 （二）对模型参数全面检查 （四） 历史拟合

附件1：关于实测压力的皮斯曼校正 附件2：关于烃类有效孔隙体积的计算

一、 数值模拟发展概况

30年代人们开始研究地下流体渗流规律并将理论用于石油开发； 50年代在模似计算的方法方面，取得较大进展；

60年代起步，人们开始用计算机解决油田开发上的一些较为简单间题，由于当时计算机的速度只有每秒几万到几十万次，实际上只能做些简单的科学运算； 70年后主要体现于计算机的快速升级带动了油藏数模的迅猛发展，大型标量机计算速度达到100--500万次，内存也高增主约16兆字节。在理论上黑油模型计算方法更趋成熟，D. W. Peaceman的<油藏数值模似基础>以及K. Aziz和A. Settari的<油藏模似>等主要著作都是在这个阶段出版的，但仍受到计算机速度和内存的限制，使用的方法一般仅限于IMPES及半隐式等，只能解决中小型油藏的模拟应用问题；

80年代则是油藏数值模似技术飞跃发展的年代，解决不同类型油藏的数模计算方法及软件相应问世，同时超级向量机的诞生，使计算机速度达到亿次，甚至几十亿次，内存高达10—20亿字节。

90年代特别是后期，油藏模似软件各模块功能也有了惊人的发展，主要体现为向一体化方面发展；即集地震、测井、油藏工程（数模）、工艺及地面集输、经济评价等为一体的大型软件方面发展。

目前油藏数值模似软件基本上形成了一套能处理各种类型油气藏和各种不同开采方式的软件系列。

?黑油模型已被广泛用于各种常规油气藏的模拟；

?裂缝模型可用来解决除砂岩以外的灰岩,花岗岩,凝灰岩和变质岩的裂缝性油气藏开发问题；

?组分模型用于凝析气藏,轻质油和挥发油藏的开发设计和混相驱的研究； ?热采模型用于稠(重)油油藏蒸气吞吐,蒸汽驱和就地燃烧的设计；

?化学驱模型用于在注入水中添加聚合物,表面活性剂,碱等各种化学剂进行三次采油提高采收率的计算和设计。 油藏数值模拟方法的新突破

随着计算机运算速度的提高,向量算法的出现和应用是软件设计上一个划时代的发展。预处理共轭梯度法更快速,有效地解各种更为复杂和困难的大型稀疏线性方程组。网格化方面不局限于静态和动态的局部网格加密技术,不规则网格、PEBI网格的出现更好的解决了在边界、断层插值计算以及面与面垂直正交的新型数模计算方法，更快速收敛。此外多重网格法、混合有限元法、流线法等都在逐步完善和发展。同时，并行处理技术给大中型油田数模工作带来了生机。

二、数值模拟基本原理

以渗流力学为基础建立数值模型，即通过一组方程组，在一定的假设条件下，考

2

虑油藏构造形态、断层位置、砂体分布、储层孔渗饱等参数的变化；流体高压物性变化；不同岩石类型；不同渗流驱替特征曲线（相渗）；井筒垂直管流等描述油藏真实的物理过程。主要包括：运动方程、状态方程和连续方程。 1．运动方程

引入单相（或多相）流的达西定律方程： 单相流达西定律： 多相流达西定律： 其中：K:空气渗透率；达西 A:截面积；厘米2 L:长度；厘米 △P:压差；大气压

Kro、Krw:油、水相对渗透率 μo:原油粘度；厘泊 转变为运动方程： 2、连续性方程

研究流入流出单元体中质量的变化方程。根据物质平衡原理，流入单元体中的流体流量减去流出单元体流体流量等于单元体流体质量变化。

单元体中流体质量发生变化的速率 油 水

将达西定律代入连续性方程得出油、水的流动方程： 油 水

3、状态方程

为求解上述两个方程中的两个未知数P、Sw，引入状态方程： （1）

（2）h、 、 为常数 （3）

（4） 、 为 的函数 So=1-Sw

由连续性方程推导出油、水的流动方程，从而引入压力方程：

饱和度方程引入含水百分数：

最后把两维两相流动方程写成如下形式：

解上述方程通过把微分方程离散化变代数形式来求解。

三、数值模拟的目的 （一）、为什么开展油藏数值模拟工作

研究和开发一个油田是一个复杂的综合性的科技问题，高精度的地震资料的处理

3

解释提供研究区域的构造、断层、边界及其走向，但地震纵向分辨率受到限制，不能很好的反映一个同相轴（地震道）中沉积砂体的物性变化特征；测井可较好的反映到小于1米以下沉积砂体的物性特征，提供可靠的地层对比结果。但作为新老油田开发方案的研究及剩余油分布的研究，是地震、地质、测井理论方法都无法做到的。地质上仅定性或半定量分析，测井用于生产监测不能以点带面。惟独油藏数值模拟工作可再现生产历史，定量分析剩余油潜力；并做到室内研究投入少、时间短，还可进行开发方案优选及经济评价工作。所以总公司强调开发方案的部署一定要开展数值模拟工作。值得强调的是油藏数值模拟工作提倡一体化，注重前期的地震解释和测井解释即油藏描述工作。 （二）、油藏数值模拟的目的 在进行油藏数值模拟工作前，首先应根据油田开发过程中存在难以解决的实际问题，提出开展此项工作的目的及意义，即最终所要达到解决问题的目标是什么？一般通过油藏数值模拟可进行以下研究工作： 1. 初期开发方案的模拟

1). 评价开发方式；如：枯竭开采、注水开发等。 2).选择合理井网、开发层系、确定井位； 3).选择合理的注采方式、注采比； 4).对油藏和流体性质敏感性研究。 2.对已开发油田历史模拟

1).核实地质储量，确定基本的驱替机理（如：是天然驱，还是注水开发。）； 2).确定产液量和生产周期； 3).确定油藏和流体特性； 4).提出问题、潜力所在区域。 3.动态预测

1).开发指标预测及经济评价

2).评价提高采收率的方法（如：一次采油、注水、注气、化学驱等） 3).剩余油饱和度分布规律的研究,再现生产历史动态 诸如：? 研究剩余油饱和度分布范围和类型；

? 单井调整：改变液流方向、注采井别、注水层位； ? 扩大水驱油效率和波及系数； 4）.潜力评价和提高采收率的方向 诸如：? 确定井位、加密井的位置；

? 确定油田开发最大产液量、产量对采收率的影响； ? 确定地面和井的设备。 5）．专题和机理问题的研究

诸如：? 对比注水、注气和天然枯竭开采动态； ? 研究各种注水方式的效果；

? 研究井距、井网对油藏动态的影响； ? 研究不同开发层系对油藏动态的影响； ? 研究注水速度对产油量和采收率的影响；

? 研究油藏平面性质和层间非均质性对油藏动态的影响； ? 验证油藏的面积和地质储量； ? 校验油藏数据；

? 为谈判和开发提供必要的数据。

4

注意

无论是对油藏进行初期开发方案、已开发油田历史模拟，还是动态预测的数值模拟工作，都要求油藏工程师要有针对性的拟定出能解决油田开发实际问题的数值模拟工作详细计划，及其开展此项工作的目的和应达到的目标是什么？

四、选择适当的数值模型及相类

对一般油藏而言，可具有两相（油＋水、油＋气、气＋水）和三相（油、气、水）的模型，而油藏维数可具有一维平面、垂直模型、两维平面（剖面）、两维径相（锥进）模型和三维平面等模型。但实际油藏其地下流体渗流机理、岩石及流体性质等地质特征不同；生产过程中的开采方式、机理不一等复杂问题，促使非常规油藏模型孕育而生。如：

?天然裂缝油气藏的模拟应用双重介质的裂缝模型； ?凝析油气田开发以及注气油田开采应用组份模型；

?视油、气田含富气或平气大小可采用组份或黑油模型； ?稠油开发用热采注蒸汽模型；

?注各种化学剂的三次采油应用化学模型；

?而对一般油藏，模拟可选用常规油气田开发的黑油模型。 注意

如何确定多相问题？需了解油藏流体性质、生产特征，以及采油机理。如： ?气藏是否带油环和下伏水区；

?有无明显相态变化的混相驱油藏； ?有无表现出凝析现象的相态特征；

?气顶、及水锥效应是否占优势的油藏等；

由此，在模拟区已确定的同时，需要认真考虑选择模型及相类的问题。 五、数据录取的准备工作

在数据录取工作中，建立精细的油藏地质模型是模拟工作成败的关键；而对生产数据的正确分析、合理取舍、根据生产情况合理分层和划分时间阶段步长，即录入生产数据卡工作也极为重要。 油藏模拟工作要研究油田开发中的问题很多，不可能一律对待，也不可能在一个模型里都予以满足，需根据研究任务和客观许可条件建立相应的地质模型。 （一）建立油藏地质模型

1.在油藏描述的基础上建立油藏地质模型

开展油藏描述工作，对油藏的地质、油层非均质特征，沉积相的详细描述和研究，根据油藏沉积相研究建立该油藏特征的沉积模式。 油藏描述分析的目的是综合所有的测井、岩心和生产测试等资料来得出一个与全油田一致的储集层模型。对各种未知的基本参数例如：对顶面深度、砂厚、孔、渗、饱等空间分布的评价中最大限度地发挥现有测井资料的作用，同时将这些参数结合所需储集层的几何特性参量进行计算，并结合地质沉积相分析提供出更为精细、完善的油田地质模型。

确定一个油藏地质模型所需的许多参数，在油藏勘探试采阶段或初期刚投入开发阶段用有限几口井的资料进行计算、解释及建模其精度是不高的，这些参数初始误差越大，则通过历史拟合达到令人满意的油藏特征描述所需的时间越长。所以，无论在对老区或新区进行数值模拟时，应对所选区块选用所有的井（特别是―关键井‖）开展油藏描述工作，进行全面的分析研究。

5

φ----孔隙度

2、如果油井测得压力数据是一个压力点，该点位于压力恢复曲线直线段上， 则：首先将质量总流率（产量）Q换算成体积总流率q ,

1/r 0)+1/r 0]?又: q = Q [fw(1 ,然后用下式计算出校正的实际地层压力。

式中：P0----校正的实际地层压力 Pwsd----折算到Dd处的压力 q-----关井前稳定总产率 ----总流度

H ----射开厚度 -----孔隙度 ------关井时间

、 ------网格x、y方向的步长

附件2：

关于烃类有效孔隙体积的计算 在渗流过程中，地层的孔隙体积实际上是发生变化的，它与地层压力和压实作用有关：

式中： ----压缩因子 ----节点当前平均压力 ----原始孔隙体积 ----岩石压缩系数 ----节点当前压力 ----基准面压力

----变化后的孔隙体积

当 时， ，说明孔隙体积没有变化，在简化三维三相黑油模型中，也常常设 ，即岩石是不可压缩的，但这对模拟结果可能带来某种影响，在空间某一微小的单元里，虽然孔隙体积变化很微弱，但它可能影响到烃类有效孔隙体积的大小，也就是说，对剩余油的分布有一定影响。因此，需进行烃类有效孔隙体积的计算。 为什么？ 通常情况下，数值模拟提供各时间阶段的剩余油分布来，并不意味着剩余油饱和度大的地方，就是挖潜或调整的重点，因为剩余油饱和度与油层厚度没有直接关系，因此，其分布状况往往掩盖了厚油层与薄油层的区别，为了使剩余油分布从定性向定量转化，可用烃类（包括油和水）有效孔隙体积PV来表示： ?（Sor ?Bo?r0/r0r+Sgr?Bg? rg/rgr）?PV = DX ? DY ? H ? 其中：DX 、 DY –---为网格节点X、Y方向步长， H---油层有效厚度，

Sor 、Sgr ---为剩余油、气饱和度， Bo、Bg—为油、气体积系数，

r0、r0r和rg、rgr---为油、气地面和地下密度。

参考书：

1． 油藏数值模拟----------范江，石油工业出版社

2． 油藏数值模拟原理--------李允、杜志敏、钟义贵，成都科技大学出版社

11

3． 油藏数值模拟技术应用成果集------中国石油天然气总公司开发生产局， 石油工业出版社

4． 黑油和组分模型的应用------------李福垲，科技出版社 5．油气藏多学科综合研究-------张一伟 刘洛 欧阳健 傅全全 王伟华，严俊君，石油工业出版社

12

本文来源：https://www.bwwdw.com/article/jhbf.html