新疆石油地质 ›› 2005, Vol. 26 ›› Issue (6): 695-697.

• 应用技术 • 上一篇    下一篇

黑油模型的预处理正交极小化解法

李中锋, 何顺利, 门成全   

  1. 中国石油大学 石油天然气工程学院,北京 102249
  • 收稿日期:2004-12-14 出版日期:2005-12-01 发布日期:2020-11-23
  • 作者简介:李中锋(1976-),男,河南遂平人,在读博士研究生,油藏工程,(Tel)010-89734268 (E-mail)lizhogfeng@sina.com.

MILU-ORTHOMIN Method for Black Oil Model

LI Zhong-feng, HE Shun-li, MEN Cheng-quan   

  1. College of Oil-Gas Engineering, China University of Petroleum, Beijing 102249, China
  • Received:2004-12-14 Online:2005-12-01 Published:2020-11-23

摘要: 数值模拟最终都要归结到大型线性方程组的求解,一个好的数值模拟软件应具备适用范围广、收敛速度快的优点。用 IMPES(隐压显饱)法求解黑油模型,引入井底压力为未知数时,压力方程的系数矩阵变成了加边的七对角阵。选择适合求解大型不对称稀疏矩阵的 MILU-ORTHOMIN(预处理正交极小化)方法,是预处理共轭梯度法的一种,是MILU(修正不完全因子分解)法和 ORTHOMIN(正交极小化)法的有机组合; 该方法中 MILU 法保持了原有稀疏矩阵的稀疏性,节省了存储空间,ORTHOMIN 法则大大提高了迭代收敛速度。根据该思想编制了水平井单井模拟器软件,使用IMPES 法并隐式处理水平井,取相同的迭代精度,在求解同一个问题时,其迭代收敛速度是 Gauss-Seidel(异步迭代)法的100多倍。

关键词: 数值模拟, 黑油模型, 方法, 水平井

Abstract: A numerical simulation at last comes down to calculation of a number of large linear equation groups, and an excellent simulator should be provided with advantages of wide application and rapid convergence.As implicit pressure and explicit saturation (IMPES) method is selected to calculate black oil model and the bottom-hole pressure is introduced, the coefficient matrix of pressure equation becomes bordered seven-symmetrical matrix.The modified incomplete factorization (MILU) -orthogonal minimization (ORTHOMIN) method can be used to calculate large asymmetric sparse matrix as one of the pretreatment conjugate gradient methods.MILU method keeps the sparse character of sparse matrix and saves the storage space, while ORTHOMIN method accelerates the rate of iteration convergence obviously.According to the MILU-ORTHOMIN method, the horizontal single well simulator is programmed.The IMPES and the implicit wellbore pressure method are employed in this paper.As the same iteration accuracy rate is selected to solve the same problem, it is found that the MILU-ORTHOMIN method makes the convergence speed increased by 100 times compared to that of normal Gauss-Seidel method (asynchronous iteration).

Key words: numerical simulation, black-oil model, method, horizontal well

中图分类号: