CN105427052A - 一种基于参考线的并联水库确定性优化调度方法 - Google Patents

Abstract

本发明公开了一种基于参考线的并联水库确定性优化调度方法，属于水利电力技术领域；本发明包括以下步骤：(1)确定调度期；(2)建立水库群确定性长期优化调度发电量最大模型；(3)确定水库群初始参考线和初始调度方案；(4)判断初始调度方案是否可行。本发明将水库最小决策限制线和防弃水线结合起来使用，可准确确定水库增发单位出力引起的不蓄电能损失值(K值)计算的起始和终止时刻，不要求各水库使用统一的调度起始和终止时刻，避免了各水库蓄、供水期划分的不同对调度结果的影响；将专家经验和水库调度理论相结合，避免了水库群调度“维数障碍”问题，适用于大规模并联水库群优化调度。

Description

一种基于参考线的并联水库确定性优化调度方法

技术领域

本发明属于水利电力技术领域，更具体地，涉及一种基于参考线的并联水库确定性优化调度方法。

背景技术

水库群调度是一个典型的多目标、多阶段、多约束、非线性优化问题，广泛应用的水库供蓄水判别系数法的最主要的缺陷是需人为地将调度期划分为供水期、蓄水期、不蓄不供期等时期，由于蓄、供水期的划分与水库调度周期内的来水过程、水库发电运行方式等因素密切相关，其起止时刻具有不确定性，难以准确确定；传统的优化方法为动态规划及其改进算法，其中改进的算法主要包括逐步优化法、动态规划逐次逼近法、状态逐密动态规划法和增量动态规划法等，与全状态枚举法相比，这些方法虽然大大减小了计算工作量，但却未从根本上克服水库群调度“维数障碍”。

发明内容

针对现有技术中存在的问题，本申请提供的是一种基于参考线的并联水库确定性优化调度方法，其中通过对并联水库的参考线进行研究和设计，在保证调度方案最优性的前提下，极大地缩短了计算时间，因而尤其适用于大规模并联水库群的调度应用场合。

为实现上述目的，本发明提供了1、一种基于参考线的并联水库确定性优化调度方法，其特征在于，所述方法包括以下步骤：

步骤(1)，设置调度期，将所述调度期均匀划分为T个等间距的时段,其中T为正整数；设置所述并联水库调度期内的初始条件和约束条件；所述初始条件包括水库调度期期初蓄水状态、期末蓄水状态和调度期内各时段入库流量；所述约束条件包括水库蓄水量约束、出力约束、出库流量约束以及水量平衡方程；

步骤(2)，根据所述并联水库中各水库的库容特性曲线、尾水关系曲线、水库优化调度的目标函数和所述约束条件，建立水库群长期优化调度发电量最大模型：

$E=max\underset{m=1}{\overset{M}{\Σ}}\underset{t=1}{\overset{T}{\Σ}}{N}_{m,t}\·\mathrm{\Δ}t$

N_m,t＝K_mQ_m,tH_m,t

其中，E表示水库群全时期最优发电量；N_m,t、Q_m,t和H_m,t分别表示t时段水库m的出力、发电流量和发电水头；K_m表示水库m的综合出力系数；

步骤(3)，确定所述各水库的初始参考线，包括确定所述各水库的防弃水线和最小决策限制线；

步骤(4)，以所述各水库的初始参考线为其初始调度线，根据入库流量信息和水量平衡方程从前往后逐时段递推至最后一个时段，计算所述各水库各时段的出库流量和出力值，直至所有水库计算完毕,统计所述并联水库各时段的总出力值；

步骤(5)，从前往后逐时段检查所述并联水库各时段的总出力值是否大于保证出力，是，则转下一时段；否，则计算不足出力并按所述并联水库不蓄电能损失最小原则调整所述各水库的时段出力值，调整完毕后根据调整后的所述各水库的蓄水状态按所述步骤(3)所述方法调整后续时段所述各水库的最小决策限制线，调整完毕后转下一时段；上述过程一直持续至调度期末；根据等微增率原理，调整后所述各水库增发单位出力引起的不蓄电能损失值，即K值，应大致相当；以待调整时段t₀所述各水库调整出力之后的K值为搜索变量，以t₀时段各水库调整之后的出力总和等于保证出力为搜索条件采用一维搜索方法调整所述各水库的初始参考线；

