一种面向风电消纳的两阶段源荷储优化调度方法与流程

本发明涉及电热联合系统经济调度技术领域，具体涉及一种面向风电消纳的两阶段源荷储优化调度方法。

背景技术：

随着我国风电装机容量的快速增长，风电占系统总发电比例逐渐增加，但目前风电弃风问题突出，在冬季供暖季，大量热电联产机组(chp)工作在“以热定电”的模式下以满足热负荷需求，其调节能力受到限制，制约了风电的消纳，同时其出力的随机性、不确定性和间歇性导致系统的功率波动进一步扩大，常规电源的调节能力难以有效应对，给电网传统调度运行方式带来了新的挑战，风电预测的不确定性对调度计划的制定与实施产生了很大的影响，且预测误差与预测时间跨度呈正相关，根据日前预测数据制定调度计划与实际参考计划偏差较大，无法满足系统的调度要求。

在实际中需要常规机组提供备用和频繁调峰，系统整体的运行成本上升，单独调度发电侧资源，会受到设备投资和爬坡等因素的限制，导致调度灵活性不足，不能充分消纳风电，但从负荷侧调控，虽然实现了“削峰填谷”，但受用户满意度和高额的调度成本，其调节范围有限，仅调度储能设备，受到设备容量、充放电成本，其消纳能力有限，源荷储协调调度通过对系统中电源侧资源、储能设备和需求响应等可调度资源进行综合调度，提高风电消纳量，减小系统运行成本。

因此需要制定优化调度策略，充分调用系统源荷储资源，通过制定电价、热价提升负荷响应的积极性，同时通过储能将能量进行时间上的平移，起到“削峰填谷”的作用，提升系统运行调度的灵活性，平滑风电波动。随预测时间的临近，风电的预测的准确性越高，用精确度高的预测值对精确度低的进行修正，提高风电预测的精度，从而制定更为合理的调度策略，减小弃风量。

技术实现要素：

针对现有的技术的不足，本发明提供一种面向风电消纳的两阶段源荷储优化调度方法，该方法通过建立两阶段调度模型，综合调用源荷储资源，进行削峰填谷，平滑风电波动，促进风电消纳，保证电热联合系统运行的经济性，并采用改进粒子群算法对两阶段调度模型进行求解。

一种面向风电消纳的两阶段源荷储优化调度方法，包括以下步骤：

步骤1：获取电热联合系统中日前的风电出力预测值、电、热负荷的预测功率值，所述电热联合系统包括风电机组、常规机组、热电联产机组三大机组；

步骤2：根据电负荷的预测功率制定峰、谷、平的分时电价；

步骤3：根据热负荷的预测功率制定峰、谷、平的分时热价；

步骤4：根据电、热负荷对分时价格的响应，建立电、热负荷需求响应模型；

步骤5：根据日前风电出力预测值和电负荷的需求响应模型、热负荷的需求响应模型，以电热联合系统的总运行费用—发电费用和分时价格用户的费用最小为目标，建立日前调度模型；

步骤6：根据实时采集的环境温度和pmv指标，建立激励型负荷、采暖负荷的日内调度响应模型；

步骤7：以弃风量最小为目标，建立电热联合系统的日内调度模型，更新风电预测功率，调用电储能、激励型负荷和采暖负荷，对日前调度进行微调；

步骤8：分别采用改进的粒子群算法求解两阶段源荷储优化调度模型，得到在风电消纳情况下n个常规机组的出力值、m个热电联产机组的出力值、蓄热式电锅炉出力值、调用的激励型负荷值pl1,t、采暖负荷的热功率值qt、电储能出力值以及风电机组出力值，所述两阶段源荷储优化调度模型包括日前调度模型和日内调度模型。

所述的步骤2具体步骤如下：

2.1)利用公式(1)～公式(2)对一天内的电负荷进行峰、谷、平分时段的划分：

式中，pav表示电负荷的平均功率，pl,t⁰表示t时段电负荷的预测功率值，tg,p表示电负荷在t时段被划分为谷时段，tp,p表示电负荷在t时段被划分为平时段，tf,p表示电负荷在t时段被划分为峰时段，αp、βp分别表示电负荷在划分时段时的比例系数；

