CN105974271B - 一种配电线路故障端口的节点阻抗矩阵计算方法 - Google Patents
一种配电线路故障端口的节点阻抗矩阵计算方法 Download PDFInfo
- Publication number
- CN105974271B CN105974271B CN201610590130.4A CN201610590130A CN105974271B CN 105974271 B CN105974271 B CN 105974271B CN 201610590130 A CN201610590130 A CN 201610590130A CN 105974271 B CN105974271 B CN 105974271B
- Authority
- CN
- China
- Prior art keywords
- node
- matrix
- parallel branch
- impedance
- load
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/081—Locating faults in cables, transmission lines, or networks according to type of conductors
- G01R31/086—Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution networks, i.e. with interconnected conductors
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/088—Aspects of digital computing
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种配电线路故障端口的节点阻抗矩阵计算方法,它包括以下过程:忽略配电线路中并联支路和故障端口之后串联支路的线路阻抗,将并联支路的负荷阻抗添加到对应节点上,将故障端口之后串联支路的负荷阻抗合并到故障端口节点上,得到故障主回路的等效电路;故障主回路的等效电路采用线路回路阻抗矩阵与并联支路节点导纳矩阵来表示;计算各个节点的并联支路节点导纳矩阵;根据故障主回路的等效电路计算故障端口节点阻抗矩阵。本发明提出了一种故障端口节点阻抗矩阵计算的新方法,它考虑到并联支路的影响,通过简单的矩阵运算即可求出故障补偿端口节点阻抗矩阵,避免了复杂的高阶矩阵的求逆运算,并且计算结果精确。
Description
技术领域
本发明涉及一种配电线路故障端口的节点阻抗矩阵计算方法,属于配电网故障计算技术领域。
背景技术
传统的电力系统故障计算通常采用节点导纳矩阵解法,这种方法没有考虑到配电网的特点,矩阵分解因子表的形成比较复杂,计算量大。配电网的故障计算还可以采用基于相分量的端口补偿法,该方法充分利用配电网一般是辐射状网络的特点,可以使用前推回代法计算,得到故障端口开路电压,而端口阻抗矩阵可以根据端口节点阻抗矩阵求出,因此极大地减少了短路故障计算工作量。
目前有关端口补偿法计算配电网故障的研究,在计算节点阻抗时,都不考虑线路中并联的无功补偿电容器与负荷(简称并联支路)的影响,因此,节点自阻抗就是节点到根节点之间回路的串联阻抗之和,两个节点之间的互阻抗是两个节点到根节点的回路的公共串联阻抗之和。这种节点阻抗的计算方法比较简单,但是,由于忽略了并联支路的影响,无法提高故障电流计算的精度。考虑并联支路后,通常可根据配电网结构与参数建立节点导纳矩阵,再将其求逆即可获得节点阻抗矩阵。因为涉及矩阵求逆问题,如果节点比较多,矩阵维数就比较高,计算过程复杂,计算量大,达不到简化计算的目的。
发明内容
针对现有技术的不足,本发明提出了一种配电线路故障端口的节点阻抗矩阵计算方法,其不仅计算方法简单,而且计算精度高。
本发明解决其技术问题采取的技术方案是:一种配电线路故障端口的节点阻抗矩阵计算方法,配电线路故障端口的短路电流采用端口补偿法计算,其特征是,所述节点阻抗矩阵计算方法包括以下过程:
忽略配电线路中并联支路和故障端口之后串联支路的线路阻抗,将并联支路的负荷阻抗添加到对应节点上,将故障端口之后串联支路的负荷阻抗合并到故障端口节点上,得到故障主回路的等效电路;
故障主回路的等效电路采用线路回路阻抗矩阵与并联支路节点导纳矩阵来表示,线路回路阻抗矩阵表示为Z0、Z1、Z2……Zi,各个节点的并联支路节点导纳矩阵表示为Y0、Y1、Y2……Yi,Z0为电源的线路回路阻抗矩阵,Z1为节点1到电源之间线路段的线路回路阻抗矩阵,Zi为故障端口节点i与节点i-1之间线路段的线路回路阻抗矩阵,Y0为根节点的并联支路节点导纳矩阵,Y1为节点1的并联支路节点导纳矩阵,Yi为故障端口节点i的并联支路节点导纳矩阵;
计算各个节点的并联支路节点导纳矩阵;
根据故障主回路的等效电路计算故障端口节点阻抗矩阵。
进一步地,所述并联支路节点导纳矩阵的计算过程包括以下步骤:
把每一相负荷等效为一个恒阻抗,阻抗值根据负荷的有功功率与无功功率求出,负荷阻抗ZL为:
负荷导纳的计算公式为:
式中,UN为负荷两端的标称电压,PL为负荷的有功功率,QL为负荷额无功功率,j为虚数单位;
当负荷接线方式采用Y0接线方式时,并联支路线路电压与并联支路电流关系可直接写为:
式(3)简写为:
中性点直接接地时的并联支路节点导纳矩阵YP为:
式中,YP为中性点直接接地时并联支路节点导纳矩阵,其中三个对角线的元素为每一相支路的导纳;是支路电流矩阵向量;是支路电压矩阵向量;
当负荷接线方式采用Y接线方式时,并联支路线路电压与并联支路电流关系可直接写为:
式中,Un是并联支路中性点对地电压;
因为中性点不接地,三相并联支路电流之和等于零,Un可根据电压以及并联支路阻抗参数求出,Un的表达式为:
将Un代入式(6)中得到:
式(8)简写为:
中性点不接地时的并联支路节点导纳矩阵Ype为:
式中,Ype为中性点不接地时的并联支路节点导纳矩阵。
式中,Ype为中性点不接地时的并联支路节点导纳矩
优选地,当负荷采用星型接线时UN为负荷两端的额定相电压,当采用三角形接线时UN为负荷两端的额定线电压。
优选地,将负荷接线方式采用三角形接线方式的等效电路转换成星型中性点不接地接线方式的等效电路,三角形接线方式的并联支路节点导纳矩阵采用中性点不接地时的并联支路节点导纳矩阵来求解。
进一步地,故障端口节点阻抗矩阵的计算过程包括以下步骤:
1)计算从根节点向电源侧看进去的电源侧节点阻抗矩阵,将电源的中性点视为直接接地,则根节点的电源侧节点阻抗矩阵就是电源回路阻抗矩阵Z0,即:
ZS0=Z0;
2)计算根节点加入并联支路后的节点阻抗矩阵,设根节点电压矩阵向量为流入系统与并联支路的电流矩阵向量为则有:
增加并联支路后的根节点的节点阻抗矩阵为:
3)计算节点1的电源侧节点阻抗矩阵,由于节点1电压与根节点电压满足以下关系:
节点1的电源侧节点阻抗矩阵为:
ZS1=Zn0+Z1
4)计算节点1加入并联支路后的节点阻抗矩阵,增加并联支路后的节点1的节点阻抗矩阵为:
5)依次计算出节点2至故障端口节点i的电源侧节点阻抗矩阵和节点阻抗矩阵,故障端口节点i的电源侧节点矩阵为:
ZSi=Zn(i-1)+Zi
增加并联支路后的故障端口节点i的节点阻抗矩阵为:
其中,Z0为电源的线路回路阻抗矩阵,Z1为节点1到电源之间线路段的线路回路阻抗矩阵,Zi为故障端口节点i与节点i-1之间线路段的线路回路阻抗矩阵,ZS0为根节点的电源侧节点阻抗矩阵,ZS1为节点1的电源侧节点阻抗矩阵,ZSi为故障端口节点i的电源侧节点阻抗矩阵,Zn0为根节点并联支路的中性点与地之间的阻抗,Zn1为节点1并联支路的中性点与地之间的阻抗,Zni为故障端口节点i并联支路的中性点与地之间的阻抗,Y0为根节点的并联支路节点导纳矩阵,Y1为节点1的并联支路节点导纳矩阵,Yi为故障端口节点i的并联支路节点导纳矩阵。
优选地,并联支路采用Y0接线时,并联支路的中性点与地之间的阻抗Zn=0,否则Zn取无穷大。
优选地,所述并联支路包括无功补偿器并联支路和/或负荷并联支路。
本发明的有益效果如下:
本发明分别求解出中性点接地和中性点不接地时的并联支路节点导纳矩阵;考虑到实际中一般负荷阻抗远大于线路阻抗,忽略故障主回路之外的分支线路阻抗(包括故障端口后的线路阻抗),将故障后负荷阻抗叠加到故障点;同时根据故障主回路的配电线路回路阻抗矩阵以及并联支路节点导纳矩阵,通过简单的阻抗串并联计算实现端口节点阻抗矩阵计算,避免了复杂的导纳矩阵求逆运算。
本发明提出了一种故障端口节点阻抗矩阵计算的新方法,它考虑到并联支路的影响,通过简单的矩阵运算即可求出故障补偿端口节点阻抗矩阵,避免了复杂的高阶矩阵的求逆运算,并且计算结果精确。
附图说明
图1为本发明的方法流程图;
图2(a)为中性点接地的星型接线方式并联支路的等效电路;
图2(b)为中性点不接地的星型接线方式并联支路的等效电路;
图2(c)为三角型接线方式并联支路中性点接地时的等效电路;
图3为5节点配电线路的计算故障端口节点阻抗矩阵简化线路示意图;
图4为图3所示配电线路的故障主回路等效电路;
图5为图3所示配电线路的补偿端口主回路等效电路。
具体实施方式
为能清楚说明本方案的技术特点,下面通过具体实施方式,并结合其附图,对本发明进行详细阐述。下文的公开提供了许多不同的实施例或例子用来实现本发明的不同结构。为了简化本发明的公开,下文中对特定例子的部件和设置进行描述。此外,本发明可以在不同例子中重复参考数字和/或字母。这种重复是为了简化和清楚的目的,其本身不指示所讨论各种实施例和/或设置之间的关系。应当注意,在附图中所图示的部件不一定按比例绘制。本发明省略了对公知组件和处理技术及工艺的描述以避免不必要地限制本发明。
如图1所示,本发明的一种配电线路故障端口的节点阻抗矩阵计算方法,配电线路故障端口的短路电流采用端口补偿法计算,所述节点阻抗矩阵计算方法包括以下过程:
忽略配电线路中并联支路和故障端口之后串联支路的线路阻抗,将并联支路的负荷阻抗添加到对应节点上,将故障端口之后串联支路的负荷阻抗合并到故障端口节点上,得到故障主回路的等效电路;
故障主回路的等效电路采用线路回路阻抗矩阵与并联支路节点导纳矩阵来表示,线路回路阻抗矩阵表示为Z0、Z1、Z2……Zi,各个节点的并联支路节点导纳矩阵表示为Y0、Y1、Y2……Yi,Z0为电源的线路回路阻抗矩阵,Z1为节点1到电源之间线路段的线路回路阻抗矩阵,Zi为故障端口节点i与节点i-1之间线路段的线路回路阻抗矩阵,Y0为根节点的并联支路节点导纳矩阵,Y1为节点1的并联支路节点导纳矩阵,Yi为故障端口节点i的并联支路节点导纳矩阵;
计算各个节点的并联支路节点导纳矩阵;
根据故障主回路的等效电路计算故障端口节点阻抗矩阵。
进一步地,所述并联支路节点导纳矩阵的计算过程包括以下步骤:
把每一相负荷等效为一个恒阻抗,阻抗值根据负荷的有功功率与无功功率求出,负荷阻抗ZL为:
负荷导纳的计算公式为:
式中,UN为负荷两端的标称电压,PL为负荷的有功功率,QL为负荷额无功功率,j为虚数单位;当负荷采用星型接线时UN为负荷两端的额定相电压,当采用三角形接线时UN为负荷两端的额定线电压。
当负荷接线方式采用Y0接线方式时,并联支路线路电压与并联支路电流关系可直接写为:
式(3)简写为:
中性点直接接地时的并联支路节点导纳矩阵YP为:
式中,YP为中性点直接接地时并联支路节点导纳矩阵,其中三个对角线的元素为每一相支路的导纳;是支路电流矩阵向量;是支路电压矩阵向量;
当负荷接线方式采用Y接线方式时,并联支路线路电压与并联支路电流关系可直接写为:
式中,Un是并联支路中性点对地电压;
因为中性点不接地,三相并联支路电流之和等于零,Un可根据电压以及并联支路阻抗参数求出,Un的表达式为:
将Un代入式(6)中得到:
式(8)简写为:
中性点不接地时的并联支路节点导纳矩阵Ype为:
式中,Ype为中性点不接地时的并联支路节点导纳矩阵。
优选地,当负荷采用星型接线时UN为负荷两端的额定相电压,当采用三角形接线时UN为负荷两端的额定线电压。
进一步地,故障端口节点阻抗矩阵的计算过程包括以下步骤:
1)计算从根节点向电源侧看进去的电源侧节点阻抗矩阵,将电源的中性点视为直接接地,则根节点的电源侧节点阻抗矩阵就是电源回路阻抗矩阵Z0,即:
ZS0=Z0;
2)计算根节点加入并联支路后的节点阻抗矩阵,设根节点电压矩阵向量为流入系统与并联支路的电流矩阵向量为则有:
增加并联支路后的根节点的节点阻抗矩阵为:
3)计算节点1的电源侧节点阻抗矩阵,由于节点1电压与根节点电压满足以下关系:
节点1的电源侧节点阻抗矩阵为:
ZS1=Zn0+Z1
4)计算节点1加入并联支路后的节点阻抗矩阵,增加并联支路后的节点1的节点阻抗矩阵为:
5)依次计算出节点2至故障端口节点i的电源侧节点阻抗矩阵和节点阻抗矩阵,故障端口节点i的电源侧节点矩阵为:
ZSi=Zn(i-1)+Zi
增加并联支路后的故障端口节点i的节点阻抗矩阵为:
其中,Z0为电源的线路回路阻抗矩阵,Z1为节点1到电源之间线路段的线路回路阻抗矩阵,Zi为故障端口节点i与节点i-1之间线路段的线路回路阻抗矩阵,ZS0为根节点的电源侧节点阻抗矩阵,ZS1为节点1的电源侧节点阻抗矩阵,ZSi为故障端口节点i的电源侧节点阻抗矩阵,Zn0为根节点并联支路的中性点与地之间的阻抗,Zn1为节点1并联支路的中性点与地之间的阻抗,Zni为故障端口节点i并联支路的中性点与地之间的阻抗,Y0为根节点的并联支路节点导纳矩阵,Y1为节点1的并联支路节点导纳矩阵,Yi为故障端口节点i的并联支路节点导纳矩阵。
下面对本发明涉及到的并联支路统一建模和故障端口节点阻抗矩计算进行详细说明。
在计算放射式配电线路故障补偿端口(简称故障端口)节点阻抗时,把故障端口到根节点之间的串联配电线路称为故障主回路。由于一般负荷阻抗远大于线路阻抗,因此,在计算故障端口节点矩阵时,本发明忽略故障主回路之外的分支线路阻抗(包括故障端口后的线路阻抗),只考虑其所带的负荷阻抗。这样,计算故障端口节点阻抗的等效电路就是一个不带分支线路的从故障端口到根节点的简单放射式线路。根据故障主回路的配电线路回路阻抗矩阵以及并联支路节点导纳矩阵,可计算出故障端口的节点阻抗矩阵。
1、并联支路统一建模
中压配电网故障计算涉及的并联支路包括无功补偿器和负荷。因采用恒阻抗模型后,负荷模型与无功补偿电容器模型是统一的。所以在下面的叙述中,不再区分负荷与无功补偿电容器,统称为负荷。
在进行配电网故障计算时,一般忽略负荷的非线性现象以及三相负荷的耦合,把每一相负荷等效为一个恒阻抗,阻抗值根据负荷的有功功率与无功功率求出,即
式中,UN为负荷两端的标称电压,PL为负荷的有功功率,QL为负荷额无功功率,j为虚数单位;当负荷采用星型接线时UN选为额定相电压,当采用三角形接线是UN选为额定线电压。
负荷导纳的计算公式为
如图2(a)、图2(b)和图2(c)所示,负荷接线方式有三种情况。由于图2(c)所示的三角形接线方式的等效电路可以转换成图2(b)所示星型中性点不接地接线方式的等效电路,因此,下面只讨论前两种接线方式的模型。
对于Y0接线方式,如图2(a)所示,并联支路线路电压与并联支路电流关系可直接写为
式(4)简写如下
中性点直接接地时的并联支路节点导纳矩阵YP为:
式中,YP是并联支路导纳矩阵,因为中性点直接接地,它也是并联支路节点导纳矩阵,其中三个对角线的元素为每一相支路的导纳;是支路电流矩阵向量;是支路电压矩阵向量。
对于Y接线方式,如图2(b)所示,并联支路线路电压与并联支路电流关系可直接写为:
式中,Un是并联支路中性点对地电压。
因为中性点不接地,三相并联支路电流之和等于零,Un可根据电压以及并联支路阻抗参数求出,表达式如下:
将Un代入式(6)中得到:
式(8)简写为:
式中,Ype为考虑中性点不接地的因素后的并联支路节点导纳矩阵,即:
由式(5)与(10)可知:在负荷中性点接地时,节点导纳矩阵等于支路导纳矩阵;在中性点不接地时,节点导纳矩阵需要根据支路导纳由式(10)求出。
2、故障端口节点阻抗矩计算
以图3所示的5节点放射式配电线路为例,故障发生在节点3和节点4之间,采用端口补偿法计算节点3短路电流,计算故障端口节点3的节点阻抗矩阵。忽略支路15与支路34的线路阻抗,将节点5的负荷阻抗接在节点1上,节点4的负荷阻抗与节点3的负荷阻抗合并,得到故障主回路如图4所示。用线路回路阻抗矩阵与并联支路节点导纳矩阵表示的故障主回路的等效电路如图5所示,其中并联支路的中性点与地之间的阻抗用Zn表示,如果并联支路采用Y0接线时,Zn=0,否则Zn取无穷大;Z0、Z1、Z2、Z3是电源、线路段01、12、23的回路阻抗矩阵,Y0、Y1、Y2、Y3是根节点、节点1、2、3的并联支路节点导纳矩阵。
根据故障主回路等效电流,计算故障端口节点阻抗矩阵步骤为:
1)计算从根节点(即节点0)向电源侧看进去的电源侧节点阻抗矩阵,将电源的中性点视为直接接地,则根节点的电源侧节点阻抗矩阵就是电源回路阻抗矩阵Z0,即:
ZS0=Z0;
2)计算根节点加入并联支路后的节点阻抗矩阵,设根节点电压矩阵向量为流入系统与并联支路的电流矩阵向量为则有
增加并联支路后的根节点的节点阻抗矩阵为:
3)计算节点1的电源侧节点阻抗矩阵,由于节点1电压与根节点电压满足以下关系:
节点1的电源侧节点阻抗矩阵为:
ZS1=Zn0+Z1
4)与根节点的计算方法相同,得到加入后并联支路,节点1的节点阻抗矩阵为
5)同理,可计算出节点2的电源侧电源侧节点阻抗矩阵为
ZS2=Zn1+Z2
6)同理,节点2的节点阻抗矩阵为
7)同理,计算出节点3的电源侧节点矩阵为
ZS3=Zn2+Z3
8)加入节点3的并联支路,得到节点3的节点阻抗矩阵为
其中,Z0为电源的线路回路阻抗矩阵,Z1为节点1到电源之间线路段的线路回路阻抗矩阵,Z2为节点2到节点1之间线路段的线路回路阻抗矩阵,Z3为故障端口节点3与节点2之间线路段的线路回路阻抗矩阵,ZS0为根节点的电源侧节点阻抗矩阵,ZS1为节点1的电源侧节点阻抗矩阵,ZS2为节点2的电源侧节点阻抗矩阵,ZS3为故障端口节点3的电源侧节点阻抗矩阵,Zn0为根节点并联支路的中性点与地之间的阻抗,Zn1为节点1并联支路的中性点与地之间的阻抗,Zn2为节点2并联支路的中性点与地之间的阻抗,Zn3为故障端口节点3并联支路的中性点与地之间的阻抗,Y0为根节点的并联支路节点导纳矩阵,Y1为节点1的并联支路节点导纳矩阵,Y2为节点2的并联支路节点导纳矩阵,Y3为故障端口节点3的并联支路节点导纳矩阵。
本发明分别求解出中性点接地和中性点不接地时的并联支路节点导纳矩阵;考虑到实际中一般负荷阻抗远大于线路阻抗,忽略故障主回路之外的分支线路阻抗(包括故障端口后的线路阻抗),将故障后负荷阻抗叠加到故障点;同时根据故障主回路的配电线路回路阻抗矩阵以及并联支路节点导纳矩阵,通过简单的阻抗串并联计算实现端口节点阻抗矩阵计算,避免了复杂的导纳矩阵求逆运算。
以上所述只是本发明的优选实施方式,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也被视为本发明的保护范围。
Claims (5)
1.一种配电线路故障端口的节点阻抗矩阵计算方法,配电线路故障端口的短路电流采用端口补偿法计算,其特征是,所述节点阻抗矩阵计算方法包括以下过程:
忽略配电线路中并联支路和故障端口之后串联支路的线路阻抗,将并联支路的负荷阻抗添加到对应节点上,将故障端口之后串联支路的负荷阻抗合并到故障端口节点上,得到故障主回路的等效电路;
故障主回路的等效电路采用线路回路阻抗矩阵与并联支路节点导纳矩阵来表示,线路回路阻抗矩阵表示为Z0、Z1、Z2……Zi,各个节点的并联支路节点导纳矩阵表示为Y0、Y1、Y2……Yi,Z0为电源的线路回路阻抗矩阵,Z1为节点1到电源之间线路段的线路回路阻抗矩阵,Zi为故障端口节点i与节点i-1之间线路段的线路回路阻抗矩阵,Y0为根节点的并联支路节点导纳矩阵,Y1为节点1的并联支路节点导纳矩阵,Yi为故障端口节点i的并联支路节点导纳矩阵;
计算各个节点的并联支路节点导纳矩阵;
根据故障主回路的等效电路计算故障端口节点阻抗矩阵;
所述并联支路节点导纳矩阵的计算过程包括以下步骤:
把每一相负荷等效为一个恒阻抗,阻抗值根据负荷的有功功率与无功功率求出,负荷阻抗ZL为:
负荷导纳的计算公式为:
式中,UN为负荷两端的标称电压,PL为负荷的有功功率,QL为负荷额无功功率,j为虚数单位;
当负荷接线方式采用Y0接线方式时,并联支路线路电压与并联支路电流关系可直接写为:
式(3)简写为:
中性点直接接地时的并联支路节点导纳矩阵YP为:
式中,YP为中性点直接接地时并联支路节点导纳矩阵,其中三个对角线的元素为每一相支路的导纳;是支路电流矩阵向量;是支路电压矩阵向量;
当负荷接线方式采用Y接线方式时,并联支路线路电压与并联支路电流关系可直接写为:
式中,Un是并联支路中性点对地电压;
因为中性点不接地,三相并联支路电流之和等于零,Un可根据电压以及并联支路阻抗参数求出,Un的表达式为:
将Un代入式(6)中得到:
式(8)简写为:
中性点不接地时的并联支路节点导纳矩阵Ype为:
式中,Ype为中性点不接地时的并联支路节点导纳矩阵;
故障端口节点阻抗矩阵的计算过程包括以下步骤:
1)计算从根节点向电源侧看进去的电源侧节点阻抗矩阵,将电源的中性点视为直接接地,则根节点的电源侧节点阻抗矩阵就是电源回路阻抗矩阵Z0,即:
ZS0=Z0;
2)计算根节点加入并联支路后的节点阻抗矩阵,设根节点电压矩阵向量为流入系统与并联支路的电流矩阵向量为则有:
增加并联支路后的根节点的节点阻抗矩阵为:
3)计算节点1的电源侧节点阻抗矩阵,由于节点1电压与根节点电压满足以下关系:
节点1的电源侧节点阻抗矩阵为:
ZS1=Zn0+Z1
4)计算节点1加入并联支路后的节点阻抗矩阵,增加并联支路后的节点1的节点阻抗矩阵为:
5)依次计算出节点2至故障端口节点i的电源侧节点阻抗矩阵和节点阻抗矩阵,故障端口节点i的电源侧节点矩阵为:
ZSi=Zn(i-1)+Zi
增加并联支路后的故障端口节点i的节点阻抗矩阵为:
其中,Z0为电源的线路回路阻抗矩阵,Z1为节点1到电源之间线路段的线路回路阻抗矩阵,Zi为故障端口节点i与节点i-1之间线路段的线路回路阻抗矩阵,ZS0为根节点的电源侧节点阻抗矩阵,ZS1为节点1的电源侧节点阻抗矩阵,ZSi为故障端口节点i的电源侧节点阻抗矩阵,Zn0为根节点并联支路的中性点与地之间的阻抗,Zn1为节点1并联支路的中性点与地之间的阻抗,Zni为故障端口节点i并联支路的中性点与地之间的阻抗,Y0为根节点的并联支路节点导纳矩阵,Y1为节点1的并联支路节点导纳矩阵,Yi为故障端口节点i的并联支路节点导纳矩阵。
2.如权利要求1所述的一种配电线路故障端口的节点阻抗矩阵计算方法,其特征是,当负荷采用星型接线时UN为负荷两端的额定相电压,当采用三角形接线时UN为负荷两端的额定线电压。
3.如权利要求1所述的一种配电线路故障端口的节点阻抗矩阵计算方法,其特征是,将负荷接线方式采用三角形接线方式的等效电路转换成星型中性点不接地接线方式的等效电路,三角形接线方式的并联支路节点导纳矩阵采用中性点不接地时的并联支路节点导纳矩阵来求解。
4.如权利要求1所述的一种配电线路故障端口的节点阻抗矩阵计算方法,其特征是,并联支路采用Y0接线时,并联支路的中性点与地之间的阻抗Zn=0,否则Zn取无穷大。
5.如权利要求1至4任意一项所述的一种配电线路故障端口的节点阻抗矩阵计算方法,其特征是,所述并联支路包括无功补偿器并联支路和/或负荷并联支路。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610590130.4A CN105974271B (zh) | 2016-07-25 | 2016-07-25 | 一种配电线路故障端口的节点阻抗矩阵计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610590130.4A CN105974271B (zh) | 2016-07-25 | 2016-07-25 | 一种配电线路故障端口的节点阻抗矩阵计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105974271A CN105974271A (zh) | 2016-09-28 |
CN105974271B true CN105974271B (zh) | 2018-09-25 |
Family
ID=56950739
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610590130.4A Active CN105974271B (zh) | 2016-07-25 | 2016-07-25 | 一种配电线路故障端口的节点阻抗矩阵计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105974271B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107230976A (zh) * | 2016-12-16 | 2017-10-03 | 中国电力科学研究院 | 一种电力系统电气距离可视化方法和装置 |
CN108879666A (zh) * | 2018-07-04 | 2018-11-23 | 国网山东省电力公司济南供电公司 | 一种基于配电网故障的分布式电源接入容量控制方法 |
CN109142967B (zh) * | 2018-07-17 | 2021-12-17 | 中国电力科学研究院有限公司 | 一种用于电网线路的变时限距离保护测距方法和装置 |
CN109599870A (zh) * | 2018-11-14 | 2019-04-09 | 北京电力经济技术研究院有限公司 | 一种抑制短路电流的联络线电抗快速计算方法 |
CN117805555B (zh) * | 2024-02-28 | 2024-05-03 | 昆明理工大学 | 双向推演的四端柔性直流输电线路故障测距方法及系统 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103927459A (zh) * | 2014-05-04 | 2014-07-16 | 华北电力大学(保定) | 一种用于含分布式电源配电网故障定位的方法 |
CN104049175A (zh) * | 2012-08-14 | 2014-09-17 | Abb研究有限公司 | 具有多个多相母线的电网中的故障分析 |
JP2014206532A (ja) * | 2013-04-12 | 2014-10-30 | 三菱電機株式会社 | 非接地配電系統の障害を解析する方法およびそのシステム |
CN104331619A (zh) * | 2014-10-30 | 2015-02-04 | 国网上海市电力公司 | 一种含逆变器接口电源电力系统的短路电流计算方法 |
CN104698346A (zh) * | 2015-03-31 | 2015-06-10 | 国网内蒙古东部电力有限公司通辽供电公司 | 含源配电网故障分析方法及装置 |
CN105548803A (zh) * | 2015-12-04 | 2016-05-04 | 国网山东省电力公司电力科学研究院 | 基于加权最小二乘算法的特高压输电线路故障定位方法 |
CN105759173A (zh) * | 2016-03-31 | 2016-07-13 | 中国科学院电工研究所 | 含分布式电源的配电网的自适应故障区段定位方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2051085A1 (en) * | 2007-10-19 | 2009-04-22 | ABB Research Ltd. | Method for fault location on series compensated power transmission lines with two-end unsychronized measurement |
-
2016
- 2016-07-25 CN CN201610590130.4A patent/CN105974271B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104049175A (zh) * | 2012-08-14 | 2014-09-17 | Abb研究有限公司 | 具有多个多相母线的电网中的故障分析 |
JP2014206532A (ja) * | 2013-04-12 | 2014-10-30 | 三菱電機株式会社 | 非接地配電系統の障害を解析する方法およびそのシステム |
CN103927459A (zh) * | 2014-05-04 | 2014-07-16 | 华北电力大学(保定) | 一种用于含分布式电源配电网故障定位的方法 |
CN104331619A (zh) * | 2014-10-30 | 2015-02-04 | 国网上海市电力公司 | 一种含逆变器接口电源电力系统的短路电流计算方法 |
CN104698346A (zh) * | 2015-03-31 | 2015-06-10 | 国网内蒙古东部电力有限公司通辽供电公司 | 含源配电网故障分析方法及装置 |
CN105548803A (zh) * | 2015-12-04 | 2016-05-04 | 国网山东省电力公司电力科学研究院 | 基于加权最小二乘算法的特高压输电线路故障定位方法 |
CN105759173A (zh) * | 2016-03-31 | 2016-07-13 | 中国科学院电工研究所 | 含分布式电源的配电网的自适应故障区段定位方法 |
Non-Patent Citations (1)
Title |
---|
基于节点阻抗矩阵的配电网故障测距算法;郑涛 等;《电网技术》;20131130;第37卷(第11期);正文第3233至3240页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105974271A (zh) | 2016-09-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105974271B (zh) | 一种配电线路故障端口的节点阻抗矩阵计算方法 | |
CN106549384B (zh) | 一种含upfc电力系统的通用潮流计算方法 | |
Ju et al. | An extension of FBS three-phase power flow for handling PV nodes in active distribution networks | |
CN111061986B (zh) | 一种多运行模式的热电综合能源系统潮流计算方法 | |
CN107681682B (zh) | 一种基于ward等值的交直流系统等值方法 | |
CN103094900B (zh) | 考虑相序混合法的分布式发电配电网三相潮流计算方法 | |
CN103971026B (zh) | 主动配电网潮流计算通用方法 | |
CN105260949B (zh) | 一种基于矩阵变换的配电网络短路计算方法 | |
WO2015106502A1 (zh) | 一种光伏发电站并网的闪变传递系数获取的方法及系统 | |
CN111668843B (zh) | 基于相序混合法的低压配电网三相潮流计算方法 | |
CN110880764B (zh) | 一种含逆变型分布式电源不平衡配电网的故障处理方法 | |
CN103984822A (zh) | 一种基于全数字实时仿真装置的三相潮流实现方法 | |
CN107576886A (zh) | 含逆变型分布式电源小电阻接地系统的接地故障分析方法 | |
CN113514731B (zh) | 一种含逆变型电源不平衡配电网的短路电流确定方法 | |
CN107769213A (zh) | 一种多换流器并联的交直流配电网潮流计算方法 | |
CN104934973B (zh) | 一种考虑温度的电力线路损耗计算方法 | |
JP5933427B2 (ja) | 非接地配電系統の三相電力潮流解析のための方法およびシステム | |
TW200522468A (en) | The power flow model of the multiterminal voltage-source converter-based high voltage DC (M-VSC-HVDC) transmission system for large-scale power systems | |
CN106485395A (zh) | 一种考虑配电网网损和扩容成本的分布式电源配置方法 | |
Teng | Implementation of a voltage sweep power flow method and comparison with other power flow techniques | |
CN106682363B (zh) | 一种类奔德斯分解的下垂控制孤岛微电网潮流计算方法 | |
CN206004320U (zh) | 模块化多电平换流器的仿真装置 | |
CN107171329A (zh) | 一种含新能源并网的配电网潮流计算方法 | |
CN107800137B (zh) | 一种基于mmc的upfc机电暂态仿真方法 | |
CN105608256A (zh) | 一种吸纳不同步长接口延时的并行计算分网方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information |
Address after: 255086 Zibo high tech Development Zone, Shandong province high Chong Park A block, room 313 Applicant after: Shandong University of Technology Address before: 255000 Zibo, Zhangdian District West Village Road, No. 266 Applicant before: Shandong University of Technology |
|
CB02 | Change of applicant information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |