CN102967779A - 一种输电线路分布参数的辨识方法 - Google Patents
一种输电线路分布参数的辨识方法 Download PDFInfo
- Publication number
- CN102967779A CN102967779A CN2012104588244A CN201210458824A CN102967779A CN 102967779 A CN102967779 A CN 102967779A CN 2012104588244 A CN2012104588244 A CN 2012104588244A CN 201210458824 A CN201210458824 A CN 201210458824A CN 102967779 A CN102967779 A CN 102967779A
- Authority
- CN
- China
- Prior art keywords
- value
- partiald
- dtri
- transmission line
- prime
- 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
Images
Landscapes
- Emergency Protection Circuit Devices (AREA)
Abstract
一种输电线路分布参数的辨识方法,它包括有如下的步骤:对于发生外部故障输电线路,采集首末两端的电信号瞬时值、计算出正序分量和零序分量;计算正序分布参数,即电阻R1,电感L1,电容C1;再计算零序分布参数,即电阻R0,电感L0,电容C0。本发明只需要采集两端电压、电流同步瞬时值即可求出线路分布参数,无需知道两侧参数及其余部分的运行状态,且该方法辨识精度不受故障条件的影响,具有操作简便和计算准确的优点,可为电网提供准确的线路参数,保证后续电力***计算如状态估计、潮流计算、故障测距、继电保护整定等计算结果的可靠性。
Description
技术领域
本发明涉及一种输电线路参数的辨识方法,特别是一种基于故障数据的输电线路分布参数的辨识方法。
背景技术
随着我国电网建设的高速发展,电网规模也逐步扩大,分布地域更广、元件更多,网络结构更复杂,电网稳定安全运行面临严峻挑战。而准确的电网模型参数是状态估计、潮流计算、故障测距、继电保护整定等电力***计算的重要基础,提高其准确性对于线路和***的安全运行具有重要意义。但由于地形和气候等各种因素的影响将导致线路参数的改变,电力部门掌握的线路参数往往存在一些不精确,因此,对于可疑线路,有必要根据已有量测信息重新估计其参数。并且鉴于目前全国各地特高压输电迅速发展,对于电压等级高、输电距离长的线路来说,分析相关问题时必须要考虑其分布参数特性的影响。
目前线路在线参数辨识的途径主要有两种,一种是基于稳态数据的参数辨识方法,这类方法通过线路两端的监控与数据采集***(SCADA)和相量采集装置(PMU)的在线实测数据估计线路参数,主要有增广状态估计法、残差灵敏度分析法等,这类方法涉及到的参数和量测都较多,各种误差容易相互影响,并且无法辨识线路的零序参数。
另一种是基于暂态数据的方法,当线路发生外部故障时,通过分析线路两端的故障录波数据信息,分离出序分量,从而达到辨识出该线路的零序与正序参数的目的。由于故障数据中含有非周期分量,利用频域模型计算需要把故障数据转化为向量,在此过程中必然引入误差导致精度不高,故要实现高精度的参数辨识通常需要线路时域模型。但对于线路的分布参数模型来说,目前常用的双曲函数模型为频域模型,只适用于稳态情况,通过时域模型实现对线路分布参数的精确辨识是一个亟待解决的技术问题。
发明内容
本发明的目的是提供一种输电线路分布参数的辨识方法,无需知道线路运行状态与双端量测信息电信号的幅值相位信息,仅在时域内对采集的故障录波数据进行计算,可精确地辨识出线路零序与正序的分布参数,且辨识精度不受故障条件的影响,为电力***计算的准确性奠定基础。
本发明的目的是通过这样的技术方案实现的,它包括有如下的步骤:
(1)、采集发生外部故障输电线路首末两端的电信号瞬时值:
在T(m-1)时间段内,每隔T时间,实时采集输电线路首端的三相电流瞬时信号序列值和三相电压瞬时信号序列值,即:
A相首端电流瞬时信号序列值是:iMA(t1),iMA(t2),…iMA(tm);
B相首端电流瞬时信号序列值是:iMB(t1),iMB(t2),…iMB(tm);
C相首端电流瞬时信号序列值是:iMC(t1),iMC(t2),…iMC(tm);
A相首端电压瞬时信号序列值是:uMA(t1),uMA(t2),…uMA(tm);
B相首端电压瞬时信号序列值是:uMB(t1),uMB(t2),…uMB(tm);
C相首端电压瞬时信号序列值是:uMC(t1),uMC(t2),…uMC(tm);
同时,在T(m-1)时间段内,每隔T时间,实时采集输电线路末端的三相电流瞬时信号序列值和三相电压瞬时信号序列值,即:
A相末端电流瞬时信号序列值是:iNA(t1),iNA(t2),…iNA(tm);
B相末端电流瞬时信号序列值是:iNB(t1),iNB(t2),…iNB(tm);
C相末端电流瞬时信号序列值是:iNC(t1),iNC(t2),…iNC(tm);
A相末端电压瞬时信号序列值是:uNA(t1),uNA(t2),…uNA(tm);
B相末端电压瞬时信号序列值是:uNB(t1),uNB(t2),…uNB(tm);
C相末端电压瞬时信号序列值是:uNC(t1),uNC(t2),…uNC(tm);
其中,5ms≤T≤10.05ms,m是采集电信号的总个数,t1、t2、…tm分别表示各个采集电信号的时刻;
(2)、计算出正序分量和零序分量:
用步骤(1)获得的电信号序列值通过瞬变正弦信号序分量的获取方法分别算出各个时刻的正序分量和零序分量,即:
输电线路首端的电流正序分量是:iM1(t1),iM1(t2),…iM1(tm);
输电线路首端的电流零序分量是:iM0(t1),iM0(t2),…iM0(tm);
输电线路首端的电压正序分量是:uM1(t1),uM1(t2),…uM1(tm);
输电线路首端的电压零序分量是:uM0(t1),uM0(t2),…uM0(tm);
输电线路末端的电流正序分量是:iN1(t1),iN1(t2),…iN1(tm);
输电线路末端的电流零序分量是:iN0(t1),iN0(t2),…iN0(tm);
输电线路末端的电压正序分量是:uN1(t1),uN1(t2),…uN1(tm);
输电线路末端的电压零序分量是:uN0(t1),uN0(t2),…uN0(tm);
(3)、计算发生外部故障输电线路的正序分布参数,即电阻R1,电感L1,电容C1:
①、首先设定线路正序分布参数初始值为R1C,L1C,C1C;
②、用下列公式计算出各个时刻的残差向量eM=[e1、e2、…em]T
式中:
l是输电线路的长度;
j是计算单元的个数,可以取1→∞的正整数;
i是与j值紧密相关的总和循环值,为正整数;
n是正整数,取值为1→m
uM1(tn)是tn时刻线路首端正序电压的瞬时值
uN1(tn)是tn时刻线路末端正序电压的瞬时值
iM1(tn)是tn时刻线路首端正序电流的瞬时值
iN1(tn)是tn时刻线路末端正序电流的瞬时值
iM1‘(tn)是iM1‘(tn)的一阶求导值;
iN1‘(tn)是iN1‘(tn)的一阶求导值;
uN1 (2j-i)(tn)是uN1(tn)的2j-i阶求导值;
iM1 (2j-i+1)(tn)是iM1(tn)的2j-i+1阶求导值;
iN1 (2j-i-1)(tn)是iN1(tn)的2j-i-1阶求导值;
uM1 (2j-i)(t)是uM1(t)的2j-i阶求导值;
③、计算出总残差ε
将ε与所设定的残差阈值ε0相比较:
当ε>ε0时,计算正序分布参数的增量
其中
将R1C+ΔR1c、L1C+ΔL1c、C1C+ΔC1c分别作为新的正序分布参数R1C、L1C、C1C,执行步骤(3)的第②步;
当ε<ε0时,说明此时参入计算的正序分布参数R1C、L1C、C1C就是发生外部故障输电线路的正序分布参数,即电阻R1=R1C,电感L1=L1C,电容C1=C1C;
(4)计算发生外部故障输电线路的零序分布参数,即电阻R0,电感L0,电容C0:
①、首先设定线路零序分布参数初始值为R0C,L0C,C0C;
②、用下列公式计算出各个时刻的残差向量e’M=[e’1、e’2、…e’m]T
式中:
l是输电线路的长度;
j是计算单元的个数,可以取1→∞的正整数;
i是与j值紧密相关的总和循环值,为正整数;
n是正整数,取值为1→m
uM0(tn)是tn时刻线路首端零序电压的瞬时值
uN0(tn)是tn时刻线路末端零序电压的瞬时值
iM0(tn)是tn时刻线路首端零序电流的瞬时值
iN0(tn)是tn时刻线路末端零序电流的瞬时值
iM0‘(tn)是iM0‘(tn)的一阶求导值;
iN0‘(tn)是iN0‘(tn)的一阶求导值;
uN0 (2j-i)(tn)是uN0(tn)的2j-i阶求导值;
iM0 (2j-i+1)(tn)是iM0(tn)的2j-i+1阶求导值;
iN0 (2j-i-1)(tn)是iN0(tn)的2j-i-1阶求导值;
uM0 (2j-i)(t)是uM0(t)的2j-i阶求导值;
③、计算出总残差ε
将ε与所设定的残差阈值ε0相比较:
当ε>ε0时,计算零序分布参数的增量
其中
将R0C+ΔR0c、L0C+ΔL0c、C0C+ΔC0c分别作为新的零序分布参数R0C、L0C、C0C,执行步骤(4)的第②步;
当ε<ε0时,说明此时参入计算的零序分布参数R0C、L0C、C0C就是发生外部故障输电线路的零序分布参数,即电阻R0=R0C,电感L0=L0C,电容C0=C0C。
本发明中,步骤(2)中所述的“瞬变正弦信号序分量的获取方法”属于现有技术,该现有技术源于专利申请号为201010199340.3、名称为“一种电力***瞬变正弦信号序分量的获取方法”的专利申请文件之中,它是先从采集的电流或电压序列值来获取输电线路的零序电压、电流瞬时函数表达式,即:
①、将步骤(1)获取的输电线路首、末端的三相电压瞬时信号uMA(tn)、uMB(tn)、uMC(tn)、uNA(tn)、uNB(tn)、uNC(tn)序列值来得到首、末端的零序电压瞬时信号uM0(tn)和uN0(tn)的序列值,将其分别依次输入到正弦逼近处理器中,经正弦逼近处理器逼近处理后,分别对应地输出随时间变化的第一电压参数系数AM1(tn)、AN1(tn)、AM0(tn)、AN0(tn)和第二电压参数系数BM1(tn)、BN1(tn)、BM0(tn)、BN0(tn),第一电压参数系数AM1(tn)、AN1(tn)、AM0(tn)、AN0(tn)和第二电压参数系数BM1(tn)、BN1(tn)、BM0(tn)、BN0(tn)必定满足瞬时电压信号的瞬变正弦函数表达式:
uM1(tn)=AM1(tn)cosωt+BM1(tn)sinωt;
uN1(tn)=AN1(tn)cosωt+B N1(tn)sinωt;
uM0(tn)=AM0(tn)cosωt+BM0(tn)sinωt;
uN0(tn)=AN0(tn)cosωt+BN0(tn)sinωt;
其中,n为从1到m的正整数,ω为输电线路上电信号的角频率;
②、将步骤(1)获取的输电线路首、末端的三相电流瞬时信号iNA(tn)、iNB(tn)、iNC(tn)、iMA(tn)、iMB(tn)、iMC(tn)序列值来得到首、末端的零序电压瞬时信号iM0(tn)和iN0(tn)的序列值,将其分别依次输入到正弦逼近处理器中,经正弦逼近处理器逼近处理后,分别对应地输出随时间变化的第一电流参数系数AM1(tn)、AN1(tn)、AM0(tn)、AN0(tn)和第二电流参数系数BM1(tn)、BN0(tn)、BM0(tn)、BN0(tn),第一电流参数系数AM1(tn)、AN1(tn)、AM0(tn)、AN0(tn)和第二电流参数系数BM1(tn)、BN0(tn)、BM0(tn)、BN0(tn)必定满足瞬时电流信号的瞬变正弦函数表达式:
iM1(tn)=AM1(tn)cosωt+BM1(tn)sinωt;
iN1(tn)=AN1(tn)cosωt+BN1(tn)sinωt;
iM0(tn)=AM0(tn)cosωt+BM0(tn)sinωt;
iN0(tn)=AN0(tn)cosωt+BN0(tn)sinωt;
其中,n为从1到m的正整数,ω为输电线路上电信号的角频率。
本发明就是将输电线路完全等效为由无穷多个计算单元彼此串联而成的电路模型,如图1所示。每个计算单元是由电阻、电感和电容构成,如图2所示,其中,电阻与电感串联后,一端为单元的输入端,另一端为单元的输出端,且与电容的一端连接,电容的另一端接地。基本思想是将输电线路无穷个计算单元上的电压和电流级联叠加,推导出考虑多阶距离无穷小的输电线路分布参数数学模型,数学模型为线路距离的函数,即由线路一端任意时间的电压和电流可依据数学模型计算出线路线路上任何一节点上的电压,其节点上计算的电压应与真实的电压接近相同。
因此,对于每个计算单元建立的微分方程如下:
un(t)=un-1(t)-RΔxin-1(t)-LΔxin-1‘(t)
in(t)=in-1(t)-CΔxun(t)
上式中:
un(t)表示每个单元输出端的电压;
un-1(t)表示每个单元输入端的电压;
Δx表示每单元的长度;
in-1(t)表示每个单元输入端的电流;
in-1‘(t)表示in-1(t)的一阶求导;
in(t)表示每个单元输出端的电流;
t表示电压或电流进入该单元输入端的时刻;
从输电线路的输入端起,第一个单元的输入端的电信号可以准确采集,然后由上面的两个方程可以解出第一个单元输出端的电压和电流值,并以此作为第二个单元的输入值,同样代入上面的两个方程,又可以解出第二个单元输出端的电压和电流值,以此类推,反复叠加推算,最终得出输电线路分布参数数学模型如下:
其中:
考虑输电线路的多阶距离无穷小的时域分布参数数学模型,与传统的传输线模型相比,基于微分方程的时域表达式不仅更适用于暂态过程,而且模型更加精确。对线路上任意一点来说,可以从线路两端的推算t时刻该点的电压电流值,在从线路端点到该点的一段上,从首端或末端的第一个单元计算开始,每增加一个单元,计算的结果与该点实际值就越接近,而且其接近的量值是随着单元数的增加而显著减少,直到单元个数计算到无穷大时,计算出的该点的电压应当与实际值完全相同。根据对线路中点的两端推算值相当的原理,构建出包含线路参数的函数,可辨识出线路参数,并且辨识精度随单元数的增加而提高。在实际计算中,基于无穷大单元个数的计算是不现实的,只要计算单元的个数能够满足绝对误差和相对误差的精度要求,就可以不用继续增加单元个数的计算,以实现本发明的最终目的,通常取计算单元的个数j=20即可满足参数辨识的精度要求。
由于采用了上述技术方案,本发明只需要采集两端电压、电流同步瞬时值即可求出线路分布参数,无需知道两侧参数及其余部分的运行状态,且该方法辨识精度不受故障条件的影响,具有操作简便和计算准确的优点,可为电网提供准确的线路参数,保证后续电力***计算如状态估计、潮流计算、故障测距、继电保护整定等计算结果的可靠性。
附图说明
图1是输电线路的分布参数等值电路
图2是图1中一个单元的电路图。
具体实施方式
本发明包括有如下的步骤:
(1)、采集发生外部故障输电线路首末两端的电信号瞬时值:
在T(m-1)时间段内,每隔T时间,实时采集输电线路首端的三相电流瞬时信号序列值和三相电压瞬时信号序列值,即:
A相首端电流瞬时信号序列值是:iMA(t1),iMA(t2),…iMA(tm);
B相首端电流瞬时信号序列值是:iMB(t1),iMB(t2),…iMB(tm);
C相首端电流瞬时信号序列值是:iMC(t1),iMC(t2),…iMC(tm);
A相首端电压瞬时信号序列值是:uMA(t1),uMA(t2),…uMA(tm);
B相首端电压瞬时信号序列值是:uMB(t1),uMB(t2),…uMB(tm);
C相首端电压瞬时信号序列值是:uMC(t1),uMC(t2),…uMC(tm);
同时,在T(m-1)时间段内,每隔T时间,实时采集输电线路末端的三相电流瞬时信号序列值和三相电压瞬时信号序列值,即:
A相末端电流瞬时信号序列值是:iNA(t1),iNA(t2),…iNA(tm);
B相末端电流瞬时信号序列值是:iNB(t1),iNB(t2),…iNB(tm);
C相末端电流瞬时信号序列值是:iNC(t1),iNC(t2),…iNC(tm);
A相末端电压瞬时信号序列值是:uNA(t1),uNA(t2),…uNA(tm);
B相末端电压瞬时信号序列值是:uNB(t1),uNB(t2),…uNB(tm);
C相末端电压瞬时信号序列值是:uNC(t1),uNC(t2),…uNC(tm);
其中,5ms≤T≤10.05ms,m是采集电信号的总个数,t1、t2、…tm分别表示各个采集电信号的时刻;
(2)、计算出正序分量和零序分量:
用步骤(1)获得的电信号序列值通过瞬变正弦信号序分量的获取方法分别算出各个时刻的正序分量和零序分量,即:
输电线路首端的电流正序分量是:iM1(t1),iM1(t2),…iM1(tm);
输电线路首端的电流零序分量是:iM0(t1),iM0(t2),…iM0(tm);
输电线路首端的电压正序分量是:uM1(t1),uM1(t2),…uM1(tm);
输电线路首端的电压零序分量是:uM0(t1),uM0(t2),…uM0(tm);
输电线路末端的电流正序分量是:iN1(t1),iN1(t2),…iN1(tm);
输电线路末端的电流零序分量是:iN0(t1),iN0(t2),…iN0(tm);
输电线路末端的电压正序分量是:uN1(t1),uN1(t2),…uN1(tm);
输电线路末端的电压零序分量是:uN0(t1),uN0(t2),…uN0(tm);
(3)、计算发生外部故障输电线路的正序分布参数,即电阻R1,电感L1,电容C1:
①、首先设定线路正序分布参数初始值为R1C,L1C,C1C;
②、用下列公式计算出各个时刻的残差向量eM=[e1、e2、…em]T
式中:
l是输电线路的长度;
j是计算单元的个数,可以取1→∞的正整数;
i是与j值紧密相关的总和循环值,为正整数;
n是正整数,取值为1→m
uM1(tn)是tn时刻线路首端正序电压的瞬时值
uN1(tn)是tn时刻线路末端正序电压的瞬时值
iM1(tn)是tn时刻线路首端正序电流的瞬时值
iN1(tn)是tn时刻线路末端正序电流的瞬时值
iM1‘(tn)是iM1‘(tn)的一阶求导值;
iN1‘(tn)是iN1‘(tn)的一阶求导值;
uN1 (2j-i)(tn)是uN1(tn)的2j-i阶求导值;
iM1 (2j-i+1)(tn)是iM1(tn)的2j-i+1阶求导值;
iN1 (2j-i-1)(tn)是iN1(tn)的2j-i-1阶求导值;
uM1 (2j-i)(t)是uM1(t)的2j-i阶求导值;
③、计算出总残差ε
将ε与所设定的残差阈值ε0相比较:
当ε>ε0时,计算正序分布参数的增量
其中
将R1C+ΔR1c、L1C+ΔL1c、C1C+ΔC1c分别作为新的正序分布参数R1C、L1C、C1C,执行步骤(3)的第②步;
当ε<ε0时,说明此时参入计算的正序分布参数R1C、L1C、C1C就是发生外部故障输电线路的正序分布参数,即电阻R1=R1C,电感L1=L1C,电容C1=C1C;
(4)计算发生外部故障输电线路的零序分布参数,即电阻R0,电感L0,电容C0:
①、首先设定线路零序分布参数初始值为R0C,L0C,C0C;
②、用下列公式计算出各个时刻的残差向量e’M=[e’1、e’2、…e’m]T
式中:
l是输电线路的长度;
j是计算单元的个数,可以取1→∞的正整数;
i是与j值紧密相关的总和循环值,为正整数;
n是正整数,取值为1→m
uM0(tn)是tn时刻线路首端零序电压的瞬时值
uN0(tn)是tn时刻线路末端零序电压的瞬时值
iM0(tn)是tn时刻线路首端零序电流的瞬时值
iN0(tn)是tn时刻线路末端零序电流的瞬时值
iM0‘(tn)是iM0‘(tn)的一阶求导值;
iN0‘(tn)是iN0‘(tn)的一阶求导值;
uN0 (2j-i)(tn)是uN0(tn)的2j-i阶求导值;
iM0 (2j-i+1)(tn)是iM0(tn)的2j-i+1阶求导值;
iN0 (2j-i-1)(tn)是iN0(tn)的2j-i-1阶求导值;
uM0 (2j-i)(t)是uM0(t)的2j-i阶求导值;
③、计算出总残差ε
将ε与所设定的残差阈值ε0相比较:
当ε>ε0时,计算零序分布参数的增量
其中
将R0C+ΔR0c、L0C+ΔL0c、C0C+ΔC0c分别作为新的零序分布参数R0C、L0C、C0C,执行步骤(4)的第②步;
当ε<ε0时,说明此时参入计算的零序分布参数R0C、L0C、C0C就是发生外部故障输电线路的零序分布参数,即电阻R0=R0C,电感L0=L0C,电容C0=C0C。
现结合实验例对本发明作进一步说明:
本实验例所针对的是全长l分别为100km,300km的线路,且电压等级220kv,该输电线路真实参数分别是
表1线路单位长度真实参数
R1[Ω/km] | 0.01273 | L1[H/km] | 9.337e-4 | C1[F/km] | 1.274e-8 |
R0[Ω/km] | 0.3864 | L0[H/km] | 4.126e-3 | C0[F/km] | 7.751e-9 |
实验例1:计算单元个数j对参数辨识精度的影响的检测
由于高压长线路分布电容较大,必须在实际运用中充分考虑线路的分布特性。随着上述线路分布参数数学模型中计算单元个数j的增大,该数学模型分布特性越来越好,也越来越精确。
考虑到线路单位长度的等效电阻电感电容的物理意义设定迭代初值如表2所示:
表2设定线路参数迭代初值
R1c[Ω/km] | 0.01 | L1c[H/km] | 1e-3 | C1c[F/km] | 1e-8 |
R0c[Ω/km] | 0.5 | L0c[H/km] | 5e-3 | C0c[F/km] | 8e-9 |
以单相接地为例,假设长度l为200km与300km的输电线路发生单相接地故障,设定阈值ε0=1*10-4,利用采集到的故障电压电流数据辨识参数,得到结果并求出相对误差,计算单元个数j对辨识精度的影响如表3所示:
表3200km线路发生单相接地,计算单元个数j对辨识精度影响仿真结果比较
表4300km线路发生单相接地,计算单元个数j对辨识精度影响仿真结果比较
表3,表4表明,当j逐渐增大时,该线路分布参数数学模型表示的分布特性越来越好,也越来越接近真实的线路,所以辨识精度越来越高。在实际使用中,计算单元通常取j=20即可满足精度要求。上述的相对误差的定义:相对误差=(计算的线路参数-实际线路参数)/实际线路参数x 100%。
实验例2:在发生不同类型故障,且量测数据受到干扰的情况下,参数辨识精度变化的检测
线路长度l取300km,假设线路分别发生单相接地,两相接地,以及单相断线故障,并在量测数据上叠加标准差为5%的高斯噪声模拟测量噪声,再进行辨识计算,迭代初值如表5所示:
表5设定线路参数迭代初值
R1c[Ω/km] | 0.01 | L1c[H/km] | 1e-3 | C1c[F/km] | 1e-8 |
R0c[Ω/km] | 0.5 | L0c[H/km] | 5e-3 | C0c[F/km] | 8e-9 |
设定阈值ε0=1*10-4,取j=20,计算结果如下表6:
表6线路全长300km,不同类型故障且叠加干扰后仿真结果比较
从表6中可以看出,本发明对各种故障类型均具有较好的适应性,并具有一定的抗干扰能力,辨识精度满足电网参数的精度要求。此外,在辨识结果中,电感的相对误差最小,电容相对误差稍大,电阻的相对误差最大,这是和参数对线路两端的量测量影响是相对应的。
由试验例可知,本发明提出的在时域内利用暂态数据对长线路分布参数进行辨识的方法可精确辨识出线路的正序与零序参数,对不同类型的故障具有良好的适应性,并且在数据受干扰情况下的依然能保持较高的辨识精度。该方法可以有效避免频域法在向量转换过程中引入误差的问题,从而有效提高了参数辨识结果的准确度,并且无需电压电流的幅值相位信息,降低了对量测数据的要求。
Claims (1)
1.一种输电线路分布参数的辨识方法,它包括有如下的步骤:
(1)、采集发生外部故障输电线路首末两端的电信号瞬时值:
在T(m-1)时间段内,每隔T时间,实时采集输电线路首端的三相电流瞬时信号序列值和三相电压瞬时信号序列值,即:
A相首端电流瞬时信号序列值是:iMA(t1),iMA(t2),…iMA(tm);
B相首端电流瞬时信号序列值是:iMB(t1),iMB(t2),…iMB(tm);
C相首端电流瞬时信号序列值是:iMC(t1),iMC(t2),…iMC(tm);
A相首端电压瞬时信号序列值是:uMA(t1),uMA(t2),…uMA(tm);
B相首端电压瞬时信号序列值是:uMB(t1),uMB(t2),…uMB(tm);
C相首端电压瞬时信号序列值是:uMC(t1),uMC(t2),…uMC(tm);
同时,在T(m-1)时间段内,每隔T时间,实时采集输电线路末端的三相电流瞬时信号序列值和三相电压瞬时信号序列值,即:
A相末端电流瞬时信号序列值是:iNA(t1),iNA(t2),…iNA(tm);
B相末端电流瞬时信号序列值是:iNB(t1),iNB(t2),…iNB(tm);
C相末端电流瞬时信号序列值是:iNC(t1),iNC(t2),…iNC(tm);
A相末端电压瞬时信号序列值是:uNA(t1),uNA(t2),…uNA(tm);
B相末端电压瞬时信号序列值是:uNB(t1),uNB(t2),…uNB(tm);
C相末端电压瞬时信号序列值是:uNC(t1),uNC(t2),…uNC(tm);
其中,5ms≤T≤10.05ms,m是采集电信号的总个数,t1、t2、…tm分别表示各个采集电信号的时刻;
(2)、计算出正序分量和零序分量:
用步骤(1)获得的电信号序列值通过瞬变正弦信号序分量的获取方法分别算出各个时刻的正序分量和零序分量,即:
输电线路首端的电流正序分量是:iM1(t1),iM1(t2),…iM1(tm);
输电线路首端的电流零序分量是:iM0(t1),iM0(t2),…iM0(tm);
输电线路首端的电压正序分量是:uM1(t1),uM1(t2),…uM1(tm);
输电线路首端的电压零序分量是:uM0(t1),uM0(t2),…uM0(tm);
输电线路末端的电流正序分量是:iN1(t1),iN1(t2),…iN1(tm);
输电线路末端的电流零序分量是:iN0(t1),iN0(t2),…iN0(tm);
输电线路末端的电压正序分量是:uN1(t1),uN1(t2),…uN1(tm);
输电线路末端的电压零序分量是:uN0(t1),uN0(t2),…uN0(tm);
(3)、计算发生外部故障输电线路的正序分布参数,即电阻R1,电感L1,电容C1:
①、首先设定线路正序分布参数初始值为R1C,L1C,C1C;
②、用下列公式计算出各个时刻的残差向量eM=[e1、e2、…em]T
式中:
l是输电线路的长度;
j是计算单元的个数,可以取1→∞的正整数;
i是与j值紧密相关的总和循环值,为正整数;
n是正整数,取值为1→m
uM1(tn)是tn时刻线路首端正序电压的瞬时值
uN1(tn)是tn时刻线路末端正序电压的瞬时值
iM1(tn)是tn时刻线路首端正序电流的瞬时值
iN1(tn)是tn时刻线路末端正序电流的瞬时值
iM1‘(tn)是iM1‘(tn)的一阶求导值;
iN1‘(tn)是iN1‘(tn)的一阶求导值;
uN1 (2j-i)(tn)是uN1(tn)的2j-i阶求导值;
iM1 (2j-i+1)(tn)是iM1(tn)的2j-i+1阶求导值;
iN1 (2j-i-1)(tn)是iN1(tn)的2j-i-1阶求导值;
uM1 (2j-i)(t)是uM1(t)的2j-i阶求导值;
③、计算出总残差ε
将ε与所设定的残差阈值ε0相比较:
当ε>ε0时,计算正序分布参数的增量
其中
将R1C+ΔR1c、L1C+ΔL1c、C1C+ΔC1c分别作为新的正序分布参数R1C、L1C、C1C,执行步骤(3)的第②步;
当ε<ε0时,说明此时参入计算的正序分布参数R1C、L1C、C1C就是发生外部故障输电线路的正序分布参数,即电阻R1=R1C,电感L1=L1C,电容C1=C1C;
(4)计算发生外部故障输电线路的零序分布参数,即电阻R0,电感L0,电容C0:
①、首先设定线路零序分布参数初始值为R0C,L0C,C0C;
②、用下列公式计算出各个时刻的残差向量e’M=[e’1、e’2、…e’m]T
式中:
l是输电线路的长度;
j是计算单元的个数,可以取1→∞的正整数;
i是与j值紧密相关的总和循环值,为正整数;
n是正整数,取值为1→m
uM0(tn)是tn时刻线路首端零序电压的瞬时值
uN0(tn)是tn时刻线路末端零序电压的瞬时值
iM0(tn)是tn时刻线路首端零序电流的瞬时值
iN0(tn)是tn时刻线路末端零序电流的瞬时值
iM0‘(tn)是iM0‘(tn)的一阶求导值;
iN0‘(tn)是iN0‘(tn)的一阶求导值;
uN0 (2j-i)(tn)是uN0(tn)的2j-i阶求导值;
iM0 (2j-i+1)(tn)是iM0(tn)的2j-i+1阶求导值;
iN0 (2j-i-1)(tn)是iN0(tn)的2j-i-1阶求导值;
uM0 (2j-i)(t)是uM0(t)的2j-i阶求导值;
③、计算出总残差ε
将ε与所设定的残差阈值ε0相比较:
当ε>ε0时,计算零序分布参数的增量
其中
将R0C+ΔR0c、L0C+ΔL0c、C0C+ΔC0c分别作为新的零序分布参数R0C、L0C、C0C,执行步骤(4)的第②步;
当ε<ε0时,说明此时参入计算的零序分布参数R0C、L0C、C0C就是发生外部故障输电线路的零序分布参数,即电阻R0=R0C,电感L0=L0C,电容C0=C0C。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210458824.4A CN102967779B (zh) | 2012-11-14 | 2012-11-14 | 一种输电线路分布参数的辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210458824.4A CN102967779B (zh) | 2012-11-14 | 2012-11-14 | 一种输电线路分布参数的辨识方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102967779A true CN102967779A (zh) | 2013-03-13 |
CN102967779B CN102967779B (zh) | 2014-11-12 |
Family
ID=47798059
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210458824.4A Active CN102967779B (zh) | 2012-11-14 | 2012-11-14 | 一种输电线路分布参数的辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102967779B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103837800A (zh) * | 2014-03-18 | 2014-06-04 | 国家电网公司 | 直流输电线路单极接地时故障位置定位方法 |
CN107123348A (zh) * | 2017-06-01 | 2017-09-01 | 国网湖北省电力公司电力科学研究院 | 一种高精度中压配电线路真型等效模型及其参数计算方法 |
CN108226710A (zh) * | 2017-12-29 | 2018-06-29 | 中国电力科学研究院有限公司 | 一种输电线路参数辨识方法及*** |
CN109001589A (zh) * | 2018-05-03 | 2018-12-14 | 西安理工大学 | 基于非全相运行故障录波数据的线路参数时域计算方法 |
CN109375048A (zh) * | 2018-07-26 | 2019-02-22 | 广东电网有限责任公司 | 一种基于故障录波数据的输电线路参数辨识方法及装置 |
CN111242459A (zh) * | 2020-01-07 | 2020-06-05 | 中国南方电网有限责任公司 | 全网设备参数异常值辨识方法及*** |
CN112083237A (zh) * | 2020-07-31 | 2020-12-15 | 西安交通大学 | 一种用于大尺度电气设备宽频特性时域测量方法及*** |
CN113702762A (zh) * | 2021-08-27 | 2021-11-26 | 山东科汇电力自动化股份有限公司 | 一种利用零序信息量的配电网单相接地故障测距方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6072317A (en) * | 1998-03-27 | 2000-06-06 | Eaton Corporation | Plug-in multifunction tester for AC electrical distribution system |
JP2005033982A (ja) * | 2003-07-11 | 2005-02-03 | Sungkyunkwan Univ | 伝送線路で高調波と状態遷移図を用いた故障検出方法 |
CN101881792A (zh) * | 2010-06-12 | 2010-11-10 | 重庆电力科学试验研究院 | 一种电力***瞬变正弦信号的序分量获取方法 |
CN102082420A (zh) * | 2010-12-08 | 2011-06-01 | 重庆大学 | 输电线路的纵差保护方法 |
-
2012
- 2012-11-14 CN CN201210458824.4A patent/CN102967779B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6072317A (en) * | 1998-03-27 | 2000-06-06 | Eaton Corporation | Plug-in multifunction tester for AC electrical distribution system |
JP2005033982A (ja) * | 2003-07-11 | 2005-02-03 | Sungkyunkwan Univ | 伝送線路で高調波と状態遷移図を用いた故障検出方法 |
CN101881792A (zh) * | 2010-06-12 | 2010-11-10 | 重庆电力科学试验研究院 | 一种电力***瞬变正弦信号的序分量获取方法 |
CN102082420A (zh) * | 2010-12-08 | 2011-06-01 | 重庆大学 | 输电线路的纵差保护方法 |
Non-Patent Citations (2)
Title |
---|
GERARDUS C. PAAP: "Symmetrical Components in the Time Domain and Their Application to Power Network Calculations", 《TRANSACTIONS ON POWER SYSTEMS》 * |
王增平等: "基于分布参数模型的T型输电线路电流差动保护新原理", 《电网技术》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103837800A (zh) * | 2014-03-18 | 2014-06-04 | 国家电网公司 | 直流输电线路单极接地时故障位置定位方法 |
CN107123348A (zh) * | 2017-06-01 | 2017-09-01 | 国网湖北省电力公司电力科学研究院 | 一种高精度中压配电线路真型等效模型及其参数计算方法 |
CN108226710A (zh) * | 2017-12-29 | 2018-06-29 | 中国电力科学研究院有限公司 | 一种输电线路参数辨识方法及*** |
CN109001589A (zh) * | 2018-05-03 | 2018-12-14 | 西安理工大学 | 基于非全相运行故障录波数据的线路参数时域计算方法 |
CN109001589B (zh) * | 2018-05-03 | 2021-01-15 | 西安理工大学 | 基于非全相运行故障录波数据的线路参数时域计算方法 |
CN109375048A (zh) * | 2018-07-26 | 2019-02-22 | 广东电网有限责任公司 | 一种基于故障录波数据的输电线路参数辨识方法及装置 |
CN109375048B (zh) * | 2018-07-26 | 2021-04-23 | 广东电网有限责任公司 | 一种基于故障录波数据的输电线路参数辨识方法及装置 |
CN111242459A (zh) * | 2020-01-07 | 2020-06-05 | 中国南方电网有限责任公司 | 全网设备参数异常值辨识方法及*** |
CN112083237A (zh) * | 2020-07-31 | 2020-12-15 | 西安交通大学 | 一种用于大尺度电气设备宽频特性时域测量方法及*** |
CN112083237B (zh) * | 2020-07-31 | 2021-06-29 | 西安交通大学 | 一种用于大尺度电气设备宽频特性时域测量方法及*** |
CN113702762A (zh) * | 2021-08-27 | 2021-11-26 | 山东科汇电力自动化股份有限公司 | 一种利用零序信息量的配电网单相接地故障测距方法 |
CN113702762B (zh) * | 2021-08-27 | 2024-04-16 | 山东科汇电力自动化股份有限公司 | 一种利用零序信息量的配电网单相接地故障测距方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102967779B (zh) | 2014-11-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102967779B (zh) | 一种输电线路分布参数的辨识方法 | |
CN102331548B (zh) | 一种配网电缆单相接地故障的故障定位方法 | |
CN102967842B (zh) | 一种电子式电流互感器的渐变性故障在线诊断方法 | |
CN103941147B (zh) | 利用暂态主频分量的配网电缆单相接地故障测距方法 | |
CN107632236B (zh) | 一种基于对端母线反射波识别的单出线输电线路单端故障测距方法 | |
CN105738769B (zh) | 基于分布参数模型的串联补偿双回线故障定位方法 | |
CN103792465A (zh) | 一种基于零序电压的配电网单相接地故障测距的方法 | |
CN102570428B (zh) | 基于电子式互感器微分输出的故障测距与距离保护方法 | |
CN104898021A (zh) | 一种基于k-means聚类分析的配电网故障选线方法 | |
CN105067948A (zh) | 一种小电流接地选线装置及单相接地检测方法 | |
CN111141995B (zh) | 一种基于比幅原理的线路双端稳态量测距方法和*** | |
CN102621388A (zh) | 基于同步时域信号的输电线路集中参数在线确定方法 | |
CN103207308A (zh) | 避雷器阻性电流和容性电流暂态值的测量方法 | |
CN107621591A (zh) | 一种基于零模行波波速变化特性的输电线路迭代测距方法 | |
CN103427405B (zh) | 基于高阶累积量的输电线路差动保护方法 | |
CN116338525B (zh) | 一种风电交流送出线路故障测距方法及*** | |
CN107329045A (zh) | 配网故障最小二乘定位算法 | |
CN105138843A (zh) | 电力***采样飞点检测及其修复方法 | |
CN104297628A (zh) | 含dg的配电网的区段故障检测与定位方法 | |
CN113253052A (zh) | 一种基于改进smmg的高压直流输电线路故障测距方法 | |
CN107247215B (zh) | 基于多点同步测量数据的配网故障粒子群定位算法 | |
CN104198889B (zh) | 基于相继动作的高压线路瞬时性接地故障单端测距方法 | |
CN110161375B (zh) | 一种基于分布电阻参数的高压直流输电线路计算模型 | |
CN115015687A (zh) | 一种四端环状柔性直流电网故障测距方法及*** | |
CN102082420B (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |