CN106226597A - 次同步谐振监测方法 - Google Patents

次同步谐振监测方法 Download PDF

Info

Publication number
CN106226597A
CN106226597A CN201610673636.1A CN201610673636A CN106226597A CN 106226597 A CN106226597 A CN 106226597A CN 201610673636 A CN201610673636 A CN 201610673636A CN 106226597 A CN106226597 A CN 106226597A
Authority
CN
China
Prior art keywords
subsynchronous resonance
low
value
snap shot
sampled value
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
Application number
CN201610673636.1A
Other languages
English (en)
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.)
State Grid Corp of China SGCC
Xuji Group Co Ltd
XJ Electric Co Ltd
Xuchang XJ Software Technology Co Ltd
Original Assignee
State Grid Corp of China SGCC
Xuji Group Co Ltd
XJ Electric Co Ltd
Xuchang XJ Software Technology Co Ltd
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 State Grid Corp of China SGCC, Xuji Group Co Ltd, XJ Electric Co Ltd, Xuchang XJ Software Technology Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201610673636.1A priority Critical patent/CN106226597A/zh
Publication of CN106226597A publication Critical patent/CN106226597A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/02Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Measurement Of Resistance Or Impedance (AREA)

Abstract

本发明涉及次同步谐振监测方法,该方法包括如下步骤:1)实时采集串补系统中线路电流,并对线路电流进行低通滤波处理;2)对低通滤波处理后的线路电流进行抽点,所述抽点是指每隔M个采样值抽取1个采样值,从而将高采样率采样值转换为低采样率采样值,其中M>1;3)对抽点后的采样值进行带通滤波,生成次同步谐振电流;4)计算次同步谐振电流的有效值ISSR,并判断ISSR是否大于次同步谐振电流定值,若大于,则旁路串补系统。本发明的方法消除了次同步谐振产生的回路参数条件,保护发电机轴系安全,提高系统稳定性和可靠性,通过采用抽点、长窗有效值算法可分段计算分段更新然后累加等措施,计算量较小,方便继电保护实现。

Description

