专利名称：利用极波电压s变换能量比值的特高压直流输电线路故障边界元件的方法
技术领域：
本发明涉及电力系统继电保护技术领域，具体地说是一种利用极波电压S变换能 量比值的特高压直流输电线路故障边界元件的方法。
背景技术：
现有的直流线路保护主要由ABB和SIEMENS两家公司提供的，其快速的单端量保 护是依据电压行波的突变量、电压行波变化率和电流变化梯度作为判据的。虽然行波保护 具有超高速动作特性，且不受电流互感器饱和、系统振荡和长线分布电容等的影响。但根据 国内外实际运行的相关资料显示，目前使用的直流线路的行波保护受以下的因素的影响(1)电压行波受线路末端平波电抗器和直流滤波器的影响，而导致电压变化率减 小。(2)方向行波变化率主要受过渡电阻的影响。3 电压变量和地模波受故障距离的影响。因此以变化率判据为基础直流线路保护易受上述因素的影响，而产生误动。直流 输电系统包括直流输电线路和直流线路两端的平波电抗器和直流滤波器，它们构成了直流 输电线路高频暂态量的“天然”边界，并且边界频率特性分析表明线路边界内、外故障信号 的高频分量存在显著的特征差异，据此可以构成单端暂态量保护元件。本发明是基于利用高频分量特征的直流输电线路单端超高速保护原理(其原理 框图如图1)上提出的一种边界元件的方法。
发明内容
本发明的目的是提供一种利用极波电压S变换能量比值的特高压直流输电线路 故障边界元件的方法。本发明的利用极波电压S变换能量比值的特高压直流输电线路故障边界元件的 方法按以下步骤进行(1)直流线路发生故障后，启动元件立即启动，根据保护安装处测得的两极直流电 压、直流电流，求出正极线极波P1GO,负极线极波P2(k)为P1GO = ZpXi1GO-U1GO⑴P2 (k) = Zp X i2 (k)-U2 (k)(2)式中，Zp为直流输电线路极波阻抗，U1 (k)为正极直流电压，ii(k)为正极直流电 流，u2(k)为负极直流电压，i2(k)为负极直流电流，k= 1、2、3....N，N为采样序列长度；(2)对极波电压进行S变换，其采样频率为100kHz，采样序列长度为200，经S变换 得到101X200的复矩阵 式中x[k]为采集到的N个离散的极波电压信号点，k = 0、l、2....N-l，N为采样 序列长度。X[n]为x[k]的傅里叶变换。S[m，η]为复时频矩阵，其列向量表示信号在某一 时刻的幅频特性，其行向量表示信号在某一频率的时域分布；(3)对得到的复矩阵中的各个元素求模，得到模矩阵，根据下式计算高频能量和低
频能量 式中|S[m，η] |为将S矩阵中的各个元素求模后得到的模矩阵。E1(Hi)为高频能 量，E2(Hi)为低频能量，m= 1、2、3...^^为时间窗长度；(4)求出高频能量与低频能量的最大比值k 式中，E1(Hi)为高频能量，E2(Hi)为低频能量，k为在所取的时间窗内高频能量与低 频能量的比值的最大值；(5)区内外故障的甄别判据，当k彡1时为区外故障，当k> 1时为区内故障。以下是本发明的设计原理1.直流输电线路单端超高速保护原理现有的直流线路保护主要由ABB和SIEMENS两家公司提供的，是根据电压行波突 变量、电压行波变化率和电流变化梯度构成快速的单端量保护的判据。但实际运行经验表 明当直流线路经高阻故障时，行波保护的电压变化率减小，导致行波保护拒动。直流输电 系统包括直流输电线路和直流线路两端的平波电抗器和直流滤波器，它们构成了直流输电 线路高频暂态量的“天然”边界，并且边界频率特性分析表明线路边界内、外故障信号的高 频分量存在显著的特征差异，据此可以构成单端暂态量保护元件。利用高频分量特征的直 流输电线路单端超高速保护原理图如图1所示。2.边界元件特高压直流输电系统结构图如图2所示。图2中，送电容量为5000MW，整流侧和逆 变侧的无功补偿容量分别为3000Mvar和3040Mvar ；每极换流单元由2个12脉冲换流器串 联组成，直流输电线路为六分裂导线，全长取为1500km，采用J. R. Marti频率相关模型；线 路两侧装有400mH的平波电抗器；M点为保护安装处。本发明创造性地提出了以平波电抗器和直流滤波器构成其物理边界，并对其幅频 特性进行了分析。如图3所示，其中U1为区外暂态电压，U2为U1经边界传变至直流线路保 护安装处的电压；Bi、B2、B3、B4为直流滤波器避雷器,D1为平波电抗器避雷器、D2为直流母线 避雷器，本发明将避雷器I、B2、B3、B4、Dp D2统称为边界避雷器。Bp B2额定电压为150kV， B3、B4额定电压为75kV，D^ D2额定电压为824kV。
图 3 中，L = 400mH、L1 = 39. 09mH、L2 = 26. 06mH、L3 = 19. 545mH、L4 = 34. 75mH、 C1 = 0. 9 μ F、C2 = 0. 9 μ F、C3 = 1. 8 μ F、C4 = 0. 675 μ F。现定义边界元件的传递函数H(j ω)为
(1)Z1(Jw)为直流滤波器阻抗，Z2(j )为平波电抗器阻抗。边界元件传递函数 Η(」ω)的幅频特性如图4。从图4可知当f < IOOOHz时，H(jco) ^ 1 ；当1000Hz < f
<2000Hz 时，H(jco)频谱曲线有振荡，特别是 f = 600Hz, f = 1200Hz, f = 1800Hz 时， H(jco) << 1 ；当f > 2000Hz时，H(jco) < 1，因此本发明选择f彡2000Hz为高频分量，f
<2000Hz为低频分量。3. S变换的基本理论S变换是一种可逆的局部时频分析方法，其基本思想是对连续小波变换和短时傅 里叶变换的发展。信号x(t)的S变换S( τ，f)定义如下 式⑵中 式⑵和式(3)中，w( τ -t，f)为高斯窗口 ； τ为控制高斯窗口在时间轴t的位 置参数；f为频率；j为虚数单位。对式(1)右边先作传统的傅里叶变换，再作傅里叶反变换，最后进行变量代换将S 变换转换成信号x(t)的傅里叶变换X(f)的函数，即 式(4)中，f兴0。这样，S变换就可以利用FFT实现快速计算。由式⑷可以得 到信号X(k)的S变换的离散表示形式S[m，η]为 于是对采集到的N个离散信号点X[k](k = 0，l，...，N-l)采用式(4)、(5)进行S 变换，变换结果为一复时频矩阵，记作S矩阵，其行对应采样时间点，列对应频率，相邻行之 间的频率差Af为 式(8)中，fs为采样频率，N为采样点数。第η行对应的频率4为
(9)4.极波的定义
直流线路发生故障后，启动元件立即启动。根据保护安装处测得的两极直流电压、 直流电流，求出正极线的极波P1GO,负极线的极波P2(k)为 式中，Zp为直流输电线路极波阻抗，U1 (k)为正极直流电压，ii(k)为正极直流电 流，u2(k)为负极直流电压，i2(k)为负极直流电流，k= 1、2、3....N，N为时间窗长度。5.基于S变换的区内、区外故障的甄别 图1所示系统在距离保护安装处IOOkm处，0. 505s发生正极接地故障，极波波形如 图5 (a)所示；故障过渡电阻为0. 1 Ω，时间窗长度选取故障后2ms，采样频率为100kHz。
本发明中，对极波电压进行S变换，采样数据长度200，经过S变换得到101 X 200 的复矩阵。η = 1对应直流分量，η = 101对应频率50kHz的极波电压的高频分量，相邻频 率的间隔为500Hz，S[m, η]矩阵中的列向量表示信号在某一时刻的幅频特性，其行向量表 示信号在某一频率下的时域分布。根据边界元件的幅频特性，选择η = 5为2000Hz，作为高频分量和低频分量的分界 频率。η ^ 5为高频分量，η < 5为低频分量，将式(12)中的各个元素求模后，再求出高频 能量和低频能量。 式(12)、(13)中，E1为高频能量，E2为低频能量。从图5可以看出区外故障时，由于边界元件对高频的衰减作用，保护安装处测到 的正极线极波的高频含量远小于低频含量。区内故障时，高频分量没有经过边界元件，因此 高频含量大于低频含量。定义高频能量与低频能量的最大比值k 因此，提出区内、区外故障判据k彡1，为区外故障(15a)k>l，为区内故障(15b)本发明与现有技术相比具有如下优点1、本方法采样频率为100kHz，时间窗为2ms，避开了直流控制系统的暂态响应过 程，结论更为精确。2、本方法不受噪声干扰和故障距离的影响，且耐受过渡电阻的性能很强，有较强 的实用性。3、本方法采用极波电压进行分析，既可以作为区内外故障判据的特征量，又可以 作为故障选极判据的特征量。
图1为直流系统输电线路单端超高速保护原理框图，图中U1U1为保护安装处获得 的直流电压和直流电流。图2为云广士SOOkV直流输电系统结构图，图中F2、F3为区外故障，FpF4为区内 故障，M为保护安装处。图3为平波电抗器和直流滤波器构成的边界元件，U1为区外的暂态电压，U2为U1 经边界传变至直流线路保护安装处的电压；Bp B2, B3、B4为直流滤波器避雷器山工为平波电 抗器避雷器、D2为直流母线避雷器；Lp L2, L3> L4为电感元件；Cp C2、C3> C4为电容元件。图4为边界元件的频谱特性图，f为频率，Hz为频率的单位，H(j ω)为频谱的幅值。图5为极波波形图，图中t/s为时间/秒，u/kV为电压/千伏。(a)为区内故障时 的正、负电压；(b)为区外故障时的正、负极波。图6为本发明的正极极波高频分量与低频分量比值的分布图，图中k为高频分量 与低频分量的比值，t/s为时间/秒。(a)为区内故障时，正极波高频分量与低频分量的比 值；(b)为区外故障时，正极波高频分量与低频分量的比值。
具体实施例方式仿真模型如图1所示，距离保护安装处IOOkm处分别发生正极接地故障，故障过渡 电阻为0. 1Ω，时间窗长度取2ms，采样频率为100kHz。(1)直流线路发生故障后，启动元件立即启动，根据公式
( ο)
)(11)求出正极线极波和负极线极波，极波波形如图(4)所示；(2)根据公式 对正极线极波进行S变换，得到101 X 200的复矩阵；(3)对复矩阵中的各元素求模，根据边界元件的幅频特性，选择η = 5为2000Hz，
作为高频分量和低频分量的分界频率。η ^ 5为高频分量，η < 5为低频分量，根据公式 101(13) 求出高频能量和低频能量。(4)求出高频能量与低频能量的比值k = 4. 7，根据k > 1 (公式(15b))，为区内故 障，判断为区内故障。本发明中对不同的故障距离、不同的接地电阻进行了仿真验证，得到正极波电压 高频能量与低频能量的比值k，结果如下表所示。
权利要求
一种利用极波电压S变换能量比值的特高压直流输电线路故障边界元件的方法，其特征在于按以下步骤进行(1)直流线路发生故障后，启动元件立即启动，根据保护安装处测得的两极直流电压、直流电流，求出正极线极波电压P1(k)，负极线极波电压P2(k)为P1(k)＝Zp×i1(k)-u1(k) (1)P2(k)＝Zp×i2(k)-u2(k) (2)式中，Zp为直流输电线路极波阻抗，u1(k)为正极直流电压，i1(k)为正极直流电流，u2(k)为负极直流电压，i2(k)为负极直流电流，k＝1、2、3....N，N为采样序列长度；(2)对极波电压进行S变换，其采样频率为100kHz，采样序列长度为200，经S变换得到101×200的复矩阵 <mrow><mi>X</mi><mo>[</mo><mi>n</mi><mo>]</mo><mo>=</mo><mfrac> <mn>1</mn> <mi>N</mi></mfrac><munderover> <mi>Σ</mi> <mrow><mi>k</mi><mo>=</mo><mn>0</mn> </mrow> <mrow><mi>N</mi><mo>-</mo><mn>1</mn> </mrow></munderover><mi>x</mi><mo>[</mo><mi>k</mi><mo>]</mo><msup> <mi>e</mi> <mrow><mo>-</mo><mi>j</mi><mn>2</mn><mi>πkn</mi><mo>/</mo><mi>N</mi> </mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo></mrow> </mrow> <mrow><mi>S</mi><mo>[</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>]</mo><mo>=</mo><munderover> <mi>Σ</mi> <mrow><mi>k</mi><mo>=</mo><mn>0</mn> </mrow> <mrow><mi>N</mi><mo>-</mo><mn>1</mn> </mrow></munderover><mi>X</mi><mo>[</mo><mi>n</mi><mo>+</mo><mi>k</mi><mo>]</mo><msup> <mi>e</mi> <mrow><mo>-</mo><mn>2</mn><msup> <mi>π</mi> <mn>2</mn></msup><msup> <mi>k</mi> <mn>2</mn></msup><mo>/</mo><msup> <mi>n</mi> <mn>2</mn></msup> </mrow></msup><msup> <mi>e</mi> <mrow><mi>j</mi><mn>2</mn><mi>πkm</mi><mo>/</mo><mi>N</mi> </mrow></msup><mo>,</mo><mi>n</mi><mo>≠</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo></mrow> </mrow> <mrow><mi>S</mi><mo>[</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>]</mo><mo>=</mo><mfrac> <mn>1</mn> <mi>N</mi></mfrac><munderover> <mi>Σ</mi> <mrow><mi>k</mi><mo>=</mo><mn>0</mn> </mrow> <mrow><mi>N</mi><mo>-</mo><mn>1</mn> </mrow></munderover><mi>x</mi><mo>[</mo><mi>k</mi><mo>]</mo><mo>,</mo><mi>n</mi><mo>=</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo></mrow> </mrow>式中x[k]为采集到的N个离散的正极波信号点，k＝0、1、1....N-1，N为采样序列长度。X[n]为x[k]的傅里叶变换。S[m，n]为复时频矩阵，其列向量表示信号在某一时刻的幅频特性，其行向量表示信号在某一频率的时域分布；(3)对得到的n+1行m列的复时频矩阵中的各个元素求模，得到模矩阵，根据下式计算高频能量和低频能量 <mrow><msub> <mi>E</mi> <mn>1</mn></msub><mrow> <mo>(</mo> <mi>m</mi> <mo>)</mo></mrow><mo>=</mo><munderover> <mi>Σ</mi> <mrow><mi>n</mi><mo>=</mo><mn>5</mn> </mrow> <mn>101</mn></munderover><mo>|</mo><mi>S</mi><mo>[</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>]</mo><mo>|</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo></mrow> </mrow> <mrow><msub> <mi>E</mi> <mn>2</mn></msub><mrow> <mo>(</mo> <mi>m</mi> <mo>)</mo></mrow><mo>=</mo><munderover> <mi>Σ</mi> <mrow><mi>n</mi><mo>=</mo><mn>1</mn> </mrow> <mn>4</mn></munderover><mo>|</mo><mi>S</mi><mo>[</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>]</mo><mo>|</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo></mrow> </mrow>式中|S[m，n]|为将S矩阵中的各个元素求模后得到的模矩阵。E1(m)为高频能量，E2(m)为低频能量，m＝1、2、3....N，N为时间窗长度；(4)求出高频能量与低频能量的最大比值k <mrow><mi>k</mi><mo>=</mo><mi>max</mi><mrow> <mo>(</mo> <mfrac><mrow> <msub><mi>E</mi><mn>1</mn> </msub> <mrow><mo>(</mo><mi>m</mi><mo>)</mo> </mrow></mrow><mrow> <msub><mi>E</mi><mn>2</mn> </msub> <mrow><mo>(</mo><mi>m</mi><mo>)</mo> </mrow></mrow> </mfrac> <mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo></mrow> </mrow>式中，E1(m)为高频能量，E2(m)为低频能量，k为在所取的时间窗内高频能量与低频能量的比值的最大值；(5)区内外故障的甄别判据，当k≤1时为区外故障，当k＞1时为区内故障。
全文摘要
本发明是一种利用极波电压S变换能量比值的特高压直流输电线路故障边界元件的方法。属电力系统继电保护技术领域。本发明为当直流线路发生故障后，启动元件启动，根据保护安装处测得的两极直流电压、直流电流，求出正负极线极波电压。选取故障后采样序列长度为200的离散极波电压信号，进行S变换，变换结果为一101×200的时频复矩阵，并对此复矩阵中的各个元素求模。根据得到的模矩阵，提取极波电压的高频分量和低频分量后，求出极波电压高频能量与低频能量的比值。根据高频能量与低频能量的比值的大小，区分区内、区外故障。大量仿真结果表明，本发明效果良好。
文档编号H02H7/26GK101860020SQ20101019190
公开日2010年10月13日 申请日期2010年6月4日 优先权日2010年6月4日
发明者彭仕欣, 束洪春, 田鑫萃 申请人:昆明理工大学