2.2)利用公式(3)～公式(5)，根据电负荷峰谷电量的比值，制定电负荷峰、谷、平的分时电价；

其中，当kp取值范围为1<kp<k1,p时，电负荷峰、谷、平的分时电价ct,p表示为：

其中，当kp取值范围为k1,p≤kp时，电负荷峰、谷、平的分时电价c't,p表示为：

式中，kp表示电负荷峰谷功率比，k1,p表示预设电负荷峰谷功率值，其中kp根据一天中的电负荷功率最大值pmax和电负荷功率最小值pmin确定，pmax表示一天中的电负荷功率最大值，pmin表示一天中的电负荷功率最小值，kp越大表示电负荷峰谷电量差越大，c0,p表示未实行分时电价时的统一电价，m1,p、m2,p、m3,p、m4,p分别表示电负荷不同分时段的比例系数，且满足m3,p<m1,p、m2,p<m4,p、m1,p<m2,p、m3,p<m4,p。

所述的步骤3具体表述为：

3.1)利用公式(6)～公式(7)对一天内的热负荷进行峰、谷、平分时段的划分：

式中，qav表示热负荷的平均功率，ql,t⁰表示t时段热负荷的预测功率值，tg,q表示热负荷在t时段被划分为谷时段，tp,q表示热负荷在t时段被划分为平时段，tf,q表示热负荷在t时段被划分为峰时段，αq、βq分别表示热负荷在划分时段时的比例系数；

3.2)利用公式(8)～公式(10)，根据热负荷峰谷电量的比值，制定热负荷峰、谷、平的分时热价；

其中，当kq取值范围为1<kq<k1,q时，热负荷峰、谷、平的分时热价ct,q表示为：

其中，当kq取值范围为k1,q≤kq时，热负荷峰、谷、平的分时热价c't,q表示为：

式中，kq表示热负荷峰谷功率比，k1,q表示预设热负荷峰谷功率值，其中kq根据一天中的热负荷功率最大值qmax和热负荷功率最小值qmin确定，qmax表示一天中的热负荷功率最大值，qmin表示一天中的热负荷功率最小值，kq越大表示热负荷峰谷电量差越大，c0,q表示未实行分时热价时的统一热价，m1,q、m2,q、m3,q、m4,q分别表示热负荷不同分时段的比例系数，且满足m3,q<m1,q、m2,q<m4,q、m1,q<m2,q、m3,q<m4,q，分时电/热价格引导电/热负荷削峰填谷，电/热峰谷差越大，采用的电/热分时价格差越大，越能调动电/热负荷响应的积极性。

所述的步骤4中根据电负荷对分时价格的响应，建立的电负荷需求响应模型为：

式中，pl,f⁰、pl,g⁰、pl,p⁰分别表示未实施分时电价时电负荷峰、谷、平的分时用电负荷，pl,f、pl,g、pl,p分别表示实施分时电价后电负荷峰、谷、平的分时用电负荷，e表示电负荷的价格弹性需求矩阵，cf,p⁰、cg,p⁰、cp,p⁰分别表示实施分时电价前电负荷峰、谷、平的用电价格，△cf,p、△cg,p、△cp,p分别表示实施分时电价后的电负荷峰、谷、平的电价变化；

所述的步骤4中根据热负荷对分时价格的响应，建立的热负荷需求响应模型为：

式中，ql,f⁰、ql,g⁰、ql,p⁰分别表示未实施分时热价时热负荷峰、谷、平的分时用热负荷，ql,f、ql,g、ql,p分别表示实施分时热价后热负荷峰、谷、平的分时用热负荷，e'表示热负荷的价格弹性需求矩阵，cf,q⁰、cg,q⁰、cp,q⁰分别表示实施分时热价前热负荷峰、谷、平的用热价格，△cf,q、△cg,q、△cp,q分别表示实施分时热价后的热负荷峰、谷、平的热价变化。

所述的步骤5具体步骤如下：

5.1)根据常规机组的经济成本系数建立常规机组一天内的成本函数：

