热采数值模拟中的区域分解算法.pdfVIP

维普资讯 石 油 勘 探 与 开 发 110 2002年 4月 PE1、ROLEUM EXPLoRA1、IoN AND I3EVELOPMENT Vo}．29 No．2 文章编号：1000—0747(2002)02一叭1O—o3 热采数值模拟中的区域分解算法 吴淑红，刘翔鹗，郭尚平 【中国石油勘探开发研究院) 摘要 ：区域分解算法在浩藏数值模拟 中的应用 已得翌 了国珂外学者的共识 ，然而 ，谊方法在热采数值模型中的应1田极少 根据热采模型的多相 、多组舟 的特点，构造 了适奇垂晦式计算方法的区域分解算法，并给 出详鳍的实现步骤和噩域收敛 的 判断条件。实例计算对 比表 明：区域舟解算法适用干热采模型．带书区域舟解算法的热采模型与常规热呆模 型的计算结 果吻合很好 ，计算结果是可信的；使用区域分解算法可进一步提高热秉模型的计算运度 ．节省计算时间。圈4参 4(吴淑生 摘) 关键词 ：区域分辑算法；热采模型 ；奎 隐式 中圈分娄号 ：TE35744 文献标识码 ：A 区域分解算法在数学领域研究较早，起源于十九 由 I=l，2，…，l2的网格组成 ，记为Q1= U B ，Q2由 世纪七十年代。然而 ，将 区域分解算法用于油藏数值 I=4，5，…，l5的网格组成 ，记为 = U B ，r1是 Q1 模拟 中则 为时较 晚。二 十世 纪九十 年代 初 ，国外 的离散边界 ，rl= U B，， 是 的离散边界，r2= Gosselino1【和NacaulNCl2等作 了一些试验性工作。 U B 。(1)式解法如下 国内章寒松_3研究 了区域分散算法在 IMPES方 法的 黑油模型 中的应用 。 本文在热采数值模型基础上 ，根据全 隐式求解方 法，构造适合全隐式热采数值模型的区域分解算法，并 进行区域 收敛判断等研究。 1热采模型的区域分解算法 11求解 问题描述 在 区域Q上 ，某一时间步内，求解如下 问题： 喜『(等 h= (嘻SiXu) l (鲁 ) (2) I [嘻PJUJ+(1 T )](1) 1 Il∑ = 0 【∑ =0 {：n (3) 将上述 问题简记为：Au=f，其 中列 向量 “=(SI，S2， “i = “ (I12) S3，X1，X2，P，T) ，A 为方程组 的系数，f为右端项 。 1．2区域分解算法 ③重复上述过程 ，直至区域收敛。 将区域 n分解为两个子域 n 和 (见 图 1)。QI 由上述步骤可 以看出，在求解非线性方程组时 ，使 维普资讯 2002年4月 吴趣红 等：热栗数值摸担中的区域分解算法 用区域分解算法实质有双重迭代过程 (见围2)。首先 求解线性方程 (7)式，得到解 她 ； 是内迭代 (牛顿迭代)，求得非线性方程组的解；外迭代 判断牛顿迭代是否收敛，若不收敛 ，令 L+1 L，