华 中 科 技 大 学 博 士 学 位 论 文
摘 要
有源电力滤波技术是一项涉及电力电子、检测技术、控制理论、信号处理及计算机技术等多学科的研究领域。随着工业现代化发展,电力系统中非线性负载的大量增加,日趋严重的谐波污染对电能的生产、利用以及设备的使用造成了很大的危害。有源电力滤波器是目前用来治理谐波污染最具前景的方法之一,它包含了多门学科的理论与技术,成为研究与开发的热点。目前有源电力滤波的研究主要集中在以下几个方面:高精度实时谐波检测和信号处理、电流快速跟踪智能控制、高功率因数变流技术以及主电路拓扑结构和非线性参数设计等。
瞬时无功功率理论是谐波检测中使用最广泛的方法,本文详细分析了基于瞬时无功功率理论的谐波实时检测性能,针对其在动态响应性能上存在的不足,研究一种基于超前网络补偿的谐波检测方法。理论分析和仿真实验结果都表明,采用校正检测系统频率响应零点的方法能有效改善瞬时无功功率理论谐波检测的动态响应性能。
基于瞬时无功功率理论谐波检测仅限于三相电路,而且受电网电压影响,小波变换在信号处理方面具有特殊的应用价值。本文对小波变换用于谐波电流实时检测的原理、小波基的选择、检测精度和动态性能进行了深入的分析,指出小波变换实时算法中的边界扩展是影响其检测性能的主要原因。通过分析电流幅值和小波系数之间的关系,提出一种基于小波系数的谐波检测实时补偿方法,使检测性能得到明显提高。
建立有源电力滤波器的数学模型,是实现其控制的理论基础。本文在导出数学模型基础上,推导出开关函数和交流输出电压的数学表达式及频谱特征。分别从满足不同性能指标的角度对有源电力滤波器主电路各项参数进行了分析,给出了它们的取值范围和设计原则。
针对谐波电流具有时变和变化率高的特点,对谐波电流跟踪预测控制进行了研究。从预测控制的三个基本特征出发,对有源电力滤波器电流预测控制的结构、原理、控制规律等方面进行研究,分析了电流预测控制参数对系统稳定性和动态响应性能的影响。采用电流预测控制有效的提高了电流控制精度和动态性能,对预测控
I
华 中 科 技 大 学 博 士 学 位 论 文
制在快速响应系统中的应用作了有益尝试。
本文以TMS320F2812控制器为核心,建立了有源电力滤波器实验平台,在该平台上对谐波检测和电压、电流控制进行了实验。对实验结果进行分析,给出了实验结论。
关键词:有源电力滤波器; 谐波; 谐波检测; 瞬时无功功率理论; 小波变换;
电流预测控制
II
华 中 科 技 大 学 博 士 学 位 论 文
Abstract
The technology of active power filter (APF) is a research field concerning power electronic, detecting technique, control theory, signal processing and computer science. With the increasing of nonlinear loads in power system, serious harmonic current pollution has threathen the running of power system and affected the use of electrical equipments. APF is one of the most prospect methods to suppress harmonic contamination. It contains several theories and technologies, becoming the research hotspot. Today, research on APF mainly concentrates on real time high precision harmonic detectiong and signal processing, current controlling, high power factor conversion technique, main circuit topology structure and nonlinear parameters design.
The instantaneous reactive power theory is widely used in harmonic detecting, and its capability is analyzed in this paper. To solve the shortcoming of instantaneous reactive power theory in dynamic performance, a new harmonic detecting method based on lead network compensation is proposed. The method improves the capability of harmonic detecting by correcting its zeros in frequency response. The theoretical analysis and computer simulations show that the method can improve the dynamic performance.
Instantanous reactive power theory based harmonic detecting is limited in three-phase circuit and is immune to the electric voltage. Wavelet transform, for its special application value in signal processing, is used in real time harmonic detecting. The principle, wavelet basis, precision and dynamic response of harmonic realtime detecting based on wavelet transform are analyzed. This method has low performance because of boundary extension in wavelet arithmetic. Through analyzing the relationship between amplitude of current and wavelet coefficients, a compensation method is proposed. This compensation method can markedly improve the performance of real time harmonic detecting based on wavelet transform.
The accurate mathematics model of active power filter is the theory base for control realization. After establishing the mathematics model, the spectrum of switching function and ac side voltage are analyzed. Studies on main circuit of APF, the calculating methods for main parameters are deduced to meet the performance.
Real time and accurate harmonic current control is the key of APF. A harmonic current control, based on predictive control, is presented in this paper. The principle, structure, character and control rule of APF with current predictive control are analyzed.
III
华 中 科 技 大 学 博 士 学 位 论 文
Especially, the stability and dynamic response of control system are discussed and the effects of parameters for performance are analyzed. Predictive control improves the precision and dynamic performance of current control and a useful attempt is made for predictive control to be applied in fast respone system.
The APF experimental platform based on TMS320F2812 is presented. Harmonic detecting and current controlling experiments are done on the platform. The results of experiment are analyzed and the conclusions are drawn.
Keywords: active power filter; harmonic; harmonic detecting; instantaneous reactive
power theory; wavelet transform; current predictive control
IV
独创性声明
本人声明所呈交的学位论文是我个人在导师指导下进行的研究工作及取得的研究成果。尽我所知,除文中已经标明引用的内容外,本论文不包含任何其他个人或集体已经发表或撰写过的研究成果。对本文的研究做出贡献的个人和集体,均已在文中以明确方式标明。本人完全意识到本声明的法律结果由本人承担。
学位论文作者签名:
日期: 年 月 日
学位论文版权使用授权书
本学位论文作者完全了解学校有关保留、使用学位论文的规定,即:学校有权保留并向国家有关部门或机构送交论文的复印件和电子版,允许论文被查阅和借阅。本人授权华中科技大学可以将本学位论文的全部或部分内容编入有关数据库进行检索,可以采用影印、缩印或扫描等复制手段保存和汇编本学位论文。 本论文属于
保密□ ,在_____年解密后适用本授权书。 不保密□。
学位论文作者签名: 指导教师签名: 日期: 年 月 日 日期: 年 月 日
华 中 科 技 大 学 博 士 学 位 论 文
1 绪论
1.1 引言
本章首先介绍电能质量及谐波产生的原因和引起的危害,然后从拓扑结构、谐波检测、控制策略等方面介绍了有源电力滤波器的基本原理及目前发展情况,最后对本文所做的研究内容和取得研究成果进行了介绍。
1.2 谐波污染及治理
1.2.1 电能质量及谐波定义
电力电子装置在电力系统、工业、交通以及家庭中的广泛应用,使得电力系统中非线性负荷大量增加。一些非线性和时变性电子装置如整流器、逆变器和各种开关电源等大规模地应用,以及工业用电弧炉、电力机车等大功率非线性和时变性电子装置的使用,使得电网中电流波形严重畸变,含有大量的谐波和无功分量,给电能的生产和传输带来危害,也影响电子设备的应用。谐波所造成的危害日趋严重,同时也造成电能极大的损耗[1][2]。在全球能源日趋紧张和环境污染空前加剧的情况下,对电力系统中谐波污染的治理和电能的合理利用,提高电能质量,成为电力电子技术领域十分紧迫和必须面对的问题[3][4][5]。近年来,由于电力电子技术发展,以及计算机技术和现代信号处理的新成果在电力电子领域的应用,极大的促进了无功功率补偿和谐波抑制的研究和应用,并且取得了一些突破性的进展。
面对谐波所带来的危害,世界各国和专家都对谐波问题予以充分的关注,特别是近几十年来,许多国家都制定了关于谐波的标准和规定,一些国际学术组织如国际电工组织(IEC)和国际电气电子工程师协会(IEEE)等经常召开国际学术会议,探讨和研究有关谐波的问题。我国先后于1984年和1993年制定《电力系统谐波管理暂行规定》(SD126-84)[6]和《电能质量-公用电网谐波》(GB/T14549-93)[7]限制电力系统中谐波的规定和国家标准。
文献[7]中对谐波测量点(harmonic measurement points)、基波(分量)(fundamental component)和谐波(分量)(harmonic component)等都分别进行了定义:
谐波(分量):对周期性交流量进行傅立叶级数分解,得到频率为基波频率大于
1
华 中 科 技 大 学 博 士 学 位 论 文
1整数倍的分量。
谐波含有率(harmonic ratio, HR):周期性交流量中含有的第n次谐波分量的方均根值与基波分量的方均根值之比(用百分数表示)。
总谐波畸变率(total harmonic distortion, THD):周期性交流量中的谐波含量的方均根值与其基波分量的方均根值之比(用百分数表示)。 1.2.2 谐波产生的原因和危害
在我国供电系统中,电网供给电压通常是频率为50Hz的正弦交流电压,当它给线性负载(如电阻、电感和电容)供电时,电流波形不发生变化,仍然是交流电。但是当正弦交流电压给非线性负载(如二极管整流器、交流调压)供电时,电流波形就会发生畸变,除了含有基波分量外,还含有谐波分量。
产生谐波的负载装置同时也大多是消耗基波无功功率的装置,谐波和无功功率既是非线性负载正常工作的要求,但也给电力系统和电力用户带来了危害。它们的影响主要表现在以下几个方面[1][8]:
(1)谐波大大增加了电力系统中发生谐振的可能性。为了提高负载的功率因数,补偿负载的无功功率,电力网中广泛使用补偿功率因数的电容器,电网中还有分布电容,同时电力系统中还存在感性部分。电力系统中的电容和感性部分组合,可能存在串联或并联的谐振条件,对一定频率的谐波造成很高的过电流或过电压而烧毁设备,引发事故;
(2)影响电气设备(如旋转电机、电容器、变压器等)的运行,使设备工作不正常,加速绝缘老化,从而缩短设备使用寿命。由于频率升高而产生的集肤效应、磁滞及涡流等,谐波电压或电流在电机绕组、转子回路及定子与转子的铁心中产生附加损耗,导致整个电机的温度升高和局部过热,缩短电机的使用寿命,情况严重时甚至会损坏电机。对变压器的影响主要体现在谐波电压可以使变压器的磁滞及涡流损耗增加,使绝缘材料承受的电气应力增大,谐波电流使变压器的涡耗增加,导致变压器局部过热,噪声增大;
(3)增加附加损耗,降低发电、输电及用电设备的效率和设备利用率; (4)使继电保护、自动装置、计算机系统以及许多电气设备出现错误动作,设备不能正常动作或操作;
(5)使测量和计量仪器、仪表不能正确指示或计量;
2
华 中 科 技 大 学 博 士 学 位 论 文
(6)干扰通信系统,降低信号的传输质量,破坏信号的正常传递,甚至损坏通信设备。
1.2.2 谐波抑制和无功补偿方法
从治理谐波污染的方式看可以分为两种。一种是从改造电力电子装置入手,使其不产生谐波,如PWM整流器[9][10][11]和多重化技术[12];另外一种方式是对电网中的谐波进行集中补偿,使电源只输出有功电流,LC无源滤波器(Passive Power Filter, PPF)[13][14][15]和有源电力滤波器(Active Power Filter, APF)[16][17][18]都属于这种方式。有源电力滤波器由于其具有可以实现动态补偿、受电网阻抗影响不大、能跟踪电网频率变化等众多优点,成为电力电子技术领域的研究热点。有源电力滤波器融合了电力电子技术、检测技术、控制技术等多种技术,特别是近几年来,许多新的理论和智能控制技术应用于有源电力滤波器中,使得有源电力滤波器成为多个学科的交叉,取得了许多成果,并且有力的推动了有源电力滤波器的发展[19][20][21]。
对产生谐波的装置进行改造,使之不产生谐波,可以从根本上抑制谐波的产生。目前该类方法主要有增加整流器相数、采用多重化连接的方式[22],或者采用PWM整流器,通过先进的控制技术、特定谐波消去法等,减少谐波的产生。这类方式可以从根本上减少进入电网中的谐波,但是并不能完全消除谐波,在装置工作时还是会产生部分谐波进入电网中,同时这种方式也增加了设备成本,大量的现有设备也无法完全进行改造。
在电网中装设滤波装置,滤除电网中的谐波,补偿电网中的无功功率。这类方法中主要有LC无源滤波器、静止无功补偿装置(Static Var Compensator, SVC)和有源电力滤波器等。下面对这类补偿装置进行简单介绍。
(1)LC无源滤波器。传统的谐波补偿装置是LC无源滤波器,它是由电容器、电抗器和电阻器适当组合而成的滤波装置,与谐波源并联,除了可以对谐波进行滤波外,还可以进行无功补偿。LC无源滤波器具有结构简单、成本低、技术成熟、维护方便、装设容易等优点,因此得到了广泛的应用,至今也是应用最多的谐波补偿方法。但是,它存在只能补偿固定次谐波、补偿特性受电网参数和运行状态的影响、可能和系统局部发生并联或串联谐振等缺点。由于无源滤波器的上述缺点,使它很难跟上现代电力系统的发展,满足电力系统的要求。
(2)静止无功补偿装置。静止无功补偿装置现在广泛应用于电力系统,它主要
3
华 中 科 技 大 学 博 士 学 位 论 文
用来对无功功率进行动态补偿,提高功率因数。它包括晶闸管控制电抗器(Thyristor Controlled Reactor, TCR)、晶闸管投切电容器(Thyristor Switched Capacitor, TSC)以及这两者的混合装置(TCR+TSC)、晶闸管控制电抗器与固定电容器(Fixed Capacitor, FC)等。80年代以来,出现了比SVC更加先进的静止无功发生器(Static Var Generator, SVG),和SVC需要大容量的电抗器、电容器等储能元件不同,SVG只在直流侧需要小容量的电容器维持电压,通过控制,SVG既可以发出无功功率,呈容性,也可以吸收无功功率,呈感性。静止无功补偿装置主要补偿电力系统中的无功功率,功能较单一。
(3)有源电力滤波器。有源电力滤波器是一种动态抑制谐波和补偿无功功率的电力电子装置,它的主电路主要由三相六桥臂的变流器组成,按照接入电网的方式可以分为并联型和串联型。与无源滤波器相比较,有源电力滤波器的优点主要表现在:它能对频率和大小都变化的谐波以及变化的无功功率进行补偿,对补偿对象的变化有极快的响应;可以同时对谐波和无功功率进行补偿,且补偿的大小可以做到连续调节;补偿无功功率和谐波所需的贮能元件容量不大;即使补偿对象电流过大,APF也不会发生过载,并能正常发挥补偿作用;APF受电网阻抗的影响不大,不容易和电网发生谐振;能够跟踪电网频率的变化,补偿性能不受电网频率变化的波动的影响;APF既可以对一个谐波和无功源单独补偿,也可以对多个谐波和无功源进行集中补偿等。正是APF的这些特点,使它成为补偿电力系统中谐波和无功的有效手段。但是有源电力滤波器也存在着价格昂贵、运行成本高和难以构造大容量补偿装置的缺点[23][24][25]。
1.3 有源电力滤波器的发展及结构类型
1.3.1 有源电力滤波器的发展
1969年,B.M.Bird和J.F.Marsh发表的论文[26]中描述了通过向交流电网中注入三次谐波电流来减少电源电流中的谐波成分,从而改善电源电流波形的新方法。尽管它没有描述一个完整的有源电力滤波器,但这种方法可以认为是有源电力滤波器基本思想的萌芽。1971年,H.Sasaki和T.Machida发表的论文[27]中首次完整的描述了有源电力滤波器的基本原理,它采用线性放大的方法产生一个与从被补偿电流中检测出来的谐波大小相等、方向相反的补偿电流,从而使电源电流只含基波成分。但
4
华 中 科 技 大 学 博 士 学 位 论 文
由于当时电力电子器件和技术发展条件的限制,该方法功率损耗大、成本高,因此仅在实验室中研究,未能在工业中实用。随着电力电子器件的发展,1976年L.Gyugyi等人提出了采用PWM控制变流器构成有源电力滤波器,确立了有源电力滤波器主电路的基本拓扑结构和控制方法[28],形成了有源电力滤波器的基本概念。在这段时期,由于电力电子器件和控制技术的限制,对有源电力滤波器的研究主要局限在实验阶段,没有获得实际应用。
1983年,H.Akagi等人提出了三相瞬时无功功率理论(Instantaneous Reactive Power Theory,又称pq理论)[29][30],并且研制出了7KVA的有源电力滤波器实验装置,用于补偿20KVA三相全控整流器产生的谐波和无功功率,获得成功。三相瞬时无功功率理论极大的促进了有源电力滤波技术的发展,并被认为是有源电力滤波器的主要理论基础之一。进入八十年代后,各种新型电力半导体器件的相继问世,电力电子器件容量和性能的不断提高,以及PWM控制技术的日趋成熟,有源电力滤波器取得了长足的进展,并在实际应用中获得了成功的应用。
九十年代后,由于计算机技术、智能控制技术和现代信号处理技术的发展,有源电力滤波器成为各种新技术的有机融合,并且极大的促进了有源电力滤波器性能的提高,每年都有大量的研究论文和成果,使得有源电力滤波器的研究不断深入。
目前,国际上有源电力滤波器的研究和应用主要以日本和欧美为代表,已经进入工业实用阶段,对改善电力系统的供电质量起到了重要的作用。在我国有源电力滤波器的研究起步较晚,并且由于各种条件的限制,目前主要以科研机构和大专院校理论实验研究为主[31],但也有部分在工业中得到了成功的应用[32]。目前,在我国电力系统谐波的治理主要以无源滤波器为主,随着能源的日趋紧张和电力污染的日趋严重,有源电力滤波器在我国有着良好的市场前景和经济效益。 1.3.2 有源电力滤波器的结构
按照有源电力滤波器接入电网的方式,可以分为并联有源电力滤波器(Shunt Active Power Filter, SAPF)和串联有源电力滤波器(Series Active Power Filter),每一类又可以单独使用或与LC滤波器混合使用等方式。
(1) 单独使用的并联有源电力滤波器
并联型有源电力滤波器主电路与谐波源并联,具有安装、维护容易的优点。单独使用的并联APF如图1-1所示,变流器与电感和直流侧电容共同组成APF主电路,
5
华 中 科 技 大 学 博 士 学 位 论 文
电感与电网相连接。APF与电源及谐波源并联,可以看成一个电流源,向电网中注入与负载中的谐波电流大小相等、方向相反的电流,使电源只输出正弦电流。图中的高频滤波器HPF与APF并联,主要用来滤除APF产生的开关频率附近的高次谐波电流[1][33][34]。
在这种结构中,APF的主电路与负载并联接入电网,其补偿的谐波电流基本上由APF提供,它可以实现对谐波、无功功率和供电点电压波动的补偿。但是,由于主电路直接连接在电网中,交流电源的电压直接加在变流器上,并且补偿电流基本上由变流器提供,因此要求变流器具有较大的容量,这一点限制了单独使用的APF在大容量场合的应用。
e非线性负载HPFAPF
图1-1 并联型有源电力滤波器
(2) 并联混合有源电力滤波器(Shunt Hybrid Active Power Filter)
与LC无源滤波器混合使用的并联型APF如图1-2所示,它主要是利用LC无源滤波器来承担部分谐波补偿任务,由APF和LC滤波器共同对谐波进行补偿,因而可以有效的减低APF的容量,充分利用APF补偿性能好和LC滤波器结构简单、成本低的优点。根据APF和LC滤波器连接方式的不同又分为图1-2中(a)(b)两种。
(a)中APF与LC滤波器并联,LC滤波器包括多组单调谐滤波器和高通滤波器,承担大部分补偿谐波和无功电流,APF主要改善整个系统的性能[35][36]。这种结构中,由于APF只补偿部分谐波,因此可以有效减少APF容量。(b)中APF与LC滤波器串联,谐波和无功电流主要由LC滤波器补偿,APF的作用是改善LC滤波器滤波特性,克服LC滤波器容易受电网阻抗的影响、易与电网阻抗发生谐振的缺点。这种结构中,APF不承受交流电源的基波电压,因此容量较小。
6
华 中 科 技 大 学 博 士 学 位 论 文
ee非线性负载APFLC滤波器非线性负载LC滤波器APF(a) APF与LC滤波器并联
(b) APF与LC滤波器串联
图1-2 并联混合有源电力滤波器
与单独使用的并联有源电力滤波器相比,混合型有源电力滤波器可以大大的降低APF的容量[37][38]。
(3) 单独使用的串联有源电力滤波器
图1-3示出了单独使用的串联型有源电力滤波器,串联型APF通过变压器与电源及谐波源串联。串联型APF主要用于补偿谐波电压,使被补偿点电压成正弦波。单独使用的串联型APF流过的是负载电流,因此,对主电路容量要求较高,同时与并联型APF相比,它的安装、投切、保护等也较复杂。
e非线性负载HPFAPF
图1-3 串联有源电力滤波器
(4) 串联混合有源电力滤波器(Series Hybrid Active Power Filter)
混合型串联有源电力滤波器[39][40]结构如图1-4所示。混合型串联有源电力滤波器中,谐波和无功主要由LC滤波器补偿,APF主要用来改善LC滤波器的滤波特性,混合型串联有源电力滤波器同样可以大大的降低APF的容量[41][42]。
7
华 中 科 技 大 学 博 士 学 位 论 文
图1-4 串联混合有源电力滤波器
1.3.3 有源电力滤波器的原理
图1-5是并联有源电力滤波器的基本原理图。
eisilicRC谐波检测i*cic+-非线性负载L驱动APF电流控制
图1-5 三相并联有源电力滤波器原理图
图1-5中,三相电压源e给非线性负载供电,当不用APF对谐波进行补偿时,电源电流is与负载电流il相等,is中除了含有与电压源同频率的基波电流外,还含有大量的谐波电流,即
is=il=if+ih (1.1) 其中,if为基波电流,ih为谐波电流。
8
华 中 科 技 大 学 博 士 学 位 论 文
采用APF对谐波进行补偿的目的是控制APF,向电网中输入与负载谐波ih大小相等的电流,最终使三相电源只输出基波电流,即在图1-5所示电流方向下有
ic=ih (1.2) is=if (1.3) 为了控制APF输出补偿电网中的谐波电流,首先必须检测出三相电网中的谐波
*
作为APF电流控制的输入给定,通过电流控制,使电流,将检测出来的谐波电流ic
*
APF输出电流跟踪检测出来的谐波电流。如果检测出来的ic中除了包含谐波电流外
还含有无功电流,则通过电流控制,APF输出的电流ic中也包含谐波和无功电流,这样,APF可以同时补偿谐波电流和无功电流[1][32][43]。
1.4 有源电力滤波器的主要技术
APF是对实时性要求非常高的电力电子装置,电网中谐波产生后,从检测谐波到电流跟踪控制,最后APF输出补偿电流,不可避免的有时间延时。如何减少各个环节的延时,提高个环节的精度,最终提高有源电力滤波器的补偿精度和动态响应性能,是有源电力滤波器研究中需要解决的问题[44][45][46]。涉及到有源电力滤波器性能的关键技术主要有以下几个方面。
1.4.1 高精度实时谐波电流检测
(1)模拟滤波器。它是最早采用检测谐波电流的方法。模拟滤波器有两种,一种是滤除基波电流来得到谐波电流分量;另外一种是将滤波器设计成带通滤波器得到基波分量,再与被检测电流相减得到谐波电流分量。模拟滤波器原理和结构简单,成本低,能够滤出一些固有频率的谐波,但缺点是很难设计截止性能好的模拟滤波器,检测精度差,误差较大,同时对电力系统及电路元件参数十分敏感。随着计算机技术、数字控制技术的发展,尤其是DSP的迅速发展,模拟滤波器已经很少采用。
(2)瞬时无功功率理论。日本学者H.Akagi1983年提出瞬时无功功率理论[29],以瞬时实功率P和瞬时虚功率Q的定义为基础,即PQ理论,后又发展为以瞬时有功电流ip和瞬时无功电流iq为基础的理论体系[47][48]。PQ检测方法在检测过程中需要电网电压信号瞬时值,当电网电压畸变时,谐波电压参与运算,最终会导致检测结果
9
华 中 科 技 大 学 博 士 学 位 论 文
有误差。ipiq检测方式在检测过程中只取电压信号的频率和相位信息,畸变电压的谐波成分在运算过程中不出现,因而检测结果不受电压波形畸变的影响,检测结果是准确的。因此,ipiq方法适应的范围更宽,能适应不对称电网和电压波形畸变时的情况。
瞬时无功功率理论能准确的检测三相电路中的谐波电流,实时性较好。但是在瞬时无功功率理论对谐波电流进行检测时采用了数字低通滤波器,因此,检测出的结果会有一定时间的延时。瞬时无功功率理论是目前有源电力滤波器中采用最多的一种谐波检测方法。
(3)基于傅立叶变换(Fourier)的检测方法。该方法是采用傅立叶变换(FFT)对电网电流进行计算,得到电网电流中的谐波分量[49][50][51]。傅立叶变换是对一段时间的信号进行分析,得到其频率成分,它是一种纯频域的分析方法,它在频域的定位是完全准确的,而在时域无任何定位性,即不能提供任何局部时间段上的频率信息。并且该方法计算量大、因而实时性不好。对非整数次谐波的检测时,由于信号负频分量引入的长范围泄露以及窗扇形损失引入的短范围泄露,导致频谱泄漏和栅栏现象等缺点,从而检测出的谐波幅值、相角和频率有误差。
(4)基于自适应的检测方法。基于自适应的谐波电流检测方法利用了自适应信号处理中的噪声对消法:当原始输入是信号和与信号不相关的加性噪声组成时,参考噪声信号通过自适应滤波器处理后与原始输入相减,经过系统自适应调整使最小均方误差最小,实现抵消干扰噪声、检测出信号的目的。当自适应方法用于谐波检测时,将基波分量看作噪声,通过最小均方误差算法来调整权重,得到被测电流中的基波分量,从而检测出谐波。该方法具有较高的检测精度,能根据系统的变化进行调节,但缺点是当电网电流发生变化时,自适应的参数调整需要一定的时间,影响了其响应速度[52][53][54]。
(5)基于小波变换(Wavelet)的检测方法。小波变换发展了短时傅立叶变换(Short Time Fourier Transform, STFT)局部化的思想,而且它可以随着频率的变化调整窗口大小,可以对信号在时间和频域都精确定位。基于小波变换的谐波检测方法克服了傅立叶分析在频域完全局部化而在时域无局部化的缺点,它能够准确计算出一定时间的频率分布,因此可以准确检测出谐波电流。目前大多数对于小波变换谐波检测的研究[55][56][57]都是对一段采样数据进行小波分析,得到其中的谐波成分,而没有涉及到实时检测。在小波变换谐波实时检测过程中,分析当前时刻的谐波成分时,需要
10
华 中 科 技 大 学 博 士 学 位 论 文
对数据进行边界扩展,对于电网电流这种周期性比较强的信号,往往采取周期扩展的方法。如果电网电流是稳定的,周期扩展能够取得很好的检测精度,但是当电网电流发生变化时,会产生误差,检测精度会下降,导致谐波电流检测动态性能较差。
检测精度和实时性是判断谐波检测方法的重要指标,各种检测方法都有优点,但每种检测方法都存在局限性。在各种谐波电流检测方法中,神经网络、小波变换等现代数字信号处理进行谐波电流实时检测是极具潜力的新型谐波检测方法。
本文对基于小波变换的谐波电流实时检测进行了深入的研究,对造成实时检测性能降低的原因进行了分析,并对如何提高小波变换实时检测性能进行了深入的研究。
1.4.2 主电路结构及参数设计
目前,有源电力滤波器主电路主要采用PWM变流器的形式,当采用单个变流器不能满足系统容量要求时,可以采用多重化或多电平(Multi-level)的主电路结构形式[1][58][59]。
(1)单个PWM变流器的主电路结构
根据主电路直流侧储能元件的不同,可以分为电压型和电流型两种,图1-6(a)、(b)分别示出了三相电压型和电流型单个PWM变流器的主电路形式。
LVaVbVaVcC
VbVcL
图1-6 (a)三相电压型PWM变流器 (b)三相电流型PWM变流器
电压型PWM变流器直流侧接有大电容,正常工作时需要通过控制保持直流侧电压恒定,可以看作电压源,交流侧通过电感与被补偿点相接。在实际工作过程中,直流侧电压会发生波动,交流侧电感对输出电流有阻碍作用,因此,电容和电感的大小对APF装置的成本和性能有很大的影响。电压型APF直流侧电容损耗较小、体积较小,因此适合构成大容量由源滤波器。
电流型PWM变流器直流侧接有大电感,正常工作时需要通过控制保持直流侧
11
华 中 科 技 大 学 博 士 学 位 论 文
电流恒定,可以看作电流源。在实际工作过程中,直流侧电流会发生波动,因此,电感的大小对APF装置的成本和性能有较大影响。电流型APF不存在直流测短路的问题,可靠性较高,动态性能好,滤波器精度较高,但由于直流侧电感中始终有电流流过,因此损耗较大。
(2)多重化主电路结构形式
当实现大容量的APF装置时,电力电子器件随着容量的增大开关频率降低,较低的开关频率会导致APF性能的降低。有源电力滤波器多重化主电路结构是将多个PWM变流器串联或并联起来,通过变压器或各相的电感将多个输出电流或电压组合在一起,使整体容量变大,可以较容易的在不改变单个器件容量的基础上实现APF装置的大容量。采用多重化方式可以利用各变换器现有的控制技术,组合成大容量变换器,还可以有效降低单个电力电子器件的开关频率,减少器件的开关损耗,以较低的器件开关频率获得较高的等效开关频率,改善电流补偿性能。
(3)多电平主电路结构形式
多电平变流器不使用变压器和电抗器,而是通过增加电力电子器件,设计多电平主电路拓扑结构,将变流器的输出由传统的两电平输出变为多电平输出。由于不采用输出变压器和动态均压电路,开关频率低,实现了开关器件所承受的电压应力小、效率高的目标,从而提高APF容量。多电平变流器也可以以较低的器件开关频率获得较高的等效开关频率,改善电流补偿性能,但多电平主电路控制方法较复杂。
(4)参数设计
根据APF接入电网的形式和与LC无源滤波器组合的形式,有源电力滤波器有不同的拓扑结构,主电路中各种元器件的参数和LC无源滤波器直接影响有源电力滤波器的性能。由于有源电力滤波器需要产生的谐波电流是时变的,并且变化速率快,主电路中电感、高频滤波器、电容、无源滤波器等各项参数对其谐波电流的跟踪性能影响很大,因此选择合适的主电路结构和参数对有源电力滤波器的性能有重要的影响。由于APF结构多样,补偿的谐波源也多种多样,对APF的容量和谐波补偿的性能指标也有不同的要求,目前,对于APF主电路各项参数的设计没有统一的理论,主要根据被补偿的谐波源选择主电路结构形式,根据主电路结构、开关器件的开关频率、电压等级、电流等级选择合适的开关器件,根据APF装置的性能和成本设计直流侧电压、电容值和交流侧电感值以及LC无源滤波器,因此对APF主电路参数有不同的设计方法。
12
华 中 科 技 大 学 博 士 学 位 论 文
文献[60]建立了单相APF的数学传递函数,通过传递函数分析了补偿电流的大小、直流侧电压的大小以及主电路参数对系统性能的影响,并通过实验进行了验证。文献[61][62]建立了并联APF频率域模型,从频率分析的角度分析了参数对系统性能的影响。文献[63][64]从内模结构的角度分析了APF的模型,并建立了内模控制下APF的离散数学模型。文献[65]将滤波器的设计看作一个多目标优化问题,在满足约束条件下建立多目标优化模型,并选用多目标模糊优化的方法进行求解并联混合型APF的系统参数。文献 [66]提出了一种串联型APF中低通滤波器的设计方法,并利用仿真软件采用在滤波范围内扫描的方法对参数进行优化。
有源电力滤波器是一个复杂系统,只要是满足APF装置性能指标的参数设计都是合理的。建立有源电力滤波器精确的数学模型,在其指导下进行主电路各项参数设计具有重要的意义。
1.4.3 快速高精度电流跟踪控制策略
当准确的检测出电网中的谐波电流后,控制APF主电路,使APF输出电流跟踪谐波电流变化。由于谐波电流具有时变和高变化率(di/dt)的特点,这就要求APF电流控制器具有较快动态响应性能和较高的控制精度,除了具有较快的动态响应性能,电流控制器的稳定性也是必须要考虑的因素。目前对APF电流控制方法的研究比较多,主要有比例积分控制、滞环控制、单周控制、滑模控制等[67]。
(1)比例积分(Proportion integral, PI)控制
有源电力滤波器电流控制传统采用的是PI调节器,它采用比例积分对三相电流误差计算,或通过坐标变换在两相坐标中进行解耦控制,控制量与三角载波进行比较或通过矢量算法生成PWM调制波,这样组成的控制系统是基于把电流误差信号控制为最小的原则来进行设计的[68][69]。PI控制原理简单,易于模拟和数字实现,器件开关频率固定,较容易设计滤波器消除有开关频率引起的高次谐波电流。但是正如前面所提到的,由于控制器的给定谐波信号是时变的,采用PI控制不能完全消除稳态误差,并且为了保证系统的稳定性,PI控制器的阶数不能太高,放大器的增益也是有限的,因此PI控制方法的跟随误差较大,不能取得很好的控制效果[70][71]。
(2)电流滞环控制(Hysteresis Control, HC)
电流滞环控制是对电流瞬时值进行比较的一种控制方法,是目前电力电子装置中广泛应用的一种非线性闭环电流控制方法。APF滞环电流控制的基本原理是:将
13
华 中 科 技 大 学 博 士 学 位 论 文
检测到的谐波电流作为输入给定,与APF的受控输出电流进行比较,它们的误差作为滞环比较器的输入,通过滞环比较器产生控制变流器开关器件通断的PWM控制信号,使APF的输出电流跟踪给定谐波电流。滞环电流控制是基于电流暂态的控制,具有动态响应速度快、鲁棒性好的优点,并且实现简单[72][73][74]。但缺点是在设定好滞环控制器的滞环宽度后,电流跟随误差是固定的,控制器输出的控制信号频率是变化的,具有较宽的频谱,滤波较困难,给滤波器参数的设计带来了一定的难度,并且可能使开关频率过高,甚至超过器件允许的最高工作频率而导致器件损坏。为了解决开关器件频率变化的缺点,许多研究提出了各种改进措施[75][76]。文献[75][76]提出一种定频滞环电流控制方法,这种方法可以使开关器件的工作频率固定,但是在这种控制方法中补偿电流的跟随误差是不固定的。
(3)滑模控制(Sliding Mode Control, SMC)
滑模控制又称为变结构控制(variable structure),它最早由前苏联学者在20世纪50年代提出,由于滑模控制系统具有较强的鲁棒性,对被控对象的模型误差、对象参数的变化以及外部干扰有极佳的不敏感性,同时也具有响应快速、动态品质好的优点,引起了广大学者的关注。由于电力电子装置主电路具有变结构特性,近年来,滑模控制也逐渐应用到电力电子及有源电力滤波器的研究中。
有源电力滤波器滑模控制的基本思想是设计滑模函数,利用控制的不连续性,依靠高频开关将系统的闭环相轨迹引导到一个滑动面上,使控制系统在滑模面上的运动渐进稳定,动态性能良好,即存在性;同时使所有的运动轨迹于有限的时间内到达滑模面,即可达性。当系统的状态向量进入滑模面后,就被约束在限定空间做小幅度、高频率的滑模运动。此时,系统的动态品质由滑模面的参数决定,而与系统参数、外部扰动无关,因此滑模控制对被控对象的模型误差、对象参数的变化以及外部扰动具有较强的不敏感性,具有较强的鲁棒性。在滑模控制中,合理的选择滑模面对控制系统的动态性能极为重要。滑模控制的不连续开关特性容易引起系统的颤动,这种颤动可能将系统中存在的高频成分激励起来,甚至使系统不稳定,这是有源电力滤波器滑模控制中需要注意的[77][78][79]。
(4)单周控制(One Cycle Control, OCC)
单周控制理论是美国学者K.M.Smedley博士在1991年提出的,它是一种通用型非线性大信号控制方法,特别适合于开关电路的控制。
单周控制的基本思想是在每个控制周期内强迫开关变量的平均值与控制参考电
14
华 中 科 技 大 学 博 士 学 位 论 文
压相等或成一定的比例,从而消除稳态或瞬态误差[80][81]。和其它定频率的控制方法相比较,单周控制方法控制电路简单,响应速度快,控制稳定性能好,具有很好的鲁棒性等特点。单周控制在控制回路中不进行误差综合,它在一个控制周期中可以自动消除稳态和瞬态误差,因此,上一个周期的误差不会影响到下一周期。单周控制的缺点是抗干扰性能较差,容易受到外界条件的影响,即使在器件误差很小的情 况下,噪声干扰也会导致控制电路的条件发生变化,从而影响电路的正常工作[82][83]。
(5)预测控制(Predictive Control, PC)[84][85]
现代控制理论自60年代初期发展以来,对自动控制技术的发展起到了积极的推动作用,它以状态空间为基础,需要对控制对象建立精确的数学模型,对于一些确定性的、线性的、时不变和易于建立数学模型的控制对象,现代控制理论能够取得较好的控制效果。实际的工业过程中,控制对象经常具有非线性、时变性和不确定性,且大多数工业过程是多变量的,难于建立精确的数学模型。采用自适应、自校正控制技术虽然能在一定程度上解决不确定性问题,但需要在线辨识对象模型,所以算法复杂、计算量大,且抗干扰能力和自适应能力差,系统的鲁棒性问题也尚待进一步解决,这些都在一定程度上限制了现代控制理论在工业控制领域的应用范围。为了解决这些问题,从70年代以来,人们除了加强对生产过程的建模、系统辨识、自适应控制、鲁棒控制等的研究外,开始打破传统控制思想的束缚,试图面向工业过程,寻找各种对模型要求低、在线计算方便、控制综合效果好的算法。计算机的发展也为这类新算法提供了物质基础。预测控制就是在这种情况下发展起来的一类新型计算机控制算法。
预测控制自产生以来已经出现了多种不同的预测控制算法。70年代至80年代初,Richalet、Mehra等提出的建立在非参数模型响应基础上的模型预测启发控制(Model Predictive Heuristic Control, MPHC)和模型算法控制(Model Algorithmic Control, MAC),以及Culter等提出的建立在非参数模型阶跃响应基础上的动态矩阵控制(Dynamic Matrix Control, DMC)等。这类算法不需要建立精确的控制对象数学模型,采用易于从生产现场测得的脉冲响应或阶跃响应来描述控制对象的动态行为,并不需要过程模型的结构和参数有关先验知识,也不需要通过复杂的系统辨识来建立过程的数学模型,它根据某一优化指标来设计控制系统,在控制过程中采用不断的在线滚动优化,且在优化过程中进行反馈校正,在一定程度上克服了系统的不确定性,使系统的鲁棒性得到增强。80年代后期,出现了基于离散参数模型的预测控制算法,
15
华 中 科 技 大 学 博 士 学 位 论 文
主要有Clarke的广义预测控制(Generalized Predictive Control, GPC)、Lelic提出的广义预测极点配置控制(Generalized Poles Placements Control, GPP)等。与基于非参数模型的MAC和DMC等预测控制算法相比较,基于参数模型的预测控制算法仍然保留了它们的基本特征,但是在预测模型上采用的是具有一定结构和参数的离散受控自回归积分滑动平均模型(CARIMA)或受控自回归滑动平均模型(CARMA)。由于参数模型是最小化模型,需要已知模型结构,但需要知道模型参数远较非参数模型少,减少了预测控制算法的计算量。
针对于有源电力滤波器这种响应速度快、非线性强的电力电子装置,本文采用预测控制对其电流跟踪控制进行了研究,对电流控制系统结构、稳定性、动态响应性能以及参数与系统性能的关系都进行了分析和探讨。
此外,直流侧电压控制也是APF关键技术之一。对于图1-6(a)所示的电压型有源电力滤波器,它的直流侧必须有一个稳定的直流电压,APF才能正常工作,目前,电压型有源电力滤波器直流侧通常是并接一个电容器,来维持直流侧电压。直流侧电压的恒定也是APF获得较好的性能指标的重要因素,因此,如何控制直流侧电压的恒定也是APF控制中的一个重要方面。目前比较常用的方式是将电压环的控制加入到有功电流环节的控制中,通过电流的控制来控制直流侧电压的恒定,所以在APF系统中有电流和电压两个闭环控制环节。
1.5 本文所做的工作
针对于上述有源电力滤波器的主要技术,本文以校企合作项目“有源电力滤波器(APF)的研究与开发”为背景,以提高有源电力滤波器性能指标为目的,对有源电力滤波器的一些关键技术进行了研究和探索,分以下几个方面。
(1)在查阅国内外文献、资料的基础上,阐述了有源电力滤波器的研究状况,对各种技术进行分析和总结。
(2)详细分析了基于瞬时无功功率理论的谐波检测的精度和动态性能,针对瞬时无功功率理论对谐波检测动态性能不好的缺点,提出一种基于超前网络补偿的谐波电流检测方法。
(3)将现代信号处理的发展成果--小波变换应用于谐波实时检测,对基于小波变换的谐波实时检测性能进行了分析,指出边界扩展是影响其实时检测精度和动态
16
华 中 科 技 大 学 博 士 学 位 论 文
响应性能的主要原因。针对小波变换实时信号检测的局限性,在理论分析和数学推导的基础上,提出一种实时补偿方法,有效提高了小波变换实时信号检测的动态响应性能。
(4)建立了并联有源电力滤波器的数学模型,在此基础上对开关函数和交流侧输出电压进行理论分析,得到了开关函数和网侧电压的数学表达式,对有源电力滤波器主电路各项参数的设计和电流控制提供了理论基础。
(5)对有源电力滤波器系统的功率流动进行了分析,在基于允许最大电压波动下对直流侧电容、满足给定误差和动态性能指标的交流侧电感和高频滤波器电容进行设计,得出了它们的取值范围。
(6)针对有源滤波器这种电流信号具有时变、高变化率的特点和滞后对系统动、静态响应有严重的负面影响的系统,提出采用预测控制解决谐波电流跟踪的控制方案,对电流预测控制(Current Predictive Control, CPC)的预测模型、滚动优化和反馈校正等方面进行研究,并对电流预测控制下滤波系统的稳定性和动态响应性能进行分析,研究预测控制参数对系统稳定性和动态响应性能的影响,给出了它们的取值方法和范围。
(7)在对有源电力滤波器理论研究的基础上,研制了以TMS320F2812为控制核心的有源电力滤波器实验装置,并进行了实验研究和分析,对本文提出的理论分析和设计方法进行了实验验证。
17
华 中 科 技 大 学 博 士 学 位 论 文 2 瞬时无功功率理论基础和实时谐波检测方法的改进
2.1 引言
本章介绍基于瞬时无功功率理论的谐波检测方法,分析其检测精度和动态响应性能,分析影响检测性能的因素,提出一种基于超前网络补偿的谐波检测方法。
2.2 瞬时无功功率理论基础
早期的瞬时无功功率理论是以功率的定义为基础的pq理论,后来发展到以瞬时有功电流ip和瞬时无功电流iq为基础的理论体系。与p、q检测方式相比较,ip、iq检测方式不需要电压信号的瞬时值,因此检测结果不受电压波形畸变的影响,同时检测结果也不受电网电压不对称的影响,因而应用更加广泛[1][86]。
在三相电路中,设三相电流的瞬时值为ia、ib、ic,它们由基波和谐波成分组成,表达式为
∞
⎧
ia=2I1sin(ωt+ϕ1)+∑2Insin(nωt+ϕn)⎪
n=2
⎪∞
2π2π⎪
=−++−iIωtϕInωt2sin()2sin[()+ϕn] (2.1) ⎨b∑n11
33n=2⎪∞
4π⎪i=2Isin(ωt−4π+ϕ)+−2Isin[n(ωt)+ϕn]∑cn11⎪33n=2⎩
其中ω为电网基波电流角频率,I1和ϕ1分别为基波电流的有效值和相位,In和
ϕn分别为各次谐波电流的有效值和初相角,n=3k±1,k=1,2,3...。
将电网电流变换到两相正交静止的α-β坐标系,可以得到α、β坐标系下的瞬时电流iα、iβ
⎡ia⎤
⎡iα⎤⎥ (2.2) i=⎢⎥=C32⎢ib⎢⎥⎣iβ⎦⎢⎣ic⎥⎦
其中C32=
2⎡1
3⎢⎣0−12−12⎤
。
32−32⎥⎦
18
华 中 科 技 大 学 博 士 学 位 论 文
在α-β坐标系下将电流矢量分别向电压矢量及其法线上进行投影,就可以得到有功电流ip和无功电流分量iq。
⎡ip⎤⎡iα⎤
⎢⎥=C⎢⎥ (2.3)
⎣iβ⎦⎣iq⎦
⎡cosωt
其中C=⎢
⎣sinωtsinωt⎤
,图2-1为有功电流、无功电流分解示意图。 ⎥−cosωt⎦
βipeiβ0iqiiαα
图2-1 有功电流、无功电流分解示意图
将电流表达式(2.1)代入式(2.2)和式(2.3)可以计算出瞬时有功电流和无功电流
∞
⎧
ip=3I1sinϕ1±3∑Insin[(n±1)ωt+ϕn]⎪⎪n=2
(2.4) ⎨∞
⎪iq=3I1cosϕ1−3∑Incos[(n∓1)ωt+ϕn]⎪n=2⎩
分析式(2.4),ip和iq都是由直流分量和交流分量组成,直流分量为
⎧⎪i=3I1sinϕ1
⎨p (2.5)
⎪⎩iq=3I1cosϕ1直流分量ip、iq是由基波电流产生的,如果将ip、iq从ip、iq中分离出来,经过合成就可以得到α-β坐标系中的基波电流iαf、iβf
⎡ip⎤⎡3I1sin(ωt+ϕ1)⎤⎡iαf⎤
⎢⎥=C⎢⎥=⎢⎥ (2.6)
ii⎣βf⎦⎣q⎦⎣−3I1cos(ωt+ϕ1)⎦再经过两相静止坐标到三相坐标的变换可以得到三相坐标系中的基波电流iaf、
19
华 中 科 技 大 学 博 士 学 位 论 文
ibf、icf
⎡⎤⎢2I1sin(ωt+ϕ1)⎥⎡iaf⎤
⎢⎥i2⎢⎥T⎡αf⎤ ⎢ibf⎥=C32⎢⎥=⎢2I1sin(ωt−π+ϕ1)⎥ (2.7)
3⎥⎣iβf⎦⎢⎢icf⎥4⎣⎦⎢2Isin(ωt−π+ϕ)⎥
11
⎢⎥3⎣⎦它等于原电流中的基波成分,可见,ipiq检测方法准确的检测出了ia、ib、ic中的基波电流,从原电流中减去基波电流就可以得到谐波电流总和iah、ibh、ich。ipiq谐波检测的原理可以用图2-2表示[87][88]。
sinea
PLL-cosiaibiciαC32iβCipiqLPFLPF_ip_iqCiαfiβfC23iafibficf--+-iah
++ibhich
图2-2 ipiq运算方式原理图
图2-2中,通过一个锁相环(PLL)和一个正、余弦信号发生器产生与ea同相位的正弦信号sinωt和余弦信号−cosωt。
低通滤波器(LPF)的目的是从有功电流ip和无功电流iq中分离出直流分量ip和
iq。如果需要对基波无功功率也进行补偿,那么在计算出ip、iq后,只需用ip进行合成得到基波有功电流,最后可以得到原电流中的谐波和无功电流的总和。
2.3 基于瞬时无功功率理论的谐波检测性能分析
2.3.1 ipiq检测方式数学模型描述
基于瞬时无功功率理论的谐波检测是在时域中对信号进行处理,该方法对当前时刻采样点进行计算,检测谐波电流时有较好的实时性,但是在谐波检测中采用了低通滤波器,它对信号的处理会有一定时间的延迟,频率响应特性会影响检测精度。实际检测过程中采样、计算存在时延,对检测性能也有一定的影响。
20
华 中 科 技 大 学 博 士 学 位 论 文
巴特沃斯(Butterworth)滤波器[89][90]的幅度响应在通带内具有最平坦的特性,且在通带和阻带内幅值是单调变化的,它的幅频特性比较适合瞬时无功功率理论的谐波检测的要求,下面以采用二阶巴特沃斯低通滤波器的瞬时无功功率理论检测谐波电流来分析ipiq检测方式的数学描述和性能。
设图2-2ipiq检测方式中二阶巴特沃斯低通滤波器的传递函数为G0(s),它具有如下的形式
ω02
(2.8) G0(s)=2
2
s+2ξω0s+ω0
其中ω0为截止频率,ξ为阻尼系数。
被检测三相电流的拉氏变换为I(s),检测得到的谐波电流拉氏变换为Ih(s),由图2-2可以得到
Ih(s)=[E−C23CCC32G0(s)]I(s) (2.9)
⎡100⎤
⎥为三阶单位矩阵。可以得到三相检测系统传递函数为 其中E=⎢010⎢⎥
⎢⎣001⎥⎦
G(s)=E−C23CCC32G0(s) (2.10)
代入C23、C和C32各矩阵值可以求得各相电流检测的传递函数
⎧Ga(s)=Iah(s)Ia(s)=1−G0(s)⎪
⎨Gb(s)=Ibh(s)Ib(s)=1−G0(s) (2.11)
⎪G(s)=I(s)I(s)=1−G(s)
chc0⎩c
将式(2.8)代入上式可以求得ipiq谐波电流检测方式中各相电流检测的传递函数
s2+2ζω0s
,k=a,b,c (2.12) Gk(s)=2
2
s+2ζω0s+ω0
2.3.2 ipiq谐波检测静态特性分析
下面以二阶巴特沃斯低通滤波器的截止频率为30Hz为例具体分析检测性能。此时检测系统的传递函数为
21
华 中 科 技 大 学 博 士 学 位 论 文
s2+266.6s
Gk(s)=2,k=a,b,c (2.13) 4
s+266.6s+3.553*10通过式(2.13)可以画出瞬时无功功率ipiq检测方式的幅频特性、相频特性曲线,分别如图2-3、图2-4所示。
幅值/dB
图2-3 谐波电流检测幅频特性曲线
相角/Deg
图2-4 谐波电流检测相频特性曲线
分析图2-3和2-4可以得到以下结论 (1)瞬时无功功率检测系统具有高通性质;
(2)幅频特性图上对应各次谐波的幅值都是0dB,没有衰减;
(3)相频特性图上对应各次谐波都有一定的相移,最大相移为2次谐波,约20º,随着谐波频率的增高,相移逐渐减少,5次谐波有10º相移,7次谐波有7º相移,高于1kHz的谐波基本无相移。
综上所述,通过瞬时无功功率理论检测的谐波幅值上没有衰减,低次谐波有较少的相移。对于主要含有5次、7次以上谐波的三相对称电路,具有较高的检测精度。
22
华 中 科 技 大 学 博 士 学 位 论 文
2.3.3 ipiq谐波检测动态特性分析
通过式(2.13)可以画出瞬时无功功率ipiq检测方式的群延迟(Group Delay)曲线图2-5和阶跃响应曲线图2-6。
群延迟/s频率/Hz
图2-5 谐波电流检测群延迟曲线
幅值/A
图2-6 谐波电流检测阶跃响应曲线
群延迟代表波形通过滤波器之后的波形失真,在图2-2中,实际使用的是低通滤波器,图2-5显示对直流信号的延迟有7.5ms,因此,瞬时无功功率理论对谐波检测有较大的时间延迟。图2-6描述检测系统在单位阶跃函数作用下的动态过程,它的上升时间为12ms,峰值时间为22ms。从动态特性图分析可以得到以下结论
(1)检测系统有较大的群延迟,对检测波形有一定的影响; (2)阶跃响应上升时间有半个基波周期的时间,快速响应性能较差; (3)动态响应超调量为4.3%。
23
华 中 科 技 大 学 博 士 学 位 论 文
2.3.4 仿真分析
(1)负载恒定情况下的谐波电流检测
对负载恒定的电网电流检测可以验证检测精度。三相整流桥带阻感负载是比较常见的非线性负载,采用瞬时无功功率ipiq检测方式对该电网电流进行谐波检测,电路结构图如下
图2-7 谐波检测电路结构图
在saber仿真环境中建立图2-7的仿真模型,利用2.2节ipiq谐波检测方式对电网电流ila、ilb和ilc进行谐波检测。其中,负载电阻R=10Ω,负载电感L=1mh,电流采样频率fs=20khz,低通滤波器采用数字式二阶巴特沃斯低通滤波器,截止频率为
30hz,其数字离散函数表达式为
2.206*10−5Z2+4.412*10−5Z+2.206*10−5
(2.14) G0(z)=2
Z−1.987Z+0.9868图2-8(a)、(b)、(c)分别是负载恒定时a相电网电流、检测出来的基波电流和谐波电流波形,图2-9(a)、(b)、(c)分别是它们的频谱图,表2-1列出了它们的幅值及检测误差。
24
华 中 科 技 大 学 博 士 学 位 论 文
60.030.00.0-30.0-60.060.030.00.0-30.0-60.060.030.00.0-30.0-60.00.10.110.120.130.140.150.160.170.180.19时间/s
图2-8 a相电网电流、基波电流和谐波电流检测波形
(b)(c)(a)谐波电流/A基波电流/A电网电流/A60.040.020.00.060.040.020.00.060.040.020.00.001002003004005006007008009001k1.1k1.2k1.3k1.4k1.5k1.6k1.7k1.8k1.9k2k频率/Hz
图2-9 a相电网电流、基波电流和谐波电流频谱
表2-1 基波和主要谐波的幅值及检测误差
50Hz 5次
7次 11次 13次
0
ila(A) 56.851 12.543 6.716 4.879 3.997 ifa(A) 56.588 0.048 0.0364 0 误差(%)
12.535 6.782 4.824 4.041 iha(A) 0.865
1.52 -0.06 0.98 -1.12 1.10
表2-1中检测出的基波电流检测误差为-0.46%,检测出的基波电流谐波总畸变率为0.89%,对主要的几次谐波检测误差都在1.5%以内。仿真结果可以看出ipiq检测方式对稳定电网电流检测具有较高的检测精度,这与2.3.2节静态性能的分析是相符合的。
(2)负载变化情况下的谐波电流检测
25
华 中 科 技 大 学 博 士 学 位 论 文
在实际应用中,负载往往是变化的,对负载变化情况下的电网电流谐波检测可以观察检测方法的动态性能。在其它条件不变的情况下,将图2-7中负载电阻从10Ω变化到5Ω,再变化到10Ω。图2-10分别是电网电流及ipiq谐波检测方式中ip、iq的变化曲线。
(a)电网电流/A(c)iq/A(b)ip/A
图2-10 电网电流、ip、iq波形
图2-2中低通滤波器主要是对直流信号产生影响。从图2-10(b)可以看出在当前负载变化情况下,通过低通滤波器的直流信号ip经过约30ms才达到新的稳态,最终也会反映到谐波检测的动态响应。图2-11分别是负载变化情况下电网电流、检测到的基波电流和谐波电流波形
谐波电流/A(b)电网、基波电流/A(a)
图2-11 电网电流、基波电流及谐波电流检测波形
为了清楚的看出谐波检测系统的动态性能,图2-11(a)中将电网电流和检测出的基波电流放在同一坐标轴中。电网电流在0.1s和0.2s发生变化,当负载电阻从10Ω变化到5Ω时,正如2.3.3分析的,由于低通滤波器的延时,ip、iq的变化并不能立
26
华 中 科 技 大 学 博 士 学 位 论 文
即发生突变,而是有约30ms的响应时间,才达到新的稳定状态,因此检测到的基波电流和谐波电流在这段时间都不准确。通过仿真结果可以看出,ipiq检测方式的动态响应性能不高,与2.3.3节的分析是一致的。
2.4 基于超前网络补偿的瞬时无功功率谐波检测研究
2.4.1 影响检测系统动态性能的原因分析
瞬时无功功率谐波检测由于采用低通滤波器,对检测系统动态性能有较大的影响。分析图2-2可以知道,低通滤波器主要是影响ip和iq中的直流信号,导致直流滤波输出ip和iq响应较慢。当被检测电流突变时,ip和iq向新的稳态值变化,需经过一段时间后才能达到新的稳态。
仍然以式(2.13)为例分析影响检测系统动态响应性能的原因。式(2.13)中,单相检测传递函数有两个共轭极点p1=−133.3+133.3j和p1=−133.3−133.3j,有两个零点z1=0和z2=−266.6。它们的分布关系如图2-12所示。
jωp1133.3z2-266.6-133.3z10-133.3p2σ
图2-12 谐波检测系统零极点分布图
共轭极点决定了检测系统的频率响应性能,零点对系统的动态性能有很大的影响[91][92]。式(2.13)的一个零点位于原点,另外一个零点远离原点,对于一个具有高通性质的系统来说,它的零点分布不合理,这是造成系统动态响应性能较低的原因。如果在不改变检测系统频率响应性能的情况下通过校正的方法使零点都位于原点附近,这样可以提高检测系统的动态响应性能。根据这样的思想,本文提出一种基于超前网络补偿的瞬时无功功率谐波检测方法。
27
华 中 科 技 大 学 博 士 学 位 论 文
2.4.2 基于超前网络补偿的瞬时无功功率谐波检测
通过上一节的分析,瞬时无功谐波检测系统的一个零点分布远离原点,导致检测系统动态响应变慢。如果在不改变系统极点的条件下,校正系统使零点都位于原点附近,将使检测系统具有较好的动态响应性能。超前网络环节具有微分性质,可以用来校正系统的零点。图2-2中低通滤波器主要是导致ip和iq的响应速度变慢的原因,超前网络环节具有微分的性质,能加快ip和iq的响应速度。基于超前网络补偿的ipiq检测方式原理图如图2-13所示。
kkτs+1T0s+1τs+1T0s+1图2-13 基于超前网络补偿的瞬时无功功率谐波检测原理图
图2-13中,超前校正环节为
Gc(s)=k
τs+1
T0s+1
(2.15)
为了加快ip、iq的响应速度,减少调节时间,在低通滤波器的输出端加入超前网络环节对ip、iq进行补偿,使ip、iq尽快达到新的稳定状态,从而达到改善检测系统的动态响应性能的目的 。
为了将检测系统远离原点的零点校正到原点附近,Gc(s)应具有微分的性质,对滞后的补偿主要由时间常数τ来补偿,它的取值可以根据需要校正的零点来选择。
1
可以消除微分的引入而产生的噪音干扰,为了保证环节Gc(s)具有超前补偿的T0s+1
作用,时间常数T0的选择应满足
T0<<τ (2.16) 增益k用来在实际系统中补偿超前环节Gc(s)对信号的衰减。超前环节Gc(s)的零极点分布关系及Bode图为图2-14
28
华 中 科 技 大 学 博 士 学 位 论 文
1τ−1T01T0−1τ90000
图2-14 超前环节零极点分布和bode图
它在1/τ~1/T0频率段对信号具有微分作用。
仍然以式(2.13)为例,分析瞬时无功功率谐波检测校正方法。在式(2.13)中主要是把远离原点的零点z2=−266.6校正到原点附近,设计校正环节
Gc(s)=
0.0075s+1
(2.17)
0.000075s+1
校正环节的零点为z=133.3,极点为p=13333.3,由图2-13可以得到校正后检测系统的传递函数为
Gk(s)=
(s+13560)(s+9.534)(s−7.442)
,k=a,b,c (2.18) 24
(s+13300)(s+266.6s+3.553*10)
它的零极点分布图为
图2-15 校正后检测系统零极点图
通过校正后,检测系统分别多了一个零点和极点。在图2-15中,有一对远离原点的零极点z3=−13560和p3=−13300,它们远离原点,并且相隔较近,是一对对偶零极点,对系统性能没有大的影响,因此可以忽略。系统的主导极点p1、p2没有变化,所以检测系统的频率响应特性没有变化。通过校正零点z2从-266.6移到-9.35,离原点较近,但是另外一个零点z1从原点移到了7.44,也在原点附近。经过校正后,
29
华 中 科 技 大 学 博 士 学 位 论 文
检测系统是一个非最小相角系统,对系统的相频特性有较大的影响,但是因为实际系统中校正环节处理的是直流信号ip、iq,所以不会对检测系统产生影响。图2-16是检测系统校正前和校正后的频率特性曲线对比。
700校正前-10幅值/dB-20-30-40-50-60-700.11101001k10k校正后频率/Hz图2-16 校正前后频率特性对比
幅值/A图2-17 校正前后阶跃响应对比
图2-16显示,经过超前网络校正后,系统的截止频率没有发生变化,但是,在截止频率附近,幅度响应曲线更加陡峭,检测系统具有更加理想的高通性能,表明经过校正后,系统有更好的频率特性。图2-17是校正前和校正后检测系统的阶跃响应性能对比。
通过校正后检测系统的阶跃响应上升时间从12ms减少到5ms,系统的动态响应性能明显提高,具有更好的动态响应性能。 2.4.3 仿真实验验证
在和2.3.4节相同条件的情况下对图2-7的电路中负载变化进行基于超前网络补偿的谐波检测仿真。图2-18分别是补偿前和补偿后ip、iq变化曲线对比。
图2-18 对变化电流检测的ip、iq及检测波形
30
华 中 科 技 大 学 博 士 学 位 论 文
从图中可以看到,在0.1s和0.2s时刻,负载发生变化。当负载电阻从10Ω变化到5Ω时,与未补偿前相比,通过补偿后,ip的变化明显加快,从30ms的动态响应时间提高到约15ms,大大缩短了响应时间,最终会提高谐波检测的响应时间。图2-19是电网电流、基于超前网络补偿后检测到的基波电流和谐波电流波形。
100.050.00.0-50.0-100.0100.050.00.0-50.0-100.060.0m80.0m0.10.120.140.160.180.20.220.24时间/s
图2-19 电网电流、基波电流及谐波电流检测波形
对比图2-10和图2-11,当负载电阻从10Ω变化到5Ω时,采用基于超前网络补偿的瞬时无功功率谐波检测,分解得到的直流信号ip、iq的调节时间从30ms减少到只有15ms,基波电流的响应明显加快,提高了检测系统的动态响应速度。同时,从图2-19中可以看到,改进的ipiq检测方式没有影响到稳定电流的检测精度。图2-10对电流变化的周期进行了局部放大,从中可以更加清楚的看到检测到的基波、谐波电流变化过程。
电网、基波电流/A谐波电流/A100.0(a)(b)0.0-100.0100.00.0-100.082.5m85.0m87.5m90.0m92.5m95.0m97.5m0.10.10250.1050.10750.110.11250.1150.1175时间/s
图2-20 电网电流、检测基波和谐波电流局部放大图
31
华 中 科 技 大 学 博 士 学 位 论 文
2.5 本章小结
本章主要研究了基于瞬时无功功率理论对电网谐波电流的实时检测。首先介绍了瞬时无功功率理论中ipiq谐波电流检测方式的原理,详细分析推导了瞬时无功功率谐波检测的传递函数,分析了其静态性能和动态响应性能。针对瞬时无功功率理论谐波检测动态性能较差的缺点,提出了一种基于超前网络补偿的瞬时无功功率实时谐波检测方法,通过校正检测系统零点,在保证检测系统截止频率不变的条件下改善其动态响应性能,分析和仿真结果证明该方法对检测系统的幅频特性和动态响应性能都有较大的改善。
在分析中主要考虑的是针对低通滤波器引起的延迟,在实际系统中,还有采样、计算时间等延迟,对系统的性能也会造成一定的影响,如果在通过超前网络校正检测系统的动态响应性能时,将实际系统的采样延迟和计算时间等因数也考虑进去,统一进行校正和补偿,可以有效的改善检测系统的动态响应性能。
瞬时无功功率理论是基于三相电路提出来的,它对三相电路的谐波检测可以取得较好的效果,但不适用于单相电路的谐波检测,如果要用于单相检测,必须对单相进行重新构造,增加了检测系统的复杂性[93][94]。同时,瞬时无功功率理论检测谐波时,需要用到电网电压信息,受电网的影响较大。现代数字信号处理技术新的理论和成果,提供了解决问题的方法和途径。下一章将对小波变换理论应用于实时谐波检测进行研究。
32
华 中 科 技 大 学 博 士 学 位 论 文 3 基于小波变换的谐波实时检测与补偿策略研究
3.1 引言
本章首先简要介绍小波变换的基本理论,包括连续小波变换、离散小波变换、多分辨率分析以及快速小波算法,然后在此基础上应用小波变换对电网电流进行谐波检测的原理及方法进行研究,并应用Mallat算法进行仿真实验,分析小波变换在实时检测中的不足之处以及产生误差的原因,提出一种基于小波系数的谐波电流实时检测和补偿方法,仿真结果表明该方法有效的减少检测误差,提高了检测的实时性。
3.2 小波变换理论基础
3.2.1 连续小波变换
设φ(t)为一平方可积函数,也即φ(t)∈L2(R),若其傅立叶变换ψ(ω)满足条件
∫
R
|ψ(ω)|2
ωdω<∞ (3.1)
则称φ(t)为一个基本小波或小波母函数,并称式(3.1)为小波函数的可容许性条件[95][96]。将小波母函数φ(t)进行伸缩和平移,就可以得到连续小波基函数
φa,b(t)=a−12φ(
t−b
),a∈R,a>0;b∈R (3.2) a
其中a为尺度因子,b为平移因子,它们都取连续变化的值。 连续小波变换(Continuous Wavelet Transform, CWT)的定义: 将任意L2(R)空间中的函数x(t)在小波基下展开,其表达式为
WTf(a,b)= −∞ +∞ t−b )dt (3.3) a 在容许条件3.1下,连续小波变换的逆变换公式为 x(t)= 1Cφ∫ +∞ 0 1t−bda+∞ (,)φ()db (3.4) CWTabx aa2∫−∞a 33 华 中 科 技 大 学 博 士 学 位 论 文 其中Cφ=∫ +∞ 0 |ψ(aω)|2 da a 从小波变换的定义可以得到以下结论 (1) 小波变换是信号x(t)与小波φa,b(t)的内积,是信号x(t)在基函数系φa,b(t)上的分解或投影; (2) 小波基具有尺度a、平移b两个参数; (3) 将信号x(t)在小波基下展开,就意味着将一个时间函数投影到二维的时间 -尺度相平面上; (4) 小波变换是一种变分辨率的时频联合分析方法,信号x(t)在某一尺度a、平移点b上的小波变换系数,实质上表征的是在位置b处,时间段aΔt上包含在中心频率为 ω0 a 、带宽为 Δω频率窗内的频率分量大小。 a 因此,小波变换是对信号在不同时间点上进行频率分解,然后将所需频率成分进行重构,从而得到原信号中需要成分的一种信号处理方法。 3.2.2 离散小波变换 连续小波变换是在连续变化的尺度因子a和平移因子b下进行的变换,变换基函数φa,b(t)具有很大的相关性,经过变换得到的连续小波系数CWTx(a,b)含有很多的冗余信息。离散小波变换(Discrete Wavelet Transform, DWT)是对尺度和平移参数进行取样,从而只用时间-尺度平面上的一些适当选取的离散点的小波变换值来描述信号[97],并由它们恢复原信号。实际应用中,一种常用的方法是对尺度a按幂级数作离散化,即取a=a0j,j=0,1,2,...,位移b可以在尺度a0下以间隔b0作均匀采样(b0应适当选择使信息仍能覆盖全b轴而不丢失),第j尺度下沿b轴以a0jb0为间隔作均匀采样,即 b=ka0jb0,k=0,1,2,...这样连续小波变换下的基函数系φa,b(t)就成为 −j a0φ[a(t−kab)]=a0φ[a0−kb0] (3.5) −j 2 −j0 j00 −j2 记作φj,k(t),j,k∈Z。取a0=2,b0=1,有 φj,k(t)=2φ(2−jt−k) (3.6) 所以,任意函数x(t)的离散小波变换为 − j2 34 华 中 科 技 大 学 博 士 学 位 论 文 WTx(j,k)=∫x(t)φj,k(t)dt (3.7) −∞+∞ 若离散小波序列{φj,k}j,k∋Z构成一个框架,其上、下边界分别为A和B,则由离散小波系数重建原信号x(t)的近似公式为 x(t)= 2 ∑WTx(j,k)⋅φj,k(t) (3.8) A+Bj,k A−B ||x||。为了减小逼近误差,使离散小波系数A+B 逼近误差的范数为||R||⋅||x||≤ 当A=B=1时,框架φj,k(t)成为一个紧框架,并且成为正交基,精确恢复原信号x(t), 此时经过框架变换后的信息没有冗余。在紧框架A=B=1离散小波变换的逆变换公式为 x(t)=∑ j,k j,k 3.2.3 多分辨率分析 多分辨率分析(Multi-resolution Analysis, MRA),也称多尺度分析是建立在函数空间概念上的理论,它是S.mallat在研究图像处理问题时引入小波分析中的,将标准正交小波基作为多尺度图像分析的一个工具,用它来描述当信号从一个较粗的逼近到一个较精细的逼近时所需要的信息增量。它的基本思想是将L2(R)空间分解为一系列具有不同分辨率的子空间,然后将L2(R)中的函数x(t)用各个子空间上的函数进行描述,即在各个子空间上进行投影,得到信号x(t)在不同分辨率子空间上的特征,从而可以在不同的尺度上由粗及精地观察信号[95]。 定义1 多分辨率分析是指满足下述性质的一系列闭子空间{Vj},j∋Z (1) 一致单调性:Vj⊂Vj−1,j∈Z (2) 渐进完全性:∩Vj={0};∪Vj=L2(R),j∈Z (3) 伸缩性:x(t)∈Vj⇔x(2jt)∈V0,j∈Z (4) 平移不变性:x(t)∈V0⇒x(t−n)∈V0,n∈Z (5) 正交基存在性:存在φ∈V0,使得{φ(t−n)}n∈Z是V0的正交基。 定义2 尺度空间{Vj}j∈Z的补空间 35 华 中 科 技 大 学 博 士 学 位 论 文 设Wj为Vj在Vj−1中的补空间,即 Vj−1=Vj⊕Wj,Wj⊥Vj ,Vj∈Vj−1,Wj∈Vj−1 (3.10)则 L2(R)=⊕Wj (3.11) j∈Z {Wj}j∈Z构成了L2(R)的一系列正交子空间,且任意子空间Wj之间是正交的。 由定义1和定义2可以知道,多分辨率分析的一系列尺度空间和补空间都是分别由同一尺度函数在不同的尺度下张成的,即一个多分辨率分析{Vj}j∈Z对应一个尺度函数(Scaling Function)φ(t),其补空间{Wj}j∈Z对应一个小波函数(Wavelet Function)ϕ(t)。 φj,k(t)=2φ(2−jt−k),j,k∈Z (3.12) ϕj,k(t)=2ϕ(2−jt−k),k∈Z (3.13) 从式(3.10)可知,如果φ(t)、ϕ(t)分别为尺度空间V0和小波空间W0的一个标准正交基函数,又由于V0⊂V−1,W0⊂V−1,所以φ(t)、ϕ(t)也必然属于V−1空间,那么φ(t)、ϕ(t)可用V−1空间的正交基φ−1,n(t)线性展开 −j 2 − j2 φ(t)= ϕ(t)= n=−∞+∞ ∑h(n)φ+∞ −1,n (t)=2∑h(n)φ(2t−n) (3.14) n=−∞+∞ +∞ n=−∞ ∑g(n)φ−1,n (t)=2∑g(n)φ(2t−n) (3.15) n=−∞ 式(3.14)、式(3.15)描述的是相邻二尺度空间基函数之间的关系,称为二尺度方程。系数h(n)=<φ,φ−1,n>,g(n)=<ϕ,φ−1,n>称为滤波器组系数,其中h(n)具有低通的性质,g(n)具有高通的性质。 3.2.4 Mallat算法 任意信号x(t)∈L2(R)在Vj−1空间的展开式为 x(t)= k=−∞ ∑c +∞ j−1,k 2(−j+1)2φ(2−j+1t−k) (3.16) 由式(3.10)有Vj−1=Vj⊕Wj,将x(t)分解一次投影到Vj和Wj空间,得到信号x(t)在Vj下的粗像和在Wj下的细节 36 华 中 科 技 大 学 博 士 学 位 论 文 x(t)= k=−∞ ∑c +∞ j,k 2 −j/2 φ(2t−k)+ −j k=−∞ ∑d +∞ j,k 2−j/2ϕ(2−jt−k) (3.17) 其中cj,k为j尺度上的剩余系数,dj,k为j尺度上的小波系数。式(3.17)中第一项为x(t)在j尺度空间的近似部分,第二项为x(t)在j尺度上的细节部分。根据式(3.14)和式(3.15)可以推得相邻尺度空间剩余系数和小波系数之间的关系为[96][98] cj,k= dj,k= m=−∞+∞ ∑h(m−2k)c +∞ j−1,m (3.18) (3.19) m=−∞ ∑g(m−2k)c j−1,m 式(3.18)和式(3.19)说明j尺度空间的剩余系数cj,k和小波系数dj,k可以通过 j−1尺度空间的剩余系数cj−1,k经过滤波器系数组h(n)和g(n)进行加权求和得到。如果将Vj空间的剩余系数cj,k进一步分解下去,可以得到Vj+1和Wj+1空间的剩余系数 cj+1,k和小波系数dj+1,k,继续分解下去,可以得到任意尺度J空间下的剩余系数和小波系数, 图3-1为Mallat分解快速算法的示意图。 图3-1 Mallat小波分解快速算法示意图 经过分解,可以得到信号x(t)的在各尺度空间的分解表达式为 x(t)= j=−∞K=−∞ ∑∑d J+∞ j,k (t)ϕj,k(t)+ k=−∞ ∑c +∞ J,k φJ,k(t) (3.20) 与小波分解类似,可以由各尺度空间的剩余系数和小波系数对剩余系数进行重构,小波变换系数的重建公式为 cj−1,m= k=−∞ ∑c +∞ j,k h(m−2k)+ k=−∞ ∑d +∞ j,k g(m−2k) (3.21) 图3-2是Mallat小波重构快速算法示意图。 图3-2 Mallat小波重构快速算法示意图 37 华 中 科 技 大 学 博 士 学 位 论 文 利用Mallat分解和重构算法进行信号处理时,通常假定原始信号x(t)的采样属于V0空间,作为c0,k。同时在实际应用Mallat算法时,由于实际信号和滤波器组都是有限长度的,因此计算时存在如何处理边界的问题,比较常用的方法是采用周期扩展或对称扩展。 3.2.5 正交小波基的性质分析 由多分辨率分析可以知道,尺度函数φ(t)和小波函数ϕ(t)满足二尺度方程式(3.14)和式(3.15),利用二尺度方程就可以构造出小波母函数,通过伸缩平移可以得到整个L2(R)空间的基。小波基具有以下数学特性: (1) 正交性 正交性(Orthonormality)是空间被正交分解的基础,从小波基的正交性出发,它所对应的低通滤波器h(n)和高通滤波器g(n)具有如下的关系 g(n)=(−1)n⋅h−n+1 (3.22) (2) 消失矩 小波不仅应当是零均值的,而且还必须具有高阶消失矩,如果ϕ(t)满足 ∫ +∞ −∞ ϕ(t)tmdt=0,(m=0,1,...M−1) (3.23) 则称小波ϕ(t)具有M阶消失拒(Cancellations),它使函数在小波展开时消去了其高阶平滑部分,使我们能够研究函数的高阶变化和某些高阶倒数中可能的奇异性。 (3) 正则性 正则性(Regularity)是函数光滑程度的一种描述,进行小波变换时,如果一段光滑的输入信号经过正则性较差的滤波器,其输出将出现不连续性,会导致出现误差,因此实际分析中希望小波具有一定的正则性。 (4) 紧支性 紧支性(Compact support)描述为如果小波函数ϕ(t)在区间[a,b]外恒等于零,那么 ϕ(t)在这个区间上紧支撑,称为紧支撑小波。支集越小的小波,其局部化的能力就 越强,应用精度越高。 (5) 对称性 如果小波函数ϕ(t)满足 ϕ(a+t)=ϕ(a−t)或ϕ(a+t)=−ϕ(a−t) (3.24) 38 华 中 科 技 大 学 博 士 学 位 论 文 则称ϕ(t)具有对称性(Symmetry)或反对称性。具有对称或反对称性的小波函数可以构造紧支的小波基,使其具有线性相位。对称性滤波器组在小波变换中可以避免信号在分解与重构中失真。 在应用小波变换进行信号处理时,可以根据实际需要结合以上的性质来合适选择合适的小波基。目前常用的小波函数主要有Haar小波、Daubechies小波、Symlets小波、Meyer小波、Morlet小波和Mexican Hat小波等,这些小波在紧支性、对称性、消失矩、正则性等方面均有不同的特点。 3.3 基于小波变换的电网谐波电流实时检测 3.3.1 谐波电流的特点 当电网给非线性负载供电时,电网电流中除含有基波电流外,还含有大量的谐波电流,谐波电流检测的目的就是要把谐波电流从电网电流中不失真的分离出来。按照谐波的定义,当基波频率为50Hz时,谐波频率有2次、3次…,n次谐波,考虑到普遍的情况,电网中最低谐波频率为100Hz。在实际应用中,电网中谐波含量主要为低次谐波,为20次谐波以下,高于20次的谐波幅值较小。因此,电网中谐波检测主要是2次到20次的谐波,即100Hz~1kHz的频率范围。在有源电力滤波器中,谐波检测的目的就是将50Hz的基波和各次谐波分离开。 从上一节小波理论分析中可以知道,小波变换可以将信号在各个小波空间进行投影,由于各个小波空间是相互正交的,并且每个小波空间代表着不同的频域,通过小波分解将信号影射到各个小波空间意味着将信号所包含的频率成分根据小波空间所划分的频域进行分离。如果事先将基波电流和谐波电流划分到不同的小波空间,通过小波变换就可以将基波电流和谐波电流分离开。基于小波变换的电网谐波电流检测主要分以下几个步骤 (1)划分小波空间,将基波频率成分划分在一个独立的小波空间中; (2)将电网电流经过式(3.20)的小波变换,得到电网电流中基波的小波系数和各频率空间的小波系数; (3)将基波小波系数按照式(3.21)单独进行重构,得到了电网中的基波电流; (4)用原电网电流减去基波电流可以得到电网中谐波总和,达到谐波电流检测的目的。 39 华 中 科 技 大 学 博 士 学 位 论 文 3.3.2 小波基的选择出发点 为了将基波电流和谐波电流完全分离,避免小波变换时发生混叠,同时在分解和重构过程中减少失真,必须选择合适的小波基。结合电网中谐波电流的特点和谐波检测的目的,基于小波变换的谐波电流检测所选择的小波基必须具有以下特征 (1)正交性。用来保证基波和谐波的完全分离。 (2)紧支性。用来保证检测的精度。 (3)对称性。用来避免信号分解与重构时的失真。 (4)高阶消失矩。用来保证有限的计算长度和计算量。 综合以上的特征,本文选择Daubechies小波db10,其长度为20,具有10阶消失矩。 3.3.3 电网电流的小波分解 式(3.18)、式(3.19)描述的系数分解可以看作一个输入离散序列cj−1,m经过双通道滤波器h(n)和g(n)的过程,其中h(n)具有低通的性质,离散序列经过h(n)后输出离散信号的低频概貌cj,m,g(n)具有高通性质,离散序列经过g(n)后输出离散信号的高频细节dj,m。因此,原离散信号经过双通道滤波器组后,原始信号的频带被等分为低频和高频两部分,带宽分别只有原信号的一半,由采样定理可知,可以将采样频率降低一半而不丢失信息,所以可以对通过滤波器组h(n)、g(n)的信号进行二抽取,使它们的数据长度减半,即对应cj,m和dj,m的数据长度分别只有原序列cj−1,m的一半,但两个滤波器输出序列总的长度与输入长度保持一致。信号n次分解的示意图为 g(-n)C0,k h(-n)2C1,k h(-n)2C2,k h(-n)2Cn,k 2d1,k g(-n)2d2,k g(-n)2dn,k 图3-3 信号n次分解示意图 图3-3中,信号每经过滤波器组h(n)、g(n)都进行2抽取,保证数据总长度与输入信号长度相等。 40 华 中 科 技 大 学 博 士 学 位 论 文 设电网电流信号x(t)经过采样后得到的离散序列xN,它的最高频率为f,序列将离散序列xN作为最高分辨率空间V0中的概貌c0,k,经过滤波器组h(n)和长度为N。 g(n)一次分解后得到频带为0~f/2的低频概貌c1,k和频带为f/2~f的高频细节d1,k,c1,k和d1,k的长度都为N/2。将低频概貌c1,k经过滤波器组h(n)和g(n)第二次分解,又可以得到频带为0~f/4的低频概貌c2,k和频带为f/4~f/2的高频细节d2,k, c2,k和d2,k的长度都为N/4,这样经过n次分解可以得到频带为0~f/2n的低频概貌cn,k和频带为f/2n~f/2n−1的高频细节dn,k,它们的长度都为N/2n,直到0~f/2n的低频概貌中只包含基波电流成分,cn,k就是基波的小波分解系数。 3.3.4 基波电流的重构 电网电流信号x(t)经过图3-3所示的分解,得到了不同频带内的小波系数 d1,k,d2,k,⋅⋅⋅dn,k和cn,k,其中cn,k代表x(t)中的基波分解系数,将cn,k单独进行重构,可以得到x(t)中的基波电流成分。图3-4是由基波系数cn,k重建基波电流的示意图。 00 0C2,kCn,k2h(-n)2h(-n)C1,k2h(-n)C0,k 图3-4 信号n次重建示意图 与图3-3相对应,在重构的过程中为了保证恢复原信号的数据长度,将小波系数进行了插2处理。 3.3.5 仿真分析 在saber仿真环境中建立图2-7的仿真模型,其中,负载电阻R=10Ω,负载电感L=1mh,采用小波变换进行谐波检测。仿真实验中采样频率4.8khz,分解层数为 5层,采用db10小波,边界扩展采用周期扩展。根据Nyquist采样定理,原始所包含的信号频带为0~2.4khz,最高谐波为48次谐波。 根据Mallat算法的特点,仿真算法分为以下几步 (a)确定小波空间。根据谐波成分应将采样信号ila进行5层小波分解。基波电流系数位于第5层低频概貌中。 41 华 中 科 技 大 学 博 士 学 位 论 文 (b)进行5层小波分解。将ila作为V0空间的低频概貌c0,k,按照图3-3进行5层分解,得到不同频率段的小波变换系数,它们所代表的频率如表3-1所示 表3-1 a相电网电流5层分解的小波系数和对应频带 小波系数 频带(Hz) d1,k 1.2k~2.4k d2,k 600~1.2k d3,k 300~600 d4,k 150~300 d5,k 75~150 c5,k 0~75 图3-5是电网电流采样信号ila波形和 表3-1中小波系数代表的频率成分波形图。 幅值/A 图3-5 电网电流及各分解后频率成分 图3-5中c5、d5…d1波形分别对应表3-1中c5,k、d5,k…d1,k的频率成分。c5即为分解出电网电流中的基波电流成分。 (c)将基波小波系数c5,k单独重构,得到基波电流c5。 (d)用电网电流ila减去基波电流c5得到谐波电流。 小波变换谐波电流实时检测,当对电网电流进行小波分解和基波电流重构算法时,需要对当前时刻数据进行边界扩展。对于电网电流这种周期性的信号,适合采用周期扩展的方法,即用上一个周期的信号来代替下一段时间的信号(长度由小波 42 华 中 科 技 大 学 博 士 学 位 论 文 基的长度决定)。如果电网电流是稳定的,周期扩展可以取得很好的检测精度,但是当电网电流发生变化时,周期扩展会带来检测误差。 (1)恒定负载下的谐波电流检测 对负载恒定的电网电流检测可以验证检测精度。保持图2-7中负载恒定,采用小波变换对电网电流进行谐波实时检测。图3-6(a)、(b)、(c)分别是电网电流、检测出的基波电流和谐波电流波形,图3-7分别它们的频谱图,表3-2列出了它们的幅值及检测误差。 谐波电流/A基波电流/A电网电流/A(c)(b)(a) 图3-6 小波变换检测到的基波和谐波电流 (c)(b)(a)谐波电流/A基波电流/A电网电流/A60.030.00.060.030.00.060.030.00.001002003004005006007008009001k1.1k1.2k1.3k1.4k1.5k1.6k1.7k1.8k1.9k2k频率/Hz 图3-7 小波变换谐波检测频谱图 表3-2 基波和主要谐波的幅值及检测误差 50hz 5次 电网电流(A)基波电流(A)谐波电流(A)误差(%) 56.588 0.048 7次 11次 13次 0 0 0 56.666 12.575 6.318 4.733 3.997 0 12.625 6.388 4.763 4.041 0 0.39 1.1 0.62 1.10 43 华 中 科 技 大 学 博 士 学 位 论 文 表3-2中检测得到的基波电流THD为0.92%。从表中可以看到,检测出来的基波电流基本上不包含谐波电流,表明从电网电流中将基波电流完全分离出来了。对主要的几次谐波检测误差最大为7次谐波,误差为1.1%,仿真结果表明基于小波变换的谐波实时检测具有很高的检测精度。 (2)负载变化情况下谐波电流检测 在电网实际运行过程中,负载往往是变化的,对变化的电网电流谐波检测可以观察检测动态性能。在其它条件不变的情况下,将图2-7中负载电阻从10Ω变化到 5Ω,再变化到10Ω,对变化的电网电流进行小波变换谐波实时检测。图3-8是谐波 检测的波形图。 电网、基波电流/A谐波电流/A100.050.00.0-50.0-100.0100.050.00.0-50.0-100.060.0m80.0m0.10.120.140.160.180.20.220.24(b)(a)时间/s 图3-8 小波变换对突变电流谐波检测 在图3-8中,为了清楚的研究小波变换谐波实时检测的动态性能,将电网电流和检测到的基波电流放在同一坐标轴中。电网电流在0.1s和0.2s时发生变化,从图中可以看到,在电流变化后的一个周期中检测出的基波电流都有误差,说明基于小波变换的谐波实时检测的动态响应性能较差。造成动态响应性能较差的原因是电网电流小波分解和基波电流重构算法时需要进行周期扩展。为了改善由于边界延拓引起检测过程的动态性能变坏,下面研究一种补偿策略。 3.4 基于小波系数的电网谐波电流检测实时补偿方法研究 3.4.1 小波分解层数与小波系数长度的关系 电网电流由基波和谐波成分组成,设其时域表达式为 44 华 中 科 技 大 学 博 士 学 位 论 文 x(t)=x0(t)+xh(t)=Asin(2πf0t+ϕ)+∑Aksin(2kπf0t+ϕk) (3.25) k=2n 其中x0(t)和xh(t)分别表示电网电流中基波电流和谐波电流成分,A、f0和ϕ分别为电网电流中基波电流的幅值、频率和相位,Ak和ϕk分别为第k次谐波的幅值和相位。 用采样频率为fs对信号x(t)进行采样,得到c0,k,则在一个基波周期中c0,k的长度为n0=fs/f0。设经过n层小波分解可以将基波成分和谐波成分分离开,第n层概貌cn,k中只包含基波电流系数。从图3-3可以知道,信号每经过一层小波分解,分解系数都要进行2抽取,小波系数长度减半,因此对x(t)进行n层分解得到第n层的剩余系数cn,k的数据长度为 nn=n0/2n=fs/(2n⋅f0) (3.26) 采样频率越高,信号x(t)中包含的最高频率成分也越高,为了将信号f(t)中的基波和其它频率成分分离开,分解的层数也必须越多。设经过n层分解,第n层小波分解系数cn,k只包含基波成分,由式(3.26)可以知道,cn,k的数据长度主要由采样频率、基波频率及分解层数确定,而分解层数又和采样频率是对应的,因此在某一个采样频率下,一个基波周期中的数据c0,k经过n层分解后得到的基波的小波系数 cn,k的长度nn是固定的。表3-2列出了在不同的采样频率下采样频率、一个周期中数据采样点、分解层数和一个周期中基波小波分解系数长度nn之间的对应关系。 表3-3 采样频率、小波分解层数和一个周期基波分解系数长度 采样频率(hz) 一个周期采样点数n0 分解层数n 一个周期基波分解系数长度nn 4.8k 9.6k 19.2k 96 192 384 5 6 7 3 3 3 从表3-3可以看出,在不同的采样频率下,为了将基波频率和谐波分离,进行小波分解的层数也不同,最终得到的一个周期中基波变换系数的长度都是3。进一步分析可以得到,由于对信号进行n层小波分解时,由于每次抽2都是均匀抽取,因此得到的第n层基波系数的3个数据点在一个基波周期中是均匀分布的,即三个点之间后一点和前一点相位差120º。 45 华 中 科 技 大 学 博 士 学 位 论 文 3.4.2 基波电流变化与小波系数的关系 小波变换是一种时频联合的分析方法,它可以精确分析某个时刻的信号频率成分,当信号发生变化时,在小波系数上也能够精确的反映这种变化,利用小波变换的这个特征对谐波实时检测进行补偿。下面首先分析基波电流幅值的变化与小波系数的关系。 设对电网电流x(t)进行n层小波分解,可以将基波电流和谐波电流分离,n层小波分解表达式为 x(t)=∑cn,kφn,k(t)+∑∑dj,kϕj,k (3.27) k=1 j=1k=1 3nnn 其中第一项为基波电流的小波分解,第二项为谐波电流成分。为了分析电流变化时基波电流变化在小波系数上的特征,对电网电流中一个周期基波进行单独的分析,由式(3.25)和式(3.27)可以得到基波电流表达式