一种考虑不同洪水阶段的水文模型参数敏感性分析方法与流程

本发明涉及水文，尤其是涉及一种考虑不同洪水阶段的水文模型参数敏感性分析方法。

背景技术：

1、近年来，受气候变化影响，由局地强降水造成的突发性洪水频繁发生，已成为造成人员伤亡的主要灾种。据统计，我国河流洪涝灾害和山洪地质灾害损失约占全国洪涝灾害经济损失的70％～80％，死亡人数占2/3左右，对人民群众生命财产安全构成了严重威胁，因此对流域突发性洪水进行准确可靠地预报成为亟待解决的重要问题。

2、水文模型是认识水文科学规律、分析水文过程及研究水文循环机理的重要科学工具。水文模型模拟结果的不确定分析是提高模型可靠性、进行有效水情预报的一个重要研究内容。参数不确定性是影响水文模型模拟结果不确定性的关键因素之一，开展模型参数不确定性及其影响因素分析对水文预报具有重要现实意义。要进行可靠的洪水预报，客观全面地分析了解洪水的特征，了解不同阶段的洪水过程中模型参数的敏感性是首要的工作。在长期进行的洪水预报研究工作中，通常将一场洪水当成一个整体过程，并采用水文模型对这些场次洪水进行模拟，计算确定性系数，洪峰误差，洪量误差等指标评价模拟结果，进而对模型参数进行率定调整，得到能够合理模拟洪水的水文模型。然而，在这个过程中通常将整个场次洪水与模拟结果相对比，忽略了洪水具有涨洪，洪峰，落洪等不同阶段，从而不利于客观全面地认识和了解在不同阶段洪水阶段中模型参数的敏感性，不利于洪水过程的精细化模拟。

3、因此，实有必要设计一种考虑不同洪水阶段的水文模型参数敏感性分析方法，以克服上述问题。

技术实现思路

1、为了避免上述问题，提供了一种考虑不同洪水阶段的水文模型参数敏感性分析方法，数据来源稳定可靠，变量之间的函数关系明确，计算效率高、结果客观合理，有利于不同洪水阶段的水文模型参数敏感性的快速合理分析，可以进一步促进数字水文学以及流域洪水防治研究的深入发展，值得推广。

2、本发明提供的一种考虑不同洪水阶段的水文模型参数敏感性分析方法，包括如下步骤：

3、步骤1，将时间t作为自变量，河道断面的流量q作为因变量，计算q随t的变化率q′；

4、步骤2，对q随t的变化率q′进行滤波，得到变化趋势线q′sth；

5、步骤3，根据变化趋势线q′sth，将洪水过程划分为涨前阶段、上涨阶段、包含洪峰阶段、持续消退阶段和回归稳定阶段；

6、步骤4，设置参数值域，在不同阶段内，按照一定的步长使水文模型参数由大到小逐渐变化，模拟得到洪水过程，并与实测洪水比较，计算不同洪水阶段的确定性系数方差σ作为评价指标，量化分析不同洪水阶段的水文模型参数敏感性。

7、优选地，步骤1中，将洪水过程河道断面的流量q视作时间t的函数q＝f(t)，计算函数q的一阶导数

8、优选地，步骤2具体包括如下子步骤：

9、2.1设置滤波窗口，其宽度为2w+1个数据，窗口内的数据点为s(t)，t的取值为[-w,…0,…,w]，构造一个n阶多项式拟合滤波窗口内的数据序列，多项式如下：

10、

11、式中：cn0、cn1、…、cnn为多项式的待定系数；k为多项式中单项的阶数，从0到n，n为多项式阶数；t为时间；f(t)为拟合的数据序列；

12、2.2计算拟合的数据序列f(t)与原数据点s(t)的残差平方和e：

13、

14、2.3令e对各系数的导数为0，使残差平方和最小，求解多项式的系数，采用拟合得到的多项式求解滤波窗口内中心点的估计值，然后移动滤波窗口，求解得到变化趋势线q′sth，具体公式如下：

15、

16、求解得到：

17、令则

18、对于待拟合的数据[s(-w),…,s(0),…,s(w)]，求出fr和pk+r，得到多项式的待定系数cn0、cn1、…、cnn，从而确定多项式；其中，每当移动窗口时，该多项式在窗口中心的值就是滤波后的输出结果。

19、优选地，步骤3具体包括如下子步骤：

20、3.1按照时间顺序，在变化趋势线q′sth上依次设置前方区间[xs,xs+σ]和后方区间[xs-σ,xs]，其区间的步长为σ，计算前方区间方差和后方区间方差之间的比值ωs，公式如下：

21、

22、式中：m为前方区间数据点的均值；xs+j为前方区间内的数据点，s+j为数据点下标，其中，s为在变化趋势线q′sth上的当前位置，j从0变化到σ；p为后方区间数据点的均值，ωs为前方区间方差和后方区间方差之间的比值；

23、3.2确定ωs最大值所在的位置，将其作为变化趋势线q′sth由平稳变上涨的转折点α；确定ωs最小值所在的位置，将其作为变化趋势线q′sth变回平稳的转折点β，并分别确定变化趋势线q′sth的最高点γ和最低点δ；

24、3.2将α点之前的部分作为涨前阶段，α点和γ点之间的部分作为上涨阶段，γ点和δ点之间的部分作为包含洪峰阶段，δ点和β点之间的部分作为持续消退阶段，β点之后的部分作为回归稳定阶段，实现洪水阶段划分。

25、优选地，步骤4具体包括如下子步骤：

26、4.1设置参数值域，在划分的不同洪水阶段中，按照一定的步长δp使水文模型参数由大到小逐渐变化，离散得到参数变化数组listpara：

27、listpara＝{δp,2δp,3δp…u×δp…m×δp}；

28、式中：δp为参数步长，u为步长下标，从1到m，m×δp为参数值域上限；

29、4.2将参数应用于水文模型中，得到洪水过程模拟结果，并与实测洪水比较，计算确定性系数dc，并统计不同洪水阶段的确定性系数方差σ作为评价指标，量化分析不同洪水阶段的水文模型参数敏感性，具体公式如下：

30、

31、式中：dc为确定性系数；qest,i和qobs,i分别为i时刻模拟和观测的径流量，m3/s；

32、

33、式中：dcu为使用u×δp为参数的水文模型模拟得到的洪水过程的确定性系数，为确定性系数的均值，σ为确定性系数的方差；

34、σ的值越大，水文模型参数敏感性越高，反之水文模型参数敏感性较低。

35、与现有技术相比，本发明具有以下有益效果：本发明将洪水过程视作时间t的函数q＝f(t)，然后计算q随t的变化率并滤波平滑处理，得到变化趋势线q′sth；在变化趋势线q′sth上寻找由平稳变上涨α和变回平稳的转折点β，以及最高点γ和最低点δ，将洪水过程划分为涨前阶段，上涨阶段，包含洪峰阶段，持续消退阶段和回归稳定阶段；再设置参数值域，在不同阶段内，按照一定的步长使水文模型参数由大到小逐渐变化，模拟得到洪水过程，并与实测洪水比较，计算确定性系数的方差作为评价指标，量化分析不同洪水阶段的水文模型参数敏感性。该分析方法数据来源稳定可靠，变量之间的函数关系明确，有利于不同洪水阶段的水文模型参数敏感性的快速合理分析，可以进一步促进数字水文学以及流域洪水防治研究的深入发展。