黑油模型半隐式全隐式和IMPIMS方法.PPTVIP

第六章 黑油模型（半隐式、全隐式和IMPIMS方法） 半隐式方法 全隐式方法 隐式压力隐式饱和度方法 第一节 半隐式方法 一、数学模型 1. 假设条件 1) 符合达西渗流定律 2) 等温渗流 3) 油气水三相和油气水三组分，油气相之间有质量交换，即气组分不仅存在于气相中，而且存在于油相；但油相和水相之间没有质量交换 4) 油相和气相随压力变化而发生相态变化 5) 岩石和流体均可压缩 6) 油藏非均质和各向异性 7) 考虑毛管力和重力 2. 组分质量守恒方程 油组分 气组分 水组分 3. 辅助方程 1) 2) 3) 4. 初始条件 5. 边界条件 外边界 1) 定压外边界 二、差分方程组的建立 先不讨论源、汇项，方程（1）、（2）、（3）的差分方程可写为 油组分 气组分 水组分 （7）、（8）、（9）的未知数为Po、Sw、Sg 1．考虑流动系数T的半隐式油组分差分方程（7）的左端项 将（10）式代入（7）式的左端项 上式中， 为已知，变量为 忽略二阶小量得 ２.考虑毛管压力项Ｐc和流动系数项T的半隐式气组分和水组分方程的左端项 气组分方程（8）的左端项 代入方程（8）的左端项后，略去二阶小量后得： 自由气 式中 溶解气 式中 在水组分方程（９）的左端项 代入方程（8）的左端项，略去二阶小量后得： 式中 3. 方程（7）、（8）、（9）的右端项 油组分方程（7） 式中 气组分方程（8） 式中 水组分方程（9） 式中 4. 由方程右端项和左端项得： 油组分方程 （11）＝（15） 气组分方程 （12）＋（13）＝（16） 水组分方程 （14）＝（17） 油组分方程（18）可写为： 气组分方程（19）可写为： 水组分方程（20）可写为： 5. 一维情况下，（21）～（23）式的展开式 油组分 水组分 对 按上游权的原则，如： 其中各导数可在上一时间步的值附近按割线斜率法计算，如： 式中dP是预先给定P的一个小的增量。 将（24）、（25）、（26）按三个变量排列得： 代表油，水，气组分方程 b,c,d 代表前一节点，本节点，后一节点的系数 1,2,3 代表 的系数 写成矩阵方程 式中 A为三对角块状矩阵 i=1,2,3……N bi,ci,di为3×3块子矩阵。 三对角块状矩阵可用追赶法同时解出 。 二维情况下为五对角块状矩阵方程； 三维情况下为七对角块状矩阵方程； 两相情况下为2×2块子矩阵。 6. 考虑产量项的半隐式处理 1) 定产油量 将产油量直接代入油组分的差分方程组中，其它各相的产量可按流动系数计算。 令 则 2) 定总液量 令各相流度比 则 3) 定井底流压 第二节 全隐式方法 一、数学模型 1. 假设条件 1) 符合达西渗流定律 2) 等温渗流 3) 油气水三相和油气水三组分流动，油气相之间有质量交换，即气组分不仅存在于气相中，而且存在于油相；但油相和水相之间没有质量交换。 4) 油相和气相随压力变化而发生相态变化。 5) 岩石和流体均可压缩 6) 油藏非均质和各向异性 7) 考虑毛管力和重力 2. 组分质量守恒方程 油组分 气组分 水组分 3. 辅助方程 1) 2) 3) 4. 初始条件 5. 边界条件 外边界 1) 定压外边界 二、差分方程组的建立 1. 预备知识—牛顿迭代法 1) 在任一个从tn到tn+1的时间步长内，设第l+1和l两次迭代结果未知函数x之差为 即 式中x即为P、Sw、Sg。 2) 令l=0时 x0=xn 求解方程组内流动系数和毛管压力的值，然后解方程组得出第一次迭代值x1，然后再求解流动系数和毛管压力的值，解方程组得第二次迭代值x2，依次类推。 3) 当