针对三相黑油模型模拟时大规模对角稀疏线性方程组存储占用大量内存、求解耗时的问题,研究大规模稀疏系数矩阵的压缩存储及求解方法.提出了有效节点压缩和块压缩存储结合的压缩存储方法,以油藏节点为单元,利用3个实数组分别存储稀疏矩阵块主对角元素、下三角及上三角的非零块元素；由于模拟中形成的系数矩阵具有位置对称性,利用3个整型数组存储下三角非零块元素的地址信息,减少了内存损耗及运算中非零元素的搜索次数.在此基... 展开>>针对三相黑油模型模拟时大规模对角稀疏线性方程组存储占用大量内存、求解耗时的问题,研究大规模稀疏系数矩阵的压缩存储及求解方法.提出了有效节点压缩和块压缩存储结合的压缩存储方法,以油藏节点为单元,利用3个实数组分别存储稀疏矩阵块主对角元素、下三角及上三角的非零块元素；由于模拟中形成的系数矩阵具有位置对称性,利用3个整型数组存储下三角非零块元素的地址信息,减少了内存损耗及运算中非零元素的搜索次数.在此基础上,采用Krylov子空间块广义极小余量GMRES迭代法结合块不完全LU预处理方法求解上述线性方程组.结果表明,解法收敛速度快、稳定性好,可以提高大型稀疏线性方程组的求解速度.经过黑油模型标准考题的测试证实,该方法能够有效求解大规模油藏数值模拟问题. 收起<<

王宝华  吴淑红  韩大匡  桓冠仁  李巧云  李小波  李华  周久宁 

中国石油勘探开发研究院;提高石油采收率国家重点实验室/中国石油大学 北京

石油勘探与开发 ISTIC EI PKU 2013年4期

数值模拟  黑油模型  有效节点压缩  块压缩存储  块ILU分解  块广义极小残量法  reservoir simulation  black oil model  active nodes compression  block compressed storage  block ILU factorization  block generalized minimal residual method 

北京市"陆相沉积油藏数值模拟系统研制" 中国石油天然气股份有限公司"新一代油藏数值模拟软件" 中国石油勘探开发研究院"油藏数值模拟GPU并行求解技术"