基于年防洪标准的水库分期优化调度研究

2018-03-16蔡振华刘英豪于东平

水资源与水工程学报 2018年1期

张涛，蔡振华，刘英豪，于东平

(1.山东省水利勘测设计院，山东济南 250013； 2.山东省水文水资源勘测局，山东济南 250014)

1 研究背景

汛限水位分期调度考虑汛期洪水的季节变化特性，在不同的汛期分期执行不同的汛限水位方案，更具有科学性、合理性。相对于汛限水位动态控制须借助于水文气象预报等一系列技术设施辅助手段，分期控制操作相对简单，易于执行[1]，目前我国已有多座水库采用了分期调度，获得较好的兴利效益。分期控制常规计算方法与全汛期固定控制的计算方法基本相同[2]，包括汛期划分、分期洪水设计及洪水调算等。采用分期调度是否会降低年防洪标准，目前存在一定争议[3-5]，主要是分期设计洪水与年设计洪水标准关系[5-7]，以及分期防洪风险组合与年防洪风险标准关系的一致性问题[8]。王善序等[6]认为现行分期洪水计算频率与重现期的倒数存在一定差异，现行汛限水位分期调度确定方法存在一定的防洪安全隐患；陈守煜等[9]采用模糊统计法划分汛期及确定分期汛限水位；曹升乐[10]根据丰枯季节暴雨特性确定分期汛限水位过程线，但分期防洪风险组合与年防洪标准的关系具有不确定性；文献[11-12]通过对分期防洪概率事件进行概率组合，在满足年防洪标准的条件下可获得一定的兴利效益。本文从分期洪水组合概率与年设计洪水概率关系一致性，以及分期防洪风险率组合与年防洪风险标准关系一致性两个方面进行探讨，以棉花滩水库为例，在不降低年防洪标准的前提下，考虑调洪、发电规则及约束条件进行汛限水位分期优化设计，寻求最大兴利效益，并确保水库防洪安全。

2 分期设计洪水计算

2.1 洪水计算方法

分期设计洪水是水库分期调洪计算的基本数据，汛期划分则直接关系到分期洪水取样及分期调度的合理性。根据洪水的季节变化规律合理划分汛期[9]，并考虑到运行管理方便，常对汛期实行硬性划分，如以整月或整旬为单元划分为汛期。采用传统水文频率分析法推求分期设计洪水，洪水样本选样为相应分期最大样本选样。考虑流域洪水的时空分布不确定性因素，适当进行跨期取样[7]，跨期通常取5～10 d。采用P-Ⅲ型等频率曲线进行频率适线，按同倍比或同频率法对典型洪水进行放大从而得到设计洪水过程线。

2.2 分期设计洪水概率组合

假定汛期分为前汛期和后汛期，年最大洪水均发生在汛期内。同一年份中，前汛期、后汛期最大洪水样本分别为Y1、Y2，全汛期最大洪水样本为Y0，即Y0=max(Y1,Y2)，其中max()为取大函数。假定前、后汛期洪水事件相互独立，其概率组合关系[2]为：

P0(Y0≥y)=P1(Y1≥y)∪P2(Y2≥y)=P1(Y1≥y)+P2(Y2≥y)-P1(Y1≥y)P2(Y2≥y)

(1)

式中：P0、P1、P2分别为全汛期、前汛期、后汛期洪水样本变量不小于y的上区间概率。

设计洪水频率计算是以洪水样本序列为依据进行频率适线，分期设计洪水组合概率应与年设计洪水概率一致[7]，即满足上述概率组合公式(1)。由于受分期跨期取样、样本容量、频率参数选取等影响，分期洪水组合与年洪水并非完全符合概率组合公式(1)，通常存在一定的偏差。

3 分期优化调度设计模型

3.1 防洪风险率计算

