黑油模型解剖.docxVIP

 MACROBUTTON MTEditEquationSection2 Equation Chapter 1 Section 1 SEQ MTEqn \r \h \* MERGEFORMAT  SEQ MTSec \r 1 \h \* MERGEFORMAT  SEQ MTChap \r 1 \h \* MERGEFORMAT 1 黑油模型理论基础 1.1 基本假设 （1）油藏中渗流是等温渗流； （2）油藏中有油、气、水三相，各相流体的渗流均符合达西定律； （3）模型考虑油组分、气组分、水组分三组分； （4）气组分在油气相、水气相之间发生质量交换； （5）相平衡瞬间完成； （6）水组分只存在于水相中，与油气相之间没有质量交换； （7）油藏岩石微可压缩，各向异性； （8）油藏流体可压缩，且考虑渗流过程中重力、毛管力的影响。 1.2 数学模型  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 1) 辅助方程：  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 2) 初始条件：  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 3) 边界条件：  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 4) 2 黑油模型程序整体结构图 3 组员及分工 姓名负责程序张敏主程序，CGIN, MF1, INTERP, INTPVT高博主程序，TRAN1, SOLMAT, MATBAL, PRTPS郝延征主程序，TABLE, UINIT1, CODES, NODES徐杨主程序，GRID1, PARM1, QRATE4 主程序 4.1 主程序主要功能 （1）定义运算所需数组； （2）需要调用和生成的文件的打开和关闭； （3）通过调用子程序给模型赋基础数据和初始数据； （4）通过调用子程序给模型的运行确定各种控制； （5）在运行过程中反复读入井点数据（包括产量和注入量、井底流压、流动指数等）以及打印输出控制码； （6）分层计算油气水地质储量； （7）进行井点产量项处理和形成压力矩阵； （8）通过调用子程序求解压力方程； （9）显式计算饱和度； （10）进行过泡点处理； （11）根据压力和饱和度增量控制，自动调节计算时间步长； （12）未饱和网格块饱和度计算； （13）变泡点处理； （14）在每一运算时间步末进行物质平衡检验； （15）打印油井、水井、气井的分层报告和分井报告； （16）通过调用子程序进行计算结果打印输出和形成文件； （17）重启动运行方式的选择和重启动文件的生成和调用； （18）运算终止的控制及错误信息的提示。 4.2 主程序流程图 5 子程序CGIN 5.1 功能介绍 用不完全LU分解预处理共轭梯度法求解压力方程隐式差分得到的大型稀疏线性方程组。 5.2 理论基础 实践表明，预处理共轭梯度型方法无论是适应性还是计算速度，都远远超过了前面的方法。这类算法的优点是迭代收敛速度不依赖于迭代因子的选取，收敛速度极快，应用范围极广，包括化学驱、混相驱、热力驱等复杂模型和各种难以求解的系数矩阵。例如，可以解决传统的共轭梯度法解决不了的非对称矩阵问题及一般算法难以求解的病态系数矩阵。因此，目前这种类型的算法已成为油藏数值模拟中最先进的求解大型线性方程组的方法。这里，黑油模型的子程序CGIN中运用的是将系数矩阵的不完全LU分解和ORTHOMIN加速法结合的预处理共轭梯度法。 5.2.1 迭代求解的基本原理 设有方程组：  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 5) 并设为非奇异矩阵，则可以构造迭代公式：  MACROBUTTON MTPlaceRef \* M