第4章：补给-渗流面边界（RSF4）条件包

4.1概述

MODFLOW代码使用recharge（RCH1）软件包适应地下水补给边界条件。但是，RCH1包在表示未受限制的系统方面有一定的局限性。如果无侧限含水层饱和到地表，含水层吸收补给的能力就会降低，剩余的水会作为地表径流排
出。RCH1软件包无法处理这种情况。相反，它继续向含水层提供补给（如在密闭系统中），水头不断上升，远远高于地表。MODFLOW-SURFACT的新RSF4软件包旨在消除这种物理上不现实的条件。如果地下水位低于用户规定的水池（积水）高程，RSF4包允许向地下水系统补给水。水池高程与无积水时的地面高程相对应。如果在任何时候，地下水位达到水池高程，模拟只允许进行尽可能多的补给，以保持规定的水池条件。剩余的充值不被接受进入模拟域。模拟的输出反映了由于系统饱和到积水高度而导致的系统体积补给通量（L3/T）的减少。请注意，RSF4包与流路由包（STR1或SFR1）断开连接，因此，该剩余充值不会以任何方式路由到规定的流中。
此外，对于MODHMS模拟，无侧限渗流选项与陆上流动域相冲突，因为陆上流动区域处于活动状态，因此，如果RSF4软件包中使用无侧限渗漏选项，则应在这些位置设置非常高的积水高程，以防该软件包干扰陆上流动。
新RSF4软件包的额外功能允许方便地输入瞬态补给/降水，而不依赖于MODFLOW的“应力期”概念。当使用此选项时，可以通过单独的文件输入充值时间序列，该文件确定域上的充值如何随时间变化，而不考虑压力期设置，以在充值快速且经常在模拟期内变化时提供输入便利性和通常较小的输入文件大小。
除了地表补给排水条件外，RSF4软件包还可用于通过规定渗流面边界节点的高程来模拟渗流面边界条件。因此，零流量出现在边界处，直到地下水位上升到渗流面高程。如果地下水位达到渗流面高程，含水层排水以维持渗流面高程（即大气压）条件。如果在模拟过程中地下水位下降到渗流面高程条件以下，则在渗流面边界处保持零流量。
RSF4包的无侧限渗流选项仅用于BCF4包的严格可变饱和模拟选项（即输入LAYCON=40或
43）。RSF4软件包的功能介绍如下，包括对公式的简要讨论和演示示例问题。

 

除了地表补排条件外，RSF4包还可以通过规定渗流面边界节点的标高来模拟渗流面的边界条件。

因此，在地下水位上升到渗流面高度之前，边界处的通量为零。

当地下水位达到渗流面标高时，含水层排水以维持渗流面标高(即大气压)的条件。

在模拟过程中，如果地下水位低于渗流面标高条件，则渗流面边界处保持零通量。

RSF4包的潜水渗流选项仅用于BCF4包的严格可变饱和模拟选项(即输入LAYCON = 40或43)。

下面介绍RSF4包的功能，包括对公式和演示示例问题的简要讨论。

4.2 rsf4边界条件包的制定

RSF4套件可适用于任何涉及补给和大气渗流面边界条件的情况。

该包所需的输入数据包括补给率，该数据存储在二维(2-D)数组(RECH)中，每个水平单元位置有一个元素，与之前的MODFLOW补给包一样。此外，为了允许潜水的补给条件，RSF4包需要池高度(PNDEL)作为输入。

只要一个单元的地下水位高度(HNEW)保持在其池高度以下，整个补给速率(RECH)就应用于边界。

如果在非线性迭代过程中，地下水位超过池的水位，则保持池的水位，系统的补给减少到在池高程条件下系统可接受的最大补给率。

如果在模拟过程中的任何时候，地下水位低于规定的池高程，则对系统施加满补给率。

此外，如果池高程高于陆地表面，则使用改进的Picard方法(Celia et al, 1990)，池水的额外体积可以延长储存期。