根据调洪规则对设计洪水进行调洪计算，当调洪高水位等于或超过防洪高水位则认为是风险事件。给定汛限水位，当对应分期设计洪水频率为p的调洪高水位达到防洪高水位时，则认为频率p即为该汛限水位下的防洪风险率。

可调控洪水下，洪水量级与防洪库容为单调函数关系，当洪水频率越小、洪水量级越大时，则需要的防洪库容则越大。为便于快速计算各汛限水位的防洪风险率，根据分期洪水的防洪库容推求防洪风险率。假定可调控洪水的调洪规则为：洪水初期洪水来多少泄多少，大于安全泄流时则按安全泄流控泄，确保调洪高水位不超过防洪高水位，设计洪水频率p下的防洪库容Wp计算公式为：

(2)

式中：Qt,p、qt,p分别为洪水频率p下的t时刻设计入库流量、相应出库流量，m3/s；Δt为时段时长，s；T为洪水离散时段数；min()为取小函数；qt,an为安全泄流量，m3/s。

对频率微分，建立洪水频率pi与相应洪水防洪库容Wpi的两变量序列。给定汛限水位h，计算该水位与防洪高水位之间的防洪库容Wh，由已知Wpi～pi两变量序列采用线性插值求得对应库容Wh的洪水频率ph，即为该汛限水位下的防洪风险率。为简化计算，本文未考虑洪水预报、调度操作等因素[13]对防洪风险的不确定性影响。

3.2 年防洪标准下的分期防洪风险组合

假定前汛期、后汛期的风险事件相互独立，根据全概率公式，全汛期的防洪风险概率组合公式为：

Pf[hmax(h)≥H]=Pf,1[hmax(h1)≥H]+

Pf,2[hmax(h2)≥H]-Pf,1Pf,2

(3)

