本文介绍了MODFLOW中的ETmax处理方法，特别是在处理瞬态ETmax时，如何通过独立于压力周期的变化、EVT和ET2包的输入，以及使用.ETS文件存储时间序列数据来实现精确的模拟。文章详细解释了相关输入参数和文件结构，强调了自适应时间步进方案的应用。
G.1 一般
| ETmax，蒸散蒸发的最大速率，在MODFLOW的一个压力周期内被认为是恒定的。 |
| 对于ETmax的快速变化，将一个短时间周期分解为模拟期间的几个压力周期将变得繁琐。 |
因此，提供了一个选项，可以在模拟期间独立于压力周期内的变化ETmax。
瞬态ETmax是以区域为单位应用的，并且MXZETS是由EVT或ET2包读取的，指示应用瞬态ETmax的最大区域数。
如果MXZETS为零，则模拟不会调用瞬态ETmax包。
然而，如果存在瞬态ETmax，EVT或ET2包将读取一个额外的整数数组，标识每个区域的面积位置。
值为零的区域表示该位置上不存在瞬态ETmax区域，并且将对该位置应用ETmax的默认值，如EVTR或ETSR数组中为压力周期读取的那样。
在值不为零的位置，EVTR或ETSR数组将使用在单独文件中提供的ETmax时间序列数据更新该区域的适当值。
此文件的默认扩展名为.ETS，其中包含为所有区域制表的ETmax时间序列。
文件的每一行包含起始时间（TSTART），结束时间（TEND），乘法因子（用于单位转换），然后是每个区域的一个列的ETmax速率。
如果后续的TSTART值与前一行的TEND值不匹配，则对于中间期间的每个区域都将应用零ETmax。
因此，不需要在.ETS文件中单独包含长时间不发生ET的行。
事件之间或事件内的时间段不需要固定，因为此模块设计用于自适应时间步进方案。
自适应时间步进方案确保模拟准确地命中每个TSTART和TEND，以确保准确地表示时间变化。
在周期内或周期间，自适应时间步进方案试图最大化计算效率。
EVT和ET2包的详细信息可以在相应的美国地质调查报告中找到，这些报告记录了这些包。
请注意，ET2包在原始MODFLOW文档中被标记为ETS1（Banta，2000）。
当前的文档描述了将这些包扩展为包括ETmax的时序输入。
G.2 蒸散蒸发（EVT）包输入
蒸散蒸发（EVT）包的输入是从IUNIT（5）中指定的单元中读取的。对于每个模拟
1.
| 数据： | NEVTOP | IEVTCB | MXZETS2 |
| 格式： | I10 | I10 | I10 |
对于每个压力周期
EVT1RP2
2.
| 数据： | INSURF | INEVTR | INEXDP | INIEVT |
| 格式： | I10 | I10 | I10 | I10 |
** 仅在INSURF≥0时输入项目3。
3.
| 数据： | SURF |
| 模块： | U2DREL |
** 仅在INEVTR≥0且MXZETS2=0时输入项目4a。
4a.
| 数据： | EVTR |
| 模块： | U2DREL |
** 仅在INEVTR≥0且MXZETS2…
4b.
| 数据： | IZNETS2 |
| 模块： | U2DREL |
请注意，根据是否请求最大ET通量（EVTR）的时间序列，将读取项目4a或4b。
** 仅在INEXDP≥0时输入项目5。
5.
| 数据： | EXDP |
| 模块： | U2DREL |
** 仅在NEVTOP=2且INIEVT≥0时输入项目6。
6.
| 数据 | IEVT |
| 模块 | U2DINT |
输入指令中使用的字段的说明
| NEVTOP | 是蒸散蒸发（ET）选项代码。 ET参数（ET表面，最大ET速率和消光深度）在二维数组SURF，EVTR和EXDP中指定，每个垂直列中有一个值。因此，ET计算了模型网格中每个垂直列中的一个单元。选项代码确定在列内为哪个单元将计算ET。 |
| 1 - 仅为顶部网格层的单元计算ET。 | |
| 2 - 用户在IEVT数组中指定每个垂直列的单元。 | |
| 3 - 为最高活动单元计算ET。 |
| IEVTCB -- 是标志和单元编号。 |
| 如果IEVTCB>0，则它是单元编号，ICBCFL（见输出控制）设置时将记录单元对单元的流量。 |
| 如果IEVTCB≤0，则不会打印或记录单元对单元的流量。 |
| MXZETS2 -- 是由EVT或ET2包读取的指示瞬态ETmax通量（EVTR）的区域时间序列的索引。 |
| =0如果不是，则按照压力周期读取最大ET通量（EVTR）。 |
| =n如果是，则使用单独的文件读取最大ET通量（EVTR）的时间序列，n是模拟中使用的总ET区域数。 |
注意：与ETmax的区域输入相关联的ET时间序列（ETS）将按照第G.4节中提供的说明读取。
| INSURF -- 是ET表面（SURF）的读取标志。 |
| 如果INSURF≥0，则将从第5项中读取包含ET表面高程的数组。 |
| 如果INSURF<0，则将重用前一压力周期的ET表面。 |
| INEVTR -- 是最大ET速率（EVTR）的读取标志。 |
| 如果INEVTR≥0，则将从第6项中读取包含最大ET速率的数组。 |
| 如果INEVTR<0，则将重用前一压力周期的最大ET速率。 |
| IZNETS2 -- 是定义每个节点的区域值的整数数组，用于ET时间序列的区域输入。 |
请注意，区域值为零的单元中应用的ET将为零。
| INEXDP -- 是消光深度（EXDP）的读取标志。 |
| 如果INEXDP≥0，则将从第7项中读取包含消光深度（EXDP）的数组。 |
| 如果INEXDP<0，则将重用前一压力周期的消光深度。 |
| INIEVT -- 是层指示器（IEVT）的读取标志。仅当ET选项（NEVTOP）等于二时使用。 |
| 如果INIEVT≥0，则将从第8项中读取包含层指示器（IEVT）的数组。 |
| 如果INIEVT<0，则将重用前一压力周期中使用的层指示器。 |
| SURF -- 是ET表面的高程（L）。 |
| EVTR -- 是最大ET速率（每单位面积的体积流速，LT-1）。 |
| EXDP -- 是ET消光深度（L）。 |
| IEVT -- 是层指示器数组。对于每个水平位置，它指示从中移除ET的层。仅当ET选项等于二时才需要它。 |
G.3. 蒸散蒸发分段（ET2）包的输入
| ET2包的输入是从名称文件中类型为“ET2”的文件中读取的，或从IUNIT（41）中指定的单元中读取的。 |
| 如果在基本包输入文件中指定了选项“FREE”，则所有单一值变量都是自由格式；否则，变量具有10个字符的字段。 |
| 由数组读取，由U2DREL或U2DINT（Harbaugh等人，2000）指示。在可选变量中显示在方括号中。 |
对于每个模拟
| 0. #文本 |
| 0号项目是可选的，'#'必须在第1列中。项目0可以重复多次。 |
| 1. NETSOP IETSCB NPETS NETSEG MXZETS1 |
| 对于每个压力周期 |
| 2. INETSS INETSR INETSX [INIETS [INSGDF]]] |
| 3. ETSS（NCOL，NROW）- U2DREL - 如果INETSS≥0 |
| 4a. ESTR（NCOL，NROW）- U2DREL - 如果INETSR≥0且MXZETS1=0 |
| 4b. IZNETS1（NCOL，NROW）- U2DREL - 如果INETSR≥0且MXZETS1=0 |
| 请注意，根据是否请求最大ET通量（ETSR）的时间序列，将读取项目4a或4b。 |
| 5. ETSX（NCOL，NROW）- U2DREL - 如果INETSX≥0 |
| 6. IETS（NCOL，NROW）- U2DINT - 如果NETSOP=2且INIETS≥0 |
| 7. PXDP（NCOL，NROW）- U2DREL - 如果NETSEG>1且INSGDF≥0 |
| 8. PETM（NCOL，NROW）- U2DREL - 如果NETSEG>1且INSGDF≥0 |
如果NETSEG> 1，则为上部段交点读取（NETSEG - 1）次项目7和8。
如果NETSEG> 2，则首先为上部段交点读取项目7和8，然后连续的重复项目7和8，以成功较低的交点。
G.4 瞬态潜在蒸散蒸发时间序列文件
当应用瞬态ET条件时，瞬态ETmax时间序列文件与EVT或ET2包一起使用，当它们各自的输入中的MXZETS标志不为零时。
此文件包含每个MXZETS区域的ETmax的时间序列信息。
请注意，区域划分在相应的EVT或ET2输入文件中定义。
在应用瞬态ET条件时，瞬态潜在蒸散蒸发时间序列文件在FORTRAN单元97上打开。
该文件采用自由格式，每个模拟内的每个ET周期都包含一行。
文件的默认文件名是fname.ETS，其中fname是模拟的文件名。
每行的数据条目如下：
| 数据 | TSTARTE TENDE TMULTE EZ(1) EZ(2) ...... EZ(MXZETS) |
| 格式 | 自由 |
输入中使用的字段的说明
| TSTARTE - 该事件的时间输入记录的开始时间 |
| TENDE - 该事件的时间输入记录的结束时间 |
| 请注意，TSTARTE和TENDE是相对于总模拟时间的，而不是相对于压力周期时间的，当模拟中发生多个压力周期时。 |
| TMULTE - 用于给定事件的所有区域的ET速率的乘法因子。此缩放对于单位转换很方便。 |
| EZ(1) - 该事件的区域1中的ET速率 |
| EZ(2) - 该事件的区域2中的ET速率 |
| ••• |
| EZ(MXZETS) - 该事件的区域MXZETS中的ET速率 |
一个示例的ETS文件如表G.1所示。
在此示例中，域被划分为8个ET区域。
在模拟的前2.375天内不发生ET，因此模拟将快速进行，根据收敛行为调整其时间步长。
在2.375天时，ET发生在一些区域，为期一小时，直到2.417天。
在2.417天时，ET的强度（和/或区域）发生变化，这些条目出现在文件的第二行中，依此类推。
请注意，在2.5到4天之间不发生ET，不需要在文件中为此情况添加条目。
fname.ETS样本文件：
# 起始时间 结束时间 乘法因子 区域1 区域2 区域3 区域4 区域5 区域6 区域7 区域8
0.0 2.375 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
2.375 2.417 1.0 0.1 0.2 0.0 0.3 0.0 0.0 0.4 0.0
2.417 3.0 2.0 0.0 0.5 0.0 0.6 0.0 0.0 0.0 0.7
3.0 4.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0在这个示例中，第一行表示模拟开始后的2.375天内没有ET。
第二行表示从2.375天到2.417天，ET发生在区域1、2、4和8，各自的速率由乘法因子调整。
第三行表示从2.417天到3.0天，ET的强度和/或区域发生了变化，每个区域的速率由乘法因子调整。第四行表示从3.0天到4.0天，没有ET发生。
请注意，这是一个示例文件，实际的ETS文件可能根据模型的要求和ETmax的时间变化进行调整。
扫一扫
1594
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