除非打开IPNDPOR标志，否则假定池将发生在单元格的整个区域上。

在这种情况下，网格块内的小湖泊和池塘可能会将其体积(在瞬态情况下)作为网格块大小的一小部分来计算，这些部分会在数组PNDPOR中额外读取。

充电模块的逐个单元流文件还包括用于瞬态模拟的池存储期的预算值。

应该进一步指出的是，在逐单元流动文件或质量预算项中列出的“补给”的负值表示在渗流面条件下从“补给”边界节点发生的放电。

RSF4软件包还可以模拟侧向渗流面边界条件。

在水位标高高于渗流面标高之前，对于划定为可能的渗流面节点的区域内的任何节点，该规范规定系统进出(或流出)的通量为零，此时在渗流面标高条件下，水从系统排出。

4.3 迅速、瞬态补给条件

MODFLOW中的压力边界条件使用边界包（Boundary Package）来表示。

边界包可用于模拟不同类型的边界条件，如指定流量、指定总流量、指定压力等。

边界包是MODFLOW输入文件的一部分，用于描述模型的外部边界及其特性。

在MODFLOW中，有几种类型的边界包，其中一些常见的包括：

1. **General-Head Boundary (GHBD) Package:** 用于指定模型中的一般水头边界条件。它允许用户指定每个单元的水头值。

2. **River (RIV) Package:** 用于表示模型中的河流或渠道。用户可以指定河流的位置、水头和其他相关属性。

3. **Well (WEL) Package:** 用于表示井，可以指定井的位置、水头和抽水/注入速率。

4. **Drain (DRN) Package:** 用于表示水下排水沟或排水管。它可以模拟地下水位高于排水沟口时的排水。

5. **Recharge (RCH) Package:** 用于指定地下水补给的区域和速率。

6. **Evapotranspiration (EVT) Package:** 用于表示蒸发和植物蒸腾的区域和速率。

7. **Lake (LAK) Package:** 用于表示湖泊或水库，可以指定湖泊的位置、水头和湖泊与地下水之间的交互。

这些包的选择取决于模型的特定要求和研究问题。用户可以根据模型的需要选择适当的边界包，并提供相关的参数和数据。

4.4 验证与应用示例

提供了两个示例问题，用于演示对RSF4包进行验证和应用。

第一个问题涉及一个非均质含水层，该含水层接收来自地表补给的水，并排入长平行排水渠。

第二个问题涉及一个堤坝沿垂直渗流面排水的情况。

这些问题将与分别针对相应情况的解析和数值解进行比较。

4.4.1 Problem 1—Flow to Parallel Drains in an Unconfined Aquifer问题1 -潜水含水层中并联排水的流动

这个问题涉及到在各向同性和均质无限不受限制含水层中平行排水的流动，受降水的垂直补给，如图4.1所示。

由于问题的二维性质，只能分析一个垂直切片。排水之间的距离为1,000英尺，含水层的深度为75英尺。含水层的初始水头为50英尺。

该区域受到恒定的降水补给，为0.01英尺/天。

排水维持50英尺的恒定水头。含水层参数包括：水力导数（K）= 1英尺/天，比水容（Sy）= 0.2。

建模方法如前所述，问题是二维的，只选择了一个切片（NROW = 1），宽度为1英尺进行分析。此外，在排水之间存在对称线（在x = 500英尺处），因此只对域的第一半部分进行离散化，如图4.2所示。

域被离散成10个厚度为7.5英尺的层和20个长度为25英尺的列。

含水层的底部位于基准面（0英尺），顶部位于75英尺。每个层次的VCONT值（垂直导数除以两个相邻层之间的层间距）为0.1333 d–1。

每个层次的初始水头值和排水节点处的边界水头均为50英尺。

所有层次的LAYCON值为43，以利用新的可变饱和流公式，并选择SIP包来解决方程。

模拟了两种情况：Case 1在域上施加0.01英尺/天的补给，Case 2在域上施加0.02英尺/天的补给。

