针对油藏压力方程混合存储结构的多重网格求解方法技术

本发明专利技术提供一种针对油藏压力方程混合存储结构的多重网格求解方法，该针对油藏压力方程混合存储结构的多重网格求解方法包括：步骤1，计算得到最细网格层上的初始解向量和相应残向量；步骤2，将细网格层上的残向量限制到粗网格层上；步骤3，计算得到相应粗网格上的初始解向量残向量；步骤4，将粗网格层上的解向量提升到相应细网格层上；步骤5，计算得到最细网格层上的解向量和相应残向量完成一次完整的AMG‑V‑Cycle循环。该针对油藏压力方程混合存储结构的多重网格求解方法既保持了多重网格法的普适性和收敛速度快等特性，又充分发挥了结构矩阵的存储与运算优势，并通过数值实验验证新方法的正确性和高效性。

【技术实现步骤摘要】

本专利技术涉及油藏数值模拟算法研究、技术应用，特别是涉及到一种针对油藏压力方程混合存储结构的多重网格求解方法。

技术介绍

1、油藏数值模拟研究最早可追溯到上世纪50年代，在70年代后期逐渐发展，至今已建立各种油藏数学模型，如黑油模型、组分模型、热采模型和化学驱模型等，其中黑油模型是描述油藏模拟问题的基础模型，化学驱模型通常是在黑油模型的基础上扩展得到。

2、压力方程求解是预条件过程的关键环节，主要选取amg来进行求解。amg法是目前国际上求解大规模pdes离散代数系统的极为有效的方法，于上世纪八十年代初由brandt等人提出，近年来得到了极大的发展和应用。目前出现了各种各样的amg法，如基于强弱连通的rs-amg法，基于聚集法的amg法(agamg法)，自适应amg法，基于falgout粗化算法、hmis粗化算法和长距离插值算法的amg法等。其中，由于agamg法具有算法复杂度低、更适合于异构体系等特点，目前已成为amg法研究领域的热点。这些agamg法对某些简单标量型偏微分方程离散代数系统具有较好的收敛效果。然而对于复杂pdes离散系统，如化本文档来自技高网...

【技术保护点】

1.针对油藏压力方程混合存储结构的多重网格求解方法，其特征在于，该针对油藏压力方程混合存储结构的多重网格求解方法包括：

2.根据权利要求1所述的针对油藏压力方程混合存储结构的多重网格求解方法，其特征在于，该针对油藏压力方程混合存储结构的多重网格求解方法还包括，在步骤1之前，针对化学驱油藏模型，基于三维结构网格采用7点中心差分格式进行离散得到代数系统，其中油藏矩阵由C、E、W、S、N、T和B这7条带构成，C表示单元中心，E、W分别代表其东西方向，S、N分别代表其南北方向，T、B分别代表其上下方向。

3.根据权利要求2所述的针对油藏压力方程混合存储结构的多重网格求解方...

【技术特征摘要】

1.针对油藏压力方程混合存储结构的多重网格求解方法，其特征在于，该针对油藏压力方程混合存储结构的多重网格求解方法包括：

2.根据权利要求1所述的针对油藏压力方程混合存储结构的多重网格求解方法，其特征在于，该针对油藏压力方程混合存储结构的多重网格求解方法还包括，在步骤1之前，针对化学驱油藏模型，基于三维结构网格采用7点中心差分格式进行离散得到代数系统，其中油藏矩阵由c、e、w、s、n、t和b这7条带构成，c表示单元中心，e、w分别代表其东西方向，s、n分别代表其南北方向，t、b分别代表其上下方向。

3.根据权利要求2所述的针对油藏压力方程混合存储结构的多重网格求解方法，其特征在于，在步骤1，对最细网格即第0层系数矩阵采用结构矩阵str方式进行存储，并使用基于str方式的高斯赛德尔或雅克比磨光运算和矩阵向量乘积线性运算得到最细网格层上的初始解向量和相应残向量。

4.根据权利要求3所述的针对油藏压力方程混合存储结构的多重网格求解方法，其特征在于，在步骤1，由于油藏矩阵具有7条带状结构，隐式井对应的镶边矩阵具有对称性和规模小的特点，针对以上两种矩阵分别使用结构矩阵和二维数组进行存储。

5.根据权利要求4所述的针对油藏压力方程混合存储结构的多重网格求解方法，其特征在于，在步骤1，采用n表示有效单元数、nwell表井的总个数、nxy＝nx*ny，这里nx、ny分别为x、y方向的单元个数；o表示非零元，m为井的最大射孔数；为保证能够形成油藏与井镶边矩阵的稀疏结构，需要对二维数...

【专利技术属性】
技术研发人员：曹伟东，于金彪，曹绪龙，曹小朋，段敏，孟薇，马泮光，张波，史敬华，郑乃元，
申请(专利权)人：中国石油化工股份有限公司，
类型：发明
国别省市：