[发明专利]一种流域水文极值的非平稳性影响指标的确定方法及装置

权利要求书

1.一种流域水文极值的非平稳性影响指标的确定方法，其特征在于，包括：

获取在预设时间段内的流域的水文极值，并对所述水文极值建立多个待检测非平稳极值模型；

确定待检测低频气候指标；

对多个所述待检测非平稳极值模型的时变参数分别进行参数值估计，并对进行参数值估计后的待检测非平稳极值模型进行拟合检测，确定多个所述待检测非平稳极值模型中能够作为对所述水文极值的非平稳性指标进行确定的可用非平稳极值模型；

对所述可用非平稳极值模型进行似然比检验，确定出对所述水文极值的非平稳性影响指标进行确定的最优非平稳极值模型；

将所述待检测低频气候指标作为协变量集添加至所述最优非平稳极值模型的时变参数中，获得进行时变参数添加后的第一非平稳极值模型；其中，所述添加是指将待检测低频气候指标的数值以加法的形式加入至时变参数的表达式中，时变参数包括形状参数、位置参数和尺度参数；

若最优非平稳极值模型的位置参数的表达式为μ(t)＝f(t)，待检测低频气候指标为m个，则在对协变量添加后的位置参数的表达式变换为：

μ(t)＝f(t)+β₁ID₁(t)+β₂ID₂(t)+…+β_mID_m(t)

ID₁(t)、ID₂(t)、ID₃(t)、…、ID_m(t)为m个协变量，β₁、β₂、β₃、…、β_m为系数，m为正整数；

根据所述最优非平稳极值模型和所述第一非平稳极值模型，确定所述待检测低频气候指标是否为所述水文极值的非平稳性影响指标；

对所述可用非平稳极值模型进行似然比检验，确定出对所述水文极值的非平稳性进行确定的最优非平稳极值模型的步骤，包括：

以时变参数的复杂等级对所述可用非平稳极值模型进行排序，确定出时变参数的复杂等级为最低等级的第一可用非平稳极值模型；

获得所述第一可用非平稳极值模型的对数似然函数的第一最大值和剩余的可用非平稳极值模型的对数似然函数的第二最大值；

根据所述第一最大值和所述第二最大值之间的差值，确定出所述最优非平稳极值模型；

根据所述第一最大值和所述第二最大值之间的差值，确定出所述最优非平稳极值模型，包括：

若所述第一最大值和所述第二最大值之间的差值大于预设值，则将所述第二最大值对应的可用非平稳极值模型确定为所述最优非平稳极值模型，其中，所述预设值是自由度为k的卡方分布χ_k²的1-α的分位点，α是指显著性水平，置信度1-α的数值通常为95％，k是进行比较时的两个可用非平稳极值模型的维数差。

2.根据权利要求1所述的流域水文极值的非平稳性影响指标的确定方法，其特征在于，对多个待检测非平稳极值模型的时变参数分别进行参数值估计的步骤，包括：

采用极大似然法对多个待检测非平稳极值模型的时变参数分别进行参数值估计。

3.根据权利要求2所述的流域水文极值的非平稳性影响指标的确定方法，其特征在于，对进行参数值估计后的待检测非平稳极值模型进行拟合检测，确定多个所述待检测非平稳极值模型中能够作为对所述水文极值的非平稳性指标进行确定的可用非平稳极值模型的步骤，包括:

对进行参数估计后的待检测非平稳极值模型进行数据标准变换处理；

获得进行数据标准变换处理后的待检测非平稳极值模型的水文极值的次序统计量；

根据获得的水文极值的次序统计量，获得待检测非平稳极值模型的概率图和分位数图；

在所述分位数图的曲线与理想曲线之间的差值小于预设数值时，确定所述待检测非平稳极值模型为可用非平稳极值模型。