[发明专利]山洪风险预测方法、装置、设备及计算机存储介质

权利要求书

1.一种山洪风险预测方法，其特征在于，所述方法包括：

根据数字高程模型将目标流域划分成若干个均匀的网格单元；

根据新安江模型计算每个所述网格单元的蒸散发量和每个所述网格单元的产流量，和每个所述网格单元的地面径流，和每个所述网格单元的壤中流，和每个所述网格单元的地下径流；

根据有限体积法求解每个所述网格单元的二维浅水方程，以计算出每个所述网格单元中水的深度，和每个所述网格单元中的水的速度，和每个所述网格单元与相邻网格单元的水量交换值；

基于所述目标流域的预设雨量，根据每个所述网格单元的蒸散发量，和所述产流量，和所述地面径流，和所述壤中流，和所述地下径流，和所述水的深度，和所述水量交换值构建所述目标流域的山洪过程，以对所述山洪风险进行预测。

2.如权利要求1所述的方法，其特征在于，所述根据数字高程模型将目标流域划分成若干个均匀的网格单元，包括：

读入所述数字高程模型的数据文件，其中，所述数据文件包括空间节点的数据；

以四个所述空间节点为单位，根据所述空间节点的数据将所述目标流域划分为若干个均匀的网格单元。

3.如权利要求1所述的方法，其特征在于，所述根据新安江模型计算每个所述网格单元的蒸散发量，包括：

将每个所述网格单元的实测水面蒸发值作为所述新安江模型的输入，以使得所述新安江模型根据三层蒸散发模式和所述实测水面蒸发值输出每个所述网格单元的上层蓄水量和每个所述网格单元的下层蓄水量；

根据每个所述网格单元的所述上层蓄水量和每个所述网格单元的下层蓄水量计算每个所述网格单元的三个土层的流蒸散发量；

根据每个所述网格单元的三个土层的流蒸散发量计算每个所述网格单元的总流蒸散发量。

4.如权利要求3所述的方法，其特征在于，所述将每个所述网格单元的实测水面蒸发值作为所述新安江模型的输入，以使得所述新安江模型根据三层蒸散发模式和所述实测水面蒸发值输出每个所述网格单元的上层蓄水量和每个所述网格单元的下层蓄水量，包括：

当P+WU≥EP时，EU=EP,EL=0,ED=0；

当P+WUEP,WL≥C×WLM时，

EU=WU+P,EL=(EP-EU)×WL/WLM,ED=0；

当P+WUEP,C×(EP-EU)≤WLWLM时，

EU=WU+P,EL=C×(EP-EU),ED=0；

当P+WUEP,WLC×(EP-EU)时，

EU=WU+P,EL=WL,ED=C×(EP-EU)-EL；

其中，WU表示所述上层蓄水量；WL表示所述下层蓄水量，EP表示预设的流域蒸散发能力，K表示蒸散发折算系数、EM表示所述实测水面蒸发值；WLM表示下层的张力水蓄水容量、C表示深层蒸散发数、所述，P表示降雨量,所述EU，所述EL，所述ED分别表示所述三层蒸散发模式中的三个土层的流域散发量。

5.如权利要求1所述的方法，其特征在于，所述根据新安江模型计算每个所述网格单元的产流量的计算式为

当时，流域局部产流，相应的流域总产流量R为：

当，即全流域产流时：

其中，WMM表示流域单点最大张力水容量、WM表示流域平均张开力水容量、PE表示扣除雨期内蒸散发后的降雨量、B表示张力水容量曲线的方次，W0表示流域初始土壤含水量，W0=WM，。