CN109980651A - 一种环网型配电系统潮流计算方法 - Google Patents
一种环网型配电系统潮流计算方法 Download PDFInfo
- Publication number
- CN109980651A CN109980651A CN201910425652.2A CN201910425652A CN109980651A CN 109980651 A CN109980651 A CN 109980651A CN 201910425652 A CN201910425652 A CN 201910425652A CN 109980651 A CN109980651 A CN 109980651A
- Authority
- CN
- China
- Prior art keywords
- matrix
- node
- branch
- looped network
- row
- 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.)
- Granted
Links
- 238000009826 distribution Methods 0.000 title claims abstract description 124
- 238000004364 calculation method Methods 0.000 title claims abstract description 24
- 239000011159 matrix material Substances 0.000 claims abstract description 228
- 238000013461 design Methods 0.000 claims abstract description 8
- 239000008186 active pharmaceutical agent Substances 0.000 claims description 11
- 238000010276 construction Methods 0.000 claims description 6
- CUZMQPZYCDIHQL-VCTVXEGHSA-L calcium;(2s)-1-[(2s)-3-[(2r)-2-(cyclohexanecarbonylamino)propanoyl]sulfanyl-2-methylpropanoyl]pyrrolidine-2-carboxylate Chemical group [Ca+2].N([C@H](C)C(=O)SC[C@@H](C)C(=O)N1[C@@H](CCC1)C([O-])=O)C(=O)C1CCCCC1.N([C@H](C)C(=O)SC[C@@H](C)C(=O)N1[C@@H](CCC1)C([O-])=O)C(=O)C1CCCCC1 CUZMQPZYCDIHQL-VCTVXEGHSA-L 0.000 claims 1
- 230000005855 radiation Effects 0.000 claims 1
- 238000000034 method Methods 0.000 abstract description 19
- 230000008569 process Effects 0.000 abstract description 9
- 230000005611 electricity Effects 0.000 description 5
- 239000000243 solution Substances 0.000 description 3
- 240000002853 Nelumbo nucifera Species 0.000 description 2
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 2
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 238000007726 management method Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000000547 structure data Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000013024 troubleshooting Methods 0.000 description 1
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
- 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
技术领域
本发明涉及低压配电技术领域，尤其涉及一种环网型配电系统潮流计算方法。
背景技术
目前，配电网潮流计算是配电网系统分析、经济运行的重要基础。配电网的网络重构、故障处理、无功优化和状态估计、线损分析等都需要用到配电网潮流计算的结果。一种性能优良的潮流计算方法是配电系统管理的关键。随着电力部门对配电网管理的重视程度不断加深，专门针对配电网的潮流计算研究广泛开展起来。由于配电网是输电网和电能用户的中间环节，电压等级较输电网低，稳态运行时网络结构呈树形、多分支的辐射状结构，而在负荷转移或故障发生时会出现弱环网运行状态。线路R/X值较高，多数情况大于1；低压配电变压器多位于负荷中心，自配电变压器低压侧向多个方向供电，可将配电系统上级变电站出口母线视为无穷大电源，多个方向的供电可分开进行计算，低压配电网的潮流从配变配电变压器低压侧流向负荷。
国内外学者根据配电网的特点提出了各种配电网潮流算法，如牛顿法、改进PQ解耦法、回路阻抗法和前推回代法等。但牛顿法需要形成导纳矩阵，并且雅可比矩阵的对角优势不复存在，很难收敛；改进PQ解耦法对R/X值较大的线路引入补偿技术，这种算法复杂化，丧失了快速解耦原有计算量小、收敛可靠的优点；回路阻抗法需要复杂的节点和支路编号，比较耗时。比较而言，前推回代潮流算法该方法物理概念明晰，但由于，配电网络结构庞大，分支较多，计算过程中对网络结构数据的搜索，影响了潮流的计算速度，数据每次前推时都计算支路的功率损耗，涉及直接对数据矩阵的计算，占用空间比较大又比较费时，且只适用于网络呈辐射状的结构的配电网络。
发明内容
本发明的目的是提供一种环网型配电系统潮流计算方法，能够快速、准确的进行潮流计算，且计算的工作量大大降低，收敛可靠，耗时少。
本发明采用的技术方案为：
一种环网型配电系统潮流计算方法,具体包括以下几个步骤：
步骤一，根据环网型配电系统的特点，建立环网型配电系统电路模型，将配电变压器低压侧作为电势节点，等值为电压幅值和相角为恒定已知量并假定三相电压对称；馈线上所有负荷等效为末端集中的节点负荷，且为恒功率即PQ负荷，馈线支路采用集中参数模型，故环网型配电系统中只有一个电势节点，其余节点都是负荷节点；环网型配电系统中部分支路连接形成回路，根据图论知识，确定环网型配电系统中回路的连支，并定义为环网状支路，环网状配电系统中除环网支路外的其他支路定义为辐射状支路；按潮流流动的方向，设定各回路电流的绕行方向，环网状支路电流的参考方向取回路电流的绕行方向；环网型配电系统中每一条支路与系统中的两个节点相连，辐射状支路潮流从电势节点出发流向负荷节点，将辐射状支路潮流流出节点定义为始端节点，另一端的节点定义为末端节点，将环网状支路的参考电流流出节点定义为始端节点，另一端的节点定义为末端节点；
步骤二，对环网型配电系统的电路模型中拓扑结构的支路和节点进行编号，电势节点编号为0，负荷节点依次编号为1,3,4……n-1，n为环网型配电系统的节点数，包括一个电势节点，n-1个负荷节点；支路从与电势节点相连支路出发，先对辐射状支路进行编号，使辐射状支路编号和与其相连的末端节点编号相同，依次编号为1,2,3……,b，b为环网型配电系统中辐射状支路的支路数；再对环网状支路进行编号，依次继续编号为b+1,b+2,……,b+m，m为环网型配电系统中环网状支路的支路数；则环网型配电系统的支路数为(b+m)，包括b条辐射状支路，m条环网状支路；
并通过初始数据矩阵DSH来表述环网型配电系统电路模型的拓扑结构和具体参数，初始数据矩阵DSH为(b+m)行，5列的矩阵，其中第i行为：DSH(i,:)＝[x,NS(x),NR(x),Z(x),S(NR(；x))其]中，i＝1,2,......,b+m；x为第i条支路编号，NS(x)为第i条支路始端节点的编号，NR(x)为第i条支路末端节点的编号，Z(x)为第i条支路阻抗，S(NR(x))为节点NR(x)的节点负荷复功率；
步骤三，通过初始数据矩阵DSH来构造环网型配电系统的节点支路关联矩阵NBH，节点支路关联矩阵NBH为n行，(b+m)列的矩阵，第j行第i列元素NBH(j,i)为：
其中，j＝1,2,......,n，i＝1,2,......,b+m；
步骤四，根据节点支路关联矩阵NBH构造环网型配电系统电路模型的回路路径矩阵LH，回路路径矩阵LH为m行，(b+m)列的矩阵，第s行第i列元素LH(s,i)为：
其中，s＝1,2,......,m，i＝1,2,......,b+m；
根据回路路径矩阵LH和初始数据矩阵DSH，构造回路阻抗矩阵ZLH，回路阻抗矩阵ZLH为m行，m列的矩阵，第s1行第s2列元素ZLH(s1,s2)为：
当s1＝s2＝s时，ZLH(s,s)表示回路s的自阻抗，等于回路s中所有支路阻抗之和，当s1≠s2时，ZLH(s1,s2)表示回路s1与回路s2的互阻抗，等于回路s1和s2共有支路阻抗代数和，s1＝1,2,......,m；s2＝1,2,......,m；s＝1,2,......,m；
其中，回路s表示环网状支路(b+s)所在的回路；
步骤五，将环网型配电系统的环网状支路断开，将环网型配电系统电路模型等效为对应的辐射型电路模型，并将注入末端节点的负荷电流用等效电流的形式来替代；
用等效数据矩阵DS来表述等效后辐射型电路模型的拓扑结构和具体参数，则等效数据矩阵DS为b行，5列的矩阵，第i行元素DS(i,:)＝DSH(i,:)；构造等效后辐射型电路模型的等效节点支路关联矩阵NB，等效节点支路关联矩阵NB为n行，b列的矩阵，第j行第i列元素NB(j,i)＝NBH(j,i)；
其中，i＝1,2,......,b；j＝1,2,......,n，i＝1,2,......,b；
步骤六，根据等效节点支路关联矩阵NB和等效数据矩阵DS构造等效后辐射型电路模型从电势节点到负荷节点的潮流流经的路径矩阵P，路径矩阵P为(n-1)行，b列的矩阵，第j行第i列元素P(j,i)为：
其中，j＝1,2,......,n-1，i＝1,2,......,b；
等效后辐射型电路模型中从电势节点0到负荷节点j的潮流流经支路的判定包括以下几步：
1)遍历NB矩阵中第(j+1)行元素NB(j,:)若NB(j,y)＝1，可确定支路y末端节点为j；
2)遍历NB矩阵中第y列元素NB(:,y)，若NB(z,y)＝-1，可确定支路y始端节点为(z-1)；
3)遍历NB矩阵中第z行元素NB(z,:)，若NB(z,t)＝1，可确定支路t末端节点为(z-1)；
4)重复以上步骤2,3，直到支路始端节点为电势节点0；
则从电势节点0到负荷节点j的潮流经路径为电势节点0……支路t，节点(z-1)，支路y，负荷节点j；
步骤七，根据等效后辐射型电路模型的路径矩阵P和等效数据矩阵DS，构造等效后辐射型电路模型从电势节点到负荷节点的潮流流经的支路阻抗矩阵ZP和节点负荷的复功率矩阵LB；支路阻抗矩阵ZP为(n-1)行，b列的矩阵，第j行第i列元素ZP(j,i)为：ZP(j,i)＝P(j,i)*Z(i)；节点负荷的复功率矩阵SB为b行，(n-1)列的矩阵，第i行第j列元素SB(i,j)为：SB(i,j)＝P(i,j)*S(j)；
其中，Z(i)＝DS(i,4)为第i条支路的阻抗；S(j)＝DS(i,5)为j节点的节点负荷复功率，支路i末端节点为j；j＝1,2,......,n-1，i＝1,2,......,b；
步骤八，计算环网型配电系统的各支路潮流和各负荷节点电压，具体包括：
1)设定各个负荷节点的电压初始值依次为V1，V2，…Vj，…Vn-1,其中，Vj表示负荷节点j的电压初始值，构造节点电压矩阵V，节点电压矩阵V为1行，(n-1)列的矩阵，其中第j列元素：
V(1,j)＝Vj
其中，j＝1,2,......,n-1；
2)计算可得回路压降矩阵VL，回路压降矩阵VL为m行，1列的矩阵，其中第s行第1列元素为：
VL(s,1)＝V(1,y)-V(1,z)；由ZLH*IL＝VL计算可得回路电流矩阵IL，回路电流矩阵IL为m行，1列的矩阵，第s回路的回路电流IL(s+b)＝IL(s,1)；
其中，环网状支路(b+s)的始端节点为y，末端节点为z，s＝1,2,......,m；
3)构造节点电压共轭矩阵VC，节点电压共轭矩阵为1行，n-1列的矩阵，其中第j列元素：
VC(1,j)＝V*(1,j)
节点负荷复功率共轭矩阵为SC，节点负荷复功率共轭矩阵SC为b行，(n-1)列的矩阵，其中第i行第j列元素为：
SC(i,j)＝SB*(i,j)
计算可得节点电流矩阵LC，节点电流矩阵LC为b行，(n-1)列的矩阵，第i行第j列元素为：
LC(i,j)＝SC(i,j)/VC(1,j)；其中，i＝1,2,......,b；j＝1,......,n-1；
4)计算可得等效节点电流矩阵LCT，等效节点电流矩阵LCT为b行，(n-1)列的矩阵，其中第i行第j列元素为：
若LC(i,j)＝0，则LCT(i,j)＝0
若节点j不与环网状支路相连，则LCT(i,j)＝LC(i,j)
若节点j与环网状支路相连，则:
A)若节点j与环网状支路k相连，且为首端节点，则令LC(i,j)＝LC(i,j)+IL(k)
B)若节点j与环网状支路k相连，且为末端节点，则令LC(i,j)＝LC(i,j)-IL(k)
C)令k＝k+1，若k≤b+m回到步骤A)，否则，到步骤D)
D)LCT(i,j)＝LC(i,j)
其中，i＝1,2,......,b；j＝1,......,n-1，k＝b+1；
5)计算可得支路电流矩阵FC，支路电流矩阵FC为b行，1列的矩阵，其中第i行元素为：
其中，i＝1,2,......,b；
6)计算可得压降矩阵D，总压降矩阵M和电压矩阵T；压降矩阵D，总压降矩阵M和电压矩阵T为n-1行，b列的矩阵；
D(j,i)＝ZP(j,i)*FC(i,1)
T(j,i)＝V0-M(j,i)
计算可得新节点电压矩阵Vnew，新节点电压矩阵Vnew为1行，(n-1)列的矩阵，第j列元素：Vnew(1,j)＝min(T(j,1),T(j,2),...T(j,(n-1)))
其中，V0为辐射型配电系统低压侧变压器出口处电压测量值，i＝1,2,......,b；j＝1,......,n-1；l＝1,2,......,i；
7)确定收敛条件，选择系统的误差精度为δ，判定|Vnew(1,j)-V(1,j)|≤δ是否成立，如果不成立，则V(1,j)＝Vnew(1,j)并回到步骤2)；
否则，可得环网型配电系统中各负荷节点电压和各支路潮流：
第j节点电压为Uj＝V(1,j)
第i支路电流为Ili＝FC(i,1)
第i支路潮流为
其中，i＝1,2,......,b；j＝1,......,n-1；
第k支路电流为Ilk＝IL(k)
第k支路潮流为
其中，k＝b+1,b+2,......,b+m；j＝1,......,n-1。
本发明根据环网型配电系统的特点建立电路模型，对电路模型中的节点和支路进行编号，并利用矩阵形式来表示电路的初始数据，尤其是电路拓扑结构的数据化描述，使电路模型中图形连接的逻辑关系清晰明确，关联矩阵的引入实现了对路径搜索过程的程序化设计，并用初始数据矩阵作为唯一数据输入，使整个程序设计直接简单；计算过程中利用图论的知识，对环网型配电系统回路进行合理展开，等效为辐射型配电系统进行计算，大大简化了计算过程，并且计算过程中主要利用矩阵元素的代数运算求解配电网的各支路潮流和各负荷节点电压，因而程序计算用时较少，收敛速度较快。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明的流程图；
图2为本发明实施例所述环网型配电系统的电路模型示意图；
图3为本发明实施例所述环网型配电系统的电路模型支路节点编号示意图；
图4为本发明实施例所述等效辐射型配电系统的电路模型示意图；
图5为本发明实施例所述等效辐射型配电系统从电势节点到负荷节点9潮流流经支路的判定示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有付出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
如图1、2和3所示，本发明具体包括以下几个步骤：
步骤一，根据环网型配电系统的特点，建立环网型配电系统电路模型，将配电变压器低压侧作为电势节点，等值为电压幅值和相角为恒定已知量并假定三相电压对称；馈线上所有负荷等效为末端集中的节点负荷，且为恒功率即PQ负荷，馈线支路采用集中参数模型，故环网型配电系统中只有一个电势节点，其余节点都是负荷节点；环网型配电系统中部分支路连接形成回路，根据图论知识，确定环网型配电系统中回路的连支，并定义为环网状支路，环网状配电系统中除环网支路外的其他支路定义为辐射状支路；按潮流流动的方向，设定各回路电流的绕行方向，环网状支路电流的参考方向取回路电流的绕行方向；环网型配电系统中每一条支路与系统中的两个节点相连，辐射状支路潮流从电势节点出发流向负荷节点，将辐射状支路潮流流出节点定义为始端节点，另一端的节点定义为末端节点，将环网状支路的参考电流流出节点定义为始端节点，另一端的节点定义为末端节点；
步骤二，对环网型配电系统的电路模型中拓扑结构的支路和节点进行编号，电势节点编号为0，负荷节点依次编号为1,3,4……n-1，n为环网型配电系统的节点数，包括一个电势节点，n-1个负荷节点；支路从与电势节点相连支路出发，先对辐射状支路进行编号，使辐射状支路编号和与其相连的末端节点编号相同，依次编号为1,2,3……,b，b为环网型配电系统中辐射状支路的支路数；再对环网状支路进行编号，依次继续编号为b+1,b+2,……,b+m，m为环网型配电系统中环网状支路的支路数；则环网型配电系统的支路数为(b+m)，包括b条辐射状支路，m条环网状支路；
并通过初始数据矩阵DSH来表述环网型配电系统电路模型的拓扑结构和具体参数，初始数据矩阵DSH为(b+m)行，5列的矩阵，其中第i行为：DSH(i,:)＝[x,NS(x),NR(x),Z(x),S(NR(；x))其]中，i＝1,2,......,b+m；x为第i条支路编号，NS(x)为第i条支路始端节点的编号，NR(x)为第i条支路末端节点的编号，Z(x)为第i条支路阻抗，S(NR(x))为节点NR(x)的节点负荷复功率；
步骤三，通过初始数据矩阵DSH来构造环网型配电系统的节点支路关联矩阵NBH，节点支路关联矩阵NBH为n行，(b+m)列的矩阵，第j行第i列元素NBH(j,i)为：
其中，j＝1,2,......,n，i＝1,2,......,b+m；
步骤四，根据节点支路关联矩阵NBH构造环网型配电系统电路模型的回路路径矩阵LH，回路路径矩阵LH为m行，(b+m)列的矩阵，第s行第i列元素LH(s,i)为：
其中，s＝1,2,......,m，i＝1,2,......,b+m；
根据回路路径矩阵LH和初始数据矩阵DSH，构造回路阻抗矩阵ZLH，回路阻抗矩阵ZLH为m行，m列的矩阵，第s1行第s2列元素ZLH(s1,s2)为：
当s1＝s2＝s时，ZLH(s,s)表示回路s的自阻抗，等于回路s中所有支路阻抗之和，
当s1≠s2时，ZLH(s1,s2)表示回路s1与回路s2的互阻抗，等于回路s1和s2共有支路阻抗代数和，s1＝1,2,......,m；s2＝1,2,......,m；s＝1,2,......,m；
其中，回路s表示环网状支路(b+s)所在的回路；
步骤五，将环网型配电系统的环网状支路断开，将环网型配电系统电路模型等效为对应的辐射型电路模型，并将注入末端节点的负荷电流用等效电流的形式来替代；
用等效数据矩阵DS来表述等效后辐射型电路模型的拓扑结构和具体参数，则等效数据矩阵DS为b行，5列的矩阵，第i行元素DS(i,:)＝DSH(i,:)；构造等效后辐射型电路模型的等效节点支路关联矩阵NB，等效节点支路关联矩阵NB为n行，b列的矩阵，第j行第i列元素NB(j,i)＝NBH(j,i)；
其中，i＝1,2,......,b；j＝1,2,......,n，i＝1,2,......,b；
步骤六，根据等效节点支路关联矩阵NB和等效数据矩阵DS构造等效后辐射型电路模型从电势节点到负荷节点的潮流流经的路径矩阵P，路径矩阵P为(n-1)行，b列的矩阵，第j行第i列元素P(j,i)为：
其中，j＝1,2,......,n-1，i＝1,2,......,b；
等效后辐射型电路模型中从电势节点0到负荷节点j的潮流流经支路的判定包括以下几步：
5)遍历NB矩阵中第(j+1)行元素NB(j,:)若NB(j,y)＝1，可确定支路y末端节点为j；
6)遍历NB矩阵中第y列元素NB(:,y)，若NB(z,y)＝-1，可确定支路y始端节点为(z-1)；
7)遍历NB矩阵中第z行元素NB(z,:)，若NB(z,t)＝1，可确定支路t末端节点为(z-1)；
8)重复以上步骤2,3，直到支路始端节点为电势节点0；
则从电势节点0到负荷节点j的潮流经路径为电势节点0……支路t，节点(z-1)，支路y，负荷节点j；
步骤七，根据等效后辐射型电路模型的路径矩阵P和等效数据矩阵DS，构造等效后辐射型电路模型从电势节点到负荷节点的潮流流经的支路阻抗矩阵ZP和节点负荷的复功率矩阵LB；支路阻抗矩阵ZP为(n-1)行，b列的矩阵，第j行第i列元素ZP(j,i)为：ZP(j,i)＝P(j,i)*Z(i)；节点负荷的复功率矩阵SB为b行，(n-1)列的矩阵，第i行第j列元素SB(i,j)为：SB(i,j)＝P(i,j)*S(j)；
其中，Z(i)＝DS(i,4)为第i条支路的阻抗；S(j)＝DS(i,5)为j节点的节点负荷复功率，支路i末端节点为j；j＝1,2,......,n-1，i＝1,2,......,b；
步骤八，计算环网型配电系统的各支路潮流和各负荷节点电压，具体包括：
1)设定各个负荷节点的电压初始值依次为V1，V2，…Vj，…Vn-1,其中，Vj表示负荷节点j的电压初始值，构造节点电压矩阵V，节点电压矩阵V为1行，(n-1)列的矩阵，其中第j列元素：
V(1,j)＝Vj
其中，j＝1,2,......,n-1；
2)计算可得回路压降矩阵VL，回路压降矩阵VL为m行，1列的矩阵，其中第s行第1列元素为：
VL(s,1)＝V(1,y)-V(1,z)；由ZLH*IL＝VL计算可得回路电流矩阵IL，回路电流矩阵IL为m行，1列的矩阵，第s回路的回路电流IL(s+b)＝IL(s,1)；
其中，环网状支路(b+s)的始端节点为y，末端节点为z，s＝1,2,......,m；
3)构造节点电压共轭矩阵VC，节点电压共轭矩阵为1行，n-1列的矩阵，其中第j列元素：
VC(1,j)＝V*(1,j)
节点负荷复功率共轭矩阵为SC，节点负荷复功率共轭矩阵SC为b行，(n-1)列的矩阵，其中第i行第j列元素为：
SC(i,j)＝SB*(i,j)
计算可得节点电流矩阵LC，节点电流矩阵LC为b行，(n-1)列的矩阵，第i行第j列元素为：
LC(i,j)＝SC(i,j)/VC(1,j)；其中，i＝1,2,......,b；j＝1,......,n-1；
4)计算可得等效节点电流矩阵LCT，等效节点电流矩阵LCT为b行，(n-1)列的矩阵，其中第i行第j列元素为：
若LC(i,j)＝0，则LCT(i,j)＝0
若节点j不与环网状支路相连，则LCT(i,j)＝LC(i,j)
若节点j与环网状支路相连，则:
A)若节点j与环网状支路k相连，且为首端节点，则令LC(i,j)＝LC(i,j)+IL(k)
B)若节点j与环网状支路k相连，且为末端节点，则令LC(i,j)＝LC(i,j)-IL(k)
C)令k＝k+1，若k≤b+m回到步骤A)，否则，到步骤D)
D)LCT(i,j)＝LC(i,j)
其中，i＝1,2,......,b；j＝1,......,n-1，k＝b+1；
5)计算可得支路电流矩阵FC，支路电流矩阵FC为b行，1列的矩阵，其中第i行元素为：
其中，i＝1,2,......,b；
6)计算可得压降矩阵D，总压降矩阵M和电压矩阵T；压降矩阵D，总压降矩阵M和电压矩阵T为n-1行，b列的矩阵；
D(j,i)＝ZP(j,i)*FC(i,1)
T(j,i)＝V0-M(j,i)
计算可得新节点电压矩阵Vnew，新节点电压矩阵Vnew为1行，(n-1)列的矩阵，第j列元素：Vnew(1,j)＝min(T(j,1),T(j,2),...T(j,(n-1)))
其中，V0为辐射型配电系统低压侧变压器出口处电压测量值，i＝1,2,......,b；j＝1,......,n-1；l＝1,2,......,i；
7)确定收敛条件，选择系统的误差精度为δ，判定|Vnew(1,j)-V(1,j)|≤δ是否成立，如果不成立，则V(1,j)＝Vnew(1,j)并回到步骤2)；
否则，可得环网型配电系统中各节点负荷电压和各支路潮流：
第j节点负荷电压为Uj＝V(1,j)
第i支路电流为Ili＝FC(i,1)
第i支路潮流为
其中，i＝1,2,......,b；j＝1,......,n-1；
第k支路电流为Ilk＝IL(k)
第k支路潮流为
其中，k＝b+1,b+2,......,b+m；j＝1,......,n-1；
参见图1，为了更好的解释本发明的方法步骤，以下以一个具体实施例进行解释说明，具体包括以下几个步骤：
步骤一，参见图2，根据环网型配电系统的特点，建立电路模型，将配电变压器低压侧作为电源节点1，等值为电压幅值和相角为恒定已知量并假定三相电压对称；馈线上所有负荷等效为节点负荷2,故配电系统中只有一个电势节点1，其余节点都是负荷节点2，配电系统中环网电路部分6中支路连接形成回路，根据图论知识，每个回路仅含一条连支，且这个连支不出现在其他回路中，确定回路的连支，定义为环网状支路7，并设定环网状支路所在回路的回路电流参考方向8，标示环网状支路电流的参考方向8与回路的绕行方向一致，环网型配电系统中除环网状支路外的其他支路称为辐射状支路4；配电系统中每一条支路与系统中的两个节点相连，辐射状支路潮流从电势节点出发流向负荷节点，将辐射状支路潮流流出节点定义为始端节点3，另一端的节点定义为末端节点5，将环网状支路的参考电流流出节点定义为始端节点，另一端的节点定义为末端节点；
步骤二，参见图3，对环网型的电路模型拓扑结构中的支路和节点进行编号，其中，配电系统的节点数n为10，包括一个电势节点，9个负荷节点；配电系统的支路数(b+m)为(9+2)，包括9个辐射状支路和2个环网状支路。源节点编号为0，负荷节点依次编号为1,2,3……9；从与源节点相连支路出发，对辐射状支路进行编号，使支路编号与相连的末端节点编号相同，依次编号为[1],[2],[3]……[9]；再对形成回路的环网状支路的进行编号，依次继续编号为[10],[11]；
并通过初始数据矩阵DS来表述电路模型的拓扑关系和具体参数如下：
其中，初始数据矩阵DS为11行，5列的矩阵，Z(i)为第i条支路阻抗，S(j)为第i条支路末端节点为j的节点负荷复功率，i＝1,2,......,9，j＝1,2,......,11；
步骤三，通过初始数据矩阵DSH可构造配电系统的节点支路关联矩阵NBH如下：
步骤四，根据节点支路关联矩阵NBH构造环网型电路模型回路路径矩阵LH如下：
构造环网型电路模型回路阻抗矩阵ZLH如下
步骤五，参见图4，将配电系统的环网型支路断开，将环网状电路模型等效为对应的辐射型电路模型如图3所示，注入末端节点负荷电流用环网状支路等效电流的形式来替代；
则用等效数据矩阵DS来表述等效后辐射型电路模型的拓扑结构和具体参数，
等效节点支路关联矩阵NB为：
步骤六，参见图5，根据等效节点支路关联矩阵NB可得由电势节点到负荷节点的潮流流经的路径矩阵P如下：
步骤七，根据路径矩阵P和等效初始数据矩阵DS，构造从电势节点到负荷节点的潮流电流流经的支路阻抗矩阵ZP和节点的负荷复功率矩阵SB如下：
步骤八，计算配电网的各支路电流和各负荷节点电压，具体包括
1)设定其他各个负荷节点电压初始值依次为V1，V2…Vj…Vn-1,其中，Vj表示第j节点电压初始值，构造这里令所有负荷节点的电压都等于电势节点电压V0，则电压矩阵V如下：
V＝[V0 V0 V0 V0 V0 V0 V0 V0 V0]
2)计算可得回路压降矩阵VL如下：
由ZLH*IL＝VL可得
回路电流矩阵IL为m行，1列的矩阵，定义IL(s+b)＝IL(s,1)；即为第s回路的回路电流；
其中，s＝1,2,......,m；
3)进一步构造节点电压共轭矩阵VC，节点电压共轭矩阵为1行，9列的矩阵，其中第j列元素：
VC(1,j)＝V*(1,j)
节点负荷复功率共轭矩阵为SC，节点负荷复功率共轭矩阵SC为9行，9列的矩阵，其中第i行第j列元素为：
SC(i,j)＝SB*(i,j)
其中，i＝1,2,......,b；j＝1,......,n-1；
4)计算可得节点电流矩阵LC，支路电流矩阵FC，
LC(i,j)＝SC(i,j)/VC(1,j)
其中，i＝1,2,......,b；j＝1,......,n-1；
为书写方便，这里不妨设定：
5)计算可得等效电流矩阵LCT，
由可得
6)计算可得压降矩阵D，总压降矩阵M和电压矩阵T；压降矩阵D；
D(j,i)＝ZP(j,i)*FC(i,1)
T(j,i)＝V0-M(j,i)
计算可得新节点电压矩阵Vnew，新节点电压矩阵Vnew为1行，9列的矩阵，第j列元素：Vnew(1,j)＝min(T(j,1),T(j,2),...T(j,(n-1)))
其中，
M(1)＝Z(1)*Il(1)
M(2)＝Z(1)*Il(1)+Z(2)*Il(2)
M(3)＝Z(1)*Il(1)+Z(2)*Il(2)+Z(3)*Il(3)
M(4)＝Z(1)*Il(1)+Z(2)*Il(2)+Z(3)*Il(3)+Z(4)*Il(4)
M(5)＝Z(1)*Il(1)+Z(5)*Il(5)
M(6)＝Z(1)*Il(1)+Z(5)*Il(5)+Z(6)*Il(6)
M(7)＝Z(1)*Il(1)+Z(7)*Il(7)
M(8)＝Z(1)*Il(1)+Z(2)*Il(2)+Z(8)*Il(8)
M(9)＝Z(1)*Il(1)+Z(2)*Il(2)+Z(3)*Il(3)+Z(9)*Il(9)
其中，V0为配电网低压侧变压器出口处电压测量值，i＝1,2,......,9；j＝1,......,9；
l＝1,2,......,i；
Vnew(1,j)
＝[V0-M(1) V0-M(2) V0-M(3) V0-M(4) V0-M(5) V0-M(6) V0-M(7) V0-M(8) V0-M(9)]
7)确定收敛条件，选择系统的误差精度为δ，判定|Vnew(1,j)-V(1,j)|≤δ是否成立，
如果不成立，则V(1,j)＝Vnew(1,j)并回到步骤2)；
否则，可得环网状电路模型中各节点负荷电压和各支路潮流如下：
第j节点电压为Uj＝V(1,j)
第i支路电流为Ili＝FC(i,1)
第i支路潮流为
其中，i＝1,2,......,b；j＝1,......,n-1；
第k支路电流为Ilk＝IL(s+b)
第k支路潮流为
其中，k＝b+1,b+2,......,b+m；s＝1,2,......,m。
本发明提供了一种环网型配电系统的潮流计算方法，针对配电网台区一个方向的配电系统进行潮流计算，首先根据环网型配电系统的特点，建立环网型配电系统电路模型，并利用矩阵形式存储配电网的原始数据，包括拓扑结构关系和电网参数，并将其作为唯一的数据输入，建立节点支路关联矩阵，并将环网型电路模型等效为对应的辐射型电路模型，进一步计算出节点负荷电流和回路电流，并用等效电流的形式来替代等效辐射型电路注入末端节点的负荷电流；最后直接利用矩阵元素的代数运算求解环网型配电系统的各支路潮流和各负荷节点电压。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (8)
1.一种环网型配电系统潮流计算方法，其特征在于：环网型配电系统潮流计算方法,其特征在于，具体包括以下几个步骤：
步骤一，根据环网型配电系统的特点，建立环网型配电系统电路模型：将配电变压器低压侧作为电势节点，等值为电压幅值和相角为恒定已知量并假定三相电压对称；馈线上所有负荷等效为末端集中的节点负荷，馈线支路采用集中参数模型，故环网型配电系统中只有一个电势节点，其余节点都是负荷节点；环网型配电系统中部分支路连接形成回路，确定环网型配电系统中回路的连支，并定义为环网状支路，环网状配电系统中除环网状支路外的其他支路定义为辐射状支路；按潮流流动的方向，设定各回路电流的绕行方向，环网状支路电流的参考方向取回路电流的绕行方向；环网型配电系统中每一条支路与系统中的两个节点相连，辐射状支路潮流从电势节点出发流向负荷节点，将辐射状支路潮流流出节点定义为始端节点，另一端的节点定义为末端节点，将环网状支路的参考电流流出节点定义为始端节点，另一端的节点定义为末端节点；
步骤二，对环网型配电系统的电路模型中拓扑结构的支路和节点进行编号，电势节点编号为0，负荷节点依次编号为1,3,4……n-1，n为环网型配电系统的节点数，包括一个电势节点，n-1个负荷节点；支路从与电势节点相连支路出发，先对辐射状支路进行编号，使辐射状支路编号和与其相连的末端节点编号相同，依次编号为1,2,3……,b，b为环网型配电系统中辐射状支路的支路数；再对环网状支路进行编号，依次继续编号为b+1,b+2,……,b+m，m为环网型配电系统中环网状支路的支路数；则环网型配电系统的支路数为(b+m)，包括b条辐射状支路，m条环网状支路；并通过初始数据矩阵DSH来表述环网型配电系统电路模型的拓扑结构和具体参数；
步骤三，通过初始数据矩阵DSH来构造环网型配电系统的节点支路关联矩阵NBH；
步骤四，根据节点支路关联矩阵NBH构造环网型配电系统电路模型的回路路径矩阵LH；
步骤五，将环网型配电系统的环网状支路断开，将环网型配电系统电路模型等效为对应的辐射型配电系统电路模型，并将注入末端节点的负荷电流用等效电流的形式来替代；
步骤六，根据等效节点支路关联矩阵NB和等效数据矩阵DS构造等效后辐射型电路模型从电势节点到负荷节点的潮流流经的路径矩阵P；
步骤七，根据等效后辐射型电路模型的路径矩阵P和等效数据矩阵DS，构造等效后辐射型电路模型从电势节点到负荷节点的潮流流经的支路阻抗矩阵ZP和节点负荷的复功率矩阵LB；
步骤八，计算环网型配电系统的各支路潮流和各负荷节点电压。
2.根据权利要求1所述的环网型配电系统潮流计算方法，其特征在于：所述的步骤二中，初始数据矩阵DSH为(b+m)行，5列的矩阵，其中第i行为：DSH(i,:)＝[x,NS(x),NR(x),Z(x),S(NR(x))]；其中，i＝1,2,......,b+m；x为第i条支路编号，NS(x)为第i条支路始端节点的编号，NR(x)为第i条支路末端节点的编号，Z(x)为第i条支路阻抗，S(NR(x))为节点NR(x)的节点负荷复功率。
3.根据权利要求1所述的环网型配电系统潮流计算方法，其特征在于：所述的步骤三中，节点支路关联矩阵NBH为n行，(b+m)列的矩阵，第j行第i列元素NBH(j,i)为：
其中，j＝1,2,......,n，i＝1,2,......,b+m。
4.根据权利要求1所述的环网型配电系统潮流计算方法，其特征在于：所述的步骤四中，回路路径矩阵LH为m行，(b+m)列的矩阵，第s行第i列元素LH(s,i)为：
其中，s＝1,2,......,m，i＝1,2,......,b+m；根据回路路径矩阵LH和初始数据矩阵DSH，构造回路阻抗矩阵ZLH，回路阻抗矩阵ZLH为m行，m列的矩阵，第s1行第s2列元素ZLH(s1,s2)为：
当s1＝s2＝s时，ZLH(s,s)表示回路s的自阻抗，等于回路s中所有支路阻抗之和，
当s1≠s2时，ZLH(s1,s2)表示回路s1与回路s2的互阻抗，等于回路s1和s2共有支路阻抗代数和，s1＝1,2,......,m；s2＝1,2,......,m；s＝1,2,......,m；其中，回路s表示环网状支路(b+s)所在的回路。
5.根据权利要求1所述的环网型配电系统潮流计算方法，其特征在于：所述的步骤五中，用等效数据矩阵DS来表述等效后辐射型电路模型的拓扑结构和具体参数，则等效数据矩阵DS为b行，5列的矩阵，第i行元素DS(i,:)＝DSH(i,:)；构造等效后辐射型电路模型的等效节点支路关联矩阵NB，等效节点支路关联矩阵NB为n行，b列的矩阵，第j行第i列元素NB(j,i)＝NBH(j,i)；其中，i＝1,2,......,b；j＝1,2,......,n，i＝1,2,......,b。
6.根据权利要求1所述的环网型配电系统潮流计算方法，其特征在于：所述的步骤六中，路径矩阵P为(n-1)行，b列的矩阵，第j行第i列元素P(j,i)为：
其中，j＝1,2,......,n-1，i＝1,2,......,b；
等效后辐射型电路模型中从电势节点0到负荷节点j的潮流流经支路的判定包括以下几步：
遍历NB矩阵中第(j+1)行元素NB(j,:)若NB(j,y)＝1，可确定支路y末端节点为j；
遍历NB矩阵中第y列元素NB(:,y)，若NB(z,y)＝-1，可确定支路y始端节点为(z-1)；
遍历NB矩阵中第z行元素NB(z,:)，若NB(z,t)＝1，可确定支路t末端节点为(z-1)；
重复以上步骤，直到支路始端节点为电势节点0。
7.根据权利要求1所述的环网型配电系统潮流计算方法，其特征在于：所述的步骤七中，支路阻抗矩阵ZP为(n-1)行，b列的矩阵，第j行第i列元素ZP(j,i)为：ZP(j,i)＝P(j,i)*Z(i)；节点负荷的复功率矩阵SB为b行，(n-1)列的矩阵，第i行第j列元素SB(i,j)为：SB(i,j)＝P(i,j)*S(j)；
其中，Z(i)＝DS(i,4)为第i条支路的阻抗；S(j)＝DS(i,5)为j节点的节点负荷复功率，支路i末端节点为j；j＝1,2,......,n-1，i＝1,2,......,b。
8.根据权利要求1所述的环网型配电系统潮流计算方法，其特征在于：所述的步骤八具体包括如下步骤：
8-1，设定各个负荷节点的电压初始值依次为V1，V2，…Vj，…Vn-1,其中，Vj表示负荷节点j的电压初始值，构造节点电压矩阵V，节点电压矩阵V为1行，(n-1)列的矩阵，其中第j列元素：V(1,j)＝Vj其中，j＝1,2,......,n-1；
8-2，计算可得回路压降矩阵VL，回路压降矩阵VL为m行，1列的矩阵，其中第s行第1列元素为：
VL(s,1)＝V(1,y)-V(1,z)；由ZLH*IL＝VL计算可得回路电流矩阵IL，回路电流矩阵IL为m行，1列的矩阵，第s回路的回路电流IL(s+b)＝IL(s,1)；
其中，环网状支路(b+s)的始端节点为y，末端节点为z，s＝1,2,......,m；
8-3，构造节点电压共轭矩阵VC，节点电压共轭矩阵为1行，n-1列的矩阵，其中第j列元素：VC(1,j)＝V*(1,j)
节点负荷复功率共轭矩阵为SC，节点负荷复功率共轭矩阵SC为b行，(n-1)列的矩阵，其中第i行第j列元素为：SC(i,j)＝SB*(i,j)
计算可得节点电流矩阵LC，节点电流矩阵LC为b行，(n-1)列的矩阵，第i行第j列元素为：LC(i,j)＝SC(i,j)/VC(1,j)；其中，i＝1,2,......,b；j＝1,......,n-1；
8-4，计算可得等效节点电流矩阵LCT，等效节点电流矩阵LCT为b行，(n-1)列的矩阵，其中第i行第j列元素为：
若LC(i,j)＝0，则LCT(i,j)＝0
若节点j不与环网状支路相连，则LCT(i,j)＝LC(i,j)
若节点j与环网状支路相连，则:
A)若节点j与环网状支路k相连，且为首端节点，则令LC(i,j)＝LC(i,j)+IL(k)
B)若节点j与环网状支路k相连，且为末端节点，则令LC(i,j)＝LC(i,j)-IL(k)
C)令k＝k+1，若k≤b+m回到步骤A)，否则，到步骤D)
D)LCT(i,j)＝LC(i,j)其中，i＝1,2,......,b，j＝1,......,n-1，k＝b+1；
8-5，计算可得支路电流矩阵FC，支路电流矩阵FC为b行，1列的矩阵，其中第i行元素为：
其中，i＝1,2,......,b，j＝1,......,n-1；
8-6，计算可得压降矩阵D，总压降矩阵M和电压矩阵T；压降矩阵D，总压降矩阵M和电压矩阵T为n-1行，b列的矩阵；
D(j,i)＝ZP(j,i)*FC(i,1)
T(j,i)＝V0-M(j,i)
计算可得新节点电压矩阵Vnew，新节点电压矩阵Vnew为1行，(n-1)列的矩阵，第j列元素：Vnew(1,j)＝min(T(j,1),T(j,2),...T(j,(n-1)))
其中，V0为辐射型配电系统低压侧变压器出口处电压测量值，i＝1,2,......,b；j＝1,......,n-1；l＝1,2,......,i；
8-7，确定收敛条件，选择系统的误差精度为δ，判定|Vnew(1,j)-V(1,j)|≤δ是否成立，如果不成立，则V(1,j)＝Vnew(1,j)并回到步骤8-2)；
否则，可得环网状电路模型中各节点电压和各支路电流，潮流如下：
第j节点电压为Uj＝V(1,j)，第i支路电流为Ili＝FC(i,1)，第i支路潮流为其中，i＝1,2,......,b；j＝1,......,n-1；第k支路电流为Ilk＝IL(s+b)，第k支路潮流为
其中，k＝b+1,b+2,......,b+m；s＝1,2,......,m。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910425652.2A CN109980651B (zh) | 2019-05-21 | 2019-05-21 | 一种环网型配电系统潮流计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910425652.2A CN109980651B (zh) | 2019-05-21 | 2019-05-21 | 一种环网型配电系统潮流计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109980651A true CN109980651A (zh) | 2019-07-05 |
| CN109980651B CN109980651B (zh) | 2022-11-15 |
Family
ID=67073873
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910425652.2A Expired - Fee Related CN109980651B (zh) | 2019-05-21 | 2019-05-21 | 一种环网型配电系统潮流计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109980651B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110445156A (zh) * | 2019-08-22 | 2019-11-12 | 浙江大学 | 含环网的三相不平衡的中低压主动配电网潮流计算方法 |
| CN110728032A (zh) * | 2019-09-20 | 2020-01-24 | 河海大学 | 一种计及环网的电-热互联综合能源系统快速潮流计算方法 |
| CN113890042A (zh) * | 2021-12-06 | 2022-01-04 | 国网浙江省电力有限公司金华供电公司 | 基于有源配电网的潮流计算方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102891485A (zh) * | 2012-10-15 | 2013-01-23 | 河海大学 | 基于序分量法的弱环配电网三相解耦潮流计算方法 |
| EP2660738A1 (en) * | 2012-04-23 | 2013-11-06 | Institute of Nuclear Energy Research Atomic Energy Council | Distribution power flow analysis system and method |
| CN103795057A (zh) * | 2014-01-23 | 2014-05-14 | 华北电力大学（保定） | 配电网潮流计算中基于搜索的配电网拓扑编号生成方法 |
| CN103793590A (zh) * | 2012-11-01 | 2014-05-14 | 同济大学 | 一种基于gpu的快速求解配电网潮流的计算方法 |
| CN106446458A (zh) * | 2016-11-01 | 2017-02-22 | 国网技术学院 | 考虑分布式电源的弱环配电网潮流计算方法 |
| CN106786602A (zh) * | 2017-02-04 | 2017-05-31 | 中南大学 | 一种配电网潮流计算方法 |
-
2019
- 2019-05-21 CN CN201910425652.2A patent/CN109980651B/zh not_active Expired - Fee Related
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2660738A1 (en) * | 2012-04-23 | 2013-11-06 | Institute of Nuclear Energy Research Atomic Energy Council | Distribution power flow analysis system and method |
| CN102891485A (zh) * | 2012-10-15 | 2013-01-23 | 河海大学 | 基于序分量法的弱环配电网三相解耦潮流计算方法 |
| CN103793590A (zh) * | 2012-11-01 | 2014-05-14 | 同济大学 | 一种基于gpu的快速求解配电网潮流的计算方法 |
| CN103795057A (zh) * | 2014-01-23 | 2014-05-14 | 华北电力大学（保定） | 配电网潮流计算中基于搜索的配电网拓扑编号生成方法 |
| CN106446458A (zh) * | 2016-11-01 | 2017-02-22 | 国网技术学院 | 考虑分布式电源的弱环配电网潮流计算方法 |
| CN106786602A (zh) * | 2017-02-04 | 2017-05-31 | 中南大学 | 一种配电网潮流计算方法 |
Non-Patent Citations (2)
| Title |
|---|
| 初壮 等: "基于回路电流与负荷节点电压方程的配电网潮流计算", 《电工技术学报》 * |
| 李宏伟 等: "考虑分布式电源的弱环配电网潮流计算方法", 《电力系统及其自动化学报》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110445156A (zh) * | 2019-08-22 | 2019-11-12 | 浙江大学 | 含环网的三相不平衡的中低压主动配电网潮流计算方法 |
| CN110728032A (zh) * | 2019-09-20 | 2020-01-24 | 河海大学 | 一种计及环网的电-热互联综合能源系统快速潮流计算方法 |
| CN110728032B (zh) * | 2019-09-20 | 2021-04-27 | 河海大学 | 一种计及环网的电-热互联综合能源系统快速潮流计算方法 |
| CN113890042A (zh) * | 2021-12-06 | 2022-01-04 | 国网浙江省电力有限公司金华供电公司 | 基于有源配电网的潮流计算方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109980651B (zh) | 2022-11-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103310065B (zh) | 计及分布式发电和储能单元的智能配网重构方法 | |
| CN109980651A (zh) | 一种环网型配电系统潮流计算方法 | |
| CN106532710B (zh) | 计及电压稳定约束的微电网潮流优化方法 | |
| CN107979092A (zh) | 一种考虑分布式电源与软开关接入的配电网动态重构方法 | |
| CN106026187B (zh) | 一种含分布式电源的配电网重构的方法及系统 | |
| CN109980650A (zh) | 一种辐射型配电系统潮流计算方法 | |
| CN111490542B (zh) | 一种多端柔性多状态开关的选址定容方法 | |
| CN106159944A (zh) | 低碳环境下基于二层规划模型的多阶段输电网规划方法 | |
| CN109494719A (zh) | 一种中低压混合配电网层次化阻抗分析方法 | |
| CN107546773A (zh) | 一种基于图论的区域多微电网动态组网方法 | |
| CN109995059B (zh) | 基于遗传算法的低压配电系统三相负荷优化平衡配置方法 | |
| CN109378819A (zh) | 基于谱聚类及模块度校验的电力系统电压分区方法及系统 | |
| CN116361603A (zh) | 一种电力系统碳排放流计算方法 | |
| CN115545290A (zh) | 一种含光伏发电配电网中的分布式储能经济优化配置方法 | |
| CN109377020B (zh) | 一种考虑配电网负荷转供能力的输电网规划方法 | |
| Ranjan et al. | Optimal conductor selection of radial distribution feeders using evolutionary programming | |
| CN110120673A (zh) | 基于戴维南等值参数辨识的分布式输配协同无功优化方法及系统 | |
| CN105914752A (zh) | 基于快速搜索与密度峰值聚类的主导节点选择方法 | |
| CN109829599A (zh) | 基于高比例可再生能源的配电网的集群划分方法及装置 | |
| CN111834996B (zh) | 电网线损计算方法和装置 | |
| CN108170952A (zh) | 基于电力电子变压器的微电网优化配置方法及装置 | |
| CN109995039B (zh) | 一种三相不平衡低压配电系统潮流计算方法 | |
| Mao et al. | Community-detection-based approach to distribution network partition | |
| CN104009469B (zh) | 一种低压配网台区网络拓扑等值图归集方法 | |
| CN105896547A (zh) | 一种风电接入下的大电网分级电压控制方法 |
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 | ||
| 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: 20221115 |