为充分利用短路后电流变化率剧增的故障信息,同时剔除噪声干扰影响,在提取故障特征(如求导)之前,应对信号进行平滑和消噪处理。
文献[2,4]介绍了低压系统早期故障检测的小波分析算法,对低压配电系统短路电流信号进行了研究,证实了小波分解算法对短路故障早期检测的有效性。
由参考文献[7]知低压线路的短路故障信号比较复杂,有可能会含有非平稳信号,虽然应用传统的傅里叶变换分析来提取故障特征量时,滤波效果好,算法简单,但直流衰减分量的存在,会受到很大的局限性,因而产生较大的误差,而且,由参考文献[19]知道采用传统的全周波傅氏算法虽然能实现故障检测,但检测时间还是比较长。
然而,小波变换则具有优良的时域局部化性能,适合短路故障信号的分析。
小波分析作为新的信号处理方法的出现,在这几年故障诊断中取得很大的进展,给故障诊断技术带来了新的生机和活力,因此,在本文中,将对一种经过形态学滤波后三次B样条小波变换,得到第四尺度细节分量的故障检测模型的有效性,通过MATLAB进行仿真论证。
1.4本文的主要工作
本文介绍了近年来兴起的小波变换这一数学工具,利用MATLAB实现小波变换在低压系统中短路故障的快速检测的仿真分析。
本论文的主要工作:
1)利用仿真方法对低压短路故障暂态过程进行分析研究,并探索出短路故障特征。
2)引入形态学小波概念,即将形态学滤波和三次B样条小波结合起来,并在MATLAB上仿真和验证形态学滤波后小波分解,得到第四尺度细节分量的模型,在实现短路故障快速检测的可靠性和可行性,并介绍了门限值的确定问题。
第2章MATLAB和小波介绍
2.1MATLAB/Simulink的介绍
MATLAB(MatrixLaboratory)是美国MathWorks公司开发的一款公认的最好的科学计算类的软件。
它具备高效率的数值计算和图形交互等仿真功能。
MATLAB具有简洁灵活的言语风格,不论是在命令窗口还是M文件,对于每一行的命令都能实时地在工作窗口看到即时的结果。
对M文件的编程语法要求不高,变量并不用像单片机编程要定义变量或声明等步骤。
在MATLAB中存在大量的数学及工程之类的函数可供调用,因此方便数值的计算。
还有其面向图形架构让使用者可以执行视觉数据分析,并划出比较高品质的图片。
强大的工具箱,包含了功能性工具箱和学科性工具箱,更是让各行各业的工作者所青睐。
完备的帮助系统,提供在线帮助,这些帮助主要是嵌在M文件中,即时性强,反应迅速,可以help命令来实现。
除此之外,还提供了演示软件demo。
总之MATLAB给人带来了最直观和最简洁的开发环境。
本文中对于课题的研究,主要用到了MATLAB的M文件编程和Simulink中的电力系统工具箱进行建模。
MathWorks公司MATLAB的5.2版本开始,就增加了电力系统工具箱,该工具箱是基于Simulink环境。
Simulink是MATLAB软件的扩展,它是Simu一词和表示仿真和Link一词表示连接组成,实现动态系统建模和仿真的一个软件包,可直接建立电路仿真模型,随意改变仿真参数,并得到及时的仿真效果,直观性强,并进一步省去编程的步骤,这也是它与MATLAB语言有比较大的不同。
实质上,Simulink的重点是提供一个直观模型的建立,其与用户交互接口是基于Windows的模型化图形输入,各个模块的内部编程都是自带的编好的,因此,利用编程,也能实现仿真图出现的效果。
在Simulink的仿真工具箱中有众多的元器件库,电力系统工具箱模块库是仿真电路的元件库。
它包括连接元件库(Connectors),电源库(ElectricalSources),基本元件库(Elements),元件库(ExtraLibrary),电机元件库(Machines),测量元件库(Measurements)和电力电子元件库(PowerElectronics)。
这些模块库包含了大多数常用电力系统元件的模块。
利用这些库模块及其它库模块,用户可方便、直观地建立各种系统模型并进行仿真。
下表2-1,对该模块中的主要组成和功能进行了比较详细的介绍。
各种模块的用法,以及参数设计都比较简单。
每个模块还提供了比较详实的帮助和相应器件的用法演示实例。
表2-1电力系统工具箱模块库的组成和相应的用途
电力系统工具箱模块库
主要构成
主要用途
电路元件模型
断路器(Breaker)、分布参数线(DistributeParameterLine)、线性变压器(LinearTransformer)、并联RLC负荷(ParallelRLCLoad),II型线路参数(IISectionLine)、饱和变压器(SaturableTransformer)、串联RLC支路(SeriesRLCBranch)、串联RLC负荷(SeriesRLCload)、过电压自动装置(SurgeArrester)。
这部分可以仿真交流输电线装置
电力电子设备模型
二极管(Diode)、GT0、理想开关(IdealSwitch)、MOS管(Mosfet)、可控晶闸管(Thyristor)的仿真模型。
这些设备模型不仅可以单独进行仿真而且可以组合在一起仿真整流电路等直流输变电的电力电子设备
电机设备模型
异步电动机(AsynchronousMachine)、励磁系统(ExcitationSystem)、水轮电机及其监测系统(HydraulicTurbineandGovernor(HTG))、永磁同步电机(PermanentMagnetSynchronousMachine)、简化的同步电机(SimplifiedSynchronousMachine)、同步电机(SynchronousMachine)。
这些模型可以仿真电力系统中发电机设备,电力拖动设备等。
接线设备模型
测量设备模型
这一部分包括接地设备、输电线母线等和用来采集线路的电压或电流值的电压表和电流表。
这些模型可以仿真电力系统的接线设备和测量设备
2.2Simulink/PowerSystem系统的仿真
在Simulink环境下,编辑模型的一般过程:
打开一个空白的编辑窗口,然后将模块库中需要的元件拉入编辑窗口,并按需要改变模块参数,再将各个模块连接起来,之后便能对整个模块进行仿真。
仿真开始最主要的就是进行设置仿真参数和选择解法器。
设置仿真参数和选择解法器,选择Simulation菜单下的Parameters命令,就会弹出一个仿真参数对话框,它主要用三个页面来管理仿真的参数,即Solver页,WorkspaceI/O页,Diagnostics页。
MATLAB中在Solver页提供了两种不同的仿真步长的选取,即变步长和固定步长。
变步长模式可以在仿真的过程中改变步长,提供误差控制和过零检测。
固定步长模式在仿真过程中提供固定的步长,不提供误差控制和过零检测。
用户还可以在第二个下拉选项框中选择对应模式下仿真所采用的算法。
表2-2中采用列表格的形式说明各个变步长的解法器。
表2-2MATLAB变步长解算器(Variable-stepSolver)
变长解法器
主要相应的数学方法
求解问题
精度阶
应用场合说明
ode45
Dormand-Prince法
非刚性
中
大部分问题采用这种算法,
ode23
Bogacki-Shampine法
非刚性
低
精度较低场合,但效率高
ode113
AdamsPECE法
非刚性
变阶:
低-高
要求速度快速的大部分场合
ode15s
NDF法
刚性
变阶:
低-中
当模型为刚性无法用ode45求解时或求解效率很低时,采用ode15s。
ode23s
二阶改进Rosenbrock法
刚性
低
当对精度要求不高时,可解决某些ode15s无法解决的刚性问题。
ode23t
自由插值的梯形法
适度刚性
低
当模型为适度刚性问题时采用ode23t。
ode23tb
隐式Runge-Kutta法
刚性
低
与ode23s类似,当对精度要求不高时,可用此法代替ode15s。
2.3小波分析的基本理论
2.3.1傅里叶分析到小波分析
傅里叶分析,包含傅里叶级数和傅里叶变换。
傅里叶变换是传统的信号分析的基础,它一直被认为是信号处理领域的完美分析方法。
它变换对应的基本函数是确定的,主要是三角函数。
傅里叶变换在频域的定位是准确的,具有最高的频域的分辨率,它实质是将信号分解到了{}的正交空间。
然而它是全局的分析方法,仅能在纯时域或完全频域,无法表述信号的时域局部性质,在时域没有任何定位性,也就是毫无分辨能力。
简单的说,傅里叶变换所反映的是整个信号全部时间下的整体频域特征,而不能提供局部时间段上的频率信息。
具体而言,即对一个信号进行Fourier变换的结果,从频谱图上可以看到这个信号所包含有哪些频率分量的信号,从其频谱图上却无法了解到这些不同频率的信号分量在时间上的关系。
为了解决这个时频无法局部化的难题,Gabor先提出了高斯变换,即短时傅里叶变换,引入了一个高斯时间窗,使得信号得以实现时频的局部化。
可惜,美中不足的是,这个时间窗是固定不变的,即高斯函数确定后,窗口的大小就唯一确定,在分析高频时,便无法得到较高的时间分辨率。
低频时,就无法得到较高的频率分辨率,而这些对于分析非平稳信号而言,都是致命的缺陷。
小波分析属于时频分析的一种,它被赋予信号分析中的显微镜的美称。
小波分解对应的小波基函数的选择是不固定的,可以是不唯一,而且对相同信号,通过不同小波基分解的结果相差会很大。
但是小波对信号所加的窗口形状是可以改变的,适合非平稳信号的分析与检测。
因此,从这几个方面说来,小波分析在分析局部特征信息时,代替傅里叶分析是一种必然的趋势。
2.3.2小波概念
小波分析方法是一种窗口大小固定但形状可变,时间窗和频率窗都可以改变的时频局部化分析方法。
1)小波的定义:
如果函数ψ(t)满足容许性条件,即
(3-1-1)
式中:
。
或(3-1-2)
将母小波经过伸缩和平移后,就可以得到一个小波序列。
2)连续小波变换(CWT)
基本小波经过伸缩和平移后得:
(3-1-3)
称其为一个小波序列,亦称为连续小波或者分析小波。
其中为伸缩因子,为平移因子。
连续小波具有以下重要性质:
线性性:
一个多分量信号