式中，pi,t表示第i个常规机组在t时刻的出力值，ai、bi、ci分别表示第i个常规机组的经济成本系数，n表示常规机组的总个数；

5.2)根据热电联产机组的经济成本系数建立热电联产机组一天内的成本函数：

式中，pj,t表示第j个热电联产机组在t时刻的输出电功率，hj,t表示第j个热电联产机组在t时刻的输出热功率，aj、bj、cj、dj、ej、fj分别表示第j个热电联产机组的经济成本系数，m表示热电联产机组的总个数；

5.3)存在弃风时蓄热式电锅炉开启，利用弃风供热，利用的是弃风电量，不存在成本；

5.4)建立电负荷一天内的分时电价成本fpl为：

式中，pl,t表示分时电价负荷；

5.5)建立热负荷一天内的分时热价成本fql为：

式中，pl,t表示分时热价负荷；

5.6)以电热联合系统的总运行费用—发电费用和分时价格用户的费用最小为目标，建立日前调度模型为：

所述的步骤6具体表述为：

6.1)由于工业用户的可调度容量大，可以应对风电的波动，促进风电的消纳，对激励型负荷增加或削减的负荷，由电网公司予以补偿，则激励型负荷的日内调度响应约束及调用成本为：

plmin≤pl1,t≤plmax(18)

式中：pl1,t为t时刻调用的激励型负荷值，plmin、plmax分别表示可调用的激励型负荷的下限值、上限值，cdr表示激励型负荷的调用成本，ρ表示激励型负荷的电网公司补偿成本。

6.2)根据pmv指标以及室外建筑物的一阶热力学模型，建立采暖负荷的日内调度响应模型为：

qt,min≤qt≤qt,max(20)

其中，

式中，qt,min表示采暖负荷的最小值，qt,max表示采暖负荷的最大值，it^pmv表示pmv指标，θtⁱⁿ表示t时刻的室内温度，θsk表示人体舒适状态的皮肤温度值，mt表示人体代谢率，r表示服装热阻，qt表示采暖负荷的热功率，θt+1ⁱⁿ表示t+1时刻的室内温度值，r表示建筑物的等效热阻，θt^out表示t时刻的室外温度值，h为常数，且满足其中对于特定的建筑物在△t时刻cair为常数。

所述的步骤7具体表述为：

7.1)以弃风电量最小为目标，建立电热联合系统的日内调度模型为：

min(pw,rel,t-pw,t)²(22)

式中，pw,rel,t表示t时段消纳的风电量，pw,t表示t时段风电出力值；

7.2)根据电热联合系统的日内调度模型，对电热联合系统的日前调度进行微调时，由于日内的风电出力预测值比日前的风电出力预测值准确，所以通过电热联合系统的日内调度模型对电热联合系统的日前调度模型进行修正，具体表述为：若产生弃风现象，首先调用激励型负荷的需求响应和采暖负荷的需求响应消纳风电负荷，若激励型负荷和采暖负荷的需求响应仍不能完全消纳风电负荷，调用电储能负荷的需求响应来消纳风电负荷。

所述的步骤8中的采用改进的粒子群算法求解日前调度模型的具体步骤为：

8.1.1)设置求解日前调度模型时的粒子的最大速度值vu,max、最大迭代次数kmax、种群规模f1，以及初始化粒子的位置xu⁰和速度vu⁰，为了缩短收敛时间，各粒子的初始位置选择时应满足公式(23)给出的取值范围，

xu,min≤xu⁰≤xu,max(u＝1,2...u)(23)

式中，xu⁰表示粒子u的初始位置，u表示求解日前调度模型时的粒子总个数，xu,min表示粒子u初始位置的最小值，xu,max表示粒子u初始位置的最大值，其中xu,min和xu,max根据电热联合系统中三大机组出力的上限值和下限值、蓄热式电锅炉出力的上限值和下限值确定；

8.1.2)计算每个粒子的适应值，并与粒子最优值和全局最优值作比较，更新粒子u的个体极值pbestu和粒子u的全局极值gbestu；

