计及三相不平衡的有源配电网重构_邱革非.pdf

2023年第41卷第2期计及三相不平衡的有源配电网重构邱革非，杨浩宇，何超，骆钊，刘铠铭，何虹辉（昆明理工大学电力工程学院，昆明650031）0引言大量光伏和风力等分布式电源（DistributedGeneration，DG）接入使得传统配电网逐步转变为多能源多路径系统，提高了配电网供电的可靠性1-2，但其灵活的接入特点以及出力的不确定性与波动性，会对配电网可靠性产生不利影响3-4。低压配电网中三相不平衡问题普遍存在5-8，其主要原因是负荷分配不均以及单相DG的接入。配电 网 重 构（Distribution Network Reconfiguration，DNR）是指通过对分段开关和联络开关的断开与闭合来改变网络拓扑结构，从而减小网损、提高电能质量、均衡网络有功和无功以及提高对DG的接纳能力9-11，进而缓解三相不平衡问题的方法。国内外学者对DNR进行了广泛而深入的研究，文献12分析了重构过程中DG出力不确定性的解决方法和各种重构模型的优缺点。文献13提出一种可以改进现有直角坐标的线性潮流计算无法处理PV节点三相不平衡线性潮流的计算方法，此方法无须迭代可直接计算配电网辐射状、弱环状和含PV摘要：针对单相光伏的大量接入引起配电网运行的三相不平衡、线路有功损耗增大问题，提出以有功损耗和三相不平衡度为目标函数，构建基于三相电流注入法的配电网动态重构模型，采用大M法对节点电压方程进行松弛，避免断开线路的约束条件对模型的影响，并采用多边形内近似法将非线性约束线性化。通过修改的IEEE 33节点配电系统的分析计算，验证了算法的有效性，提高了算法的求解效率和求解精度。关键词：三相电流注入法；分布式电源；配电网重构；有功损耗；三相不平衡度文献标志码：A中图分类号：TM744文章编号：1008-6218（2023）02-0014-09doi:10.19929/ki.nmgdljs.2023.0020基金项目 云南省应用基础研究计划资助项目“多能源互补综合能源系统的耦合机理与市场运行优化机制研究”（2101AT070080）；云南电网有限责任公司科技项目“以电能替代为核心的综合能源系统指标体系与配置优化研究”（YNKJXM20190087）Abstract：Aiming at the proldem of threephase unbalance and increased active power loss of distribution network caused bylarge access of single phase photovoltaic,taking active power loss and three phase imbalance as objective functions,propose a dynamic distribution network reconstruction model based on threephase current injection method.The large Mmethod is used to relax the node voltage equation to avoid the influence of the constraint condition of disconnected lineson the model.The nonlinear constraint is linearized by polygonal inner approximation.Through the analysis and calculationof the modified IEEE 33 node distribution system,the effectiveness of the proposed algorithm is verified,and the efficiencyand accuracy of the proposed algorithm are improved.Key words：threephase current injection method;distributed generation;distribution network reconfiguration;active loss;threephase voltage imbalanceActive Distribution Network Reconfiguration Considering ThreePhase ImbalanceQIU Gefei,YANG Haoyu,HE Chao,LUO Zhao,LIU Kaiming,HE Honghui（College of Electric Power Engineering,Kunming University of Science and Technology,Kunming650031,China）内 蒙 古 电 力 技 术INNER MONGOLIA ELECTRIC POWER142023年第41卷第2期邱革非，等：计及三相不平衡的有源配电网重构节点3种情况的三相不平衡潮流。文献14针对含不确定性的三相不平衡运行配电网重构问题，以有功网损最小为优化目标函数，考虑DG和负荷注入功率的不确定性进行建模，对一天内的时段进行合理划分，采用基于三相Distflow潮流方程建立含不确定性预算的配电网两阶段鲁棒动态重构模型，但鲁棒性是考虑了模型最差情况下的安全运行，结果相对保守，经济性较差。文献15建立了考虑系统有功损耗、节点最低电压幅值及开关操作次数的配电网重构多目标优化模型，结合粒子群算法和Pareto支配方法来求解，即多目标粒子群算法。虽然粒子群算法实现简单、收敛速度快，但很容易陷入局部最优从而导致早熟收敛。文献16提出了一种含分布式电源的主动配电网重构策略，以最小化网损为目标函数，将改进的教与学优化算法应用于配电网供电路径的快速优化。该方法收敛速度快、寻优时间短、运行稳定，但没有考虑配电网三相不平衡运行情况。大M法是运筹优化建模中常用到的方法，目的是将非线性问题转化成线性问题，使得问题能够使用运筹经典模型来求解。文献17从数学优化的角度研究了含分布式电源配电网重构的二阶锥规划模型。首先建立以配电网有功损耗最小化为目标的优化模型，其次采用凸松弛的Distflow模型描述配电网的潮流方程，并基于生成树的配电网连通性约束，将Distflow中的节点电压用大M约束进行松弛，但没有考虑分布式电源的波动性和负荷的时变性对配电网的影响。本文在大M法松弛模型约束的前提下，提出以降低有功损耗和三相不平衡度为目标，基于三相电流注入法（ThreePhase Current Injection Emethod，TCIM）18的配电网重构方法。该方法使用单相节点导纳矩阵描述三相网络，充分考虑大量单相分布式电源的波动性和单相负荷的时变性对网络的影响；使用多边形内近似法将所提模型转为混合整数化线 性 规 划（Mixed IntegerLinearProgramming，MILP）；加入各相充放电功率可以单独控制19的三相储能设备（Energy Storage，ES），提高单相资源的时空调用能力；加入单相静止无功发生器（StaticVar Generator，SVG）补偿无功，缓解网络无功不足和光伏逆变器牺牲部分有功出力来提高无功出力而造成的有功不足，提高网络单相的调节能力，配合重构改变供电路径，降低网络有功损耗，提高三相平衡度。最后，通过改进IEEE 33节点配电系统的计算与分析，验证本文方法的有效性。1优化模型1.1目标函数配电网运行时，由于单相负荷和单相分布式电源的时空分布不均，三相不平衡会导致相电压偏移变大，线路损耗增加20，影响电网的运行及用户用电。因此，本文以全天有功损耗最小和三相不平衡度最小为目标函数，其中三相不平衡度可用式（1）表示：|sk,t=|Vsk,t-Vk,tVk,tVk,t=sVsk,t3，（1）式中：集合L1，L2，L3，U，V，W，s；Vk,tt时刻节点k三相电压平均值；Vsk,tt时刻节点k的s相电压幅值；sk,tt时刻节点k的s相三相不平衡度。对于整个系统，三相不平衡度可用式（2）表示：t=1NkNssk,t3，（2）式中：tt时刻网络三相不平衡度；N网络节点集合。网络的有功损耗Ploss,t计算见式（3）：Ploss,t=()k,i Es()Pssd,ki,t-Psrsv,ki,t，（3）式中：E网络的线路集合；Pssd,ki,tt时刻线路ki上s相的送端功率；Psrsv,ki,tt时刻线路ki上s相的受端功率，其中：Pssd,ki,t=Vsk,t()reIski,t()re+Vsk,t()imIski,t()im，（4）Psrsv,ki,t=Vsi,t()reIski,t()re+Vsi,t()imIski,t()im，（5）式中：Iski,t()ret时刻线路ki上s相电流的实部；Iski,t()imt时刻线路ki上s相电流的虚部；Vsk,t()ret时刻节点k上s相电压实部；Vsk,t()imt时刻节点k上s相电压虚部；Vsi,t()ret时刻节点i上s相电压实部；Vsi,t()imt时刻节点i上s相电压虚部。本文考虑有功损耗和三相不平衡度的目标函数（F）如式（6）所示：152023年第41卷第2期内 蒙 古 电 力 技 术minF=tT1Ploss,tP*loss,t+2t*t，（6）式中：Tt的集合；1、2有功损耗和三相不平衡度的权重系数，且1+2=1；P*loss,t未优化前t时刻有功损耗；*t未优化前t时刻网络三相不平衡度。1.2约束条件决策变量在满足目标函数的同时，还需要满足潮流约束、电流注入约束、线路连接状态约束、网络辐射状态运行约束、开关动作次数约束、电压和电流幅值约束、DG出力约束、ES和SVG运行约束。1.2.1潮流约束在三相网络中，各相之间存在耦合，由于相间耦合关系的非线性会造成模型非凸，给求解各相的物理量带来很大困难，本文使用相间互导纳来近似地表示三相之间的耦合关系，相间互导纳为常量，可将相间非线性耦合线性化，并基于电路理论分别求解三相网络中各相的物理量。相间互导纳如式（7）所示：Y=|GL1L1kk+jBL1L1kkGL1L2kk+jBL1L2kkGL1L3kk+jBL1L3kkGL2L1kk+jBL2L1kkGL2L2kk+jBL2L2kkGL2L3kk+jBL1L3kkGL1L1kk+jBL1L1kkGL3L2kk+jBL3L2kkGL1L3kk+jBL3L3kkGL1L1ik+jBL1L1ikGL2L1ik+jBL2L1ikGL3L1ik+jBL3L1ikGL2L1ik+jBL2L1ikGL2L2ik+jBL2L2ikGL2L3ik+jBL2L3ikGL3L1ik+jBL3L1ikGL3L2ik+jBL3L2ikGL3L3ik+jBL3L3ik|GL1L1ki+jBL1L1kiGL1L2ki+jBL1L2kiGL1L3ki+jBL1L3kiGL2L1ki+jBL1L1kiGL2L2ki+jBL2L2kiGL2L3ki+jBL2L3kiGL3L1ki+jBL3L1kiGL3L2ki+jBL3L2kiGL3L3ki+jBL3L3kiGL1L1ii+jBL1L1iiGL1L2ii+jBL1L2iiGL1L3ii+jBL1L3iiGL2L1ii+jBL2L1kiGL2L2ii+jBL2L2iiGL2L3ii+jBL2L3iiGL3L1ii+jBL3L1iiGL3L2ii+jBL3L2iiGL3L3ii+jBL3L3ii，（7）式中：Y线路kk、ki、ii和sq相之间构成的66导纳矩阵；Y=G+jB,q，s，G为导纳Y的实部，B为导纳Y的虚部。网络的节点电压方程如式（8）所示，对于三相不含中性线的潮流模型，可根据式（8）推导出余下潮流约束：I=VY，（8）Iscalc,k,t()re=iN()iq()GsqkiVqi,t()re-BsqkiVqi,t()im，（9）Iscalc,k,t()im=iN()iq()GsqkiVqi,t()im+BsqkiVqi,t()re，（10）Iski,t()im=qBsqki()Vqi,t()re-Vqk,t()re+Gsqki()Vqi,t()im-Vqk,t()im，（11）Iski,t()re=qGsqki()Vqi,t()re-Vqk,t()re-Bsqki()Vqi,t()im-Vqk,t()im，（12）式中：N（i）与节点k相关联的节点集合，规定k节点到i节点为正方向；I网络节点电流向量；V网络节点电压向量；Iscalc,k,t()ret时刻节点k上s相电流实部；Iscalc,k,t()imt时刻节点k上s相电流虚部；Vqi,t()ret时刻节点i上q相电压实部；Vqi,t()imt时刻节点i上q相电压虚部；Vqk,t()ret时刻节点k上q相电压实部；Vqk,t()imt时刻节点k上q相电压虚部；Gsqki线路ki的电导；Bsqki线路ki上sq相的电纳。1.2.2电流注入约束如式（13）所示，节点的注入电流大小由注入功率确定，各节点注入有功功率()Psinj,k,t和无功功率()Qsinj,k,t大小如式（14）和（15）所示，将式（14）和（15）带入式（13）即可得出三相不含中性线模型，模型中各节点各相电流注入方程如式（16）和（17）所示。Sinj=IinjV，（13）Psinj,k,t=Pspv,k,t-Psload,k,t+PsES,k,t，（14）Qsinj,k,t=Qspv,k,t-Qsload,k,t+Qssvg,k,t，（15）Psinj,k,t=Vsk,t()reIsinj,k,t()re+Vsk,t()imIsinj,k,t()im，（16）Qsinj,k,t=Vsk,t()imIsinj,k,t()re-Vsk,t()reIsinj,k,t()im，（17）式中：Sinj网络的总注入功率；Iinj网络的总注入电流；Qspv,k,tt时刻节点k上s相光伏的无功功率；Qsload,k,tt时刻节点k上s相负荷的无功功率；Qssvg,k,tt时刻节点k上s相SVG的无功功率；Pspv,k,tt时刻节点k上s相光伏的有功功率；Psload,k,tt时刻节点k上s相负荷的有功功率；PsES,k,tt时刻节点k上s相ES的有功功率；Isinj,k,t()ret时刻节点k上s相注入电流的实部；Isinj,k,t()imt时刻节点k上s相注入电流的虚部。根据欧姆定律，其他节点流入该节点的电流大小等于该节点自身注入电流大小，据此可以得到基于注入电流不平衡量的潮流方程见式（18）和（19）。Isk,t()re=Iscalc,k,t()re-Isinj,k,t()re=0，（18）Isk,t()im=Iscalc,k,t()im-Isinj,k,t()im=0，（19）式中：Isk,t()ret时刻节点k上s相注入电流实部的不平衡量；162023年第41卷第2期邱革非，等：计及三相不平衡的有源配电网重构Isk,t()imt时刻节点k上s相注入电流虚部的不平衡量。1.2.3线路连接状态约束Iscalc,k,t()reiN()iq-()1-ki,tM()re+GsqkiVqi,t()re-BsqkiVqi,t()im，（20）Iscalc,k,t()reiN()iq()1-ki,tM()re+GsqkiVqi,t()re-BsqkiVqi,t()im，（21）Iscalc,k,t()imiN()iq-()1-ki,tM()im+GsqkiVqi,t()im+BsqkiVqi,t()re，（22）Iscalc,k,t()imiN()iq()1-ki,tM()im+GsqkiVqi,t()im+BsqkiVqi,t()re，（23）式中：M()re、M()im正实数；ki,tt时刻线路ki的连接状态。ki,t=1表示在t时刻线路ki闭合，式（20）（23）分别等价于式（9）和式（10），此时线路两端电压幅值存在相关性；ki,t=0表示在t时刻线路断开，由于较大正数M()re和M()im的存在，导致式（20）（23）中的约束不起作用，线路两端电压幅值不存在相关性，但断开线路ki两端电压幅值与网络其他节点存在相关性，因此断开线路两端电压幅值由其他节点的约束关系决定。1.2.4网络辐射状态运行约束()k,i Eki,t=n-1，（24）iN()iki,t=1，（25）ki,t+ik,t=ki,t，（26）1i=0，（27）0ki,t1，（28）式中：ki,tt时刻节点k与节点i之间的关系变量。式（24）表明线路与节点之间须满足的数量关系，否则会出现孤立的节点；ki,t=1时，节点k是节点i的父节点，反之节点k不是节点i的父节点；式（25）和（26）表明线路ki连接时，有且只有一个父节点；式（27）表示平衡节点一定不是父节点。1.2.5开关动作次数约束tTkiE|ki,t-ki,t-1Ns，（29）式中：Ns最大动作开关次数。1.2.6电压幅值约束|Vsk,t=()Vsk,t()re2+()Vsk,t()im2Vk,t()minVsk,tVk,t()max，（30）式中：Vk,t()mint时刻节点k上电压幅值的下限；Vk,t()maxt时刻节点k上电压幅值的上限。1.2.7电流幅值约束|Iski,t()mag=()Iski,t()re2+()Iski,t()im2()Iski,t()re2+()Iski,t()im2ki,tIki()max，（31）式中：Iski,t()magt时刻线路ki上s相电流幅值大小；Iki()max线路ki流经电流幅值的上限。1.2.8DG出力约束在约束式（14）和（15）中，光伏的有功和无功分别用节点注入的有功和无功表示，并将配电网出口母线认为是平衡节点，其他节点为PQ节点。光伏逆变器发出的无功功率不是无限的，而是受到多方面因素制约，本文考虑以可控的光伏逆变器为模型建模，其约束见式（32）（34）。Qspv,k,t=Pspv,k,ttan,tantanmin，（32）()Qspv,k,t2+()Pspv,k,t2Spv,k()max，（33）Pspv,k,t()Qspv,k,t2+()Pspv,k,t2tanmin，（34）式中：Qspv,k,tt时刻节点k上s相光伏逆变器发出的无功；Spv,k()max节点k上光伏的最大容量；tan光伏逆变器的功率因数；tanmin光伏逆变器的最小功率因数。1.2.9ES运行约束|EsES,k,t+1=EsES,k,t+()+Ps,+ES,k,t-Ps,-ES,k,t0Ps,+ES,k,tP+ES,k()maxBs,+ES,k,t0Ps,-ES,k,tP-ES,k()maxBs,-ES,k,tBs,-ES,k,t+Bs,+ES,k,t1EES,k()minEsES,k,tEES,k()max，（35）式中：Ps,+ES,k,tt时刻节点k上s相ES的充电功率；Ps,-ES,k,tt时刻节点k上s相ES的放电功率；172023年第41卷第2期内 蒙 古 电 力 技 术Bs,+ES,k,tt时刻节点k上s相充电功率的01变量；Bs,-ES,k,tt时刻节点k上s相放电功率的01变量；EsES,k,tt时刻节点k上s相ES的电量；+充电效率；-放电效率；EES,k()minES电量上限；EES,k()maxES电量下限；P+ES,k()max节点k的ES充电功率上限；P-ES,k()max节点k的ES放电功率上限。1.2.10SVG运行约束Qsvg,t()minQssvg,tQsvg,t()max，（36）式中：Qsvg,t()minSVG补偿无功的上限；Qsvg,t()maxSVG补偿无功的下限。2模型处理与求解根据优化模型，式（30）、（31）、（33）、（34）都为混合整数非线性问题，尽管现有的求解器能对其求解，但难以保证其全局最优性。本文将上述非线性约束线性化处理，从而将模型转化成线性问题，以便快速求解，并保证其最优性。2.1圆形约束线性化式（30）电压约束的可行域为小圆之外、大圆之内，见图1；式（31）和（33）电流和功率约束的可行域为圆的内部，见图2，可采用多边形内近似法对圆形约束进行线性化21。为了考虑转化后的精度，采用圆内接正12边形来近似表示上述圆形约束，多边形内近似法系数见表1。式（30）、（31）和（33）可以分别线性化表示为：|rc0Vsk,t()re+rc1Vsk,t()im+rc2Vk,t()max0,c1,12rc0Vsk,t()re+rc1Vsk,t()im+rc2Vk,t()min0,c1,12，（37）rc0Iski,t()re+rc1Iski,t()im+ki,trc2Iki()max0,c1,12，（38）rc0Pspv,k,t+rc1Qspv,k,t+rc2Spv,k()max0,c1,12，（39）式中：rc0、rc1、rc2线性化约束对应的系数；c多边形的边数，且cN。2.2功率因数线性化式（34）为功率因数取值范围约束，将功率因数的 取 值 设 置 在 滞 后 0.85 与 超 前 0.85 之 间，即tanmin=0.85。当功率因数为定值时，可使用不等式的变化方法，将式（35）转化为DG输出的无功功率约束。图1电压约束的多边形内近似法示意图Fig.1 Schematic diagram of the polygon innerapproximation for voltage constraintsVimVre图2电流和功率约束的多边形内近似法示意图Fig.2 Schematic diagram of the polygon innerapproximation for current and power constraintsQ/IimP/Ire0图中：Vim、Vre电压的虚部和实部。图中：Q无功功率；P有功功率；Iim、Ire电流的虚部和实部。表1多边形内近似法系数Tab.1 Coefficient of the polygon inner approximation methodc123456789101112rc0110.267 8-0.267 9-1-1-1-1-0.267 90.267 911rc10.267 911110.267 9-0.267 9-1-1-1-1-0.267 9rc2-1-1.366-1-1-1.366-1-1-1.366-1-1-1.366-1182023年第41卷第2期邱革非，等：计及三相不平衡的有源配电网重构-Pspv,k,t1-tan2mintanminQspv,k,tPspv,k,t1-tan2mintanmin。（40）3算例分析采用改进IEEE 33节点配电测试系统验证方法的有效性，任何网络连接关系都可以用节点导纳矩阵表示。本算例使用Opendss软件建立改进IEEE33导纳矩阵，使用基于python平台的pyomo建模语言建模，采用IPOPT求解器进行求解。计算机配置英特尔Core-i5-10600KF，主频4.10GHz，内存16GB。改进IEEE 33节点测试系统见图3，选择容量为500 kVA 的 10 kV/400 V 配电变压器，电压上下限为5%，在节点8-21、9-15、12-22、18-33、25-29设置了联络开关，初始状态均为断开，并在节点4、8、12、16、20、30接入单相光伏，光伏的接入情况以及额定容量设置如表2所示。节点10接入三相ES，三相容量均为30 kVA，最大充放电功率为15 kW，节点12接入单相SVG，补偿区间为025 kvar。单次潮流的优化只可以决定某一个时刻优化结果，为了直观地对同一时间不同条件下的算例进行分析评估，算例中以1天24 h为优化时长，步长为60 min，对配电网运行进行优化重构。本文采用文献22中提供的100个住宅的负荷使用情况，所提供的负荷数据以恒阻抗、恒电流、恒功率表示，即时变的ZIP负荷模型。为了充分体现网络的三相不平衡情况，在各节点各相中随机分配负荷。3.1原始网络分析为了表现本文所提方案对配电网的积极影响，先对原始系统进行分析，原始系统1天24 h内三相电压变化以及光伏出力和负荷有功功率变化见图4和图5。通过图4和图5可以看出，00:0006:00，光伏出力和负荷消耗功率基本保持一致，系统各相电压波动很小。在09:0017:00时，光伏出力大于负荷功率需求，但由于光伏的容量以及并网位置不同，三相电压的波动情况也不同。线路L1相光伏并网位置靠近线路末端，因此L1相末端节点电压提升幅度较大；线路L2相光伏并网位置靠近线路前端，因此L2相前端节点电压提升幅度较大；线路L3相的光伏并网位置靠近线路中端，因此L3相前端和末端电压提升幅度都较大。19:0024:00，三相负荷功率需求大于光伏出力，各相电压均下降。光伏并网位置不同以及其出力的波动性和负荷的时变性，导致各节点出现电压波动大、电能质量差等情况。3.2配电网动态重构分析表3为考虑ES和SVG配电网的33节点动态重构结果，图6为重构后网络三相电压情况，与重构前的电压情况对比，电压质量得到明显改善，各相电压均衡，没有出现电压越限的情况。为了对比本文所提方法与其他算法结果的差异，选取了改进遗传算法23和二阶锥规划算法24与本文算法以及遍历算法结果对比，本文仅选取了波动较大的11:0012:00进行对比，优化结果见表4。根据对比结果可知，二阶锥规划算法和本文的算法可以获得正确的重构策略，即遍历算法得到的结果；而改进遗传算法虽然能够降低网络有功损耗，但是陷入了局部最优解。从计算的精度来看，改进遗传算法编码不能全面将优化问题的约束表示出来，因此需要考虑对不可行解采用阈值，进而增加了工作量和求解时间，而二阶锥规划算法属于二次规划范畴，求解效率还是低于本文模型。3.3ES和SVG运行分析ES不平衡输入、输出和SVG的无功补偿可以控制三相网络某相的电压波动，从而有效提高网络三相平衡度。本文选取了11:0012:00，此时网络不平衡度以及有功损耗对比情况见表5。由表5可以看出，ES和SVG对网络三相不平衡有一定的改善，同时也能够降低有功损耗，但由于ES和SVG的出力有限以及供电路径的限制，对网络图3IEEE 33节点网络图Fig.3 IEEE 33 node network diagram12345678 9 10 111213 14 15 16 17 182324 25262728 29 30 31 32 33PVSVGPVPVPVPVESPV19 20 21 22表2光伏接入节点额定容量Tab.2 Rated capacity of photovoltaic access nodekVA接入节点4812162030L1相3015L2相2025L3相2520192023年第41卷第2期内 蒙 古 电 力 技 术的调节能力有限；重构可以改变网络供电路径，从而减小网络有功损耗。因有功损耗主要取决于线路两端的电压差，各节点电压越接近平衡节点电压（电压额定值），则有功损耗越小，所以网络三相不平衡度也会随之减小，但联络开关的开、断是针对三相线路，不能仅断开三相中某一相，对单相调节能力有限。在考虑ES和SVG的配电网重构中，在改变供电路径的同时能够根据实际情况合理分配单相资源，进一步改善三相不平衡情况和减小有功损耗。4结束语分布式电源灵活的接入特点以及出力的不确定性和负荷的时变性给配电网带来诸多挑战，同时在实际运行中，配电网多为三相不平衡状态，为了解决三相不平衡问题，本文提出以有功最小和三相不平衡度最小为目标函数的三相不平衡配电网重构模型。结果表明，三相电流注入法使用单相节点导纳矩阵描述三相网络，充分考虑了大量单相分布式电源的波动性和单相负荷的时变性对网络的影响，因此在计及三相不平衡的有源配电网重构模型中，本文算法的求解精度要高于其他算法。本文使用大M法对断开线路约束的松弛，避免了断开线路的约束对模型的影响，使用圆内接正12边形来近似表示上述圆形约束，将非线性约束线性化，提高了算法的求解效率；同时配电网动态重构缓解了配电网的三相不平衡情况，减少了配电网的有功损耗，并随着SVG和ES的并网，重构对配电网的优化效果将得到进一步提高。图5光伏和负荷有功功率曲线Fig.5 PV output and load active power curve140120100806040200-2004812162024t/h有功功率/kW光伏负荷1.101.051.000.950.9048121620243020101.101.091.081.071.061.051.041.031.021.011.000.980.970.95U（p.u.）t/h节点（a）L1相1.101.051.000.950.9048121620243020101.101.091.081.071.061.051.041.031.021.011.000.980.970.95U（p.u.）t/h节点（b）L2相1.101.051.000.950.9048121620243020101.101.091.081.071.061.051.041.031.021.011.000.980.970.95U（p.u.）t/h节点（c）L3相图4重构前三相电压变化Fig.4 Change of threephase voltagebefore reconstruction表333节点动态重构结果Tab.3 Dynamic reconstruction results of 33 nodes时段00：00-04：0004：00-06：0006：00-18：0018：00-24：00断开线路9-10,8-21,12-22,25-29,18-337-8,9-10,12-22,25-29,18-337-8,9-10,14-15,25-29,32-337-8,9-10,14-15,27-28,32-33202023年第41卷第2期邱革非，等：计及三相不平衡的有源配电网重构参考文献：1 于达，张玮，王辉.布式电源并网控制研究综述J.齐鲁工业大学学报，2022，36（1）：5965.YU Da,ZHANG Wei,WANG Hui.Survey on gridconnectedcontrol of distributed generationJ.Journal of Shandong Instituteof Light Industry,2022,36(1):5965.2 肖永江，于永进，张桂林.基于改进乌燕鸥算法的分布式电源优化配置J.电力系统保护与控制，2022，50（3）：148155.XIAO Yongjiang,YU Yongjin,ZHANG Guilin.Optimal configurationof distributed power generation based on an improved sootyternoptimizationalgorithmJ.PowerSystemProtectionandControl,2022,50(3):148155.3 赵永生，王吉文，凌松，等.分布式电源接入对配电网供电可靠性影响分析J.电气自动化，2020，42（6）：3234，111.ZHAO Yongsheng,WANG Jiwen,LING Song,et al.Analysison the influence of distributed power supply accessupon powersupplyreliabilityofdistributionnetworksJ.ElectricalAutomation,2020,42(6):3234,111.4 吴斌，单来支，王庆，等.基于三相不平衡的分布式光伏接入配电网的优化配置研究J.电力系统保护与控制，2019，47（24）：108114.WU Bin,SHAN Laizhi,WANG Qing,et al.Research on optimalconfiguration of distributed photovoltaic access distribution networkbased on three phase unbalanceJ.Power System Protectionand Control,2019,47(24):108114.5 冯亚明.低压配电网三相不平衡问题及处理对策探究J.科技创新与应用，2021，11（25）：116118.6 高永键.低压配电网三相不平衡问题及对策J.通信电源技术，2018，35（10）：251252.GAO Yongjian.Analysis of Threephase Unbalance of LowvoltageDistributionNetworkandCountermeasuresJ.TelecomPowerTechnologies,2018,35(10):251252.7 杨超.三相不平衡对低压配电网的安全经济性评估D.济南：济南大学，2019.8 张铧文，马莹.低压配电网三相负荷不平衡优化问题分析J.现代工业经济和信息化，2020，10（9）：2223.ZHANG Huawen,MA Ying.Analysis of Three phase LoadImbalance Optimization in Lowvoltage Distribution NetworkJ.Modern Industrial Economy and Informationization,2020,10(9):2223.9 秦昕，李修华，要航，等.含有分布式电源配电网重构的研究现状及发展趋势J.电气开关，2015，53（4）：7478.QIN Xin,LI Xiuhua,YAO Hang,et al.The Status Quo ofthe Current Research on Distribution Network ReconstructionofDistributedPowerSupplyandItsTendencyofFutureDevelopmentJ.Electric Switcher,2015,53(4):7478.10 邓桂秀，江修波，蔡金锭.浅析配电网网络重构研究现状和发展J.机电技术，2013，36（4）：8487.11 刘路宁，彭春华，温泽之，等.基于配电网动态重构的分布式光伏表4算法对比结果Tab.4 Algorithm comparison results算法改进遗传算法二阶锥规划算法本文算法遍历算法断开支路7-8,9-10,12-13,25-29,32-33重构前网损/kWh23.623.623.623.6网损/kWh16.414.814.614.2计算时间/s10236245372354表5不同方案结果对比Tab.5 Comparison of results of different schemes运行情况原始网络ES和SVG重构ES、SVG和重构网络最大三相不平衡度/%3.012.322.421.36有功损耗/kWh22.620.119.214.6图6重构后三相电压变化Fig.6 Change of threephase voltage after reconstruction1.101.051.000.950.9048121620243020101.101.091.081.071.061.051.041.031.021.011.000.980.970.95t/h节点U（p.u.）（a）L1相1.101.051.000.950.9048121620243020101.101.091.081.071.061.051.041.031.021.011.000.980.970.95t/h节点（b）L2相U（p.u.）1.101.051.000.950.9048