小波变换在电力系统扰动检测中的应用 - 范文中心

小波变换在电力系统扰动检测中的应用

08/03

《现代电子技术》2003年第17期总第160期仿真与测试

小波变换在电力系统扰动检测中的应用

王润怀, 辛云宏

1

2

(1. 陕西秦岭(秦华) 发电有限责任公司 陕西华阴 714206; 2. 西安电子科技大学电子工程学院 陕西西安 710071)

摘 要:利用小波变换的时频局域化特征, 将电力系统中常见的4种扰动(电压突降、电压突升、瞬间间断、瞬时振荡) 信号波形的时间和频率信息以可视的形式表现出来。计算机仿真结果表明:所提出的方法能同时表示出4种扰动的时间和频率信息, 从而能有效的检测到这些变化的起始和终止时刻, 可广泛应用于多种扰动的定位及分析。

关键词:小波变换; 电力系统; 扰动检测; 频率

中图分类号:T M 712     文献标识码:B      文章编号:1004373X (2003) 1702903

Application of Wavelet Transform to Detect Power System Disturbances

WA N G R unhuai 1, XIN Y unho ng 2

(1. Shaanx i Qinli ng (Q i nhua ) Electri c Power Co. Lt d. , Hua yin, 714206, China ; 2. Colleg e of El ect ri cs Eng ineering , Xi dian Uni v ersit y, Xi ′an, 710071, China)

Abstract :A w av elet transfor m appro ach is pro po sed t o detect pow er system disturbances in this paper . W ith the t ime

fr e-

quency localizatio n char acter istics embedded in wav elet s, t he time and fr equency infor mation o f t he w avefo rm o f fo ur po wer sy stem disturbances including vo lta ge sag , vo ltage sw ell, momentar y int er r uption and o scillato ry is pr esented as a visualized scheme. R e-sults of co mputer simula tio n illust rate t he pr acticality and advantag es of the pro po sed metho d fo r t he a pplicatio n .

Keywords :w avelet tr ansfor m; pow er sy stem; distur bance detect; fr equency

1 引 言

随着高压非线性电子控制设备在供电线路上的广泛应用, 线路电压波形失真度的大小在今天看来更为重要。该失真对供电质量的影响越来越受到供电部门及用户的关注。从经济利益考虑, 有效地监视这些电力系统扰动, 对于供电部门和用户都是很有价值的。

传统的处理电力系统扰动信号的方法是基于离散傅里叶变换的方法, 这种方法的主要缺点是:

(1) 不能同时进行时间和频率分析。(2) 不能应用于非平稳信号的分析。

为了克服以上缺陷, 人们从很多方面对傅里叶变换做了改进, 其中最有名的就是窗口傅里叶变换, 窗口傅里叶变换虽然在很多方面改变了傅里叶变换的缺陷, 但是他对于信号来说相当于一个带宽固定的带通滤波器, 不能根据不同的需要对信号做多尺度的分析, 而小波变换的方法则可以通过伸缩和平移等运算功能对函数或信号进行多尺度细化分析, 解决了傅里叶变换不能解决的许多难题。将小波变换方法应用于检测电力系统的扰动, 具有许多优点:

(1) 在检测信号随时间的动态变化时非常有效, 在信号变化快速的时间间隔内, 该方法能放大感兴趣的区域, 以获得一个更清晰的信号特征。

(2) 能同时表示出4种扰动的时间和频率信息。(3) 可广泛应用于多种扰动的定位, 在那些对供电质量要求较高的场合, 本方法是传统的扰动检测方法的更新换代方法。2 小波变换的简单描述2. 1 小波变换的基本概念

具有有限能量的函数f (t ) [f (t ) ∈L 2(R ) ]的小波变换定义为以函数族7a , b (t ) =分核的积分变换, 如式(1) 所示:

W f (a , b ) =

7() 为积

a a

∫f (t ) 7(t ) d t

=∫f (t )

7() d t

a a

+∞

a , b

-∞+∞

-∞

(1)

其中:a 是尺度参数; b 是位移参数;

函数7a , b (t ) 称为小波。

根据Parseval 定理W f (a , b ) 可以写成如下形式:  W f (a , b ) =

收稿日期:20030617

∫f d (X ) 7-∞+∞

a , b

(X ) d X (2)

王润怀等:小波变换在电力系统扰动检测中的应用

其中:f d (X ) , 7d (X ) 分别为f (t ) , 7a , b (t ) 的傅里叶变换。定义时间域的质心或一阶原点距为:

   t 0

t û7=

∫û7

-∞-∞+∞-∞+∞

a , b

恒Q 带通滤波器。

