一种水电站群期末蓄能最大短期优化调度方法技术

本发明专利技术公开了一种水电站群期末蓄能最大短期优化调度方法，属于水电调度领域。首先将原模型中各时段总出力与总负荷需求的一致性约束，代替为各时段总出力与总负荷需求最大差值为零的约束；而后借助凝聚函数法将简化后的最大差值转换为便于处理的形式。在此基础上，采用拉格朗日松弛法建立对偶优化模型，并以次梯度法求解。本发明专利技术方法将存在规模为96的拉格朗日乘子向量的期末蓄能最大优化调度问题简化为只有一个拉格朗日乘子的问题，降低了问题的复杂性，采用次梯度法和逐次逼近搜索算法相结合的方法可有效求解，具有很强的实用性，应用广泛。

【技术实现步骤摘要】

本专利技术属于水电调度领域，是一种涉及复杂约束条件下，水电站群期末蓄能最大短期优化调度方法。
技术介绍
水电站短期优化调度一般指以15分钟为时段，未来几天内的优化调度问题，可以大致分为两类，第一类是给定各水库末水位，在调度期内按照调度目标调节水量时空分布，此类模型偏重于在满足各电站的中长期控制方式前提下，提高水能利用效率，是广义的以水定电方法，目标通常为发电量最大、发电效益最大等；第二类是给定各时段水电系统负荷需求，在各电站间分配出力，是广义的以电定水方法，主要采用蓄能最大等模型，适于梯级和跨流域水电站群负荷分配问题。蓄能最大模型求解受到负荷需求约束条件限制，较困难。目前流行的拉格朗日松弛算法求解效率较低，一方面由于短期调度每天包括96个时段，拉格朗日乘子向量规模较大，需进行很多次向量更新；另一方面，由于对偶问题包括出力持续性、爬坡速度等复杂约束，因而求解耗时较多。
技术实现思路
针对以拉格朗日松弛法求解水电站群短期蓄能最大模型时，乘子向量规模较大对计算效率的限制问题，提出:首先以各时段总出力与其目标间最大差值为零代替各时段总出力目标约束；而后借助凝聚函数法将简化后的最大差值进行转换。在此基础上，采用拉格朗日松弛法建立对偶模型，并采用次梯度法和逐次逼近搜索算法相结合的方法。本专利技术的技术方案如下:一种水电站群期末蓄能最大优化调度方法，以15分钟为时段，进行未来I天内的梯级和跨流域水电站群负荷分配；首先将原模型中各时段总出力与总负荷需求的一致性约束，代替为各时段总出力与总负荷需求最大差值为零的约束；而后借助凝聚函数法将简化后的最大差值进行转换，问题复杂度较原问题大大降低。最后模型采用拉格朗日松弛和逐次逼近搜索算法求解。本专利技术首先将原模型中各时段总出力与总负荷需求的一致性约束，代替为各时段总出力与总负荷需求最大差值为零的约束，而后借助凝聚函数法将简化后的最大差值转换为便于处理的形式。该方法大大降低了问题的复杂性，采用次梯度法和逐次逼近搜索算法相结合的方法可有效求解，具有很强的实用性，应用广泛。附图说明图1是初始发电流量过程；图2是扩展调整范围过程示意；图3是流量变化速度修正示意；图4是流量波动修正示意；图5是出力升降时段修正示意；图6是各电站出力过程；图7是系统总出力过程具体实施例方式以下结合技术方案详细叙述本专利技术的具体实施例。(一)建立短期蓄能最大优化调度模型1.1目标函数本文档来自技高网...

