如何提高数值模拟计算性能-LearnCMG.PDFVIP

第十一期：如何提高数值模拟计算性能（1）                    如何提高数值模拟计算性能（1） 本期内容是如何诊断较差的数值模拟计算性能，以及相关的建议和改进方法。 1、如何读计算记录信息？ 除了常规写入输出结果文件（*.OUT文件）的模拟结果信息外，STARS将每个模拟 步长的计算信息写到屏幕上（采用前台运行，即用“Run Immediately” 的方式提交任务） 如下图；或者将计算信息直接写入日志文件（*.log）中（采用后台运行，即用“Submit to Scheduler”的方式提交任务）。 整个记录信息分为：Time Step (时间步) 、Time （时间）、Production （生产）、Injection （注入）、Mat Bal Err （物质平衡误差）以及Maximum Changes （最大变化）六大部分。 1) Time Step (时间步) ：本部分共有四列，分别是No. (timestep number 时间步数）, Size Days （时间步长）, IT （迭代次数），CUT （截断数-收敛失效的次数）； 2) Time （时间）：主要是写明本时间步的时间及日期； 3) Production （生产）：显示油、气、水的产量以及气油比（GOR ）和含水率（Wat Cut ）指标； 4) Injection （注入）：显示注气量、注水量； 5) Mat Bal Err （物质平衡误差）：以百分数形式表示本步长内物质平衡误差； 6) Maximum Changes （最大变化量）显示本步长内压力、饱和度（油、气、水三 相）及温度的最大变化值。 2、数值模拟计算是怎么一回事？ 数值模拟计算是用计算机数学模型来表示、表达及模拟地下油藏的物理模型。数学 1   第十一期：如何提高数值模拟计算性能（1）                    模型模拟物理模型，不但需要精准的数学计算能力，还需要对物理过程的准确描述。这 个过程就是将是物理模型用数学模型进行时空离散，也就是说用网格去离散空间，用时 间步去离散时间。一个成熟的商业化数值模拟软件是经过多年的开发、测试及改进逐渐 形成的，就STARS而言，内部有几百万行的源程序。对于使用者来说，我们无需知道软 件内部复杂的程式，而只需要了解它如何为我们工作的就足够了。 首先，达西定律、控制方程及边界条件是模拟软件的基础，这些方程都是非线性化 的偏微分方程，模拟器需要在每一个时间步做三件事情，1）确定时间步控制，确定步 长；2 ）用牛顿迭代将非线性方程组线性化，实际也就是去微分的过程，我们称之为外 循环，形成线性的大规模矩阵；3 ）利用核心算法器（SOLVER ）求解矩阵，例如STARS 的普通版算法器是AIMSOL ，并行版的算法器是Parasol ，算法器的速度及稳定性是衡量 一个数模软件性能的重要指标。 3 、每一个时间步长（Timestep Size ）是如何确定的？ 每一个时间步的步长大小是由模拟器自动确定的，通常由受下几个因素影响： a ）本时间步与上一个时间步压力、饱和度及温度等值的最大变化量 如果在两个时间步内，压力饱、和度剧烈变化，会导致步长变小。这可以通过很多 经验证明，例如，我们通常在注入井开井时要把时间步缩小到DTWELL （最小步长）， 是因为此时井所在网格内的压力、饱和度会剧烈变化。另外，这对建模的模型质量提出 了要求，就是说相邻网格的空隙体积（PV ）尽量不要变化太大，保持和谐、均匀的变 化。如果一个网格的PV是1000立方米，附件的另