基于多维度的暴雨洪水相似过程判定方法与流程

本发明涉及水文预报领域，特别涉及一种基于多维度的暴雨洪水相似过程判定方法。

背景技术：

1、洪水预报预测是防洪工程措施的重要组成部分。目前，基于暴雨洪水知识相似性分析的预报方法有很多，如主成分分析法、灰色关联分析法、相关系数法和加权距离系数法等。尽管大量关于降雨相似分析模型的研究已经开展，但目前相似性方法只是根据某一维的相似进行判断，不够全面。为此，本文提出了一种基于时间、空间维度和水文状态综合的暴雨洪水相似过程判定方法，它具有从多个维度考虑相似性的优点，极大的提高相似性查找的准确性。

技术实现思路

1、本发明所要解决的技术问题是提供一种基于多维度的暴雨洪水相似过程判定方法，考虑暴雨洪水多个维度的水文特征，提高相似查找的准确率。

2、为解决上述技术问题，本发明所采用的技术方案是：一种基于多维度的暴雨洪水相似过程判定方法，包括如下步骤：

3、s1、收集已知流域的暴雨洪水数据，并根据水文学原理划分场次；

4、s2、划分出数据集中的历史查找样本和测试样本；

5、s3、构建水文状态相似计算模型，获取历史查找样本的水文状态相似度；

6、s4、以水文状态相似为基准，构建孪生网络提取图像特征开展时间维度和空间维度相似分析，得到历史查找样本在时间维度和空间维度的相似度；

7、s5、基于lasso回归进行多维度相似度的线性拟合，计算水文状态相似度、时间维度相似度和空间维度相似度的权重参数；

8、s6、基于拟合的权重参数，融合水文状态相似、时间维度相似和空间维度相似构建场次降雨的多维度相似过程查找模型；

9、s7、查找历史相似暴雨洪水场次。

10、优选的方案中，所述步骤s1包括如下步骤：

11、s101、收集已知流域的场次降雨径流数据；

12、s102、以统一降雨径流的时间步长为要求划分场次，划分场次要保证大、中、小场次过程尽量分布均值。

13、优选的方案中，所述步骤s2中，按照一定比例将数据集划分为历史查找样本和测试样本。

14、优选的方案中，所述s3包括如下步骤：

15、s301、随机在测试样本中指定一个样本；

16、s302、将指定的测试样本输入到程序，使用归一化将各项特征数据转换到一定数值范围内，假设xji为某一项特征值，则归一化后得到的数值yji如下式所示：

17、

18、其中，j表示某一项水文特征，i表示某一项特征的某一个数值，xjmin表示在j特征列中的最小值，xjmax表示在j特征列中的最大值；

19、s303、将指定的测试样本c与每一个历史查找样本l的同一特征的值进行相减，返回计算结果绝对值如下式所示：

20、

21、s304、计算每项水文特征与洪峰的相关系数wj，逐一将历史查找样本每项水文特征的差值乘以对应相关系数并求和，得到测试样本和每一个历史查找样本的水文状态相似度sl：

22、

23、s305、将s304中历史查找样本按相似度从小到大排序，选取最相似的数个样本。

24、优选的方案中，所述s4中包括如下步骤：

25、s401、将s305中查找的结果场次样本作为时间、空间维度的历史相似基准，再查找其与测试样本的时间维度和空间维度相似度；

26、s402、以降雨量为横坐标轴，降雨时间为纵坐标绘出各个场次的柱状图；

27、s403、提取降雨观测区域的经纬度和其观测的降雨量，绘出各个场次的色斑图；

28、s404、对柱状图和色斑图分别进行裁剪，对裁剪的图像进行数据增强，将图像数据集扩展；

29、s405、以resnet152模型为图像特征提取器，其中每个残差块表示为：

30、yi＝h(xi)+f(xi,ωi)

31、xi+1＝f(yi)

32、其中，f为残差函数，f为relu函数，ωi为权值矩阵，xi和yi分别为第i层的输入和输出，函数h(xi)表示为：

33、h(xi)＝xi

34、残差函数f(xi,ωi)的定义为：

35、f(xi,ωi)＝ωi·σ(b(ω′i)·σ(b(xi)))

36、其中，b(xi)是批处理规范化，“·”表示卷积，σ函数表示的含义为：σ(x)＝max(x,0)；

37、s406、构建resnet152模型并迁移学习在imagenet数据集上已经训练好的权重，在此基础上分别使用s404中的柱状图和色斑图数据集进行训练并更新保存模型参数，分别作为查找时间维度和空间维度相似度的特征提取权重参数；

38、s407、建立孪生网络相似性分析框架，分别载入s406中的模型参数，输入历史查找样本和测试样本对应柱状图或色斑图，计算由resnet152模型提取的特征向量间欧氏距离，作为s401中历史相似基准的历史查找样本和测试样本在时间维度和空间维度的相似度，其中欧氏距离dw的计算方法为：

39、

40、其中，dw代表两个样本特征x1和x2的欧式距离，p表示样本的特征维数；

41、优选的方案中，所述s5包括如下步骤：

42、s501、将s2划分的历史查找样本中随机选取部分样本作为线性拟合的训练样本，未选样本作为训练样本的历史查找样本；

43、s502、计算训练样本与历史查找样本的洪峰参数绝对差值并作为因变量y，历史查找样本对应训练样本的水文状态相似度、时间维度相似度和空间维度相似度作为自变量(x1,x2,x3)，构建基于lasso回归的多元线性方程，获取多维度相似权重参数，具体方程为：

44、y＝w1x1+w2x2+w3x3+b

45、其中，b表示回归常数，w1为水文状态相似的权重参数，w2为时间维度相似的权重参数，w3为空间维度相似的权重参数；

46、优选的方案中，所述步骤s6中，基于s502得到的权重参数，将多维相似度与对应的权重系数相乘再求和，得到融合相似度sr，构建场次降雨的多维度相似过程查找模型，其中sr计算方法为：

47、sr＝s1·w1+s2·w2+s3·w3

48、其中，s1为水文状态相似度，s2为时间维度相似度，s3为空间维度相似度。

49、优选的方案中，所述步骤s7中，按照融合相似度大小依次排列历史查找样本，得到融合相似度最小的样本，作为最相似历史暴雨洪水场次。

50、本发明提供的一种基于多维度的暴雨洪水相似过程判定方法，具有以下有益效果：

51、1、本发明创造性的基于状态之间的相关性来挖掘场次降雨间的水文相似性，可以更加高效查找水文状态特征最相似的历史场次，查找过程具有可解释性高的优势。

52、2、本发明利用高层数的resnet152模型作为图像特征提取网络，基于resnet152模型并构建孪生网络来计算具备场次降雨时间和空间维度特征的图像相似度，深度挖掘每个历史场次与预测场次之间存在的时间维度特征和空间维度特征的相似度。

53、3、本发明基于lasso回归多元线性拟合的方法，可以在权衡模型拟合数据和模型复杂度之间找到一个平衡点，生成一个更简化的模型，并将多维相似度作为自变量，对应场的洪峰预测绝对误差作为因变量，有效的分析出多维相似度在预测精准性的权重分布，得到多维相似度各自的权重参数。

54、4、本发明更全面考虑的暴雨洪水的多个水文特征，基于多元线性回归融合场次降雨的水文状态特征、时间维度特征和空间维度特征，查找在融合相似度上最相似历史相似暴雨洪水场次，可以有效的提高相似查找的准确率。