次同步谐振监测方法
技术领域
本发明属于电力系统继电保护技术领域,具体涉及一种次同步谐振监测方法。
背景技术
继电保护装置要求在系统发生故障时能可靠地切除故障设备,保证系统的安全稳定运行,但保护装置的不正确动作(包括拒动和误动)造成的危害也是巨大的。
近年来,多条特高压线路投入运行,为提高线路输送容量,特高压线路均采用了串联电容器方案。另外,对于500kV以上电压等级的长线路,为避免新建线路所需的巨大投资,也多采用增加串补电容来扩容的方法。
线路增加串补电容后,虽然扩大了线路功率输送容量,但同时也增加了系统出现次同步谐振的风险。如在系统出现扰动时,会导致同步电机功角摇摆,若摇摆频率与电气系统谐振同频并与电机轴系固有振荡频率发生互补(如电气谐振频率与轴系固有频率之和为50Hz),则可能发生次同步谐振现象。频率低于50Hz的谐振通常认为是SSR(次同步谐振)。出现SSR时,电容器组应被旁路,否则会严重损坏谐振的发电机轴系。
专利文献《串补系统中次同步谐振的检测方法》(申请号:200910029127.5)给出了一种次同步谐振SSR的检测方法,即对于采集到的线路电流采用低通、带阻级联滤波器进行滤波,得到次同步谐振电流,然后求次同步谐振电流的有效值,当次同步谐振电流的有效值大于次同步谐振电流定值时,旁路串补系统。其中对线路电流处理获取的次同步电流的过程计算复杂,计算量大,所以,需要一种简单易实现的方法实现次同步谐振的监测。
发明内容
本发明提供了次同步谐振监测方法,解决现有的次同步谐振检测方法处理过程复杂、计算量大的缺陷。
为解决上述技术问题,本发明的次同步谐振监测方法包括:
1)实时采集串补系统中线路电流,并对线路电流进行低通滤波处理;
2)对低通滤波处理后的线路电流进行抽点,所述抽点是指每隔M个采样值抽取1个采样值,从而将高采样率采样值转换为低采样率采样值,其中M>1;
3)对抽点后的线路电流采样值进行带通滤波,得到次同步谐振电流;
4)计算次同步谐振电流的有效值ISSR,并判断ISSR是否大于次同步谐振电流定值,若大于,则旁路串补系统。
所述步骤1)中对线路电流进行低通滤波处理时的低通滤波器公式为:
y ( n ) = Σ k = 0 11 b ( k ) x ( n - k )
其中,y(n)为低通滤波后采样值,x(n)为原始采样值,b(k)为低通滤波器系数。
所述步骤2)中对低通滤波处理后的线路电流进行抽点的公式为:
y(n)=x(4k)
其中y(n)为抽点后采样值;x(4k)为经低通滤波后的线路电流采样值;n、k为自然数序列。
所述步骤3中对抽点后的采样值进行带通滤波所采用的带通滤波器由6个二阶滤波器级联而成,各二阶滤波器的公式为:
x 1 ( n ) = b 01 x ( n ) + b 11 x ( n - 1 ) + b 21 x ( n - 2 ) - [ a 11 x 1 ( n - 1 ) + a 21 x 1 ( n - 2 ) ] x k ( n ) = b 0 k x k - 1 ( n ) + b 1 k x k - 1 ( n - 1 ) + b 2 k x k - 1 ( n - 2 ) - [ a 1 k x k ( n - 1 ) + a 2 k x k ( n - 2 ) ] y ( n ) = A x 6 ( n )
其中,k=2,3,...,6,x(n)为带通滤波前的采样值,y(n)为带通滤波后采样值,分别为带通滤波系数,A为增益系数。
所述步骤4)计算次同步谐振电流有效值ISSR采用长窗有效值算法,公式为:
I S S R = Σ k = 0 k = N - 1 i k 2 / N
其中,ik为次同步谐振电流,N为时间窗大小。
N=300。
所述步骤3)中带通滤波的通带频率为7~42Hz。
本发明的有益效果:本发明的方法通过对串补线路电流采样值进行低通滤波、抽点运算、带通滤波,仅保留有效的次同步谐振分量电流,并通过长窗有效值算法计算出次同步谐振分量电流,低通滤波与抽点算法用于降低后续算法计算量,带通滤波算法用于提取出线路电流中的次同步谐振分量,长窗有效值算法用于计算线路电流低频分量有效值大小,当次同步谐振电流分量大于定值门槛时,经延时确认后旁路线路串补系统,从而消除谐振产生的回路参数条件,保护发电机轴系安全,提高系统稳定性和可靠性。本发明通过采用抽点、长窗有效值算法可分段计算分段更新然后累加等措施,计算量较小,方便继电保护实现。
附图说明
图1为本实施例的串补系统结构图;
图2为本实施例的次同步谐振监测方法流程图;
图3为本实施例的带通滤波器结构图。
具体实施方式
下面结合附图,对本发明的技术方案作进一步详细的说明。
次同步谐振监测方法主要应用于超高压及特高压串补系统继电保护装置中,作为一个独立的功能模块,集成在串补系统继电保护装置中,该功能模块始终投入,其流程如图2所示。具体实现过程为:
1)实时采集串补系统中线路电流,并对线路电流进行低通滤波处理;
2)对低通滤波处理后的线路电流进行抽点,所述抽点是指每隔M个采样值抽取1个采样值,从而将高采样率采样值转换为低采样率采样值,其中M>1;
3)对抽点后的线路电流采样值进行带通滤波,得到次同步谐振电流;
4)计算次同步谐振电流的有效值ISSR,并判断ISSR是否大于次同步谐振电流定值,若大于,则旁路串补系统。
下面对上述步骤进行详细阐述:
1、首先实时采集串补系统三相线路电流,并对各相电流进行低通滤波,滤波公式如下:
y ( n ) = Σ k = 0 11 b ( k ) x ( n - k )
其中,y(n)为滤波后采样值,x(n)为原始采样值,b(k)为滤波器系数,b(k)可优选为b=[-0.014626139573560 -0.025006738161872 -0.0028330908309590.075858546295221 0.190287862792706 0.276488181623457 0.2764881816234570.190287862792706 0.075858546295221 -0.002833090830959 -0.025006738161872 -0.014626139573560]。
2、对经低通滤波后的线路电流采样值进行抽点,将高采样率转换为低采样率,降低计算量。抽点算法公式为:
y(n)=x(4k)
其中y(n)为抽点后采样值;x(4k)为经低通滤波后的线路电流采样值;n、k为1、2、3……自然数序列;本实施例抽点前的采样率为每秒4000次、抽点后采样率为每秒1000次。
3.对抽点后的低采样率线路电流采样值进行带通滤波计算,以滤除次同步谐振频率之外的成分,仅保留次同步谐振频率成分。带通滤波器由6个二阶滤波器级联而成,其结构如图3所示,各滤波器具体公式为:
x 1 ( n ) = b 01 x ( n ) + b 11 x ( n - 1 ) + b 21 x ( n - 2 ) - [ a 11 x 1 ( n - 1 ) + a 21 x 1 ( n - 2 ) ] x k ( n ) = b 0 k x k - 1 ( n ) + b 1 k x k - 1 ( n - 1 ) + b 2 k x k - 1 ( n - 2 ) - [ a 1 k x k ( n - 1 ) + a 2 k x k ( n - 2 ) ] y ( n ) = A x 6 ( n )
其中,k=2,3,...,6,x(n)为带通滤波前的采样值,y(n)为带通滤波后采样值,分别为带通滤波系数,A为增益系数。可优选为[1.0,-1.0712198125161931,0.25375374337006074,1.0,-1.8896824277577728,0.95055847065015886;
1.0,0.042103342323017712,0.0011041087284387788,1.0,-1.9727728412407597,0.97442666890251062;
1.0,-1.0865673927717474,0.50822439050164991,1.0,-1.8752648911479859,0.91819128979385733;
1.0,-1.1677101765793967,0.16771043714493428,1.0,-1.8761705220985851,0.8991747756194588;
1.0,-1.8757065301390308,0.98309920183577482,1.0,-1.903765665799370,0.91154560677330532;
1.0,-1.2072499843350848,0.20728012362303741,1.0,-1.9129919796696275,0.98410952762209791];
增益系数A优选为0.0000026264783383456857。
本实施例中带通滤波的通带频率为7~42Hz。
4.计算带通滤波后的次同步谐振电流有效值,具体算法如下:
I S S R = Σ k = 0 k = N - 1 i k 2 / N
其中,ik为次同步谐振电流,N为时间窗大小。
为降低有效低频分量的计算误差,可采用300ms时间窗计算ISSR有效值,300ms时间窗对应的N为300点。
5.判别次同步谐振电流有效值ISSR是否大于次同步谐振定值,若大于次同步谐振定值,则旁路串补系统,否则,不做任何操作。
本实施例的次同步谐振监测方法通过对线路电流进行低通滤波、抽点、带通滤波、长窗有效值计算,得出次同步谐振电流分量,并判别该分量大小,若该分量大于次同步谐振旁路定值,则旁路串补系统,从而改变系统参数,消除谐振条件,保护发电机轴系的安全。次同步谐振保护电流取线路电流,如图1所示CT2。该方法的特点为通过采用抽点、长窗有效值算法可分段计算分段更新然后累加等措施,计算量较小,方便继电保护实现。
以上给出了具体的实施方式,但本发明不局限于所描述的实施方式。本发明的基本思路在于上述基本方案,对本领域普通技术人员而言,根据本发明的教导,设计出各种变形的模型、公式、参数并不需要花费创造性劳动。在不脱离本发明的原理和精神的情况下对实施方式进行的变化、修改、替换和变型仍落入本发明的保护范围内。

