[发明专利]一种考虑快速启停设备的多能虚拟电厂调节能力计算方法

权利要求书

1.一种考虑快速启停设备的多能虚拟电厂聚合方法，其特征在于，该方法包括以下步骤：

(1)记多能虚拟电厂调节能力的变量集合为其中，t为调度时段的编号，Υ为调度时段的集合，为多能虚拟电厂在调度时段t的联络线有功功率的基准值，和分别为多能虚拟电厂在调度时段t的联络线有功功率最大向上调节量和最大向下调节量；

(2)建立多能虚拟电厂调节能力计算模型的目标函数：

S_SD＝S_G∪S_CHP∪S_GB∪S_AC,S_QD＝S_EB∪S_EC

式中，Δt为相邻调度时段之间的时间间隔，i为多能虚拟电厂中任一设备的编号，S_SD为多能虚拟电厂内部所有传统火电机组、热电联产机组、燃气锅炉和吸收式机组共同构成的集合，S_G为多能虚拟电厂内部所有传统火电机组共同构成的集合，S_CHP为多能虚拟电厂内部所有热电联产机组共同构成的集合，S_GB为多能虚拟电厂内部所有燃气锅炉共同构成的集合，S_AC为多能虚拟电厂内部所有吸收式机组共同构成的集合，S_QD为多能虚拟电厂内部所有电锅炉和电制冷机共同构成的集合，S_EB为多能虚拟电厂内部所有电锅炉共同构成的集合，S_EC为多能虚拟电厂内部所有电制冷机共同构成的集合，为调度时段t时的购电价格，为多能虚拟电厂在调度时段t参与辅助服务市场向上调节单位有功功率的备用容量价格，为多能虚拟电厂在调度时段t参与辅助服务市场向下调节单位有功功率的备用容量价格，c^s,i为设备i的开机成本，为设备i的固定费用为表征设备i在调度时段t开机动作的0-1变量，设备i在调度时段t由关机状态转变为开机状态，则取值为1，其他情况下取值为0，为表征设备i在调度时段t运行状态的0-1变量，设备i在调度时段t处于开机状态时的取值为1，设备i在调度时段t处于关机状态时的取值为0，x_U为第一阶段决策变量和共同构成集合，i∈S_SD，即：

为调度时段t的上级电网可能向多能虚拟电厂下发的关口调度计划的有功功率值，满足：Ω₁为所有调度时刻的构成的集合，即cⁱ为传统火电机组、燃气锅炉和吸收式机组的运行成本系数，i∈S_G∪S_GB∪S_AC；对于传统火电机组，为传统火电机组i在调度时段t产生的有功功率；对于燃气锅炉，为燃气锅炉i在调度时段t产生的热功率；对于吸收式机组，为吸收式机组i在调度时段t产生的热功率；和分别为多能虚拟电厂在调度时段t参与辅助服务市场向上和向下调节单位有功功率的能量价格，和分别为多能虚拟电厂在调度时段t相对联络线有功功率的基准值向上调节和向下调节的有功功率值；c^E,i和c^H,i为热电联产机组i的成本系数，i∈S_CHP；和分别为热电联产机组i在调度时段t产生的有功功率和热功率，i∈S_CHP；x_L为第二阶段决策变量和共同构成集合，即：

其中，中i∈S_CHP，中i∈S_HL∪S_CL，中i∈S_EB∪S_EC，中i∈S_EB∪S_EC；

为热电联产机组i在调度时段t第k个组合系数，，i∈S_CHP；和分别为蓄电池i在调度时段t的充电功率和放电功率，为蓄电池i在调度时段t的电量，为向热负荷i消耗的热功率，为建筑i在调度时段t的室内温度；

(3)建立多能虚拟电厂调节能力计算模型的约束条件：

(3-1)建立设备运行状态的约束条件：

(3-2)建立设备运行的约束条件：

式中，对于传统火电机组，P^i,min和P^i,max分别为传统火电机组i在产生的有功功率的下限和上限；对于燃气锅炉，P^i,min和P^i,max分别为燃气锅炉i产生的热功率的下限和上限；对于吸收式机组，P^i,min和P^i,max分别为吸收式机组i产生的热功率的下限和上限；对于电锅炉，P^i,min和P^i,max分别为电锅炉i消耗的电功率的下限和上限；对于电制冷机，P^i,min和P^i,max分别为电制冷机i消耗的电功率的下限和上限；对于传统火电机组，为传统火电机组i在调度时段t-1产生的有功功率；对于燃气锅炉，为燃气锅炉i在调度时段t-1产生的热功率；对于吸收式机组，为吸收式机组i在调度时段t-1产生的热功率；对于电锅炉，为电锅炉i在调度时段t-1消耗的电功率；对于电制冷机，为电制冷机i在调度时段t-1消耗的电功率；R^i,up和R^i,down分别为设备i的向上爬坡速率和向下爬坡速率，R^i,su和R^i,sd分别为设备i的开机速率和关机速率；

