《油藏数值模拟基础培训(第六讲)剖析》由会员分享，可在线阅读，更多相关《油藏数值模拟基础培训(第六讲)剖析（25页珍藏版）》请在金锄头文库上搜索。
1、油藏数值模拟基础培训 （第六讲） 张世明，2009年4月，职工大学 油藏数值模拟基础培训 求实 创新 超越 REALITY，INNOVATION,TRANSCENDENCY REGION/SOLUTION要点分析 REGION部分要点分析 求实 创新 超越 REALITY，INNOVATION,TRANSCENDENCY 油藏不同部位可能有不同的流体属性，比如不同断块的油密度，粘度不 同，或你的油藏岩性在纵向或垂向有变化，那你就需要在你的模型中设置 流体或岩性分区。 分区设置REGION： 储量分区FIPNUM：如果你想输出模型不同部位的储量，你需要设置储量分区 流体分区PVTNUM：如果你的
2、模型不同部位流体PVT属性不同，你需要设置流体分 区。 岩性分区SATNUM：如果你的模型不同部位岩性不同，需要用不同的相渗曲线和 毛管压力曲线，你需要设置岩性分区。 平衡区分区EQLNUM：如果你的模型有不同的油水或油气界面，你需要设平衡区 分区。 REGION部分要点分析 求实 创新 超越 REALITY，INNOVATION,TRANSCENDENCY 分区设置REGION： 我们假设你的模型有东西两个断块，两断块被封闭断层分割。东断块的油比 西断块的油密度重，在你提供油藏流体PVT表时你需要提供两个表，一个是密度 重的PVT表，另一个是密度轻的PVT表。在你的分区部分将东断块的流体分区
3、值设 为1，将西断块流体分区值设为2。软件在计算东断块的流体流动时将自动用第一 个PVT表（密度重的PVT表），在计算西断块流体流动时用第二个PVT表（密度轻 的PVT表）。 REGION部分要点分析 求实 创新 超越 REALITY，INNOVATION,TRANSCENDENCY 分区设置REGION： REGION部分关键字 REGION部分要点分析 求实 创新 超越 REALITY，INNOVATION,TRANSCENDENCY 储量分区： REGION部分 FIPNUM 1 2 3 / 流体分区： PROPS部分 PVTO table1 / table2 / table3 / PV
4、DG table1 / table2 / table3 / PVTW table1 / table2 / table3 / PROPS部分 PVTW table1 / table2 / table3 / DENSITY table1 / table2 / table3 / ROCK table1 / table2 / table3 / REGION部分 PVTNUM 1 2 3 / REGION部分要点分析 求实 创新 超越 REALITY，INNOVATION,TRANSCENDENCY 岩性分区： PROPS部分 SWOF table1 / table2 / table3 / SGOF t
5、able1 / table2 / table3 / REGION部分 SATNUM 1 2 3 / 平衡区分区： SOLUTION部分 EQUIL line1 / line2 / line3 / RSVD table1 / table2 / table3 / REGION部分 EQLNUM 1 2 3 / REGION部分要点分析 求实 创新 超越 REALITY，INNOVATION,TRANSCENDENCY 分区设置REGION： 同一个流体区内可以有不同的平衡区，但用一个平衡区 内不容许存在不同的流体区。其余的分区可以交叉。 SOLUTION部分要点分析 求实 创新 超越 REALIT
6、Y，INNOVATION,TRANSCENDENCY （1）计算过渡带高度。由油水界面和油气界面深度以及相渗曲线提供的最大毛管压 力计算。 （2）计算每一个网格初始的油相，水相，气相压力分布。首先将在流体属性部分提 供的油，气，水地面密度折算为地下密度。基于参考点的深度和对应压力以及油水 界面，油气界面深度，过渡带高度，结合油，气，水地下密度计算其他深度处的油 ，气，水相压力。 （3) 由每个网格的油，气，水压力计算油水和油气毛管压力 （4) 计算饱和度分布。这部分计算主要用你提供的相渗曲线端点值。将油水界面以 下的含水饱和度设为你在油水相渗曲线中提供的最大含水饱和度，通常为1。将油气 界面以
7、上的含气饱和度设为你提供的油气相渗曲线的最大值。油气界面以上的含水 饱和度为束缚水饱和度。在油区的含油饱和度为1减束缚水饱和度。在过渡带的含油 和含水饱和度由你提供的毛管压力曲线得到。 （5）计算初始溶解油气比或泡点压力的分布。 模型初始化过程： SOLUTION部分要点分析 求实 创新 超越 REALITY，INNOVATION,TRANSCENDENCY 模型初始化过程： 平衡初始化 根据已知的油藏压力系数、油水及油气界面、界面处的毛管压力和平 衡系统的相对渗透率曲线等，考虑重力及毛管压力影响，建立油藏原始 压力及饱和度分布。 油区： SWI 气顶： SWI 油气过渡带： PCOG + S
8、WI 油水过渡带： PCOW + SWI 水区： 100%SW SOLUTION部分要点分析 求实 创新 超越 REALITY，INNOVATION,TRANSCENDENCY SOLUTION部分要点分析 求实 创新 超越 REALITY，INNOVATION,TRANSCENDENCY SOLUTION部分要点分析 求实 创新 超越 REALITY，INNOVATION,TRANSCENDENCY 各网格的初始压力和相初始饱和度； 各网格块的初始溶解气油比和（或）溶解油气比； 油藏流体属性对深度的依赖关系，包括API、饱和溶解气油比、 泡点压力、饱和溶解油气比OGR和露点压力随深度的变化；
9、 解析水体的初始状况。 平衡初始化 SOLUTION部分要点分析 求实 创新 超越 REALITY，INNOVATION,TRANSCENDENCY 平衡法指定界面的深度、一个参考深度及此参考深度处的压力； 假设满足流体静力学平衡； 平衡法适用于有水平界面的初始未钻开油藏； 这里可能需要输入不同的毛管力； 用黑油的状态方程来计算各处各相的静压力； 从PROPS部分的饱和度方程中得到各区域的相饱和度； 过渡区中的Sw和Sg通过反查饱和度函数中的毛管力曲线得到。 平衡初始化 SOLUTION部分要点分析 求实 创新 超越 REALITY，INNOVATION,TRANSCENDENCY 平衡初始化
10、 平衡法计算初始压力和饱和度是基于饱和度函数、流体界面、参考深 度和参考压力的。如果已知的是处于油区的一个参考深度处的压力，则 会循环求解油相的黑油状态方程和油相的静压力方程。 黑油状态方程： 油相静压力方程： 通过以上方程得到各处的油相压力。其中（r）和（s）分别代表油藏 和地面条件，h代表深度。给定油水界面，用同样的方法求解气相和水相 的状态方程，最终可以得到油藏中各个位置处的各相静压力。 SOLUTION部分要点分析 求实 创新 超越 REALITY，INNOVATION,TRANSCENDENCY 相初始饱和度 计算出初始压力之后，ECLIPSE会为每一个区域赋予相饱和度值。 但这不包
11、括过渡区。相饱和度是从PROPS部分定义的饱和度函数的饱 和度端点值中得到，而且油相饱和度总是通过水相和气相饱和度计算得 到的。其中： 在气区，Sg为最大气相饱和度SGU，是输入的气相饱和度函数表中 的最大含气饱和度值。Sw是最小值SWL，是输入的水相饱和度函数表 中的最小含水饱和度值，即原生水或束缚水饱和度Swco。而油相饱和 度为So=1-SGU-SWL。 SOLUTION部分要点分析 求实 创新 超越 REALITY，INNOVATION,TRANSCENDENCY 相初始饱和度 在水区，Sw为最大水相饱和度SWU，是输入的水相饱和度函数表 中的最大含水饱和度值，此值常取值为1。Sg为最
12、小值SGL，是输入的 气相饱和度函数表中的最小含气饱和度值，即原生气饱和度Sgco，此 值常取为0。而油相饱和度为So=1-SWU-SGL。 在油区，含水和含气饱和度都处于最小值，Sg=SGL=Sgco， Sw=SWL=Swco。则含油饱和度为So=1-SWL-SGL。通常，So=1- Swco。 SOLUTION部分要点分析 求实 创新 超越 REALITY，INNOVATION,TRANSCENDENCY 过渡区内的相初始饱和度 在过渡区，相的饱和度是由毛管力控制的。毛管力的定义为： Pcow=Po-Pw和Pcog=Pg-Po。毛管力可以通过计算任意深度处相 间的静压力差得到。然后，ECL
13、IPSE会反查输入的饱和度函数中的 毛管压力曲线来确定过渡区的含水和含气饱和度。理论上，在Sw=0 时，水相的毛管力曲线Pcow趋于无穷大。而实际上，毛管力曲线应 该在Swco处终止，以保证不会在饱和度函数的最低含水饱和度处出 现不连续的情况。 让我们考虑一个穿越油水界面OWC的倾斜网格。通过计算网格的 中心深度，及前面的方法所计算出来的初始饱和度Soi、Swi、Sgi， 就可以得到此网格块的饱和度。这种方法提供了一种完全稳定的初 始化方法，其相饱和度与静压力是一致的。这种方法为块中心平衡 。 SOLUTION部分要点分析 求实 创新 超越 REALITY，INNOVATION,TRANSCE
14、NDENCY SOLUTION = EQUIL - -Depth Pi WOC Pcow|woc GOR Pcog|Goc RSVD RVVD ACC 1500 150.0 1580 0.0 1500 0.0 1* 1* 10 SOLUTION = EQUIL - -Depth Pi WOC Pcow|woc GOR Pcog|Goc RSVD RVVD ACC 1500 150.0 1580 0.0 0.0 0.0 1 1* 10 RSVD 2000.0 96.11 3000.0 96.11 / SOLUTION部分要点分析 求实 创新 超越 REALITY，INNOVATION,TRANS
15、CENDENCY 关于EQUIL关键字 在单相、油水和气水两相模型中忽略GOC； 在单相和油气模型中忽略OWC； 在对初始高于露点压力的凝析气藏的模拟运算中，给定的油水界面OWC 与油气界面GOC的深度应该相同； 如果是仅含有挥发油的自由气，则油气GOC界面可以低于储层底部； 如果初始时没有可动水，则油水界面WOC可以低于储层底部； 第9项对精度的选择默认为-5，N0分别对应水平块和倾斜块的初 始化平衡计算。 不同油藏对EQUIL关键字的参数项的解释是不同的。 SOLUTION部分要点分析 求实 创新 超越 REALITY，INNOVATION,TRANSCENDENCY 重启计算 历史拟合结束后需要进行产量预测，在进行产量预测计算时，不需要 再从历史拟合开始时进行计算，可以直接从历史拟合结束的时间接着往下 算。这种应用上一次计算的输出作为下一次计算的初始输入计算就叫重启 计算。 要进行重启计算，首先要定义重启时间步的输出。可以用RPTRST来 定义输出每时间步，每月，每年或每隔几月几年重启时间步文件。如果采 用多文件格式输出，则文件后缀为：.X0000, .X0001等，如果是单文件输 出，则输出文件为.UNRST. 重启文件记录了每时间步模型压力分布，饱和度分布，溶解油气比分 布，同时也记录所有井的井位，射孔位置，产量控制。不
