本发明公开了一种含环网的三相不平衡的中低压主动配电网潮流计算方法。针对含环网的三相不平衡的中低压主动配电网，获得已知基础网络数据，通过虚拟电流补偿法，将环网解列为解环网，处理获得潮流迭代后的电网中各节点的潮流后电压，并且不断潮流处理迭代，每次潮流迭代后，计算潮流迭代前后电网各节点电压之差的最大值：并判断直至达到要求，实现完成了潮流处理。本发明方法能处理三相不平衡的中低压主动配电网中的环状不平衡网络，更快速地进行计算处理，占用内存小，潮流结果准确性好，对于三相不平衡的中低压主动配电网的构建能优化运行，向电网调度提供准确数据。
1.一种含环网的三相不平衡的中低压主动配电网潮流计算方法，所述的含环网的三相不平衡的中低压主动配电网包括三相不平衡的中压主动配电网、三相不平衡的低压主动配电网和中低压配电变压器的三个部分，中低压配电变压器连接三相不平衡的中压主动配电网和三相不平衡的低压主动配电网之间；主动配电网是指含分布式电源的配电网，环网是指由多条供电线路构成的“口”字型电网结构，环网处于三相不平衡的中压配电网，环网线路是指组成环网的多条供电线路中的距离源节点最远的一条供电线路；所述的中低压主动配电网中存在节点，节点包含用电负荷和分布式电源，各节点之间相连的线路为支路，以中低压主动配电网输入端所连接的上级供电变压器成为源节点，上级供电变压器与分布式电源共同向中低压主动配电网输送提供电能；其特征在于：方法包括以下几个步骤：
1)通过中低压主动配电网中的传感器采集获得以下中低压主动配电网的已知基础网络数据，包括：中低压主动配电网：总节点数n、线路长度Lij，其中i,j均表示节点的序数，i,j∈n；线路的单位阻抗Z、额定电压矩阵VN、中压配电网额定电压VN.M、低压配电网额定电压VN.L；
中压主动配电网节点的用电负荷的电能消耗有功功率PLoad.i.M、中压主动配电网节点的用电负荷的电能消耗无功功率QLoad.i.M、中压主动配电网节点的分布式电源的电能产生的有功功率PGen.i.M；
低压主动配电网节点的用电负荷的电能消耗有功功率PLoad.i.L、低压主动配电网节点的用电负荷的电能消耗无功功率QLoad.i.L、低压主动配电网节点的分布式电源的电能产生的有功功率PGen.i.L；
中低压配电变压器：第i个配电变压器的变比nt.i、第i个配电变压器的变比矩阵Ti、第i个配电变压器的阻抗Zt.i、第i个配电变压器的三相阻抗矩阵Zt.i、配电变压器电流电压修正矩阵AT；
先根据中压主动配电网节点的用电负荷的电能消耗有功功率PLoad.i.M、中压主动配电网节点的用电负荷的电能消耗无功功率QLoad.i.M、中压主动配电网节点的分布式电源的电能产生的有功功率PGen.i.M与中压主动配电网的各节点的额定电压VN.M计算得到中压主动配电网中各节点的综合流出电流Ii.M：其中，()*表示矩阵的共轭计算；j为虚数单位；
同时根据低压主动配电网节点的用电负荷的电能消耗有功功率PLoad.i.L、低压主动配电网节点的用电负荷的电能消耗无功功率QLoad.i.L、低压主动配电网节点的分布式电源的电能产生的有功功率PGen.i.L与低压主动配电网的各节点的额定电压VN.L计算得到低压主动配电网中各节点的综合流出电流Ii.L：最后根据中压配电网中各节点的流出电流Ii.M与低压配电网中各节点的流出电流Ii.L，共同获得中低压主动配电网中各节点的综合流出电流Ii：针对第二次迭代及以后的每次迭代时：
利用上次迭代后的各节点的电压Vi′以及主动配电网节点的用电负荷的电能消耗有功功率PLoad.i、主动配电网节点的用电负荷的电能消耗无功功率QLoad.i.、主动配电网节点的分布式电源的电能产生的有功功率PGen.i，计算获得中低压主动配电网中各节点的综合流出电流Ii：其中，Vi′表示上一次迭代后的电网各节点的电压；
3)将含环网的三相不平衡的中低压主动配电网中的环网线路去除获得解环网的三相不平衡的中低压主动配电网，由各节点的综合流出电流Ii处理获得解环网的节点的综合流出电流矩阵IH，并且获得解环网的三相不平衡的中低压主动配电网的节点支路关联矩阵A和支路电压关联矩阵C以及解环网的节点电压降落矩阵ΔVH和带环网的节点电压降落矩阵ΔV；
4)根据额定电压矩阵VN与带环网的节点电压降落矩阵ΔV相减得到迭代后的电网各节点的电压矩阵Vi′，完成本次潮流迭代：Vi′＝VN-ΔV
5)通过以下公式计算潮流迭代前后电网各节点电压之差的最大值ΔVi.max：ΔVi.max＝max(|Vi′-Vi|)
其中，Vi表示上一次迭代后节点i的电压，Vi′表示本次迭代后节点i的电压；
6)不断迭代重复步骤(2)～(5)进行潮流迭代，每次潮流迭代后，采用以下方式进行判断潮流迭代是否收敛；
若本次潮流迭代后收敛，则输出本次潮流迭代后的处理结果，获得本次潮流迭代后的电网中各节点的潮流后电压Vi′；
若本次潮流迭代后未收敛，将本次潮流迭代后的处理结果作为下次潮流迭代前的潮流前数值，将本次潮流迭代后的电网中各节点的潮流后电压Vi′作为下次潮流迭代时的电网中各节点的潮流前电压Vi，返回步骤2)进行下一次潮流迭代处理。
2.根据权利要求1所述的一种含环网的三相不平衡的中低压主动配电网潮流计算方法，其特征在于：所述1)中的低压主动配电网额定电压VN.L通过中压主动配电网额定电压VN.M与配电变压器的变比矩阵Ti处理计算得到：VN.M＝Ti·VN.L。
3.根据权利要求1所述的一种含环网的三相不平衡的中低压主动配电网潮流计算方法，其特征在于：所述1)中的额定电压矩阵VN采用以下方式处理获得：根据中压配电网额定电压VN.M与低压配电网额定电压VN.L共同组成了中低压主动配电网中各节点的额定电压VN：各个节点的额定电压VN组成的列矩阵作为额定电压矩阵VN。
4.根据权利要求1所述的一种含环网的三相不平衡的中低压主动配电网潮流计算方法，其特征在于：所述1)中，根据中压主动配电网节点的用电负荷的电能消耗有功功率PLoad.i.M和低压主动配电网节点的用电负荷的电能消耗有功功率PLoad.i.L和中压主动配电网节点的分布式电源的电能产生的有功功率PGen.i.M采用以下方式处理获得电网节点的用电负荷的电能消耗有功功率PLoad.i、电网节点的用电负荷的电能消耗无功功率QLoad.i和电网节点的分布式电源的电能产生的有功功率PGen.i：根据中压主动配电网节点的用电负荷的电能消耗有功功率PLoad.i.M与低压配电网节点的用电负荷的电能消耗有功功率PLoad.i.L共同组成了电网节点的用电负荷的电能消耗有功功率PLoad.i：根据中压主动配电网节点的用电负荷的电能消耗无功功率QLoad.i.M与低压主动配电网节点的用电负荷的电能消耗无功功率QLoad.i.L共同组成了电网节点的用电负荷的电能消耗无功功率QLoad.i：根据中压主动配电网节点的分布式电源的电能产生的有功功率PGen.i.M与低压主动配电网节点的分布式电源的电能产生的有功功率PGen.i.L共同组成了电网节点的分布式电源的电能产生的有功功率PGen.i：
5.根据权利要求1所述的一种含环网的三相不平衡的中低压主动配电网潮流计算方法，其特征在于：所述步骤3)具体如下：
3.1)由各节点的综合流出电流Ii组成节点的带环网的综合流出电流矩阵I，表示为：
3.2)由环网线路上流经的电流Bij.H拓展连接在节点的综合流出电流矩阵I之后组成解环网的节点的综合流出电流矩阵IH，表示为：
3.3)将含环网的三相不平衡的中低压主动配电网中的环网线路去除，得到一个解环网的三相不平衡的中低压主动配电网，并获得解环网的三相不平衡的中低压主动配电网的节点支路关联矩阵A，其中基于节点i的综合流出电流Ii对电网中节点i和节点j之间支路上的电流Bij的流经路径，节点支路关联矩阵A第i行j列元素取值如下：其中，Ti表示第i个配电变压器的变比矩阵，表示为：
3.4)通过以下公式，根据解环网的节点的综合流出电流矩阵IH与解环网的三相不平衡的中低压主动配电网的节点支路关联矩阵A获得解环网的电网支路的电流矩阵BH：其中，XA表示节点支路关联矩阵A的第一修正矩阵，YA表示节点支路关联矩阵A的第二修正矩阵；Bij为电网中节点i和节点j之间支路上的电流；B表示带环网的电网支路的电流矩阵，计算为：
3.5)采用以下方式生成解环网的三相不平衡的中低压主动配电网的支路电压关联矩阵C，其中基于电网中节点i和节点j之间支路上的电流Bij及配电变压器的位置，支路电压关联矩阵C中对应元素取值如下：其中，Zij表示线路阻抗矩阵，AT为配电变压器的电流电压修正矩阵，表示为：线路阻抗矩阵Zij根据线路单位阻抗Z、线路长度Lij和三阶单位矩阵1得到：Zij＝Z×Lij·1
Zt.i表示第i个配电变压器的三相阻抗矩阵，表示为：
3.6)通过以下公式，根据解环网的电网支路的电流矩阵BH与解环网的三相不平衡的中低压主动配电网的支路电压关联矩阵C获得解环网的节点电压降落矩阵ΔVH：其中，XC表示支路电压关联矩阵C的第一修正矩阵，YC表示支路电压关联矩阵C的第二修正矩阵，ZC表示支路电压关联矩阵C的第三修正矩阵；节点i的电压降落ΔVi为电网中节点i相对于额定电压VN的差；ΔV表示带环网的节点电压降落矩阵，表示为：
3.7)根据上述解环网的节点电压降落矩阵ΔVH的公式和解环网的电网支路的电流矩阵BH的公式，获得以下公式的解环网的节点电压降落矩阵ΔVH：接着通过库朗Kron简化处理得到带环网的节点电压降落矩阵ΔV为：
K＝C·A+XC·YA；L＝C·XA+XC；M＝YC·A+ZC·YA；N＝YC·XA+ZC。
6.根据权利要求1所述的一种含环网的三相不平衡的中低压主动配电网潮流计算方法，其特征在于：所述步骤6)中，潮流迭代是否收敛的判敛依据为：潮流迭代前后电网节点电压之差的最大值ΔVi.max是否小于等于10-5：ΔVi.max≤10-5
如果满足上述公式，则本次潮流迭代后收敛；否则本次潮流迭代后尚未收敛。
7.根据权利要求1所述的一种含环网的三相不平衡的中低压主动配电网潮流计算方法，其特征在于：所述步骤6)中，潮流迭代是否收敛的判敛依据为：潮流迭代计算次数k超过
100次，则判断该含环网的三相不平衡的中低压主动配电网潮流计算不收敛，结束含环网的三相不平衡的中低压主动配电网潮流计算。
含环网的三相不平衡的中低压主动配电网潮流计算方法
技术领域
[0001] 本发明属于配电网潮流处理领域，涉及含环网的三相不平衡的中低压主动配电网潮流计算方法。
背景技术
[0002] 潮流处理是对配电网进行管理与规划，并实现安全、稳定、经济运行的基础。近年来，随着能源电力、“云大物移智”等技术的不断进步与应用，分布式电源与泛在电力物联网发展结合度越来越高。随着风、光等分布式电源的快速大量接入，传统无源配电网正向有源主动配电网转变。在带来显著环境、经济和社会效益的同时，对配电网潮流计算也带来了新的挑战。首先，配电网由于线路不对称布置和转置缺失、负荷不均匀分布等因素影响本身具有显著的三相不平衡特征，而分布式电源(通常为单相)的随机接入进一步加剧了网络不平衡性；其次，分布式电源不仅直接影响其通常所接入的低压配网，同时间接影响其上游中压配网；另外，中压配电网由于供电可靠性的要求，其环网拓扑获得了更多重视和应用。
[0003] 截至目前，尚无潮流处理全面考虑了三相不平衡、多电压等级、分布式电源和弱环网等配电网新特征。
发明内容
[0004] 针对上述背景技术中的问题，本发明提出了一种含环网的三相不平衡的中低压主动配电网潮流计算方法，适用于多电压等级和环网的三相不平衡主动配电网综合直接潮流处理，应用该含环网的三相不平衡的中低压主动配电网潮流计算方法，保留了直接潮流的精准性与高效性，同时将其应用范围扩展至含多电压等级、环网的三相不平衡主动配电网，并可依据分布式电源的出力进行潮流计算，为新形势下配电网的潮流计算提供了有效解决方案。
[0005] 本发明的技术方案是：
[0006] 所述的含环网的三相不平衡的中低压主动配电网(以下文中部分地方简称电网)包括三相不平衡的中压主动配电网、三相不平衡的低压主动配电网和中低压配电变压器的三个部分，中低压配电变压器连接三相不平衡的中压主动配电网和三相不平衡的低压主动配电网之间；主动配电网是指含分布式电源的配电网。环网是指由多条供电线路构成的“口”字型电网结构，环网处于三相不平衡的中压配电网，环网线路是指组成环网的多条供电线路中的距离源节点最远的一条供电线路；所述的中低压主动配电网中存在节点，节点包含用电负荷和分布式电源，各节点之间相连的线路为支路，以中低压主动配电网输入端所连接的上级供电变压器成为源节点，上级供电变压器与分布式电源共同向中低压主动配电网输送提供电能；
[0007] 所述的节点为将汇集、分配和传送能量(电能)的设备。例如为住宅区用电中转站。
[0008] 方法包括以下几个步骤：
[0009] 1)通过中低压主动配电网中的传感器或者通过工具采集获得以下中低压主动配电网的已知基础网络数据，包括：
[0010] 中低压主动配电网：总节点数n、线路长度Lij，其中i,j均表示节点的序数，i,j∈n；线路的单位阻抗Z、额定电压矩阵VN、中压配电网额定电压VN.M、低压配电网额定电压VN.L；
[0011] 中压主动配电网节点的用电负荷的电能消耗有功功率PLoad.i.M、中压主动配电网节点的用电负荷的电能消耗无功功率QLoad.i.M、中压主动配电网节点的分布式电源的电能产生的有功功率PGen.i.M；
[0012] 低压主动配电网节点的用电负荷的电能消耗有功功率PLoad.i.L、低压主动配电网节点的用电负荷的电能消耗无功功率QLoad.i.L、低压主动配电网节点的分布式电源的电能产生的有功功率PGen.i.L；
[0013] 中低压配电变压器：第i个配电变压器的变比nt.i、第i个配电变压器的变比矩阵Ti、第i个配电变压器的阻抗Zt.i、第i个配电变压器的三相阻抗矩阵Zt.i、配电变压器电流电压修正矩阵AT；
[0014] 2)处理获得电网中各节点的综合流出电流Ii：
[0015] 针对第一次迭代时：
[0016] 先根据中压主动配电网节点的用电负荷的电能消耗有功功率PLoad.i.M、中压主动配电网节点的用电负荷的电能消耗无功功率QLoad.i.M、中压主动配电网节点的分布式电源的电能产生的有功功率PGen.i.M与中压主动配电网的各节点的额定电压VN.M计算得到中压主动配电网中各节点的综合流出电流Ii.M：
[0017]
[0018] 其中，()*表示矩阵的共轭计算；j为虚数单位。
[0019] 同时根据低压主动配电网节点的用电负荷的电能消耗有功功率PLoad.i.L、低压主动配电网节点的用电负荷的电能消耗无功功率QLoad.i.L、低压主动配电网节点的分布式电源的电能产生的有功功率PGen.i.L与低压主动配电网的各节点的额定电压VN.L计算得到低压主动配电网中各节点的综合流出电流Ii.L：
[0020]
[0021] 最后根据中压配电网中各节点的流出电流Ii.M与低压配电网中各节点的流出电流Ii.L，共同获得中低压主动配电网中各节点的综合流出电流Ii：
[0022]
[0023] 针对第二次迭代及以后的每次迭代时：
[0024] 利用上次迭代后的各节点的电压Vi′以及主动配电网节点的用电负荷的电能消耗有功功率PLoad.i、主动配电网节点的用电负荷的电能消耗无功功率QLoad.i.、主动配电网节点的分布式电源的电能产生的有功功率PGen.i，计算获得中低压主动配电网中各节点的综合流出电流Ii：
[0025]
[0026] 其中，Vi′表示上一次迭代后的电网各节点的电压，由步骤4)求得；
[0027] 3)将含环网的三相不平衡的中低压主动配电网中的环网线路去除获得解环网的三相不平衡的中低压主动配电网，由各节点的综合流出电流Ii处理获得解环网的节点的综合流出电流矩阵IH，并且获得解环网的三相不平衡的中低压主动配电网的节点支路关联矩阵A和支路电压关联矩阵C以及解环网的节点电压降落矩阵ΔVH和带环网的节点电压降落矩阵ΔV；
[0028] 4)根据额定电压矩阵VN与带环网的节点电压降落矩阵ΔV相减得到迭代后的电网各节点的电压矩阵Vi′，完成本次潮流迭代：
[0029] Vi′＝VN-△V
[0030] 5)通过以下公式计算潮流迭代前后电网各节点电压之差的最大值△Vi.max：
[0031] △Vi.max＝max(|Vi′-Vi|)
[0032] 其中，Vi表示上一次迭代后节点i的电压，Vi′表示本次迭代后节点i的电压；
[0033] 需要说明的是，对于第一次迭代，上一次迭代后的节点i的电压Vi为各节点的额定电压VN；第二次迭代后(包括第二次迭代)，上一次迭代后的节点i的电压Vi为上一次迭代中步骤5)所计算的中低压主动配电网各节点的电压。
[0034] 6)不断迭代重复步骤(2)～(5)进行潮流迭代，每次潮流迭代后，采用以下方式进行判断潮流迭代是否收敛；
[0035] 若本次潮流迭代后收敛，则输出本次潮流迭代后的处理结果，获得本次潮流迭代后的电网中各节点的潮流后电压Vi′；
[0036] 若本次潮流迭代后未收敛，将本次潮流迭代后的处理结果作为下次潮流迭代前的潮流前数值，将本次潮流迭代后的电网中各节点的潮流后电压Vi′作为下次潮流迭代时的电网中各节点的潮流前电压Vi，返回步骤2)进行下一次潮流迭代处理。
[0037] 本发明采用潮流计算的结果进行电网调度，优化电网中元件的优化配置，进行电网状态的监测。
[0038] 所述1)中的低压主动配电网额定电压VN.L通过中压主动配电网额定电压VN.M与配电变压器的变比矩阵Ti处理计算得到：
[0039] VN.M＝Ti·VN.L。
[0040] 所述1)中的额定电压矩阵VN采用以下方式处理获得：
[0041] 根据中压配电网额定电压VN.M与低压配电网额定电压VN.L共同组成了中低压主动配电网中各节点的额定电压VN：
[0042]
[0043] 各个节点的额定电压VN组成的列矩阵作为额定电压矩阵VN。
[0044] 所述1)中，根据中压主动配电网节点的用电负荷的电能消耗有功功率PLoad.i.M和低压主动配电网节点的用电负荷的电能消耗有功功率PLoad.i.L和中压主动配电网节点的分布式电源的电能产生的有功功率PGen.i.M采用以下方式处理获得电网节点的用电负荷的电能消耗有功功率PLoad.i、电网节点的用电负荷的电能消耗无功功率QLoad.i和电网节点的分布式电源的电能产生的有功功率PGen.i：
[0045] 根据中压主动配电网节点的用电负荷的电能消耗有功功率PLoad.i.M与低压配电网节点的用电负荷的电能消耗有功功率PLoad.i.L共同组成了电网节点的用电负荷的电能消耗有功功率PLoad.i：
[0046]
[0047] 根据中压主动配电网节点的用电负荷的电能消耗无功功率QLoad.i.M与低压主动配电网节点的用电负荷的电能消耗无功功率QLoad.i.L共同组成了电网节点的用电负荷的电能消耗无功功率QLoad.i：
[0048]
[0049] 根据中压主动配电网节点的分布式电源的电能产生的有功功率PGen.i.M与低压主动配电网节点的分布式电源的电能产生的有功功率PGen.i.L共同组成了电网节点的分布式电源的电能产生的有功功率PGen.i：
[0050]
[0051] 所述步骤3)具体如下：
[0052] 3.1)由各节点的综合流出电流Ii组成节点的带环网的综合流出电流矩阵I，表示为：
[0053]
[0054] 3.2)由环网线路上流经的电流Bij.H拓展连接在节点的综合流出电流矩阵I之后组成解环网的节点的综合流出电流矩阵IH，表示为：
[0055]
[0056] 3.3)将含环网的三相不平衡的中低压主动配电网中的环网线路去除/解除，得到一个解环网的三相不平衡的中低压主动配电网，并获得解环网的三相不平衡的中低压主动配电网的节点支路关联矩阵A，其中基于节点i的综合流出电流Ii对电网中节点i和节点j之间支路上的电流Bij的流经路径，节点支路关联矩阵A第i行j列元素取值如下：
[0057]
[0058] 其中，Ti表示第i个配电变压器的变比矩阵，表示为：
[0059]
[0060] 上述第i个配电变压器的变比矩阵Ti是根据变压器接线方式决定，采用Dyn11型接线。
[0061] 其中，1为三阶单位矩阵，0为三阶零矩阵；
[0062]
[0063]
[0064] 3.4)通过以下公式，根据解环网的节点的综合流出电流矩阵IH与解环网的三相不平衡的中低压主动配电网的节点支路关联矩阵A获得解环网的电网支路的电流矩阵BH：
[0065]
[0066] 其中，XA表示节点支路关联矩阵A的第一修正矩阵，YA表示节点支路关联矩阵A的第二修正矩阵；Bij为电网中节点i和节点j之间支路上的电流；B表示带环网的电网支路的电流矩阵，计算为：
[0067]
[0068] 所述的第一修正矩阵XA由n×1个维度为3×3的子矩阵组成的列矩阵，其中，根据环网线路所对应的连接节点：第i节点与第j节点，在修正矩阵XA中，第i个子矩阵为单位矩阵1，第j个子矩阵为-1，其余子矩阵为零矩阵0；第二修正矩阵YΑ由1×n个维度为3×3的零矩阵0组成的横矩阵。
[0069] 3.5)采用以下方式生成解环网的三相不平衡的中低压主动配电网的支路电压关联矩阵C，其中基于电网中节点i和节点j之间支路上的电流Bij及配电变压器的位置，支路电压关联矩阵C中对应元素取值如下：
[0070]
[0071] 其中，Zij表示线路阻抗矩阵，AT为配电变压器的电流电压修正矩阵，表示为：
[0072]
[0073] 上述配电变压器的电流电压修正矩阵AT是根据变压器接线方式决定，采用Dyn11型接线。
[0074] 线路阻抗矩阵Zij根据线路单位阻抗Z、线路长度Lij和三阶单位矩阵1得到：
[0075] Zij＝Z×Lij·1
[0076] 其中，Lij表示节点i和节点j之间的线路的长度；
[0077] Zt.i表示第i个配电变压器的三相变比矩阵，表示为：
[0078]
[0079] 3.6)通过以下公式，根据解环网的电网支路的电流矩阵BH与解环网的三相不平衡的中低压主动配电网的支路电压关联矩阵C获得解环网的节点电压降落矩阵ΔVH：
[0080]
[0081] 其中，XC表示支路电压关联矩阵C的第一修正矩阵，YC表示支路电压关联矩阵C的第二修正矩阵，ZC表示支路电压关联矩阵C的第三修正矩阵；节点i的电压降落△Vi为电网中节点i相对于额定电压VN的差；ΔV表示带环网的节点电压降落矩阵，表示为：
[0082]
[0083] 第一修正矩阵XC由n×1个维度为3×3的零矩阵0组成的列矩阵；第二修正矩阵YC由1×n个维度为3×3的子矩阵组成的横矩阵，其中，根据环网线路所对应的连接节点：第i节点与第j节点，在修正矩阵YC中，第i个子矩阵为支路电压关联矩阵C的第i个子矩阵，第j个子矩阵为负的支路电压关联矩阵C的第j个子矩阵，其余子矩阵为零矩阵0；第三修正矩阵ZC为环网线路的线路阻抗矩阵Zij.H。
[0084] 3.7)根据上述解环网的节点电压降落矩阵ΔVH的公式和解环网的电网支路的电流矩阵BH的公式，获得以下公式的解环网的节点电压降落矩阵ΔVH：
[0085]
[0086] 接着通过库朗(Kron)简化处理得到带环网的节点电压降落矩阵ΔV为：
[0087] ΔV＝(K-L·N-1·M)·I
[0088] K＝C·A+XC·YA；L＝C·XA+XC；M＝YC·A+ZC·YA；N＝YC·XA+ZC。
[0089] 其中矩阵L与矩阵M互为转置矩阵，即：L＝MT。
[0090] 所述步骤6)中，潮流迭代是否收敛的判敛依据为：潮流迭代前后电网节点电压之差的最大值△Vi.max是否小于等于10-5：
[0091] △Vi.max≤10-5
[0092] 如果满足上述公式，则本次潮流迭代后收敛；否则本次潮流迭代后尚未收敛。
[0093] 所述步骤6)中，潮流迭代是否收敛的判敛依据为：潮流迭代计算次数k超过100次，则判断该含环网的三相不平衡的中低压主动配电网潮流计算不收敛，结束含环网的三相不平衡的中低压主动配电网潮流计算。
[0094] 本发明的有益效果：
[0095] 本发明方法能处理三相不平衡的中低压主动配电网中的环状不平衡网络，并且能更快速地进行计算处理，占用内存小，潮流结果准确性好，对于含环网的三相不平衡的中低压主动配电网的构建能优化运行，向电网调度提供准确数据，提供准确的网络状态。
[0096] 本发明可以兼容包括光伏、风电、储能、微型燃气轮机、电动汽车等分布式电源的配电网络。
[0097] 本发明可以计算各种接线类型与变比的配电变压器网络。
[0098] 本发明可以准确计算配电变压器两侧三相不平衡电压的幅值与相角。
[0099] 本发明对大量节点的配电网仍具备较好的鲁棒性。
[0100] 本发明中低压主动配电网潮流计算提升了直接潮流的适用性，且在提升适用性的基础上依然保持了原电力系统直接潮流的高效性与高鲁棒性。
附图说明
[0101] 图1为本发明实施例的含环网的三相不平衡的中低压主动配电网示例图。
具体实施方式
[0102] 下面结合附图和实施例对本发明作进一步说明。
[0103] 按照本发明发明内容完整方法实施的实施例及其实施过程如下：
[0104] 1)通过电网中的传感器或者通过工具采集获得以下电网的已知基础网络数据，包括：
[0105] 电网：总节点数n、线路长度Lij，其中i,j均表示节点的序数，i,j∈n；线路单位阻抗Z、中压配电网额定电压VN.M；中压配电网节点的电能消耗有功功率PLoad.i.M、中压配电网节点的电能消耗无功功率QLoad.i.M、中压主动配电网节点的分布式电源的电能产生的有功功率PGen.i.M；低压配电网节点的电能消耗有功功率PLoad.i.L、低压配电网节点的电能消耗无功功率QLoad.i.L、低压主动配电网节点的分布式电源的电能产生的有功功率PGen.i.L；
[0106] 中低压配电变压器：第i个配电变压器的变比Zt.i；第i个配电变压器的三相变比矩阵Zt.i；配电变压器电流电压修正矩阵AT；
[0107] 具体实施例如如下电-热联合系统：
[0108] 如图1所示的含环网的三相不平衡的中低压主动配电网由六个节点组成。其中，存在六个节点，节点1代表上级供电变压器，为源节点；节点2、节点3和节点6均为中压配电网节点；节点3与节点6之间通过环网连接；节点3与节点4之间存在一个配电变压器，其连接方式为Dyn11；节点4与节点5为低压配电网节点；
[0109] 并且获得电网与热网的已知基础网络数据为：
[0110] 电力节点数：5；
[0111] 线路参数：线路单位长度电阻3.7Ω/km，线路单位长度电抗0.369Ω/km[0112] 线路长度：线路12：0.01km；线路23：0.01km；线路26：0.01km；线路45：0.01km[0113] 中压配电网电力节点额定电压：10kV；低压配电网电力节点额定电压：220V；
[0114] 节点用电设备的A相功率：节点2：有功功率500kW、无功功率0kVar；节点3：有功功率5000kW、无功功率0kVar；节点4：有功功率5000kW、无功功率0kVar；节点4：有功功率400kW、无功功率0kVar；节点5：有功功率400kW、无功功率0kVar。B相功率为A相的1.02倍，C相功率为A相的1.05倍。
[0115] 分布式电源：节点2：A、B、C相功率均为50kW；节点5：A、B、C相功率均为50kW；
[0116] 配电变压器变比：10000/220。
[0117] 1、根据步骤3.1)～3.4)，解环网的电网支路的电流矩阵BH为：
[0118]
[0119] 2、根据步骤3.5)～3.6)，解环网的节点电压降落矩阵ΔVH：
[0120]
[0121] 3、根据步骤3)至步骤5)，完成含环网的三相不平衡的中低压主动配电网潮流计算。
[0122] 实施例进行含环网的三相不平衡的中低压主动配电网潮流处理后，总迭代次数进行5次，计算时间为0.0027s。获得最后结果为：节点2潮流后电压依次为：A相9.97kV、B相9.96kV、C相9.94kV；节点3潮流后电压依次为：A相9.95kV、B相9.92kV、C相9.9kV；节点4潮流后电压依次为：A相218V、B相216V、C相213V；节点5潮流后电压依次为：A相215.6V、B相
211.8V、C相209V；节点6潮流后电压依次为：A相9.95kV、B相9.93kV、C相9.92kV。
[0123] 由此，可见本发明所提含环网的三相不平衡的中低压主动配电网潮流计算方法可以计算环网的三相不平衡的中低压主动配电网，同时，具有高效性与高鲁棒性。在配电变压器所连两侧网络均为三相不平衡的配电网时，本发明所提含环网的三相不平衡的中低压主动配电网潮流计算方法可以准确计算出两侧网络的不平衡的三相电压。在三相不平衡的中压主动配电网还有环网时，本发明所提含环网的三相不平衡的中低压主动配电网潮流计算方法可以计算出环网线路两端节点的不平衡的三相电压。在三相不平衡的中压主动配电网与三相不平衡的低压主动配电网均具有分布式电源时，可以准确计算出节点的不平衡的三相电压。