步骤(6)，输出所述各水库调整后的参考线和各时段出力值。

优选地，所述步骤(3)包括以下子步骤：

(3-1)以减少弃水能为原则确定所述各水库的防弃水线，其包括以下子步骤：

(3-1-1)已知所述并联水库中的某一水库一个时段的时段末蓄水状态和入库流量信息，令所述时段出库流量等于水轮机最大过流，由水量平衡方程求得所述某一水库所述时段初的蓄水状态；

(3-1-2)根据所述某一水库所述时段初的蓄水状态、时段末的蓄水状态、入库流量和最大容许出力，计算求得所述某一水库的弃水流量，若所述弃水流量大于零，则减小所述某一水库的出库流量使得所述弃水流量为零，并根据所述水量平衡方程重新计算所述时段初的蓄水状态，若所述弃水流量等于零，则所述时段初的蓄水状态取所述步骤(3-1-1)所求得的值，将计算得到的所述时段初的蓄水状态作为上一时段末的蓄水状态；

(3-1-3)按所述步骤(3-1-1)和(3-1-2)从调度期末递推计算至调度期初，得到所述某一水库的一条蓄水过程线；

(3-1-4)对(3-1-3)所述蓄水过程线从调度期初至调度期末逐时段检查是否有弃水发生，若各时段均无弃水，则所述蓄水过程线即为所述某一水库的防弃水线；若某时段有弃水发生，则进行相应调整，使弃水集中在最后一个或若干个时段，调整后的蓄水过程线即为所述某一水库的防弃水线；

(3-1-5)按所述步骤(3-1-1)～(3-1-4)的方法求得所述各水库的防弃水线；

(3-2)根据保持高水位运行以减少不蓄电能损失的原则，确定所述各水库的最小决策限制线，其包括以下子步骤：

(3-2-1)根据所述并联水库中的一个水库某一时段初的蓄水状态和入库流量信息，计算求得所述一个水库所述某一时段的最小出力对应的出库流量，并与最小下泄流量比较，取两者中较大者作为所述一个水库所述某一时段的出库流量，根据水量平衡方程获得所述一个水库所述某一时段末的蓄水状态，并将其作为所述一个水库所述某一时段的下一时段初的蓄水状态；

(3-2-2)按所述步骤(3-2-1)从调度期初递推计算至调度期末，得到一条蓄水过程线，即为所述一个水库的最小决策限制线；

(3-2-3)按所述步骤(3-2-1)和(3-2-2)的方法求得所述各水库的最小决策限制线；

(3-3)对所述一个水库相同时刻的所述防弃水线和所述最小决策限制线的蓄水状态值取小，获得一条蓄水过程线，作为所述一个水库的初始参考线；

(3-4)按步骤(3-3)的方法求得所述各水库的初始参考线。

优选地，所述步骤(5)中所述各水库K值计算包括以下子步骤：

(5-1)确定所述各水库K值计算的起始时段和终止时段；所述各水库K值计算的起始时段取待调整时段t₀；所述各水库K值计算的终止时段的确定包括如下步骤：

(5-1-1)计算所述各水库t₀时段增发发单位出力后的时段末蓄水状态；根据所述各水库t₀时段初的蓄水状态、入库流量和出力值以及水量平衡方程计算所述各水库t₀时段末的蓄水状态。

(5-1-2)根据步骤(5-1-1)所得所述各水库t₀时段末的蓄水状态，按步骤(3-2)的方法重新计算所述各水库从t₀+1时段末至调度期末的最小决策限制线，设其与所述各水库的防弃水线相交于时段t₁，则时段t₁即为所述各水库K值计算的终止时段；

(5-2)计算所述各水库t₀时段增发单位出力后在t₀时段至t₁时段内总的不蓄电能损失值，即为所述各水库的K值；所述各水库中任一水库某一时段的不蓄电能损失值为所述任一水库所述某一时段的水头减少量与入库流量的乘积。

总体而言，按照本发明的上述技术构思与现有技术相比，主要具备以下的技术优点：

1、本发明将水库最小决策限制线和防弃水线结合使用，可达到在尽量不弃水的情况下尽可能抬高水电站发电水头的效果，提高水电站能量转换效率；