(3-3)建立可再生能源机组运行的约束条件：

式中，为可再生能源机组i在调度时段t产生的有功功率的预测值，S_W为多能虚拟电厂内部所有可再生能源机组共同构成的集合；

(3-4)建立热电联产机组运行的约束条件：

式中，EPⁱ为热电联产机组i的可行域端点的集合，P^i,k和H^i,k分别为热电联产机组i的可行域的第k个端点的有功功率值和热功率值；

(3-5)建立蓄电池运行的约束条件：

式中，P^c,i,max和P^dc,i,max分别为蓄电池i充电功率的最大值和放电功率的最大值，E^i,min和E^i,max分别为蓄电池i电量的最小值和最大值，η^c,i和η^dc,i分别为蓄电池i充电效率和放电效率，e_t-1^ES,i为蓄电池i在调度时段t-1的电量；

(3-6)建立冷负荷和热负荷的约束条件：

式中，为冷负荷i或热负荷i的热容，为热负荷i或冷负荷i在调度时段t-1的室内温度，Uⁱ为热负荷i或冷负荷i的热导，为调度时段t的环境温度；τ^i,min和τ^i,max分别为热负荷i或冷负荷i室内温度的最小值和最大值，S_HL为多能虚拟电厂内部所有热负荷构成的集合，S_CL为能虚拟电厂内部所有冷负荷构成的集合；

(3-7)建立能量平衡的约束条件：

式中，η^EB,i为电锅炉i的制热系数，COPⁱ为电制冷机i的制冷系数，P_t^L,i为电负荷i在调度时段t消耗的有功功率，S_L为多能虚拟电厂内部所有电负荷构成的集合；

(3-8)建立不确定性的约束条件：

式中，为多能虚拟电厂在调度时段t参与辅助服务获得收益，为上级电网在调度时段t采用的多能虚拟电厂的联络线有功功率向上调节量的比例，为上级电网在调度时段t采用的多能虚拟电厂的联络线有功功率向下调节量的比例，为表征上级电网在调度时段t要求多能虚拟电厂的联络线有功功率向上调节和向下调节的0-1变量，当上级电网在调度时段t要求多能虚拟电厂的联络线有功功率向上调节时，的取值为1，当上级电网在调度时段t要求多能虚拟电厂的联络线有功功率向下调节时，的取值为0；

(4)建立一个由步骤(2)中建立的多能虚拟电厂调节能力计算模型的目标函数和步骤(3)中建立的多能虚拟电厂调节能力计算模型的约束条件共同构成多能虚拟电厂调节能力计算模型，该多能虚拟电厂调节能力计算模型中的决策变量包括：调度时段t的上级电网可能向多能虚拟电厂下发的关口调度计划的有功功率值多能虚拟电厂在调度时段t的联络线有功功率调节量的上限和下限传统火电机组i在调度时段t产生的有功功率开机动作变量和运行状态变量燃气锅炉i在调度时段t产生的热功率开机动作变量和运行状态变量吸收式机组i在调度时段t产生的热功率开机动作变量和运行状态变量热电联产机组i在调度时段t第k个组合系数开机动作变量和运行状态变量电锅炉i在调度时段t的开机动作变量电制冷机i在调度时段t的开机动作变量蓄电池i在调度时段t的充电功率和放电功率蓄电池i在调度时段t的电量在调度时段t向热负荷和冷负荷提供热功率热负荷和冷负荷在调度时段t的室内温度将该模型写为矩阵形式，即：

s.t.Ax+By≥q

Ω＝{u|Hu≥r}

Θ＝{v,z|Fv+Ez+Dx+Cy+Gu≥s,Tz≥p}

式中，y、x、u、z、v分别为由决策变量构成的列向量，其中，y由和共同构成，i∈S_SD,t∈Υ；x由和共同构成，t∈Υ；u由构成，t∈Υ；z由构成，i∈S_QD,t∈Υ；v由共同构成，中i∈S_QD,t∈Υ，中i∈S_CHP,t∈Υ；

cy+dx的具体表达式为：

其中，c为目标函数中y对应的系数矩阵，d为目标函数中x对应的系数矩阵；

hu+ez+fv的具体表达式为：

其中，h为目标函数中u对应的系数矩阵，e为目标函数中z对应的系数矩阵，f为目标函数中v对应的系数矩阵；