结果与讨论模拟的结果分别在图4.3和4.4中显示，对于Case 1和Case 2，分别显示了不同时间值。

对于Case 1模拟，水位保持在地表以下，并且整个补给在整个模拟过程中都被应用到域中。

该案例的稳态解与Bear（1979）的解析解进行了比较，如图4.3所示，显示出良好的匹配。对于Case 2，水位在流动域的某些部分达到地表。在这些地方，水位然后维持在池的高程（地表）处，地下水系统的补给减少，如模拟输出所反映。

4.4.2 Problem 2—Seepage Through a Square Embankment  问题2 - 正方形筑坝的渗流

这个问题涉及到正方形土堤的稳定渗透。

选择这个问题是为了测试RSF4 Package模拟渗透面的能力。

Cooley（1983）提供了这个问题的严格（可变饱和流）数值解。

几何形状、边界条件和各向同性土壤参数如图4.5所示。

一个边长为10米的各向同性土堤，具有水平导数性质K = 0.01 m/d，在左侧保持10米的水位，右侧维持2米的水头。

然而，沿着这个边界可能会发展出一个渗漏面。由于这个问题的二维性质，只分析问题的垂直剖面。

图4.1 示意图，描述不受限制含水层中流向平行排水的流动。

图4.2 用于模拟平行排水流动的离散化方案。

图4.3 平行排水问题的数值和解析解比较（补给 = 0.01英尺/天）。

图4.4 平行排水问题水位的瞬时变化（补给 = 0.02英尺/天）。

图4.5 正方形土堤的稳定渗透。

Modeling Approach

问题的配置是二维的。

因此，只选择了一个切片（NROW=1），宽度为1米（DELC=1米）进行分析。

在垂直方向上，域被离散成12列和10个1米厚度的层。

在左右边界处，使用0.1米的列宽（DELR）以准确捕捉边界效应。

在左边界施加10米的水头，在右边界的底部两个单元格上施加2米的水头。

右边界沿着剩余单元格可能有渗漏面，并且通过RSF4 Package将其纳入模拟。

MODFLOW所需的VCONT值为0.01 d–1。

每个层次的输入LAYCON值为43，以使用BCF4 Package的可变饱和流公式，SIP Package用于解决矩阵方程。

模拟所使用的输入数据文件列在附录D中。

结果和讨论进行了以初始水头为10米的稳态模拟。

图4.6显示了模拟的水位高程和域底部的压力水头。

这些结果与Cooley（1983）和Huyakorn等人（1986）使用有限差分和有限元模型获得的结果非常相符。

4.3 Rapid, Transient Recharge Conditions

MODFLOW中的一个应力周期内，补给和降水被视为恒定的。

对于短时期的强降雨事件，将一个短时间段分割成多个应力周期在模拟中可能变得繁琐。

因此，提供了一个选项，在模拟期间独立于应力周期内变化补给或降水。

瞬时补给被区域性地应用，RSF4包通过读取一个指示瞬时补给应用区域数量的索引MXZRCH。

如果MXZRCH为零，则瞬时补给包在模拟中未被调用。

然而，如果存在瞬时补给，RSF4包会读取一个额外的整数数组，用于标识每个区域的面积位置。

区域值为零表示在该位置没有瞬时补给区域，并且将默认的补给值应用于该位置，就像在RECH数组的应力周期中读取的那样。

在区域值不为零的位置，将在.RECH数组的应力周期中为该区域的适当值添加默认补给值，并使用单独文件中提供的瞬时补给时间序列数据。

此文件的默认扩展名为.RTS，包含了所有区域的补给时间序列表格化的数据。

文件的每一行包含一个事件的开始时间（TSTART）、结束时间（TEND）、一个乘法因子（用于单位转换），然后是每个区域的补给/降水速率的一列。

如果后续的TSTART值与前一行的TEND值不匹配，则在临时期间对每个区域应用零补给。

因此，不需要将长时间间隔的休息期作为.RTS文件中的单独行包含。

