考虑水库调度影响的洪水频率分析方法

技术领域

本发明属于洪水频率分析技术领域，具体涉及考虑水库调度影响的洪水频率分析方法。

背景技术

洪水频率分析法可用于推求设计洪水，广泛应用于防洪工程实践中。传统的研究方法主要包括：美国1968年撰写了设计洪水计算指南；爱尔兰国立大学Cunnane教授受世界气象组织(WMO)委托，于1989年撰写出版了《洪水频率分析的统计分布》；Hosking和Wallis 1997年提出了以线性矩为基础的区域洪水频率分析方法；中国自20世纪50年代开始开展洪水频率分析工作，在频率曲线线型、经验频率公式、参数估计方法、设计洪水过程线放大等方面进行了较为深入的研究。

Eagleson首先提出了基于物理机制的洪水频率曲线公式推导法，这一方法主要包含降水和径流模型的设定以及径流峰值与降雨转换关系。Sivapalan将气候变化和土地利用等因素纳入洪水频率分析中，推导了考虑季节影响下的洪水频率公式。Botter假定降水服从泊松分布，通过体现气候变化和土地利用影响的模型参数，推导了基于物理机制的洪水频率曲线解析模型。

传统的研究侧重于考虑气候变化及土地利用对洪水频率分析的影响，却忽视了水库调度对自然状态下的水资源时空特征进行重新分配带来的影响。少量研究表明水库建设前后会对洪水洪峰洪量有一定程度的削减作用，但水库对洪水频率分析的具体影响尚未有解析公式。

发明内容

本发明是为了解决上述问题而进行的，目的在于提供一种考虑水库调度影响的洪水频率分析方法。本发明为了实现上述目的，采用了以下方案：

本发明提供一种考虑水库调度影响的洪水频率分析方法，其特征在于，包括以下步骤：步骤1：根据水库库容与水库面积、水位关系，以及出库流量与水位关系得到水库库容与出库流量的相关关系；步骤2：入库出库洪水过程采用简化三角形方式处理，入库出库三角形面积差为入库洪水形成的水库最大调蓄水量得到水量平衡关系式；步骤3：联立求解步骤1和2得到的关系式，推求出库流量与入库流量的关系式，并将该关系式代入无水库调度影响下的洪水频率公式，推导出水库调度影响下的洪水频率公式。

本发明提供的考虑水库调度影响的洪水频率分析方法，还可以具有以下特征：在步骤1中，得到的水库库容与出库流量的关系式为：式中，V为水库最大调蓄水量，a为水库特征参数，Q_out为水库调节下的洪峰出库流量。

本发明提供的考虑水库调度影响的洪水频率分析方法，还可以具有以下特征：在步骤2中，简化三角形方式为线性处理洪水过程线。

本发明提供的考虑水库调度影响的洪水频率分析方法，还可以具有以下特征：在步骤2中，水量平衡关系式为：式中，T为涨水退水历时，Q_in为无水库调节下的洪峰入库流量。

本发明提供的考虑水库调度影响的洪水频率分析方法，还可以具有以下特征：在步骤3中，推求出的出库流量与入库流量的关系式为：

本发明提供的考虑水库调度影响的洪水频率分析方法，还可以具有以下特征：在步骤3中，推导出的水库调度影响下的洪水频率公式为：

式中，β,α,δ依次为尺度参数、形状参数、位置参数。

发明的作用与效果

与现有技术相比，本发明的有益效果在于：

(1)本发明充分考虑了水库调度影响下的洪峰流量变化，简化了水库入库出库流量的计算方法，推导出了水库调度影响下的洪水频率公式。

(2)本发明能够将人类活动的影响(水库建设)纳入洪水频率分析考虑中，有助于了解水库调度下水资源时空特征变化情况。

(3)本发明可广泛应用于水库调度影响下洪水频率分析中，为面对人类活动影响下的水文特征变化提供科学的技术支撑。

附图说明

图1为本发明实施例中所涉及的考虑水库调度影响的洪水频率分析方法的流程图；

图2为本发明实施例中所涉及的简化三角形处理方法的示意图；以及

图3为本发明实施例中所推求出的水布垭水库的洪水频率曲线图。

具体实施方式

以下结合附图对本发明涉及的考虑水库调度影响的洪水频率分析方法、以及该方法应用于非一致性条件下防洪调度方案优化问题的具体实施方案进行详细地说明。

