[发明专利]一种针对含有模糊信息历史洪水的L‑moment估参方法

摘要

本发明公开了一种针对含有模糊信息历史洪水的L‑moment估参方法，通过对含有模糊信息的历史洪水，在合理范围内对它可能的量值进行离散，然后求得含有模糊信息的历史洪水在不同离散值下，频率估计结果的离差平方和，通过随机动态规划法求解得到使离差平方和最小时的离散值组合，则在该离散值组合下计算得到的频率分布线型的参数为最优解，该方法既不增加冗余信息，又不提高历史洪水的误差，还能够得到使估计结果最优的参数，拓展了L‑moment方法的适用度，为水库的洪水频率分布函数的参数估计提供重要的参考依据。

主权项

一种针对含有模糊信息历史洪水的L‑moment估参方法，其特征在于包括如下步骤：步骤(1)，收集历史洪水和实测洪水资料，并区分含有模糊信息的历史洪水；所述收集到的历史洪水，设其最大调查期为N，调查的历史洪水个数为m，实测系列的年数为n，实测系列中含特大洪水个数为l；X1,X2,...,Xn‑l,Xn‑l+1,...,Xn‑l+m是样本容量为n‑l+m的样本系列，并按从小到大的顺序排列，用Xj表示系列中的第j个值；对其中认定为可以参与频率分析计算的历史洪水样本，称为较精确的历史洪水；对其中只能估计排位、但不能估计其量值的历史洪水样本，称为含有模糊信息的历史洪水；设样本系列中共有a场历史洪水为含有模糊信息的历史洪水，按照从小到大的顺序排列，分别记为其中n‑l+1＜s1＜s2＜...＜sa＜n‑l+m；步骤(2)，对含有模糊信息的历史洪水的量值在合理区间进行离散；步骤(3)，将不同的离散值代入可考虑历史洪水的L‑moment方法计算公式中，计算频率分布函数的参数，得到频率分布函数；对于第i场含有模糊信息的历史洪水记为将它的第j个离散状态记为其中1≤j≤Max；设取第k1个离散状态，记为设取第ki个离散状态，记为其中1≤ki≤Max；将代入可考虑历史洪水的L‑moment方法计算公式中，求得a场含有模糊信息的历史洪水分别在离散状态k1,k2,…,ka取值时，计算频率分布函数的参数，从而可以得到相应的频率分布函数，记为F(x)；步骤(4)，将具有较精确量值的历史洪水和实测洪水序列代入步骤(3)中得到的频率分布函数，计算相应的理论频率，并通过经验频率公式计算前述具有较精确量值的历史洪水和实测洪水序列的经验频率；步骤(5)，根据步骤(4)中得到的理论频率和经验频率，计算在不同的离散值下频率估计结果的离差平方和；步骤(6)，通过随机动态规划法求解得到使离差平方和最小时的离散值组合，在该离散值组合下计算得到频率分布线型的参数。