[发明专利]一种绘制P‑III型分布频率曲线的优化方法

说明书

技术领域

本发明涉及一种频率曲线的优化方法，特别涉及一种绘制P-III型分布频率曲线的优化方法。

背景技术

洪水频率分析是估计水利水电工程设计洪水的主要内容，也是水利工程建设时计算设计标准的重要依据。通过构建洪水频率分析模型来进行设计洪水计算，该模型由洪水频率分布函数和参数估计方法构成。20世纪60年代，我国根据大量长期洪水系列分析结果和多年来设计工作的实际经验，规定P-III型分布作为我国水文分析计算的分布函数，采用矩法估计参数的初值，采用数学期望公式计算经验频率，将各项洪水的数值和相应的经验频率点绘在格纸上，以统计参数的初估值绘制频率曲线，如果曲线与经验频率点据拟合不够好，可调整统计参数，直至曲线与点据拟合较好、较满意为止，由此确定采用的统计参数，这种参数估计方法被称为目估适线法(CF)，中国的洪水频率分析模型被称为P-III/CF模型。

按照P-III/CF模型的适线方法，在频率曲线定线时，一般要尽可能地照顾历史洪水和特大洪水点据，即绘制的频率曲线上尾部要尽量靠近它。由于不同工程设计人员的经验、视差有差异，采用该方法得到的洪水频率曲线存在较大的任意性，在工程实际中经常出现不同的设计部门会计算得到不同设计结果的情况，难以制定统一的工程设计标准。而且，采用人工多次调整参数的方法，也耗费了大量的时间和精力。

刘力等(刘力、周建中、杨俊杰等，粒子群优化适线法在水文频率分析中的应用.水文，2009,29(2)，21—23.)提出了基于粒子群优化算法的优化适线法，并将该方法应用到年最大流量频率分析中，该方法虽然在一定程度上避免了人为因素的不确定性，但是它以样本数据的实际值和理论曲线估计出来的在同频率下的估计值的离差绝对值的和最小为目标函数，没有照顾到历史洪水和特大洪水，不符合中国工程设计的实际要求。李宏伟等(李宏伟，宋松柏.蚁群算法在水文频率曲线参数计算中的应用.人民黄河，2009,31(4)，38—40.)提出了一种可以应用于水文频率参数计算的基于网格划分策略的蚁群算法，该方法不仅参数众多，操作繁琐，同样没有照顾到历史洪水和特大洪水。

发明内容

针对现有技术存在的不足，本发明提供了一种绘制P-III型分布频率曲线的优化方法。

本发明通过建立P-III型分布频率曲线优化模型，对理论频率和经验频率的相对误差赋予不同的权重，以相对误差权重的平方和最小为目标函数，并根据流域概况和水库的实际情况添加参数的约束条件，以矩法估计的参数为初值，采用改进的非常快速模拟退火算法(MVFSA)求解模型，得到参数的最优解，从而绘制出P-III型分布频率曲线。

为解决上述技术问题，本发明采用如下的技术方案：

一种绘制P-III型分布频率曲线的优化方法，包括如下步骤：

步骤1，收集历史洪水和实测洪水资料；

对于收集到的历史洪水，设其最大调查期为N，调查的历史洪水个数为m，实测系列的年数为n，实测系列中含特大洪水个数为l；X₁,X₂,...,X_n-l,X_n-l+1,...,X_n-l+m是样本容量为n-l+m的样本系列，并按从小到大的顺序排列，用X_j表示系列中的第j个值；

步骤2，计算洪水序列的经验频率，并用矩法初步估计分布函数的参数；

对于步骤1中得到的样本系列X_j，将它的经验频率记作P₀(X_j)；

步骤3，建立P-III型分布频率曲线优化模型；

将步骤1中的洪水数据代入步骤2得到的频率分布函数，可以求得不同的洪水点据对应的理论频率，记作P₁(X_j),j＝1,2,...,n-l+m，从而可以得到理论频率与经验频率的相对误差RE(X_j)，对不同的洪水点据，赋予相应的权重α_j；采用相对误差权重的平方和最小为目标函数，并根据流域概况和水库的实际情况添加参数的约束条件，建立P-III型分布频率曲线优化模型；

步骤4，采用改进的非常快速模拟退火算法(MVFSA)求解步骤3中的模型，从而得到P-III型分布函数的各项参数，从而绘制频率曲线。

所述改进的非常快速模拟退火算法(MVFSA)为现有技术，具体参见同济大学学报(自然科学版)2006年第34卷第8期，1121—1125页，作者为陈华根，李丽华，许惠平，陈冰。