【技术保护点】
一种水电站群期末蓄能最大优化调度方法，其特征包括以下步骤，(一)建立短期蓄能最大优化调度模型1.1?目标函数MaxF=Σm=1MESm+Σm=1MElm---(1)其中：M为跨流域水电站群参与计算电站的总数，1≤m≤M；t为时段标号，T为计算时段总数，控制期内1≤t≤T；ESm表示m号电站及其全部上游电站死水位以上水量在m号电站可产生的电量，采用式（2）计算；Elm表示m号电站上游电站在m号电站产生的滞时电量，采用式（3）计算；ESm={[VmT+WT(m)]/ηm}WT(m)=Σk=1Km{VUm[k]T+WT(Um[k])}---(2)其中：表示水库m第t时段初的库容；ηm表示m号电站平均耗水率；WT(m)表示m号电站水库及其全部上游电站水库调度期末死水位以上蓄水量；Um表示m号电站直接上游电站标号数组；Km为m号水电站水库的直接上游电站数目；Elm=WL(m)/ηm---(3)其中：，WL(m)表示m号电站全部上游电站滞时水量计算函数；，表示有直接上下游关系的上游Um[k]号水电站水库和下游m号水电站水库间的最大最小滞时时段数；Δt表示第t时段的小时数；表示Um?[k]号电站在n时段出库流量为时，在下游m号电站T?1时段以后产生的入库流量，按式（4）计算：为Um?[k]号电站水库出库流量为时的滞时时段数；1.2?约束条件除常规水量平衡、?发电流量约束、电站出力约束、库水位约束、出库流量约束等条件外，该模型尚需满足如下约束：水电站出力爬坡限制：?（MW）表示m号水电站相邻时段最大出力升降限制；水电站出力波动限制：????????????????????????（6）tvm是m号水电站出力升降最小间隔时段数，即在一轮出力升降过程最高最低点需持续最少tvm个时段；最小出力升降时段数限制：m号水电站出力上升开始至下降开始，或由下降开始至上升开始的时间间隔不少于tpm个时段；水电总负荷约束：?其中：Nt表示第t时段的水电负荷；（二）短期蓄能最大优化调度方法2.1模型简化将约束，t=1,2,…,T，转换为：；而后采用凝聚函数法将转换为，；P为参数，取为0.01；2.2建立对偶问题针对水电系统总负荷约束采用拉格朗日松弛法建立对偶优化问题；在不采用2.1模型时，建立以式（8）为目标函数，以水量平衡、发电流量约束、电站出力约束、库水位约束、出库流量约束、水电站出力爬坡限制、最小出力升降时段数限制为约束条件的优化问题；minBmaxAF1=F+Σt=1T[λt(Σm=1Mpmt-Nt)]---(8)λt为t时段拉格朗日乘子向量分量，A=p11,p12,...,p1Tp21,p22,...,p2T...pM1,pM2,...,pmT，B={λ1,λ2,…,λT}；而在采用简化方式时，对偶优化问题约束条件相同，目标函数为：minγmaxAF2=F+γ[1Pln{Σt=1TeP[f(t)-max1≤t0≤Tf(t0)]}+max1≤t≤Tf(t)],f(t)=|Σm=1Mpmt-Nt|---(9)γ为拉格朗日乘子；2.3优化求解计算流程如下：（1）初设i=1和拉格朗日乘子γi；（2）采用逐次逼近的搜索算法求解以式（9）为目标函数，以水量平衡、发电流量约束、电站出力约束、库水位约束、出库流量约束、水电站出力爬坡限制、最小出力升降时段数限制为约束条件的优化问题，此时；（3）采用次梯度法更新乘子，即γi+1=γi+β*F2，?β为乘子更新系数，设i=i+1；（4）若|γi+1?γi|＜ε，ε为收敛精度，转步骤（5），否则返回步骤（2）；（5）在优化结果的基础上，以为目标消除对偶间隙，仍采用逐次逼近的算法求解；流程中采用的逐次逼近的搜索算法需要使得搜索模式满足出力爬坡、波动和升降时段约束，本专利技术采用如下方法：（1）在t时段发电流量由变为后，依次向前后方向扩展，以满足出力波动约束；计原具体方法是依次向前调整t?1，t?2，t?tvm/2时段的出力，直到；向后搜索的方法类似；（2）在扩展范围边缘修正出力变化幅度，以满足出力爬坡约束，设扩展范围的范围为[t1,t2]，则首先搜索最大的满足，（t0＜t1）的时段t0，由时段依次按照爬坡出力值调整t1?1，t1?2，…t0+1的发电流量；而后搜索最小的满足，（t3＞t2）的时段t3，而后由时段依次按照爬坡出力值调整t2+1，t2+2，…t3?1的发电流量；（3...

【技术特征摘要】
1.一种水电站群期末蓄能最大优化调度方法，其特征包括以下步...

【专利技术属性】
技术研发人员：何光宏，戴建炜，武新宇，林成，廖胜利，汪明清，苏华英，王俊莉，
申请(专利权)人：贵州乌江水电开发有限责任公司，大连理工大学，
类型：发明
国别省市：