规模化空调负荷参与电网调峰优化调度方法与流程

技术特征：

1.规模化空调负荷参与电网调峰优化调度方法，其特征在于：包括以下步骤，

步骤1，计及空调负荷本体运行特性，建立适用于电网调峰的空调负荷虚拟调峰机组模型；

步骤2，进行空调虚拟调峰机组的聚合；

步骤3，建立空调负荷虚拟调峰机组聚合层可调度容量及出力成本评估模型；

步骤4，空调负荷虚拟调峰机组聚合层间调峰容量分解；

步骤5，空调负荷虚拟调峰机组聚合层内优化调度。

2.根据权利要求1所述的规模化空调负荷参与电网调峰优化调度方法，其特征在于：用于电网调峰的空调负荷虚拟调峰机组模型为，

ACVPG＝{C(t),C(t,max),C(t,min),ρ(C(t)),D^min,D^max,IN}

其中，ACVPG表示空调负荷虚拟调峰机组，C(t)表示t时段空调负荷虚拟调峰机组出力，C(t,max)表示t时段空调负荷虚拟调峰机组出力上限，C(t,min)表示t时段空调负荷虚拟调峰机组出力下限，ρ(C(t))表示t时段空调负荷虚拟调峰机组出力成本，D^min表示空调负荷虚拟调峰机组最小调度持续时长，D^max表示空调负荷虚拟调峰机组最大调度持续时长，IN表示t时段空调负荷虚拟调峰机组最小调度间隔。

3.根据权利要求2所述的规模化空调负荷参与电网调峰优化调度方法，其特征在于：步骤1具体包括，

1.1)建立空调负荷等效热参数模型；

${\stackrel{\·}{T}}_a\left(t\right)+\frac{G_a}{C_a}\·T_a\left(t\right)=\frac{1}{C_a}\·\left(G_a{T}_o\right(t)+Q_h(t)+Q_m(t)+Q_s(t\left)\right)$

其中，C_a表示建筑物等效热容，G_a表示建筑物等效热导，T_a(t)表示t时段室内温度，T_o(t)表示t时段室外温度，Q_h(t)表示t时段空调制冷量，Q_m(t)表示t时段室内人员及设备所散发的热量，Q_s(t)表示t时段进入房间的太阳辐射热量；

1.2)建立空调负荷虚拟调峰机组出力模型；

$C\left(t\right)=P_{air}^0\left(t\right)-P_{air}\left(t\right)$

其中，表示空调负荷调控前基线功率，P_air(t)表示空调负荷调控后实际功率；

1.3)建立空调负荷虚拟调峰机组出力上下限模型；

空调负荷温度设定为约定最高值，最大降负荷量为空调负荷虚拟调峰机组出力上限，出力下限为0；

$C(t,max)=P_{air}^0\left(t\right)-P_{air}^{\min }\left(t\right)$

式中：表示空调负荷温度设定最高值并对空调运行参数优化后的总功耗；

1.4)建立空调负荷虚拟调峰机组出力成本模型；

$\mathrm{\ρ}\left(C\right(t\left)\right)=\underset{i^{\′}=0}{\overset{2}{\Σ}}{k}_{i^{\′}}{\[C\left(t\right)\]}^{i^{\′}}$

其中，k_i′表示空调负荷虚拟调峰机组出力成本系数，i′表示C(t)的次数；

1.5)建立空调负荷虚拟调峰机组最小调度持续时长模型；

D^min＝t_on+t_g

其中，t_on为空调负荷最小在线时长，t_g为用户附加最小调度时长；

1.6)建立空调负荷虚拟调峰机组最大调度持续时长模型；

对空调负荷进行调控后，室温上升到温度区间最大值的时长与人体对最高室温的忍耐时长之和，即为空调负荷虚拟调峰机组最大可调度持续时长；

D^max＝Δt+Δt′

$\mathrm{\Δ}t=\frac{C_a}{G_a}ln\frac{T_a^{set}-T_o\left(t\right)+\frac{1}{G_a}\left(Q_h\right(t)+Q_m(t)+Q_s(t\left)\right)}{T_a^{\max }-T_o\left(t\right)+\frac{1}{G_a}\left(Q_h\right(t)+Q_m(t)+Q_s(t\left)\right)}$

