本发明涉及一种基于上库弃水期望的梯级弃水概率求解方法。
背景技术：
：水资源缺乏问题是制约国民经济发展的瓶颈，随着水库的修建以及调度水平的提高，洪水资源已经从过去的“猛兽”转变为可合理利用的水资源。如何处理洪水资源利用的风险与防汛安全之间的平衡，在确保安全度汛的前提下充分利用洪水资源已是当前调度工作者面临的严峻挑战。技术实现要素：本发明的目的在于克服现有技术中存在的上述不足，而提供一种可根据上游水库当前状态、区间来水情况、本库当前状态以及未来运行方式求解出本库可能发生弃水的概率，用以指导水库调度工作的基于上库弃水期望的梯级弃水概率求解方法。本发明解决上述问题所采用的技术方案是：一种基于上库弃水期望的梯级弃水概率求解方法，其特征在于：步骤如下：步骤1)：通过收集梯级各站区间洪水峰量资料，进行排频计算，依据绘制上库分期的特征洪量频率曲线，其中p为洪量频率，m为排序号，n为总年数；步骤2)：根据上库分期的特征洪量频率曲线，结合上库当前状态及运行方式，采用面积法求解上库在当前状态及运行方式下的弃水期望e弃以及上库下泄水量w下泄；其中，pi为第i个洪量频率，pi+1为第i+1个洪量频率，w弃i为pi对应的弃水量，w弃i+1为pi+1对应的弃水量；w下泄＝e弃+w发电，w发电为发电用水量；若求解龙头水库的弃水概率，可直接跳过本步骤2)至步骤3)；步骤3)：根据本库当前库水位、运行方式及w下泄，求解本库的可用库容v可用，v可用＝v最大可用-v当前+w发电-w下泄，v最大可用为水库的最大可用库容，v当前为当前水位对应的库容；本步骤3)中，若求解龙头水库弃水概率则w下泄为0。步骤4)：根据本库区间洪量频率曲线，查出v可用所对应的洪水频率，v可用所对应的洪水频率即为弃水概率p弃。步骤1)中，采用梯级各站的区间洪水峰量资料，避免了洪水资料受上库调节作用的影响。步骤2)中，采用面积法求解弃水期望值，利用有限的频率值和对应的弃水量值求解出水库在当前水位和运行方式下的弃水量，解决了弃水期望难以求解的问题。本发明提供了一种直观、可表征防洪风险的指标，本发明在考虑梯级水库间水力联系的基础上，提出基于上游水库弃水期望的梯级弃水概率，即用于指导水库调度工作，对提升洪水资源利用水平，提高水库运行经济效益具有重要意义。本发明所述梯级各站区间洪水峰量资料为流域主汛期内的洪水资料。使得计算资料更加精确。本发明所述特征洪量频率曲线的特征洪量为3d洪量或5d洪量。本发明所述特征洪量频率曲线为皮尔逊ⅲ型曲线。本发明v最大可用为汛限水位与死水位之间的库容，或者，v最大可用为正常蓄水位与死水位之间的库容。本发明步骤2)中，采用面积法求解弃水期望e弃，利用有限的洪量频率和洪量频率对应的弃水量求解上库在当前水位和运行方式下的弃水量，用于解决弃水期望难以求解的问题。本发明与现有技术相比，具有以下优点和效果：考虑了梯级水库间水力联系，水库当前运行状态及本库未来运行方式等信息，解决了水库弃水期望难以量化计算的难题，实现了水库调度中实时面临的弃水概率准确、便捷量化计算，结果更加符合实际，为水库调度工作提供重要指导作用。附图说明图1是本发明实施例基于上库弃水期望的梯级弃水概率求解方法的流程示意图。具体实施方式下面结合附图并通过实施例对本发明作进一步的详细说明，以下实施例是对本发明的解释而本发明并不局限于以下实施例。实施例。参见图1。本实施例为一种基于上库弃水期望的梯级弃水概率求解方法，步骤如下：步骤1)：通过收集梯级各站区间洪水峰量资料，进行排频计算，依据绘制上库分期的特征洪量频率曲线，其中p为洪量频率，m为排序号，n为总年数。作为优选，为了获取可靠性高的梯级各站区间洪水峰量资料，梯级各站区间洪水峰量资料为流域主汛期内的洪水资料，使得计算资料更加精确。采用梯级各站的区间洪水峰量资料，避免了洪水资料受上库调节作用的影响。作为优选，特征洪量频率曲线的特征洪量为3d洪量或5d洪量。作为优选，特征洪量频率曲线为皮尔逊ⅲ型曲线。以西南地区某水库为例，选取可靠性较高的区间洪水峰量资料，绘制分期的特征洪量频率曲线，如表1、表2所示，表1为本库主汛期洪水三日、五日洪量频率曲线，表2为上库主汛期洪水三日、五日洪量频率曲线。表1单位：亿m3项目均值cvcs/cv0.01％0.02％0.05％0.10％0.20％0.50％三天洪量2.770.46412.81210.9109.218.1五天洪量3.770.47417.816.715.11412.811.21％2％5％10％20％50％75％90％95％99％7.266.425.34.453.582.411.861.601.491.41108.867.296.094.893.262.502.162.021.92表2单位：亿m3项目均值cvcs/cv0.01％0.02％0.05％0.10％0.20％0.50％三天洪量1.610.623.57.597.076.395.875.364.67五天洪量2.070.643.51110.29.228.457.686.661％2％5％10％20％50％75％90％95％99％4.143.622.936.271.881.170.900.810.740.695.895.114.18.672.561.551.151.080.970.89步骤2)：根据上库分期的特征洪量频率曲线(以三日洪量为例)，结合上库当前状态及运行方式，采用面积法求解上库在当前状态及运行方式下的弃水期望e弃以及上库下泄水量w下泄。其中，pi为第i个洪量频率，pi+1为第i+1个洪量频率，w弃i为pi对应的弃水量，w弃i+1为pi+1对应的弃水量；w下泄＝e弃+w发电，w发电为发电用水量。采用面积法求解弃水期望值，利用有限的频率值和对应的弃水量值求解出水库在当前水位和运行方式下的弃水量，解决了弃水期望难以求解的问题。本实施例弃水期望e弃以及上库下泄水量w下泄如表3所示。表3为上库主汛期三日洪量弃水期望求解表。表3单位：m、亿m3注：上库三日满发用水量为2.58亿m3。若求解龙头水库的弃水概率，可直接跳过本步骤2)至步骤3)。步骤3)：根据本库当前库水位、运行方式及w下泄，求解本库的可用库容v可用，v可用＝v最大可用-v当前+w发电-w下泄，v最大可用为水库的最大可用库容，v当前为当前水位对应的库容。作为优选，v最大可用为汛限水位与死水位之间的库容，或者，v最大可用为正常蓄水位与死水位之间的库容。本库可用库容求解表参见表4。表4单位：m、亿m3本步骤3)中，若求解龙头水库弃水概率则w下泄为0。步骤4)：根据本库区间洪量频率曲线，查出v可用所对应的洪水频率，v可用所对应的洪水频率即为弃水概率p弃。本库弃水概率求解表如表5所示。表5单位：m、亿m3本实施例基于上库弃水期望的梯级弃水概率求解方法不同于传统方法，通过引入上库弃水期望，根据上游水库当前状态、区间来水情况、本库当前状态以及未来运行方式求解出本库可能发生弃水的概率。本实施例的基于上库弃水期望的梯级弃水概率求解方法更加符合实际情况，实用化程度更高的水库弃水概率求解方法。本实施例通过直观、可表征防洪风险的指标，在考虑梯级水库间水力联系的基础上，提出基于上游水库弃水期望的梯级弃水概率，即用于指导水库调度工作，对提升洪水资源利用水平，提高水库运行经济效益具有重要意义。此外，需要说明的是，本说明书中所描述的具体实施例，其零、部件的形状、所取名称等可以不同，本说明书中所描述的以上内容仅仅是对本发明结构所作的举例说明。凡依据本发明专利构思所述的构造、特征及原理所做的等效变化或者简单变化，均包括于本发明专利的保护范围内。本发明所属
技术领域：
的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代，只要不偏离本发明的结构或者超越本权利要求书所定义的范围，均应属于本发明的保护范围。当前第1页12