8.1.3)为了避免粒子群算法陷入全局最优或局部最优，利用公式(24)对ω^g进行动态调整，以保证算法的收敛性，

式中，ω^g表示第g次迭代的惯性权重，ωmax表示最大惯性权重，ωmin表示最小惯性权重，kmax表示求解日前调度模型时的迭代总次数；

8.1.4)利用公式(25)将第g次迭代的惯性权重ω^g代入速度更新公式中，更新微粒速度；

vud^g＝ω^gvud^g-1+c1rand(0,1)(pbestu-xu^g)+c2rand(0,1)(gbestu-xu^g)(25)

式中，vud^g表示粒子u第g次迭代的速度，vud^g-1表示粒子u第g-1次迭代的速度，c1、c2分别表示学习因子，pbestu表示粒子u的个体极值，gbestu表示粒子u的全局极值，xu^g表示粒子u第g次迭代的位置；

8.1.5)利用公式(26)迭代更新粒子的位置，

xu^g＝xu^g-1+vud^g(26)

式中，xu^g-1表示粒子u第g-1次迭代的位置，xu^g表示粒子u第g次迭代的位置，vud^g表示粒子u第g次迭代的速度；

8.1.6)重复执行步骤8.1.2)～步骤8.1.5)，直至达到求解日前调度模型时的迭代总次数kmax或全局最优位置满足最小界限，算法执行结束，输出在风电消纳情况下的n个常规机组的出力值、m个热电联产机组的出力值，蓄热式电锅炉的出力值以及风电机组的出力值；

所述的步骤8中的采用改进的粒子群算法求解日内调度模型的具体步骤如下：

8.2.1)设置求解日内调度模型时的粒子的最大速度值vv,max、最大迭代次数k'max、种群规模f2，以及初始化粒子的位置xv⁰和速度vv⁰，为了缩短收敛时间，各粒子的初始位置选择时应满足公式(27)给出的取值范围，

xv,min≤xv⁰≤xv,max(v＝1,2...v)(27)

式中，xv⁰表示粒子v的初始位置，v表示求解日内调度模型时的粒子总个数，xv,min表示粒子v初始位置的最小值，xv,max表示粒子v初始位置的最大值，其中xv,min和xv,max根据电热联合系统中风电机组出力的上限值和下限值、可调用的激励型负荷的上限值plmax、下限值plmin，采暖负荷的最小值qt,min、最大值qt,max，电储能出力的上限值和下限值确定；

8.2.2)计算每个粒子的适应值，并与粒子最优值和全局最优值作比较，更新粒子v的个体极值pbestv和粒子v的全局极值gbestv；

8.2.3)为了避免粒子群算法陷入全局最优或局部最优，利用公式(28)对ω^h进行动态调整，以保证算法的收敛性，

式中，ω^h表示第h次迭代的惯性权重，ωmax表示最大惯性权重，ωmin表示最小惯性权重，k'max表示求解日内调度模型时的迭代总次数；

8.2.4)利用公式(29)将第h次迭代的惯性权重ω^h代入速度更新公式中，更新微粒速度；

vvd^h＝ω^hvvd^h-1+c3rand(0,1)(pbestv-xv^h)+c4rand(0,1)(gbestv-xv^h)(29)

式中，vvd^h表示粒子v第h次迭代的速度，vvd^h表示粒子v第h-1次迭代的速度，c3、c4分别表示学习因子，pbestv表示粒子v的个体极值，gbestv表示粒子v的全局极值，xv^h表示粒子v第h次迭代的位置；

8.2.5)利用公式(30)迭代更新粒子的位置，

xv^h＝xv^h-1+vvd^h(30)

式中，xv^h-1表示粒子v第h-1次迭代的位置，xv^h表示粒子v第h次迭代的位置，vvd^h表示粒子v第h次迭代的速度；

8.2.6)重复执行步骤8.2.2)～步骤8.2.5)，直至达到求解日内调度模型时的迭代总次数k'max或全局最优位置满足最小界限，算法执行结束，输出在风电消纳情况下调用的激励型负荷值pl1,t、采暖负荷的热功率值qt、电储能出力值以及风电机组的出力值。