Θ为z和v的可行域；A、B、q、H、r、F、E、D、C、G、T、p是以(2)为目标函数以(3)为约束条件的模型转换为矩阵形式后，对应形式的约束条件分别对应的系数矩阵；

(5)利用改进的嵌套列约束生成算法对步骤(4)中建立的矩阵形式的多能虚拟电厂调节能力计算模型进行求解，具体步骤如下：

(5-1)记外层列约束生成算法的目标函数的下界和上界分别为LB_out和UB_out，设置LB_out的初始值为-∞，设置UB_out的初始值为∞；

(5-2)记外层列约束生成算法的迭代次数为l，设置l的初始值为0；

(5-3)利用分支定界算法求解外层主问题优化模型，外层主问题优化模型形式如下：

s.t.Ax+By≥q,

Fv^k+Ez^k+Dx+Cy+Gu^k≥s,k＝1,2,...,l

η≥hu^k+ez^k+fv^k,Tz^k≥p,k＝1,2,...,l

u^k∈Ω＝{u¹,u²,...u^NU},k＝1,2,...,l.

其中，η为外层列约束生成算法引入的辅助变量，记外层主问题优化模型的最优解为η^*、x^*和y^*，其中，η^*为η的最优解，x^*为x的最优解，y^*为y的最优解；

令LB_out＝cy^*+dx^*+η^*，对LB_out和UB_out之间的关系进行判断：如果不满足UB_out-LB_out≤ε_out，则进行步骤(5-4)；如果满足UB_out-LB_out≤ε_out，则计算结束，本次迭代计算获得的x^*和y^*即为多能虚拟电厂调节能力计算模型的最优解，获得多能虚拟电厂调节能力计算结果，结束计算；

(5-4)记内层列约束生成算法目标函数的下界和上届分别为LB_in和UB_in，设置LB_in的初始值为-∞，设置UB_in的初始值为∞，记内层列约束生成算法的迭代次数为m，设置m的初始值为0；

(5-5)利用分支定界算法求解对偶后的内层主问题优化模型，内层主问题优化模型如下：

s.t.Hu≥r,

θ≤ez^j+(s-Ez^j-Dx^*-Cy^*-Gu)^Tω^j,j＝1,2,...,NZ

ω^j≥0,j＝1,2,...,NZ

其中ω为引入的对偶变量构成的列向量，记内层主问题优化模型的最优解为和其中，为Θ的最优解，为u的最优解；

令对LB_in和UB_in之间的关系进行判断：如果满足UB_in-LB_in≤ε_in，则令并进行步骤(5-8)；如果不满足UB_in-LB_in≤ε_in，则进行步骤(5-6)；

(5-6)利用分支定界算法求解内层子问题优化模型，内层子问题优化模型如下：

Tz≥p,δ≥0,

其中δ为内层子问题中引入的松弛变量构成的列向量；J为松弛变量的系数矩阵，为一个对角矩阵，其取值规则为：对应约束(3-7)的对角元取值为1，其他值均为0；记内层子问题优化模型的最优解为z^*、v^*和δ^*，其中，z^*为z的最优解，v^*为v的最优解，δ^*为δ的最优解；

令LB_in＝max(LB_in,ez^*+fv^*+σ1^Tδ^*)，对LB_in和UB_in之间的关系进行判断：如果满足UB_in-LB_in≤ε_in，则令并进行步骤(5-8)；如果不满足UB_in-LB_in≤ε_in，则进行步骤(5-7)；

(5-7)建立新变量ω^m+1，ω^m+1为列约束生成算法添加的第m+1组对偶变量，建立下列约束条件，并添加至内层主问题优化模型：

θ≤(s-Ez^j-Dx^*-Cy^*-Gu)^Tω^m+1,

令m＝m+1，然后重新返回步骤(5-5)；

(5-8)令UB_out＝min(UB_out,cy^*+dx^*+UB_in)，对LB_out和UB_out之间的关系进行判断：如果满足UB_out-LB_out≤ε_out，则计算结束，本次迭代计算获得的x^*和y^*即为多能虚拟电厂调节能力计算模型的最优解，获得多能虚拟电厂调节能力计算结果，结束计算；

如果不满足UB_out-LB_out≤ε_out，则建立新变量z^l+1和v^l+1，向外层主问题优化模型增加如下列约束：

Fv^l+1+Ez^l+1+Dx+Cy+Gu^l+1≥s,

η≥hu^l+1+ez^l+1+fv^l+1,Tz^l+1≥p,

u^l+1＝u*.

令l＝l+1，然后重新返回步骤(5-3)。