第7章. 稳态和瞬态模拟

时间，那黑暗而狭窄的地峡，位于两个永恒之间。

——查尔斯·卡莱布·科尔顿

而时间……必然会有一个终点。

——《亨利四世》第一部

在前六章中，我们描述了一个数值模型的设计，并提及了稳态和瞬态（时间相关）条件。几乎所有的模型都以稳态模拟为起点，而一些建模目标也需要进行瞬态模拟。在本章中，我们讨论了稳态模拟的特征性质，其中一些在前面已经介绍过（第4至6章的问题都被构建为稳态问题），但本章的大部分内容涉及瞬态模拟。

7.1 稳态模拟

在稳态模型中，方程（Eqn (3.12)）中的 是0，水文参数（第5.4节）和计算的水头和通量在时间上是恒定的（时间不变）。仅仅使用稳态解通常就足以满足许多建模目标，比如分析平均地下水流动模式和流速；估算从流失河流中的平均年漏水量；计算区域性的水位梯度；模拟受长期抽水影响的流向。稳态模型还可以用来预测时间平均应力的影响，如长期干旱或预测的长期抽水。此外，由于大多数瞬态模型使用稳态解作为初始条件（第7.4节），稳态模型通常是瞬态建模的第一步。在将输入输入到执行地下水流动代码的图形用户界面（GUI）（第3.6节）时，稳态模拟通常是默认的模拟选项。

7.1.1 初始水头

数值模型通常是通过迭代求解的（第3.5和3.7节），在稳态模拟开始时，需要在所有节点上设定初始水头。（如第3.5节所述的直接数值解通常不需要起始水头。）初始水头影响解收敛所需的迭代次数，但在大多数情况下，只要解已经收敛到足够的精度（第3.5节），最终的稳态解通常不依赖于起始水头的值。

起始水头可以在整个网格/网格上保持恒定，或者可以指定任意值。任意的起始水头通常基于地表高程或指定的内部和周边边界条件。

当起始水头相对接近最终水头解时，解需要的迭代次数较少。选择不当的起始水头可能会导致迭代求解器在初始迭代过程中超出解，从而导致干燥节点（第4.5节），影响收敛和最终水头解。

因此，始终应在解的范围内指定起始水头。

起始水头通常基于指定的指定初始水头边界条件，并且几乎总是设置在模型基底底部的底部高程之上，在三维模型中通常设置在顶层的底部高程之上。

一些图形用户界面（GUI）通过在内部和周边边界处使用指定的水头值或使用最上层的顶高来自动生成默认的起始水头。

7.1.2 边界条件

在稳态模型中，内部和周边边界条件（第4.2和4.3节）是时间不变的，材料特性和水文参数的值也是时间不变的（第5.4节）。

在指定边界条件时必须小心，因为稳态解受其边界强烈影响。

这是因为稳态解没有由储存贡献的水，因此所有的水都来自内部和周边边界（以及分布在区域内的源）。

或者换句话说，稳态状态表示无限时间内的水头，这意味着水文应力达到周边边界，边界处的条件在整个模型中完全传播。

7.1.3 表征稳态条件

设计为通用或简单筛选模型的解释性稳态模型可能不需要校准（第1.3节），但大多数稳态模型都经过校准（第9章），以适应地下水系统的实测观测数据，通常是水头和流量。

在需要校准时，建模者必须整理一组实测观测数据，以与模拟的稳态水头（图9.3）和流量进行比较。

然而，需要进行判断，因为地下水系统的实测数据很少处于真正的稳态；相反，水头和流量会不断变化，以响应自然和人工充沛和排泄的变化（图7.1）。

因此，在大多数情况下，建模者必须从固有的瞬态实测数据中提取代表性的稳态校准目标。有两种组装稳态校准数据集的选项。

1. 代表平均年、季节、月份或模型目标重要的其他选定时间段的水头的时均值，可以从观测井的水头测量中计算得出（图7.2）。

平均基流可以从流量测量中计算得出（图9.4(b)）。

判断一个时间段是否代表稳态条件的标准之一是评估该时段初期的场地测量值是否类似于该时段末尾的测量值。时段开始和结束时水头的相似值表明在该时段内储存的变化可以忽略不计，并且该时段内的平均水头可以使用（例如，Dripps等人，2006年）。

图7.1：箱形图显示美国俄勒冈州温带气候下水文地质单元观测井中水头的均值、中位数和范围（Snyder, 2008）。

图7.2：新西兰监测井的水文图。使用1990年至2000年的时期得出了用于稳态模型的平均水头目标（以虚线表示）（Scott和Thorley，2009）。

2. 当水头随时间变化很小时，可以直接使用来自代表性时间段的实测数据作为稳态校准数据。理想情况下，这些测量反映的是水位和流量几乎恒定的时间段，即准稳态或伪稳态。例如，某个代表性时间测得的水头可能被认为代表了长期平均条件。另一种选择是使用季节性的准稳态水头。例如，在一些北纬较高的地区，山谷地区的地下水系统主要受到雪融充沛的影响，因此三月中旬至晚期可能是一个伪稳态期，即在冬季结束之前，春季充沛之前水文图形趋于平稳。在其他温带气候中，秋季中旬可能是一个准稳态的时期（图7.3）。界定和连续稳态模拟（第7.2节）校准到准稳态水头。

稳态模型的输出是水头的空间数组，通常由代码输出为各层的水头的二维数组，并伴随着内部和周边边界的相关流量。图形用户界面（GUI）和/或可视化软件以面积（图2.5和图9.3）和剖面（图4.22(c)）视图显示各层的模拟水头。模拟的稳态水量平衡（第3.6节）报告的单位与用于参数值输入的单位一致（例如，体积流量以m³/d、ft³/d、m³/s或其他单位报告）。

图7.3：得克萨斯州特立尼地下水层监测井的水文图（2009年7月至2012年7月），显示2011年夏末抽水率较低、秋雨降临之前的伪稳态条件（改编自中德克萨斯地下水保护区，The Hydro Blog，2012年8月，http://www.centraltexasgcd.org/the-hydro-blog/）。

7.2 稳态还是瞬态？

在Box 4.1中，我们讨论了模型设计中的一个基本决策，即是否需要一个三维模型。在本节中，我们讨论另一个基本决策——是否需要一个瞬态模型。稳态模型在创建、执行和后处理方面本质上比瞬态模型更简单，仅仅因为这个原因，通常情况下如果它们能够足够地满足建模目标，就更受欢迎。稳态模型只需要一组校准数据并产生一组结果，而瞬态模型则需要多组校准数据并产生比稳态模型更多的输出。瞬态模型还需要额外的输入，即储存参数（第5.4节）和初始条件（第7.4节）的值。

然而，在许多情况下，瞬态模型显然是必需的，例如，计算水位如何随着日常或季节性抽水率的变化而波动，季节性的地下水提取导致