<实施例>

如图1所示，本实施例所提供的考虑水库调度影响的洪水频率分析方法包括以下步骤：

步骤1.根据水库库容与水库面积、水位关系，以及出库流量与水位关系得到水库库容与出库流量的相关关系；

由于库容及库水位存在三次方关系，自由溢流水库存在出库流量与库水位有二次方关系，由此得出水库库容与出流流量的关系式：

式中，V为水库最大调蓄水量，a为水库特征参数，Q_out为水库调节下的洪峰出库流量；

步骤2.采用简化三角形方式处理入库出库洪水过程，并将水库出流过程按照自由溢流方式处理，如图2所示，当水位处于起调水位时，出库流量为0；出库流量按照水库库容～出流流量关系计算，出库流量过程与入库过程相交于点C，过后出库流量随库水位降低而逐步减小，因此可将出库过程的涨水段简化为直线，图2中三角形ΔDAB为入库过程，ΔCAB为出库过程，涨水退水历时T的计算：已知水库集水面积、河长、坡降等参数时，采用设计洪水计算方法推算，由水量平衡关系式得到水库最大调蓄水量V为入库出库三角形面积差：

式中，T为涨水退水历时，Q_in为无水库调节下的洪峰入库流量；

步骤3.联立求解上述步骤1和2得到的关系式，推求出库流量与入库流量的关系式：

假设无水库调度影响下的洪峰序列服从皮尔逊Ⅲ型频率曲线公式，即：

式中，β,α,δ依次为尺度参数、形状参数、位置参数；

求Q_out关于Q_in的偏导并将其代入原洪水频率分析公式中可推导出水库调度影响下的洪水频率公式：

上述公式5便是基于水库调度影响下的洪水频率分析公式，由该公式可以看出，水库调度改变了原有的洪水频率分析公式。

另外，为了进行误差估计，可以在上述步骤1的基础上，进一步进行以下步骤2’和3’：

步骤2’.入库洪水过程采用线性处理，出库洪水过程采用实际曲线积分计算，得到水量平衡关系式：

σdT＝δV^0.5dT+dV，(公式6)

式中，σ为入库洪水系数，δ为出库洪水系数；

求解微分方程可得V和T的关系式：

式中，函数P(x)为关于x的牛顿迭代解析式；

步骤3’.根据步骤2’中水库最大调蓄水量V和涨水退水历时T的关系式，及步骤1中水库最大调蓄水量V和出库流量Q_out的关系式，得到出库流量Q_out和涨水退水历时T的关系式，通过牛顿迭代公式可求得出库流量Q_out数值解：

将公式8代入原洪水频率分析公式4中，推导出水库调度影响下的洪水频率公式：

然后，采用牛顿迭代编程计算出库流量Q_out的误差，将其与实际流量和公式5一起比较进行误差估计。

下面以清江流域中的水布垭水库为例，对本实施例所提供的洪水频率分析方法进行验证：

清江流域横贯湖北省西南，干流全长423km，总落差1430m。流域呈南北窄，东西长的狭长形，面积约为17000km²。水布垭水库坝址位于清江流域中游，距恩施约117km，距渔峡口水文站17.6km，控制面积为10860km²，比渔峡口集水面积(11906km²)小8.8％。渔峡口水文站是水布垭的设计代表站。

将水布垭水库作为实例研究，根据其62年(1951-2012)实测径流资料，绘制如图3所示的频率曲线，比较原频率曲线与推导得到的频率曲线公式，可见推导得到的频率曲线线型和原皮尔逊Ⅲ型曲线线型不一致。由于水库调节的消峰影响，出库流量明显比入库流量小。由图3可知，所推导的频率曲线公式拟合出库流量的效果比原皮尔逊Ⅲ型曲线好。推导的频率曲线公式与实际出库流量经验频率的均方根误差为0.09，小于原皮尔逊Ⅲ型曲线的均方根误差0.31。由此可验证结论：水库调度改变了原有的洪水频率分析公式。

以上实施例仅仅是对本发明技术方案所做的举例说明。本发明所涉及的考虑水库调度影响的洪水频率分析方法并不仅仅限定于在以上实施例中所描述的内容，而是以权利要求所限定的范围为准。本发明所属领域技术人员在该实施例的基础上所做的任何修改或补充或等效替换，都在本发明的权利要求所要求保护的范围内。