[发明专利]考虑发电调节次数约束的在线滚动调度方法

权利要求书

1.一种考虑发电调节次数约束的在线滚动调度方法，其特征在于,该方法包括以下步骤：

1）设置滚动调度的时间长度为T_span，T_span的取值范围为1～4小时；设置滚动调度的时间长度内相邻的滚动调度的时间段的时间间隔长度为T_space，T_space的取值范围为5～15分钟；滚动调度的时段个数T如式（1）所示：

T=TspanTspace---(1)]]>

设置发电机的调整阈值系数r，r的取值范围一般为1~10；

设置电力系统中各台发电机在T个时段内包含的子时段个数S_i以及各个子时段的末尾时段其中S_i为第i台发电机在T个时段内包含的子时段个数，为第i台发电机在第j个子时段的末尾时段；设置障碍参数σ，其取值范围为0~1；设置迭代精度ε，其取值范围为1e-5至1e-3；

设置电力系统中所有发电机的原始发电计划值其中N_g为发电机的台数；为第i台发电机在第t个时段发出的有功功率的原始计划值；

2）计算所有的发电机的单位发电调整价格其中为第i台发电机在第t个时段的发电调整价格；在点处对第i台发电机在第t个时段的发电成本函数做泰勒展开，取值为泰勒展开式的一次项系数；

3）建立基于指数函数的考虑发电调节的次数约束以及时间间隔约束的在线滚动调度模型如下：

minxf(x)]]>

s.t.h(x)=0    （2）

g(x)≤g‾]]>

式（2）模型中，x为决策变量，包括：第i台发电机在第t个时段的向上发电调整量第i台发电机在第t个时段的向下发电调整量第i台发电机在第t个时段的调整量占爬坡速率的比值其中：

i＝1,2,...,N_g,t＝1,2,...,T；

f(x)为所有发电机在所有时段的发电调整价格之和的目标函数，如式（3）为：

f(x)=Σi=1NgΣt=1Tait(puit+pdit)---(3)]]>

等式约束h(x)=0包括：各个时段的发电负荷平衡约束，如式（4）所示：

Σi=1Ngpit-dt=0,∀t=1,2,...,T---(4)]]>

式（4）中，d_t为在第t个时段的系统负荷预测值；

调整后各台发电机在各个时间段发出的有功功率约束，如式（5）所示：

p0,it+puit-pdit-pit=0,∀t=1,2,...,T,i=1,2,...,Ng---(5)]]>

式（5）中，为调整后第i台发电机在第t个时段发出的有功功率；

式（2）模型中，不等式约束包括：

各台发电机的发电调整总次数约束：

Σt=1T(1-e-ruit)≤Ci,∀i=1,2,...,Ng---(6)]]>

式（6）中，C_i为第i台发电机在T个时段内调节发电次数的最大值；

各台发电机在各个子时段的发电调整次数约束，如式（7）所示：

Σt=sij-1+1sij(1-e-ruit)≤nij,∀i=1,2,...,Ng,j=1,2,...,Si---(7)]]>

式（7）中，为第i台发电机在第j个子时段内调节发电次数的最大值；

各台发电机在各个时段的爬坡速率约束，如式（8）所示：

-pit+pit-1-uitRDit-1≤0pit-pit-1-uitRUit-1≤0,∀t=1,2,...,T,i=1,2,...,Ng]]>

式（8）中，为第i台发电机在第t个时段的最大向下调节量，为第i台发电机在第t个时段的最大向上调节量；

各台发电机在各个时段发出的有功功率约束，如式（9）所示：

-pit≤-Pmin,itpit≤Pmax,it,∀t=1,2,...,T,i=1,2,...,Ng---(9)]]>

式（9）中，为第i台发电机在第t个时段发出的有功功率的最小值，为第i台发电机在第t个时段发出的有功功率的最大值；

4）求解考虑发电调节的次数约束以及时间间隔约束的在线滚动调度模型，具体过程如下：