2、本发明所涉及的方法体现了专家调度经验，在保证跳读方案最优性的同时可大大减小计算工作量，且所得结果可理解性强，便于调度人员在实际运用中动态决策；

3、本发明与水库蓄供水判别系数法相比，其优势在于可准确确定水库增发单位出力引起的不蓄电能损失值(K值)计算的终止时刻；

4、本发明与传统的优化方法(动态规划及其改进算法)相比，其优势在于不存在水库群调度“维数障碍”问题，适用于大规模并联水库群优化调度。

附图说明

图1是本发明基于参考线的并联水库确定性优化调度的总体技术路线；

图2是本发明基于参考线的并联水库确定性优化调度流程图；

图3是采用0.618法以不蓄电能损失最小原则将不足出力分配至各水库的流程图；

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。

水库发电过程中涉及到的能量有水库初始蓄能、输入能、损耗能和输出能，根据能量守恒定律，在水库初始蓄能和输入能一定的情况下，损耗能越少，则输出的电能越多。基于上述原理，可获得水库发电优化调度须遵循的两条基本原则：(1)尽量避免弃水，以减少弃水能损耗；(2)尽量抬高水库运行水位，以减少不蓄电能损失。根据水库初始蓄水位和计划期内入库流量过程，运用上述原则可确定水库最优运行策略。

图2所示为基于参考线的并联水库确定性优化调度方法的流程图，具体包括以下步骤：

步骤1设置调度期，将调度期(1年)均匀划分为T(T为正整数，在本发明实施例中，T＝12或36)个等间距的时段；设置水库调度期内的约束条件，包括蓄水量约束：出力约束：出库流量约束：水量平衡方程：V_m,t+1＝V_m,t+(I_m,t-R_m,t)·Δt和水库期初、期末水位限制：Z_m,1＝Z_ms，Z_m,T+1＝Z_ms。其中，V _m,t、V_m,t、 N _m,t、N_m,t、 R _m,t、R_m,t、分别表示t(t＝1,2，…，T)时段水库m的最小蓄水量、蓄水量、最大蓄水量、最小出力、出力、最大出力、最小下泄流量、下泄流量和最大下泄流量；I_m,t表示m水库t时段的入库流量；Δt表示时段步长；Z_ms、Z_mend分别为水库m期初和期末水位，为给定值。

步骤2根据水库的库容特性曲线、尾水关系曲线、以及水库优化调度的目标函数和上述约束条件，以水库各时段入库流量作为输入，建立水库群长期优化调度发电量最大模型如下：

$E=max\underset{m=1}{\overset{M}{\Σ}}\underset{t=1}{\overset{T}{\Σ}}{N}_{m,t}\·\mathrm{\Δ}t$

N_m,t＝K_mQ_m,tH_m,t

其中，E表示水库群全时期最优发电量；N_m,t、Q_m,t和H_m,t分别表示t时段水库m的出力、发电流量和发电水头；K_m表示水库m的综合出力系数。

步骤3确定水库群初始参考线，包括确定各水库防弃水线和最小决策限制线，优选地，步骤3包括以下子步骤：

(3-1)确定水库防弃水线。水库运行过程中，减少弃水有利于提高发电量。水库用于发电的下泄水量在以下两种情况下会形成弃水：一是下泄流量大于水库水轮机最大过流量PlantsMaxQ_m；二是水库按PlantsMaxQ_m下泄，其出力值大于时段最大出力限制具体包括以下子步骤：

(3-1-1)已知水库m计划期末水位Z_m,T+1和时段入库流量I_m,T,由库容特性线性插值得其时段末库容V_m,T+1，令所述时段其出库流量R_m,T等于水库水轮机最大过流PlantsMaxQ_m，由V_m,T+1、R_m,T、I_m,T和水量平衡方程，可求得水库时段初库容V_m,T，若则取若V_m,T＜V _m,T，则取V_m,T＝V _m,T，求得时段初库容V_m,T后由库容特性线性插值获得时段初水位Z_m,T；

(3-1-2)由V_m,T、V_m,T+1、I_m,T、PlantsMaxQ_m、和水量平衡方程可求得弃水流量若则调整出库流量为根据水量平衡方程重新计算V_m,T，若则取若V_m,T＜V _m,T，则取V_m,T＝V _m,T，求得V_m,T后由库容特性线性插值获得时段初水位Z_m,T；若则所述时段初水位Z_m,T取上述步骤(3-1-1)所获得的值；将Z_m,T作为上一时段的时段末水位；

