一种基于流量重现期的流域山洪风险分布式表达方法

1.一种基于流量重现期的流域山洪风险分布式表达方法，其特征在于：包括以下步骤：S1、构建流域分布式新安江模型，得到分布式网格流量序列；S2、基于数字高程模型DEM进行水文分析计算，得到流域汇流累积量的空间分布，对所述流域汇流累积量进行归一化处理，生产空间个体网格对流域出口流量贡献度空间分布；将流域汇流累积量进行归一化处理，具体为：计算网格(i,j)的归一化流域汇流累积量，即流域出口流量贡献度ACC i,j ，无量纲；i,j分别代表计算网格所在行数和列数；若ACC max ＝ACC min ，则ACC i,j ＝1；若ACC max ≠ACC min ，则ACC i,j ＝(ACC i,j -ACC min )/(ACC max -ACC min )；遍历流域范围内的所有网格，得到空间个体网格对流域出口流量贡献度空间分布ACC con ；其中，m为研究区网格化后的行数；n为研究区网格化后的列数；S3、将流域出口控制站的设计流量值与流域出口流量贡献度空间分布进行空间矩阵计算，得到流域空间网格设计流量；根据流域出口流量贡献度空间分布和流域出口控制站设计流量值，通过矩阵运算，求得各个流域空间网格设计流量QGrid，单位：立方米/秒；其中，Q P＝0.01 、Q P＝0.02 、Q P＝0.05 、Q P＝0.1 表示出口控制站100年一遇、50年一遇、20年一遇、10年一遇的设计流量值；单位：立方米/秒；S4、将分布式模型计算所得的分布式网格流量序列与流域空间网格设计流量进行差分、叠置、融合计算，得到流域山洪风险时空图分布图集；S401、设置开始时间和结束时间，输入驱动数据，所述驱动数据包括：降雨、蒸发能力，运行已率定好的分布式新安江模型，得到分布式网格流量序列：其中，Qt为t时刻流域内的分布式网格流量序列，单位：立方米/秒；S402、将某时刻分布式网格流量序列Qt与各重现期的流域空间网格设计流量QGrid进行差分、叠置、融合计算，得到网格差分流量Qt dis ：其中，Qt dis 为t时刻流域内的网络差分流量，单位：立方米/秒；S403、对所述网络差分流量进行条件运算，得到各网格是否超过给定设计流量标准的矩阵，即流域山洪风险时空图集Rt，无量纲；若Qt dis ＞＝0,Rt＝1；有给定重现期的山洪风险；若Qt dis ＜0,Rt＝0；无给定重现期的山洪风险；S404、重复S402和S403，获得时序风险矩阵R：R＝[R1,R2,...,Rn]；S405、对所述时序风险矩阵R进行颜色渲染，生成流域山洪风险时空间分布图集。 2.根据权利要求1所述的一种基于流量重现期的流域山洪风险分布式表达方法，其特征在于，所述S3的具体步骤为：S301、获取流域出口控制站的设计流量值；若流域出口控制站已有设计流量成果，直接使用；若流域出口控制站无已有设计流量成果，存在长时间序列流量观测值，采用P-III频率曲线进行频率分析，获得相应频率的流域出口控制站设计流量值；若流域出口控制站无设计流量成果，同时也无长时间序列流量观测值，采用相似流域缩放的方式获得相应频率的流域出口控制站设计流量值。