式中，Δt表示室温上升到温度区间最大值的时长，Δt′表示人体对温度区间最大值的忍耐时长，表示空调调控前的设定温度；

1.7)建立空调负荷虚拟调峰机组最小调度间隔模型；

调度后室温从温度区间最大值下降到温度区间最小值的时长为空调负荷虚拟调峰机组最小调度间隔；

$IN=\frac{C_a}{G_a}ln\frac{T_a^{max}-T_o\left(t\right)+\frac{1}{G_a}\left(Q_h\right(t)+Q_m(t)+Q_s(t\left)\right)}{T_a^{min}-T_o\left(t\right)+\frac{1}{G_a}\left(Q_h\right(t)+Q_m(t)+Q_s(t\left)\right)}.$

4.根据权利要求1所述的规模化空调负荷参与电网调峰优化调度方法，其特征在于：基于改进模糊C均值聚类算法进行空调负荷虚拟调峰机组聚合，即通过对传统模糊C均值聚类算法进行改进后，利用改进的模糊C均值聚类算法对空调负荷虚拟调峰机组进行聚类分组；

具体过程为，

2.1)改进模糊C均值聚类算法初始化，即设定聚类中心搜索终止条件δ、迭代终止条件ξ，权重系数m′；

2.2)计算空调负荷虚拟调峰机组数据集X＝{x₁,x₂,…,x_n}中所有数据点的密度指标；

$M_i^k\left(x_i\right)=\underset{l=1}{\overset{n}{\Σ}}\exp (-\frac{\left|\right|x_i-x_l|{|}^2}{{(r_a/2)}^2})$

其中，k表示密度指标计算次数，表示第k次计算中数据点x_i的密度指标，x_i∈X，r_a表示数据点x_i的邻域半径，l为整数，l∈[1,n]，l≠i，n为数据集中数据点的个数；

2.3)选出密度指标最大的数据点作为一个聚类中心，剔除该数据点后，计算剩余数据点的密度指标；

$M_i^k\left(x_i\right)=M_i^{k-1}\left(x_i\right)-M_c^{k-1}\underset{l=1}{\overset{n}{\Σ}}\exp (-\frac{\left|\right|x_i-x_c^{k-1}|{|}^2}{{(r_b/2)}^2})$

其中，表示第k-1次计算中数据点x_i的密度指标，c表示聚类中心，表示第k-1次计算中选出的聚类中心的密度指标，表示第k-1次计算中选作聚类中心的数据点，r_b表示密度指标减小的邻域半径；

2.4)判断是否成立，如果成立，则转至步骤2.5，如果不成立，则转至步骤2.3；

其中，表示第k次计算得到的聚类中心的密度指标，表示第1次计算得到的聚类中心的密度指标；

2.5)输出c个聚类中心；

2.6)计算空调负荷虚拟调峰机组调度特性指标初始隶属度矩阵U；

$u_{ij}=\frac{1}{\underset{k=1}{\overset{c}{\Σ}}{(d_{ij}/d_{kj})}^{2/(m^{\′}-1)}}$

其中，u_ij为隶属度矩阵U中的元素，表示第j个数据属于第i个类的隶属度，j为整数，j∈[1,n]，d_ij＝||x_j-v_i||表示数据点x_j与聚类中心v_i间的欧式距离；

2.7)更新空调负荷虚拟调峰机组调度特性指标聚类中心，即根据数据点隶属度矩阵U更新聚类中心矢量V；

$v_i=\frac{\underset{j=1}{\overset{n}{\Σ}}{\left(u_{ij}\right)}^{m^{\′}}{x}_j}{\underset{j=1}{\overset{n}{\Σ}}{\left(u_{ij}\right)}^{m^{\′}}},\∀i\∈c$

式中，v_i为聚类中心矢量V元素，x_j表示数据集中的数据点；

2.8)更新空调负荷虚拟调峰机组调度特性指标隶属度矩阵，即以更新后的聚类中心对隶属度矩阵进行更新；

2.9)判断终止条件||u_ij(k)-u_ij(k-1)||≤ξ是否成立，如果成立，则结束并输出结果；否则，依次重复步骤2.7和2.8；

