CN111799799B - 一种基于区间泰勒展开法的交直流混合配电网区间潮流计算方法 - Google Patents
一种基于区间泰勒展开法的交直流混合配电网区间潮流计算方法 Download PDFInfo
- Publication number
- CN111799799B CN111799799B CN202010673187.7A CN202010673187A CN111799799B CN 111799799 B CN111799799 B CN 111799799B CN 202010673187 A CN202010673187 A CN 202010673187A CN 111799799 B CN111799799 B CN 111799799B
- Authority
- CN
- China
- Prior art keywords
- power
- node
- distribution network
- current
- fluctuation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 72
- 238000000034 method Methods 0.000 title claims abstract description 62
- 239000011159 matrix material Substances 0.000 claims abstract description 156
- 238000011217 control strategy Methods 0.000 claims abstract description 27
- 230000003190 augmentative effect Effects 0.000 claims abstract description 19
- 238000002347 injection Methods 0.000 claims description 115
- 239000007924 injection Substances 0.000 claims description 115
- 230000005540 biological transmission Effects 0.000 claims description 37
- 230000005283 ground state Effects 0.000 claims description 17
- 238000010248 power generation Methods 0.000 claims description 17
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 239000011541 reaction mixture Substances 0.000 claims description 3
- 238000012423 maintenance Methods 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 12
- 238000000342 Monte Carlo simulation Methods 0.000 description 8
- 230000005622 photoelectricity Effects 0.000 description 5
- 238000005457 optimization Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006386 neutralization reaction Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/02—Circuit arrangements for ac mains or ac distribution networks using a single network for simultaneous distribution of power at different frequencies; using a single network for simultaneous distribution of ac power and of dc power
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及一种基于区间泰勒展开法的交直流混合配电网区间潮流计算方法,将区间泰勒展开法应用到交直流混合配电网中,结合增广直角坐标潮流模型,提出了一种适用于求解交直流混合配电网中存在间歇性电源以及波动性负荷的区间潮流计算方法。该方法考虑了主从控制策略与对等控制策略下直流配电网中节点电压和功率的波动对电压源换流器和交流配电网内潮流与调制比的影响,采用区间泰勒展开法避免了区间迭代与非线性规划,基于增广直角坐标潮流模型的雅克比矩阵与海森矩阵稀疏度极高,在IEEE33节点扩展的交直流混合配电网***中进行区间潮流计算,结果表明,本发明在增广直角坐标下用时0.0814秒,在直角坐标下用时0.4757秒,可以显著提高计算效率。
Description
技术领域
本发明涉及电力***技术领域,特别是一种基于区间泰勒展开法的交直流混合配电网区间潮流计算方法。
背景技术
随着世界各国对可再生能源的重视程度不断加强,风力发电,光伏发电技术逐渐成熟,并成功接入大电网参与***运行。风电,光电受天气因素影响较大,风力发电功率出力主要取决于风速,光伏发电功率出力主要取决于光照强度。由于天气随机性较强,造成发电功率出力存在不确定性,实际发电出力与预测值存在偏差,当偏差超出一定值时,功率输入大电网会引起电压波动,降低电能质量,进而影响设备运行。此外,波动性负荷接入电网同样会引起电压波动,因此为提前获取电压波动信息,对电网进行不确定性潮流运算是非常必要的。
随着交直流混合网络的发展,不确定性潮流算法的研究重点由交流网络转向交直流混合网络,限于电压源换流器多种控制方式的约束,当前关于交直流混合网络不确定性算法的研究并不多,主要包括非线性优化法与确定性等效法,其中非线性优化方法计算时间长,而确定性等效无法独立完成计算,在电压源换流器处需要结合非线性优化方法,同时无法获取节点电压相角的边界值。
增广直角坐标潮流模型中节点电压与节点注入电流皆为状态量,因此功率方程由全耦合的非线性方程转化为线性耦合的二次方程,基于这一特性,雅可比矩阵中的元素多为常数,在迭代计算中无需反复计算,同时也大大提高了雅可比矩阵与海森矩阵的稀疏率,提高潮流计算效率。
发明内容
有鉴于此,本发明的目的是提出一种基于区间泰勒展开法的交直流混合配电网改进区间潮流算法,避免了区间迭代与非线性优化,采用增广直角坐标潮流模型,雅可比矩阵与海森矩阵的稀疏率极高,在IEEE33节点扩展的交直流混合配电网***中进行区间潮流计算,结果表明,在增广直角坐标下用时0.0814秒,在直角坐标下用时0.4757秒,可以显著提高计算效率。
本发明采用以下方案实现:一种基于区间泰勒展开法的交直流混合配电网区间潮流计算方法,包括以下步骤:
步骤S1:针对交直流混合配电网的主从控制策略与对等控制策略分别确定每台VSC对应的控制方式,将VSC的多种控制方式与区间潮流算法相结合,分析区间潮流中不同控制方式下直流配电网节点电压,节点注入电流以及直流配电网与VSC间有功功率传输波动状况;
步骤S2:提出区间潮流中VSC调制比与VSC内部节点电压以及直流配电网节点电压应满足的电压关系式,确定直流配电网中的电压波动变量,分析电压波动变量对VSC内部节点电压和交流配电网节点电压以及调制比的影响;(直流侧还是交流侧)
步骤S3:根据交直流混合配电网中的节点电压方程,功率方程以及电压关系式列写增广直角坐标下的雅可比矩阵与海森矩阵;
步骤S4:选定负荷功率波动节点与光伏发电、风力发电接入节点,采用区间算法,确定节点功率波动量与区间中心值,建立波动变量矩阵;
步骤S5:采用交直流解耦潮流算法,首先进行直流配电网的区间潮流计算,获取直流配电网节点电压,节点注入电流区间变量对波动变量的一阶导数与二阶导数,计算直流配电网节点电压、节点电流与节点注入功率的区间值;
步骤S6:获取直流配电网节点电压与节点注入功率区间值,计及直流配电网输入交流配电网的功率值波动与直流配电网的电压波动状况对交流配电网的影响,选取交流配电网负荷功率波动节点,进行交流配电网的区间潮流计算。
进一步的,所述步骤S1的具体内容为:
交直流混合配电网的控制策略主要有主从控制策略与对等控制策略,其中,主从控制策略下,一台VSC作为主站,采用定直流电压控制方式,稳定直流电压值,其他VSC采用定功率控制方式,控制有功功率传输值;对等控制策略下,所有换流站皆采用下垂控制方式,同时调节直流配电网节点电压与节点注入电流有功功率值;
对于单台VSC,控制方式包括定直流电压、定交流电压、定有功或无功功率控制以及下垂控制方式;
定直流电压控制方式:
控制节点电压值为定值,区间潮流计算中功率波动引起该节点注入电流发生波动,因此节点注入功率为区间值;
上式中为直流配电网节点注入功率,为节点电压,为节点注入电流,为节点注入功率波动量,为节点注入电流波动量;无其他功率注入时,直流配电网注入VSC的功率为功率波动过程中,当波动量幅值较高时,功率传输方向会发生变化,发生整流与逆变状态的转换。直流配电网节点注入功率的波动引起VSC内部节点电压幅值与相角的波动,进而引起交流配电网内各状态量的波动;
定功率控制方式:
定功率控制方式下,其他节点功率注入值发生波动时,会引起该节点电压与节点注入电流发生微小波动,二者乘积计算后的功率波动值低于收敛精度,因此,在允许误差范围内认为节点注入功率为定值;由于采用交直流解耦潮流算法,定有功功率传输值改为定Pc,Pc为VSC与直流配电网间的有功功率传输值;无功功率传输值Qs维持定值,波动量设为0;
下垂控制方式:
直流配电网中连接VSC的节点出现电压上升或下降时,VSC根据电压偏差值决定注入节点的功率值;
根据区间算法,节点注入功率为:
定交流电压控制方式:
增广直角坐标下定交流电压控制方式维持节点i电压幅值恒为1,采用如下控制方式:
ei=1 fi=0 (5)
上式中ei为节点电压实部,fi为节点电压虚部;定交流电压方程中不存在功率波动变量,且雅可比矩阵中仅是关于ei的修正系数,因此,区间潮流计算中能够保持节点电压幅值恒为1,即电压幅值Ui=1。
进一步地,所述步骤S2的具体内容为:
结合区间泰勒展开法与交直流混合配电网的控制方式得到交直流混合配电网中采用区间泰勒展开法的潮流模型,VSC两侧需要满足电压关系式如下
上式中,Mk为第k台VSC电压调制比,采用标幺值,基准值其中,ei,fi为第k台VSC内部节点电压的实部与虚部,每台VSC包含一个内部节点,为直流配电网中与第k台VSC相连的节点电压值;控制方式为定直流电压控制时,节点电压为定值;控制方式为下垂控制与定有功功率控制时,节点电压为区间值,直流配电网在基态下运行时,节点电压在标幺值1附近,因此
上式中直流配电网中与VSC相连接节点的电压波动量[ΔVlowj,ΔVhj]对VSC内部电压ei,fi以及电压调制比Mk皆会产生影响,使VSC内部电压ei,fi以及电压调制比Mk发生波动,进而使交流配电网中连接VSC的节点,其节点电压,节点注入电流发生波动。
进一步地,所述步骤S3具体包括以下步骤:
步骤S31:建立交直流混合配电网区间潮流中的功率方程,节点电压方程以及VSC电压关系式,进行区间泰勒展开,建立不确定性潮流下的三个确定性等式;
考虑风电,光电等分布式电源的波动性以及负荷波动性的不确定潮流方程中,交流配电网中节点i功率波动表示为区间值:
[Pi]=Pic+[-ΔPi,ΔPi] (9)
[Qi]=Qic+[-ΔQi,ΔQi] (10)
Pic,Qic分别为有功,无功功率区间中心值,ΔPi,ΔQi分别为有功,无功功率波动量;
考虑光电等分布式电源的波动性以及负荷波动性的不确定潮流方程中,直流配电网中节点i功率波动表示为区间值;
考虑直流配电网中节点电压波动的不确定潮流方程中,直流配电网中与第k台VSC连接的节点j电压波动表示为区间值;
此外,直流配电网与VSC间的有功功率传输量受电压波动的影响也会发生波动,对VSC以及交流配电网中的潮流与调制比产生影响;因此,交直流混合配电网不确定性潮流中的波动变量包括分布式电源注入节点的有功,无功功率,波动性负荷以及直流配电网中与VSC相连的节点电压波动量与有功功率传输波动量,皆包含在矩阵ε内
上式中,εPi=[-ΔPi,ΔPi]表示节点有功功率波动量,有功功率波动量通过节点注入有功功率与节点有功功率负荷进行区间相加后得到,εQi=[-ΔQi,ΔQi]表示节点无功功率波动量,无功功率波动量通过节点注入无功功率与节点无功功率负荷进行区间相加后得到,i=1,2,…,n;表示直流配电网中与VSC相连节点的节点电压波动量,表示直流配电网与VSC间的有功功率传输波动量,k=1,2…L;交流配电网中节点电压方程与节点功率方程区间表达式为:
ΔPi=(Pic+εPi)-ei(ε)ai(ε)-fi(ε)bi(ε)=0 (15)
ΔQi=(Qic+εQi)-fi(ε)ai(ε)+ei(ε)bi(ε)=0 (16)
i=1,2,…,n,PV节点采用式(17)替代式(16);直流配电网中节点电压方程与节点功率方程区间表达式为:
i=1,2,…,r;定电压节点采用式(20)替代式(19);下垂控制方式中,直流配电网中与第k台VSC连接的节点j区间表达式如下:
第k台VSC电压关系式区间表达式如下:
联合公式(14)-(22),交直流混合配电网不确定潮流方程表示为矩阵方程:
不确定潮流矩阵方程ΔI(W(ε),ε)=0包含公式(14)和(18),ΔP(W(ε),ε)=0包含公式(15),(19)和(21),ΔQ(W(ε),ε)=0包含公式(16),ΔU(W(ε),ε)=0包含公式(17),(20)和(22)。其中W表示为:
W=[e1(ε),f1(ε),a1(ε),b1(ε),e2(ε),f2(ε),a2(ε),b2(ε)…M1(ε),M2(ε)…V1 d(ε),I1 d(ε),V2 d(ε),I2 d(ε)…]
(24)
直流配电网含有2r个状态量,交流配电网与VSC含有4n+L个状态量;
功率波动前***运行在基态值下,波动变量ε=εc=0,因此,采用区间泰勒展开法对功率不确定潮流方程在***基态运行值处进行区间泰勒展开,由于功率方程仅含节点电压,节点注入电流等状态量的二次项与波动变量的一次项,因此,将方程展开到二阶后,其高阶无穷小为0,节点电压方程仅包含状态量的一次项,其二阶导数为0,函数F(W(ε),ε)=0的泰勒展开式为:
式中Δε=[-Δε,Δε],要满足上述等式恒为0,只要满足:
F(W(εc),εc)=0 (26)
将函数的一阶导数公式(27),二阶导数公式(28)分别展开如下所示
所以,得到以下三个确定性方程
F(W(εc),εc)=0 (31)
直流配电网中,x=2r,y表示有功功率波动变量的数目;交流配电网与VSC中,x=4n+L,y表示有功,无功与电压波动变量的数目;
步骤S32:将三个确定性等式进行展开,计算功率函数对功率波动变量的导数以及VSC中电压关系式对直流配电网中电压波动量的导数;式(31)为波动变量ε=εc=0时的基态潮流计算,采用牛拉法进行潮流计算;
式(32)展开式如下:
上式中z为矩阵方程F(W(ε),ε)=0的数目,直流配电网中,z=x=2r,矩阵方程中F1-Fr为节点电压方程,Fr+1-F2r为有功功率方程;交流配电网与VSC中,z=x=4n+L,矩阵方程中F1-F2n为节点电压方程,F2n+1-F4n为有功,无功功率方程,F4n+1-F4n+L为VSC电压关系式方程;
矩阵方程对波动量的导数中,节点不发生功率波动时该导数为0,发生功率波动时该导数为波动量系数1。直流配电网中与第k台VSC相连的节点电压发生波动时,VSC电压关系式方程(22)对电压波动量的导数的值如下:
k=1,2…L,第k台VSC电压关系式是交流配电网矩阵方程F(W(ε),ε)=0的第4n+k个方程,Mkc为第k台VSC的电压调制比基态运行值;
式(33)表达式如下:
由于需要满足电压关系式(22),相较于纯交流配电网,连接直流配电网的交流配电网矩阵方程对波动量的二阶导数增加了矩阵当εi为直流配电网中节点i电压波动量且Wa同时为与节点i相连的VSC调制比时,交流配电网中第4n+k个方程,即第k台VSC电压关系式对波动量的二阶导数为非零元素,具体计算如下:
Wc为状态量在基态下的运行值,具体如下:
步骤S33:计算增广直角坐标潮流模型下的雅可比矩阵与海森矩阵,分析采用增广直角坐标潮流模型的雅可比矩阵与海森矩阵的稀疏度;
交流配电网中:
式(40)中矩阵各个元素如下:
式(41)-(43)中,an,bn分别为交流配电网节点n注入电流的实部与虚部,en,fn为节点n电压的实部与虚部,Gnn,Bnn为节点n自导纳的实部与虚部,G1n=Gn1为节点1与节点n互导纳的实部,二者相等,B1n=Bn1节点1与节点n互导纳的虚部,二者相等。
式(44)中,分别为第1台VSC内部节点电压的实部与虚部,下标i为该节点在交流网络中的节点标号,上标(1,2i-1)表示元素ei在该矩阵中的位置。分别为第L台VSC内部节点电压的实部与虚部,下标j为该节点在交流网络中的节点标号,上标同样表示位置;
直流配电网中:
直流配电网中含有r个节点;
上述矩阵AI,AV,DI,DV,I中非对角元素皆为0,因此雅可比矩阵中稀疏度很高;式(36)中关于增广直角坐标下的海森矩阵H(F),具体计算如下:
交流配电网中:
节点电压方程最高阶为1阶,海森矩阵H(F1)-H(F2n)皆为4n+L阶的0矩阵;
节点功率方程的海森矩阵H(F2n+1)-H(F4n)为稀疏矩阵,节点i的有功功率方程海森矩阵H(F2n+2i-1)为4n+L阶矩阵,其中矩阵中值为1,其他元素皆为0,无功功率方程海森矩阵稀疏度相同,i最高取值为n;
VSC电压关系式方程的海森矩阵H(F4n+1)-H(F4n+L)皆为4n+L阶矩阵,其中为矩阵中的三个非零元素,其他皆为0;Mj为第j台换流站调制比,ecj,fcj为第j台换流站交流侧电压实虚部,j最高取值为L;
直流配电网中:
节点电压方程海森矩阵H(F1)-H(Fr)为2r阶0矩阵;
进一步地,所述步骤S4的具体内容为:
交流配电网中考虑风电,光电等分布式发电波动性的不确定潮流方程中区间变量为功率值,其中节点i电源注入功率区间值为
考虑负荷功率波动性的不确定潮流方程中区间变量为功率值,其中节点i负荷功率区间值为
对节点有功,无功功率进行区间运算:
式(47)-(56)中,为节点i电源注入有功功率区间中心值,为节点i负荷有功功率区间中心值,节点i电源注入有功功率区间波动量,为节点i负荷有功功率区间波动量,节点i负荷无功功率区间中心值,为节点i电源注入无功功率区间中心值,为节点i负荷无功功率区间波动量,为节点i电源注入无功功率区间波动量;
直流配电网中考虑光电等分布式发电波动性的不确定潮流方程中区间变量为功率值,其中节点i电源注入功率区间值为
考虑负荷功率波动性的不确定潮流方程中区间变量为功率值,其中节点i负荷功率区间值为
对节点有功功率进行区间运算
将计算得到的交直流配电网中有功功率,无功功率波动变量代入式(13),节点注入有功功率,无功功率波动幅度为±20%,在功率基值Pic,Qic的基础上进行计算,得到节点有功功率,无功功率波动量,对波动变量εPi以及εQi赋值,得到波动变量矩阵ε,进而对交直流混合配电网进行区间潮流计算。
进一步地,所述步骤S5的具体内容为:
选取直流配电网中有功功率电源波动节点与负荷波动节点,计算节点有功功率波动幅值,节点注入功率波动幅度为±20%,在功率基值的基础上计算功率波动量;节点注入有功功率的基值为节点注入有功功率发生波动前的有功功率运行值;求得有功功率波动量后,对波动变量赋值,得到波动变量矩阵ε,此时波动变量矩阵ε仅包括直流配电网中的有功功率波动变量,为1×r阶;进行直流配电网的区间潮流计算,通过式(46)首先计算直流配电网的雅可比矩阵,进而根据式(34)求取状态量对波动变量的一阶导数,通过式(30)计算直流配电网的海森矩阵,进而根据式(36)求取状态量Vi d,Ii d对波动变量的二阶导数,,通过式(38)获取直流配电网节点电压,节点注入流的区间值。
进一步地,所述步骤S6的具体内容为:
步骤S5已获取直流配电网节点注入功率区间值,通过直流配电网区间潮流计算,得到了直流配电网中节点电压与节点注入电流的区间值;因此,通过式(3)获取直流配电网与VSC间有功功率传输区间值,根据式(9)计算得到有功功率传输波动量通过式(7),(8)获取直流配电网中与VSC相连接的节点电压波动量因此计及直流配电网输入交流配电网的有功功率值波动与直流配电网的电压波动状况对交流配电网的影响,选取交流配电网功率波动节点,节点有功,无功注入功率波动幅度皆为±20%,在功率基值Pic,Qic的基础上计算有功功率波动量εPi与无功功率波动量εQi;因此,直流配电网与VSC间的有功功率传输波动量直流配电网节点电压波动量交流配电网节点注入有功功率波动量εPi以及节点注入无功功率波动量εQi组成交流配电网的波动变量矩阵ε,为1×(2n+2L)阶矩阵;通过式(40)首先计算交流配电网的雅克比矩阵,进而根据式(34)求解交流配电网节点电压,节点注入电流以及电压调制比对波动变量以及εQi的一阶导数;通过式(30)首先计算海森矩阵,进而根据式(36)求解交流配电网节点电压,节点注入电流以及电压调制比对波动变量εPi以及εQi的二阶导数,最后通过式(38)求解交流配电网节点电压,节点注入电流以及电压调制比的区间值。与现有技术相比,本发明有以下有益效果:
本发明将区间泰勒展开法应用到交直流混合配电网中,结合增广直角坐标潮流模型,提出了一种适用于求解交直流混合配电网中存在间歇性电源以及波动性负荷的区间潮流计算方法。该方法实现了交直流混合配电网在主从控制策略与对等控制策略下的区间潮流计算,同时考虑到VSC对有功,无功功率的传输控制,以及VSC两侧应满足的电压关系式,能以较高的计算效率完成大规模交直流配电网的区间潮流计算。
附图说明
图1为本发明实施例的交直流配电网结构图。
图2为本发明实施例的主从控制下直流配电网节点电压幅值分布图。
图3为本发明实施例的主从控制下直流配电网节点注入电流幅值分布图。
图4为本发明实施例的主从控制下直流配电网节点注入功率分布图。
图5为本发明实施例的主从控制下交流配电网节点电压幅值分布图。
图6为本发明实施例的主从控制下交流配电网节点电压相角分布图。
图7为本发明实施例的对等控制下直流配电网电压幅值分布图。
图8为本发明实施例的对等控制下直流配电网节点注入功率分布图。
图9为本发明实施例的对等控制下交流配电网节点电压幅值分布图。
图10为本发明实施例的对等控制下交流配电网节点电压相角分布图。
具体实施方式
下面结合附图及实施例对本发明做进一步说明。
应该指出,以下详细说明都是例示性的,旨在对本申请提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
本实施例提供一种基于区间泰勒展开法的交直流混合配电网区间潮流计算方法,包括以下步骤:
步骤S1:针对交直流混合配电网的主从控制策略与对等控制策略分别确定每台VSC对应的控制方式,将VSC的多种控制方式与区间潮流算法相结合,分析区间潮流中不同控制方式下直流配电网节点电压,节点注入电流以及直流配电网与VSC间有功功率传输波动状况;
步骤S2:提出区间潮流中VSC调制比与VSC内部节点电压以及直流配电网节点电压应满足的电压关系式,确定直流配电网中的电压波动变量,分析电压波动变量对VSC内部节点电压和交流配电网节点电压以及调制比的影响;(直流侧还是交流侧)
步骤S3:根据交直流混合配电网中的节点电压方程,功率方程以及电压关系式列写增广直角坐标下的雅可比矩阵与海森矩阵;
步骤S4:选定负荷功率波动节点与光伏发电、风力发电接入节点,采用区间算法,确定节点功率波动量与区间中心值,建立波动变量矩阵;
步骤S5:采用交直流解耦潮流算法,首先进行直流配电网的区间潮流计算,获取直流配电网节点电压,节点注入电流区间变量对波动变量的一阶导数与二阶导数,计算直流配电网节点电压、节点电流与节点注入功率的区间值;
步骤S6:获取直流配电网节点电压与节点注入功率区间值,计及直流配电网输入交流配电网的功率值波动与直流配电网的电压波动状况对交流配电网的影响,选取交流配电网负荷功率波动节点,进行交流配电网的区间潮流计算。
在本实施例中,所述步骤S1的具体内容为:
交直流混合配电网的控制策略主要有主从控制策略与对等控制策略,其中,主从控制策略下,一台VSC作为主站,采用定直流电压控制方式,稳定直流电压值,其他VSC采用定功率控制方式,控制有功功率传输值;对等控制策略下,所有换流站皆采用下垂控制方式,同时调节直流配电网节点电压与节点注入电流有功功率值;
对于单台VSC,控制方式包括定直流电压、定交流电压、定有功或无功功率控制以及下垂控制方式;
定直流电压控制方式:
控制节点电压值为定值,区间潮流计算中功率波动引起该节点注入电流发生波动,因此节点注入功率为区间值;
上式中为直流配电网节点注入功率,为节点电压,为节点注入电流,为节点注入功率波动量,为节点注入电流波动量;无其他功率注入时,直流配电网注入VSC的功率为功率波动过程中,当波动量幅值较高时,功率传输方向会发生变化,发生整流与逆变状态的转换。直流配电网节点注入功率的波动引起VSC内部节点电压幅值与相角的波动,进而引起交流配电网内各状态量的波动;
定功率控制方式:
定功率控制方式下,其他节点功率注入值发生波动时,会引起该节点电压与节点注入电流发生微小波动,二者乘积计算后的功率波动值低于收敛精度,因此,在允许误差范围内认为节点注入功率为定值;由于采用交直流解耦潮流算法,定有功功率传输值改为定Pc,Pc为VSC与直流配电网间的有功功率传输值;无功功率传输值Qs维持定值,波动量设为0;
下垂控制方式:
直流配电网中连接VSC的节点出现电压上升或下降时,VSC根据电压偏差值决定注入节点的功率值;
根据区间算法,节点注入功率为:
定交流电压控制方式:
增广直角坐标下定交流电压控制方式维持节点i电压幅值恒为1,采用如下控制方式:
ei=1 fi=0 (5)
上式中ei为节点电压实部,fi为节点电压虚部;定交流电压方程中不存在功率波动变量,且雅可比矩阵中仅是关于ei的修正系数,因此,区间潮流计算中能够保持节点电压幅值恒为1,即电压幅值Ui=1。
在本实施例中,所述步骤S2的具体内容为:
结合区间泰勒展开法与交直流混合配电网的控制方式得到交直流混合配电网中采用区间泰勒展开法的潮流模型,VSC两侧需要满足电压关系式如下
上式中,Mk为第k台VSC电压调制比,采用标幺值,基准值其中,ei,fi为第k台VSC内部节点电压的实部与虚部,每台VSC包含一个内部节点,为直流配电网中与第k台VSC相连的节点电压值;控制方式为定直流电压控制时,节点电压为定值;控制方式为下垂控制与定有功功率控制时,节点电压为区间值,直流配电网在基态下运行时,节点电压在标幺值1附近,因此
上式中直流配电网中与VSC相连接节点的电压波动量[ΔVlowj,ΔVhj]对VSC内部电压ei,fi以及电压调制比Mk皆会产生影响,使VSC内部电压ei,fi以及电压调制比Mk发生波动,进而使交流配电网中连接VSC的节点,其节点电压,节点注入电流发生波动。
上述所提到的使VSC内部电压ei,fi以及电压调制比Mk发生波动,具体为在不考虑功率波动量的前提下,电压波动量与调制比呈反比,所以电压波动较高时,电压调制比Mk波动减小。VSC内部电压ei,fi与电压波动量呈正比或反比取决于ei,fi的正负,因此其波动大小不确定。电压波动量对交流配电网节点电压,节点注入电流的波动是间接影响,不存在绝对的正比或反比关系。
实际运行中,功率波动量与电压波动量同时作用,因此,VSC内部电压ei,fi以及电压调制比Mk、交流配电网节点电压,节点注入电流的波动是多个波动变量共同作用的结果。
在本实施例中,所述步骤S3具体包括以下步骤:
步骤S31:建立交直流混合配电网区间潮流中的功率方程,节点电压方程以及VSC电压关系式,进行区间泰勒展开,建立不确定性潮流下的三个确定性等式;
考虑风电,光电等分布式电源的波动性以及负荷波动性的不确定潮流方程中,交流配电网中节点i功率波动表示为区间值:
[Pi]=Pic+[-ΔPi,ΔPi] (9)
[Qi]=Qic+[-ΔQi,ΔQi] (10)
Pic,Qic分别为有功,无功功率区间中心值,ΔPi,ΔQi分别为有功,无功功率波动量;
考虑光电等分布式电源的波动性以及负荷波动性的不确定潮流方程中,直流配电网中节点i功率波动表示为区间值;
考虑直流配电网中节点电压波动的不确定潮流方程中,直流配电网中与第k台VSC连接的节点j电压波动表示为区间值;
此外,直流配电网与VSC间的有功功率传输量受电压波动的影响也会发生波动,对VSC以及交流配电网中的潮流与调制比产生影响;因此,交直流混合配电网不确定性潮流中的波动变量包括分布式电源注入节点的有功,无功功率,波动性负荷以及直流配电网中与VSC相连的节点电压波动量与有功功率传输波动量,皆包含在矩阵ε内
上式中,εPi=[-ΔPi,ΔPi]表示节点有功功率波动量,有功功率波动量通过节点注入有功功率与节点有功功率负荷进行区间相加后得到,εQi=[-ΔQi,ΔQi]表示节点无功功率波动量,无功功率波动量通过节点注入无功功率与节点无功功率负荷进行区间相加后得到,表示直流配电网中与VSC相连节点的节点电压波动量,表示直流配电网与VSC间的有功功率传输波动量;
k=1,2…L交流配电网中节点电压方程与节点功率方程区间表达式为:
ΔPi=(Pic+εPi)-ei(ε)ai(ε)-fi(ε)bi(ε)=0 (15)
ΔQi=(Qic+εQi)-fi(ε)ai(ε)+ei(ε)bi(ε)=0 (16)
i=1,2,…,n,PV节点采用式(17)替代式(16);直流配电网中节点电压方程与节点功率方程区间表达式为:
i=1,2,…,r;定电压节点采用式(20)替代式(19);下垂控制方式中,直流配电网中与第k台VSC连接的节点j区间表达式如下:
第k台VSC电压关系式区间表达式如下:
联合公式(14)-(22),交直流混合配电网不确定潮流方程表示为矩阵方程:
不确定潮流矩阵方程ΔI(W(ε),ε)=0包含公式(14)和(18),ΔP(W(ε),ε)=0包含公式(15),(19)和(21),ΔQ(W(ε),ε)=0包含公式(16),ΔU(W(ε),ε)=0包含公式(17),(20)和(22)。其中W表示为:
W=[e1(ε),f1(ε),a1(ε),b1(ε),e2(ε),f2(ε),a2(ε),b2(ε)…M1(ε),M2(ε)…V1 d(ε),I1 d(ε),V2 d(ε),I2 d(ε)…]
(24)
直流配电网含有2r个状态量,交流配电网与VSC含有4n+L个状态量;
功率波动前***运行在基态值下,波动变量ε=εc=0,因此,采用区间泰勒展开法对功率不确定潮流方程在***基态运行值处进行区间泰勒展开,由于功率方程仅含节点电压,节点注入电流等状态量的二次项与波动变量的一次项,因此,将方程展开到二阶后,其高阶无穷小为0,节点电压方程仅包含状态量的一次项,其二阶导数为0,函数F(W(ε),ε)=0的泰勒展开式为:
式中Δε=[-Δε,Δε],要满足上述等式恒为0,只要满足:
F(W(εc),εc)=0 (26)
将函数的一阶导数公式(27),二阶导数公式(28)分别展开如下所示
所以,得到以下三个确定性方程
F(W(εc),εc)=0 (31)
直流配电网中,x=2r,y表示有功功率波动变量的数目;交流配电网与VSC中,x=4n+L,y表示有功,无功与电压波动变量的数目;
步骤S32:将三个确定性等式进行展开,计算功率函数对功率波动变量的导数以及VSC中电压关系式对直流配电网中电压波动量的导数;式(31)为波动变量ε=εc=0时的基态潮流计算,采用牛拉法进行潮流计算;
式(32)展开式如下:
上式中z为矩阵方程F(W(ε),ε)=0的数目,直流配电网中,z=x=2r,矩阵方程中F1-Fr为节点电压方程,Fr+1-F2r为有功功率方程;交流配电网与VSC中,z=x=4n+L,矩阵方程中F1-F2n为节点电压方程,F2n+1-F4n为有功,无功功率方程,F4n+1-F4n+L为VSC电压关系式方程;
矩阵方程对波动量的导数中,节点不发生功率波动时该导数为0,发生功率波动时该导数为波动量系数1。直流配电网中与第k台VSC相连的节点电压发生波动时,VSC电压关系式方程(22)对电压波动量的导数的值如下:
k=1,2…L,第k台VSC电压关系式是交流配电网矩阵方程F(W(ε),ε)=0的第4n+k个方程,Mkc为第k台VSC的电压调制比基态运行值;
式(33)表达式如下:
由于需要满足电压关系式(22),相较于纯交流配电网,连接直流配电网的交流配电网矩阵方程对波动量的二阶导数增加了矩阵当εi为直流配电网中节点i电压波动量且Wa同时为与节点i相连的VSC调制比时,交流配电网中第4n+k个方程,即第k台VSC电压关系式对波动量的二阶导数为非零元素,具体计算如下:
Wc为状态量在基态下的运行值,具体如下:
步骤S33:计算增广直角坐标潮流模型下的雅可比矩阵与海森矩阵,分析采用增广直角坐标潮流模型的雅可比矩阵与海森矩阵的稀疏度;
交流配电网中:
式(40)中矩阵各个元素如下:
式(41)-(43)中,an,bn分别为交流配电网节点n注入电流的实部与虚部,en,fn为节点n电压的实部与虚部,Gnn,Bnn为节点n自导纳的实部与虚部,G1n=Gn1为节点1与节点n互导纳的实部,二者相等,B1n=Bn1节点1与节点n互导纳的虚部,二者相等。
式(44)中,别为第1台VSC内部节点电压的实部与虚部,下标i为该节点在交流网络中的节点标号,上标(1,2i-1)表示元素ei在该矩阵中的位置。分别为第L台VSC内部节点电压的实部与虚部,下标j为该节点在交流网络中的节点标号,上标同样表示位置;
直流配电网中:
直流配电网中含有r个节点;
上述矩阵AI,AV,DI,DV,I中非对角元素皆为0,因此雅可比矩阵中稀疏度很高;式(36)中关于增广直角坐标下的海森矩阵H(F),具体计算如下:
交流配电网中:
节点电压方程最高阶为1阶,海森矩阵H(F1)-H(F2n)皆为4n+L阶的0矩阵;
节点功率方程的海森矩阵H(F2n+1)-H(F4n)为稀疏矩阵,节点i的有功功率方程海森矩阵H(F2n+2i-1)为4n+L阶矩阵,其中矩阵中值为1,其他元素皆为0,无功功率方程海森矩阵稀疏度相同,i最高取值为n;
VSC电压关系式方程的海森矩阵H(F4n+1)-H(F4n+L)皆为4n+L阶矩阵,其中为矩阵中的三个非零元素,其他皆为0;Mj为第j台换流站调制比,ecj,fcj为第j台换流站交流侧电压实虚部,j最高取值为L;
直流配电网中:
节点电压方程海森矩阵H(F1)-H(Fr)为2r阶0矩阵;
在本实施例中,所述步骤S4的具体内容为:
交流配电网中考虑风电,光电等分布式发电波动性的不确定潮流方程中区间变量为功率值,其中节点i电源注入功率区间值为
考虑负荷功率波动性的不确定潮流方程中区间变量为功率值,其中节点i负荷功率区间值为
对节点有功,无功功率进行区间运算:
式(47)-(56)中,为节点i电源注入有功功率区间中心值,为节点i负荷有功功率区间中心值,ΔPi S为节点i电源注入有功功率区间波动量,ΔPi L为节点i负荷有功功率区间波动量,为节点i负荷无功功率区间中心值,为节点i电源注入无功功率区间中心值,为节点i负荷无功功率区间波动量,为节点i电源注入无功功率区间波动量;
直流配电网中考虑光电等分布式发电波动性的不确定潮流方程中区间变量为功率值,其中节点i电源注入功率区间值为
考虑负荷功率波动性的不确定潮流方程中区间变量为功率值,其中节点i负荷功率区间值为
对节点有功功率进行区间运算
将计算得到的交直流配电网中有功功率,无功功率波动变量代入式(13),节点注入有功功率,无功功率波动幅度为±20%,在功率基值Pic,Qic的基础上进行计算,得到节点有功功率,无功功率波动量,对波动变量εPi以及εQi赋值,得到波动变量矩阵ε,进而对交直流混合配电网进行区间潮流计算,具体计算过程在步骤S5,步骤S6中进行详细介绍。
在本实施例中,所述步骤S5的具体内容为:
选取直流配电网中有功功率电源波动节点与负荷波动节点,计算节点有功功率波动幅值,节点注入功率波动幅度为±20%,在功率基值的基础上计算功率波动量;节点注入有功功率的基值为节点注入有功功率发生波动前的有功功率运行值;求得有功功率波动量后,对波动变量赋值,得到波动变量矩阵ε,此时波动变量矩阵ε仅包括直流配电网中的有功功率波动变量,为1×r阶;进行直流配电网的区间潮流计算,通过式(46)首先计算直流配电网的雅可比矩阵,进而根据式(34)求取状态量Vi d,对波动变的一阶导数,通过式(30)计算直流配电网的海森矩阵,进而根据式(36)求取状态量对波动变量的二阶导数,,通过式(38)获取直流配电网节点电压,节点注入电流的区间值。
在本实施例中,所述步骤S6的具体内容为:
步骤S5已获取直流配电网节点注入功率区间值,通过直流配电网区间潮流计算,得到了直流配电网中节点电压与节点注入电流的区间值;因此,通过式(3)获取直流配电网与VSC间有功功率传输区间值,根据式(9)计算得到有功功率传输波动量通过式(7),(8)获取直流配电网中与VSC相连接的节点电压波动量因此计及直流配电网输入交流配电网的有功功率值波动与直流配电网的电压波动状况对交流配电网的影响,选取交流配电网功率波动节点,节点有功,无功注入功率波动幅度皆为±20%,在功率基值Pic,Qic的基础上计算有功功率波动量εPi与无功功率波动量εQi;因此,直流配电网与VSC间的有功功率传输波动量直流配电网节点电压波动量交流配电网节点注入有功功率波动量εPi以及节点注入无功功率波动量εQi组成交流配电网的波动变量矩阵ε,为1×(2n+2L)阶矩阵;通过式(40)首先计算交流配电网的雅克比矩阵,进而根据式(34)求解交流配电网节点电压,节点注入电流以及电压调制比对波动变量以及εQi的一阶导数;通过式(30)首先计算海森矩阵,进而根据式(36)求解交流配电网节点电压,节点注入电流以及电压调制比对波动变量εPi以及εQi的二阶导数,最后通过式(38)求解交流配电网节点电压,节点注入电流以及电压调制比的区间值。较佳的,在本实施例中,一具体实例如下:本实施例对IEEE33节点扩展的交直流混合配电网进行了区间潮流计算,在处理器Intel(R)Core(TM)i5-3230CPU下采用Matlab2014a编程。***基准容量SB=100MVA,交流侧与直流侧基准电压为UB=10kV。交流配电网中接入风力发电,直流配电网中接入光伏发电,具体接入位置如图1所示。功率波动幅值设定±20%,不确定功率包括风力发电机发出的有功功率与无功功率,光伏发电发出的有功功率,以及负荷功率波动,其中交流配电网中负荷功率波动选定节点2、3、4、5、9、10、11、16、17、19、22、27、29、31,直流配电网负荷功率不发生波动。采用本实施例提出的方法与蒙特卡罗法进行比较,其中蒙特卡罗法区间内默认为均匀分布,计算次数为10000次。
电压源换流器内主要参数包括阻抗值Zc与电纳值Bf,其具体数值在表1中给出。同时给出了交直流混合配电网采用主从控制策略与对等控制策略时电压源换流器采取的控制方式。
表1电压源换流器参数及控制方式
1主从控制
1)直流配电网区间潮流计算结果分析
IEEE33节点扩展的交直流配电网的直流配电网包含6个节点,节点2与节点4连接电压源换流器,节点5与节点6连接光伏发电,波动范围皆为±20%。电压源换流器1为定控制方式,电压源换流器2为定Pc控制方式,在主从控制方式场景下分别采用区间泰勒展开法与蒙特卡罗法进行直流配电网的区间潮流计算。
由图2可以看出采用区间泰勒展开法得到的节点电压幅值计算结果对蒙特卡罗法的计算结果具有完全的包裹性,可以证明该方法的有效性。同时对节点注入电流进行了区间潮流计算,计算结果如图3所示,图4所示为根据节点电压与节点注入电流的区间值得到的直流配电网节点注入功率的波动值,电压源换流器1向节点2输入的有功功率在[0.0071,0.0111]范围内波动,当直流配电网中波动幅度较大时,功率传输方向会发生改变。节点4为定有功功率控制方式,节点4电压在[0.9917,0.9947]范围内波动,与电流值相乘后获得的节点注入功率在[0.000037,0.000058]范围内,上下限功率差值极小,低于收敛精度0.0001。
2)交流配电网区间潮流计算结果分析
交流配电网中节点26处接入风力发电机,同时选取12个节点的有功功率与无功功率负荷,波动范围皆为±20%,两台电压源换流器直流侧皆受直流配电网波动量的影响,因此,交流配电网部分含有28个波动变量,分别采用区间泰勒展开法与蒙特卡罗法进行区间潮流计算。
由图5、6可以看出,区间泰勒展开法的节点电压幅值计算结果对蒙特卡罗法的计算结果具有完全的包裹性,可以证明该方法的有效性。节点1为平衡节点,主从控制方式场景中,节点14连接电压源换流器2,控制方式为定Us,因此节点1与节点14电压幅值始终不变,节点1与节点14附近节点电压幅值波动范围较小,PV节点可有效抑制电压波动量,远离PV节点的节点电压幅值波动较剧烈。
电压源换流器调制比M1,M2的调制比标幺值皆为1,区间潮流计算中,随着电压值的波动,调制比为满***直流配电网的电压关系,也随之波动,计算结果如下
表2电压调制比分布
2对等控制
对等控制方式场景中,下垂控制参数kdroop=-50,直流配电网与交流配电网中的波动变量不变,分别采用区间泰勒展开法与蒙特卡罗法进行区间潮流计算。
1)直流配电网区间潮流计算
由图7可以看出,区间泰勒展开法的计算结果对蒙特卡罗法计算结果同样具有完全的包裹性,波动量相同时,相比于主从控制方式,对等控制方式的节点4电压波动幅度有效减小。
节点注入功率与节点电压差值呈线性关系,当电压值低于电压标准值时,配电网根据下垂控制方式增加节点功率注入量,抬高电压值。
图8为节点注入功率的波动状况,由节点2与节点4可以看出,电压达到上界时,对应节点注入功率的下界,因此,采用区间泰勒展开法获得的结果完全合理。
2)交流配电网区间潮流计算
采用对等控制方式时,直流配电网中节点2与节点4的电压值与注入功率皆为波动量,输入到交流配电网与电压源换流器中,因此,进行交流配电网的区间潮流计算时,波动变量在原基础上增加了2个,为30个波动变量。
与主从控制方式相比,参照电压波动图图3与图8,功率波动图图5与图9可以看出,对等控制方式下节点4电压与节点2注入功率的波动范围均缩小,同时由于对等控制方式下节点电压与节点注入功率皆发生了波动,因此交流配电网中连接电压源换流器1的节点7与节点34,其节点电压幅值波动范围缩小,电压源换流器2采用定交流电压控制,因此电压幅值波动不受影响。
如图10电压源换流器2相较于主从控制中的定功率控制方式,对等控制方式下功率波动范围增加,因此节点35功角波动范围增大。
表3电压调制比分布
相较于主从控制中直流配电网节点2的定电压控制方式,对等控制方式下节点2电压波动范围增大,相较于主从控制中直流配电网节点4的定功率控制方式,对等控制方式下节点4电压波动范围减小。
以上所述仅为本发明的较佳实施例,凡依本发明申请专利范围所做的均等变化与修饰,皆应属本发明的涵盖范围。
Claims (6)
1.一种基于区间泰勒展开法的交直流混合配电网区间潮流计算方法,其特征在于,包括以下步骤:
步骤S1:针对交直流混合配电网的主从控制策略与对等控制策略分别确定每台VSC对应的控制方式,将VSC的多种控制方式与区间潮流算法相结合,分析区间潮流中不同控制方式下直流配电网节点电压,节点注入电流以及直流配电网与VSC间有功功率传输波动状况;
步骤S2:提出区间潮流中VSC调制比与VSC内部节点电压以及直流配电网节点电压应满足的电压关系式,确定直流配电网中的电压波动变量,分析电压波动变量对VSC内部节点电压和交流配电网节点电压以及调制比的影响;
步骤S3:根据交直流混合配电网中的节点电压方程,功率方程以及电压关系式列写增广直角坐标下的雅可比矩阵与海森矩阵;
步骤S4:选定负荷功率波动节点与光伏发电、风力发电接入节点,采用区间算法,确定节点功率波动量与区间中心值,建立波动变量矩阵;
步骤S5:采用交直流解耦潮流算法,首先进行直流配电网的区间潮流计算,获取直流配电网节点电压,节点注入电流区间变量对波动变量的一阶导数与二阶导数,计算直流配电网节点电压、节点电流与节点注入功率的区间值;
步骤S6:获取直流配电网节点电压与节点注入功率区间值,计及直流配电网输入交流配电网的功率值波动与直流配电网的电压波动状况对交流配电网的影响,选取交流配电网负荷功率波动节点,进行交流配电网的区间潮流计算;
其中,所述步骤S3具体包括以下步骤:
步骤S31:建立交直流混合配电网区间潮流中的功率方程,节点电压方程以及VSC电压关系式,进行区间泰勒展开,建立不确定性潮流下的三个确定性等式;
考虑分布式电源的波动性以及负荷波动性的不确定潮流方程中,交流配电网中节点i功率波动表示为区间值:
[Pi]=Pic+[-△Pi,△Pi] (9)
[Qi]=Qic+[-△Qi,△Qi] (10)
Pic,Qic分别为有功,无功功率区间中心值,△Pi,△Qi分别为有功,无功功率波动量;
考虑分布式电源的波动性以及负荷波动性的不确定潮流方程中,直流配电网中节点i功率波动表示为区间值;
考虑直流配电网中节点电压波动的不确定潮流方程中,直流配电网中与第k台VSC连接的节点j电压波动表示为区间值;
此外,直流配电网与VSC间的有功功率传输量受电压波动的影响也会发生波动,对VSC以及交流配电网中的潮流与调制比产生影响;因此,交直流混合配电网不确定性潮流中的波动变量包括分布式电源注入节点的有功,无功功率,波动性负荷以及直流配电网中与VSC相连的节点电压波动量与有功功率传输波动量,皆包含在矩阵ε内;
上式中,εPi=[-△Pi,△Pi]表示节点有功功率波动量,有功功率波动量通过节点注入有功功率与节点有功功率负荷进行区间相加后得到,εQi=[-△Qi,△Qi]表示节点无功功率波动量,无功功率波动量通过节点注入无功功率与节点无功功率负荷进行区间相加后得到,i=1,2,…,n;表示直流配电网中与VSC相连节点的节点电压波动量,表示直流配电网与VSC间的有功功率传输波动量,k=1,2…L;交流配电网中节点电压方程与节点功率方程区间表达式为:
△Pi=(Pic+εPi)-ei(ε)ai(ε)-fi(ε)bi(ε)=0 (15)
△Qi=(Qic+εQi)-fi(ε)ai(ε)+ei(ε)bi(ε)=0 (16)
i=1,2,…,n,PV节点采用式(17)替代式(16);直流配电网中节点电压方程与节点功率方程区间表达式为:
i=1,2,…,r;定电压节点采用式(20)替代式(19);下垂控制方式中,直流配电网中与第k台VSC连接的节点j区间表达式如下:
第k台VSC电压关系式区间表达式如下:
联合公式(14)-(22),交直流混合配电网不确定潮流方程表示为矩阵方程:
不确定潮流矩阵方程△I(W(ε),ε)=0包含公式(14)和(18),△P(W(ε),ε)=0包含公式(15),(19)和(21),△Q(W(ε),ε)=0包含公式(16),△U(W(ε),ε)=0包含公式(17),(20)和(22);其中W表示为:
直流配电网含有2r个状态量,交流配电网与VSC含有4n+L个状态量;
功率波动前***运行在基态值下,波动变量ε=εc=0,因此,采用区间泰勒展开法对功率不确定潮流方程在***基态运行值处进行区间泰勒展开,由于功率方程仅含节点电压,节点注入电流状态量的二次项与波动变量的一次项,因此,将方程展开到二阶后,其高阶无穷小为0,节点电压方程仅包含状态量的一次项,其二阶导数为0,函数F(W(ε),ε)=0的泰勒展开式为:
式中△ε=[-△ε,△ε],要满足上述等式恒为0,只要满足:
F(W(εc),εc)=0 (26)
将函数的一阶导数公式(27),二阶导数公式(28)分别展开如下所示:
所以,得到以下三个确定性方程:
F(W(εc),εc)=0 (31)
直流配电网中,x=2r,y表示有功功率波动变量的数目;交流配电网与VSC中,x=4n+L,y表示有功,无功与电压波动变量的数目;
步骤S32:将三个确定性等式进行展开,计算功率函数对功率波动变量的导数以及VSC中电压关系式对直流配电网中电压波动量的导数;式(31)为波动变量ε=εc=0时的基态潮流计算,采用牛拉法进行潮流计算;
式(32)展开式如下:
上式中z为矩阵方程F(W(ε),ε)=0的数目,直流配电网中,z=x=2r,矩阵方程中F1-Fr为节点电压方程,Fr+1-F2r为有功功率方程;交流配电网与VSC中,z=x=4n+L,矩阵方程中F1-F2n为节点电压方程,F2n+1-F4n为有功,无功功率方程,F4n+1-F4n+L为VSC电压关系式方程;
矩阵方程对波动量的导数中,节点不发生功率波动时该导数为0,发生功率波动时该导数为波动量系数1;直流配电网中与第k台VSC相连的节点电压发生波动时,VSC电压关系式方程(22)对电压波动量的导数的值如下:
k=1,2…L,第k台VSC电压关系式是交流配电网矩阵方程F(W(ε),ε)=0的第4n+k个方程,Mkc为第k台VSC的电压调制比基态运行值;
式(33)表达式如下:
由于需要满足电压关系式(22),相较于纯交流配电网,连接直流配电网的交流配电网矩阵方程对波动量的二阶导数增加了矩阵当εi为直流配电网中节点i电压波动量且Wa同时为与节点i相连的VSC调制比时,交流配电网中第4n+k个方程,即第k台VSC电压关系式对波动量的二阶导数为非零元素,具体计算如下:
Wc为状态量在基态下的运行值,具体如下:
步骤S33:计算增广直角坐标潮流模型下的雅可比矩阵与海森矩阵,分析采用增广直角坐标潮流模型的雅可比矩阵与海森矩阵的稀疏度;
交流配电网中:
式(40)中矩阵各个元素如下:
式(41)-(43)中,an,bn分别为交流配电网节点n注入电流的实部与虚部,en,fn为节点n电压的实部与虚部,Gnn,Bnn为节点n自导纳的实部与虚部,G1n=Gn1为节点1与节点n互导纳的实部,二者相等,B1n=Bn1节点1与节点n互导纳的虚部,二者相等;
式(44)中,fi (1,2i)分别为第1台VSC内部节点电压的实部与虚部,下标i为该节点在交流网络中的节点标号,上标(1,2i-1)表示元素ei在该矩阵中的位置;分别为第L台VSC内部节点电压的实部与虚部,下标j为该节点在交流网络中的节点标号,上标同样表示位置;
直流配电网中:
直流配电网中含有r个节点;
上述矩阵AI,AV,DI,DV,I中非对角元素皆为0,因此雅可比矩阵中稀疏度很高;式(36)中关于增广直角坐标下的海森矩阵H(F),具体计算如下:
交流配电网中:
节点电压方程最高阶为1阶,海森矩阵H(F1)-H(F2n)皆为4n+L阶的0矩阵;
节点功率方程的海森矩阵H(F2n+1)-H(F4n)为稀疏矩阵,节点i的有功功率方程海森矩阵H(F2n+2i-1)为4n+L阶矩阵,其中矩阵中值为1,其他元素皆为0,无功功率方程海森矩阵稀疏度相同,i最高取值为n;
VSC电压关系式方程的海森矩阵H(F4n+1)-H(F4n+L)皆为4n+L阶矩阵,其中为矩阵中的三个非零元素,其他皆为0;Mj为第j台换流站调制比,ecj,fcj为第j台换流站交流侧电压实虚部,j最高取值为L;
直流配电网中:
节点电压方程海森矩阵H(F1)-H(Fr)为2r阶0矩阵;
2.根据权利要求1所述的一种基于区间泰勒展开法的交直流混合配电网区间潮流计算方法,其特征在于,所述步骤S1的具体内容为:
交直流混合配电网的控制策略主要有主从控制策略与对等控制策略,其中,主从控制策略下,一台VSC作为主站,采用定直流电压控制方式,稳定直流电压值,其他VSC采用定功率控制方式,控制有功功率传输值;对等控制策略下,所有换流站皆采用下垂控制方式,同时调节直流配电网节点电压与节点注入电流有功功率值;
对于单台VSC,控制方式包括定直流电压、定交流电压、定有功或无功功率控制以及下垂控制方式;
定直流电压控制方式:
控制节点电压值为定值,区间潮流计算中功率波动引起该节点注入电流发生波动,因此节点注入功率为区间值;
上式中为直流配电网节点注入功率,为节点电压,为节点注入电流,为节点注入功率波动量,为节点注入电流波动量;无其他功率注入时,直流配电网注入VSC的功率为功率波动过程中,当波动量幅值较高时,功率传输方向会发生变化,发生整流与逆变状态的转换;直流配电网节点注入功率的波动引起VSC内部节点电压幅值与相角的波动,进而引起交流配电网内各状态量的波动;
定功率控制方式:
定功率控制方式下,其他节点功率注入值发生波动时,会引起该节点电压与节点注入电流发生微小波动,二者乘积计算后的功率波动值低于收敛精度,因此,在允许误差范围内认为节点注入功率为定值;由于采用交直流解耦潮流算法,定有功功率传输值改为定Pc,Pc为VSC与直流配电网间的有功功率传输值;无功功率传输值Qs维持定值,波动量设为0;
下垂控制方式:
直流配电网中连接VSC的节点出现电压上升或下降时,VSC根据电压偏差值决定注入节点的功率值;
根据区间算法,节点注入功率为:
定交流电压控制方式:
增广直角坐标下定交流电压控制方式维持节点i电压幅值恒为1,采用如下控制方式:
ei=1 fi=0 (5)
上式中ei为节点电压实部,fi为节点电压虚部;定交流电压方程中不存在功率波动变量,且雅可比矩阵中仅是关于ei的修正系数,因此,区间潮流计算中能够保持节点电压幅值恒为1,即电压幅值Ui=1。
3.根据权利要求1所述的一种基于区间泰勒展开法的交直流混合配电网区间潮流计算方法,其特征在于,所述步骤S2的具体内容为:
结合区间泰勒展开法与交直流混合配电网的控制方式得到交直流混合配电网中采用区间泰勒展开法的潮流模型,VSC两侧需要满足电压关系式如下
上式中,Mk为第k台VSC电压调制比,采用标幺值,基准值其中, ei,fi为第k台VSC内部节点电压的实部与虚部,每台VSC包含一个内部节点,为直流配电网中与第k台VSC相连的节点电压值;控制方式为定直流电压控制时,节点电压为定值;控制方式为下垂控制与定有功功率控制时,节点电压为区间值,直流配电网在基态下运行时,节点电压在标幺值1附近,因此
4.根据权利要求1所述的一种基于区间泰勒展开法的交直流混合配电网区间潮流计算方法,其特征在于,所述步骤S4的具体内容为:
交流配电网中考虑分布式发电波动性的不确定潮流方程中区间变量为功率值,其中节点i电源注入功率区间值为
考虑负荷功率波动性的不确定潮流方程中区间变量为功率值,其中节点i负荷功率区间值为
对节点有功,无功功率进行区间运算:
△Pi=△Pi L+△Pi S (53)
式(47)-(56)中,为节点i电源注入有功功率区间中心值,为节点i负荷有功功率区间中心值,△Pi S为节点i电源注入有功功率区间波动量,△Pi L为节点i负荷有功功率区间波动量,为节点i负荷无功功率区间中心值,为节点i电源注入无功功率区间中心值,为节点i负荷无功功率区间波动量,为节点i电源注入无功功率区间波动量;
直流配电网中考虑分布式发电波动性的不确定潮流方程中区间变量为功率值,其中节点i电源注入功率区间值为
考虑负荷功率波动性的不确定潮流方程中区间变量为功率值,其中节点i负荷功率区间值为
对节点有功功率进行区间运算
△Pi d=△Pi dL+△Pi dS (61)
5.根据权利要求1所述的一种基于区间泰勒展开法的交直流混合配电网区间潮流计算方法,其特征在于,所述步骤S5的具体内容为:
选取直流配电网中有功功率电源波动节点与负荷波动节点,计算节点有功功率波动幅值,节点注入功率波动幅度为±20%,在功率基值的基础上计算功率波动量;节点注入有功功率的基值为节点注入有功功率发生波动前的有功功率运行值;求得有功功率波动量后,对波动变量赋值,得到波动变量矩阵ε,此时波动变量矩阵ε仅包括直流配电网中的有功功率波动变量,为1×r阶;进行直流配电网的区间潮流计算,通过式(46)首先计算直流配电网的雅可比矩阵,进而根据式(34)求取状态量Vi d,对波动变量的一阶导数,通过式(30)计算直流配电网的海森矩阵,进而根据式(36)求取状态量Vi d,对波动变量的二阶导数,通过式(38)获取直流配电网节点电压,节点注入电流的区间值。
6.根据权利要求2所述的一种基于区间泰勒展开法的交直流混合配电网区间潮流计算方法,其特征在于,所述步骤S6的具体内容为:
步骤S5已获取直流配电网节点注入功率区间值,通过直流配电网区间潮流计算,得到了直流配电网中节点电压与节点注入电流的区间值;因此,通过式(3)获取直流配电网与VSC间有功功率传输区间值,根据式(9)计算得到有功功率传输波动量通过式(7),(8)获取直流配电网中与VSC相连接的节点电压波动量因此计及直流配电网输入交流配电网的有功功率值波动与直流配电网的电压波动状况对交流配电网的影响,选取交流配电网功率波动节点,节点有功,无功注入功率波动幅度皆为±20%,在功率基值Pic,Qic的基础上计算有功功率波动量εPi与无功功率波动量εQi;因此,直流配电网与VSC间的有功功率传输波动量直流配电网节点电压波动量交流配电网节点注入有功功率波动量εPi以及节点注入无功功率波动量εQi组成交流配电网的波动变量矩阵ε,为1×(2n+2L)阶矩阵;通过式(40)首先计算交流配电网的雅克比矩阵,进而根据式(34)求解交流配电网节点电压,节点注入电流以及电压调制比对波动变量 以及εQi的一阶导数;通过式(30)首先计算海森矩阵,进而根据式(36)求解交流配电网节点电压,节点注入电流以及电压调制比对波动变量εPi以及εQi的二阶导数,最后通过式(38)求解交流配电网节点电压,节点注入电流以及电压调制比的区间值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010673187.7A CN111799799B (zh) | 2020-07-13 | 2020-07-13 | 一种基于区间泰勒展开法的交直流混合配电网区间潮流计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010673187.7A CN111799799B (zh) | 2020-07-13 | 2020-07-13 | 一种基于区间泰勒展开法的交直流混合配电网区间潮流计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111799799A CN111799799A (zh) | 2020-10-20 |
CN111799799B true CN111799799B (zh) | 2022-03-08 |
Family
ID=72806811
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010673187.7A Active CN111799799B (zh) | 2020-07-13 | 2020-07-13 | 一种基于区间泰勒展开法的交直流混合配电网区间潮流计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111799799B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115588989A (zh) * | 2022-10-25 | 2023-01-10 | 国网安徽省电力有限公司合肥供电公司 | 基于交直流混合配电网交互模型的潮流计算方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107123994A (zh) * | 2017-04-28 | 2017-09-01 | 华南理工大学 | 区间无功优化模型的线性化求解方法 |
CN107204617A (zh) * | 2017-01-24 | 2017-09-26 | 华南理工大学 | 基于线性规划的直角坐标形式的区间潮流计算方法 |
CN107221935A (zh) * | 2017-06-12 | 2017-09-29 | 国网江西省电力公司经济技术研究院 | 基于置信度变换的交直流混联输电***区间潮流方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2341594A1 (en) * | 2009-12-29 | 2011-07-06 | Converteam Technology Ltd | Power collection and transmission systems |
CN107104442B (zh) * | 2017-04-20 | 2020-07-14 | 河海大学 | 计及参数模糊性的含风电场电力***概率潮流计算方法 |
CN109242365B (zh) * | 2018-11-06 | 2022-01-14 | 国电南瑞科技股份有限公司 | 一种电-热互联综合能源***的区间潮流计算方法 |
CN110445146B (zh) * | 2019-07-24 | 2023-06-23 | 广东电网有限责任公司电力调度控制中心 | 一种基于含多端柔直混联电网的潮流计算模型的潮流计算方法 |
CN110518591B (zh) * | 2019-08-22 | 2021-06-15 | 中国农业大学 | 一种不确定电力***的潮流计算方法 |
CN110932282B (zh) * | 2019-12-25 | 2021-06-22 | 福州大学 | 一种增广直角坐标下的基于vsc内部修正方程矩阵与交替迭代法的潮流计算方法 |
-
2020
- 2020-07-13 CN CN202010673187.7A patent/CN111799799B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107204617A (zh) * | 2017-01-24 | 2017-09-26 | 华南理工大学 | 基于线性规划的直角坐标形式的区间潮流计算方法 |
CN107123994A (zh) * | 2017-04-28 | 2017-09-01 | 华南理工大学 | 区间无功优化模型的线性化求解方法 |
CN107221935A (zh) * | 2017-06-12 | 2017-09-29 | 国网江西省电力公司经济技术研究院 | 基于置信度变换的交直流混联输电***区间潮流方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111799799A (zh) | 2020-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Roy et al. | Dynamic stability analysis of hybrid islanded DC microgrids using a nonlinear backstepping approach | |
CN105119280B (zh) | 基于锥优化的交直流混合结构有源配电网运行优化方法 | |
Hossain et al. | Robust control strategy for PV system integration in distribution systems | |
CN107171341B (zh) | 一种基于分布式计算的输配电网一体化无功优化方法 | |
CN105119275B (zh) | 一种计及统一潮流控制器的电力***动态最优潮流的算法 | |
CN112531790B (zh) | 一种虚拟电厂动态灵活性评估方法 | |
Mohammadi et al. | Fuzzy sliding-mode based control (FSMC) approach of hybrid micro-grid in power distribution systems | |
CN107947192A (zh) | 一种下垂控制型孤岛微电网的无功优化配置方法 | |
CN105162144B (zh) | 一种配电网电压无功优化的线性逼近求解方法 | |
CN106549392A (zh) | 一种配电网协调控制方法 | |
CN102684201A (zh) | 一种基于电压越限概率的含风电场电网无功优化方法 | |
CN114362267B (zh) | 考虑多目标优化的交直流混合配电网分散式协调优化方法 | |
CN110323979B (zh) | 一种考虑电压稳定的发电机励磁调差系数优化整定方法 | |
CN105305401B (zh) | 一种光伏高压直流并网变流器的并网控制方法 | |
CN104158198A (zh) | 配电网优化潮流控制装置和方法 | |
Ambia et al. | Power management of hybrid micro-grid system by a generic centralized supervisory control scheme | |
CN103368186A (zh) | 一种风电***无功优化方法 | |
CN108306334A (zh) | 基于粒子群优化算法的风电场内部无功优化策略 | |
Adebiyi et al. | Performance analysis of grid-tied photovoltaic system under varying weather condition and load | |
CN108376984A (zh) | 一种电力***节点电压自动调整方法 | |
CN111799799B (zh) | 一种基于区间泰勒展开法的交直流混合配电网区间潮流计算方法 | |
Bayati et al. | Accurate modeling of DC microgrid for fault and protection studies | |
CN108376997B (zh) | 一种考虑分布式电源不确定性的有源配电网孤岛划分方法 | |
CN106021754A (zh) | 考虑vsc无功越限调整策略的混联电网概率潮流算法 | |
Singh et al. | Particle swarm optimization approach for distributed generation allocation planning for voltage profile improvement |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |