CN114696329A - 一种中压不接地配电网的三相潮流统一计算方法 - Google Patents

一种中压不接地配电网的三相潮流统一计算方法 Download PDF

Info

Publication number
CN114696329A
CN114696329A CN202210256418.3A CN202210256418A CN114696329A CN 114696329 A CN114696329 A CN 114696329A CN 202210256418 A CN202210256418 A CN 202210256418A CN 114696329 A CN114696329 A CN 114696329A
Authority
CN
China
Prior art keywords
phase
matrix
ungrounded
node
voltage
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
Application number
CN202210256418.3A
Other languages
English (en)
Other versions
CN114696329B (zh
Inventor
张逸
林金榕
李传栋
刘雄飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fuzhou University
Original Assignee
Fuzhou University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Fuzhou University filed Critical Fuzhou University
Priority to CN202210256418.3A priority Critical patent/CN114696329B/zh
Publication of CN114696329A publication Critical patent/CN114696329A/zh
Application granted granted Critical
Publication of CN114696329B publication Critical patent/CN114696329B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • 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/26—Arrangements for eliminating or reducing asymmetry in polyphase 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/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)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Operations Research (AREA)
  • Computing Systems (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明涉及一种中压不接地配电网的三相潮流统一计算方法。本发明分析不接地配电网三相潮流在节点变量、变压器导纳矩阵和零序分量方面的特殊性，设计基于功率平衡方程，约束变压器不接地绕组节点零序电流等于零，利用牛拉法构建不接地配电网三相潮流统一算法。

Description

一种中压不接地配电网的三相潮流统一计算方法
技术领域
本发明涉及一种中压不接地配电网的三相潮流统一计算方法。
背景技术
大量分布式新能源将接入中低压配电网，中低压配电网的运行和控制均面临重大挑战。同时，受到环境和人为因素的影响，新能源接入呈现较大的波动性和随机性。大量新能源的剩余功率形成逆向潮流，通过低压配网馈入中压配网，引起中压配网的电压越限、三相不平衡和稳定性问题等。面对日趋复杂的电网，对有效、准确地计算潮流提出了更高的要求。三相潮流计算作为基本仿真工具，对满足未来中压配电网的分析需求具有重要意义。
现有配电网的三相潮流计算普遍考虑为中性点直接接地的运行方式。然而，我国中压配电网存在中性点不接地的运行方式。在这种运行方式下，零序分量存在自然约束。潮流计算时，若采用不恰当的零序分量约束方式可能导致结果失真乃至不收敛。现有研究针对不接地配电网的三相潮流计算均提出了有效的实现方式，但还存在如下问题：1)基于电流平衡方程，导致算法难以与主网联合仿真；2)未形成对不接地配电网三相潮流特殊性的完整认识，导致缺乏统一的不接地配电网三相潮流零序约束方式，如文献[1]约束不接地配电网各节点的零序电压和零序电流等于零；文献[2]约束不接地配电网各节点的零序电流等于零；文献[3]约束不接地配电网中变压器不接地绕组节点的零序电压等于零。
发明内容
本发明的目的在于提供一种中压不接地配电网的三相潮流统一计算方法，该方法对具有不接地绕组的变压器导纳矩阵进行处理，删除变压器不接地绕组节点的两个功率方程，代以对应节点的零序电流方程。基于功率平衡方程，利用牛顿-拉夫逊算法实现中压不接地配电网的三相潮流计算，对未来的配电网分析具有重要意义。
为实现上述目的，本发明的技术方案是：一种中压不接地配电网的三相潮流统一计算方法，包括如下步骤：
步骤S1、输入中压不接地配电网网络的网架结构、变压器和线路等设备参数信息、节点信息、负荷信息并初始化；
步骤S2、对具有不接地绕组的变压器导纳矩阵进行处理；
步骤S3、删除变压器不接地绕组节点的两个潮流方程，替换为相应节点的零序电流实部与虚部等于零的约束方程；
步骤S4、利用功率平衡方程，基于牛顿-拉夫逊算法生成功率不平衡矩阵、雅可比矩阵和电压不平衡矩阵，进行算法迭代，并判断收敛，实现中压不接地配电网的三相潮流计算。
在本发明一实施例中，所述步骤S2具体实现方式为：
变压器三相导纳矩阵的模型如式(1)所示；
Figure BDA0003548871310000021
式中：YT为变压器三相导纳矩阵，Yff、Yft、Ytf、Ytt为YT的子导纳矩阵；
Figure BDA0003548871310000022
式中，α、β为常数；
若Yff可逆，α取0；若Yff不可逆，α取10-5～10-8；若Ytt可逆，β取0；若Ytt不可逆，β取10-5～10-8；
经过上述处理，具有不接地绕组的变压器三相导纳矩阵为：
Figure BDA0003548871310000023
在本发明一实施例中，所述步骤S3具体实现方式为：
以电网中均为PQ节点为例，满足的原始潮流方程如下式所示；
Figure BDA0003548871310000024
其中，
Figure BDA0003548871310000025
Figure BDA0003548871310000026
Gp、Bp的子矩阵定义如下所示，
Figure BDA0003548871310000031
Figure BDA0003548871310000032
Figure BDA0003548871310000033
Figure BDA0003548871310000034
Figure BDA0003548871310000035
式中，diag(*)表示利用向量*组成对角矩阵；Pp表示p相有功功率；Qp表示p相无功功率；Upx表示p相节点电压实部；Upy表示p相节点电压虚部；Uabcx表示三相节点电压实部；Uabcy表示三相节点电压虚部；G1,Ω表示节点1对于节点Ω的互电导，为3×3维矩阵，其它同理；B1,Ω表示节点1对于节点Ω的互电纳，为3×3维矩阵，其它同理；gab表示a相对b相的互电导，其它同理；bab表示a相对b相的互电纳，其它同理；Y表示三相节点导纳矩阵；B表示三相节点电纳矩阵；G表示三相节点电导矩阵；上标i、j表示节点编号；Ω表示节点总数；
以删除变压器不接地绕组节点的B、C相无功方程为例，添加对应节点的零序电流的实部和虚部满足的方程；此时，满足的潮流方程如下式所示；
Figure BDA0003548871310000041
其中，
Yp＝Gp+jBp
式中：I0x为零序电流实部；I0y为零序电流虚部；Ypx表示Yp的实部；Ypy表示Yp的虚部；Uay(i1)表示按照向量i1的元素提取Uay对应行向量组成新矩阵，其它同理；i1为网络的节点编号向量，但不包括变压器不接地绕组节点的编号；i2为变压器不接地绕组节点的节点编号向量。
在本发明一实施例中，所述步骤S4具体实现方式为：
牛顿-拉夫逊潮流方法的基础方程为：
F＝JΔx
功率不平衡矩阵F为：
Figure BDA0003548871310000042
雅可比矩阵J为：
Figure BDA0003548871310000051
电压不平衡矩阵Δx为：
Figure BDA0003548871310000052
式中：Pabc表示三相有功功率；
Figure BDA0003548871310000053
表示PV节点给定的三相电压。
其中，功率不平衡矩阵F已经删除变压器不接地绕组节点B、C相的ΔQb和ΔQc，添加这些节点的ΔI0x和ΔI0y；雅可比矩阵中已经删除变压器不接地绕组节点B、C相的
Figure BDA0003548871310000054
和
Figure BDA0003548871310000055
添加这些节点的
Figure BDA0003548871310000056
和
Figure BDA0003548871310000057
Figure BDA0003548871310000058
表示按照向量i1的元素提取
Figure BDA0003548871310000059
对应行向量组成新矩阵，其它同理；
与零序分量相关的偏导矩阵元素按下式生成：
Figure BDA00035488713100000510
相较于现有技术，本发明具有以下有益效果：
本发明提出一种中压不接地配电网的三相潮流计算方法。该方法对具有不接地绕组的变压器导纳矩阵进行处理，删除变压器不接地绕组节点的2个功率方程，代以对应节点的零序电流方程。基于功率平衡方程，利用牛顿-拉夫逊算法实现中压不接地配电网的三相潮流计算，对未来的配电网分析具有重要意义。
附图说明
图1为本发明方法流程图。
图2为本发明YNy变压器结构示例。
图3为5节点系统。
具体实施方式
下面结合附图，对本发明的技术方案进行具体说明。
如图1所示，本发明一种中压不接地配电网的三相潮流统一计算方法，包括如下步骤：
步骤S1、输入中压不接地配电网网络的网架结构、变压器和线路等设备参数信息、节点信息、负荷信息并初始化；
步骤S2、对具有不接地绕组的变压器导纳矩阵进行处理；
步骤S3、删除变压器不接地绕组节点的两个潮流方程，替换为相应节点的零序电流实部与虚部等于零的约束方程；
步骤S4、利用功率平衡方程，基于牛顿-拉夫逊算法生成功率不平衡矩阵、雅可比矩阵和电压不平衡矩阵，进行算法迭代，并判断收敛，实现中压不接地配电网的三相潮流计算。
本发明的具体内容如下：
中压不接地配电网的三相潮流特殊性：
1)节点变量
传统接地网络的三相潮流，网络中各节点拥有12个变量，已知6个，待求6个。例如，PQ节点已知的是三相的有功功率和无功功率，待求的是三相电压的实部和虚部(或电压幅值和相角)，满足方程组如式(1)所示。
Figure BDA0003548871310000061
其中，
Figure BDA0003548871310000071
Figure BDA0003548871310000072
Gp、Bp的子矩阵定义如式(4)-(8)所示，
Figure BDA0003548871310000073
Figure BDA0003548871310000074
Figure BDA0003548871310000075
Figure BDA0003548871310000076
Figure BDA0003548871310000077
式中，diag(*)表示利用向量*组成对角矩阵；Pp表示p相有功功率；Qp表示三相无功功率；Upx表示节点p相电压实部；Upy表示节点p相电压虚部，p＝a,b,c；Uabcx表示节点三相电压实部；Uabcy表示节点三相电压虚部；G1,Ω表示节点1对于节点Ω的互电导，为3×3维矩阵，其它同理；B1,Ω表示节点1对于节点Ω的互电纳，为3×3维矩阵，其它同理；gab表示a相对b相的互电导，其它同理；bab表示a相对b相的互电纳，其它同理；Y表示节点三相导纳矩阵；B表示三相节点电纳矩阵；G表示三相节点电导矩阵；上标i、j表示节点编号；Ω表示节点总数。
需要指出，式(1)中的方程数等于待求变量数，数学上确定存在解。然而，在面对不接地网络时，零序分量自然为零，节点的变量之间并不完全独立，方程减少自由度。此时，若添加约束方程，将使方程数大于待求变量数。此类方程组在数学上定义为超定方程，不能保证有解。可见，若需保证方程组有解，需删除不独立的节点变量及其对应的方程。
2)变压器导纳矩阵
以YNy变压器模型为例进行说明，相应的导纳矩阵如式(9)-(11)所示。从矩阵形式可判断YNy变压器导纳矩阵和其子矩阵均是奇异的。
Figure BDA0003548871310000081
Figure BDA0003548871310000082
Figure BDA0003548871310000083
式中：y1为三相变压器的正序导纳；k表示三相变压器的变比。
可见，由于具有不接地绕组的变压器导纳矩阵是奇异的，导致无法直接应用于不接地网络的三相潮流计算。
3)零序分量
对具有不接地绕组变压器的零序分量进行分析，仍以YNy变压器为例，如图2所示。由于变压器不接地绕组在结构上不存在对地通路，导致零序电流没有回路，即ITaj+ITbj+ITcj＝3IT0j＝0，IT0表示变压器绕组的零序电流。当j侧流过零序电流时，即ILaj+ILbj+ILcj＝3IL0j≠0，此零序电流只能通过对地电纳形成回路，且不接地网络的对地电纳数值很小，从而在变压器不接地绕组侧形成较大的零序电压。
针对以上关于中压不接地配网三相潮流的特点，本发明对具有不接地绕组的变压器导纳矩阵进行处理，约束不接地网络中变压器的不接地绕组节点的零序电流为零。由于自然存在的零序分量约束，使节点变量不独立，本发明删除变压器不接地绕组节点的2个潮流方程，使传统接地网络的三相潮流适用于不接地网络。
本方法的具体过程如下：
步骤S1、输入中压不接地配电网网络的网架结构、变压器和线路等设备参数信息、节点信息、负荷信息并初始化；
步骤S2、对具有不接地绕组的变压器导纳矩阵进行处理；
步骤S3、删除变压器不接地绕组节点的两个潮流方程，替换为相应节点的零序电流实部与虚部等于零的约束方程；
步骤S4、利用功率平衡方程，基于牛顿-拉夫逊算法生成功率不平衡矩阵、雅可比矩阵和电压不平衡矩阵，进行算法迭代，并判断收敛，实现中压不接地配电网的三相潮流计算。
进一步地，所述步骤S2的具体内容为：
三相变压器导纳矩阵的模型如式(9)所示。提取式(9)的对角矩阵，依据式(12)对其进行处理：
Figure BDA0003548871310000091
式中，α、β为常数。
若Yff可逆，α取0；若Yff不可逆，α取10-5～10-8；若Ytt可逆，β取0；若Ytt不可逆，β取10-5～10-8。
经过上述处理，具有不接地绕组的变压器的三相导纳矩阵为：
Figure BDA0003548871310000092
进一步地，所述步骤S3的具体内容为：
以网络均为PQ节点为例，原始的潮流方程如式(1)所示。以删除变压器不接地绕组节点的B、C相无功潮流方程为例，添加相应节点的零序电流的实部和虚部满足的方程。此时，该节点满足的潮流方程如式(14)所示。
Figure BDA0003548871310000101
其中，
Yp＝Gp+jBp (15)
式中：I0x为零序电流实部；I0y为零序电流虚部；Ypx表示Yp的实部；Ypy表示Yp的虚部；Uay(i1)表示按照向量i1的元素提取Uay对应行向量组成新矩阵，其它同理；i1为网络的节点编号向量，但不包括变压器不接地绕组节点的编号；i2为变压器不接地绕组节点的节点编号向量。
进一步地，所述步骤S4的具体内容为：
牛顿-拉夫逊潮流方法的基础方程为：
F＝JΔx (16)
功率不平衡矩阵F为：
Figure BDA0003548871310000102
雅可比矩阵J为：
Figure BDA0003548871310000111
电压不平衡矩阵Δx为：
Figure BDA0003548871310000112
式中：Pabc表示三相有功功率；
Figure BDA0003548871310000113
表示PV节点给定的三相电压。
其中，功率不平衡矩阵F已经删除变压器不接地绕组节点B、C相的ΔQb和ΔQc，添加这些节点的ΔI0x和ΔI0y；雅可比矩阵中已经删除变压器不接地绕组节点B、C相的
Figure BDA0003548871310000114
和
Figure BDA0003548871310000115
添加这些节点的
Figure BDA0003548871310000116
和
Figure BDA0003548871310000117
Figure BDA0003548871310000118
表示按照向量i1的元素提取
Figure BDA0003548871310000119
对应行向量组成新矩阵，其它同理；
Figure BDA00035488713100001110
为了将本方法阐述地更加清晰，结合下述例子，对本方法进一步详细说明。本发明的示意性实施方式及其说明仅用于解释本发明，并不作为对本发明的限定。
假定一个如图3所示的5节点系统模型。
1)确定网络的输入参数，包括网架结构、2台变压器参数、2条线路参数、负荷参数、5个节点的类型信息，确定算法的收敛精度，启动电压；
2)根据公式(12)对YNd11和Dyn11变压器的导纳矩阵进行处理，其中，YNd11变压器的Yff可逆，Ytt不可逆，Dyn11变压器的Yff不可逆，Ytt可逆；
3)删除节点3和节点4的B、C相无功潮流方程，替换为相应节点的零序电流方程，方程如式(14)所示，式中各个参数如式(21)所示。
Figure BDA0003548871310000121
4)生成牛顿-拉夫逊潮流算法的相关矩阵。牛拉法相关矩阵的参数如式(22)-(25)所示，其中雅可比矩阵的常规偏导矩阵生成公式见文献[4]。
功率不平衡矩阵F为：
Figure BDA0003548871310000131
其中，
Figure BDA0003548871310000132
雅可比矩阵中与零序分量相关的偏导矩阵为：
Figure BDA0003548871310000133
电压不平衡矩阵Δx为：
Figure BDA0003548871310000134
参考文献：
[1]索南加乐,付伟,刘文涛,孙丹丹,王向兵,齐军.中性点不接地配网潮流算法研究[J].中国电机工程学报,2006(10):86-91.
[2]颜伟,杨焕燕,赵霞,余娟,杜跃明,王俊杰,邓亚拓.中性点不接地配网的三相潮流模型[J].中国电机工程学报,2013,33(19):135-141+10.
[3]高亚静,苗宏佳,吴文传,巨云涛,李飞.处理不接地配电网三相潮流不收敛的变压器建模新方法[J].电力系统自动化,2014,38(18):53-57.
[4]李传栋,林金榕,严翠云,杜培.一种电网三相谐波潮流相量矩阵计算方法[P].福建省：CN112670987A,2021-04-16.。
以上是本发明的较佳实施例，凡依本发明技术方案所作的改变，所产生的功能作用未超出本发明技术方案的范围时，均属于本发明的保护范围。

