龚锦霞, 郑元黎
(上海电力学院 电气工程学院, 上海 200090)
在可再生能源领域,风力发电因成本较低、技术成熟、可靠性较高,得到了大规模发展,而大规模并网风力发电的随机性给电力系统的经济调度和安全运行带来极大的挑战[1]。在大规模风电场并网的背景下,研究含随机因素的最优潮流具有重要意义[2-6]。风电功率的随机波动极有可能导致实时调度或控制时出现线路功率阻塞或节点电压过高等现象,而采用考虑随机因素最优潮流模型的调度方案对随机因素具有耐受性,可预防风电功率波动对电网的影响。
机会约束规划作为解决含随机变量的随机优化方法,近年来在电力系统领域得到应用。文献[7]引入机会约束规划,计算并网风电场的最大注入功率,并应用混合智能算法求解优化模型。文献[8]在机会约束规划的框架下确定了系统所需的最优旋转备用容量,并采用基于蒙特卡洛仿真的遗传算法进行求解。文献[9]建立了含风电场的随机动态最优潮流,采用随机模拟的混合智能技术来求解随机规划模型。文献[10]结合采样平均技术和遗传算法,求解考虑随机因素的机组组合问题。文献[11]基于机会约束模型研究了分布式电源的位置和容量优化问题,并采用Monte-Carlo抽样法和遗传算法相结合的方法进行求解。这些问题的求解过程中都需要最优潮流计算作为工具。
现有的文献中,基于机会约束规划的随机最优潮流(Stochastic Optimal Power Flow,SOPF)模型的求解一般采用智能优化算法[11]。对于含风电场的电力系统随机最优潮流,采用智能算法进行求解时,计算量较大。文献[12]中提出的改进原始对偶内点法可以求解目标函数分段可微的电力系统最优潮流计算问题,具有较快的计算速度和良好的收敛性。本文在文献[12]中内点法的基础上,加入与随机因素相关的调节因子,将随机规划问题转化为确定的非线性优化问题,迭代求解含风电的随机调度模型,并通过算例分析进行验证。
以发电机组报价之和最小为目标,则目标函数的公式为
(1)
式中:CG——可调节发电机的集合;
C2,i,C1,i,C0,i——节点i的发电成本的相关系数;
PGi——节点i的发电机有功功率。
等式约束条件为
(2)
式中:PWi,QWi——风电场的有功功率和无功功率;
PDi,QDi——负荷的有功功率和无功功率;
Ui,Uj——节点i和j的电压幅值;
Gij,Bij——节点i和j间的互导纳;
θij——节点i和j间的相角差;
QGi——节点i处发电机的无功功率。
不等式约束条件为
(3)
(4)
其中,
式中:PGi,max——机组最大调度出力;
P′Gi——前一时刻的机组出力;
Uup——有功出力上升速率;
PGi,min——机组最小调度出力;
Ddo——有功出力下降速率。
机会约束条件为
(5)
(6)
(7)
(8)
式中:P(·)——不等式成立的概率值;
βi——一组事先给定的置信水平,其大小反映了电力系统运行水平要求,βi∈[0,1]。
考虑负荷和风电场的功率波动,其公式分别为
(9)
(10)
式中:PD——负荷功率实际值;
PW——风电场的输出功率;
由于负荷和风电功率的计算是采用不同的预测方法,故认为预测误差相互独立。负荷的预测误差满足正态分布,风电功率的预测误差采用半不变量和Cornish-Fisher级数[13-14]表示为
(11)
式中:φ(x)——标准正态分布的概率密度函数;
总结式(1)~式(8),得到机会约束下的随机最优潮流模型为
(12)
式中:X=[PG,QG,θ,V]。
求解式(12)的难点在于机会约束的处理。通过增加调节因子,可将随机最优潮流模型转化为确定性优化问题。因此,式(12)中的机会约束可变为
C(X)≤a
(13)
式中:a——调节因子。
a表示随机因素对确定性优化问题的影响程度,也是随机因素和确定性表达的联系变量,a中的元素均大于等于零。由式(13)可以看出,a的取值反映了解空间的压缩程度。
根据式(12)和式(13)可得优化模型的改进拉格朗日函数为
(14)
式中:λ——等式约束的拉格朗日乘子;
μH,μC——不等式约束的拉格朗日乘子;
γ——扰动因子;
ZH,ZC——不等式约束的松弛因子;
Zm——松弛变量,Zm>0。
式(14)对应的KKT条件为
(15)
式中:fX=∂f/∂X;
GX=∂G/∂X;
HX=∂H/∂X;
CX=∂C/∂X;
μ=[μH;μC];
Z=[ZH;ZC]。
状态变量ΔX′和功率波动量的关系如下
ΔX′=J-1ΔW
(16)
式中:ΔX′=[Δθ,ΔV];
J——雅克比增广矩阵;
ΔW——负荷和风力发电机的功率预测误差。
不等式约束的波动量和ΔW的关系如下
ΔC=(CX·J-1)ΔW
(17)
式中:ΔC——C(X)的波动量。
ΔCj的第i阶半不变量的计算式为
(18)
式中:gi,ΔCj——ΔC的第j个元素的第i阶半不变量;
i——半不变量的阶数,i=1,2,3,…,N;
M=CX·J-1,其元素为Mi,j;
N——ΔW的元素个数;
gi,ΔWj——ΔW的第j个元素的第i阶半不变量。
(19)
将式(19)代入式(11),可得到Cj的概率密度函数fCj,k(c)。根据如下等价关系
(20)
式中:pj——满足P{Cj≤0}=βj的分位点。
当P{Cj≤0}=βj时:pj=0;pj>0表示概率越限,即不满足机会约束;pj≤0表示满足机会约束。
定义越界数列hβ为
(21)
(22)
式中:ρ——比例因子,取ρ=0.5。
根据式(15)和式(22),设定初始值k=0,a=0,则式(12)的求解流程如下。
步骤1k=k+1,按照下列公式计算优化变量
其中,ΔX,ΔZ,Δλ,Δμ的计算过程及符号意义参考文献[12]。加入a后,改进计算过程如下
ΔZ=-[H(X);C(X)]-Z-a-[HX;CX]ΔX
Δμ=-μ+[Z]-1(γe-[μ]ΔZ)
步骤2根据式(21)计算越界数列hβ,同时得到越界的机会约束集合Cy,k={Cj|hβ,j>0}。如果Cy,k=Ø,计算结束;否则,k=1时,采用式(22)更新调节因子;k>1时,如果Cy,k⊆Cy,k-1,采用式(16)更新调节因子,否则,ak=ρ·ak-1。
步骤3如果max(hβ)<ζ(ζ为一较小的数值),则计算结束;否则,如果k 3.1.1 算例数据 算例系统如图1所示。系统节点参数、负荷参数、发电机维修率及负荷半不变量数据均见文献[15-16]。 图1 算例系统示意 本算例选择节点9和12为风力发电场的并网点,当前风电场预测值分别为22 MW和20 MW。表1列出了以5 MW为统计区间的风电功率预测误差前四阶半不变量,只包括与本文算例相关的部分数据。功率基准值为100 MW。机会约束包括线路视在功率和节点电压幅值约束,设定置信水平均为95%。为说明线路传输功率越限,设置节点2和节点3之间的支路视在功率上限为50 MW;其他支路潮流的视在功率上限为200 MW。 表1 风电预测误差的半不变量统计(部分) 3.1.2 中间计算过程分析 采用本文所提的方法,终止条件为max(hβ)<5×10-5,经过12次迭代后可得到符合终止调节的解,迭代过程及计算量见表2。其中,Nk为第k次迭代的越界约束的个数,Ik为第k次迭代时内点法的迭代次数,Mk为第k次迭代的越界约束。 表2 迭代过程及计算量 表2中I1为第一步优化模型的迭代次数。Cy,1中含有如下7个元素 即当调节因子a=0时,支路2-3之间存在传输功率越限,节点7和节点9~12的电压大于其上限值。随着迭代过程中不断更新调节因子,越界约束的个数逐渐减小至零。a中非零元素对应于迭代中的越限机会约束,其中第1次迭代中有7个机会约束越限,迭代结果如图2所示。 图2 调节因子a的迭代结果 当某一机会约束满足后,其对应的调节因子的元素变为零,越界约束的个数减少,当迭代结束时为零。图3给出了迭代过程中7个机会约束对应的分位点的变化。从图3的结果还可以得到:若终止条件采用max(hβ)<1×10-4,则经过8次迭代可得到满意解。 图3 机会约束对应的分位点的变化 根据表2,图2和图3可以得出如下结论。整个计算过程中牛顿步迭代次数为172,是12次迭代中第1步优化模型的迭代次数之和;本算例中,C1和C7的变化对系统的总迭代次数影响较大,因此当迭代结束时,C1和C7小于零的概率接近给定的置信度。虽然最终解a中对应于C2~C6的元素为零,然而a中对应于C1和C7的元素的分量却增大,最终C2~C6满足机会约束。从表2可以看出,计算过程中Cy,k⊆Cy,k-1,即越限的机会约束数目保持不增,具有良好的收敛性。 3.1.3 最终结果分析 将本文所采用的随机最优潮流(SOPF)与确定性最优潮流方法(Deterministic Optimal Power Flow,DOPF)进行计算和比较:两种方法下的发电机出力调度方案、节点电压、支路3线路功率,以及C1~C7小于零的概率值分别如表3~6所示。 表3 发电机出力比较 表4 节点电压比较 表5 支路3线路功率比较 表6 C1~C7小于零的概率值 由表4可以看出,采用本文SOPF方法得到的节点7和节点9~12的电压幅值比DOPF方法的计算值小,可有效地预防风电场功率波动造成的节点电压越限。根据表5数据,采用DOPF计算可得到的调度方案中支路2-3的视在功率为49.988 MW,在实时调整中有可能出现线路阻塞的情况。采用SOPF方法,支路2-3的视在功率为43.433 MW,离设定的线路功率极限(50 MW)保持一定的裕度。表6表明,采用SOPF方法得到的解可有效地满足机会约束,即C1<0的概率大于等于0.95,并且C1~C7刚达到所设定的置信度,说明调节过程保证了调整量最小。 将风力发电机组接入节点9和节点26,预测功率分别为50 MW和40 MW,预测误差的概率参数如表1所示,并假设负荷的方差为其预测值的8%,约束支路8上的潮流上限设为400 MW。第1次迭代得到5个越限的机会约束数为 即支路5~8之间存在传输功率越限,节点3~5的电压大于其上限值。采用SOPF方法计算迭代10次后满足终止条件max(hβ)<5.0×10-5,调节过程如表7所示。 表7 118节点系统迭代过程及计算量 表8是两种方法下初次越界的机会约束最终的概率值。最优潮流和随机潮流的目标函数值分别为126 236.15元和126 360.57元。 由表8可以看出,118节点算例显示SOPF方法仍具有良好的收敛速度。 表8 C118,1~C118,5小于零的概率值 本文所有程序均采用Intel(R) Core(TM) i5-2410M CPU@ 2.3 GHz计算机,在MATLAB 2010a环境下运行,相关算法基于Matpower4.0中的子程序,计算时间如表9所示。 表9 算例计算时间 综上可见,与DOPF相比,通过SOPF方法得到的调度方案与电压和线路功率极限值保持了一定的距离,是更为安全的运行点,能经受和适应风电功率预测误差和负荷功率波动的影响。 随机最优潮流计算是考虑随机变量的电力系统分析的基础。本文基于机会约束规划建立含风电场的电力系统随机最优潮流模型,并采用内点法结合随机调节因子更新来求解该模型。利用分位点和机会约束概率之间的等价关系,在内点法中加入与分位点相关的调节因子,通过更新调节因子来调节机会约束满足的概率。该方法结合了内点法和半不变量理论的优点。算例分析表明,该方法可以有效求解含波动性电源和负荷的电力系统随机最优潮流。3 算例分析
3.1 4节点算例
3.2 118节点算例
4 结 语