MIKE-SHE模型运行误差:“UZ柱短于初始水深”

Question

我正在尝试使用MIKE-SHE(2020版，DHI)运行一个集水区模型，其中包括一个不饱和区(UZ)和一个饱和区(SZ)。我已经指定了地面流和SZ的“有限差分”方法，以及UZ的理查兹方程。我已经输入了UZ和SZ材料的属性数据(目前沙子的默认值用于UZ);然而，我得到以下错误，然后是终止:

误差:UZ柱短于(7 146)的初始水深到地下水的深度:0.90 m UZ柱长度:0.90 m(剖面号)。6)终止……呼叫MIKE 1D在0 timesteps后停止请求…Ok异常终止

看起来错误是由于单个网格单元格(7146)造成的，但我不知道错误的原因，也不知道我可以调整什么来解决它。

Answer 1

看起来你的模型中的UZ列在某些位置的设置中延伸得不够深。

在MIKE-SHE中，当使用UZ Richard方程和SZ有限差分耦合时，模型中的UZ ' soil '剖面和SZ material column需要重叠:使UZ剖面底部深度延伸下面深圳的顶部。此外，地下水位高度应始终高于UZ剖面的底部。虽然在SZ模块中模拟了地下水位，但如果地下水位下降到低于模型UZ剖面底部的高度，则存在一个有效的不饱和层下面模型的UZ计算域。这可能会导致模拟UZ剖面底部的压头计算出现误差。

你收到的模型误差显示地下水位位于UZ剖面的最底部，而不是在UZ剖面内，在模拟开始时指定的网格单元格(7146)处。这可能对其他细胞也适用，但在检测到第一个病例后，该模型就终止了。

应该定义UZ概要文件以进行扩展下面模拟中地下水位所能达到的最深深度。

也有人建议，最上面的计算层您的SZ柱的定义应使地下水位始终位于该SZ层顶部。这是因为模型中UZ和SZ之间的水交换仅计算了SZ顶层。正因为如此，为了安全起见，UZ剖面可以被定义为向下延伸到刚刚低于顶部SZ计算层的底部高程的深度。

为了解决您的问题，尝试进一步向下扩展您的UZ“土壤”剖面，以覆盖您期望的地下水位的全范围深度:调整“土壤剖面”表和计算单元“垂直离散化”表的相关土壤/材料剖面。注意:MIKE-SHE界面将其标记为“土壤剖面”，但它可能需要扩展到通常定义为“土壤”的范围之外，进入通常不饱和的风化岩层。

您可能还需要调整顶部SZ计算层的厚度(或顶部和底部的高度)，以便在整个模拟过程中水位也位于顶部层内。如果预测地下水位会大幅上下移动，那么在整个模拟过程中获得正确的深度可能是一个迭代过程。