(t ) û2d t

2

a , b

(t ) ûd t

1/2

(3)

时间域的标准偏差R 7a , b (即二阶中心距的平方根) 为:

  R 7a , b =

(t -t ) û7(t ) ûd t

∫û7(t ) ûd t

2

a , b

2

-∞

a , b

2

2

1, 0

(4)

图1 小波变换分辨元胞的空间表示

令:t =

∫û7

+∞-∞

+∞

t û71, 0(t ) ûd t -∞

(t ) ûd t

2*

, 则有:

(5) (6)

2. 2 小波函数

在小波理论中, 常用的小波基函数有好多种, 每一个基函数都有其特点, 本文将Daubechies (db2) 小波应用于电力系统扰动的检测, 下面对其做一简要介

绍。

Daubechies 函数是由世界著名小波分析学者In-rid Daubechies 构造的小波函数, 除了db1(即Haar 小波) 外, 其他小波没有明确的表达式, 但转换函数h 的平方是很明确的。

N -1

    t 0=at +b 同时还可推导出:

    R 7a , b =a R 71, 0

d

同样, 在频域运用相同的定义形式, 规定X 07为a , b

频域的质心, R 7d a , b 为标准偏差, 可推导出下列的关系式:

d

X 07

1, 0

(7)    X =a

R 7d 1, 0

d

   R 7a , b =(8)

a

d

由此可以得出:在此空间中, 围绕点心(t 0, X 07) a , b

d 7a , b

d X 07

1, 0