4-1)构造所述在线滚动调度模型的增广拉格朗日函数，如式（10）所示

L(x,y,w,μ,u)=f(x)-yTh(x)-wT[g(x)+u-g‾]-μΣj=1rgln(uj)---(10)]]>

式（10）中，x，u为优化模型的原始变量，y，w为优化模型的对偶变量；其中，u为不等式约束的松弛变量，y为与等式约束对应的对偶变量，w为与不等式约束对应的对偶变量，r_g为矢量w的维度，μ为壁垒系数；各变量满足约束条件u≥0,w ≤0；

4-2）设置原始变量x，u、对偶变量y，w和壁垒系数μ的初始值，分别记为x⁽⁰⁾,u⁽⁰⁾，y⁽⁰⁾,w⁽⁰⁾和μ⁽⁰⁾，x⁽⁰⁾,y⁽⁰⁾,u⁽⁰⁾,w⁽⁰⁾,μ⁽⁰⁾一般取为0；设置迭代次数k=0；

4-3）计算拉格朗日增广函数对各个原始变量和对偶变量的一阶偏导数，如式（11）所示：

▿xL=▿xf(x(k))-▿xh(x(k))y(k)-▿xg(x(k))w(k)▿yL=h(x(k))▿wL=g(x(k))+u(k)-g‾▿uL=-w(k)-μ(k)U(k)-1e---(11)]]>

式（11）中，U(k)=diag(u1(k),u2(k),...,urg(k)),e=(1,1,...,1)rgT;]]>

4-4)建立如式（12）的修正方程，求解得到第k次迭代原始变量x的修正量Δx^(k)和对偶变量y的修正量Δy^(k)：

H′(k)▿xh(x(k))▿xTh(x(k))0Δx(k)Δy(k)=▿xL-▿yL---(12)]]>

式（12）中，H′(k)=-[▿x2f(x(k))-▿x2h(x(k))y(k)-▿x2g(x(k))w(k)]+▿xg(x(k))U(k)-1W(k)▿xTg(x(k)),]]>W(k)=diag(w1(k),w2(k),...,wrg(k));]]>

4-5)根据式（13）计算第k次迭代变量u,w的修正量Δu^(k),Δw^(k)：

Δu(k)=-▿xg(x(k))TΔx(k)-[g(x(k))+u(k)-g‾]]]>

Δw(k)=U(k)-1W(k)▿xg(x(k))TΔx(k)---(13)]]>

-U(k)-1{[U(k)W(k)e+μ(k)e]-W(k)[g(x(k))+u(k)-g‾]}]]>

4-6）根据式（14）计算第k次迭代原始变量的修正步长和对偶变量的修正步长为：

stepP(k)=0.9995min{mini=1,2,...,rg(-ui(k)Δui(k):Δui(k)<0),1}]]>（14）

stepD(k)=0.9995min{mini=1,2,...,rg(-wi(k)Δwi(k):Δwi(k)>0),1}]]>

4-7）根据式（15）对第k次迭代的原始变量和对偶变量进行修正：

x(k+1)u(k+1)=x(k)u(k)+stepP(k)Δx(k)Δu(k)]]>（15）

y(k+1)w(k+1)=y(k)w(k)+stepD(k)Δy(k)Δw(k)]]>

4-8）根据式（16）计算第k次迭代的原始-对偶互补间隙ρ^(k)和壁垒系数μ^(k)：

ρ^(k)=-u^(k)Tw

μ(k)=σρ(k)2rg---(16)]]>

4-9）对原始-对偶互补间隙ρ^(k)进行判断，若ρ^(k)<ε，则进行步骤5）；若ρ^(k)≥ε，则设置k=k+1，并进行步骤4-3）；

5）根据式（17）得到调整后各台发电机在各个时段发出的有功功率，如式(17)所示:

pit=p0,it+puit-pdit,∀t=1,2,...,T,i=1,2,...,Ng---(17);]]>

6）将调整后各台发电机在各个时段发出的有功功率下发到各个发电厂并执行，进行在线滚动调度。