本发明的有益效果：

1)本发明考虑了风电的不确定性，随着时间的临近，风电的预测精度逐渐增大，通过时前风电预测对日前预测进行修正，提高了风电的预测精度；

2)本发明通过日前负荷预测，得出峰谷电量比，根据峰谷电量比制定合适的分时电价，峰谷比大时，分时电价的价格差异性大，峰谷比小时，分时电价的价格差异性小，提高负荷参与需求响应的积极性，实现负荷的削峰填谷；

3)本发明在日内调度时，调用激励型负荷和采暖负荷，激励型负荷可根据电网的调度要求适时地调整用电，可调度量大；采暖负荷的可调范围为一个区间，可以根据风电出力的情况实时调整采暖负荷的用热量，提高风电的利用率；

4)本发明将需求响应模型按照不同的特性分类，即分为电、热负荷分时响应模型、激励型负荷响应模型和采暖负荷响应模型，在日前和日内调度两个不同的时间尺度下选取相应的需求响应进行调度；

附图说明

图1为本发明中的一种面向风电消纳的两阶段源荷储优化调度方法流程图。

图2为本发明中的分时电价制定流程图。

图3为本发明中的改进粒子群算法求解流程图。

图4为本发明中的建筑物一节热力学模型图。

具体实施方式

下面结合附图和具体实施实例对发明做进一步说明。

如图1所示，一种面向风电消纳的两阶段源荷储优化调度方法，包括以下步骤：

步骤1：获取电热联合系统中日前的风电出力预测值、电热负荷预测功率值，所述电热联合系统包括风电机组、常规机组、热电联产机组三大机组，本发明中所述常规机组为纯发电的火电机组；

步骤2：根据电负荷的预测功率制定峰、谷、平的分时电价，具体步骤如下：

2.1)利用公式(1)～公式(2)对一天内的电负荷进行峰、谷、平分时段的划分：

式中，pav表示电负荷的平均功率，pl,t⁰表示t时段电负荷的预测功率值，tg,p表示电负荷在t时段被划分为谷时段，tp,p表示电负荷在t时段被划分为平时段，tf,p表示电负荷在t时段被划分为峰时段，αp、βp分别表示电负荷在划分时段时的比例系数；

2.2)利用公式(3)～公式(5)，根据电负荷峰谷电量的比值，制定电负荷峰、谷、平的分时电价，其流程图如图2所示；

其中，当kp取值范围为1<kp<k1,p时，电负荷峰、谷、平的分时电价ct,p表示为：

其中，当kp取值范围为k1,p≤kp时，电负荷峰、谷、平的分时电价c't,p表示为：

式中，kp表示电负荷峰谷功率比，k1,p表示预设电负荷峰谷功率值，其中kp根据一天中的电负荷功率最大值pmax和电负荷功率最小值pmin确定，pmax表示一天中的电负荷功率最大值，pmin表示一天中的电负荷功率最小值，kp越大表示电负荷峰谷电量差越大，c0,p表示未实行分时电价时的统一电价，m1,p、m2,p、m3,p、m4,p分别表示电负荷不同分时段的比例系数，且满足m3,p<m1,p、m2,p<m4,p、m1,p<m2,p、m3,p<m4,p；

步骤3：根据热负荷的预测功率制定峰、谷、平的分时热价，具体表述为：

3.1)利用公式(6)～公式(7)对一天内的热负荷进行峰、谷、平分时段的划分：

式中，qav表示热负荷的平均功率，ql,t⁰表示t时段热负荷的预测功率值，tg,q表示热负荷在t时段被划分为谷时段，tp,q表示热负荷在t时段被划分为平时段，tf,q表示热负荷在t时段被划分为峰时段，αq、βq分别表示热负荷在划分时段时的比例系数；

3.2)利用公式(8)～公式(10)，根据热负荷峰谷电量的比值，制定热负荷峰、谷、平的分时热价；

其中，当kq取值范围为1<kq<k1,q时，热负荷峰、谷、平的分时热价ct,q表示为：

