[发明专利]河流设计最低通航水位的确定方法和装置

权利要求书

1.一种河流设计最低通航水位的确定方法，其特征在于，包括：

确定河流的目标断面在预设历史时期内的非一致性水位参数序列，其中，所述非一致性水位参数序列包括所述预设历史时期中每一单位时间段内与预设保证率对应的水位参数值；

确定与所述非一致性水位序列对应的预设非平稳皮尔逊Ш型概率密度函数集，其中，所述预设非平稳皮尔逊Ш型概率密度集包括至少两个非平稳皮尔逊Ш型概率密度函数，且所述至少两个非平稳皮尔逊Ш型概率密度函数均包括待确定的时变参数，所述时变参数包括以下至少一项：位置参数、尺度参数和形状参数；

在对数似然函数：取最大值的情况下，分别确定所述预设非平稳皮尔逊Ш型概率密度函数集中每一个非平稳皮尔逊Ш型概率密度函数的时变参数的取值其中，所述f(z_(t)|θ_j(t))为所述预设非平稳皮尔逊Ш型概率密度函数集中的第j个皮尔逊Ш型概率密度函数，所述j为1至S之间的任一整数，所述S为所述预设非平稳皮尔逊Ш型概率密度函数集中包括的非平稳皮尔逊Ш型概率密度函数的数量，所述为所述预设非平稳皮尔逊Ш型概率密度函数集中第j个皮尔逊Ш型概率密度函数的时变参数值，所述N为所述预设历史时期内包含的单位时间段的数量，所述z_(t)为所述非一致性水位参数序列，所述t为1至N中的任一整数；

根据公式：计算出所述预设非平稳皮尔逊Ш型概率密度函数集中每一个非平稳皮尔逊Ш型概率密度函数的评价指标，其中，所述EI_j为所述预设非平稳皮尔逊Ш型概率密度函数集中的第j个非平稳皮尔逊Ш型概率密度函数的评价指标，所述p为预设惩罚函数，所述df_j为所述预设非平稳皮尔逊Ш型概率密度函数集中的第j个非平稳皮尔逊Ш型概率密度函数的整体自由度；

根据与数值最小的评价指标对应的一个非平稳皮尔逊Ш型概率密度函数，确定目标非平稳皮尔逊Ш型分布函数；

根据所述目标非平稳皮尔逊Ш型分布函数，确定所述目标断面在现状期内与预设重现期对应的设计最低通航水位，其中，所述目标非平稳皮尔逊Ш型分布函数包括：时变参数，且所述时变参数根据所述非一致性水位参数序列确定。

2.根据权利要求1所述的方法，其特征在于，所述根据目标非平稳皮尔逊Ш型分布函数，确定所述目标断面在现状期内与预设重现期对应的设计最低通航水位的步骤，包括：

根据所述非一致性水位参数序列，确定水位参数值发生突变的目标单位时间和突变后的水位参数序列，其中，所述突变后的水位参数序列包括所述非一致性水位参数序列中位于所述目标单位时间及位于所述目标单位时间之后的各单位时间内的水位参数值；

根据与所述突变后的水位参数序列对应的目标非平稳皮尔逊Ш型分布函数，确定所述目标断面在现状期内与预设重现期对应的设计最低通航水位。

3.根据权利要求2所述的方法，其特征在于，所述根据所述非一致性水位参数序列，确定水位参数值发生突变的目标单位时间和突变后的水位参数序列的步骤，包括：

采用跳跃诊断法，根据所述非一致性水位参数序列，确定水位参数值发生突变的目标单位时间；

根据所述目标单位时间和所述非一致性水位参数序列，确定突变后的水位参数序列。

4.根据权利要求3所述的方法，其特征在于，所述目标单位时间通过以下过程确定：

根据公式：在所述H取最小值的情况下，确定所述τ的取值；

其中，所述Z(t)为所述非一致性水位参数序列，所述所述所述t为所述单位时间，且所述t为1至N中的任一整数，所述N为所述预设历史时期内包括的单位时间的数量，所述τ为所述目标单位时间。