式中：Pf[hmax(h)≥H]为全汛期汛限水位h下的调洪高水位不低于H的风险率，简记为Pf；hmax(h)为对应汛限水位h的调洪高水位；Pf,1[hmax(h)≥H]为前汛期汛限水位h1下的调洪高水位不低于H的风险率，简记为Pf,1;Pf,2[(hmax(h2)≥H]为后汛期汛限水位h2下的调洪高水位不低于H的风险率，简记为Pf,2。

当全汛期风险率一定时，则存在满足条件的多种分期汛限水位概率组合[11]，通过对各分期进行风险及效益转移，寻求整体最大效益。一般情况下，主汛期大洪水发生频次高于次汛期，其防洪任务较重，因此，可通过降低前汛期的兴利效益，以提高防洪效益；降低后汛期的防洪效益，以提高兴利效益。

3.3 效益最大化目标

分期汛限水位防洪风险率组合按年防洪标准控制，拟定汛限水位组合方案，根据长期发电调度规则、调洪规则进行长系列模拟，计算主要效益指标，以综合效益最大化为目标确定最优方案。

(1)发电效益。考虑洪水的丰枯季节变化特性，在枯水期内多为低出力运行，与丰水期相比，枯水期的发电效益增加效果更加显著。另外，水电站保证出力的日保证率越大，则对电网的平衡作用越稳定。考虑电网平衡作用，不同季节时段采用不同的效益权重，发电效益为：

(4)

式中：Ec(x)为方案x下的年均发电量，kW·h；Ei,n(x)为方案x下第n年第i分期的发电量，kW·h；N为年数；I为年内分期数；Ci,n为第n年第i分期的的效益权重。

考虑权重时发电效益需进行归一化处理，归一化公式为：

(5)

(6)

式中：K为较大洪水场次；VH为防洪高水位对应的库容，m3；Vk,max(x)为方案x下第k场洪水调洪高水位对应的库容，m3。

(3)供水效益。主要包括农业、生活、工业、生态等供水效益。

(4)综合效益最大化目标。结合约束条件，对各效益指标进行经济效益定量量化，如发电效益通过发电量及电网价格衡量，供水效益通过供水量及水费价格衡量，效益最大化目标函数为：

(7)

(j=1,2,…,J)

式中：Bm(x)为方案x下的第m个效益指标的经济效益；M为效益指标个数；Rj(x)为x方案下第j个约束，如下泄流量等；Rj,min、Rj,max为第j个约束的最小、最大约束范围；J为约束个数。

对于以发电为主的水库，在满足各约束条件前提下提高发电效益是主要目标，当其它指标差异不大的情况下，可以发电效益等主目标进行方案优选。当各方案指标差异不明显、难于定量量化时，可采用多目标模糊评价法[15]进行方案评价。

3.4 约束条件

调洪及发电调度需满足库水位限制、下游防洪安全、机组出力等约束条件。

(1)库水位。包括水库死水位、限制水位、防洪高水位等约束，库水位Zt满足：

Zt,min≤Zt≤Zt,max

(8)

式中：Zt,min为t时刻水库死水位或其它运用要求的最低水位，m；Zt,max为t时刻上限水位，m，其中汛期洪水期内为防洪高水位(遭遇超标准洪水时则按调洪规则确定)，汛期非洪水期内为汛限水位，非汛期为正常蓄水位。

(2)电站出力。包括水电站发电降低出力、预想出力等约束，出力Nt满足：

Nt,low≤Nt≤Nt,exp

(9)

式中：Nt,low、Nt,exp分别为t时刻水电站的降低出力、预想出力，kW。

(3)发电流量。包括水轮机发电允许最小、最大流量等约束，发电流量根据出力、上游水位及尾水位等确定，发电流量Qt,E满足：

(10)

式中：QE,min、QE,max分别为水轮机允许的的最小、最大过流能力，m3/s；A为出力系数；Zt,wei为t时刻发电尾水位，m。

(4)下泄流量。水库控泄阶段的下游用水、最大泄流等约束，控泄流量qt满足：

qt,min≤qt≤qt,an

(11)

式中：qt,min为t时刻水库下游用水下限流量，m3/s；qt,an为t时刻下游安全泄流量，m3/s。

4 分期调度实例计算

4.1 水库基本情况

(1)主要参数。棉花滩水库是以发电为主，兼有防洪、航运等功能的大型水库，总库容20.35×108m3，机组保证出力8.8×104kW；防洪标准为50年一遇，防洪高水位为173m；前汛期为5月1日-6月30日，后汛期为7月1日-8月25日。

(2)调洪及发电基本规则。洪水调度以库水位和入库流量作为判别标准，前汛期、后汛期安全泄流分别为5 210、4 650m3/s，当库水位低于173m时，泄流按来水控制且不大于安全泄流；当库水位高于173m时，逐步加大泄流直至全部敞泄。汛期当遇到洪水则蓄到汛限水位，尽量减少弃水，加大出力；在汛限水位附近及以上时满荷预想出力，尽快回落到汛限水位。

(3)约束条件。水库运行最低水位146.0m，长系列中无超标准洪水，汛期洪水期内的上限水位取173m，汛期非洪水期内的上限水位为相应汛限水位，非汛期的上限水位为173m；电站降低出力6.8×104kW，预想出力60×104kW；水轮机允许的的最小、最大流量分别为70、280m3/s；水库下游用水下限为40m3/s，控泄阶段泄流不超过相应安全泄流。

(4)分期设计洪水计算。采用P-Ⅲ型曲线对洪水样本进行频率适线；以前汛期(1973年6月)、后汛期(1959年8月)大洪水为典型洪水，采用同频率法对洪峰、洪量进行过程放大得到设计洪水过程。设计洪水频率参数、典型频率设计成果以及典型洪水统计参数见表1。表1中Qm为洪峰；W24、W72、W120分别为最大24、72、120h洪量。

表1 设计洪水及典型洪水统计参数 m3/s， 108 m3

通过验算，前、后汛期洪水组合概率与全汛期概率基本一致。如洪峰为7 783m3/s时，对应全汛期洪峰频率为2%；对应前、后汛期洪峰频率分别为1.49%、0.58%，根据公式(1)，概率组合为2.06%，与全汛期频率2%基本一致。

4.2 防洪风险率计算

拟定汛限水位试算范围前汛期为168.40～169.85m，后汛期为169.00～170.45m，水位步长0.01m，共146组。假定洪水均为可调控，根据调洪规则，推求各频率分期设计洪水的防洪库容，频率区间[0.01%～5%]，频率步长0.01%，得到洪水频率与相应防洪库容两变量序列。分期防洪库容与风险率、库水位与相应防洪库容关系曲线见图1。根据拟定分期汛限水位与防洪高水位之间的库容，采用线性插值方法求得分期汛限水位对应的防洪风险率，见图2。

根据公式(3)，按年防洪风险率2%进行分期汛限水位防洪风险率组合，部分方案组合见表2。

表2 分期汛限水位防洪风险组合

4.3 长系列模拟计算

该水库于2001年建成，建成后的天然水文资料一致性较差。因此，主要对1959年1月1日-2001年12月31日的日流量，根据长期发电、调洪规则进

行长系列模拟，日流量系列见图3。前汛期汛限水位为168.50～169.5m，水位步长0.01m，后汛期为组合防洪风险率为2%的相应汛限水位；初始调算水位为死水位146.0m。

图3 1959-2001年入库日流量系列

计算指标包括全年、非汛期、前汛期、后汛期发电量以及发电、日保证率、洪峰≥2 000m3/s的大洪水平均余留库容(简记平均余留库容)、弃水量、超安全泄流次数。结合地区电力平衡作用，各分期的发电效益权重系数非汛期取1.25、前汛期取0.9、后汛期取1.05，通过归一化处理后求得全年(权) 发电量。部分汛限水位组合方案的计算结果见表3；平均全年发电量、全年(权)发电量以及组合防洪风险率为2%下的前、后汛期汛限水位组合见图4。

表3 长系列模拟计算指标

图4 分期汛限水位组合及发电效益图

4.4 结果分析

(1)效益分析。方案10中，汛期采用单一汛限水位169.12m，年均发电量为16.075×108kW·h，发电效益相对较小，日保证率93.10%，平均余留库容2.409×108m3。方案4中，汛限水位前汛期取168.78m、后汛期取169.58m时，年均发电量为16.087×108kW·h，相对于方案10，平均每年多发电量120×104kW·h，非汛期多发电量640×104kW·h，日保证率93.82%，平均余留库容2.475×108m3。该方案的年发电量、日保证率、平均余留库容等指标均优于方案10相应指标，其不考虑权重和考虑权重的年均发电量均为最大。

(2)上游洪水遭遇分析。根据1959-2001年洪水统计，日最大流量大于2 000m3/s的洪水共计61场，其中前汛期45场、后汛期16场；日最大流量大于3 000m3/s的洪水共计13场，其中前汛期10场、后汛期3场。前汛期发生大洪水的场次基本为后汛期的3倍，一般情况下，前汛期汛限水位越低，全汛期大洪水的平均余留库容则越大，有利于防洪安全，因此宜适当降低前汛期汛限水位，以获得较大的防洪效益。

(3)下游泄洪分析。长系列模拟的日水位均未超过防洪高水位，最大泄流均未超过安全泄流。最大场次洪水为1973年前汛期6月中旬的大洪水，最大24h洪量达5.347×108m3，方案4中的调洪高水位为170.23m，最大泄流达到安全泄流5 210m3/s。

综上分析，以发电效益作为主目标评价各方案时，方案4为最优。即汛限水位前汛期取168.78m、后汛期取169.58m，分期防洪风险率分别为1.217%、0.783%，组合风险率为1.990%，能满足年防洪标准及各约束条件；相对于单一汛限水位控制方案，多年平均全年、非汛期发电量提高比例分别为0.07%、0.74%，具一定的兴利效益。