Claims (4)

1.一种中压不接地配电网的三相潮流统一计算方法，其特征在于，包括如下步骤：
步骤S1、输入中压不接地配电网网络的网架结构，变压器和线路的设备参数信息、节点信息、负荷信息并初始化；
步骤S2、对具有不接地绕组的变压器导纳矩阵进行处理；
步骤S3、删除变压器不接地绕组节点的两个潮流方程，替换为相应节点的零序电流实部与虚部等于零的约束方程；
步骤S4、利用功率平衡方程，基于牛顿-拉夫逊算法生成功率不平衡矩阵、雅可比矩阵和电压不平衡矩阵，进行算法迭代，并判断收敛，实现中压不接地配电网的三相潮流计算。
2.根据权利要求1所述的一种中压不接地配电网的三相潮流统一计算方法，其特征在于，所述步骤S2具体实现方式为：
变压器三相导纳矩阵的模型如式(1)所示；
Figure FDA0003548871300000011
式中：YT为变压器三相导纳矩阵；Yff、Yft、Ytf、Ytt为YT的子导纳矩阵，均为3×3维；
Figure FDA0003548871300000012
式中，α、β为常数；
若Yff可逆，α取0；若Yff不可逆，α取10-5～10-8；若Ytt可逆，β取0；若Ytt不可逆，β取10-5～10-8；
经过上述处理，具有不接地绕组的变压器三相导纳矩阵为：
Figure FDA0003548871300000013
3.根据权利要求1所述的一种中压不接地配电网的三相潮流统一计算方法，其特征在于，所述步骤S3具体实现方式为：
以电网中均为PQ节点为例，满足的原始潮流方程如下式所示：
Figure FDA0003548871300000021
其中，
Figure FDA0003548871300000022
Figure FDA0003548871300000023
Gp、Bp的子矩阵定义如下所示，
Figure FDA0003548871300000024
Figure FDA0003548871300000025
Figure FDA0003548871300000026
Figure FDA0003548871300000027
Figure FDA0003548871300000028
式中：diag(*)表示利用向量*组成对角矩阵；Pp表示p相有功功率；Qp表示p相无功功率；Upx表示p相节点电压实部；Upy表示p相节点电压虚部；Uabcx表示三相节点电压实部；Uabcy表示三相节点电压虚部；G1,Ω表示节点1对于节点Ω的互电导，为3×3维矩阵，其它同理；B1,Ω表示节点1对于节点Ω的互电纳，为3×3维矩阵，其它同理；gab表示a相对b相的互电导，其它同理；bab表示a相对b相的互电纳，其它同理；Y表示三相节点导纳矩阵；B表示三相节点电纳矩阵；G表示三相节点电导矩阵；上标i、j表示节点编号；Ω表示节点总数；
以删除变压器不接地绕组节点的B、C相无功方程为例，添加对应节点的零序电流的实部和虚部满足的方程；此时，网络满足的潮流方程如下式所示；
Figure FDA0003548871300000031
其中，
Yp＝Gp+jBp
式中：I0x为零序电流实部；I0y为零序电流虚部；Ypx表示Yp的实部；Ypy表示Yp的虚部；Uay(i1)表示按照向量i1的元素提取Uay的对应行向量组成新矩阵，其它同理；i1为网络的节点编号向量，但不包括变压器不接地绕组节点的编号；i2为变压器不接地绕组节点的节点编号向量。
4.根据权利要求1所述的一种中压不接地配电网的三相潮流统一计算方法，其特征在于，所述步骤S4具体实现方式为：
牛顿-拉夫逊潮流方法的基础方程为：
F＝JΔx
功率不平衡矩阵F为：
Figure FDA0003548871300000041
雅可比矩阵J为：
Figure FDA0003548871300000042
电压不平衡矩阵Δx为：
Figure FDA0003548871300000043
式中：Pabc表示三相有功功率；
Figure FDA0003548871300000044
表示PV节点给定的三相电压；
其中，功率不平衡矩阵F已经删除变压器不接地绕组节点B、C相的ΔQb和ΔQc，添加这些节点的ΔI0x和ΔI0y；雅可比矩阵中已经删除变压器不接地绕组节点B、C相的
Figure FDA0003548871300000045
和
Figure FDA0003548871300000046
添加这些节点的
Figure FDA0003548871300000047
和
Figure FDA0003548871300000048
表示按照向量i1的元素提取
Figure FDA0003548871300000049
对应的行向量组成新矩阵，其它同理；
与零序分量相关的偏导矩阵元素按下式生成：
Figure FDA0003548871300000051
CN202210256418.3A 2022-03-16 2022-03-16 一种中压不接地配电网的三相潮流统一计算方法 Active CN114696329B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210256418.3A CN114696329B (zh) 2022-03-16 2022-03-16 一种中压不接地配电网的三相潮流统一计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210256418.3A CN114696329B (zh) 2022-03-16 2022-03-16 一种中压不接地配电网的三相潮流统一计算方法