的子波变换的分辨元胞由[(t 0±a R 71, 0) ×(±

a

d R 07

1, 0

]确定, 如图1所示。该元胞的尺度大小随尺度参a

õ假设P (y ) =式的系数, 则有:

∑C

k =0

N -1+k k

-1+k

y k , 其中C N k 为二项

  ûm 0(X ) û2=(cos 2

其中:

2N -1

N

) P (sin 2)

22

(11)

数a 而变化, 但是分辨元胞的面积却与尺度参数a 及

定位参数b 无关。换句话说, 相空间中分辨元胞的宽度(时间窗口宽度) 与高度(频率窗口的高度) 在各处是不一样的, 即他是尺度参数a 的函数, 以致使元胞的面积由Heisenberg 测不准原理限定为一常数, 从而在时间域中高频分量的时间局域化分辨率的提高是以频率局域化中由R 7d a , b 测定的不确定性增大为其代价的, 分析检测高频分量时(a 减小) , 时间窗口自动变窄, 频率窗口自动变高; 分析检测低频分量时(a 增大) , 时间窗口自动变宽, 频率窗口自动变窄, 实现了时-频窗口的自适应变化。从滤波的观点来看, 7a , b (t ) 的频谱

d

7d a , b (X ) 具有带通特性, 中心频率X 0=X 07, 带宽BW

a , b

   m 0=h k e -ik X (12) ∑2k =0

õ小波函数7和尺度函数U 的有效支撑长度为2N -1, 小波函数7的消失距阶数为N 。

õdb N 大多数不具有对称性。

õ正则性随着序号的增加而增加。õ函数具有正交性。3 扰动模拟及分析

下面我们将Daubechies (db 2) 小波应用于影响供电质量的4种情况进行研究。这4种现象包括:电压突降、电压突升、瞬间间断、瞬时振荡, 从分析结果可以看出, Daubechies 小波能有效的检测到这些变化的起始和终止时刻。3. 1 电压突降

电压突降指的是系统电压下降了10%~90%, 并且持续时间为0. 5个周期到1min 。这一现象有可能是由故障电流、重负载或是大电动机的接通所引起的, 当电压下降了30%或更多时就认为情况很严重。图2

=2R 7d a , b , 这样, 该带通滤波器的品质因数Q 则为:X 7d a , b a X

7d 1, 0X 7d 1, 0X 0

====CONST (9) Q =

BW 2R 7d a , b 2a R 7d 1, 02R 7d 1, 0

Q 为恒常数的重要意义在于:子波变换时频窗口

的自适应特性相当于恒Q 滤波技术, 子波相当于一个

《现代电子技术》2003年第17期总第160期(a ) 为运用小波变换对该情况处理的结果。图中的原始信号表示突降的电压扰动信号, a 1、a 2分别表示一次、二次分解的低频系数, d 1、d 2分别表示一次、二次分解的高频系数, 图2均采用这种表示形式。从图中可以看出, 经过二次分解, 其二次分解的高频系数在该事件的起始和结尾都有一个大的变化, 而其他地方的信号则很小, 扰动的出现很容易看出来。3. 2 电压突升

电压突升的情况与电压突降的情况刚好相反, 指的是系统电压暂短的上升, 这一现象经常出现在三相电路中单相短路时的正常相上。电压突升有可能给连接其上的设备造成永久性损伤。图2(b ) 为小波变换对该情况处理的结果。从图中可以看出, 经过二次分

仿真与测试

解, 在电压突升信号的起始和结尾都有一个大的变化, 扰动表现的很明显。3. 3 瞬间间断

瞬间间断是指供电系统上电压的瞬间丢失, 当系统电压的下降率为90%~100%, 并且持续时间为0. 5个周期到1m in 时, 则认为系统处于瞬间间断状态。这一现象常常是由电路故障引起开关动作产生的。小波变换对该情况处理的结果如图2(c ) 所示。3. 4 瞬时振荡

如果扰动的持续时间小于电压突降或突升的持续时间, 则认为扰动是瞬时的, 若瞬时扰动呈现振荡特性则称其为瞬时振荡, 瞬时振荡通常是电容器的开关引起的。小波变换对该情况处理的结果如图2(d)

所示。

图2 小波变换的仿真结果

4 结 语

从模拟测试结果可以看出, 小波变换的方法在检测和定位电力系统的这4种不同类型的扰动时是非常成功的, 将小波变换的方法运用于电力系统供电质量的监视是可行的。在模拟测试过程中, 没有考虑现场的各种干扰, 关于存在干扰分析结果如何, 还有待于进一步的研究。

参 考 文 献

[1] 宋国乡, 王世儒, 甘小冰. 数值分析[M ]. 西安:西安电子科技大学, 1998.

[2] 赵松年, 熊小芸. 子波变换与子波分析[M ].

北京:电子工业出版社, 1996.

(下转第34页)   

王西芳等:限制部分客户机上网的实现方法

当某些客户机被禁止一段时间后, 由于其他原因要求访问某个站点时, 可以在Black List 中选定该客户机的IP, 然后点击M odify (修改) , 在右边的Exter-nal IP (外部IP) 中指定Internet 网站的IP 地址, 点击OK 即可。这样该客户机即能通过Syg ate 上网。但是只能浏览你所指定的网址, 不能浏览其他网页。如果对该客户机无任何限制时, 可以将他从Black List 删除掉, 以恢复他的上网访问权。

2. 2 创建White List , 限制客户机的访问权限

White List 是局域网中被允许上网的一组IP 地址列表, 这里的IP 地址也是指在局域网中为各客户机配置的内部IP 地址。

当局域网中需要限制上网的客户机大于要上网的客户机时, 使用Black List 比较麻烦, 这时可以通过创建White List (白名单) 来实现。具体方法中:在Perm is-, sions Edito r 窗口, 单击White List 页面, 选择“Add ”出现Add BWList Item(添加黑白名单) 窗口(同图2) , 下面的操作步骤与黑名单相同, 不同之处是, 在这里所输入的Internal IP 地址, 表示允许上网的客户机的内部IP, Ex ternal IP 地址表示禁止浏览网站的IP 地址。

在创建White List 时, 计算机房管理员可根据自己实际情况, 按照创建Black List 时同样的方法选择协议类型、端口类型, 设定访问网站的时间段。同样当需要限制的客户机全部设置完毕后, 必须将Activate Black List 打上勾, 以使Black List 生效。3 结 语

通过Sy gate 的设置, 机房管理员就可以轻松、方便地管理学生机房, 同时净化了学生机房网络环境, 拒绝垃圾信息, 更好地为学生学习掌握计算机知识服务, 并对学生身心的健康发展将发挥更大的作用。

参 考 文 献

[1] 林育深. 陈德. 学生机房上网管理的几种有效

方法[J ]. 广东教育学院学报, 2002, 22(2) :

9394.

[2] 杨金胜, 李晓霞. 局域网维护技巧[J].邯郸农业高等专科学校学报, 2002, 19(1) :2526. [3] 刘军. 基于局域网的机房自动化管理系统[J ].

南京建筑工程学报, 2002, 61(2) :6367.

作者简介 王西芳 女, 1965年出生, 毕业于西安交通大学计算机科学与工程系, 计算机工程师。主要从事计算研究与应用

工作。

(上接第28页) 4 结 语

需要指出的是, M atlab 6. 5对应的神经网络工具箱还提供很多其他网络模型、算法函数, 本文主要以BP 网络模型, 利用Levenberg M arquardt 规则训练前向网络为例介绍了工具箱的应用方法, 其他网络模型和算法的应用可以此类推。基于该工具箱开发的设备状态分类器有利于复杂设备故障诊断技术水平的不

断提高。

参 考 文 献

[1] 温熙森, 胡茑庆, 邱静. 模式识别与状态监控

[M]. 长沙:国防科技大学出版社, 1997.

[2] 飞思科技产品研发中心. Matlab 6. 5辅助神经

网络分析与设计[M ]. 北京:电子工业出版社,

2003.

作者简介 丁晓伟 女, 1978年出生, 东南大学仪器科学与工程系硕士研究生。主要研究方向为设备故障诊断。

(上接第31页)

[3] 胡昌华, 等. 基于M ATLAB 的系统分析与设计

[M ]. 西安:西安电子科技大学出版社, 1999. [4] 孙才新, 李新, 扬永明. 从白噪声中提取局部

放电信号的小波变换方法研究[J]. 电工技术学报, 1999, 14(3) :4750.

[5] 陈继东. 小波分析应用与在线监测中信噪分离

的研究[J ]. 电网技术. 1999, 23(11) :54

57.

[6] 吕柏劝, 李天铎. 一种基于小波网络的故障检测

方法[J]. 电工技术学报, 1997, 12(4) :34

36.

作者简介 王润怀 男, 1966年出生, 工程师。主要从事通讯及电力电子技术的研究工作。


相关内容

  • 电气自动化技术 毕业论文
    电气自动化技术 目 录 摘要„„„„„简要介绍电气自动化技术的概念及其包括的专业知识 关键字„„„„„„„„„„„„控制.系统.检测.网络化 前言„„„„„„„„„„„„„„„„„„„„„„„„1 第一章 自动控制统„„„„„„„„„„„„ ...
  • 高次谐波及其抑制措施
    高次谐波及其抑制措施 字号显示:大 中 小 2006-05-23 09:21:47 来源:电子查询网 以前,人们基本上只根据电压的幅值和周波的稳定性来衡量电能的质量.近年来,随着工业生产的发展,尤其是冶金.化工产品的开发.电气铁道系统的不断 ...
  • 关于计量经济学模型随机扰动项的讨论
    第26卷第2期 2009年2月 统计研究 Statistical Research V01.26.No.2 Feb.2帅9 关于计量经济学模型随机扰动项的讨论' 李子奈李鲲鹏 内容提要:论文指出了计量经济学模型中源生的随机扰动项和衍生的随机 ...
  • 一种新型的电梯能量回馈并网系统
    计算机系统应用 http://www.c-S-&org.ca 2012年第2l卷第3期 一种新型的电梯能量回馈并网系缈 彭继慎,王伟伟,宋立业 (辽宁工程技术大学电气与控制工程学院,葫芦岛125105) 摘要:针对普通电梯变频器不能 ...
  • 激光功率计doc
    概 况 H14490是一种新型宽光谱响应.高灵敏.快响应.低温漂,数字直读的激 光功率测量仪器.该仪器配套IP - 200.IP - 1000系列探测器.分别适用最大量 程2W .10W .50W .200W 和分辨率l u w.10u w ...
  • 三相四桥臂逆变器控制技术研究_顾和荣
    第39卷 第24期 电力系统保护与控制 Vol.39 No.24 2011年12月16日 Power System Protection and Control Dec.16, 2011 三相四桥臂逆变器控制技术研究 顾和荣, 王德玉, 沈 ...
  • 信息科学技术概论考试大全
    第一章 1.通信系统包括哪些组成部分,能否举例说明这些部分及其作用? 答:信源:作用是把各种消息转换成原始电信号,称之为消息信号或者基带信号. 例如电话机.摄像机.电传机和计算机等各种数字终端设备.发送设备:基本功 能是将信源和信道匹配起来 ...
  • 城市轻轨再生制动能量吸收的仿真研究
    统 仿 真 学 报 V ol. 21 No. 15 2009年8月 Journal of System Simulation Aug., 2009 第21卷第15期 系 城市轻轨再生制动能量吸收的仿真研究 李 政,潘孟春,胡 楷 (国防科技 ...
  • 自动控制理论答案(孙扬声版)
    T2-1 判断下列方程式所描述的系统的性质:线性或非线性,定常或时变,动态或静态. 1dytd2yt2: (4)sinty(t)3u(t): ytu(t)ut(1)2yt3t: (3)2 dtdt 2 ( ...
  • 双闭环V-M调速系统设计
    电力电子与电机拖动 综合课程设计 题 目: 双闭环V -M 调速系统设计 专 业: ***********班 学 号: ********** 姓 名: AUT 完成日期: **** 指导教师: aut 电力电子与电机拖动综合课程设计任务书 ...