(3-1-3)按照步骤(3-1-1)和步骤(3-1-2)从调度期末递推计算至调度期初，得到一条水位过程线(t＝1,2，…，T，T+1)，相应的库容为(t＝1,2，…，T，T+1)；

(3-1-4)对步骤(3-1-3)中所述水位过程线从调度期初至调度期末逐时段检查是否有弃水发生。若各时段均无弃水发生，则所述水位过程线即为水库防弃水线；若时段t₁有弃水发生，则进行相应调整，调整原则是使弃水集中在最后一个或若干个时段，以保持水库高水位发电，减小不蓄电能损失。具体调整方法如下：调整水库t₁时段出库流量等于水库水轮机最大过流PlantsMaxQ_m，由和水量平衡方程，可求得水库时段末库容若则取 $V_{m,t_1+1}={\stackrel{\&OverBar;}{V}}_{m,t_1+1};$ 若 $V_{m,t_1+1}<{\underset{\&OverBar;}{V}}_{m,t_1+1},$ 则取 $V_{m,t_1+1}={\underset{\&OverBar;}{V}}_{m,t_1+1},$ 由PlantsMaxQ_m、和水量平衡方程可求得弃水流量若则调整出库流量为根据水量平衡方程重新计算若 $V_{m,t_1+1}>{\stackrel{\&OverBar;}{V}}_{m,t_1+1},$ 则取 $V_{m,t_1+1}={\stackrel{\&OverBar;}{V}}_{m,t_1+1};$ 若 $V_{m,t_1+1}<{\underset{\&OverBar;}{V}}_{m,t_1+1},$ 则取 $V_{m,t_1+1}={\underset{\&OverBar;}{V}}_{m,t_1+1},$ 求得后由库容特性线性插值获得时段末水位并计算所述时段较调整前增加的蓄水量后续时段水库末库容均增加若遇t₂(t₁＜t₂＜T)时段末库容大于最大库容的情况，则令计算所述时段增加的蓄水量后续时段水库末库容在的基础上均增加调整后的水位过程线即为防弃水线；

(3-1-5)按所述步骤(3-1-1)～(3-1-4)的方法求得所有水库的防弃水线；

(3-2)确定最小决策限制线。水电站运行过程应尽量抬高水位发电，以减小不蓄电能损失，增发电量，具体包括以下子步骤：

(3-2-1)已知水库调度期初水位Z_m,1和时段入库流量I_m,1，计算所述时段水库按最小出力N _m,1运行对应的出库流量并与所述时段容许的最小下泄流量R _m,1比较，取两者中取较大者作为所述时段出库流量R_m,1，即所述时段出库流量根据水量平衡方程，计算所述时段末库容V_m,2,若则取若则取V_m,2＝V _m,2，再根据库容特性线性插值得所述时段末水位Z_m,2,将其作为下一时段的时段初水位；

(3-2-2)按步骤(3-2-1)从调度期初递推计算至调度期末，得到一条水位过程线(t＝1,2，…，T,T+1)，即为水库m最小决策限制线。

(3-2-3)按所述步骤(3-2-1)和(3-2-2)的方法求得所有水库的最小决策限制线；

(3-3)对水库相同时刻的防弃水线和最小决策限制线的水位值取小，得到水库初始参考线(t＝1,2，…，T,T+1)，即 $Z_{m,t}^{refer0}=min(Z_{m,t}^{unsp},Z_{m,t}^{\min D});$

(3-4)按步骤(3-3)的方法求得所有水库初始参考线；

步骤4以各水库初始参考线为其初始调度线，根据各水库各时段入库流量I_m,t和水量平衡方程从前往后逐时段计算各水库各时段出库流量R_m,t和出力N_m,t，直至计算完调度期最后一个时段为止。若时段t遇则取所有水库计算完毕后统计水库群各时段总出力值