Publications (2)

Publication Number Publication Date
CN114696329A true CN114696329A (zh) 2022-07-01
CN114696329B CN114696329B (zh) 2024-07-02

Family

ID=82139870

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210256418.3A Active CN114696329B (zh) 2022-03-16 2022-03-16 一种中压不接地配电网的三相潮流统一计算方法

Country Status (1)

Country Link
CN (1) CN114696329B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117828242A (zh) * 2024-03-06 2024-04-05 国网上海能源互联网研究院有限公司 一种矩阵压缩存储格式的fpga潮流计算系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103701123A (zh) * 2014-01-10 2014-04-02 贵州电网公司信息通信分公司 一种用于不接地配电网的高斯-塞德尔三相潮流计算方法
US20150199301A1 (en) * 2014-01-10 2015-07-16 Tsinghua University Method and device for obtaining three-phase power flow of power distribution network with ungrounded transformers
CN105703363A (zh) * 2016-04-11 2016-06-22 重庆大学 一种基于线电压的不接地配电网三相潮流计算方法
CN106026082A (zh) * 2016-06-27 2016-10-12 重庆大学 一种不接地配电网的三相潮流计算方法
CN111355236A (zh) * 2020-04-10 2020-06-30 大连海事大学 一种计及中性点电压变量的配电网三相潮流计算方法
CN112383047A (zh) * 2020-10-20 2021-02-19 重庆大学 基于中低压网络测控特性的中压馈线电网三相潮流计算方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103701123A (zh) * 2014-01-10 2014-04-02 贵州电网公司信息通信分公司 一种用于不接地配电网的高斯-塞德尔三相潮流计算方法
US20150199301A1 (en) * 2014-01-10 2015-07-16 Tsinghua University Method and device for obtaining three-phase power flow of power distribution network with ungrounded transformers
CN105703363A (zh) * 2016-04-11 2016-06-22 重庆大学 一种基于线电压的不接地配电网三相潮流计算方法
CN106026082A (zh) * 2016-06-27 2016-10-12 重庆大学 一种不接地配电网的三相潮流计算方法
CN111355236A (zh) * 2020-04-10 2020-06-30 大连海事大学 一种计及中性点电压变量的配电网三相潮流计算方法
CN112383047A (zh) * 2020-10-20 2021-02-19 重庆大学 基于中低压网络测控特性的中压馈线电网三相潮流计算方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
颜伟;杨焕燕;赵霞;余娟;杜跃明;王俊杰;邓亚拓;: "中性点不接地配网的三相潮流模型", 中国电机工程学报, no. 19, 5 July 2013 (2013-07-05), pages 135 - 141 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117828242A (zh) * 2024-03-06 2024-04-05 国网上海能源互联网研究院有限公司 一种矩阵压缩存储格式的fpga潮流计算系统
CN117828242B (zh) * 2024-03-06 2024-05-14 国网上海能源互联网研究院有限公司 一种矩阵压缩存储格式的fpga潮流计算系统