事件之间或事件内的时间段不需要固定，因为该模块设计用于自适应时间步进方案。

自适应时间步进方案确保模拟精确地达到每个TSTART和TEND，以确保准确表示事件。

在事件内或事件之间，自适应时间步进方案试图最大化计算效率。

值得注意的是，在RECH数组中提供的补给会与RTS文件中提供的值相加，以在任何节点提供总的补给。

因此，通过RSF包的RECH数组可以按应力周期基础提供应用灌溉或其他人为补给的同时，通过RTS文件可以提供快速变化的降水事件。

4.4 Verification and Application Examples

提供两个示例问题以演示RSF4包的验证和应用。

第一个问题涉及到一个不受限制的含水层，接收来自地表的补给水，并排入长平行排水沟。

第二个问题涉及到一个沿着垂直渗流面排水的土堤。这些问题将与各自情境下的解析和数值解进行比较。

4.4.1 Problem 1—Flow to Parallel Drains in an Unconfined Aquifer

这个问题涉及到在各向同性和均质无限不受限制含水层中，受到来自降水的垂直补给的情况下，流向平行排水的流动，如图4.1所示。

由于问题的二维性质，只能分析一个垂直切片。排水之间的距离为1,000英尺，含水层的深度为75英尺。

含水层的初始水头为50英尺。

该区域受到恒定的降水补给，为0.01英尺/天。排水维持50英尺的恒定水头。

含水层参数包括：

水力导数（K）= 1英尺/天
比水容（Sy）= 0.2

Modeling Approach

正如前面所述，这个问题是二维的，只选择了一个切片（NROW = 1），宽度为1英尺（DELC = 1英尺）进行分析。

此外，在排水之间存在一个对称线（在x = 500英尺处），因此只对域的第一半部分进行了离散化，如图4.2所示。

域被离散成10个层次，每个层次的厚度为7.5英尺，以及20列，每列的长度为25英尺。

含水层的底部位于基准面（0英尺），顶部位于75英尺。

每个层次的VCONT值（垂直导数除以两个相邻层之间的层间距）为0.1333 d–1。

每个层次的初始水头值和排水节点处的边界水头均为50英尺。

对于所有层次使用LAYCON值为43，以利用新的可变饱和流公式，并选择SIP包来解决方程。

进行了两种情况的模拟：Case 1在域上施加0.01英尺/天的补给，Case 2在域上施加0.02英尺/天的补给。

Results and Discussion

模拟的结果分别在图4.3和4.4中显示了Case 1和Case 2的不同时间值。

对于Case 1模拟，水位保持在地表以下，整个补给在整个模拟过程中都被应用到域中。

该案例的稳态解与Bear（1979）的解析解进行了比较，如图4.3所示，显示出良好的匹配。

对于Case 2，水位在流动域的某些部分达到地表。

在这些地方，水位然后维持在池的高程（地表）处，地下水系统的补给减少，正如在模拟输出中所反映的。

4.4.2 Problem 2—Seepage Through a Square Embankment

这个问题涉及到一个正方形土堤的稳定渗流。

选择这个问题是为了测试RSF4 Package模拟渗透面的能力。

Cooley（1983）提供了对这个问题的严格（可变饱和流）数值解。

几何形状、边界条件和各向同性土壤参数如图4.5所示。

一个边长为10米的各向同性土堤，具有水平导数性质K = 0.01 m/d，在左侧保持10米的水位，右侧维持2米的水头。

然而，沿着这个边界可能会发展出一个渗漏面。由于这个问题的二维性质，只分析了问题的垂直截面。

Figure 4.1 Schematic Diagram Depicting Flow to Parallel Drains in an Unconfined
Aquifer.图4.1 示意图，描述在不受限制的含水层中流向平行排水的流动。

Figure 4.2 Discretization Scheme Used for the Simulation of Flow through Parallel
Drains.图4.2 用于模拟流经平行排水的离散化方案。

