一种基于小波能量比值的逆变站区内外故障识别方法与流程

1.本发明属于逆变站区内外故障识别技术领域，尤其涉及一种基于小波能量比值的逆变站 区内外故障识别方法。


背景技术：

2.换流站是连接交直流系统的重要组成部分，换流站保护用以保障站内设备免受故障威胁， 在交直流混联系统保护中起到至关重要的作用。为了防止在逆变站区内发生故障时逆变站设 备损坏，实际工程为逆变站配置了多类保护，然而现场的运行经验表明，逆变站区外交流系 统故障导致逆变站保护误动作的事件时有发生，严重时将引起直流系统闭锁，威胁电网安全 稳定，即逆变站的部分保护缺乏对于逆变站区内外故障的辨识能力。因此，为防止逆变站保 护在区外故障时误动作，亟需引入一种逆变站区内外故障识别方法，赋予逆变站保护识别逆 变站区内外故障的能力。
3.

技术实现要素：

4.为了能准确、可靠地识别出逆变站区内外故障，本发明提供了一种基于小波能量比值的 逆变站区内外故障识别方法。
5.本发明的一种基于小波能量比值的逆变站区内外故障识别方法，包括以下步骤：
6.a、数据采集
7.以采样频率f
s
＝10khz实时采集电流数据：实时采集逆变站的n个换流变压器tb1～tb
n
一次侧的三相电流信号，分别为i
a_tb1
(t)、i
b_tb1
(t)、i
c_tb1
(t)，i
a_tb2
(t)、i
b_tb2
(t)、i
c_tb2
(t)，
…
， i
a_tbn
(t)、i
b_tbn
(t)、i
c_tbn
(t)；实时采集交流滤波器出口处的三相电流信号，分别为i
a_acf
(t)、i
b_acf
(t)、 i
c_acf
(t)，其中t为采样时刻。
8.b、故障分量的模量计算
9.b1、分别计算换流变压器一次侧a相、b相以及c相电流之和i
a_tb
(t)、i
b_tb
(t)、i
c_tb
(t)， 算式为
10.b2、计算i
a_tb
(t)、i
b_tb
(t)、i
c_tb
(t)的故障分量δi
a_tb
(t)＝i
a_tb
(t)
‑
i
a_tb
(t
‑
t)、δi
b_tb
(t)＝i
b_tb
(t)
‑ꢀ
i
b_tb
(t
‑
t)、δi
c_tb
(t)＝i
c_tb
(t)
‑
i
c_tb
(t
‑
t)；计算i
a_acf
(t)、i
b_acf
(t)、i
c_acf
(t)的故障分量δi
a_acf
(t)＝i
a_acf
(t)
‑
i
a_acf
(t
‑
t)、δi
b_acf
(t)＝i
b_acf
(t)
‑
i
b_acf
(t
‑
t)、δi
c_acf
(t)＝i
c_acf
(t)
‑
i
c_acf
(t
‑
t)；其中t表示 一个交流工频周期。
11.b3、对δi
a_tb
(t)、δi
b_tb
(t)、δi
c_tb
(t)进行相模变换得到换流变压器一次侧a相、b相以及c 相总电流故障分量的模量，为δi
tb
(t)＝[2
×
δi
a_tb
(t)
‑
δi
b_tb
(t)
‑
δi
c_tb
(t)]/3；对δi
a_acf
(t)、δi
b_acf
(t)、 δi
c_acf
(t)进行相模变换得到交流滤波器出口处电流故障分量的模量，为 δi
acf
(t)＝[2
×
δi
a_acf
(t)
‑
δi
b_acf
(t)
‑
δi
c_acf
(t)]/3。
[0012]
c、故障时刻的确定
[0013]
c1、计算逆变站交流母线三相电压u
a
(t)、u
b
(t)、u
c
(t)的故障分量δu
a
(t)＝u
a
(t)
‑
u
a
(t
‑
t)、 δu
b
(t)＝u
b
(t)
‑
u
b
(t
‑
t)、δu
c
(t)＝u
c
(t)
‑
u
c
(t
‑
t)。
[0014]
c2、求δu
a
(t)、δu
b
(t)、δu
c
(t)绝对值的最大值δu
max
＝max(|δu
a
(t)|、|δu
b
(t)|、|δu
c
(t)|)；判 断δu
max
>δu
set
是否成立，若否，则返回步骤a；若是，则确定发生了故障，记录此时为故障 时刻，将此时刻记为算法启动时刻，并进入下一步；其中δu
set
为识别算法启动阈值。
[0015]
d、小波能量比值的求取
[0016]
d1、使用db1小波对信号δi
tb
(t)、δi
acf
(t)执行一层小波分解，产生的细节系数即为小波 变换系数，信号δi
tb
(t)的小波变换系数记为w
t
，信号δi
acf
(t)的小波变换系数记为w
f
；由小 波变换原理可知，信号δi
tb
(t)、δi
acf
(t)高频小波能量所对应的频率范围为2.5khz<f<5khz。
[0017]
d2、分别计算δi
tb
(t)、δi
acf
(t)的高频小波能量，其中j＝1 代表算法启动后的第一个采样点，n代表积分时间窗内采样点的个数。
[0018]
d3、计算δi
tb
(t)与δi
acf
(t)的高频小波能量比值e
f
/e
t
，定义k
e
＝e
f
/e
t
。
[0019]
e、逆变站区内外故障识别
[0020]
判断k
e
>k
set
是否成立，若是，则判断发生的故障为逆变站区外故障，若否，则判断发生 的故障为逆变站区内故障，其中k
set
为故障识别阈值。
[0021]
进一步的，识别算法启动阈值δu
set
的取值为0.01～0.1倍故障启动所用电压测点处电压的 额定值。
[0022]
进一步的，积分时间窗内采样点的个数n取值为10。
[0023]
进一步的，故障识别阈值k
set
取值为5。
[0024]
本发明与现有技术相比的有益技术效果为：
[0025]
1、本发明可实现对逆变站区内外故障的准确识别，在逆变站保护误动作之前识别出故障 发生在逆变站区内还是逆变站区外，即赋予了逆变站保护对于逆变站区内外故障的辨识能力， 为提高逆变站保护的可靠性打下理论基础。
[0026]
2、本发明仅需采集换流变压器一次侧三相电流信号和交流滤波器出口处三相电流信号， 无需额外安装测点，且测点距离逆变侧的控制保护系统较近，无需依赖于远距离通信来汇总 采集到的信号数据。
[0027]
3、本发明适用范围广，由于高压直流输电换流器工作需要大量的无功功率，则交流滤波 器须有无功补偿的作用，即交流滤波器高频阻抗很低，本发明即是利用此特点提出的逆变站 区内外故障识别方法，综上本发明所提方法适用于所有的高压直流输电工程。
[0028]
4、本发明在各种工况下具有良好的适应性，在不同故障类型、不同故障位置、不同过渡 电阻、不同故障初始角以及不同噪声干扰下，所提逆变站区内外故障识别方法仍能够可靠识 别出故障发生于逆变站区内还是逆变站区外，说明了本发明适应性强，可靠性高。
附图说明
[0029]
图1为高压直流输电系统逆变站区内外故障分布示意图。
[0030]
图2为f1处a相发生金属性接地故障时δi
tb
(t)与δi
acf
(t)的高频小波能量图。
[0031]
图3为f4‑
l1处a相发生金属性接地故障时δi
tb
(t)与δi
acf
(t)的高频小波能量图。
具体实施方式
[0032]
下面结合附图和具体实施方式对本发明做进一步的详细说明。
[0033]
本发明的一种基于小波能量比值的逆变站区内外故障识别方法，包括以下步骤：
[0034]
a、数据采集
[0035]
以采样频率f
s
＝10khz实时采集电流数据：实时采集逆变站的n个换流变压器tb1～tb
n
一次侧的三相电流信号，分别为i
a_tb1
(t)、i
b_tb1
(t)、i
c_tb1
(t)，i
a_tb2
(t)、i
b_tb2
(t)、i
c_tb2
(t)，
…
， i
a_tbn
(t)、i
b_tbn
(t)、i
c_tbn
(t)；实时采集交流滤波器出口处的三相电流信号，分别为i
a_acf
(t)、i
b_acf
(t)、 i
c_acf
(t)，其中t为采样时刻。
[0036]
b、故障分量的模量计算
[0037]
b1、分别计算换流变压器一次侧a相、b相以及c相电流之和i
a_tb
(t)、i
b_tb
(t)、i
c_tb
(t)， 则
[0038]
b2、计算i
a_tb
(t)、i
b_tb
(t)、i
c_tb
(t)的故障分量δi
a_tb
(t)＝i
a_tb
(t)
‑
i
a_tb
(t
‑
t)、δi
b_tb
(t)＝i
b_tb
(t)
‑ꢀ
i
b_tb
(t
‑
t)、δi
c_tb
(t)＝i
c_tb
(t)
‑
i
c_tb
(t
‑
t)；计算i
a_acf
(t)、i
b_acf
(t)、i
c_acf
(t)的故障分量δi
a_acf
(t)＝ i
a_acf
(t)
‑
i
a_acf
(t
‑
t)、δi
b_acf
(t)＝i
b_acf
(t)
‑
i
b_acf
(t
‑
t)、δi
c_acf
(t)＝i
c_acf
(t)
‑
i
c_acf
(t
‑
t)；其中t表示 一个交流工频周期。
[0039]
b3、对δi
a_tb
(t)、δi
b_tb
(t)、δi
c_tb
(t)进行相模变换得到换流变压器一次侧a相、b相以及c 相总电流故障分量的模量，为δi
tb
(t)＝[2
×
δi
a_tb
(t)
‑
δi
b_tb
(t)
‑
δi
c_tb
(t)]/3；对δi
a_acf
(t)、δi
b_acf
(t)、 δi
c_acf
(t)进行相模变换得到交流滤波器出口处电流故障分量的模量，为 δi
acf
(t)＝[2
×
δi
a_acf
(t)
‑
δi
b_acf
(t)
‑
δi
c_acf
(t)]/3。
[0040]
c、故障时刻的确定
[0041]
c1、计算逆变站交流母线三相电压u
a
(t)、u
b
(t)、u
c
(t)的故障分量δu
a
(t)＝u
a
(t)
‑
u
a
(t
‑
t)、 δu
b
(t)＝u
b
(t)
‑
u
b
(t
‑
t)、δu
c
(t)＝u
c
(t)
‑
u
c
(t
‑
t)。
[0042]
c2、求δu
a
(t)、δu
b
(t)、δu
c
(t)绝对值的最大值δu
max
＝max(|δu
a
(t)|、|δu
b
(t)|、|δu
c
(t)|)；判 断δu
max
>δu
set
是否成立，若否，则返回步骤a；若是，则确定发生了故障，记录此时为故障 时刻，将此时刻记为算法启动时刻，并进入下一步；识别算法启动阈值δu
set
的取值为0.01～0.1 倍故障启动所用电压测点处电压的额定值。
[0043]
d、小波能量比值的求取
[0044]
d1、使用db1小波对信号δi
tb
(t)、δi
acf
(t)执行一层小波分解，产生的细节系数即为小波 变换系数，信号δi
tb
(t)的小波变换系数记为w
t
，信号δi
acf
(t)的小波变换系数记为w
f
；由小 波变换原理可知，信号δi
tb
(t)、δi
acf
(t)高频小波能量所对应的频率范围为2.5khz<f<5khz。
[0045]
d2、分别计算δi
tb
(t)、δi
acf
(t)的高频小波能量，其中j＝1 代表算法启动后的第一个采样点，n代表积分时间窗内采样点的个数，取值为10。
[0046]
d3、计算δi
tb
(t)与δi
acf
(t)的高频小波能量比值e
f
/e
t
，定义k
e
＝e
f
/e
t
。
[0047]
e、逆变站区内外故障识别
[0048]
判断k
e
>k
set
是否成立，若是，则判断发生的故障为逆变站区外故障，若否，则判断发生 的故障为逆变站区内故障，其中k
set
为故障识别阈值，取值为5。
[0049]
仿真实验
[0050]
基于pscad/emtdc仿真平台上搭建受端连接有3条交流线路的
±
800kv交直流系统仿 真模型，模型中3条受端交流线路l1～l3的长度分别为200km、110km、80km。高压直流输 电系统逆变站区内外故障分布示意图如附图1所示，其中f1代表换流变压器一次侧故障，f2代表换流变压器二次侧故障，f3代表换流器直流侧出口故障，均为逆变站区内故障；f4‑
l1代 表交流线路l1上发生的故障，为逆变站区外故障。在模型上设置换流变压器tb1一次侧f1处 a相发生金属性接地故障，信号δi
tb
(t)与δi
acf
(t)经小波变换后的高频小波能量如附图2所示， 计算可得k
e
＝0.9862<5，可判断出故障为逆变站区内故障；在模型上设置交流线路l1上点f4‑
l1处a相发生金属性接地故障，故障距离，即距离换流母线的距离为10km，信号δi
tb
(t)与δi
acf
(t) 经小波变换后的高频小波能量如附图3所示，计算可得k
e
＝22.8290>5，可判断出故障为逆变 站区外故障。
[0051]
为验证所提故障识别方法在不同故障类型与不同故障位置下的适应性，分别在模型上设 置不同故障类型与不同故障位置的逆变站区内外故障，并利用本发明所提方法对故障进行识 别，得到的仿真结果如表1所示。表1中ag表示a相接地故障，abg代表ab两相接地故 障，abcg代表三相接地故障；表1中的故障距离代表f4‑
l1距离换流母线的距离；f1、f2、f3、 f4‑
l1的过渡电阻均为60ω，故障初始角均为0
°
。
[0052]
表1不同故障类型与不同故障位置下的仿真结果
[0053]
[0054]
根据表1中结果可知，当逆变站区内故障时，k
e
均小于5，此时判断发生的故障为逆变 站区内故障；当逆变站区外故障时，k
e
均大于5，此时判断发生的故障为逆变站区外故障。 因此可知，无论故障发生于逆变站区内还是逆变站区外，本发明均能够准确识别，与本发明 所要实现的目的一致。
[0055]
为验证所提故障识别方法在不同过渡电阻下的适应性，分别在模型上设置不同过渡电阻 的逆变站区内故障和逆变站区外故障，并利用本发明所提算法对故障进行识别，得到的仿真 结果如表2所示。以a相接地故障为例，其中f4‑
l1处发生的故障在距离换流母线10km处， 故障初始角均为0
°
。
[0056]
表2不同过渡电阻下的仿真结果
[0057][0058]
根据表2中结果可知，在逆变站区内发生不同过渡电阻下的故障时，k
e
均小于5，此时 判断发生的故障为逆变站区内故障；在逆变站区外发生不同过渡电阻下的故障时，k
e
均大于 5，此时判断发生的故障为逆变站区外故障。因此可知，无论故障发生于逆变站区内还是逆变 站区外，本发明均能够准确识别，与本发明所要实现的目的一致。
[0059]
为验证所提故障识别方法在不同故障初始角下的适应性，分别在模型上设置不同故障初 始角下的逆变站区内故障和逆变站区外故障，并利用本发明所提算法对故障进行识别，得到 的仿真结果如表3所示。以a相接地故障为例(过渡电阻设置为60ω)，f4‑
l1处发生的故障 在距离换流母线30km处。
[0060]
表3不同故障初始角下的仿真结果
[0061][0062]
根据表3中结果可知，在逆变站区内发生不同故障初始角下的故障时，k
e
均小于5，
此 时判断发生的故障为逆变站区内故障；在逆变站区外发生不同故障初始角下的故障时，k
e
均 大于5，此时判断发生的故障为逆变站区外故障。因此可知，无论故障发生于逆变站区内还 是逆变站区外，本发明均能够准确识别，与本发明所要实现的目的一致。
[0063]
为验证所提故障识别方法在噪声干扰下的适应性，在模型上设置不同程度噪声干扰下的 逆变站区内故障和逆变站区外故障，并利用本发明所提算法对故障进行识别，以a相接地故 障为例，且过渡电阻设置为60ω，f4‑
l1处发生的故障在距离换流母线60km处，故障初始角 为0
°
，得到的仿真结果如表4所示。
[0064]
表4噪声干扰下的仿真结果
[0065][0066]
根据表4中结果可知，在噪声干扰下，逆变站区内发生故障时，k
e
均小于5，此时判断 发生的故障为逆变站区内故障；逆变站区外发生故障时，k
e
均大于5，此时判断发生的故障 为逆变站区外故障。因此可知，无论故障发生于逆变站区内还是逆变站区外，本发明均能够 准确识别，与本发明所要实现的目的一致。