3.10 步骤 10 — 重新定义水平网格的离散
MIKE SHE的优点之一是您可以随意改变网格的尺寸而不需要重新定义模型参数。该优点便于您对模型进行网格离散的灵敏度分析,这在其他模型中很难实现。
关于Demo 版的说明:Demo版最多支持70x70个网格,所以您无法在Demo版本中运行前处理器或者运行网格加密的模型。 3.10.1 改变水平网格的分辨率
在Model Domain and Grid对话框 ? ? ? ?
双倍加密X 和Y 方向的网格 NX = 130 NY = 140
将Cell Size 定为 250 m
图3.30
模型网格设置界面
3.10.2 运行前处理器
? 点击,启动前处理器
查看前处理数据,确认所有的输入数据都转化为250m网格的数据 将网格尺寸改为400m: 在Model Domain and Grid 对话框 ? ? ? ? ?
改变X 和Y 方向的网格数目 NX = 100 NY = 100
将Cell Size 设为 400 m 运行前处理器
2-29
MIKE SHE分布式水文模型培训程? DHI
图3.31
模型网格尺寸设为400m
将网格尺寸恢复为500m: 在Model Domain and Grid 对话框 ? ? ? ?
改变X 和Y方向的网格数目 NX = 65 NY = 70
将Cell Size 设为 500 m
图3.32
模型网格尺寸设为500m
3.10.3 保存SHE文件
3.11 步骤 11 — 定义输出和率定项
Detailed Time Series Output可以输出模型范围内某点的时间序列图。它还可以进行模型误差统计以及模拟结果和实测纪录的比较。
详细时间序列以HTML的形式输出,可将该HTML文件并入网页以查看率定的进展。 3.11.1 结果储存
? 在Storing of Water balance data一栏打勾
? 在Storing interval for grid series output一栏,将Prec(ipitation),SZ-heads 和 SZ-fluxes的储存频率分别设为24小时
2-30 MIKE SHE分布式水文模型培训程? DHI
图3.33
结果存储设置界面
3.11.2 定义详细时间序列的输出
? ? ? ? ?
点击按键,在表中增加4行 对于每一项, 命名
选择 head elevation in saturated zone 作为 Data Type 各点的X,Y坐标如下表:
Name Obs 5 Obs 35 Obs 37 Obs 65 X-coordinate 24103 13800 17508 8866 Y-coordinate 8419 17607 20887 25549 ? 将Depth 设为 5
? 确认没有选上Incl. Obs. Data
2-31
MIKE SHE分布式水文模型培训程? DHI
图3.34
时间序列输出设置界面
3.11.3 网格输出
如果您在storing of Water Balance data一栏打勾,那么Grid series output对话框的所有water balance项都会自动被选上。
? 如果您没有在storing of Water Balance data一栏打勾,那么您应该选择以下项 precipitation rate,
depth to phreatic surface,
head-elevation in saturated zone,
groundwater flow in the X-direction, and groundwater flow in the Y-direction.
Note:饱和带的输出(即存储)时间步长必须是最大运算时间步长的整数倍,例
如最大运算时间步长为24小时,则输出的时间步长可以是24、48、72小时,以此类推。恒定流模拟同理。
图3.35
网格输出设置界面
2-32 MIKE SHE分布式水文模型培训程? DHI
3.12 步骤 12 — 运行模拟
在运行模拟之前,您应该返回到Setup菜单,对模型设定作一些修改,包括模拟时间长度、时间步长和结果储存。 3.12.1 定义时间步长参数
在Time step control 对话框:
? 将Initial time step设为24 hrs
? 将Max allow SZ time step设为 24 hrs
? 将Max infiltration amount per time step设为10 不改变其他的默认参数。
图3.36
时间步长参数设置界面
3.12.2 模拟控制
在SZ Computational Control Parameter 对话框:
? Solver Type选择Preconditioned Conjugate Gradient,Transient ? 将Maximum Head change per iteration 设为 0.0001 [m] 模拟收敛参数的默认值一般适用于大多数情况而无需更改。
2-33
MIKE SHE分布式水文模型培训程? DHI