其中，当kq取值范围为k1,q≤kq时，热负荷峰、谷、平的分时热价c't,q表示为：

式中，kq表示热负荷峰谷功率比，k1,q表示预设热负荷峰谷功率值，其中kq根据一天中的热负荷功率最大值qmax和热负荷功率最小值qmin确定，qmax表示一天中的热负荷功率最大值，qmin表示一天中的热负荷功率最小值，kq越大表示热负荷峰谷电量差越大，c0,q表示未实行分时热价时的统一热价，m1,q、m2,q、m3,q、m4,q分别表示热负荷不同分时段的比例系数，且满足m3,q<m1,q、m2,q<m4,q、m1,q<m2,q、m3,q<m4,q，分时电/热价格引导电/热负荷削峰填谷，电/热峰谷差越大，采用的电/热分时价格差越大，越能调动电/热负荷响应的积极性；

步骤4：根据电、热负荷对分时价格的响应，建立电、热负荷需求响应模型；

其中电负荷的需求响应模型为：

式中，pl,f⁰、pl,g⁰、pl,p⁰分别表示未实施分时电价时电负荷峰、谷、平的分时用电负荷，pl,f、pl,g、pl,p分别表示实施分时电价后电负荷峰、谷、平的分时用电负荷，e表示电负荷的价格弹性需求矩阵，cf,p⁰、cg,p⁰、cp,p⁰分别表示实施分时电价前电负荷峰、谷、平的用电价格，△cf,p、△cg,p、△cp,p分别表示实施分时电价后的电负荷峰、谷、平的电价变化；

其中热负荷的需求响应模型为：

式中，ql,f⁰、ql,g⁰、ql,p⁰分别表示未实施分时热价时热负荷峰、谷、平的分时用热负荷，ql,f、ql,g、ql,p分别表示实施分时热价后热负荷峰、谷、平的分时用热负荷，e'表示热负荷的价格弹性需求矩阵，cf,q⁰、cg,q⁰、cp,q⁰分别表示实施分时热价前热负荷峰、谷、平的用热价格，△cf,q、△cg,q、△cp,q分别表示实施分时热价后的热负荷峰、谷、平的热价变化；

步骤5：根据日前风电出力预测值和电负荷的需求响应模型、热负荷的需求响应模型，以电热联合系统的总运行费用—发电费用和分时价格用户的费用最小为目标，建立日前调度模型，具体步骤如下：

5.1)根据常规机组的经济成本系数建立常规机组一天内的成本函数：

式中，pi,t表示第i个常规机组在t时刻的出力值，ai、bi、ci分别表示第i个常规机组的经济成本系数，n表示常规机组的总个数；

5.2)根据热电联产机组的经济成本系数建立热电联产机组一天内的成本函数：

式中，pj,t表示第j个热电联产机组在t时刻的输出电功率，hj,t表示第j个热电联产机组在t时刻的输出热功率，aj、bj、cj、dj、ej、fj分别表示第j个热电联产机组的经济成本系数，m表示热电联产机组的总个数；

5.3)存在弃风时蓄热式电锅炉开启，利用弃风供热，利用的是弃风电量，不存在成本；

5.4)建立电负荷一天内的分时电价成本为：

式中，pl,t表示分时电价负荷；

5.5)建立热负荷一天内的分时热价成本为：

式中，pl,t表示分时热价负荷；

5.6)以电热联合系统的总运行费用—发电费用和分时价格用户的费用最小为目标，建立日前调度模型为：

步骤6：根据实时采集的环境温度和pmv指标，建立激励型负荷、采暖负荷的日内调度响应模型，具体表述为：

6.1)由于工业用户的可调度容量大，可以应对风电的波动，促进风电的消纳，对激励型负荷增加或削减的负荷，由电网公司予以补偿，则激励型负荷的日内调度响应约束及调用成本为：

plmin≤pl1,t≤plmax(18)

式中：pl1,t为t时刻调用的激励型负荷值，plmin、plmax分别表示可调用的激励型负荷的下限值、上限值，cdr表示激励型负荷的调用成本，ρ表示激励型负荷的电网公司补偿成本；

