本发明涉及洪涝灾害预测,具体地说是一种基于降雨的洪水淹没预测方法。
背景技术：
1、洪水淹没的机理是水源区和被淹没区有通道（如溃口和开闸放水等）并存在水位差，产生淹没是一个动态而致平衡的过程。然而高强度大范围的降雨导致低洼处积水的过程更趋近于无源淹没过程。现有技术中尚未有基于降雨等无源淹没过程导致的洪水淹没预测的方法。
2、因此,本技术提出了一种基于降雨的洪水淹没预测方法，利用洪水淹没区水位-水量关系，对淹没面积进行快速估测，能够在确定靶区的基础上，快速确定监测范围，为后续的洪水监测数据的下载、灾情信息的提取提供更细致的范围圈定依据。
技术实现思路
1、本发明的目的是克服现有技术的不足,提供了一种基于降雨的洪水淹没预测方法，利用洪水淹没区水位-水量关系，对淹没面积进行快速估测，能够在确定靶区的基础上，快速确定监测范围，为后续的洪水监测数据的下载、灾情信息的提取提供更细致的范围圈定依据。
2、为了达到上述目的,本发明提供一种基于降雨的洪水淹没预测方法，包括：
3、1.1湖水体积变化分析；
4、1.2历史雨量累积分析；
5、1.3湖水体积变化预测模型；
6、1.4淹没范围预测；
7、1.1湖水体积变化分析的具体步骤为：
8、s1，淹没分析：
9、有源淹没：
10、s1-1，新建空的淹没图：根据输入的dem栅格文件，生成具有相同栅格信息的空文件作为淹没图文件，淹没图栅格文件与输入dem栅格文件具有相同的行列值、四角坐标、栅格大小、坐标系，所有栅格值均为空；
11、s1-2，处理输入源点：
12、s1-2-1，根据输入的源点坐标，计算其在dem上的对应栅格，取其行号r，列号c，高程值e；
13、s1-2-2，新建源点列表，用于存储源点，源点的类型为元祖，即源点列表写作[(r，c),……]，其中r为行号，c为列号；
14、s1-2-3，将源点栅格的行列号作为第一个源点存入源点列表中；
15、s1-2-4，定义当前淹没高程e，将其赋值为e，作为淹没第一步的淹没高程；
16、s1-3，以e为淹没高程，根据源点列表，完成e水位下的有源淹没扩展蔓延算法，更新淹没图：
17、s1-3-1，新建源点删除列表：用于存放不再位于边界的，需要从源点列表中删除的源点信息；
18、s1-3-2，遍历源点列表：对于当前源点，取源点的8个邻域栅格，形成邻域列表；设置边界标志，并初始化为false；
19、s1-3-3，遍历邻域列表，对于当前邻域点，判断邻域点在淹没图中是否值为空；若值不为空，则代表该点已经淹没，不处理当前邻域点；若值为空，则进行下一步判断；判断邻域点的高程值e’，如果e’＞e，若该邻域点未被淹没，不处理该邻域点，将当前源点的边界标志设置为ture；若e’≤e则该点被淹没，将该邻域点加入源点列表，更新该点在淹没图中的栅格值为e；
20、s1-3-4，遍历所有邻域点后，如果边界点标志仍为false，说明该源点的所有邻域皆已被淹没，不再位于边界，将该源点加入源点删除列表中；
21、s1-3-5，遍历所有源点之后，从源点列表中将记录在源点删除列表中的源点删除；淹没高程增加一个步长，即e=e+estep；若源点列表不为空则重复步骤s1-3开始；否则淹没图计算完成，退出程序；
22、s1-3-6，通过淹没图算法，获取“水位-淹没面积关系”与“水位-淹没体积关系”，“水位-淹没体积关系”用于将降雨量与淹没范围联系起来的换算依据；无源淹没：
23、通过统计水位以下的dem体积，即可得到s1-3-6中的“水位-淹没体积关系”；
24、s2，实测水位数据：
25、通过收集现有的相关水文数据中各站点的水位数据，并根据不同站点间的最接近水位数据作为建模水位；
26、s3，水位-体积换算模型：
27、s3-1，采用dem模型对水位-淹没体积关系进行统计，统计面积覆盖整个水域，结合历年水位变化范围和高水汛期中的某段数据，作为换算的依据；
28、s3-2，对水位-体积统计结果进行拟合：
29、对s3-1中的结果进行拟合，通过水位-体积直接换算方法用于保留地形影响下计算目标水位e的体积，计算公式为：
30、公式一：
31、其中，n为水位e的整数部分，为查询到的水位为n时对应的淹没体积，为查询到的水位为n+1时对应的淹没体积；
32、s4，水位体积变化量：
33、基于dem的水位-淹没体积换算、计算获得水位体积变化量，针对每日的水位实测数据，通过换算获得每日的湖水体积，第n日水位体积的变化量为，其计算公式为：
34、公式二：；水位体积变量是未来24小时后水体积的变化量；
35、1.2历史雨量累积分析的具体步骤为：
36、s10：雨量数据收集：
37、通过era5收集历年陆地小时雨量数据集；
38、s20：雨量数据规整：
39、将小时降雨量规整为24小时累积降雨量，对雨量数据按距离进行规整并裁剪掉流域以外的数据；
40、s30：历史累计雨量：
41、以历史水位最高点作为统计对象，降雨数据时间以0时~24时累积，空间从主体范围按照覆盖目标区的预设的半径作为起点，按照预设的步长，计算主体中心点的环状区域的单日累积雨量数据；
42、1.3湖水体积变化预测模型的具体步骤为：
43、s100，预测模型变量选择：
44、基于降雨量数据，通过构建实测水位变化与积累雨量的统计关系模型来综合考虑流出量、蒸发、地下水交换因素的整体影响；
45、s200，变量相关性分析：
46、s200-1，通过对环状区域进行累加，构造圆形区域降雨量累计变量，进行双环状范围和圆形范围的双变量分析；
47、s200-2，考虑到时间因素，对水位水体积变化量和不同距离的累积雨量进行时序特征分析；
48、s200-3，基于圆形范围累积雨量，对比水位体积变化量与不同范围的累计雨量之间的相关系数，用于进一步确定雨量累积空间范围；
49、s300，模型初步拟合：
50、s300-1，对不同距离范围的降雨量数据进行多轮主成分分析；
51、s300-2，对各轮主成分合成的累积雨量与水位体积变量进行拟合，并与圆形区域累积雨量与水位体积变量的拟合结果进行对比；
52、s400，样本筛选：通过水位体积变化量的筛选，用以排除掉以泄出、增散发为主导因素造成的水体变化关系样本；
53、对全年样本通过水位不同进行多组筛选，并对筛选后的样本进行相关性分析；
54、对不同筛选数据集中，水位体积变化量与不同范围的雨量累积量之间的相关性与拟合进行计算；
55、s500，降雨累积时长对比：分析对比不同累积时长的累积雨量与水位体积变化量的相关性，确定模型中使用的降雨量累计时长；
56、s600，模型及拟合效果：对比仅水位筛选的拟合情况和水位＋体积变化量筛选的拟合结果；
57、s600-1，对比线性、二阶多项式、三阶多项式拟合效果；
58、s600-2，对比非线性的神经网络方法拟合效果；
59、s600-3，使用经过水位筛选的24小时降雨累积时长的、水位体积变化量10°以上的样本子集进行拟合建模，计算公式为：
60、公式三：；
61、其中为水位体积变化量，单位为m3；p为流域内圆形范围内的24小时累计雨量，单位为m3；c和b分别为常数项和系数；
62、1.4淹没范围预测的具体步骤为：
63、s1000，未来24小时累计雨量计算：
64、s1000-1，下载栅格格式的未来24小时预报雨量数据；
65、s1000-2，制作雨量累积范围文件，提取流域内部分，获得预报雨量的累积范围文件；
66、s1000-3，对累积雨量进行投影，投影结果水平方向单位为m；接着是要用雨量累积范围文件对投影后的雨量文件进行提取，提取出累积范围内的雨量；最后对提取出的雨量文件进行表面体积计算；
67、s2000，水位预测：
68、s2000-1，基于累积雨量计算水体体积变化量，将24小时预报降雨累积值代入公式三，计算获得预测水位体积变化量；
69、s2000-2，获取当前水位数据测值e，将当前水位值代入公式一获得当前水位对应的湖水体积，并根据预测的水位体积变化量获得预测湖水体积，公式为：
70、公式四：；
71、对公式一进行改写，获得对应的预测水位，公式为：
72、公式五：；
73、其中，为当水位为整数n时对应的湖水体积，为当水位为整数n+1时对应的湖水体积，在计算式，首先找到水位n，使＜＜；
74、s3000，淹没范围预测：
75、对制作好的淹没图进行插值处理，使其与实测水位数据精度一致。
76、s10中era5的平面分辨率为0.1°，换算后为10km。
77、s300-1中的多轮主成分分析包括四轮：
78、第一轮以180km以内环形区域累积降雨量为因子；
79、第二轮以全部距离的圆形区域累积雨量为因子；
80、第三轮以90～180km之间的圆形区域累积雨量为因子；
81、第四轮以120～180km之间的圆形区域累积雨量为因子。
82、s600-3中的c和b分别为2062612888.591和0.100611。
83、本发明同现有技术相比，基于dem数据进行淹没分析，统计水位-淹没面积与水位-淹没体积关系，构建淹没面积随水位变化的洪水淹没预测模型，在此基础上利用水量总体积相等原理，发展基于总降雨量数据的淹没范围快速预测方法，进一步圈定洪水监测区范围。
1.一种基于降雨的洪水淹没预测方法，其特征在于，包括：
2.根据权利要求1所述的基于降雨的洪水淹没预测方法，其特征在于，所述s10中era5的平面分辨率为0.1°，换算后为10km。
3.根据权利要求1所述的基于降雨的洪水淹没预测方法，其特征在于，所述s300-1中的多轮主成分分析包括四轮：
4.根据权利要求1所述的基于降雨的洪水淹没预测方法，其特征在于，所述s600-3中的c和b分别为2062612888.591和0.100611。