Claims (7)

1.次同步谐振监测方法,其特征在于,该方法包括如下步骤:
1)实时采集串补系统中线路电流,并对线路电流进行低通滤波处理;
2)对低通滤波处理后的线路电流进行抽点,所述抽点是指每隔M个采样值抽取1个采样值,从而将高采样率采样值转换为低采样率采样值,其中M>1;
3)对抽点后的线路电流采样值进行带通滤波,得到次同步谐振电流;
4)计算次同步谐振电流的有效值ISSR,并判断ISSR是否大于次同步谐振电流定值,若大于,则旁路串补系统。
2.根据权利要求1所述次同步谐振监测方法,其特征在于,所述步骤1)中对线路电流进行低通滤波处理时的低通滤波器公式为:
y ( n ) = Σ k = 0 11 b ( k ) x ( n - k )
其中,y(n)为低通滤波后采样值,x(n)为原始采样值,b(k)为低通滤波器系数。
3.根据权利要求1所述次同步谐振监测方法,其特征在于,所述步骤2)中对低通滤波处理后的线路电流进行抽点的公式为:
y(n)=x(4k)
其中y(n)为抽点后采样值;x(4k)为经低通滤波后的线路电流采样值;n、k为自然数序列。
4.根据权利要求1所述次同步谐振监测方法,其特征在于,所述步骤3中对抽点后的采样值进行带通滤波所采用的带通滤波器由6个二阶滤波器级联而成,各二阶滤波器的公式为:
x 1 ( n ) = b 01 x ( n ) + b 11 x ( n - 1 ) + b 21 x ( n - 2 ) - [ a 11 x 1 ( n - 1 ) + a 21 x 1 ( n - 2 ) ] x k ( n ) = b 0 k x k - 1 ( n ) + b 1 k x k - 1 ( n - 1 ) + b 2 k x k - 1 ( n - 2 ) - [ a 1 k x k ( n - 1 ) + a 2 k x k ( n - 2 ) ] y ( n ) = Ax 6 ( n )
其中,k=2,3,...,6,x(n)为带通滤波前的采样值,y(n)为带通滤波后采样值,分别为带通滤波系数,A为增益系数。
5.根据权利要求1所述次同步谐振监测方法,其特征在于,所述步骤4)计算次同步谐振电流有效值ISSR采用长窗有效值算法,公式为:
I S S R = Σ k = 0 k = N - 1 i k 2 / N
其中,ik为次同步谐振电流,N为时间窗大小。
6.根据权利要求5所述次同步谐振监测方法,其特征在于,N=300。
7.根据权利要求1所述次同步谐振监测方法,其特征在于,所述步骤3)中带通滤波的通带频率为7~42Hz。
CN201610673636.1A 2016-08-16 2016-08-16 次同步谐振监测方法 Pending CN106226597A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610673636.1A CN106226597A (zh) 2016-08-16 2016-08-16 次同步谐振监测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610673636.1A CN106226597A (zh) 2016-08-16 2016-08-16 次同步谐振监测方法