其中，u_ij(k)表示第k次计算得到的第j个数据属于第i个类的隶属度值，u_ij(k-1)表示第k-1次计算得到的第j个数据属于第i个类的隶属度值。

5.根据权利要求1所述的规模化空调负荷参与电网调峰优化调度方法，其特征在于：空调负荷虚拟调峰机组聚合层可调度容量及出力成本评估模型为，

$C(g,max)=\underset{t=1}{\overset{T}{\Σ}}\underset{m=1}{\overset{M}{\Σ}}C(t,m)\·v(t,m)$

$\mathrm{\ρ}\left(C\right(t,g\left)\right)=\frac{\underset{m=1}{\overset{M}{\Σ}}\mathrm{\ρ}\left(C\right(t,m,max\left)\right)}{C(t,g,max)}$

其中，C(g,max)表示空调负荷虚拟调峰机组聚合层g在调度时段的最大可调度容量，M表示空调负荷虚拟机组聚合层g聚合的空调负荷虚拟调峰机组台数，m为整数，m∈[1,M]，T调度总时段数，C(t,m)表示调度时段t第m台空调负荷虚拟调峰机组的出力，v(t,m)表示调度时段t第m台空调负荷虚拟调峰机组是否参与调峰标识，ρ(C(t,g))表示空调负荷虚拟调峰机组聚合层g在调度时段t的出力成本，ρ(C(t,m,max))表示调度时段t第m台空调负荷虚拟调峰机组最大出力时的出力成本，C(t,m,max)表示空调负荷虚拟调峰机组聚合层g在调度时段t的最大可调度容量。

6.根据权利要求5所述的规模化空调负荷参与电网调峰优化调度方法，其特征在于：建立空调负荷虚拟调峰机组聚合层可调度容量评估模型的过程为，

3.1)建立可调度容量评估目标函数；

$\max C(g,max)=\underset{t=1}{\overset{T}{\Σ}}\underset{m=1}{\overset{M}{\Σ}}C(t,m)\·v(t,m)$

$minmax\left(\underset{m=1}{\overset{M}{\Σ}}C\left(t,m\right)v\left(t,m\right)\right)-\min \left(\underset{m=1}{\overset{M}{\Σ}}C\left(t,m\right)v\left(t,m\right)\right)$

3.2)明确评估约束条件；

具体包括：

空调负荷虚拟调峰机组出力约束；

C(t,m)≤C(t,m,max)

空调负荷虚拟调峰机组最大调度时长时长约束；

$\underset{k1=t}{\overset{t+D_m^{\max }}{\Σ}}v(k1,m)\≤D_m^{\max },\∀m\∈M,\∀t\∈\[1,T-D_m^{\max }\]$

$\underset{t=1}{\overset{D_m^{\max }-H_m^{init}+1}{\Σ}}v(t,m)\≤D_m^{max}-H_m^{init},\∀m\∈M$

其中，为第m台空调负荷虚拟调峰机组调度初始时段已调度时长；

空调负荷虚拟调峰机组最小调度持续时长约束；

$\underset{t=1}{\overset{{\mathrm{VG}}_m^{on}}{\Σ}}\[1-v(t,m)\]=0,\∀m\∈M$

$\underset{k2=t}{\overset{t+{\mathrm{IN}}_m^{on}-1}{\Σ}}v(k2,m)\≥{\mathrm{IN}}_m^{on}\[v(t,m)-v(t-1,m)\],\∀m\∈M,\∀t\∈\[{\mathrm{VG}}_m^{on}+1,T-{\mathrm{IN}}_m^{on}+1\]$

$\underset{k3=t}{\overset{T}{\Σ}}\{v(k3,m)-\[v(t,m)-v(t-1,m)\]\}\≥0,\∀m\∈M,\∀t\∈\[T-{\mathrm{IN}}_m^{on}+2,T\]$

${\mathrm{VG}}_m^{on}=min\[T,max\left(\right({\mathrm{IN}}_m^{on}-{\mathrm{IN}}_m^{inon}\left)v\right(0,m),0)\]$

