分布式并行油藏模拟高效求解器的构建ESTABLISHMENT OF THE HIGH-EFFICIENCY SOLVER FOR DISTRIBUTED PARALLEL RESERVOIR SIMULATION
赵国忠
ZHAO Guozhong
摘要(Abstract):
在分析、借鉴近年来油藏模拟线性系统求解器研究进展与数值实验的基础上,考虑自有并行模拟器的架构特点,以简洁实用为目的采取两步预处理方法,构建了新的高效求解器。在第一步预处理过程中,把构成雅可比矩阵的累积项和流动项分开考虑,提出了不同于以往的解耦矩阵;在第二步预处理方法的选择上,兼顾了分布式并行的实现便利和处理效率,提出了采用z方向线红黑排序的块高斯-赛德尔法(ZLBGS)迭代求解方案。SPE标准算例和油藏模拟实例测试都显示,新的并行求解器计算效率高于采用ASM+BILU单步预处理的原有求解器,提升了自主并行模拟器的性能。
On the bases of analyzing and referring the recent research advances and numerical experiments of the solver in the linear system of the oil reservoir simulation,considering the infrastructural features of the in-house parallel simulator,a new high-efficiency solver was established for the goal of simplicity and practicality by the use of two-stage preconditioning method. In the process of the first step,the decoupling matrix different from ones accepted widely in the past was proposed by splitting the accumulative item and flux parts of Jacobi matrix; on the choice of the second step preconditioning method,taking the convince and efficiency into account for the distributed parallel simulation,Z-line red / black-ordering block Gauss-Seidel( ZLBGS) iterating method is proposed. The test results from SPE standard modeling and actual field simulation both demonstrate that the computing efficiency by the new parallel solver is higher than the previous one by means of single-step ASM + BILU preconditioner,thus the performances of the in-house parallel simulator have been promoted.
关键词(KeyWords):
油藏模拟;求解器;分布式并行
reservoir simulation;solver;distributed parallel
基金项目(Foundation): 中国石油天然气集团公司重大科技专项“大庆油田原油4000万吨持续稳产关键技术研究”(2011E-12)
作者(Author):
赵国忠
ZHAO Guozhong
参考文献(References):
- [1]Watt J W.A Method for Improving Line Successive Overrelaxation in Anisotropic Problems-A Theoretical Analysis[C]//Trans.Soc.Pet.Engr.,1973.
- [2]Vinsome P K W.ORTHOMIN,an Iterative Method for Solving Sparse Banded:Sets of Simultaneous Linear Equations[R].SPE5729,1976.
- [3]Wallis J R,Kendall R P,Little T E.Constrained Pressure Residual Acceleration of Conjugate Residual Methods[R].SPE 13536,1985.
- [4]Cao H,Tchelepi H,Wallis J,et al.Parallel Scalable Unstructured CPR-Type Linear Solver for Reservoir Simulation[R].SPE96809,2005.
- [5]Stiiben K,Clees T,Schneider M,et al.Algebraic Multigrid Methods(AMG)for the Efficient Solution of Fully Implicit Formulations in Reservoir Simulation[R].SPE 105832,2007.
- [6]Fung L S,Dogru A H.Parallel Unstructured Solver Methods for Simulation of Complex Giant Reservoir[R].SPE 106237,2007.
- [7]Al-Shaalan T M,Klie H,Dogru A.H,et al.Studies of Robust Two Stage Preconditioners for the Solution of Fully Implicit Multiphase Flow Problems[R].SPE 118722,2009.
- [8]Cusini M,Lukyanov A A,Natvig J,et al.A Constrained Pressure Residual Multiscale(CPR-MS)Compositional Solver[C]//ECMOR XIV-14th European Conference on the Mathematics of Oil Recovery,Catania,Sicily,Italy,2014.
- [9]Christie M A,Blunt M J.Tenth SPE Comparative Solution Project:A Comparison of Upscaling Techniques[R].SPE 72469,2001.
- [10]赵国忠,李保树.三维三相油藏模拟新解法[J].石油学报,1995,16(4):68-74.
- [11]赵国忠,尹芝林,吴邕.大庆油田PC集群大规模油藏数值模拟[J].西南石油学院学报,2003,25(6):35-39.
- [12]赵国忠,孙巍,何鑫.基于分层注水数学模型的油藏模拟[J].东北石油大学学报,2012,36(6):82-87.