6.2)根据pmv指标以及室外建筑物的一阶热力学模型，其中建筑物的一阶热力学模型如图4所示，建立采暖负荷的日内调度响应模型为：

qt,min≤qt≤qt,max(20)

其中，

式中，qt,min表示采暖负荷的最小值，qt,max表示采暖负荷的最大值，it^pmv表示pmv指标，θtⁱⁿ表示t时刻的室内温度，θsk表示人体舒适状态的皮肤温度值，θsk近似取为33.5℃，mt表示人体代谢率，mt取为60w/m²，r表示服装热阻，r取为0.11(m²·℃)/w，qt表示采暖负荷的热功率，θt+1ⁱⁿ表示t+1时刻的室内温度值，r表示建筑物的等效热阻，θt^out表示t时刻的室外温度值，h为常数，且满足其中对于特定的建筑物在△t时刻cair为常数；

步骤7：以弃风量最小为目标，建立电热联合系统的日内调度模型，更新风电预测功率，调用电储能、激励型负荷和采暖负荷，对日前调度进行微调，具体表述为：

7.1)以弃风电量最小为目标，建立电热联合系统的日内调度模型为：

min(pw,rel,t-pw,t)²(22)

式中，pw,rel,t表示t时段消纳的风电量，pw,t表示t时段风电出力值；

7.2)根据电热联合系统的日内调度模型，对电热联合系统的日前调度进行微调时，由于日内的风电出力预测值比日前的风电出力预测值准确，所以通过电热联合系统的日内调度模型对电热联合系统的日前调度模型进行修正，具体表述为：若产生弃风现象，首先调用激励型负荷的需求响应和采暖负荷的需求响应消纳风电负荷，若激励型负荷和采暖负荷的需求响应仍不能完全消纳风电负荷，调用电储能负荷的需求响应来消纳风电负荷；

仅调风电机组、chp机组、风电机组及分时负荷不能满足功率平衡条件，出现弃风时的判断条件如公式(31)所示，此时考虑调用激励型负荷和采暖负荷消纳风电：

pw,t+pj,t+pi,t＞pl,t+pdg,t(31)

式中：pw,t表示风电机组的出力值，pj,t表示表示第j个热电联产机组在t时刻的输出电功率，pi,t表示常规机组的电出力，pdg,t表示蓄热式电锅炉用电量。

调节激励型负荷仍产生弃风时的判断条件如公式(32)所述，此时考虑调用储电设备消纳风电：

pw,t+pj,t+pi,t＞pl,t+pl1,t+pdg,t(32)

式中：pl1,t表示t时刻调用的激励型负荷值；

步骤8：分别采用改进的粒子群算法求解两阶段源荷储优化调度模型，得到在风电消纳情况下n个常规机组的出力值、m个热电联产机组的出力值，蓄热式电锅炉出力值，调用的激励型负荷值pl1,t、采暖负荷的热功率值qt、电储能出力值以及风电机组出力值，所述两阶段源荷储优化调度模型包括日前调度模型和日内调度模型；

如图3所示，采用改进的粒子群算法求解日前调度模型的具体步骤为：

8.1.1)设置求解日前调度模型时的粒子的最大速度值vu,max、最大迭代次数kmax、种群规模f1，以及初始化粒子的位置xu⁰和速度vu⁰，为了缩短收敛时间，各粒子的初始位置选择时应满足公式(23)给出的取值范围，

xu,min≤xu⁰≤xu,max(u＝1,2...u)(23)

式中，xu⁰表示粒子u的初始位置，u表示求解日前调度模型时的粒子总个数，xu,min表示粒子u初始位置的最小值，xu,max表示粒子u初始位置的最大值，其中xu,min和xu,max根据电热联合系统中三大机组出力的上限值和下限值、蓄热式电锅炉出力的上限值和下限值确定；

8.1.2)计算每个粒子的适应值，并与粒子最优值和全局最优值作比较，更新粒子u的个体极值pbestu和粒子u的全局极值gbestu；

8.1.3)为了避免粒子群算法陷入全局最优或局部最优，利用公式(24)对ω^g进行动态调整，以保证算法的收敛性，