Figure 4.3 Comparison of Numerical and Analytical Solutions for the Parallel
Drain Problem (Recharge = 0.01 ft/d).图4.3 平行排水问题的数值解和解析解比较（补给 = 0.01英尺/天）。

Figure 4.4 Transient Movement of the Water Table for the Parallel Drain Problem
(Recharge = 0.02 ft/d).图4.4 平行排水问题水位的瞬时变化（补给 = 0.02英尺/天）。

Figure 4.5 Steady Seepage through a Square Embankment.图4.5 正方形土堤的稳定渗流。

Modeling Approach

问题的配置是二维的。因此，只选择了一个切片（NROW=1），宽度为1米（DELC=1米）进行分析。

在垂直方向上，域被离散成12列和10个1米厚度的层。

在左右边界处，使用0.1米的列宽（DELR）以准确捕捉边界效应。

沿左边界施加10米的水头，在右侧边界底部的两个单元格上施加2米的水头。

右边界沿着其余的单元格可能是渗漏面，并通过RSF4 Package纳入模拟。

MODFLOW所需的VCONT值为0.01 d–1。

每个层次使用LAYCON值为43，以使用BCF4 Package的可变饱和流公式，并选择SIP包来解决矩阵方程。用于模拟的输入数据文件在附录D中列出。

Results and Discussions

使用初始水头为10米进行了问题的稳态模拟。

图4.6显示了在域的基底沿着模拟的水位高程和压力水头。

这些结果与使用Cooley（1983）和Huyakorn等人（1986）的有限差分和有限元模型获得的结果非常吻合。

4.5 Input Instructions for the RSF4 Package

4.5.1 RSF4 Input

对Recharge Seepage Face（RSF4）包的输入是从IUNIT(8)中指定的单元读取的。

RSF4需要的其他输入变量显示在阴影框中。

请注意，新的RSF4包括原始MODFLOW（McDonald和Harbaugh，1988）的Recharge（RCH1）包。

瞬时补给条件的输入通过两个文件进行。第一个是在IUNIT(8)上的补给渗透面包输入文件，该文件指示模型中瞬时补给的存在以及补给区域的数量。

第二个文件包含每个区域的补给时间序列（默认扩展名为.RTS）。

此文件是无格式的，当应用瞬时补给时，它会在内部打开。

1

Data:	NRCHOP	IRCHCB	IUNCNF	MXSEEP	MXZRCH	IPNDPOR
Format:	I10	I10	I10	I10	I10	I10

**为每个应力周期输入2到6项。

2.

Data:	INRECH	INIRCH	INPNDEL	INSEEP
Format:	I10	I10	I10	I10

**仅当INRECH≥ 0时输入第3a项。

3a.

Data:	Utility
Module:	RECH(NCOL, NROW) U2DREL. Note that U2DREL is the utility module that reads two-dimensional array (i.e., RECH for item 3a).

**仅当INRECH > 0且MXZRCH > 0时输入第3b项。表格化的分区补给值将添加到上述第3a项中读取的默认补给值。

3b.

Data:	Utility
Module:	IZNRCH (NCOL, NROW) U2DINT

*仅当补给选项（NRCHOP）等于2且INIRCH≥ 0时输入第4项。

4.

Data:	Utility
Module:	IRCH(NCOL, NROW) U2DINT

**仅当使用非均质选项（IUNCNF> 0）且INPNDEL ≥0时，输入第5项（包括5a和5b）。

5a.

Data:	Utility
Module:	PNDEL(NCOL, NROW) U2DREL

**仅当淤积孔隙度标志开启（IPNDPOR=1）时输入第5b项。

5b.

Data:	Utility
Module:	PNDOR(NCOL, NROW) U2DREL

**仅当INSEEP > 0时输入第6项。

6.

Data:	Layer	Row	Column	Elevation
Format:	I10	I10	I10	F10.0

（输入项6通常包含每个渗流面边界单元的一条记录。如果INSEEP为负或零，则不读取输入项6。）输入说明字段的解释：

