CN106972539B - 一种基于锥规划的分布式电源就地电压控制策略整定方法 - Google Patents
一种基于锥规划的分布式电源就地电压控制策略整定方法 Download PDFInfo
- Publication number
- CN106972539B CN106972539B CN201710341630.9A CN201710341630A CN106972539B CN 106972539 B CN106972539 B CN 106972539B CN 201710341630 A CN201710341630 A CN 201710341630A CN 106972539 B CN106972539 B CN 106972539B
- Authority
- CN
- China
- Prior art keywords
- node
- distributed generation
- generation resource
- voltage
- active
- 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.)
- Expired - Fee Related
Links
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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/381—Dispersed generators
-
- 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
一种基于锥规划的分布式电源就地电压控制策略整定方法:根据选定的有源配电***,输入有源配电***结构及参数;依据有源配电***结构及参数,考虑分布式电源出力和负荷的时序特性,建立有源配电网分布式电源就地电压控制策略整定模型;将模型中的目标函数线性化,非线性约束进行线性化或锥转换,从而使有源配电网分布式电源就地电压控制策略整定模型转化为二阶锥模型;采用求解二阶锥规划的数学解算器对二阶锥模型进行计算求解;输出的求解结果,包括分布式电源就地电压控制策略的相关参数、分布式电源处电压分布情况、无功补偿情况和有功削减情况。本发明能够得到分布式电源就地电压控制策略,解决连续时间序列下的分布式电源就地电压控制策略的整定问题。
Description
技术领域
本发明涉及一种分布式电源的就地控制策略。特别是涉及一种基于锥规划的分布式电源就地电压控制策略整定方法。
背景技术
对能源和环境的高度关注使得配电网的发展面临着新的压力和挑战,这些压力和挑战同时也是推动传统配电网向有源配电网发展的重要机遇。近年来,包括光伏(Photovoltaic,PV)、风机等在内的分布式电源(Distributed Generation,DG)大量接入配电网。分布式电源渗透率的不断提高使有源配电网的双向潮流、无功与电压控制等问题更加突出,尤其是分布式电源的接入无论在时间还是空间上都具有很强的不确定性,在运行中往往导致馈线功率大幅波动、馈线间功率不平衡、以及电压越限等诸多问题。
目前,有源配电网主要采用集中式控制和就地控制两种不同的策略来实现***的运行控制。集中式控制策略统一调配光伏、风机等可控资源,利用全局信息进行全局优化,但分布式电源渗透率的提高使分布式电源的数据量显著增加,带来了沉重的通信和数据处理负担,增大了时间延迟,并且当中央控制***出现故障时,整个***都有失效的风险;另外,出于隐私以及安全方面的考虑,集中式控制难以获取全局的详细信息。而就地控制仅仅依靠本地测量信息,虽然无法实现全局最优,但不需要节点间的信息交流或远程量测,从而减少了通信的数据量,降低了控制变量的维度;并且,当分布式发电波动较大时,就地控制策略可以迅速响应,从而快速抑制波动。
分布式电源就地电压控制策略根据整定出的分布式电源无功控制策略和有功削减策略来实现对有源配电***的运行控制,改善电压水平,提高分布式电源消纳能力。若分布式电源用尽无功功率仍不能有效地将节点电压控制在理想区间时,加入有功削减策略可以解决这一问题,从而使***电压维持在一个理想的水平。
由于分布式电源的功率调节是连续变化的,其运行优化问题从单一的时间断面扩展为连续的时间序列,必须要以时序的有源配电网分布式电源就地电压控制策略整定模型作为优化问题的求解基础。该模型数学本质上是混合整数非线性规划问题,给计算的求解带来较大挑战。因此,需要一种能够快速求解上述混合整数非线性规划问题的模型与算法,用以求解有源配电网分布式电源就地电压控制策略整定模型,从而制定出分布式电源就地电压控制策略。
发明内容
本发明所要解决的技术问题是,提供一种满足***运行约束、建立有源配电网分布式电源就地电压控制策略整定模型、制定分布式电源就地电压控制策略的基于锥规划的分布式电源就地电压控制策略整定方法。
本发明所采用的技术方案是:一种基于锥规划的分布式电源就地电压控制策略整定方法,包括如下步骤:
1)根据选定的有源配电***,输入有源配电***结构及参数,包括:线路参数、负荷水平、网络拓扑连接关系,分布式电源的接入位置、容量及参数,分布式电源及负荷的日运行特性预测曲线,***安全运行电压约束和支路电流限制,***基准电压和基准功率的初值;
2)依据有源配电***结构及参数,考虑分布式电源出力和负荷的时序特性,建立有源配电网分布式电源就地电压控制策略整定模型,包括:选取根节点为平衡节点,设定有源配电***电压偏差之和、***损耗和有功削减量最小为目标函数,分别考虑***潮流约束、***安全运行约束和分布式电源运行约束;
3)将步骤2)所述模型中的目标函数线性化,非线性约束进行线性化或锥转换,从而使有源配电网分布式电源就地电压控制策略整定模型转化为二阶锥模型;
4)采用求解二阶锥规划的数学解算器对二阶锥模型进行计算求解;
5)输出步骤4)的求解结果,包括分布式电源就地电压控制策略的相关参数、分布式电源处电压分布情况、无功补偿情况和有功削减情况。
步骤2)中:
(1)所述的有源配电***电压偏差之和、***损耗和有功削减量最小为目标函数表示为
min f=αfV+βfL+γfP (1)
式中,α、β和γ分别为***电压偏差情况fV、***损耗fL和有功削减情况fP的权重系数,其中***电压偏差情况fV、***损耗fL和有功削减情况fP的表达式如下:
式中,NT和NN分别为时间断面数和***节点总数;Ωb为***支路集合;Vt,i为t时段节点i的电压幅值;Rij为支路ij的电阻,It,ij为t时段节点i流向节点j的电流幅值;为t时段节点i的分布式电源有功削减量;为Vt,i的期望电压区间,当Vt,i不在此区间时,目标函数fv用来减小电压偏差;Vi p为节点i的分布式电源有功削减策略的电压整定点,为Vi p的电压阈值,当Vi p没达到时,由于会出现过早进行有功削减的情况,因此此时目标函数fP用来减小有功削减;
(2)所述的***潮流约束表示为
式中,Ωb为支路的集合;Pt,ij为t时段节点i流向节点j的有功功率,Qt,ij为t时段节点i流向节点j的无功功率;Pt,ik为t时段节点i流向节点k的有功功率,Qt,ik为t时段节点i流向节点k的无功功率;Rij为支路ij的电阻,Xij为支路ij的电抗;It,ij为t时段节点i流向节点j的电流幅值;Vt,i为t时段节点i的电压幅值,Vt,j为t时段节点j的电压幅值;Pt,i为t时段节点i上注入的有功功率之和,分别为t时段节点i上分布式电源注入的有功功率、负荷消耗的有功功率,Qt,i为t时段节点i上注入的无功功率之和,分别为t时段节点i上分布式电源注入的无功功率、负荷消耗的无功功率;
(3)所述的***安全运行约束表示为
式中,Vmax和Vmin为***最大允许电压值和最小允许电压值;Vt,i为t时段节点i的电压幅值;It,ij为t时段节点i流向节点j的电流幅值;Imax为支路最大允许电流值;
(4)所述的分布式电源运行约束表示为
式中,为t时段节点i上分布式电源可提供的无功功率最大值;为节点i的分布式电源容量;为t时段节点i上分布式电源注入的有功功率;为t时段节点i上分布式电源注入的无功功率;为t时段节点i上分布式电源的有功削减量;为t时段节点i上分布式电源的有功功率预测值;和g(Vti)共同构成分布式电源就地电压控制策略的表达式;为分布式电源无功控制策略的表达式,存在调节死区[Vi q,min,Vi q,max],此时分布式电源产生的无功功率为0var;g(Vt,i)为节点i的分布式电源有功削减策略的表达式,g(Vt,i)存在调节死区[0.8,Vi p],此时分布式电源削减的有功功率为0W;下式分别构成和g(Vt,i):
步骤3)包括:
(1)采用U2,t,i和I2,t,ij替换二次项和将目标函数fV、***潮流约束和***运行约束线性化:
(Vmin)2≤U2,t,i≤(Vmax)2 (24)
I2,t,ij≤(Imax)2 (25)
式中,Ωb为支路的集合;Pt,ij为t时段节点i流向节点j的有功功率,Qt,ij为t时段节点i流向节点j的无功功率;Pt,ik为t时段节点i流向节点k的有功功率,Qt,ik为t时段节点i流向节点k的无功功率;Rij为支路ij的电阻,Xij为支路ij的电抗;Pt,i为t时段节点i上注入的有功功率之和,Qt,i为t时段节点i上注入的无功功率之和,分别为t时段节点i上分布式电源注入的无功功率、负荷消耗的无功功率;Vmax和Vmin为***最大允许电压值和最小允许电压值;Imax为支路最大允许电流值;NT和NN分别为时间断面数和***节点总数;Vt,i为t时段节点i的电压幅值;为Vt,i的最小期望电压值,为Vt,i的最大期望电压值;
(2)目标函数fV中含有绝对值项|U2,t,i-1|,引入辅助变量At,i,并增加约束:
At,i≥0 (30);
(3)目标函数fP中含有绝对值项由g(Vt,i)表达式得到,与Vi p成反比,将目标函数fP等价变换为目标函数fD:
目标函数fD中含有绝对值项|1.1-Vi p|,引入辅助变量Bi,并增加约束:
Bi≥0 (34)
式中,Vi p为节点i的分布式电源有功削减策略的电压整定点,为Vi p的电压阈值;
(4)分布式电源就地电压控制策略的表达式和g(Vt,i)为非线性表达式,采用分段线性化实现对和g(Vt,i)的精确线性化;通过引入辅助变量at,i,n,q n=1,2,…,6、dt,i,n,q n=1,2,…,5和at,i,n,p n=1,2,3,4、dt,i,n,p n=1,2,3,采用线段来近似和g(Vt,i)所定义的曲线,如下所示:
Vt,i=0.8at,i,1,q+0.9at,i,2,q+at,i,3,qVi q,min+at,i,4,qVi q,max+1.1at,i,5,q+1.2at,i,6,q(35)
at,i,1,q≤dt,i,1,q,at,i,6,q≤dt,i,5,q (37)
at,i,n,q≤dt,i,n,q+dt,i,n-1,q,n=2,3,4,5 (38)
at,i,n,q≥0,dt,i,n,q∈{0,1} (39)
式中,at,i,n,q n=1,2,…,6为连续变量,dt,i,n,q n=1,2,…,5为整数变量;
Vt,i=0.8at,i,1,p+at,i,2,pVi p+1.1at,i,3,p+1.2at,i,4,p (41)
g(Vt,i)=at,i,3,p+at,i,4,p (42)
at,i,1,p≤dt,i,1,p,at,i,4,p≤dt,i,3,p (43)
at,i,n,p≤dt,i,n,p+dt,i,n-1,p,n=2,3 (44)
at,i,n,p≥0,dt,i,n,p∈{0,1} (45)
式中,at,i,n,p n=1,2,3,4为连续变量,dt,i,n,p n=1,2,3为整数变量;
引入辅助整数变量ci,1,q、ci,2,q和ci,p分别将非线性乘积项at,i,3,qVi q,min、at,i,4,qVi q ,max和at,i,2,pVi p线性化,则at,i,3,qVi q,min、at,i,4,qVi q,max和at,i,2,pVi p可表示为:
at,i,3,qVi q,min=0.90at,i,3,q+0.01at,i,3,qci,1,q,0≤ci,1,q≤20 (47)
at,i,4,qVi q,max=0.90at,i,4,q+0.01at,i,4,qci,2,q,0≤ci,2,q≤20 (48)
ci,1,q≤ci,2,q (49)
at,i,2,pVip=0.90at,i,2,p+0.01at,i,2,pci,p,0≤ci,p≤20 (50)
ci,2,q≤ci,p (51)
其中,at,i,3,qci,1,q、at,i,4,qci,2,q和at,i,2,pci,p依旧是非线性乘积项,因此引入二进制变量li,1,m,q、li,2,m,q和li,m,p m=0,1,…,4表示at,i,3,qci,1,q、at,i,4,qci,2,q和at,i,2,pci,p;
引入辅助变量wt,i,1,m,q=at,i,3,qli,1,m,q、wt,i,2,m,q=at,i,4,qli,2,m,q和wt,i,m,p=at,i,2, pli,m,p,并取M为足够大的正实数,并增加以下约束:
at,i,3,q-(1-li,1,m,q)M≤wt,i,1,m,q≤at,i,3,q (55)
0≤wt,i,1,m,q≤li,1,m,qM (56)
at,i,4,q-(1-li,2,m,q)M≤wt,i,2,m,q≤at,i,4,q (57)
0≤wt,i,2,m,q≤li,2,m,qM (58)
at,i,2,p-(1-li,m,p)M≤wt,i,m,p≤at,i,2,p (59)
0≤wt,i,m,p≤li,m,pM (60)
(5)将约束条件式进行线性化,采用U2,t,i和I2,t,ij替换二次项和
再进一步进行凸松弛,得到二阶锥约束式:
本发明的一种基于锥规划的分布式电源就地电压控制策略整定方法,能够解决连续时间序列下的分布式电源就地电压控制策略的整定问题,充分考虑分布式电源和负荷的随机性和波动性,建立有源配电网分布式电源就地电压控制策略整定模型,采用二阶锥规划方法进行求解,得到分布式电源就地电压控制策略。
附图说明
图1是改进的PG&E69节点算例结构图;
图2是本发明的基于锥规划的分布式电源就地电压控制策略整定方法流程图;
图3是分布式电源及负荷运行特性预测曲线;
图4a是优化后得到的节点15处分布式电源就地电压控制策略;
图4b是优化后得到的节点26处分布式电源就地电压控制策略;
图4c是优化后得到的节点46处分布式电源就地电压控制策略;
图4d是优化后得到的节点54处分布式电源就地电压控制策略;
图5a是方案II和方案III下节点15处的分布式电源无功补偿情况;
图5b是方案II和方案III下节点26处的分布式电源无功补偿情况;
图5c是方案II和方案III下节点46处的分布式电源无功补偿情况;
图5d是方案II和方案III下节点54处的分布式电源无功补偿情况;
图6是方案II和方案III下节点26处的分布式电源有功削减情况;
图7a是优化前后节点15处的电压分布情况;
图7b是优化前后节点26处的电压分布情况;
图7c是优化前后节点46处的电压分布情况;
图7d是优化前后节点54处的电压分布情况。
具体实施方式
下面结合实施例和附图对本发明的一种基于锥规划的分布式电源就地电压控制策略整定方法做出详细说明。
如图2所示,本发明的一种基于锥规划的分布式电源就地电压控制策略整定方法,包括如下步骤:
1)根据选定的有源配电***,输入有源配电***结构及参数,包括:线路参数、负荷水平、网络拓扑连接关系,分布式电源的接入位置、容量及参数,分布式电源及负荷的日运行特性预测曲线,***安全运行电压约束和支路电流限制,***基准电压和基准功率的初值;
2)依据有源配电***结构及参数,考虑分布式电源出力和负荷的时序特性,建立有源配电网分布式电源就地电压控制策略整定模型,包括:选取根节点为平衡节点,设定有源配电***电压偏差之和、***损耗和有功削减量最小为目标函数,分别考虑***潮流约束、***安全运行约束和分布式电源运行约束;具体如下:
(1)所述的有源配电***电压偏差之和、***损耗和有功削减量最小为目标函数表示为
min f=αfV+βfL+γfP (1)
式中,α、β和γ分别为***电压偏差情况fv、***损耗fL和有功削减情况fP的权重系数,其中***电压偏差情况fV、***损耗fL和有功削减情况fP的表达式如下:
式中,NT和NN分别为时间断面数和***节点总数;Ωb为***支路集合;Vt,i为t时段节点i的电压幅值;Rij为支路ij的电阻,It,ij为t时段节点i流向节点j的电流幅值;为t时段节点i的分布式电源有功削减量;为Vt,i的期望电压区间,当Vt,i不在此区间时,目标函数fV用来减小电压偏差;Vi p为节点i的分布式电源有功削减策略的电压整定点,为Vi p的电压阈值,当Vi p没达到时,由于会出现过早进行有功削减的情况,因此此时目标函数fP用来减小有功削减;
(2)所述的***潮流约束表示为
式中,Ωb为支路的集合;Pt,ij为t时段节点i流向节点j的有功功率,Qt,ij为t时段节点i流向节点j的无功功率;Pt,ik为t时段节点i流向节点k的有功功率,Qt,ik为t时段节点i流向节点k的无功功率;Rij为支路ij的电阻,Xij为支路ij的电抗;It,ij为t时段节点i流向节点j的电流幅值;Vt,i为t时段节点i的电压幅值,Vt,j为t时段节点j的电压幅值;Pt,i为t时段节点i上注入的有功功率之和,分别为t时段节点i上分布式电源注入的有功功率、负荷消耗的有功功率,Qt,i为t时段节点i上注入的无功功率之和,分别为t时段节点i上分布式电源注入的无功功率、负荷消耗的无功功率;
(3)所述的***安全运行约束表示为
式中,Vmax和Vmin为***最大允许电压值和最小允许电压值;Vt,i为t时段节点i的电压幅值;It,ij为t时段节点i流向节点j的电流幅值;Imax为支路最大允许电流值;
(4)所述的分布式电源运行约束表示为
式中,为t时段节点i上分布式电源可提供的无功功率最大值;为节点i的分布式电源容量;为t时段节点i上分布式电源注入的有功功率;为t时段节点i上分布式电源注入的无功功率;为t时段节点i上分布式电源的有功削减量;为t时段节点i上分布式电源的有功功率预测值;和g(Vt,i)共同构成分布式电源就地电压控制策略的表达式;为分布式电源无功控制策略的表达式,存在调节死区此时分布式电源产生的无功功率为0var;g(Vt,i)为节点i的分布式电源有功削减策略的表达式,g(Vt,i)存在调节死区[0.8,Vi p],此时分布式电源削减的有功功率为0W;下式分别构成和g(Vt,i):
3)将步骤2)所述模型中的目标函数线性化,非线性约束进行线性化或锥转换,从而使有源配电网分布式电源就地电压控制策略整定模型转化为二阶锥模型;包括:
(1)采用U2,t,i和I2,t,ij替换二次项和将目标函数fv、***潮流约束和***运行约束线性化:
(Vmin)2≤U2,t,i≤(Vmax)2 (24)
I2,t,ij≤(Imax)2 (25)
式中,Ωb为支路的集合;Pt,ij为t时段节点i流向节点j的有功功率,Qt,ij为t时段节点i流向节点j的无功功率;Pt,ik为t时段节点i流向节点k的有功功率,Qt,ik为t时段节点i流向节点k的无功功率;Rij为支路ij的电阻,Xij为支路ij的电抗;Pt,i为t时段节点i上注入的有功功率之和,Qt,i为t时段节点i上注入的无功功率之和,分别为t时段节点i上分布式电源注入的无功功率、负荷消耗的无功功率;Vmax和Vmin为***最大允许电压值和最小允许电压值;Imax为支路最大允许电流值;NT和NN分别为时间断面数和***节点总数;Vt,i为t时段节点i的电压幅值;为Vt,i的最小期望电压值,为Vt,i的最大期望电压值;
(2)目标函数fV中含有绝对值项|U2,t,i-1|,引入辅助变量At,i,并增加约束:
At,i≥0 (30);
(3)目标函数fP中含有绝对值项由g(Vt,i)表达式得到,与Vi p成反比,将目标函数fP等价变换为目标函数fD:
目标函数fD中含有绝对值项|1.1-Vi p|,引入辅助变量Bi,并增加约束:
Bi≥0 (34)
式中,Vi p为节点i的分布式电源有功削减策略的电压整定点,为Vi p的电压阈值;
(4)分布式电源就地电压控制策略的表达式和g(Vt,i)为非线性表达式,采用分段线性化实现对和g(Vt,i)的精确线性化;通过引入辅助变量at,i,n,q=1,2,…,6、dt,i,n,qn=1,2,…,5和at,i,n,p n=1,2,3,4、dt,i,n,p n=1,2,3,采用线段来近似和g(Vt,i)所定义的曲线,如下所示:
Vt,i=0.8at,i,1,q+0.9at,i,2,q+at,i,3,qVi q,min+at,i,4,qVi q,max+1.1at,i,5,q+1.2at,i,6,q(35)
at,i,1,q≤dt,i,1,q,at,i,6,q≤dt,i,5,q (37)
at,i,n,q≤dt,i,n,q+dt,i,n-1,q,n=2,3,4,5 (38)
at,i,n,q≥0,dt,i,n,i∈{0,1} (39)
式中,at,i,n,q n=1,2,…,6为连续变量,dt,i,n,qn=1,2,…,5为整数变量;
Vt,i=0.8at,i,1,p+at,i,2,pVi p+1.1at,i,3,p+1.2at,i,4,p (41)
g(Vt,i)=at,i,3,p+at,i,4,p (42)
at,i,1,p≤dt,i,1,p,at,i,4,p≤dt,i,3,p (43)
at,i,n,p≤dt,i,n,p+dt,i,n-1,p,n=2,3 (44)
at,i,n,p≥0,dt,i,n,p∈{0,1} (45)
式中,at,i,n,p n=1,2,3,4为连续变量,dt,i,n,pn=1,2,3为整数变量;
引入辅助整数变量ci,1,q、ci,2,q和ci,p分别将非线性乘积项at,i,3,qVi q,min、at,i,4,qVi q ,max和at,i,2,pVi p线性化,则at,i,3,qVi q,min、at,i,4,qVi q,max和at,i,2,pVi p可表示为:
at,i,3,qVi q,min=0.90at,i,3,q+0.01at,i,3,qci,1,q,0≤ci,1,q≤20 (47)
at,i,4,qVi q,max=0.90at,i,4,q+0.01at,i,4,qci,2,q,0≤ci,2,q≤20 (48)
ci,1,q≤ci,2,q (49)
at,i,2,pVi p=0.90at,i,2,p+0.01at,i,2,pci,p,0≤ci,p≤20 (50)
ci,2,q≤ci,p (51)
其中,at,i,3,qci,1,q、at,i,4,qci,2,q和at,i,2,pci,p依旧是非线性乘积项,因此引入二进制变量li,1,m,q、li,2,m,q和li,m,p m=0,1,…,4表示at,i,3,qci,1,q、at,i,4,qci,2,q和at,i,2,pci,p;
引入辅助变量wt,i,1,m,q=at,i,3,qli,1,m,q、wt,i,2,m,q=at,i,4,qli,2,m,q和wt,i,m,p=at,i,2, pli,m,p,并取M为足够大的正实数,并增加以下约束:
at,i,3,q-(1-li,1,m,q)M≤wt,i,1,m,q≤at,i,3,q (55)
0≤wt,i,1,m,q≤li,1,m,qM (56)
at,i,4,q-(1-li,2,m,q)M≤wt,i,2,m,q≤at,i,4,q (57)
0≤wt,i,2,m,q≤li,2,m,qM (58)
at,i,2,p-(1-li,m,p)M≤wt,i,m,p≤at,i,2,p (59)
0≤wt,i,m,p≤li,m,pM (60)
(5)将约束条件式进行线性化,采用U2,t,i和I2,t,ij替换二次项和
再进一步进行凸松弛,得到二阶锥约束式:
4)采用求解二阶锥规划的数学解算器对二阶锥模型进行计算求解;
5)输出步骤4)的求解结果,包括分布式电源就地电压控制策略的相关参数、分布式电源处电压分布情况、无功补偿情况和有功削减情况。
本发明的一种基于锥规划的分布式电源就地电压控制策略整定方法实现了有源配电网分布式电源就地电压控制策略的求解。
对于本发明的实例,首先输入改进的PG&E69节点***中线路元件的阻抗值、负荷元件的有功功率基准值和功率因数、网络拓扑连接关系,算例结构如图1所示,详细参数见表1和表2;节点15和节点46分别接入一组光伏***,容量均为0.8MVA;节点26和节点54分别接入一组风机,容量均为1MVA;分布式电源及负荷运行特性预测曲线如图3所示;各节点电压幅值(标幺值)的安全运行上下限分别为1.10和0.90;最后设置***的基准电压为12.66kV、基准功率为1MVA。节点电压期望运行区间为0.97p.u.-1.03p.u.,就地电压有功控制曲线的电压阈值为1.05p.u.,***电压偏差情况、***损耗和有功削减情况的权重系数分别取0.22、0.5和0.28。
分别采用三种方案进行对比分析,方案I不使用控制手段,方案II采用本发明的基于锥规划的分布式电源就地电压控制策略,方案III采用集中控制策略,仿真结果见表3。
执行优化计算的计算机硬件环境为Intel(R)Xeon(R)CPU E5-1620,主频为3.70GHz,内存为32GB;软件环境为Windows 10操作***。
利用预测数据可以优化出本发明的基于锥规划的分布式电源就地电压控制策略的相关参数,见图4a~图4d,然后分布式电源可以根据基于锥规划的分布式电源就地电压控制策略来实时调节其无功补偿量和有功削减量,见图5a~图5d和图6,从而有效地减小电压偏差,降低网络损耗,见表3和图7a~图7d。由图7a~图7d可以看出,当不使用控制手段时,光伏、风机等分布式电源的接入会导致剧烈的电压波动。采用分布式电源进行就地电压调节后,当节点电压较低时,分布式电源发出无功功率起到电压支撑的作用;当节点电压较高时,分布式电源通过吸收无功功率来降低节点电压,若用尽无功功率仍不能有效地将节点电压控制在理想区间时,分布式电源削减有功功率,从而使***电压维持在一个理想的水平。
表1 PG&E69节点算例负荷接入位置及功率
表2 PG&E69节点算例线路参数
表3不同控制策略下的仿真结果比较
控制策略 | 电压最小值/p.u. | 电压最大值/p.u. | 网损/kW |
I.不使用控制策略 | 0.9326 | 1.0620 | 1815.0 |
II.就地电压控制策略 | 0.9516 | 1.0511 | 1265.4 |
III.集中控制策略 | 0.9554 | 1.0405 | 1088.4 |
Claims (2)
1.一种基于锥规划的分布式电源就地电压控制策略整定方法,其特征在于,包括如下步骤:
1)根据选定的有源配电***,输入有源配电***结构及参数,包括:线路参数、负荷水平、网络拓扑连接关系,分布式电源的接入位置、容量及参数,分布式电源及负荷的日运行特性预测曲线,***安全运行电压约束和支路电流限制,***基准电压和基准功率的初值;
2)依据有源配电***结构及参数,考虑分布式电源出力和负荷的时序特性,建立有源配电网分布式电源就地电压控制策略整定模型,包括:选取根节点为平衡节点,设定有源配电***电压偏差之和、***损耗和有功削减量最小为目标函数,分别考虑***潮流约束、***安全运行约束和分布式电源运行约束;包括:
(1)所述的有源配电***电压偏差之和、***损耗和有功削减量最小为目标函数表示为
min f=αfV+βfL+γfP (1)
式中,α、β和γ分别为***电压偏差情况fV、***损耗fL和有功削减情况fP的权重系数,其中***电压偏差情况fV、***损耗fL和有功削减情况fP的表达式如下:
式中,NT和NN分别为时间断面数和***节点总数;Ωb为***支路集合;Vt,i为t时段节点i的电压幅值;Rij为支路ij的电阻,It,ij为t时段节点i流向节点j的电流幅值;为t时段节点i的分布式电源有功削减量;为Vt,i的期望电压区间,当Vt,i不在此区间时,目标函数fV用来减小电压偏差;为节点i的分布式电源有功削减策略的电压整定点,为的电压阈值,当没达到时,由于会出现过早进行有功削减的情况,因此此时目标函数fP用来减小有功削减;
(2)所述的***潮流约束表示为
式中,Ωb为支路的集合;Pt,ij为t时段节点i流向节点j的有功功率,Qt,ij为t时段节点i流向节点j的无功功率;Pt,ik为t时段节点i流向节点k的有功功率,Qt,ik为t时段节点i流向节点k的无功功率;Rij为支路ij的电阻,Xij为支路ij的电抗;It,ij为t时段节点i流向节点j的电流幅值;Vt,i为t时段节点i的电压幅值,Vt,j为t时段节点j的电压幅值;Pt,i为t时段节点i上注入的有功功率之和,分别为t时段节点i上分布式电源注入的有功功率、负荷消耗的有功功率,Qt,i为t时段节点i上注入的无功功率之和,分别为t时段节点i上分布式电源注入的无功功率、负荷消耗的无功功率;
(3)所述的***安全运行约束表示为
式中,Vmax和Vmin为***最大允许电压值和最小允许电压值;Vt,i为t时段节点i的电压幅值;It,ij为t时段节点i流向节点j的电流幅值;Imax为支路最大允许电流值;
(4)所述的分布式电源运行约束表示为
式中,为t时段节点i上分布式电源可提供的无功功率最大值;为节点i的分布式电源容量;为t时段节点i上分布式电源注入的有功功率;为t时段节点i上分布式电源注入的无功功率;为t时段节点i上分布式电源的有功削减量;为t时段节点i上分布式电源的有功功率预测值;和g(Vt,i)共同构成分布式电源就地电压控制策略的表达式;为分布式电源无功控制策略的表达式,存在调节死区此时分布式电源产生的无功功率为0var;g(Vt,i)为节点i的分布式电源有功削减策略的表达式,g(Vt,i)存在调节死区此时分布式电源削减的有功功率为0W;下式分别构成和g(Vt,i):
3)将步骤2)所述模型中的目标函数线性化,非线性约束进行线性化或锥转换,从而使有源配电网分布式电源就地电压控制策略整定模型转化为二阶锥模型;
4)采用求解二阶锥规划的数学解算器对二阶锥模型进行计算求解;
5)输出步骤4)的求解结果,包括分布式电源就地电压控制策略的相关参数、分布式电源处电压分布情况、无功补偿情况和有功削减情况。
2.根据权利要求1所述的一种基于锥规划的分布式电源就地电压控制策略整定方法,其特征在于,步骤3)包括:
(1)采用U2,t,i和I2,t,ij替换二次项和将目标函数fV、***潮流约束和***运行约束线性化:
(Vmin)2≤U2,t,i≤(Vmax)2 (24)
I2,t,ij≤(Imax)2 (25)
式中,Ωb为支路的集合;Pt,ij为t时段节点i流向节点j的有功功率,Qt,ij为t时段节点i流向节点j的无功功率;Pt,ik为t时段节点i流向节点k的有功功率,Qt,ik为t时段节点i流向节点k的无功功率;Rij为支路ij的电阻,Xij为支路ij的电抗;Pt,i为t时段节点i上注入的有功功率之和,Qt,i为t时段节点i上注入的无功功率之和,分别为t时段节点i上分布式电源注入的无功功率、负荷消耗的无功功率;Vmax和Vmin为***最大允许电压值和最小允许电压值;Imax为支路最大允许电流值;NT和NN分别为时间断面数和***节点总数;Vt,i为t时段节点i的电压幅值;为Vt,i的最小期望电压值,为Vt,i的最大期望电压值;
(2)目标函数fV中含有绝对值项|U2,t,i-1|,引入辅助变量At,i,并增加约束:
At,i≥0 (30);
(3)目标函数fP中含有绝对值项由g(Vt,i)表达式得到,与成反比,将目标函数fP等价变换为目标函数fD:
目标函数fD中含有绝对值项引入辅助变量Bi,并增加约束:
Bi≥0 (34)
式中,为节点i的分布式电源有功削减策略的电压整定点,为的电压阈值;
(4)分布式电源就地电压控制策略的表达式和g(Vt,i)为非线性表达式,采用分段线性化实现对和g(Vt,i)的精确线性化;通过引入辅助变量at,i,n,q n=1,2,…,6、dt,i,n,q n=1,2,…,5和at,i,n,p n=1,2,3,4、dt,i,n,p n=1,2,3,采用线段来近似和g(Vt,i)所定义的曲线,如下所示:
at,i,1,q≤dt,i,1,q,at,i,6,q≤dt,i,5,q 37)
at,i,n,q≤dt,i,n,q+dt,i,n-1,q,n=2,3,4,5 38)
at,i,n,q≥0,dt,i,n,q∈{0,1} (39)
式中,at,i,n,q n=1,2,…,6为连续变量,dt,i,n,q n=1,2,…,5为整数变量;
g(Vt,i)=at,i,3,p+at,i,4,p 42)
at,i,1,p≤dt,i,1,p,at,i,4,p≤dt,i,3,p 43)
at,i,n,p≤dt,i,n,p+dt,i,n-1,p,n=2,3 44)
at,i,n,p≥0,dt,i,n,p∈{0,1} (45)
式中,at,i,n,p n=1,2,3,4为连续变量,dt,i,n,p n=1,2,3为整数变量;
引入辅助整数变量ci,1,q、ci,2,q和ci,p分别将非线性乘积项和线性化,则和可表示为:
ci,1,q≤ci,2,q 49)
ci,2,q≤ci,p 51)
其中,at,i,3,qci,1,q、at,i,4,qci,2,q和at,i,2,pci,p依旧是非线性乘积项,因此引入二进制变量li,1,m,q、li,2,m,q和li,m,p m=0,1,…,4表示at,i,3,qci,1,q、at,i,4,qci,2,q和at,i,2,pci,p;
引入辅助变量wt,i,1,m,q=at,i,3,qli,1,m,q、wt,i,2,m,q=at,i,4,qli,2,m,q和wt,i,m,p=at,i,2, pli,m,p,并取M为足够大的正实数,并增加以下约束:
at,i,3,q-(1-li,1,m,q)M≤wt,i,1,m,q≤at,i,3,q (55)
0≤wt,i,1,m,q≤li,1,m,qM (56)
at,i,4,q-(1-li,2,m,q)M≤wt,i,2,m,q≤at,i,4,q (57)
0≤wt,i,2,m,q≤li,2,m,qM (58)
at,i,2,p-(1--li,m,p)M≤wt,i,m,p≤at,i,2,p (59)
0≤wt,i,m,p≤li,mm,pM (60)
(5)将约束条件式进行线性化,采用U2,t,i和I2,t,ij替换二次项和
再进一步进行凸松弛,得到二阶锥约束式:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710341630.9A CN106972539B (zh) | 2017-05-13 | 2017-05-13 | 一种基于锥规划的分布式电源就地电压控制策略整定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710341630.9A CN106972539B (zh) | 2017-05-13 | 2017-05-13 | 一种基于锥规划的分布式电源就地电压控制策略整定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106972539A CN106972539A (zh) | 2017-07-21 |
CN106972539B true CN106972539B (zh) | 2019-05-14 |
Family
ID=59332172
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710341630.9A Expired - Fee Related CN106972539B (zh) | 2017-05-13 | 2017-05-13 | 一种基于锥规划的分布式电源就地电压控制策略整定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106972539B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107591797A (zh) * | 2017-08-18 | 2018-01-16 | 天津大学 | 一种智能软开关的集中和就地联合控制策略整定方法 |
CN109617092A (zh) * | 2018-12-19 | 2019-04-12 | 中国电力科学研究院有限公司 | 一种交直流混合电网动态无功优化方法及*** |
CN109728603A (zh) * | 2018-12-26 | 2019-05-07 | 天津大学 | 有源配电网分布式电源分区就地电压控制策略整定方法 |
CN110289646B (zh) * | 2019-06-19 | 2022-12-20 | 国网天津市电力公司 | 基于元模型的智能软开关就地控制策略优化方法 |
CN113422371B (zh) * | 2021-07-22 | 2022-05-20 | 天津大学 | 一种基于图卷积神经网络的分布式电源就地电压控制方法 |
CN115411739A (zh) * | 2022-09-22 | 2022-11-29 | 广东电网有限责任公司 | 数据-物理融合驱动的柔性配电***自适应电压控制方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104537258A (zh) * | 2015-01-12 | 2015-04-22 | 天津大学 | 分布式储能参与有源配电网运行调节的锥优化建模方法 |
CN105740973A (zh) * | 2016-01-25 | 2016-07-06 | 天津大学 | 基于混合整数锥规划的智能配电网综合电压无功优化方法 |
CN106159974A (zh) * | 2016-08-02 | 2016-11-23 | 清华大学 | 一种输配协调的分布式无功电压优化方法 |
CN106655177A (zh) * | 2017-01-18 | 2017-05-10 | 中国南方电网有限责任公司电网技术研究中心 | 基于扩展二阶锥规划的分布式电源最大接入能力计算方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9093842B2 (en) * | 2012-08-16 | 2015-07-28 | Mitsubishi Electric Research Laboratories, Inc. | Method for globally optimizing power flows in electric networks |
-
2017
- 2017-05-13 CN CN201710341630.9A patent/CN106972539B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104537258A (zh) * | 2015-01-12 | 2015-04-22 | 天津大学 | 分布式储能参与有源配电网运行调节的锥优化建模方法 |
CN105740973A (zh) * | 2016-01-25 | 2016-07-06 | 天津大学 | 基于混合整数锥规划的智能配电网综合电压无功优化方法 |
CN106159974A (zh) * | 2016-08-02 | 2016-11-23 | 清华大学 | 一种输配协调的分布式无功电压优化方法 |
CN106655177A (zh) * | 2017-01-18 | 2017-05-10 | 中国南方电网有限责任公司电网技术研究中心 | 基于扩展二阶锥规划的分布式电源最大接入能力计算方法 |
Non-Patent Citations (1)
Title |
---|
基于混合整数二阶锥规划的三相有源配电网无功优化;刘一兵等;《电力***自动化》;20140810;第38卷(第15期);第58-64页 |
Also Published As
Publication number | Publication date |
---|---|
CN106972539A (zh) | 2017-07-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106972539B (zh) | 一种基于锥规划的分布式电源就地电压控制策略整定方法 | |
CN106786631B (zh) | 基于锥规划的分布式电源就地电压无功控制策略整定方法 | |
CN106887852B (zh) | 一种间歇式分布式电源就地电压无功控制策略整定方法 | |
US20180309294A1 (en) | Reactive power optimization system and method of power grid based on the double-fish-swarm algorithm | |
CN105119280B (zh) | 基于锥优化的交直流混合结构有源配电网运行优化方法 | |
CN109728603A (zh) | 有源配电网分布式电源分区就地电压控制策略整定方法 | |
CN105574238B (zh) | 一种用于在线稳定裕度评估的断面有功控制策略分析方法 | |
CN106329523B (zh) | 考虑不确定性的有源配电网智能软开关鲁棒优化建模方法 | |
Xin et al. | Cooperative control strategy for multiple photovoltaic generators in distribution networks | |
CN110504691A (zh) | 一种计及vsc控制方式的交直流配电网最优潮流计算方法 | |
CN107591797A (zh) | 一种智能软开关的集中和就地联合控制策略整定方法 | |
CN106655227B (zh) | 一种基于智能软开关的有源配电网馈线负载平衡方法 | |
CN106655177B (zh) | 基于扩展二阶锥规划的分布式电源最大接入能力计算方法 | |
CN109149586A (zh) | 面向智能软开关的有源配电网分区分散式电压控制方法 | |
CN105375513B (zh) | 一种基于实时在线等值的110千伏风电场自动电压控制方法 | |
CN108023364A (zh) | 基于凸差规划的配电网分布式电源最大接入能力计算方法 | |
CN109390953A (zh) | 含分布式电源和电动汽车的低压配电网无功电压协调控制方法和*** | |
Chen et al. | An improved coordination method of multi-terminal MMC-HVDC system suitable for wind farm clusters integration | |
CN108233408B (zh) | 一种mmc-mtdc***自适应下垂控制方法 | |
CN103715687B (zh) | 一种主动配电网支路功率拥塞实时控制方法 | |
Zhang et al. | Distributed secondary control for island microgrids with expected dynamic performance under communication delays | |
CN110365023A (zh) | 一种考虑大电网稳定控制电压限值的自动电压控制方法 | |
CN107451743B (zh) | 一种考虑源荷运行特性的有源配电网孤岛运行方法 | |
CN103701142A (zh) | 考虑离散控制变量的主动配电网实时无功-电压控制方法 | |
CN109659934A (zh) | 基于短路电流裕度的新能源并网规划优选方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190514 |