计及风电出力不确定性的抽水蓄能-风电联合调度方法
【专利摘要】本发明公开了一种计及风电出力不确定性的抽水蓄能?风电联合调度方法，步骤包括：针对含抽蓄电站、常规机组及风电场的发电系统制定抽水蓄能?风电联合运行计划，建立基于常规机组的燃料费用和启停成本的目标函数并确定约束条件，生成常规机组所有可能出力的集合并作为搜索空间，在搜索空间中生成位置随机分布的决策变量，将决策变量作为国家采用帝国竞争算法求解结果来控制发电系统各常规机组的出力。本发明针对风电出力的不确定性引入旋转备用风险约束，精细化考虑了风电出力不确定性的影响，能够实现含风电发电系统运行经济性与可靠性的统筹协调，将抽水蓄能电站视为平抑风电出力波动性的技术手段与风电场联合运行，能够显著提高系统的经济效益。
【专利说明】
计及风电出力不确定性的抽水蓄能-风电联合调度方法
技术领域
[0001] 本发明涉及抽水蓄能-风电联合调度方法，具体涉及一种计及风电出力不确定性 的抽水蓄能-风电联合调度方法。
【背景技术】
[0002] 风电是一种无污染、绿色的可再生能源，其分布广泛、能量密度高、利用效率高，适 合大规模开发，因此，风力发电技术已受到了世界各国的高度重视。然而，由于风电出力具 有随机性与波动性的特点，其大规模接入电网将会给传统的调度方法带来极大的挑战，因 此，研究考虑风电接入的发电系统日前调度方法具有重要的理论价值与现实意义。目前很 多专家学者从不同角度研究了风电接入的发电系统日前调度问题，普遍考虑了风电出力不 确定性对系统的影响，但是从运行模式而言，均采用传统火电机组作为平抑风电出力波动 性的手段，然而，传统火电机组调峰限制因素较多，且成本较高，仅利用其与风电联合运行， 很难获得较低的运行成本和理想的平抑效果。大规模电力储能技术可望成为未来平抑风电 出力波动性和随机性的一种有效技术手段。而抽水蓄能电站作为现阶段较为成熟的大规模 储能技术之一，如何利用其与风电联合运行，缓解风电大规模并网给发电系统日前调度带 来的不利影响，已成为近年来人们的研究热点。因此，如何考虑风电出力不确定性，解决各 台发电机组之间的负荷经济分配问题，已经成为一项亟待解决的关键技术问题。
【发明内容】
[0003] 本发明要解决的技术问题:针对现有技术的上述问题，提供一种能够有效平抑风 电的波动性、降低系统的旋转备用需求、显著提高系统经济效益、实现了含风电电力系统运 行经济性与可靠性的统筹协调、调度求解快速高效的计及风电出力不确定性的抽水蓄能-风电联合调度方法。
[0004] 为了解决上述技术问题，本发明采用的技术方案为：
[0005] -种计及风电出力不确定性的抽水蓄能-风电联合调度方法，步骤包括：
[0006] 1)针对含抽蓄电站、常规机组及风电场的发电系统，将抽蓄电站作为平抑风电出 力波动性的技术手段控制抽蓄电站抽水工况与发电工况间的转换与风电场联合运行，制定 抽水蓄能-风电联合运行计划，建立基于常规机组的燃料费用和启停成本的目标函数，确定 发电系统的约束条件；
[0007] 2)根据抽水蓄能-风电联合运行计划生成发电系统中各常规机组所有可能出力的 集合；
[0008] 3)将发电系统中各常规机组所有可能出力的集合作为搜索空间，在搜索空间中生 成位置随机分布的多个指定维数的决策变量，将决策变量作为帝国竞争算法的国家、根据 目标函数计算国家的势力、采用帝国竞争算法求解得到最优的一个决策变量作为结果来控 制发电系统中各常规机组的出力。
[0009] 优选地，所述步骤1)中的目标函数如式(1)所示；
[0010]
⑴ f=l n~i
[0011] 式(1)中，/?)表示目标函数，T为发电系统的调度周期，N为发电系统中发电常规 机组数，Cn表示常规机组η的运行成本系数，K表示常规机组η在t时刻的出力，u nt表示常规 机组η在t时刻的状态，iMt-υ表示常规机组η在t-Ι时刻的状态，unt值为0表示常规机组η处于 停机状态，u nt值为1表示常规机组η处于开机状态，Snt为常规机组η在t时刻的启动费用。
[0012] 优选地，所述步骤1)中确定发电系统的约束条件包括如式(2)所示不计网损的情 况下的系统功率平衡约束；
[0013]
C2)
[0014] 式(2)中，N为发电系统中发电常规机组数，<表示常规机组η在t时刻的出力，unt表 示常规机组η在t时刻的状态，g表示在t时段风电的等效出力，/匕表示第t时段的负荷预测 值。
[0015] 优选地，所述步骤1)中确定发电系统的约束条件还包括如式(3)所示的常规机组 出力约束；
[0016] ΡΓ < P；, < P^x (3)
[0017] 式⑶中，€表示常规机组η在t时刻的出力，iTm表示常规机组η的最小出力，表 示常规机组η的最大出力。
[0018] 优选地，所述步骤1)中确定发电系统的约束条件还包括如式(4)所示的上升爬坡 速率约束、如式(5)所示的下降爬坡速率约束；
[0019] (4)
[0020] (5)
[0021] 式(4)和式(5)中，Κ表示常规机组η在t时刻的出力，/Γ1表示常规机组η在t-Ι时刻 的出力，Unt表示常规机组η在t时刻的状态，UnU-υ表示常规机组η在t-ι时刻的状态，DRn表示 常规机组η的下降爬坡速率，U Rn表示常规机组η的上升爬坡速率，△ T表示单位调度时段。
[0022] 优选地，所述步骤1)中确定发电系统的约束条件还包括如式(6)所示的抽蓄电站 库容约束、如式(7)所示的抽蓄电站输出功率约束、如式(8)所示的抽蓄电站日抽发电量约 束；
[0023]
(6)
[0024] 式(6)中，W。为抽蓄电站上水库初始水量，WmaxS抽蓄电站上水库最大水量，qs为平 均水量转换系数，T为发电系统的调度周期，4为时段t抽蓄电站的发电功率，％为平均电量 转换系数，Ps为抽蓄电站的抽水功率，6为时段t抽蓄电站处于抽水工况的机组数，W min为抽 蓄电站上水库最小水量；
[0025]
(7、
[0026] 式(7)中，舄表示时段t抽蓄电站处于发电工况的机组数，PC,min表示抽蓄电站发电 机组的最低出力，疼为时段t抽蓄电站的发电功率，P h表示抽蓄电站发电机组的额定功率， 苟为时段t抽蓄电站处于抽水工况的机组数；
[0027] (8)
[0028] JM? We/、，」拽;ti由蓄电站的发电功率，nc为平均电量转换系数，ns为平均水量 转换系数，τ为调度周期T内某一时刻，T为发电系统的调度周期，Ps为抽蓄电站的抽水功率， 端为时段t抽蓄电站处于抽水工况的机组数。
[0029] 优选地，所述步骤1)中确定发电系统的约束条件还包括如式(9)所示的正旋转备 用风险指数约束、如式（10)所示的负旋转备用风险指数约束，所述正旋转备用风险指数为 风电的实际出力少于风电平均出力导致的正旋转备用不足的概率，所述负旋转备用风险指 数为风电的实际出力大于平均出力导致的负旋转备用不足的概率；
[0030] (9)
[0031] (10)
[0032]式(9)和式（10)中，Qd为正旋转备用风险指数，Qu为负旋转备用风险指数，β为旋转 备用风险阈值。
[0033] 优选地，所述正旋转备用风险指数的函数表达式如式(11)所示；
[0034] ft = q{Pw I {Pw < P，r)&(RGp < P"--Plv+Rfp)} ( 11 )
[0035] 式（11)中，Qd为正旋转备用风险指数，q{}表示大括号{}中事件发生的概率，大括 号U中事件为风电的实际出力少于风电平均出力导致的正旋转备用不足的事件，Pw为风机 的输出功率，落为风电平均出力，R Gp为发电系统可以提供的正旋转备用，RLp为发电系统中 常规机组的正旋转备用需求。
[0036] 优选地，所述负旋转备用风险指数的函数表达式如式(12)所示；
[0037] β( = q{Pw I (P"- >P,r)&(Ra, < Pw -Pw + RLh)} (]2)
[0038] 式（12)中，Qu为负旋转备用风险指数，q{}表示大括号{}中事件发生的概率，大括 号U中事件为风电的实际出力大于平均出力导致的负旋转备用不足的事件，Pw为风机的输 出功率，添为风电平均出力，R Gp为发电系统可以提供的正旋转备用，RLp为发电系统中常规 机组的正旋转备用需求。
[0039]优选地，所述步骤3)的详细步骤包括：
[0040] 3.1)将发电系统中各常规机组所有可能出力的集合作为搜索空间，在搜索空间中 生成位置随机分布的多个指定维数的决策变量，将决策变量作为帝国竞争算法的国家 country = [XI，X2，…，xg]，其中g表示指定的维数，得到在搜索空间中的多个国家，每一个国 家对应一个目标函数值/(if h
[0041 ] 3.2)按照式（13)计算各个国家的势力，按照各个国家的势力大小将国家分为N。个 帝国主义国家，其余国家为殖民地国家；
[0042]
( 13)
[0043]式（13)中，Cn表示第η个国家的势力，fn表示第η个国家的目标函数值，fi表示第i个 国家的目标函数值，m表示国家数量；
[0044] 3.3)依照帝国主义国家间的势力大小来接管殖民地国家，使得势力越大的帝国主 义国家殖民地国家越多，一个帝国主义国家及其所接管的殖民地国家组成一个帝国；
[0045] 3.4)将殖民地国家代表的搜索空间中的位置向所属帝国的位置靠近并随机移动 一定的距离，判断是否存在帝国主义国家势力小于殖民地国家，如果存在则将殖民地国家、 所属帝国的帝国主义国家之间互换位置；
[0046] 3.5)根据式（14)计算各个帝国的总势力，在总势力最弱的帝国中挑选势力最弱的 殖民地国家按照式（15)的概率分配给其他Nc-1个有可能占有它的帝国，如果存在殖民地国 家数量为〇的帝国主义国家则消除该帝国主义国家的位置；
[0047] ..........
...... (14) Λ'。v J
[0048] 式（14)中，CAr表示帝国r的总势力，fr表示帝国r的目标函数值，〇为权重参数，又 a 表示帝国:r占有殖民地的目标函数平均值，N。表示帝国主义国家数量；
[0049]
(15)
[0050]式（15)中，Pj表示将总势力最弱的帝国中挑选势力最弱的殖民地国家分配给帝国 j的概率，CAj表示帝国j的总势力，Nc-1为可能占有总势力最弱的帝国中挑选势力最弱的殖 民地国家的帝国数量；
[0051] 3.6)判断是否只剩下唯一的帝国主义国家，如果不是则跳转执行步骤3.4)，否则， 将剩下唯一的帝国主义国家作为最优解输出。
[0052] 本发明计及风电出力不确定性的抽水蓄能-风电联合调度方法具有下述优点：
[0053] 1、本发明针对含抽蓄电站、常规机组及风电场的发电系统，将抽蓄电站作为平抑 风电出力波动性的技术手段控制抽蓄电站抽水工况与发电工况间的转换与风电场联合运 行，制定抽水蓄能-风电联合运行计划，能够有效平抑风电的波动性，降低系统的旋转备用 需求，显著提高了系统的经济效益。
[0054] 2、本发明建立基于常规机组的燃料费用和启停成本的目标函数，确定发电系统的 约束条件，精细化地考虑了风电出力不确定性的影响，实现了含风电电力系统运行经济性 与可靠性的统筹协调。
[0055] 3、本发明在目标函数的基础上，根据抽水蓄能-风电联合运行计划生成发电系统 中各常规机组所有可能出力的集合，将发电系统中各常规机组所有可能出力的集合作为搜 索空间，在搜索空间中生成位置随机分布的多个指定维数的决策变量，将决策变量作为帝 国竞争算法的国家、根据目标函数计算国家的势力、采用帝国竞争算法求解得到最优的一 个决策变量作为结果来控制发电系统中各常规机组的出力，利用帝国竞争算法解决了抽水 蓄能-风电联合运行模型的求解难题，从而能够快速有效地解决各台发电常规机组之间的 负荷经济分配问题。
【附图说明】
[0056] 图1为本发明实施例方法的整体流程图。
[0057] 图2为本发明实施例中的抽水蓄能-风电联合运行计划示意图。
[0058] 图3为本发明实施例中所使用帝国竞争算法流程图。
【具体实施方式】
[0059] 下文将以包含一个风电场、10台常规机组以及一个抽水蓄能电站的发电系统为 例，对本实施例计及风电出力不确定性的抽水蓄能-风电联合调度方法进行进一步说明如 下：
[0060] 如图1所示，本实施例计及风电出力不确定性的抽水蓄能-风电联合调度方法的步 骤包括：
[0061] 1)针对含抽蓄电站、常规机组及风电场的发电系统，将抽蓄电站作为平抑风电出 力波动性的技术手段控制抽蓄电站抽水工况与发电工况间的转换与风电场联合运行，制定 抽水蓄能-风电联合运行计划，建立基于常规机组的燃料费用和启停成本的目标函数，确定 发电系统的约束条件；
[0062] 2)根据抽水蓄能-风电联合运行计划生成发电系统中各常规机组所有可能出力的 集合；
[0063] 3)将发电系统中各常规机组所有可能出力的集合作为搜索空间，在搜索空间中生 成位置随机分布的多个指定维数的决策变量，将决策变量作为帝国竞争算法的国家、根据 目标函数计算国家的势力、采用帝国竞争算法求解得到最优的一个决策变量作为结果来控 制发电系统中各常规机组的出力。
[0064] 本实施例中将抽蓄电站作为平抑风电出力波动性的技术手段控制抽蓄电站抽水 工况与发电工况间的转换与风电场联合运行时，具体以风电场小时级的出力计划为基础， 并考虑抽水蓄能电站出力特性及运行成本等因素，将平抑风电出力波动性作为目标，制定 出抽水蓄能-风电联合运行计划，得到的抽水蓄能-风电联合运行计划如图2所示。参见图2 可知，风电出力较高（高于预设阈值)时，抽水蓄能电站工作在抽水工况;风电出力较低(低 于预设阈值）时，抽水蓄能电站则转换为发电工况。通过抽水蓄能电站与风电场的联合运 行，可以有效地平抑风电出力的波动性。
[0065] 本实施例中，步骤1)中的目标函数如式(1)所示；
[0066]
1 )
[0067] 式（1)中，./_(<)表示目标函数，T为发电系统的调度周期，N为发电系统中发电常规 机组数，Cn表示常规机组η的运行成本系数，g表示常规机组η在t时刻的出力，u nt表示常规 机组η在t时刻的状态，iMt-υ表示常规机组η在t-Ι时刻的状态，unt值为0表示常规机组η处于 停机状态，u nt值为1表示常规机组η处于开机状态，Snt为常规机组η在t时刻的启动费用。
[0068] 本实施例中，步骤1)中确定发电系统的约束条件包括如式(2)所示不计网损的情 况下的系统功率平衡约束；
[0069]
(2)
[0070] 式(2)中，N为发电系统中发电常规机组数，冗表示常规机组η在t时刻的出力，unt表 示常规机组η在t时刻的状态，巧表示在t时段风电的等效出力，匕表示第t时段的负荷预测 值。
[0071] 本实施例中，步骤1)中确定发电系统的约束条件还包括如式(3)所示的常规机组 出力约束；
[0072] P^' < P' < ΡΓ' (3；
[0073] 式⑶中，g表示常规机组η在t时刻的出力，表示常规机组η的最小出力，表 示常规机组η的最大出力。
[0074] 本实施例中，步骤1)中确定发电系统的约束条件还包括如式(4)所示的上升爬坡 速率约束、如式(5)所示的下降爬坡速率约束；
[0075] (4)
[0076] (5)
[0077] 式(4)和式(5)中，#表示常规机组η在t时刻的出力，/f1表示常规机组η在t-Ι时刻 的出力，Unt表示常规机组η在t时刻的状态，u n( t-υ表示常规机组η在t-Ι时刻的状态，DRn表示 常规机组η的下降爬坡速率，URn表示常规机组η的上升爬坡速率，△ T表示单位调度时段。
[0078] 本实施例中，步骤1)中确定发电系统的约束条件还包括如式(6)所示的抽蓄电站 库容约束、如式(7)所示的抽蓄电站输出功率约束、如式(8)所示的抽蓄电站日抽发电量约 束；
[0079]
(6)
[0080] 式(6)中，W。为抽蓄电站上水厍初始水量，Wmax为抽蓄电站上水库最大水量，qs为平 均水量转换系数，T为发电系统的调度周期，桴为时段t抽蓄电站的发电功率，％为平均电量 转换系数，Ps为抽蓄电站的抽水功率，€为时段t抽蓄电站处于抽水工况的机组数，w min为抽 蓄电站上水库最小水量；
[0081] ( 7 )
[0082] A W ;屮，，衣不H、」拐;t畑畜电站处于发电工况的机组数，Pg,min表不抽蓄电站发电 机组的最低出力，g为时段t抽蓄电站的发电功率，Ph表示抽蓄电站发电机组的额定功率， 於为时段t柚蓄电姑朴干柚7tc Τι况的机组数；
[0083]
(8)
[0084] 式(8)中，巧为时段t抽蓄电站的发电功率，％为平均电量转换系数，ns为平均水量 转换系数，τ为调度周期T内某一时刻，T为发电系统的调度周期，Ps为抽蓄电站的抽水功率， 埒为时段t抽蓄电站处于抽水工况的机组数。
[0085] 本实施例中，步骤1)中确定发电系统的约束条件还包括如式(9)所示的正旋转备 用风险指数约束、如式（10)所示的负旋转备用风险指数约束，所述正旋转备用风险指数为 风电的实际出力少于风电平均出力导致的正旋转备用不足的概率，所述负旋转备用风险指 数为风电的实际出力大于平均出力导致的负旋转备用不足的概率；
[0086] (9)
[0087] Qu^ (10)
[0088]式(9)和式（10)中，Qd为正旋转备用风险指数，Qu为负旋转备用风险指数，β为旋转 备用风险阈值。旋转备用风险阈值β通常在〇~1 〇 %之间取值。
[0089] 本实施例中，正旋转备用风险指数的函数表达式如式(11)所示；
[0090] ed = q{Pw I {Pw < Pw) & (RGi, <P"- Pw + Rlp)} (11)
[0091] 式（11)中，Qd为正旋转备用风险指数，q{}表示大括号{}中事件发生的概率，大括 号U中事件为风电的实际出力少于风电平均出力导致的正旋转备用不足的事件，Pw为风机 的输出功率，馬为风电平均出力，R Cp为发电系统可以提供的正旋转备用，RLp为发电系统中 常规机组的正旋转备用需求。
[0092] 本实施例中，负旋转备用风险指数的函数表达式如式(12)所示；
[0093] = q{Pw I (P,r >Pfr) & (Ra, < P"- - P"- + RLt,)\ i 12)
[0094] 式（12)中，Qu为负旋转备用风险指数，q{}表示大括号{}中事件发生的概率，大括 号U中事件为风电的实际出力大于平均出力导致的负旋转备用不足的事件，Pw为风机的输 出功率，硌为风电平均出力，R Gp为发电系统可以提供的正旋转备用，RLp为发电系统中常规 机组的正旋转备用需求。
[0095] 帝国竞争算法是具有全局搜索能力的智能算法，其借鉴了人类在封建帝国时代， 各帝国之间相互竞争并侵占对方殖民地以发展壮大的过程。该算法将初始种群个体定义为 国家，按照势力大小，将国家为"帝国主义国家"和"殖民地国家"，势力作为衡量国家强大与 否的指标，与优化目标函数相关，种群通过模拟帝国之间的竞争、获取殖民地的过程以寻求 最优解。该算法主要分为初始化帝国、同化殖民地国家、帝国主义国家竞争、最弱帝国灭亡 四个部分。本实施例中步骤3)利用帝国竞争算法解决了抽水蓄能-风电联合运行模型的求 解难题，从而能够快速有效地解决各台发电机组之间的负荷经济分配问题。
[0096] 如图3所示，本实施例中步骤3)的详细步骤包括：
[0097] 3.1)将发电系统中各常规机组所有可能出力的集合作为搜索空间，在搜索空间中 生成位置随机分布的多个指定维数的决策变量，将决策变量作为帝国竞争算法的国家 country = [XI，X2，…，xg]，其中g表示指定的维数，得到在搜索空间中的多个国家，每一个国 家对应一个目标函数值/(<);
[0098] 3.2)按照式（13)计算各个国家的势力，按照各个国家的势力大小将国家分为N。个 帝国主义国家，其余国家为殖民地国家；
[0099]
.(.1.3:)
[0100] 式（13)中，Cn表示第η个国家的势力，fn表示第η个国家的目标函数值，fi表示第i个 国家的目标函数值，m表示国家数量；
[0101] 3.3)依照帝国主义国家间的势力大小来接管殖民地国家，使得势力越大的帝国主 义国家殖民地国家越多，一个帝国主义国家及其所接管的殖民地国家组成一个帝国；
[0102] 3.4)将殖民地国家代表的搜索空间中的位置向所属帝国的位置靠近并随机移动 一定的距离，判断是否存在帝国主义国家势力小于殖民地国家，如果存在则将殖民地国家、 所属帝国的帝国主义国家之间互换位置；
[0103] 3.5)根据式（14)计算各个帝国的总势力，在总势力最弱的帝国中挑选势力最弱的 殖民地国家按照式（15)的概率分配给其他Nc-1个有可能占有它的帝国，如果存在殖民地国 家数量为〇的帝国主义国家则消除该帝国主义国家的位置；
[0104]
( 14 '
[0105] 式（14)中，CAr表示帝国r的总势力，fr表示帝国r的目标函数值，σ为权重参数， 表示帝国:r占有殖民地的目标函数平均值，N。表示帝国主义国家数量；
[0106]
(.15)
[0107] 式（15)中，Pj表示将总势力最弱的帝国中挑选势力最弱的殖民地国家分配给帝国 j的概率，CAj表示帝国j的总势力，Nc-1为可能占有总势力最弱的帝国中挑选势力最弱的殖 民地国家的帝国数量；
[0108] 3.6)判断是否只剩下唯一的帝国主义国家，如果不是则跳转执行步骤3.4)，否则， 将剩下唯一的帝国主义国家作为最优解输出。
[0109] 通过上述步骤，完成计及风电出力不确定性的抽水蓄能-风电联合优化调度。
[0110] 本实施例中计算过程中，常规机组参数和负荷数据以及风电场出力数据来源于科 研文献。使用了 Matlab进行了编程实现，取发电成本最小的那次作为最终计算结果。为了对 比分析本实施例计及风电出力不确定性的抽水蓄能-风电联合调度方法的有效性与正确 性，制定出抽水蓄能-风电联合运行计划时确立以下三种运行方式：
[0111] 运行方式1:不计抽水蓄能电站调度，风电场独立运行，令旋转备用风险阈值β为0， 风电的正负旋转备用需求为风电出力的最大波动幅值，此时发电系统的调度模型实质上为 传统风电确定性调度模型。运行方式2:计入抽水蓄能电站调度，风电场与抽水蓄能电站联 合运行，且令旋转备用风险阈值β为〇,此时发电系统的调度模型实质上为计及抽水蓄能-风 电联合运行的风电确定性调度模型。运行方式3:计入抽水蓄能电站调度，风电场与抽水蓄 能电站联合运行，且令旋转备用风险阈值β为〇.〇1，此时发电系统的调度模型实质上为计及 风电出力不确定性的抽水蓄能-风电联合调度模型。本实施例中，按运行方式3进行仿真计 算得出的24小时优化计算结果如表1所示，三种运行方式下的系统运行成本如表2所示。
[0112] 表1:运行方式3下10台常规机组的优化计算结果。
[0113]
[0114] 由表1可知，本申请调度方法的仿真结果为24个时刻的发电费用、启停费用，以及 各个时刻常规机组的启停状态与出力安排，为一次帝国竞争算法的输出结果。
[0115] 表2:三种运行方式下的系统运行成本比较表。
[0116]
[0117]由表2可知，运行方式3下的运行成本均低于运行方式1与运行方式2，验证了本申 请调度方法的正确性与有效性。
[0118]以上所述仅是本发明的优选实施方式，本发明的保护范围并不仅局限于上述实施 例，凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出，对于本技术领域 的普通技术人员来说，在不脱离本发明原理前提下的若干改进和润饰，这些改进和润饰也 应视为本发明的保护范围。
【主权项】
1. 一种计及风电出力不确定性的抽水蓄能-风电联合调度方法，其特征在于步骤包括： 1) 针对含抽蓄电站、常规机组及风电场的发电系统，将抽蓄电站作为平抑风电出力波 动性的技术手段控制抽蓄电站抽水工况与发电工况间的转换与风电场联合运行，制定抽水 蓄能-风电联合运行计划，建立基于常规机组的燃料费用和启停成本的目标函数，确定发电 系统的约束条件； 2) 根据抽水蓄能-风电联合运行计划生成发电系统中各常规机组所有可能出力的集 合； 3) 将发电系统中各常规机组所有可能出力的集合作为搜索空间，在搜索空间中生成位 置随机分布的多个指定维数的决策变量，将决策变量作为帝国竞争算法的国家、根据目标 函数计算国家的势力、采用帝国竞争算法求解得到最优的一个决策变量作为结果来控制发 电系统中各常规机组的出力。2. 根据权利要求1所述的计及风电出力不确定性的抽水蓄能-风电联合调度方法，其特 征在于，所述步骤1)中的目标函数如式(1)所示；式α)中，/(<)表示目标函数，T为发电系统的调度周期，N为发电系统中发电常规机组 数，Cn表示常规机组η的运行成本系数，<表示常规机组η在t时刻的出力，unt表示常规机组η 在t时刻的状态，表示常规机组η在t-1时刻的状态，unt值为O表示常规机组η处于停机 状态，u nt值为1表示常规机组η处于开机状态，Snt为常规机组η在t时刻的启动费用。3. 根据权利要求1所述的计及风电出力不确定性的抽水蓄能-风电联合调度方法，其特 征在于，所述步骤1)中确定发电系统的约束条件包括如式(2)所示不计网损的情况下的系 统功率平衡约束；式(2)中，N为发电系统中常规机组数，表示常规机组η在t时刻的出力，unt表示常规机 组η在t时刻的状态，表示在t时段风电的等效出力，/?表示第t时段的负荷预测值。4. 根据权利要求3所述的计及风电出力不确定性的抽水蓄能-风电联合调度方法，其特 征在于，所述步骤1)中确定发电系统的约束条件还包括如式(3)所示的常规机组出力约束；式(3)中，<表示常规机组η在t时刻的出力，Zfm表示常规机组η的最小出力，Cax表示常 规机组η的最大出力。5. 根据权利要求4所述的计及风电出力不确定性的抽水蓄能-风电联合调度方法，其特 征在于，所述步骤1)中确定发电系统的约束条件还包括如式(4)所示的上升爬坡速率约束、 如式(5)所示的下降爬坡速率约束；式(4)和式(5)中，<表示常规机组η在t时刻的出力，if1表示常规机组η在t-Ι时刻的出 力，Unt表示常规机组η在t时刻的状态，iMt-υ表示常规机组η在t-ι时刻的状态，DRn表示常规 机组η的下降爬坡速率，URn表示常规机组η的上升爬坡速率，ΛΤ表示单位调度时段。6. 根据权利要求5所述的计及风电出力不确定性的抽水蓄能-风电联合调度方法，其特 征在于，所述步骤1)中确定发电系统的约束条件还包括如式(6)所示的抽蓄电站库容约束、 如式(7)所示的抽蓄电站输出功率约束、如式(8)所示的抽蓄电站日抽发电量约束；式(6)中，W。为抽蓄电站上水库初始水量，WmaxS抽蓄电站上水库最大水量，ns为平均水 量转换系数，T为发电系统的调度周期，为时段t抽蓄电站的发电功率，取为平均电量转换 系数，Ps为抽蓄电站的抽水功率，g为时段t抽蓄电站处于抽水工况的机组数，W min为抽蓄电 站上水库最小水量；式(7)中，苟表示时段t抽蓄电站处于发电工况的机组数，PC,min表示抽蓄电站发电机组 的最低出力，踔为时段t抽蓄电站的发电功率，Ph表示抽蓄电站发电机组的额定功率，硌为 时段t抽蓄电站处于抽水工况的机组数；式(8)中，罐为时段t抽蓄电站的发电功率，％为平均电量转换系数，1为平均水量转换 系数，τ为调度周期T内某一时刻，T为发电系统的调度周期，Ps为抽蓄电站的抽水功率，句为 时段t抽蓄电站处于抽水工况的机组数。7. 根据权利要求6所述的计及风电出力不确定性的抽水蓄能-风电联合调度方法，其特 征在于，所述步骤1)中确定发电系统的约束条件还包括如式(9)所示的正旋转备用风险指 数约束、如式（10)所示的负旋转备用风险指数约束，所述正旋转备用风险指数为风电的实 际出力少于风电平均出力导致的正旋转备用不足的概率，所述负旋转备用风险指数为风电 的实际出力大于平均出力导致的负旋转备用不足的概率； Qd细 (9) (10) 式(9)和式（10)中，Qd为正旋转备用风险指数，Qu为负旋转备用风险指数，β为旋转备用 风险阈值。8. 根据权利要求7所述的计及风电出力不确定性的抽水蓄能-风电联合调度方法，其特 征在于，所述正旋转备用风险指数的函数表达式如式(11)所示；式(11)中，Qd为正旋转备用风险指数，q{}表示大括号{}中事件发生的概率，大括号{}中 事件为风电的实际出力少于风电平均出力导致的正旋转备用不足的事件，Pw为风机的输出 功率，晷为风电平均出力，Rcp为发电系统可以提供的正旋转备用，R lp为发电系统中常规机 组的正旋转备用需求。9. 根据权利要求7所述的计及风电出力不确定性的抽水蓄能-风电联合调度方法，其特 征在于，所述负旋转备用风险指数的函数表达式如式(12)所示； Qu=q{Pw\{P}r>P,r)&{Rai<Pir-Pw + K,)\ C12) 式(12)中，Qu为负旋转备用风险指数，q{}表示大括号{}中事件发生的概率，大括号{}中 事件为风电的实际出力大于平均出力导致的负旋转备用不足的事件，Pw为风机的输出功 率，愚为风电平均出力，RgpS发电系统可以提供的正旋转备用，R lpS发电系统中常规机组 的正旋转备用需求。10. 根据权利要求1~9中任意一项所述的计及风电出力不确定性的抽水蓄能-风电联 合调度方法，其特征在于，所述步骤3)的详细步骤包括： 3.1) 将发电系统中各常规机组所有可能出力的集合作为搜索空间，在搜索空间中生成 位置随机分布的多个指定维数的决策变量，将决策变量作为帝国竞争算法的国家country =[Xl，X2，…，Xg ]，其中g表示指定的维数，得到在搜索空间中的多个国家，每一个国家对应 一个目标函数值/(硭）; 3.2) 按照式（13)计算各个国家的势力，按照各个国家的势力大小将国家分为N。个帝国 主义国家，其余国家为殖民地国家；式（13)中，Cn表示第η个国家的势力，fn表示第η个国家的目标函数值，心表示第i个国家 的目标函数值，m表示国家数量； 3.3) 依照帝国主义国家间的势力大小来接管殖民地国家，使得势力越大的帝国主义国 家殖民地国家越多，一个帝国主义国家及其所接管的殖民地国家组成一个帝国； 3.4) 将殖民地国家代表的搜索空间中的位置向所属帝国的位置靠近并随机移动一定 的距离，判断是否存在帝国主义国家势力小于殖民地国家，如果存在则将殖民地国家、所属 帝国的帝国主义国家之间互换位置； 3.5) 根据式（14)计算各个帝国的总势力，在总势力最弱的帝国中挑选势力最弱的殖民 地国家按照式（15)的概率分配给其他Nc-I个有可能占有它的帝国，如果存在殖民地国家数 量为〇的帝国主义国家则消除该帝国主义国家的位置；式(14)中，Cai·表示帝国r的总势力，f r表示帝国r的目标函数值，σ为权重参数，Is，表示帝 国r占有殖民地的目标函数平均值，Ν。表示帝国主义国家数量；式(15)中，Pj表示将总势力最弱的帝国中挑选势力最弱的殖民地国家分配给帝国j的概 率，CAj表示帝国j的总势力，Nc-I为可能占有总势力最弱的帝国中挑选势力最弱的殖民地国 家的帝国数量； 3.6) 判断是否只剩下唯一的帝国主义国家，如果不是则跳转执行步骤3.4)，否则，将剩 下唯一的帝国主义国家作为最优解输出。
【文档编号】G06Q50/06GK105932722SQ201610340010
【公开日】2016年9月7日
【申请日】2016年5月20日
【发明人】陈道君, 张斌, 张磊, 呙虎, 沈阳武, 左剑, 崔挺
【申请人】国家电网公司, 国网湖南省电力公司, 国网湖南省电力公司电力科学研究院