第四章:补渗面边界(rsf4)条件包

在物理和工程领域中，非受限系统通常指的是一类地下水系统或含水层。让我为你解释一下：

• 非受限系统（Unconfined System）：这是一种地下水系统，其中地下水位与地表接触，没有上方的不透水层或隔水层。在非受限系统中，地下水可以自由上升或下降，受到地表的影响。

• 特点：非受限系统的特点包括：

  • 地下水位与地表高程相等。
  • 地下水可以通过渗透或补给进入系统。
  • 地下水流动受到地表的影响，例如降雨和蒸发。

• 应用：非受限系统在水资源管理、环境保护和地下水模拟中具有重要作用。例如，研究非受限系统的水位变化可以帮助我们了解地下水资源的可持续性。

在 MODFLOW-SURFACT 中，新的 RSF4 Package 旨在解决 RCH1 Package 在表示非受限系统时的局限性。让我为你解释一下：

• RCH1 Package：MODFLOW 中的 Recharge (RCH1) 包用于模拟地下水补给边界条件。然而，RCH1 Package 在表示非受限系统时存在一定的局限性。

• 非受限系统的问题：在非受限含水层饱和到地表时，含水层吸收供应的补给水的能力会减小，剩余的水会以地表径流的形式流失。

• RCH1 Package 无法处理这种情况。相反，它继续为含水层提供补给水（就像受限系统一样），使水头持续上升，远远超过地表。

• RSF4 Package：MODFLOW-SURFACT 的新 RSF4 Package 旨在消除这种物理上不现实的情况。RSF4 Package 允许将补给水引入地下水系统，前提是水位低于用户规定的池（积水）高程。

• 池高度：池高度对应于没有积水时的地表高程。

• 模拟结果：如果在任何时候，水位达到池高度，模拟只允许发生足够的补给以维持规定的池条件。剩余的补给不被接受到模拟域中。模拟的输出反映了由于系统饱和到池高度而导致的系统体积补给通量（L3 /T）的减少。请注意，RSF4 Package 与流量路由包（STR1 或 SFR1）是解耦的，因此这部分多余的补给不会以任何方式路由到指定的河流中。

RSF4边界条件包适用于涉及补给和大气渗漏面边界条件的任何情况。该包所需的输入数据包括补给速率，存储在二维（2-D）数组（RECH）中，每个水平单元位置一个元素，与之前的MODFLOW补给包相同。此外，为了允许非受限的补给条件，RSF4包需要池高度（PNDEL）作为输入。只要单元格的水位高度（HNEW）低于其池高度，整个补给速率（RECH）将应用于边界。如果在非线性迭代过程中，水位高度超过池高度，池高度将保持不变，系统的补给速率将减少到系统在池高度条件下可接受的最大补给速率。如果在模拟过程中的任何时候，水位低于规定的池高度，将全面应用补给速率到系统中。此外，如果池高度高于地表，存储项将通过使用修改后的Picard方法（Celia等，1990年）增加积水的附加体积。假定积水发生在单元格的整个区域内，除非打开了IPNDPOR标志。在这种情况下，网格块内的小湖泊和池塘的体积可以（在瞬态情况下）作为网格块大小的一部分额外读取到数组PNDPOR中。补给模块的逐单元流量文件还包括瞬态模拟的积水存储项的预算值。还应注意，流量文件或质量预算项中列出的“补给”的负值表示在渗漏面条件下发生的排放，来自“补给”边界节点。https://github.com/Umercia/rsf

4.3 快速、瞬态的补给条件

在MODFLOW中，补给和降水在一个应力周期内被视为恒定。对于快速的降雨事件，将一个短时间段分成多个应力周期在模拟过程中可能变得繁琐。因此，我们提供了一个选项，可以在模拟中独立于应力周期来改变补给或降水。瞬态补给是按区域应用的，RSF4包通过读