步骤5从前往后逐时段检查水库群时段总出力值是否大于保证出力N_b，是，则转下一时段；否，则计算不足出力N_bz并按水库群不蓄电能损失最小原则将不足出力N_bz分配至各水库，调整完毕后根据调整后的水库时段末蓄水状态按步骤(3-2)所述方法调整后续时段各水库的最小决策限制线，调整完毕后转下一时段。上述过程一直持续至调度期末。根据等微增率原理，调整后各水库增发单位出力引起的不蓄电能损失值(K值)应大致相当。采用0.618法调整水库群初始参考线，以待调整时段t₀各水库调整出力之后的K值为搜索变量，以待调整时段t₀各水库调整后的出力总和等于保证出力为搜索条件，优选地，包括以下步骤：

(5-1)确定K值的下限值。包括以下步骤：

(5-1-1)计算各水库增发单位出力引起的不蓄电能损失值K_m，包括以下步骤：

(5-1-1-1)将水库待调整时段t₀的出力增加一个单位值，计算出时段末水位，作为下一时段的时段初水位，按所述步骤(3-2)的方法计算待调整时段t₀至调度期最后一个时段T的水库最小决策限制线设其与所述水库防弃水线相交于时刻t₁，记t₁时刻所在时段为t₂。

(5-1-1-2)计算所述t₀～t₂时期内水库各时段水头减少量ΔH_m,t和不蓄电能损失值K_m。水库m时段t的水位减少量应为调整前后水库平均水位之差，即 ${\mathrm{\&Delta;H}}_{m,t}=0.5\&CenterDot;(Z_{m,t}^{\min D}+Z_{m,t+1}^{\min D})-0.5\&CenterDot;(Z_{m,t}^{\&prime;\min D}+Z_{m,t+1}^{\&prime;\min D}).$ 水库m时段t内不蓄电能损失值K_m,t应为所述时段入库流量I_m,t与水头减少量ΔH_m,t之乘积，即K_m,t＝I_m,t·ΔH_m,t。

(5-1-1-3)按所述步骤(5-1-1-1)和(5-1-1-2)的方法计算所有水库增发单位出力引起的不蓄电能损失值K_m。

(5-1-2)K值的下限值取各水库增发单位出力引起的不蓄电能损失值中的最小者，即K_min＝min(K_m),其中(m＝1,2,…,M)。

(5-2)确定K值的上限值。包括以下步骤：

(5-2-1)计算各水库减发单位出力引起的不蓄电能增加值K'_m，包括以下步骤：

(5-2-1-1)取水库待调整时段t₀末的水位为水库死水位根据水量平衡方程计算所述水库待调整时段t₀的出库流量R_m,t、出力N_m,t和弃水流量若则减小出库流量R_m,t，使弃水流量并计算时段末水库水位值Z_m,t+1和时段出力值N_m,t；否则，保持出库流量R_m,t、时段末水位Z_m,t+1和时段出力值N_m,t不变。按所述步骤(3-2)的方法计算待调整时段t₀至调度期最后一个时段T的最小决策限制线设其与所述水库防弃水线相交于时刻t₃，记t₃时刻所在时段为t₄。

(5-2-1-2)将水库待调整时段t₀的出力减小一个单位值，计算出时段末水位Z_m,t+1，作为下一时段初水位，按所述步骤(3-2)的方法计算所述t₀～T时期内的水库最小决策限制线与所述水库防弃水线相交于时刻t₅，记t₅时刻所在时段为t₆。

(5-2-1-3)计算所述t₀～t₆时期内水库各时段水头增加量ΔH'_m,t和不蓄电能增加值K'_m。水库m时段t的水位增加量应为调整后与调整前水库平均水位之差，即 ${\mathrm{\&Delta;H}}_{m,t}=0.5\&CenterDot;(Z_{m,t}^{\&prime;\&prime;\&prime;\min D}+Z_{m,t+1}^{\&prime;\&prime;\&prime;\min D})-0.5\&CenterDot;(Z_{m,t}^{\&prime;\&prime;\min D}+Z_{m,t+1}^{\&prime;\&prime;\min D}).$ 水库m时段t内不蓄电能增加值K'_m,t应为所述时段入库流量I_m,t与水头增加量ΔH'_m,t之乘积，即K'_m,t＝I_m,t·ΔH'_m,t。

(5-2-1-4)按所述步骤(5-2-1-1)和(5-2-1-3)的方法计算所有水库增发单位出力引起的不蓄电能增加值K'_m。

(5-2-2)取K值的上限值为各水库减发单位出力引起的不蓄电能增加值K'_m中的最大者，即K_max＝max(K'_m),其中(m＝1,2,…,M)。