- NRCHOP - 是补给选项代码（标志）。

补给速率在二维数组RECH中定义，每个垂直列有一个值。

因此，补给应用于每个垂直列中的一个单元，选项代码确定选择哪个单元用于补给。

  - 1 - 补给仅应用于顶部网格层。

  - 2 - 补给的垂直分布在数组IRCH中指定。

  - 3 - 补给应用于每个垂直列中最高的活动单元。

常头节点截取补给并阻止更深的渗透。

- IRCHCB - 是流动模拟的标志和单元编号。

  - 如果IRCHCB > 0，则是记录每当ICBCFL（请参见输出控制）设置时，将记录单元对单元的流动项的单元编号。

  - 如果IRCHCB ≠ 0，则不会打印或记录单元对单元的流动项。

- IUNCNF - 是不受限选项代码（标志）。

  - 如果IUNCNF > 0，则采用补给渗透面边界条件，补给（RECH）将应用到水位达到PNDEL数组中描述的集水高程的地方。

  - 如果IUNCNF ≠ 0，则总是像受限情况一样应用指定的补给（RECH）。

- MXSEEP - 是一次激活的渗漏边界单元的最大数量。

- MXZRCH - 是输入补给的分区时间序列标志。

  - 如果MXZRCH = n，则补给随时间变化，n是模拟中使用的总补给区域的数量。

  - 如果MXZRCH = 0，则MODFLOW的输入结构用于补给，只在应力周期内发生变化。

注意：
  - MXZRCH > 0表示补给的变化具有不同（可能小得多）的时间尺度，比MODFLOW的应力周期小，将为每个区域提供独立的补给时间序列（在文件*.RTS中详细描述）。

  - 区域索引数组标识位于特定区域内的网格块。索引值为零允许对于该应力周期的MODFLOW默认应用补给。

  - 补给的时间序列中的区域位置在应力周期之间可以变化。因此，这些区域本身可以根据应力周期的基础进行变化。这由标志INRECH控制，该标志还控制每个应力周期的补给率的读取。

- IPNDOR - 读取淤积孔隙度的标志。

  - 如果IPNDPOR = 0，则不读取淤积孔隙度，假定完整单元面积被淤积。

  - 如果IPNDPOR = 1，则读取淤积孔隙度数组。淤积孔隙度是由湖泊/池塘浸湿的总单元面积的比例，并用于确定总淤积储存，如果发生淤积的话。

- INRECH - 是RECH读取标志。

  - 如果INRECH ≠ 0，则读取补给速率（RECH）的数组。

  - 如果INRECH < 0，则使用前一个应力周期中的补给速率。

注意：如果MXZRCH > 0，则INRECH还控制用于分区索引（IZNRCH）的输入数组。因此，区域位置可能在每个应力周期中发生变化。

- INIRCH - 是IRCH读取标志。当NRCHOP为2时，如果INIRCH ≠ 0，则读取层编号（IRCH）的数组。如果INIRCH < 0，则重用前一个应力周期中使用的数组（IRCH）。

注意：当NRCHOP为1或3时，将忽略INIRCH。

- INPNDEL - 是当使用不受限选项时（IUNCNF > 0）的PNDEL读取标志。

  - 如果INPNDEL ≠ 0，则读取集水高程（PNDEL）的数组。

  - 如果INPNDEL < 0，则使用前一个应力周期中的集水高程。

- INSEEP - 是渗漏边界

单元的标志和计数器。

  - 如果INSEEP < 0，则将重用前一个应力周期中的渗漏面边界数据。

  - 如果INSEEP ≠ 0，则INSEEP是当前应力周期中渗漏面边界的数量。

- RECH - 是补给速率（L/T）。只有当INRECH大于或等于零时才读取。

- IZNRCH - 是标识每个节点的区域值的整数数组。

  - 默认情况下，每个应力周期中读入RECH数组的补给值将应用于具有零区域值的单元格。

