CN107248736A - 一种配电网线路正序参数的在线辨识方法 - Google Patents
一种配电网线路正序参数的在线辨识方法 Download PDFInfo
- Publication number
- CN107248736A CN107248736A CN201710379619.1A CN201710379619A CN107248736A CN 107248736 A CN107248736 A CN 107248736A CN 201710379619 A CN201710379619 A CN 201710379619A CN 107248736 A CN107248736 A CN 107248736A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- mtr
- mover
- 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
- 238000000034 method Methods 0.000 title claims abstract description 45
- 239000011159 matrix material Substances 0.000 claims abstract description 23
- 238000013178 mathematical model Methods 0.000 claims description 9
- 238000005259 measurement Methods 0.000 abstract description 13
- 238000004364 calculation method Methods 0.000 description 7
- 238000004088 simulation Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 101100499229 Mus musculus Dhrsx gene Proteins 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002847 impedance measurement Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 239000002699 waste material Substances 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R27/00—Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
- G01R27/02—Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant
- G01R27/08—Measuring resistance by measuring both voltage and current
-
- 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]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种配电网线路正序参数的在线辨识方法,首先根据单回配电线路三相等值模型,得到单回配电线路相分量模型,并以单回配电线路两端的电流电压分量建立线路阻抗参数辨识数学模型;在单回配电线路两端装设配电网相量量测装置,测量所述单回配电线路两端的三相电压、电流相量;根据所述单回配电线路两端的三相电压、电流相量和所建立的线路阻抗参数辨识数学模型,利用最小二乘法计算待辨识的线路阻抗矩阵;再采用对称分量法将待辨识的线路阻抗矩阵解耦,得到配电网线路正序阻抗参数。上述方法仅需采集线路两端的电流、电压相量信息,可操作性强、易于实施,且辨识得到的正阻抗参数精度高、更加可信。
Description
技术领域
本发明涉及电力系统控制技术领域,尤其涉及一种配电网线路正序参数的在线辨识方法。
背景技术
目前,配电网直接与日常生活中的用电负荷和工业、农业、商业的用电设备相连,配电线路在运行中容易老化、受环境腐蚀以及受施工、改造、事故等影响,势必会造成线路阻抗的变化。能够准确了解配电网的参数,有利于主动配电网的构建与控制,对于配电网的继电保护整定计算、故障分析、网损计算、短路电流以及故障定位有着重要的意义,而不恰当的线路参数会使得计算结果与实际情况不一致,从而构成系统潜在的危险或造成不必要的浪费。配电网是电力系统中分配电能的重要一环,其与输电网络有着很大的不同,首先三相不平衡是配电网的一个重要特征,潮流计算时需要分三相计算;其次伴随着分布式能源的接入与频繁的倒闸操作,配电网无时无刻不处在变化当中;配电网的大部分参数是未知的,而主网架的参数大部分已知,则主网参数辨识方法对于配网来说并不适用。
在传统的线路理论计算中通常根据Carson模型,利用线路几何均距、材料结构等物理参数,结合气温、地理位置等根据公式计算出电抗、电阻和电纳,或者从电工手册或产品目录中查得单位长度线路的参数在乘以长度得到。理论计算通常只考虑到完全对称的情况,并且无法考虑到实时温度、弧垂等实际存在的问题,必然就导致了该方法的结果和实际参数存在较大的差异,而且由于电力线路参数受运行环境影响易发生阻抗参数变化,则理论计算将会存在较大误差。另外一种方式为离线停电测量,在新建线路投运前或将已运行的线路停电后,采用外加电源,利用电压表、电流表等各种表计测量线路数据,经人工读取表值并结合相应的公式计算出各个参数,但该方法存在被测线路须停电、多回运行的平行线路无法测量互感等问题。
为了提高线路阻抗测量精度,在线测量法逐渐被应用,随着PMU(相量量测装置)和WAMS(广域量测系统)的大量应用,大量利用SCADA/WAMS提供的数据,采用参数估计实现在线线路参数的辨识,其估计方法主要包含增广状态估计法和残差灵敏度分析法。但是鉴于WAMS/PMU主要用在高压电网,无法在配电网应用,同时随着主动配电网的精益化控制需求,现有技术中还没有配电网相量量测装置的配电网参数辨识技术。
发明内容
本发明的目的是提供一种配电网线路正序参数的在线辨识方法,该方法仅需采集线路两端的电流、电压相量信息,可操作性强、易于实施,且辨识得到的正阻抗参数精度高、更加可信。
一种配电网线路正序参数的在线辨识方法,所述方法包括:
步骤1、根据单回配电线路三相等值模型,得到单回配电线路相分量模型,并以单回配电线路两端的电流电压分量建立线路阻抗参数辨识数学模型;
步骤2、在单回配电线路两端装设配电网相量量测装置,测量所述单回配电线路两端的三相电压、电流相量;
步骤3、根据所述单回配电线路两端的三相电压、电流相量和所建立的线路阻抗参数辨识数学模型,利用最小二乘法计算待辨识的线路阻抗矩阵;
步骤4、采用对称分量法将待辨识的线路阻抗矩阵解耦,得到配电网线路正序阻抗参数。
在所述步骤1中:
所得到单回配电线路相分量模型为:
其中,表示单回配电线路两端的a、b、c三相的电流相量;表示单回配电线路两端的a、b、c三相的电压相量;K为M或N,表示单回配电线路的两端;
进一步得到线路阻抗参数辨识数学模型表示为:
其中,Zaa、Zbb、Zcc分别为a、b、c三相的线路串联阻抗;
Zab为a、b相之间的互阻抗;Zbc为b、c相之间的互阻抗;Zac为a、c相之间的互阻抗。
在所述步骤3中:
将单回配电线路相分量模型进一步表示为:
其中,
由此可得,待辨识的线路阻抗矩阵
X=[X11 X12 X13 X14 X15 X16 X17 X18 X19]T
并满足如下线性方程:
AX=B
其中,A为某时刻单回配电线路两端的a、b、c三相的电压相量构成的系数矩阵,B为某时刻单回配电线路两端的a、b、c三相的电流相量组成的常数项,X为待辨识的线路阻抗矩阵;
由此可得,待辨识的线路阻抗矩阵:
X=A-1B。
在所述步骤4中,由对称分量法可知:
其中,算子a为ej120°,表示三相电压相量或是三相电流相量,表示为a相正序分量、负序分量、零序分量电压或电流分量;
将上式记为:
F=T-1Fp
单回配电线路中三相电压降与三相电流存在如下关系:
将上式简写为:
再将三相电压降和三相电流替换为序分量,得到:
进而得到:
由此,序分量的阻抗矩阵Zp表示为:
其中,Z(1)为配电网线路正序阻抗参数。
由上述本发明提供的技术方案可以看出,上述方法仅需采集线路两端的电流、电压相量信息,可操作性强、易于实施,且辨识得到的正阻抗参数精度高、更加可信。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域的普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他附图。
图1为本发明实施例所提供配电网线路正序参数的在线辨识方法流程示意图;
图2为本发明实施例所述单回配电线路的三相模型示意图;
图3为本发明实施例所举实例中10kV仿真系统示意图。
具体实施方式
下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明的保护范围。
下面将结合附图对本发明实施例作进一步地详细描述,如图1所示为本发明实施例所提供配电网线路正序参数的在线辨识方法流程示意图,所述方法包括:
步骤1、根据单回配电线路三相等值模型,得到单回配电线路相分量模型,并以单回配电线路两端的电流电压分量建立线路阻抗参数辨识数学模型;
在所述步骤1中,所得到单回配电线路相分量模型为:
如图2所示为本发明实施例所述单回配电线路的三相模型示意图,参考图2:
表示单回配电线路两端的a、b、c三相的电流相量;表示单回配电线路两端的a、b、c三相的电压相量;K为M或N,表示单回配电线路的两端。
进一步得到线路阻抗参数辨识数学模型表示为:
其中,Zaa、Zbb、Zcc分别为a、b、c三相的线路串联阻抗;
Zab为a、b相之间的互阻抗;Zbc为b、c相之间的互阻抗;Zac为a、c相之间的互阻抗。
步骤2、在单回配电线路两端装设配电网相量量测装置,测量所述单回配电线路两端的三相电压、电流相量;
步骤3、根据所述单回配电线路两端的三相电压、电流相量和所建立的线路阻抗参数辨识数学模型,利用最小二乘法计算待辨识的线路阻抗矩阵;
在所述步骤3中,将单回配电线路相分量模型进一步表示为:
其中,
由此可得,待辨识的线路阻抗矩阵
X=[X11 X12 X13 X14 X15 X16 X17 X18 X19]T
并满足如下线性方程:
AX=B
其中,A为某时刻单回配电线路两端的a、b、c三相的电压相量构成的系数矩阵,B为某时刻单回配电线路两端的a、b、c三相的电流相量组成的常数项,X为待辨识的线路阻抗矩阵;具体表示为:
由此可得,待辨识的线路阻抗矩阵:
X=A-1B。
步骤4、采用对称分量法将待辨识的线路阻抗矩阵解耦,得到配电网线路正序阻抗参数。
在所述步骤4中,由对称分量法可知:
其中,算子a为ej120°,表示三相电压相量或是三相电流相量,表示为a相正序分量、负序分量、零序分量电压或电流分量;
上式记为:
F=T-1Fp
单回配电线路中三相电压降与三相电流存在如下关系:
上式简写为:
再将三相电压降和三相电流替换为序分量,得到:
进而得到
由此,序分量的阻抗矩阵Zp表示为:
其中,Z(1)为配电网线路正序阻抗参数。
下面再以具体的实例对本发明所述方法进行论证,如图3所示为本发明实施例所举实例中10kV仿真系统示意图,利用PSCAD搭建10kV仿真系统,对单回配电线路L1的正序参数进行辨识,10kV线路L1为单回线,线路长度为10km;正序参数设计值为:
电阻R1=2.732Ω、电抗XL1=3.364Ω。
假设单回配电线路L1两端装设了配电网相量量测装置,并设置如下实验,以表明本发明所述方法的有效性:
实验一:直接利用仿真数据,不进行处理;
实验二:在理想仿真数据中叠加随机高斯噪声,模拟真实的配电网相量量测装置数据;其中,电压电流幅值的量测误差标准差为0.1%,相角误差为0.1°;
两种实验方案下,本发明所述方法得到的正序参数辨识结果如下表1所示:
表1
上表1表明:在理想仿真数据下本发明辨识结果与设计值近乎一致,表明本发明所述方法是可行的;实际运行中,从配电网相量量测装置获得的数据含一定量测噪声,实验二通过在理想仿真数据基础上叠加随机高斯噪声来模拟实际运行工况,实验结果表明在添加噪声之后,本发明方法辨识结果仍具有较高精度。
如上述实验结果所示本发明所述方法适用于含量测噪声实际配电网相量量测装置的量测数据,能有效削弱噪声对正序参数辨识的不利影响,所得正序参数值可信度高。
综上所述,本发明实施例所提供的方法仅需采集线路两端的电流、电压相量信息,可操作性强、易于实施,且辨识得到的正阻抗参数精度高、更加可信;同时该方法不受线路所带负荷影响,亦不受周围环境,线路位置等影响,受外界因素干扰少,在外因影响下也能实现准确计算,实用性高。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明披露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求书的保护范围为准。
Claims (4)
1.一种配电网线路正序参数的在线辨识方法,其特征在于,所述方法包括:
步骤1、根据单回配电线路三相等值模型,得到单回配电线路相分量模型,并以单回配电线路两端的电流电压分量建立线路阻抗参数辨识数学模型;
步骤2、在单回配电线路两端装设配电网相量量测装置,测量所述单回配电线路两端的三相电压、电流相量;
步骤3、根据所述单回配电线路两端的三相电压、电流相量和所建立的线路阻抗参数辨识数学模型,利用最小二乘法计算待辨识的线路阻抗矩阵;
步骤4、采用对称分量法将待辨识的线路阻抗矩阵解耦,得到配电网线路正序阻抗参数。
2.根据权利要求1所述配电网线路正序参数的在线辨识方法,其特征在于,在所述步骤1中:
所得到单回配电线路相分量模型为:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>Z</mi>
<msub>
<mover>
<mi>I</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>M</mi>
<mi>a</mi>
<mi>b</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mover>
<mi>U</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>M</mi>
<mi>a</mi>
<mi>b</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>U</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>N</mi>
<mi>a</mi>
<mi>b</mi>
<mi>c</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>Z</mi>
<msub>
<mover>
<mi>I</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>N</mi>
<mi>a</mi>
<mi>b</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mover>
<mi>U</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>N</mi>
<mi>a</mi>
<mi>b</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>U</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>M</mi>
<mi>a</mi>
<mi>b</mi>
<mi>c</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,表示单回配电线路两端的a、b、c三相的电流相量;表示单回配电线路两端的a、b、c三相的电压相量;K为M或N,表示单回配电线路的两端;
进一步得到线路阻抗参数辨识数学模型表示为:
<mrow>
<mi>Z</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>Z</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>Z</mi>
<mn>12</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>Z</mi>
<mn>13</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>Z</mi>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>Z</mi>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>Z</mi>
<mn>23</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>Z</mi>
<mn>31</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>Z</mi>
<mn>32</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>Z</mi>
<mn>33</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mi>a</mi>
<mi>a</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mi>a</mi>
<mi>c</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mi>b</mi>
<mi>b</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mi>b</mi>
<mi>c</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mi>a</mi>
<mi>c</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mi>b</mi>
<mi>c</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mi>c</mi>
<mi>c</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mn>11</mn>
</msub>
<mo>+</mo>
<msub>
<mi>jX</mi>
<mn>11</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mn>12</mn>
</msub>
<mo>+</mo>
<msub>
<mi>jX</mi>
<mn>12</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mn>13</mn>
</msub>
<mo>+</mo>
<msub>
<mi>jX</mi>
<mn>13</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mn>21</mn>
</msub>
<mo>+</mo>
<msub>
<mi>jX</mi>
<mn>21</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mn>22</mn>
</msub>
<mo>+</mo>
<msub>
<mi>jX</mi>
<mn>22</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mn>23</mn>
</msub>
<mo>+</mo>
<msub>
<mi>jX</mi>
<mn>23</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mn>31</mn>
</msub>
<mo>+</mo>
<msub>
<mi>jX</mi>
<mn>31</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mn>32</mn>
</msub>
<mo>+</mo>
<msub>
<mi>jX</mi>
<mn>32</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mn>33</mn>
</msub>
<mo>+</mo>
<msub>
<mi>jX</mi>
<mn>33</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,Zaa、Zbb、Zcc分别为a、b、c三相的线路串联阻抗;
Zab为a、b相之间的互阻抗;Zbc为b、c相之间的互阻抗;Zac为a、c相之间的互阻抗。
3.根据权利要求2所述配电网线路正序参数的在线辨识方法,其特征在于,在所述步骤3中:
将单回配电线路相分量模型进一步表示为:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>Z</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mover>
<mi>U</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>M</mi>
<mi>a</mi>
<mi>b</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>U</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>N</mi>
<mi>a</mi>
<mi>b</mi>
<mi>c</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mover>
<mi>I</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>M</mi>
<mi>a</mi>
<mi>b</mi>
<mi>c</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>Z</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mover>
<mi>U</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>N</mi>
<mi>a</mi>
<mi>b</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>U</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>M</mi>
<mi>a</mi>
<mi>b</mi>
<mi>c</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mover>
<mi>I</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>N</mi>
<mi>a</mi>
<mi>b</mi>
<mi>c</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,
由此可得,待辨识的线路阻抗矩阵
X=[X11 X12 X13 X14 X15 X16 X17 X18 X19]T
并满足如下线性方程:
AX=B
其中,A为某时刻单回配电线路两端的a、b、c三相的电压相量构成的系数矩阵,B为某时刻单回配电线路两端的a、b、c三相的电流相量组成的常数项,X为待辨识的线路阻抗矩阵;
由此可得,待辨识的线路阻抗矩阵:
X=A-1B。
4.根据权利要求3所述配电网线路正序参数的在线辨识方法,其特征在于,在所述步骤4中,由对称分量法可知:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mover>
<mi>F</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>a</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>F</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>b</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>F</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>c</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mi>a</mi>
</mtd>
<mtd>
<msup>
<mi>a</mi>
<mn>2</mn>
</msup>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<msup>
<mi>a</mi>
<mn>2</mn>
</msup>
</mtd>
<mtd>
<mi>a</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>a</mi>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mover>
<mi>F</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>a</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>F</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>a</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>F</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>a</mi>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,算子a为ej120°,表示三相电压相量或是三相电流相量,表示为a相正序分量、负序分量、零序分量电压或电流分量;
将上式记为:
F=T-1Fp
单回配电线路中三相电压降与三相电流存在如下关系:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>U</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>M</mi>
<mi>a</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>U</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>N</mi>
<mi>a</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>U</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>M</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>U</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>N</mi>
<mi>b</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>U</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>M</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>U</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>N</mi>
<mi>c</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mi>a</mi>
<mi>a</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mi>a</mi>
<mi>c</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mi>b</mi>
<mi>b</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mi>b</mi>
<mi>c</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mi>a</mi>
<mi>c</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mi>b</mi>
<mi>c</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mi>c</mi>
<mi>c</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mover>
<mi>I</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>M</mi>
<mi>a</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>I</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>M</mi>
<mi>b</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>I</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>M</mi>
<mi>c</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
将上式简写为:
<mrow>
<mi>&Delta;</mi>
<msub>
<mover>
<mi>U</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>a</mi>
<mi>b</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>Z</mi>
<msub>
<mover>
<mi>I</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>a</mi>
<mi>b</mi>
<mi>c</mi>
</mrow>
</msub>
</mrow>
再将三相电压降和三相电流替换为序分量,得到:
<mrow>
<mi>T</mi>
<mi>&Delta;</mi>
<msub>
<mover>
<mi>U</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
<mo>)</mo>
</mrow>
</msub>
<mo>=</mo>
<mi>Z</mi>
<mi>T</mi>
<msub>
<mover>
<mi>I</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>M</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mrow>
进而得到:
<mrow>
<mi>&Delta;</mi>
<msub>
<mover>
<mi>U</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
<mo>)</mo>
</mrow>
</msub>
<mo>=</mo>
<msup>
<mi>T</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>Z</mi>
<mi>T</mi>
<msub>
<mover>
<mi>I</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>M</mi>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>Z</mi>
<mi>p</mi>
</msub>
<msub>
<mover>
<mi>I</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>M</mi>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mrow>
由此,序分量的阻抗矩阵Zp表示为:
<mrow>
<msub>
<mi>Z</mi>
<mi>p</mi>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msub>
</mtd>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mo>*</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</msub>
</mtd>
<mtd>
<mo>*</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,Z(1)为配电网线路正序阻抗参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710379619.1A CN107248736B (zh) | 2017-05-25 | 2017-05-25 | 一种配电网线路正序参数的在线辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710379619.1A CN107248736B (zh) | 2017-05-25 | 2017-05-25 | 一种配电网线路正序参数的在线辨识方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107248736A true CN107248736A (zh) | 2017-10-13 |
CN107248736B CN107248736B (zh) | 2020-12-29 |
Family
ID=60017158
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710379619.1A Expired - Fee Related CN107248736B (zh) | 2017-05-25 | 2017-05-25 | 一种配电网线路正序参数的在线辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107248736B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108052725A (zh) * | 2017-12-08 | 2018-05-18 | 华北电力大学 | 一种基于指数型目标函数的抗差配电网线路正序参数辨识方法 |
CN108051700A (zh) * | 2017-10-19 | 2018-05-18 | 北京交通大学 | 基于μPMU的配电线路参数辨识的相分量故障测距方法 |
CN108089058A (zh) * | 2017-12-13 | 2018-05-29 | 华北电力大学 | 一种抗差配电网线路正序参数的辨识方法 |
CN109462245A (zh) * | 2018-12-29 | 2019-03-12 | 西安交通大学 | 基于最小二乘法的局域电网负序不平衡综合治理方法 |
CN111242459A (zh) * | 2020-01-07 | 2020-06-05 | 中国南方电网有限责任公司 | 全网设备参数异常值辨识方法及系统 |
CN111475929A (zh) * | 2020-03-20 | 2020-07-31 | 广西电网有限责任公司电力科学研究院 | 基于配电网实境试验平台监测数据的反演验证方法及系统 |
CN113285453A (zh) * | 2021-06-11 | 2021-08-20 | 燕山大学 | 一种配电网线路参数反向识别方法和系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11230999A (ja) * | 1998-02-17 | 1999-08-27 | Chubu Electric Power Co Inc | 高調波特性測定用計測装置 |
CN101699301A (zh) * | 2009-10-29 | 2010-04-28 | 江苏省电力试验研究院有限公司 | 感应电测量线路参数的方法 |
-
2017
- 2017-05-25 CN CN201710379619.1A patent/CN107248736B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11230999A (ja) * | 1998-02-17 | 1999-08-27 | Chubu Electric Power Co Inc | 高調波特性測定用計測装置 |
CN101699301A (zh) * | 2009-10-29 | 2010-04-28 | 江苏省电力试验研究院有限公司 | 感应电测量线路参数的方法 |
Non-Patent Citations (1)
Title |
---|
赵进全等: ""一种非对称输电线路参数的在线测量方法"", 《西安交通大学学报》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108051700A (zh) * | 2017-10-19 | 2018-05-18 | 北京交通大学 | 基于μPMU的配电线路参数辨识的相分量故障测距方法 |
CN108051700B (zh) * | 2017-10-19 | 2019-11-08 | 北京交通大学 | 基于μPMU的配电线路参数辨识的相分量故障测距方法 |
CN108052725A (zh) * | 2017-12-08 | 2018-05-18 | 华北电力大学 | 一种基于指数型目标函数的抗差配电网线路正序参数辨识方法 |
CN108089058A (zh) * | 2017-12-13 | 2018-05-29 | 华北电力大学 | 一种抗差配电网线路正序参数的辨识方法 |
CN109462245A (zh) * | 2018-12-29 | 2019-03-12 | 西安交通大学 | 基于最小二乘法的局域电网负序不平衡综合治理方法 |
CN109462245B (zh) * | 2018-12-29 | 2021-01-19 | 西安交通大学 | 基于最小二乘法的局域电网负序不平衡综合治理方法 |
CN111242459A (zh) * | 2020-01-07 | 2020-06-05 | 中国南方电网有限责任公司 | 全网设备参数异常值辨识方法及系统 |
CN111475929A (zh) * | 2020-03-20 | 2020-07-31 | 广西电网有限责任公司电力科学研究院 | 基于配电网实境试验平台监测数据的反演验证方法及系统 |
CN111475929B (zh) * | 2020-03-20 | 2022-12-20 | 广西电网有限责任公司电力科学研究院 | 基于配电网实境试验平台监测数据的反演验证方法及系统 |
CN113285453A (zh) * | 2021-06-11 | 2021-08-20 | 燕山大学 | 一种配电网线路参数反向识别方法和系统 |
CN113285453B (zh) * | 2021-06-11 | 2022-10-25 | 燕山大学 | 一种配电网线路参数反向识别方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN107248736B (zh) | 2020-12-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107248736B (zh) | 一种配电网线路正序参数的在线辨识方法 | |
CN102435851B (zh) | 一种双回输电线路零序参数测量方法 | |
CN100526895C (zh) | 抗分布电容电流和过渡电阻影响的线路单端故障测距方法 | |
CN105067950B (zh) | 基于纵向阻抗的双端量故障测距方法 | |
Al-Mohammed et al. | A fully adaptive PMU-based fault location algorithm for series-compensated lines | |
CN107328996B (zh) | 抗差输电线路正序参数辨识方法 | |
CN109031000A (zh) | 一种基于非故障扰动就地测量电网短路容量的方法及系统 | |
CN107064736B (zh) | 一种含多t接逆变型分布式电源配电网的故障定位方法 | |
CN105720594B (zh) | 纯有功三相不平衡负荷补偿容量的补偿方法 | |
CN101505061A (zh) | 一种基于时域仿真的戴维南等值参数跟踪的计算方法 | |
CN104237677B (zh) | 220kV及以上电网单回输电线路零序参数在线辨识方法 | |
CN103869184A (zh) | 基于pmu数据的单回输电线路零序参数抗差辨识方法 | |
CN105375484A (zh) | 一种基于pmu的电力系统分布式动态状态估计方法 | |
CN103869171A (zh) | 超高压同塔四回交流/双回双极直流线路零序参数测量法 | |
CN107609802A (zh) | 一种含多馈入直流的交直流混联电力系统功率/电压稳定裕度评估方法 | |
CN106054023A (zh) | 一种输电线路单端测距中估计两侧系统阻抗的方法 | |
Skok et al. | Hybrid state estimation model based on PMU and SCADA measurements | |
CN103972889B (zh) | 一种配电线路阻抗在线辨识方法 | |
CN108075469B (zh) | 一种配电变压器零序阻抗计算方法 | |
CN108089058A (zh) | 一种抗差配电网线路正序参数的辨识方法 | |
CN106159941B (zh) | 一种考虑实际量测误差传递特性的电力系统状态估计方法 | |
Vicol | On-line overhead transmission line And transformer parameters identification based on PMU measurements | |
CN105203877B (zh) | 能够消除量测误差影响的单回输电线路零序参数辨识方法 | |
Goklani et al. | Quantification of minimum unbalance required for accurate estimation of sequence parameters of transmission line using PMU data | |
CN103217585A (zh) | 一种t型输电线路正序电容测量方法 |
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: 20201229 |