(5-3)采用0.618法在闭区间[K_min,K_max]内搜索得到一个K^*，使待调整时段t₀各水库调整出力之后的K值与K^*之差的绝对值最小，即min(K-K^*)，且使水库群出力总和等于保证出力N_b，即包括以下步骤：

(5-3-1)在闭区间[K_min,K_max]内取K₁＝K_min+0.382(K_max-K_min),各水库以待调整时段t₀调整出力之后的K值与K₁之差的绝对值最小为目标,即以min(K-K₁)为目标，采用0.618法搜索得到其出力值,所有水库搜索完毕后计算待调整时段t₀水库群的总出力值，记为N₁，并保存各水库调整之后的参考线取K₂＝K_min+0.618(K_max-K_min),各水库以待调整时段t₀调整出力之后的K值与K₂之差的绝对值最小为目标，即以min(K-K₂)为目标，采用0.618法搜索得到其出力值,所有水库搜索完毕后计算待调整时段t₀水库群的总出力值，记为N₂，并保存各水库调整之后的参考线若|N₁-N_b|＜|N₂-N_b|,则令K_max＝K₂，否则令K_min＝K₁。

所述步骤(5-3-1)中采用0.618法搜索水库出力值，包括以下步骤(以水库m和K＝K₁为例)：

(5-3-1-1)确定水库m出力上限和下限下限值取初始方案中水库待调整时段t₀的出力，上限值取所述水库所述时段最大容许出力

(5-3-1-2)在闭区间内在取并令水库待调整时段t₀的出力等于N₃，即计算时段末水库水位值并按(5-1-1)所述方法计算所述水位值下水库增发单位出力引起的不蓄电能损失值K₃；取并令水库待调整时段t₀的出力等于N₄，即计算时段末水库水位值并按(5-1-1)所述方法计算所述水位值下水库增发单位出力引起的不蓄电能损失值K₄。若|K₃-K₁|＜|K₄-K₁|,则令否则令

(5-3-1-3)重复步骤(5-3-1-2)，直至满足精度要求为止，即输出水库m调整之后的出力值N'_m,t；

(5-4)重复步骤(5-3)，直至满足精度要求为止，即|K_max-K_min|＜ε'，输出K^*；

步骤(6)输出各水库调整后的出力值N'_m,t和调整后参考线

本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (3)

1.一种基于参考线的并联水库确定性优化调度方法，其特征在于，所述方法包括以下步骤：

步骤(1)，设置调度期，将所述调度期均匀划分为T个等间距的时段,其中T为正整数；设置所述并联水库调度期内的初始条件和约束条件；所述初始条件包括水库调度期期初蓄水状态、期末蓄水状态和调度期内各时段入库流量；所述约束条件包括水库蓄水量约束、出力约束、出库流量约束以及水量平衡方程；

步骤(2)，根据所述并联水库中各水库的库容特性曲线、尾水关系曲线、水库优化调度的目标函数和所述约束条件，建立水库群长期优化调度发电量最大模型：

$E=max\underset{m=1}{\overset{M}{\&Sigma;}}\underset{t=1}{\overset{T}{\&Sigma;}}{N}_{m,t}\&CenterDot;\mathrm{\&Delta;}t$

N_m,t＝K_mQ_m,tH_m,t

其中，E表示水库群全时期最优发电量；N_m,t、Q_m,t和H_m,t分别表示t时段水库m的出力、发电流量和发电水头；K_m表示水库m的综合出力系数；

步骤(3)，确定所述各水库的初始参考线，包括确定所述各水库的防弃水线和最小决策限制线；

步骤(4)，以所述各水库的初始参考线为其初始调度线，根据入库流量信息和水量平衡方程从前往后逐时段递推至最后一个时段，计算所述各水库各时段的出库流量和出力值，直至所有水库计算完毕,统计所述并联水库各时段的总出力值；

步骤(5)，从前往后逐时段检查所述并联水库各时段的总出力值是否大于保证出力，是，则转下一时段；否，则计算不足出力并按所述并联水库不蓄电能损失最小原则调整所述各水库的时段出力值，调整完毕后根据调整后的所述各水库的蓄水状态按所述步骤(3)所述方法调整后续时段所述各水库的最小决策限制线，调整完毕后转下一时段；上述过程一直持续至调度期末；根据等微增率原理，调整后所述各水库增发单位出力引起的不蓄电能损失值，即K值，应大致相当；以待调整时段t₀所述各水库调整出力之后的K值为搜索变量，以t₀时段各水库调整之后的出力总和等于保证出力为搜索条件采用一维搜索方法调整所述各水库的初始参考线；