Also Published As

Publication number Publication date
CN114696329B (zh) 2024-07-02

Similar Documents

Publication Publication Date Title
Chen et al. Three-phase cogenerator and transformer models for distribution system analysis
Penido et al. Three-phase power flow based on four-conductor current injection method for unbalanced distribution networks
Balamurugan et al. Review of power flow studies on distribution network with distributed generation
CN109063408B (zh) 一种基于直算法的电磁机电暂态仿真方法
CN106066918B (zh) 一种含分布式电源和非线性负荷的短路电流计算方法
Chen et al. Hybrid three phase load flow
CN108270248B (zh) 一种基于准稳态灵敏度分析的有功安全校正控制方法
Liu et al. Linear power flow calculation of distribution networks with distributed generation
CN107204617A (zh) 基于线性规划的直角坐标形式的区间潮流计算方法
CN113206517B (zh) 孤岛微电网频率和电压恢复控制方法、装置、设备和介质
CN114696329A (zh) 一种中压不接地配电网的三相潮流统一计算方法
CN111355236A (zh) 一种计及中性点电压变量的配电网三相潮流计算方法
Nduka et al. A robust augmented nodal analysis approach to distribution network solution
Jiang et al. Power–current hybrid rectangular formulation for interior-point optimal power flow
CN114204564B (zh) 一种含逆变型新能源电网的短路电流计算方法及装置
CN103956735A (zh) 一种分布式发电系统的谐波潮流分析方法
CN108599167A (zh) 一种辐射状配电网的线性化潮流计算方法
He et al. Fast short-circuit current calculation of unbalanced distribution networks with inverter-interfaced distributed generators
Tan et al. Automatic generation and parameter verification of large-scale EMT simulation models based on TSP projects
CN109494748A (zh) 基于节点类型和修正的雅可比矩阵的牛顿法潮流计算方法
CN106410811B (zh) 首次迭代小阻抗支路端点改变雅可比矩阵的潮流计算方法
CN114188945B (zh) 一种含光伏电源的配电网短路电流计算方法及装置
Campello et al. Reviewing the large electrical network equivalent methods under development for electromagnetic transient studies in the Brazilian Interconnected power system
CN113030643B (zh) 一种配网电压源型分布式电源的故障分析方法及系统
CN104600700B (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