第一章油藏数值模拟进展.PPT

第五章 黑油模型（IMPES方法） 数学模型 差分方程组的建立 井、初始场分布和过泡点处理 资料输入 第一节 数学模型 一 、假设条件 1. 符合达西渗流定律 2. 等温渗流 3. 油气水三相和油气水三组分，水组分存在于水相中，油组分存在于油相中，气组分不仅存在于气相中，而且存在于油相和水相中 4. 三维(x,y,z)方向流动 5) 岩石和流体均可压缩 6) 油藏非均质和各向异性 7) 考虑毛管力和重力 二、数学模型 1. 组分质量守恒方程 油组分 水组分 气组分 2. 辅助方程 3. 初始条件和边界条件 假设边界不规则的油藏中有若干口井生产或注入，求油藏中的压力和饱和度分布。 I.C 第二节 差分方程组的建立 一、方程（1）、（2）、（3）的右端项 根据（4）式，得 将（10）式代入（9）式得到只含So、Sw、Po的方程 二、组分质量守恒方程 将（7）、（8）、（11）代入（1）、（2）、（3） 的右端项得组分质量守恒方程 油组分 水组分 气组分 三、IMPES方法压力方程和饱和度方程 1. 思路 1) 乘以适当的系数，合并（12）（13）（14）以消除SW，SO得到只含变量PO，PW，Pg的压力方程。 2) 由毛管压力公式PCow=PO-PW， PCgo=Pg-Po，得到只含变量 的压力方程。 3) 达西系数项及毛管压力采用上一时间值，因此可得只含变量PO的线性 代数 方程组。 4) 解线性代数方程组后，求得 ，再求 5) 将 代入方程（1），求得 。 代入方程（2），求得 。 然后 2. 具体算法 1) 压力方程 (12)×(Bo-RsoBg)+(13)×(Bw-RswBg)+(14)×Bg 得 2) 隐式压力差分方程组的建立 对（15）式写差分方程时，由于网格节点多，因此可用简化形式。 式中 A— 传导系数， ， P —压力 对压力方程（15）进行隐式差分后，两端乘以 后，令 可得到以下隐 式压力差分方程组。 根据方程（16），对第i，j，k个网格写差分方程，由于其邻节点有六个，可形成七对角系数矩阵方程。若排列方式不同，则系数矩阵形式不同。 如标准排列，若有二层网格，先按k方向，后J方向，最后为I方向，顺序排列。 其系数矩阵如下： 如D4排列 对第i,j,k个网格，可写成如下的一般式： 3) 显式饱和度方程的建立 利用（1）式，经差分后，可得油饱和度方程 利用（2）式，经差分后，可得水饱和度方程 (18)、(19)式中的Pn+1可用隐式压力方程组(17)求得。 第三节 黑油模型特殊问题的处理 一、井的问题 它是内边界问题的处理 令 则 井的问题处理时，就Pwf的取法可分为两大类： 1. 显式井底压力 认为Pwf已知，即 。 2. 隐式井底压力 认为Pwf未知，即 。 以下讨论显式井底压力时多层情况下井的处理问题。 1. 定产量 1) 定总产油量Qo 设该井有n层，可利用各层的 将Qo分配到各层。 各层油产量 各层水产量 各层气产量 2) 定总液量Qt 首先计算各相的流度比 油相 水相 气相 则可计算油产量 再利用 （1）（2）（3）式可求得 Qok Qwk Qgk 3) 定注入井的总注水量Qwi或Qgi 设该井有n层，可利用各层的总流度 将Qwi或Qgi分配到各层。 2. 定井底流压 1) 隐式地层压力 若为生产井 PID＝PI 注入井 PID＝WI 将（7）式代入压力差分方程组（17）的Ql中去，由于将Ql变成了Pn+1和Pwf，方程（17）中的系数E， B将有变化 式中 当计算出 后，即可算出Ql，然后计算Qlk 2) 显式地层压力 计算出Ql 后，即可计算Qlk 二、初始压力场的分布 已知油气和油水接触面的压力和