步骤(6)，输出所述各水库调整后的参考线和各时段出力值。

2.如权利要求1所述的方法，其特征在于，所述步骤(3)包括以下子步骤：

(3-1)以减少弃水能为原则确定所述各水库的防弃水线，其包括以下子步骤：

(3-1-1)已知所述并联水库中的某一水库一个时段的时段末蓄水状态和入库流量信息，令所述时段出库流量等于水轮机最大过流，由水量平衡方程求得所述某一水库所述时段初的蓄水状态；

(3-1-2)根据所述某一水库所述时段初的蓄水状态、时段末的蓄水状态、入库流量和最大容许出力，计算求得所述某一水库的弃水流量，若所述弃水流量大于零，则减小所述某一水库的出库流量使得所述弃水流量为零，并根据所述水量平衡方程重新计算所述时段初的蓄水状态，若所述弃水流量等于零，则所述时段初的蓄水状态取所述步骤(3-1-1)所求得的值，将计算得到的所述时段初的蓄水状态作为上一时段末的蓄水状态；

(3-1-3)按所述步骤(3-1-1)和(3-1-2)从调度期末递推计算至调度期初，得到所述某一水库的一条蓄水过程线；

(3-1-4)对(3-1-3)所述蓄水过程线从调度期初至调度期末逐时段检查是否有弃水发生，若各时段均无弃水，则所述蓄水过程线即为所述某一水库的防弃水线；若某时段有弃水发生，则进行相应调整，使弃水集中在最后一个或若干个时段，调整后的蓄水过程线即为所述某一水库的防弃水线；

(3-1-5)按所述步骤(3-1-1)～(3-1-4)的方法求得所述各水库的防弃水线；

(3-2)根据保持高水位运行以减少不蓄电能损失的原则，确定所述各水库的最小决策限制线，其包括以下子步骤：

(3-2-1)根据所述并联水库中的一个水库某一时段初的蓄水状态和入库流量信息，计算求得所述一个水库所述某一时段的最小出力对应的出库流量，并与最小下泄流量比较，取两者中较大者作为所述一个水库所述某一时段的出库流量，根据水量平衡方程获得所述一个水库所述某一时段末的蓄水状态，并将其作为所述一个水库所述某一时段的下一时段初的蓄水状态；

(3-2-2)按所述步骤(3-2-1)从调度期初递推计算至调度期末，得到一条蓄水过程线，即为所述一个水库的最小决策限制线；

(3-2-3)按所述步骤(3-2-1)和(3-2-2)的方法求得所述各水库的最小决策限制线；

(3-3)对所述一个水库相同时刻的所述防弃水线和所述最小决策限制线的蓄水状态值取小，获得一条蓄水过程线，作为所述一个水库的初始参考线；

(3-4)按步骤(3-3)的方法求得所述各水库的初始参考线。

3.如权利要求1所述的方法，其特征在于，所述步骤(5)中所述各水库K值计算包括以下子步骤：

(5-1)确定所述各水库K值计算的起始时段和终止时段；所述各水库K值计算的起始时段取待调整时段t₀；所述各水库K值计算的终止时段的确定包括如下步骤：

(5-1-1)计算所述各水库t₀时段增发发单位出力后的时段末蓄水状态；根据所述各水库t₀时段初的蓄水状态、入库流量和出力值以及水量平衡方程计算所述各水库t₀时段末的蓄水状态。

(5-1-2)根据步骤(5-1-1)所得所述各水库t₀时段末的蓄水状态，按步骤(3-2)的方法重新计算所述各水库从t₀+1时段末至调度期末的最小决策限制线，设其与所述各水库的防弃水线相交于时段t₁，则时段t₁即为所述各水库K值计算的终止时段；

(5-2)计算所述各水库t₀时段增发单位出力后在t₀时段至t₁时段内总的不蓄电能损失值，即为所述各水库的K值；所述各水库中任一水库某一时段的不蓄电能损失值为所述任一水库所述某一时段的水头减少量与入库流量的乘积。