CN108089058A - 一种抗差配电网线路正序参数的辨识方法 - Google Patents
一种抗差配电网线路正序参数的辨识方法 Download PDFInfo
- Publication number
- CN108089058A CN108089058A CN201711329019.0A CN201711329019A CN108089058A CN 108089058 A CN108089058 A CN 108089058A CN 201711329019 A CN201711329019 A CN 201711329019A CN 108089058 A CN108089058 A CN 108089058A
- 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.)
- Pending
Links
Classifications
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种抗差配电网线路正序参数的辨识方法,首先在单回配电网线路两端装设配电网相量量测装置,测量配电网线路两端的三相电压、电流相量;根据所得到的配电网线路两端的三相电压、电流相量,基于线路等效模型建立待辨识的线路阻抗计算方程;根据Huber抗差估计理论及所建立的线路阻抗计算方程,建立基于Huber估计的抗差目标函数,并得到线路阻抗相参数;采用对称分量法将所得到的线路阻抗相参数解耦,得到配电网线路正序阻抗参数。上述方法仅需采集线路两端的电流、电压相量信息,可操作性强,易于实施,且辨识得到的正阻抗参数精度高,具有很好的抗差能力,辨识结果更加可信。
Description
技术领域
本发明涉及电力系统运行控制技术领域,尤其涉及一种抗差配电网线路正序参数的辨识方法。
背景技术
目前,配电网直接与日常生活中的用电负荷和工业、农业、商业的用电设备相连,配电线路在运行中容易老化、受环境腐蚀以及受施工、改造、事故等影响,势必会造成线路阻抗的变化。能够准确了解配电网的参数,有利于主动配电网的构建与控制,对于配电网的继电保护整定计算、故障分析、网损计算、短路电流以及故障定位有着重要的意义,而不恰当的线路参数会使得计算结果与实际情况不一致,从而构成系统潜在的危险或造成不必要的浪费。配电网是电力系统中分配电能的重要一环,其与输电网络有着很大的不同。首先,三相不平衡是配电网的一个重要特征,潮流计算时需要分三相计算;其次,伴随着分布式能源的接入与频繁的倒闸操作,配电网无时无刻不处在变化当中;配电网的大部分参数是未知的,而主网架的参数大部分已知,则主网参数辨识方法对于配网来说并不适用。
在传统的线路理论计算中通常根据Carson模型,利用线路几何均距、材料结构等物理参数,结合气温、地理位置等根据公式计算出电抗、电阻和电纳,或者从电工手册或产品目录中查得单位长度线路的参数在乘以长度得到。理论计算通常只考虑到完全对称的情况,并且无法考虑到实时温度、弧垂等实际存在的问题,必然就导致了该方法的结果和实际参数存在较大的差异。而且由于电力线路参数受运行环境影响易发生阻抗参数变化,则理论计算将会存在较大误差。另外一种方式为离线停电测量,在新建线路投运前或将已运行的线路停电后,采用外加电源,利用电压表、电流表等各种表计测量线路数据,经人工读取表值并结合相应的公式计算出各个参数,该方法存在被测线路须停电、多回运行的平行线路无法测量互感等问题。
为了提高线路阻抗测量精度,在线测量法逐渐被应用,同时随着PMU(相量量测装置)和WAMS(广域量测系统)的大量应用,大量利用SCADA/WAMS提供的数据,采用参数估计实现在线线路参数的辨识,其估计方法主要包含2类方法:增广状态估计法和残差灵敏度分析法。但是鉴于WAMS/PMU主要用在高压电网,无法在配电网应用;另一方面,随着主动配电网的精益化控制需求,符合配电网需求的配电网相量量测装置开发已见报道,但并没有配电网相量量测装置的配电网参数辨识的解决方案。
发明内容
本发明的目的是提供一种抗差配电网线路正序参数的辨识方法,该方法仅需采集线路两端的电流、电压相量信息,可操作性强,易于实施,且辨识得到的正阻抗参数精度高,具有很好的抗差能力,辨识结果更加可信。
本发明的目的是通过以下技术方案实现的:
一种抗差配电网线路正序参数的辨识方法,所述方法包括:
步骤1、在单回配电网线路两端装设配电网相量量测装置,测量配电网线路两端的三相电压、电流相量;
步骤2、根据所得到的配电网线路两端的三相电压、电流相量,基于线路等效模型建立待辨识的线路阻抗计算方程;
步骤3、根据Huber抗差估计理论及所建立的线路阻抗计算方程,建立基于Huber估计的抗差目标函数,并得到线路阻抗相参数;
步骤4、采用对称分量法将所得到的线路阻抗相参数解耦,得到配电网线路正序阻抗参数。
由上述本发明提供的技术方案可以看出,上述方法仅需采集线路两端的电流、电压相量信息,可操作性强,易于实施,且辨识得到的正阻抗参数精度高,具有很好的抗差能力,辨识结果更加可信。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域的普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他附图。
图1为本发明实施例提供的抗差配电网线路正序参数的辨识方法流程示意图;
图2为本发明实施例所举实例中单回配电线路的三相模型示意图;
图3为本发明实施例所举实例中10kV仿真系统示意图。
具体实施方式
下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明的保护范围。
下面将结合附图对本发明实施例作进一步地详细描述,如图1所示为本发明实施例提供的抗差配电网线路正序参数的辨识方法流程示意图,所述方法包括:
步骤1、在单回配电网线路两端装设配电网相量量测装置,测量配电网线路两端的三相电压、电流相量;
步骤2、根据所得到的配电网线路两端的三相电压、电流相量,基于线路等效模型建立待辨识的线路阻抗计算方程;
在该步骤2中,建立待辨识的线路阻抗计算方程的过程具体为:
如图2所示为本发明实施例所举实例中单回配电线路的三相模型示意图,参考图2,所得到的配电网线路两端的三相电压、电流相量满足如下方程:
其中,和为K端的三相电流和电压相量值,K为M或N,表示单回配电网线路的两端;
Z为待辨识的线路阻抗矩阵,进一步表示为:
其中,Zaa、Zbb、Zcc分别为a、b、c三相的线路串联阻抗;Zab为a、b相之间的互阻抗;Zbc为b、c相之间的互阻抗;Zac为a、c相之间的互阻抗;x为待辨识的参数相量;
上述方程进一步表示为:
其中,
由此,在某一时刻待辨识的参数相量
x=[x11 x12 x13 x14 x15 x16 x17 x18 x19]T
满足如下的线性方程:
Ax=B
其中,A为该时刻配电网线路两端的三相电压构成的系数矩阵;B为该时刻配电网线路两端的三相电流组成的常数项。
具体实现,针对所述线性方程Ax=B,若所述配电网相量量测装置的数据存在随机量测误差,则该方程进一步表示为:
Ax=B+v
其中,v为方程残差相量。
步骤3、根据Huber抗差估计理论及所建立的线路阻抗计算方程,建立基于Huber估计的抗差目标函数,并得到线路阻抗相参数;
在该步骤中,Huber抗差估计理论是假设实际观测数据服从Huber分布,Huber分布是污染分布的一种,其主体是正态分布,干扰部分服从Laplace分布;
具体来说,Huber分布的概率密度为:
其中,为标准正态分布密度,在区间-c≤x≤c内,观测值服从正态分布;在x>c时,观测值服从Laplace分布;c的取值在1.0~2.0之间,最常用的c值取1.5;
Huber分布的极大似然估计为Huber估计,其极值函数为:
其中,vi为余差,即vi=(yci-ymi)/s,yci和ymi分别表示输出矢量的计算值和测量值;
在|vi|≤k区间,取量测误差的标准差σ0;
在|vi|>k区间,取Kmad,基于Kmad=med|yci-ymi|。
基于上述分析,所建立的基于Huber估计的抗差目标函数表示为:
其中,i为采样数据点;N为采样总数;vi为第i个量测数据点的残差向量;参数k为可设定的固定值,在本发明实施例中取为标准差σ0。
该基于Huber估计的抗差目标函数在观测误差较小时即为最小二乘方法的目标函数,但在观测误差较大时,能够削弱其对目标函数的不利影响,自动剔除配电网相量量测装置数据中的坏数据。
步骤4、采用对称分量法将所得到的线路阻抗相参数解耦,得到配电网线路正序阻抗参数。
该步骤的具体过程为:
首先构造矩阵T
其中,算子a为ej120°;
由对称分量法可知存在如下关系:
其中,表示三相电压相量或三相电流相量;表示为a相正序分量、负序分量、零序分量电压或电流分量;
进一步的,配电网线路中三相电压降与三相电流存在如下关系:
将其简写为:
再将三相电压降和三相电流替换为相应的序分量,得到:
进一步将其表示为
Zp为序分量的阻抗矩阵,即:
其中,Z(1)即为所要得到的配电网线路正序阻抗参数。
下面再以具体的实例对上述辨识方法进行论证与说明:
本实施例利用PSCAD搭建10kV仿真系统,如图3所示为本发明实施例所举实例中10kV仿真系统示意图,对单回配电线路L1的正序参数进行辨识。10kV线路L1为单回线,线路长度为10km;正序参数设计值为:电阻R1=2.516Ω、电抗XL1=3.588Ω。
假设线路L1两端装设了配电网相量量测装置,测量到了正常运行时线路L1两端的三相不对称电压、电流相量,采样间隔为1ms。现用正常运行期间的1s数据并根据实施方式辨识线路L1的正序阻抗,并设置如下实验,以表明本发明所述方法的有效性。
实验一:直接利用仿真数据,不进行处理;
实验二:在理想仿真数据中叠加随机高斯噪声,模拟真实的配电网相量量测装置数据;其中,电压电流幅值的量测误差标准差为0.1%,相角误差为0.1°;
实验三:在实验一的基础上随机将10组电流幅值量测置零,以模拟配电网相量量测装置量测出现坏数据。
三种实验方案下,按照本发明所述方法得到的正序参数辨识结果如下表1所示:
表1
由表1可知:在理想仿真数据下本发明辨识结果与设计值近乎一致,表明本发明所述方法是基本可行的;实际运行中,从配电网相量量测装置获得的数据含一定量测噪声时所述方法略优于最小二乘法。但若在配电网相量量测装置量测数据存在了坏数据时,所述方法要明显优于最小二乘法,辨识结果的可信度。
对比实验有力地表明本发明所述方法更适用于含量测噪声甚至坏数据的实际配电网相量量测装置量测数据,能有效消弱坏数据对正序阻抗参数辨识的不利影响,所得正序阻抗参数值更可信。
值得注意的是,本发明实施例中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明披露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求书的保护范围为准。
Claims (5)
1.一种抗差配电网线路正序参数的辨识方法,其特征在于,所述方法包括:
步骤1、在单回配电网线路两端装设配电网相量量测装置,测量配电网线路两端的三相电压、电流相量;
步骤2、根据所得到的配电网线路两端的三相电压、电流相量,基于线路等效模型建立待辨识的线路阻抗计算方程;
步骤3、根据Huber抗差估计理论及所建立的线路阻抗计算方程,建立基于Huber估计的抗差目标函数,并得到线路阻抗相参数;
步骤4、采用对称分量法将所得到的线路阻抗相参数解耦,得到配电网线路正序阻抗参数。
2.根据权利要求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>
其中,和为K端的三相电流和电压相量值,K为M或N,表示单回配电网线路的两端;
Z为待辨识的线路阻抗矩阵,进一步表示为:
<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相之间的互阻抗;x为待辨识的参数相量;
上述方程进一步表示为:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>Z</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<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>
<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>
<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>
<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为该时刻配电网线路两端的三相电压构成的系数矩阵;B为该时刻配电网线路两端的三相电流组成的常数项。
3.根据权利要求2所述抗差配电网线路正序参数的辨识方法,其特征在于,针对所述线性方程Ax=B,若所述配电网相量量测装置的数据存在随机量测误差,则该方程进一步表示为:
Ax=B+v
其中,v为方程残差相量。
4.根据权利要求3所述抗差配电网线路正序参数的辨识方法,其特征在于,在所述步骤3中,所建立的基于Huber估计的抗差目标函数表示为:
<mrow>
<mi>E</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>N</mi>
</munderover>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<msup>
<msub>
<mi>v</mi>
<mi>i</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
<mn>2</mn>
</mfrac>
<msub>
<mo>|</mo>
<mrow>
<mo>|</mo>
<msub>
<mi>v</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<mo>&le;</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>|</mo>
<msub>
<mi>v</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<mo>-</mo>
<mfrac>
<msup>
<mi>k</mi>
<mn>2</mn>
</msup>
<mn>2</mn>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mo>|</mo>
<mrow>
<mo>|</mo>
<msub>
<mi>v</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<mo>></mo>
<mi>k</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
其中,i为采样数据点;N为采样总数;vi为第i个量测数据点的残差向量;参数k为可设定的固定值。
5.根据权利要求1所述抗差配电网线路正序参数的辨识方法,其特征在于,所述步骤4的具体过程为:
首先构造矩阵T
<mrow>
<mi>T</mi>
<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>
</mrow>
其中,算子a为ej120°;
由对称分量法可知存在如下关系:
<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相正序分量、负序分量、零序分量电压或电流分量;
进一步的,配电网线路中三相电压降与三相电流存在如下关系:
<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>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>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mrow>
进一步将其表示为
<mrow>
<mi>&Delta;</mi>
<msub>
<mover>
<mi>U</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
<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>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
<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>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
<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 |
---|---|---|---|
CN201711329019.0A CN108089058A (zh) | 2017-12-13 | 2017-12-13 | 一种抗差配电网线路正序参数的辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711329019.0A CN108089058A (zh) | 2017-12-13 | 2017-12-13 | 一种抗差配电网线路正序参数的辨识方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108089058A true CN108089058A (zh) | 2018-05-29 |
Family
ID=62175458
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711329019.0A Pending CN108089058A (zh) | 2017-12-13 | 2017-12-13 | 一种抗差配电网线路正序参数的辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108089058A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112713589A (zh) * | 2020-12-18 | 2021-04-27 | 国网辽宁省电力有限公司 | 一种基于电能质量分析仪测量的低压配电网线路阻抗估计方法 |
CN116805799A (zh) * | 2023-08-22 | 2023-09-26 | 国网福建省电力有限公司 | 一种考虑电压约束和线路参数的分布式光伏消纳方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101419253A (zh) * | 2008-12-17 | 2009-04-29 | 河南电力试验研究院 | 一种特高压输电线路正序、零序参数测量方法及系统 |
CN107248736A (zh) * | 2017-05-25 | 2017-10-13 | 华北电力大学 | 一种配电网线路正序参数的在线辨识方法 |
CN107328996A (zh) * | 2017-06-12 | 2017-11-07 | 华北电力大学 | 抗差输电线路正序参数辨识方法 |
-
2017
- 2017-12-13 CN CN201711329019.0A patent/CN108089058A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101419253A (zh) * | 2008-12-17 | 2009-04-29 | 河南电力试验研究院 | 一种特高压输电线路正序、零序参数测量方法及系统 |
CN107248736A (zh) * | 2017-05-25 | 2017-10-13 | 华北电力大学 | 一种配电网线路正序参数的在线辨识方法 |
CN107328996A (zh) * | 2017-06-12 | 2017-11-07 | 华北电力大学 | 抗差输电线路正序参数辨识方法 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112713589A (zh) * | 2020-12-18 | 2021-04-27 | 国网辽宁省电力有限公司 | 一种基于电能质量分析仪测量的低压配电网线路阻抗估计方法 |
CN116805799A (zh) * | 2023-08-22 | 2023-09-26 | 国网福建省电力有限公司 | 一种考虑电压约束和线路参数的分布式光伏消纳方法 |
CN116805799B (zh) * | 2023-08-22 | 2023-10-27 | 国网福建省电力有限公司 | 一种考虑电压约束和线路参数的分布式光伏消纳方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100526895C (zh) | 抗分布电容电流和过渡电阻影响的线路单端故障测距方法 | |
Al-Mohammed et al. | A fully adaptive PMU-based fault location algorithm for series-compensated lines | |
CN102435851B (zh) | 一种双回输电线路零序参数测量方法 | |
CN107248736A (zh) | 一种配电网线路正序参数的在线辨识方法 | |
CN107328996B (zh) | 抗差输电线路正序参数辨识方法 | |
CN104836223A (zh) | 电网参数错误与不良数据协同辨识与估计方法 | |
CN103838959A (zh) | 偏最小二乘回归应用于配电网谐波源定位与检测中的方法 | |
CN106054023B (zh) | 一种输电线路单端测距中估计两侧系统阻抗的方法 | |
CN103869184A (zh) | 基于pmu数据的单回输电线路零序参数抗差辨识方法 | |
CN109283430A (zh) | 一种基于电压分布原理的配电网故障测距方法 | |
CN109523165A (zh) | 一种基于有限信息的电压暂降状态估计方法 | |
CN105044649B (zh) | 一种电流互感器误差水平动态检测方法及系统 | |
CN104237677B (zh) | 220kV及以上电网单回输电线路零序参数在线辨识方法 | |
Negash et al. | Optimal placement of phasor measurement unit for system observability:(case study: southern region power system network, Ethiopia) | |
CN110333394A (zh) | 一种低压配电网线路阻抗估计方法 | |
CN102590686B (zh) | 一种双极直流输电线路无线电干扰确定方法 | |
CN108089058A (zh) | 一种抗差配电网线路正序参数的辨识方法 | |
CN107797017B (zh) | 一种电力变压器带电检测损耗特征参数的方法 | |
CN103972889B (zh) | 一种配电线路阻抗在线辨识方法 | |
CN105203877B (zh) | 能够消除量测误差影响的单回输电线路零序参数辨识方法 | |
Li et al. | Augmented state estimation method for fault location based on on-line parameter identification of PMU measurement data | |
Skok et al. | Hybrid state estimation model based on PMU and SCADA measurements | |
CN108052725A (zh) | 一种基于指数型目标函数的抗差配电网线路正序参数辨识方法 | |
CN109470988A (zh) | 一种适用于t型输电线路的故障测距系统及方法 | |
CN109241626A (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180529 |
|
RJ01 | Rejection of invention patent application after publication |