ECLPSE收敛性.doc

1. 介绍本文介绍如何改进收敛性，已经通过对数据体的简单修改来加速模拟器的计算。 某些情况下，进行油藏模拟的数据体需要非常小的时间步或耗费大量的运算时间来求解。修改已有的关键字或加入新的关键字有时能够在不改变结果的工程计算精度的基础上，非常有效地改进模拟器的计算速度。 我们将会展示具有该功能的关键字。 如果数据体不能够修改，那么某些时候可以通过调整收敛准则来改进模拟器的表现。运算的结果依赖于收敛准则以及其具体的参数设置，所以在某些方案中，需要在准确性和速度间进行平衡。 我们将会介绍如何检测这些方案以及如何选择合适的收敛准则。 大多数的建议都通用于所有的模拟器。 当在 ECLIPSE BLACKOIL, COMPOSITIONAL 和 THERMAL模拟器中，其使用的细节部分有不同点时，我们将会强调不同点并解释如何调整不同模拟器的数据体以改进表现。但我们使用 ECLIPSE 并没有特定指定 BLACKOIL, COMPOSITIONAL 或THERMAL，说明相关讨论都通用于所有的模拟器。 2. 时间步和迭代步的背景知识 一次 ECL 模拟计算 是由一个或多个报告步组成的 一个报告步 是由一个或多个时间步组成的 一个时间步 是由一个或多个非线性迭代组成的 一个非线性迭代 是由一个或多个线性迭代组成的 但你在构建ECL数据体的时候，需要指定报告步。所以你可以直接控制报告步的数量以及其间的时间跨度。 减少报告步的数量有时能够减少计算时间。 ECL具有默认的参数值，用于控制到下一个你指定的报告步所需要的时间步数量。通常情况下，默认值都能够很好的工作，但是某些时候，你需要调整某些默认值，以减少时间步数量。通常情况下，更少的时间步意味着更少的计算时间。在某些计算困难的数据体中，减少最大的时间步能够加速运算。 不同的默认值也控制着求解每个时间步所需要的非线性迭代数量。这些值通常情况下是不能修改的。在少数情况下，收敛准则的调整能够改进模拟器的表现。多数情况下，模拟器表现的巨大改进是需要通过分析导致非线性收敛问题的具体原因，再修改数据体来得到的。本文的主要部分即为解释如何避免这类型问题出现，以及如何找到并修正非线性迭代问题。 在线性迭代中出现问题时，通常情况下通过调整线性收敛控制来修正问题已经较迟了。我们将会解释在极端情况下哪些控制能够进行调整。但是最好的建议是通过控制时间步和非线性迭代来避免这样的问题。 2.1 时间步和迭代的报告 报告步、时间步和非线性迭代步的数量都可以在PRT和LOG－屏幕上出现的简要格式的输出（交互式运算）或者输出到 log 文件（批处理或者后台运算）中找到。 在Unix系统中，你可以使用 grep 命令来找到到所有必要的信息；在 PCs 上，你可以使用合适的编辑器来找到相关信息。黑油和组分模拟器（包括热采模型）采用不同的形式报告形式。 在 UNIX 系统中，某ECL 黑油 log文件叫做 BASE.LOG，使用命令 grep TSTEP BASE.LOG BASE.STEPS 将会创建 BASE.STEPS文件，该文件中每一行为一时间步的相关输出。 每行为如下格式： STEP 15 TIME= 400.00 DAYS ( +30.0 DAYS REPT 3 ITS) (4-FEB-2003) STEP 15 表示这是第15 时间步。 TIME= 400.00 DAYS 表示自模拟开始时间起已有400个模拟天。 +30.0 DAYS 表示当前的时间步为30天。 REPT 表示为什么选择30天的时间步长。REPT 表示已计算到某报告步。3 ITS 表示需要3个非线性迭代来求解30天长的时间步。 (4-FEB-2003) 为当前的模拟时间。 在UNIX系统中，组分或热采模拟器的 log 文件叫做BASE.LOG，使用命令 grep ; BASE.LOG BASE.STEPS 将会得到BASE.STEPS文件，该文件中每一行为一时间步的相关输出。 每行为如下格式： Rep ; 400.0 30.0 8.7838 .19498 1.4E05 32884. 1.2E06 4843.6 .00000 1.3E06 3 Rep 表示已计算到某报告步。 30.0 表示当前的时间步为30天。 下面的8项数据分别表示 气油比GOR, 含水WATER CUT, 油 气 水的生成速度，平均油藏压力和气水的注入速度。 最后的3表示需要3个非线性迭代来求解30天长的时间步。 如果使用了AIM选项那么每行将会有额外的数据，Rep ; 400.0 30.0 8.7838 .19498 1.4E05 32884. 1.2E06 4843.6 .00000 1.3E06 3