Publications (1)

Publication Number Publication Date
CN106226597A true CN106226597A (zh) 2016-12-14

Family

ID=57552431

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610673636.1A Pending CN106226597A (zh) 2016-08-16 2016-08-16 次同步谐振监测方法

Country Status (1)

Country Link
CN (1) CN106226597A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109038619A (zh) * 2018-08-16 2018-12-18 东北大学 一种基于旁路电流补偿的次同步振荡抑制装置及方法
CN111473023A (zh) * 2020-04-22 2020-07-31 中国飞机强度研究所 一种液压缸谐振智能监测系统及定位方法
CN113672863A (zh) * 2021-07-08 2021-11-19 南京国电南自电网自动化有限公司 一种发电机启机保护的简化相量计算方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3944846A (en) * 1974-03-19 1976-03-16 Southern California Edison Company Subsynchronous relay
CN1093171A (zh) * 1994-01-31 1994-10-05 清华大学 监测电力系统次同步振荡的方法及监测仪
CN101261293A (zh) * 2007-03-08 2008-09-10 国网南京自动化研究院 基于自适应滤波器的电力稳态信号跟踪测量法
CN101465545A (zh) * 2009-01-08 2009-06-24 华北电网有限公司张家口供电公司 串补系统中次同步谐振的检测方法
CN102753982A (zh) * 2010-02-17 2012-10-24 Abb研究有限公司 次同步谐振检测的方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3944846A (en) * 1974-03-19 1976-03-16 Southern California Edison Company Subsynchronous relay
CN1093171A (zh) * 1994-01-31 1994-10-05 清华大学 监测电力系统次同步振荡的方法及监测仪
CN101261293A (zh) * 2007-03-08 2008-09-10 国网南京自动化研究院 基于自适应滤波器的电力稳态信号跟踪测量法
CN101465545A (zh) * 2009-01-08 2009-06-24 华北电网有限公司张家口供电公司 串补系统中次同步谐振的检测方法
CN102753982A (zh) * 2010-02-17 2012-10-24 Abb研究有限公司 次同步谐振检测的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
孙妙平 等: "串补电路中次同步谐振检测和抑制的研究", 《铁道学报》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109038619A (zh) * 2018-08-16 2018-12-18 东北大学 一种基于旁路电流补偿的次同步振荡抑制装置及方法
CN111473023A (zh) * 2020-04-22 2020-07-31 中国飞机强度研究所 一种液压缸谐振智能监测系统及定位方法
CN113672863A (zh) * 2021-07-08 2021-11-19 南京国电南自电网自动化有限公司 一种发电机启机保护的简化相量计算方法及系统