- IRCH - 是定义应用补给的每个垂直列中的层的层数数组。只有当NRCHOP为2并且INIRCH大于或等于零时才读取。

- PNDEL - 是定义每个垂直列中的集水高程的集水高程数组。

通常，如果地表以上没有水池，集水高程与地表高程相一致。仅在IUNCNF大于零且INPNDEL大于或等于零时才读取。

- PNDPOR - 是定义当水位高于域顶部时由池塘/湖泊占据的单元面积的分数的淤积孔隙度数组。

在瞬时模拟中，该分数用于跟踪淤积水的体积/速率。

- Layer - 是表示渗漏面边界的单元的层号。

- Row - 是表示渗漏面边界的单元的行号。

- Column - 是表示渗漏面边界的单元的列号。

- Elevation - 是渗漏面的高程。

4.5.2 Transient Recharge Time-Series File

瞬时补给时间序列文件在应用瞬时补给条件时在FORTRAN Unit 96上被内部打开。

该文件采用自由格式，每个事件在模拟中占据一行。每行上的数据条目如下：

数据：
- TSTARTR - 降雨/降水事件的开始时间

- TENDR - 降雨/降水事件的结束时间

- TMULTR - 应用于给定事件的所有区域的补给速率的乘法因子。这种缩放对于单位转换很方便。

- RZ(1) - 该事件的区域1的补给速率

- RZ(2) - 该事件的区域2的补给速率

- ...

- RZ(MXZRCH) - 该事件的区域MXZRCH的补给速率

格式：
- 自由格式

输入中使用的字段的解释。

请注意，TSTARTR和TENDR是相对于总模拟时间而不是应力周期时间的，当在模拟中发生多个应力周期时。在示例的RTS文件中，域被划分为8个补给区域。

在模拟的前2.375天内没有发生降雨或降水，因此模拟会在第一个补给事件之前迅速进行，根据收敛行为调整其时间步长大小。

在2.375天时，某些区域发生了一个小时的降雨事件，持续到2.417天。

在2.417天时，补给事件的强度（和/或区域）发生变化，在文件的第二行中发生这些条目，依此类推。

请注意，在2.5到4天之间没有降雨事件，并且对于这种情况不需要文件中的条目。

TSTARTR(d)	TENDR(d)	TMULTR	RZ(1)ft/d	RZ (2)ft/d	RZ-3	RZ-4	RZ-5	RZ-6	RZ-7	RZ-8
2.375	2.417	1	0	0	0.225	0.177	0.033	0	0	0
2.417	2.458	1	0	0	0	0.13	0.13	0.022	0	0
2.458	2.5	1	0	0	0	0	0.088	0	0	0
4	4.042	1	0.169	0.551	0.551	0.337	0.178	0	0	0
4.042	4.083	1	0.151	0.447	0.333	0.106	0.08	0	0.691	0.633
4.083	4.125	1	0.1	0.05	0	0	0	0.15	0.379	0.3
4.125	4.167	1	0	0	0	0	0	0	0.083	0.022

注:MXZRCH =最大充值区数= 8

这个表格是一个RTS文件的样本，包含了不同时间段的补给区域和降雨强度。每一行都代表一个时间段，具体的列包括：

- `TSTARTR(d)`：补给或降雨事件开始的时间（天）
- `TENDR(d)`：补给或降雨事件结束的时间（天）
- `TMULTR`：时间步长的倍增因子
- `RZ(1)ft/d` 到 `RZ-8`：不同补给区域的补给速率（单位是英尺/天）

具体来说，对于第一行：
- 从时间2.375天开始，到2.417天结束，补给区域RZ(1)的速率是0.225英尺/天，RZ(2)到RZ-8的速率都是0。
- `TMULTR` 是1，表示时间步长在此期间不发生变化。

然后，文件中有几行描述了接下来的时间段，每一行都包括新的补给区域和速率。这种格式允许模型在模拟期间动态地改变补给条件，以更好地模拟真实情况。