其中，表示第m台空调负荷虚拟调峰机组在调度初始时段必须被调度的连续时长；表示第m台空调负荷虚拟调峰机组所允许的最小调度持续时长；表示第m台空调负荷虚拟调峰机组在调度初始时段已被连续调度时长；

空调负荷虚拟调峰机组最小调度间隔约束；

$\underset{t=1}{\overset{{\mathrm{VG}}_m^{off}}{\Σ}}v(t,m)=0,\∀m\∈M$

$\underset{k4=t}{\overset{t+{\mathrm{IN}}_m-1}{\Σ}}\[1-v(k4,m)\]\≥{\mathrm{IN}}_m\[v(t-1,m)-v(t,m)\],\∀m\∈M,\∀t\∈\[{\mathrm{VG}}_m^{off}+1,T-{\mathrm{IN}}_m+1\]$

$\underset{k3=t}{\overset{T}{\Σ}}\{1-v(k3,m)-\[v(t-1,m)-v(t,m)\]\}\≥0,\∀m\∈M,\∀t\∈\[T-{\mathrm{IN}}_m+2,T\]$

${\mathrm{VG}}_m^{off}=min\[T,max\left(\right({\mathrm{IN}}_m-{\mathrm{IN}}_m^{inoff}\left)\right(1-v\left(0,m\right)),0)\]$

其中，表示第m台空调负荷虚拟调峰机组在调度初始时段避免被调度的连续时长；表示第m台空调负荷虚拟调峰机组在调度初始时段未被调度的连续时长，k1、k2、k3、k4均为表示调度时段的中间变量，IN_m表示第m台空调负荷虚拟调峰机组最小调度间隔。

7.根据权利要求1所述的规模化空调负荷参与电网调峰优化调度方法，其特征在于：空调负荷虚拟调峰机组聚合层间调峰容量分解具体包括以下步骤，

4.1)建立空调负荷虚拟调峰机组聚合层间调峰容量分解目标函数

$\min F_p=\underset{t=1}{\overset{T}{\Σ}}\underset{g=1}{\overset{G}{\Σ}}\mathrm{\ρ}\left(C\right(t,g\left)\right)\·C(t,g)\·\mathrm{\Δ}t$

式中，F_p表示电网空调负荷调度总经济成本，g表示空调负荷虚拟调峰机组聚合层序号，G表示参与电网调峰的空调负荷虚拟调峰机组聚合层总数，C(t,g)表示t时段聚合层g的出力计划，Δt表示调度时段间隔；

4.2)明确调峰容量分解约束条件；

具体包括：

空调负荷虚拟调峰机组聚合层出力上下限约束；

0≤C(t,g)≤C(g,max)

式中，C(g,max)表示空调负荷虚拟调峰机组聚合层上报的最大可调度容量；

空调负荷虚拟调峰机组聚合层总体出力约束；

$P\left(t\right)\≤\underset{g=1}{\overset{G}{\Σ}}C(t,g)$

式中，P(t)表示时段t电网空调负荷削峰总容量需求。

8.根据权利要求1所述的规模化空调负荷参与电网调峰优化调度方法，其特征在于：空调负荷虚拟调峰机组聚合层内优化调度具体包括以下步骤，

5.1)建立空调负荷虚拟调峰机组聚合层内优化调度目标函数；

$\min F_a=\underset{t=1}{\overset{T}{\Σ}}\underset{m=1}{\overset{M}{\Σ}}\mathrm{\ρ}\left(C\right(t,m\left)\right)\·v(t,m)\·\mathrm{\Δ}t$

$min\mathrm{\Δ}C=\underset{t=1}{\overset{T}{\Σ}}|C(t,g)-\underset{m=1}{\overset{M}{\Σ}}C(t,m)|$

式中，F_a表示聚合层空调负荷调度总经济成本，C(t,m)表示调度时段t第m台空调负荷虚拟调峰机组出力，ΔC表示聚合层调峰偏差；

5.2)明确优化调度约束条件；

具体包括：

空调负荷虚拟调峰机组总体出力约束；

$C(t,g)\≤\underset{m=1}{\overset{M}{\Σ}}C(t,m)\·v(t,m),\∀t\∈T$

空调负荷虚拟调峰机组调度特性约束。