式中，ω^g表示第g次迭代的惯性权重，ωmax表示最大惯性权重，ωmax＝0.9，ωmin表示最小惯性权重，ωmin＝0.4，kmax表示求解日前调度模型时的迭代总次数；

8.1.4)利用公式(25)将第g次迭代的惯性权重ω^g代入速度更新公式中，更新微粒速度；

vud^g＝ω^gvud^g-1+c1rand(0,1)(pbestu-xu^g)+c2rand(0,1)(gbestu-xu^g)(25)

式中，vud^g表示粒子u第g次迭代的速度，vud^g-1表示粒子u第g-1次迭代的速度，c1、c2分别表示学习因子，c1＝c2＝2，pbestu表示粒子u的个体极值，gbestu表示粒子u的全局极值，xu^g表示粒子u第g次迭代的位置；

8.1.5)利用公式(26)迭代更新粒子的位置，

xu^g＝xu^g-1+vud^g(26)

式中，xu^g-1表示粒子u第g-1次迭代的位置，xu^g表示粒子u第g次迭代的位置，vud^g表示粒子u第g次迭代的速度；

8.1.6)重复执行步骤8.1.2)～步骤8.1.5)，直至达到求解日前调度模型时的迭代总次数kmax或全局最优位置满足最小界限，算法执行结束，输出在风电消纳情况下的n个常规机组的出力值、m个热电联产机组的出力值，蓄热式电锅炉的出力值以及风电机组的出力值；

其中采用改进的粒子群算法求解日内调度模型的具体步骤如下：

8.2.1)设置求解日内调度模型时的粒子的最大速度值vv,max、最大迭代次数k'max、种群规模f2，以及初始化粒子的位置xv⁰和速度vv⁰，为了缩短收敛时间，各粒子的初始位置选择时应满足公式(27)给出的取值范围，

xv,min≤xv⁰≤xv,max(v＝1,2...v)(27)

式中，xv⁰表示粒子v的初始位置，v表示求解日内调度模型时的粒子总个数，xv,min表示粒子v初始位置的最小值，xv,max表示粒子v初始位置的最大值，其中xv,min和xv,max根据电热联合系统中风电机组出力的上限值和下限值、可调用的激励型负荷的上限值plmax、下限值plmin，采暖负荷的最小值qt,min、最大值qt,max，电储能出力的上限值和下限值确定；

8.2.2)计算每个粒子的适应值，并与粒子最优值和全局最优值作比较，更新粒子v的个体极值pbestv和粒子v的全局极值gbestv；

8.2.3)为了避免粒子群算法陷入全局最优或局部最优，利用公式(28)对ω^h进行动态调整，以保证算法的收敛性，

式中，ω^h表示第h次迭代的惯性权重，ωmax表示最大惯性权重，ωmin表示最小惯性权重，k'max表示求解日内调度模型时的迭代总次数；

8.2.4)利用公式(29)将第h次迭代的惯性权重ω^h代入速度更新公式中，更新微粒速度；

vvd^h＝ω^hvvd^h-1+c3rand(0,1)(pbestv-xv^h)+c4rand(0,1)(gbestv-xv^h)(29)

式中，vvd^h表示粒子v第h次迭代的速度，vvd^h表示粒子v第h-1次迭代的速度，c3、c4分别表示学习因子，c3＝c4＝2，pbestv表示粒子v的个体极值，gbestv表示粒子v的全局极值，xv^h表示粒子v第h次迭代的位置；

8.2.5)利用公式(30)迭代更新粒子的位置，

xv^h＝xv^h-1+vvd^h(30)

式中，xv^h-1表示粒子v第h-1次迭代的位置，xv^h表示粒子v第h次迭代的位置，vvd^h表示粒子v第h次迭代的速度；

8.2.6)重复执行步骤8.2.2)～步骤8.2.5)，直至达到求解日内调度模型时的迭代总次数k'max或全局最优位置满足最小界限，算法执行结束，输出在风电消纳情况下调用的激励型负荷值pl1,t、采暖负荷的热功率值qt、电储能出力值以及风电机组的出力值。