Similar Documents

Publication Publication Date Title
WO2015165286A1 (zh) Lc并联回路失谐故障的继电保护方法和装置
CN102955060B (zh) 一种去除电力系统故障信号中衰减直流分量的方法
WO2015081826A1 (zh) 电力线路保护方法
CN104638664B (zh) 风电厂次同步振荡保护方法及保护装置
CN106226597A (zh) 次同步谐振监测方法
CN110954743B (zh) 一种分布式录波装置及小电流接地选线方法
CN108051702B (zh) 基于单相接地故障录波数据的故障线路参数计算方法
CN109188174B (zh) 一种交流线路出口故障的判别方法
CN104466920B (zh) 一种断路器失灵保护方法
CN102255292A (zh) 一种基于参数识别的高压输电线路距离保护方法
CN101051738A (zh) 特高压输电线路的负序方向纵联保护方法
CN102082420B (zh) 输电线路的纵差保护方法
CN110703134B (zh) 一种基于故障序分量的小电流接地选线选相方法
CN107656134A (zh) 一种自适应滤除直流分量的故障电流检测方法及装置
Saha et al. DWT and BPNN based fault detection, classification and estimation of location of HVAC transmission line
CN112039036B (zh) 一种交流线路瞬时值变化量距离保护方法及系统
CN107015115B (zh) 一种同塔双回输电线路的故障测距方法
CN105223403A (zh) 一种并网逆变器网侧谐振电流信息的小波包提取方法
CN110601153B (zh) 直流配电网单端量的继电保护方法
CN109001589B (zh) 基于非全相运行故障录波数据的线路参数时域计算方法
Zou et al. Mathematical morphology based phase selection scheme in digital relaying
CN107765074A (zh) 牵引变流器谐波电流信号强度的检测方法及检测系统
Bhupatiraju et al. A wavelet based protection scheme for EHV asymmetrical teed circuits
Paul et al. ANFIS based single line to ground fault location estimation for transmission lines
Sidhu et al. A fast distance relay using adaptive data window filters

Legal Events

Date Code Title Description
C06 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: 20161214

RJ01 Rejection of invention patent application after publication