CN109254202B - 一种应用于配电网的同步相量测量装置 - Google Patents
一种应用于配电网的同步相量测量装置 Download PDFInfo
- Publication number
- CN109254202B CN109254202B CN201810954532.7A CN201810954532A CN109254202B CN 109254202 B CN109254202 B CN 109254202B CN 201810954532 A CN201810954532 A CN 201810954532A CN 109254202 B CN109254202 B CN 109254202B
- Authority
- CN
- China
- Prior art keywords
- synchronous
- phasor measurement
- distribution network
- signal
- kaiser
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000001360 synchronised effect Effects 0.000 title claims abstract description 29
- 238000005259 measurement Methods 0.000 title claims abstract description 28
- 238000005070 sampling Methods 0.000 claims abstract description 12
- 238000000034 method Methods 0.000 claims description 5
- 238000001914 filtration Methods 0.000 claims description 3
- 230000010355 oscillation Effects 0.000 claims description 3
- 238000000691 measurement method Methods 0.000 claims 1
- 230000001629 suppression Effects 0.000 abstract description 3
- 238000013461 design Methods 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R25/00—Arrangements for measuring phase angle between a voltage and a current or between voltages or currents
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/22—Flexible AC transmission systems [FACTS] or power factor or reactive power compensating or correcting units
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
本发明公开了一种应用于配电网同步相量测量装置，包括：低通滤波器、模数转换器、数字信号处理器、微处理器、同步采样时钟单元、工频振荡器和输出、显示、存储单元。本发明的DSP芯片中，将嵌入基于Kaiser加权卷积自乘窗的同步相量测量算法，它具有很强的谐波、间谐波和噪声干扰抑制能力，从而更适用于配电网中；同时采用工频振荡器调整输出结果值，使输出结果准确度更高。
Description
技术领域
本发明涉及配电网运行状态监测领域，具体涉及一种应用于配电网的同步相量测量装置。
背景技术
为实现电力系统的动态安全监测，快速发现电力系统中可能存在的故障[1]，动态同步相量(Dynamic Synchrophasor)的测量起到了重要的作用。电气与电子工程师协会于2011年发布了修订版的IEEE C37.118.1标准。该标准对不同动态条件下相量测量单元(Phasor Measurement Unit,PMU)的测量准确度提出了明确的指标。
目前，PMU被大量应用在输电网中。随着大量新能源(如光伏设施和储能装置等)接入配电网中，配电网的潮流方向从原来的单向变成双向。且这些新能源装置使配电网呈现强烈的不确定性和动态性。为了能实现配电网的实时监测，未来PMU也被考虑应用到配电网中。但由于大量电力电子装置的使用，配电网中还含有丰富的谐波、间谐波分量。而且配电网受噪声影响明显。如何在谐波和噪声含量丰富时，仍能准确地测量动态同步相量，成为一个亟须解决的问题。
论文“基于DFT的电力系统相量及功率测量新算法”(王茂海等，电力系统自动化，2005)通过对稳态条件下的DFT及其改进算法测量相量的误差产生机理进行分析，从而对测量得到的幅值和相位进行修正。该算法具有较高的准确度，但是在动态条件下，该算法并不适用。论文“Dynamic Phasor and Frequency Estimates Through Maximally FlatDifferentiators”(PLATAS-GARZA等，IEEE Transactions on Instrumentation&Measurement，2010)将动态同步相量用泰勒级数表示，并基于动态同步相量与信号之间的关系近似表示实际信号，最后通过最小二乘法求解得到同步相量。该算法原理清晰，但由于该算法在信号建模时，未考虑当信号中存在丰富的谐波和间谐波时，算法的准确度下降。针对这一问题，论文“实现动态相量测量的FIR数字滤波器最优设计”(汪芙平等，中国电机工程学报，2014)在对信号建模时，将其表示成动态同步相量和指定次谐波和间谐波的叠加，以此实现该次谐波和间谐波的陷波设计。但是该方法需要提前已知谐波、间谐波的幅值、频率、相位等参数。然而，在实际情况下，这些参数均是未知的。论文“动态条件下的同步相量测量算法的研究”(麦瑞坤等，中国电机工程学报，2009)应用泰勒级数表示动态同步相量，并对该动态同步相量以及原始信号作傅里叶变换，再应用不同时刻的结果联立求解后得到同步相量。该算法在一定程度上提高了动态同步相量测量的准确度，但是该算法计算量大，且未考虑间谐波的影响。
发明内容
为此，本发明的目的在于设计一种应用于配电网的同步相量测量装置，旨在实现间谐波、谐波和噪声的强烈抑制，从而在提高配电网相量测量的准确度。
为了实现上述目的，本发明的实施例公开了一种应用于配电网的同步相量测量装置，包括以下部分：
用于滤除配电网电压电流信号中高频成分的低通滤波器；
用于将模拟信号转换成数字信号的模拟数字转换器；
用于计算同步相量、频率、幅值和相位的数字信号处理器；
用于控制整个同步相量测量装置的单片机；
用于产生同步采样时钟的同步采样时钟单元；
用于产生工频振荡信号的工频振荡器；
输出、显示和存储单元；
所述的低通滤波器、模拟数字转换器、数字信号处理器、单片机依次连接，单片机的输出端连接输出、显示和存储单元；同步采样时钟单元分别与模拟数字转换器和工频振荡器连接；所述的工频振荡器与数字信号处理器连接。
本发明所采用的Kaiser加权卷积自乘窗FIR滤波器，相对于未采用窗函数加权的测量测量滤波器器，具有优异的旁瓣衰减特性，这能在很大程度上实现谐波和间谐波的抑制，从而提高谐波和间谐波含量丰富时的动态同步相量测量准确度。另外，上述滤波器在基波频率附近呈现出一种平滑特性，即在一定的频率范围内，信号的增益均为1，这正是实现带限信号动态同步相量测量的根本原理。
附图说明
图1是本发明的一种应用于配电网的同步相量测量装置的方框图。
具体实施方式
下面通过参考附图描述的实施例是示例性的，仅用于解释本发明，而不能理解为对本发明的限制。
在本发明的描述中，需要理解的是，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性。
参照下面的描述和附图，将清楚本发明的实施例的这些和其他方面。在这些描述和附图中，具体公开了本发明的实施例中的一些特定实施方式，来表示实施本发明的实施例的原理的一些方式，但是应当理解，本发明的实施例的范围不受此限制。相反，本发明的实施例包括落入所附加权利要求书的精神和内涵范围内的所有变化、修改和等同物。
图1是本发明一个实施例的一种应用于配电网的同步相量测量装置的方框图，包括以下部分：用于滤除配电网电压电流信号中高频成分的低通滤波器；用于将模拟信号转换成数字信号的模拟数字转换器；用于计算同步相量、频率、幅值和相位的数字信号处理器；
用于控制整个同步相量测量装置的单片机；用于产生同步采样时钟的同步采样时钟单元；用于产生工频振荡信号的工频振荡器；输出、显示和存储单元；所述的低通滤波器、模拟数字转换器、数字信号处理器、单片机依次连接，单片机的输出端连接输出、显示和存储单元；同步采样时钟单元分别与模拟数字转换器和工频振荡器连接；所述的工频振荡器与数字信号处理器连接。
在本发明的一个实施例中，所述的数字信号处理器嵌入了Kaiser加权卷积自乘窗FIR滤波器，为使本领域技术人员进一步理解本发明，将通过以下实施例进行进一步地说明。
(1)基于Kaiser加权卷积自乘窗的相量测量算法(FIR滤波器)的设计原理为：
动态同步相量被定义为一个频带远小于基波频率的带限信号x(t)，其可以表示为
其中，f0为基波频率，为动态同步相量，a(t)为波动的信号幅度，为波动的相位值。在很小的一段时间T内，可将动态同步相量近似表示成κ阶泰勒级数，其表达式为
其中，ψ(κ)(t)表示ψ(t)的κ阶导数。这样，在任意一个时间段T内，带限信号x(t)就可以近似表示为
其中，对上述连续信号进行离散采样，则式(5)可表示成矩阵形式
xκ＝Eκψκ
其中，
xκ＝[xκ(-N),...,xκ(0),...,xκ(N)]T
ψκ＝[ψκ *,ψκ-1 *,...,ψ0 *,ψ0,...,ψκ+1,ψκ]T
其中，ψk＝ψ(k)(0)/(k！(Ts k))，Ts为采样间隔，k＝0,1,...,κ。
其中，n＝-N,…,0,…,N，时间窗长度Nw＝2N+1。这样，就实现了动态同步相量的近似表示，即x≈xκ＝Eκψκ。接着，可用最小二乘法实现参考时刻(n＝0)的同步相量估计。但是当信号中含有丰富的谐波和间谐波分量时，上述方法存在较大误差。这里考虑采用Kaiser加权卷积自乘窗的方法实现谐波和间谐波分量的抑制。即
Wx≈Wxκ＝WEκψκ
其中，W＝diag{w2(n)}。算子diag{w2(n)}表示矩阵的对角化运算。这样，ψκ矩阵的线性估计最优解为其最小二乘估计，即
ψκ＝(Eκ HWHWEκ)-1Eκ HWHWx
从而可得矩阵(Eκ HWHWEκ)-1Eκ HWHW中，与同步相量ψ0对应的行向量即为所设计的FIR滤波器。
(2)Kaiser加权卷积自乘窗
这里设计的Kaiser加权卷积自乘窗的表达式为
其中，wkaiser代表Kaiser窗函数，A为归一化系数，其为中的所有元素之和，/>表示卷积运算。
其中，Kaiser窗函数的表达式为：
式中，β是其形状参数；I0(β)是第1类变形零阶贝塞尔函数；
Claims (1)
1.基于Kaiser加权卷积自乘窗FIR滤波器的同步相量测量方法，该方法包括应用于配电网的同步相量测量装置，所述的同步相量测量装置包括以下部分：
用于滤除配电网电压或电流信号中高频成分的低通滤波器；
用于将模拟信号转换成数字信号的模拟数字转换器；
用于计算同步相量、频率、幅值和相位的数字信号处理器；
用于控制整个同步相量测量装置的单片机；
用于产生同步采样时钟的同步采样时钟单元；
用于产生工频振荡信号的工频振荡器；
输出、显示和存储单元；
所述的低通滤波器、模拟数字转换器、数字信号处理器、单片机依次连接，单片机的输出端连接输出、显示和存储单元；同步采样时钟单元分别与模拟数字转换器和工频振荡器连接；所述的工频振荡器与数字信号处理器连接；
所述的数字信号处理器嵌入了Kaiser加权卷积自乘窗FIR滤波器；
其特征在于同步相量测量方法包括下述步骤：
S01：将动态电压或电流信号用有限项泰勒级数表示；
S02：对表示后的信号和原信号均施加Kaiser加权卷积自乘窗函数；
S03：表示后的信号与原信号近似相等，从而生成同步相量测量的频率滤波器；
S04：应用上述滤波器测算同步相量和频率；
所述的Kaiser加权卷积自乘窗，其表达式为：
其中，wkaiser代表Kaiser窗函数，A为归一化系数，其为中的所有元素之和，/>表示卷积运算；
步骤S01所述的将动态电压或电流信号用有限项泰勒级数表示，通用表达式为：
其中，ψ(k)(t)表示ψ(t)的k阶导数；
所述的工频振荡器，通用表达式为：
其中，f0表示电力系统频率。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810954532.7A CN109254202B (zh) | 2018-08-21 | 2018-08-21 | 一种应用于配电网的同步相量测量装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810954532.7A CN109254202B (zh) | 2018-08-21 | 2018-08-21 | 一种应用于配电网的同步相量测量装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109254202A CN109254202A (zh) | 2019-01-22 |
| CN109254202B true CN109254202B (zh) | 2024-03-26 |
Family
ID=65049628
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810954532.7A Active CN109254202B (zh) | 2018-08-21 | 2018-08-21 | 一种应用于配电网的同步相量测量装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109254202B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111830348B (zh) * | 2020-07-21 | 2023-05-02 | 国网四川省电力公司电力科学研究院 | 一种低频振荡下的动态同步相量测量方法和装置 |
| CN112433093B (zh) * | 2020-11-06 | 2023-08-22 | 国网四川省电力公司电力科学研究院 | 一种基于Kaiser窗的动态相量估计方法及系统 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5222189A (en) * | 1989-01-27 | 1993-06-22 | Dolby Laboratories Licensing Corporation | Low time-delay transform coder, decoder, and encoder/decoder for high-quality audio |
| CN101510687A (zh) * | 2009-03-18 | 2009-08-19 | 天津大学 | 电网中利用窗函数实现多采样率信号的频率转换方法 |
| CN105182073A (zh) * | 2015-09-06 | 2015-12-23 | 中国电力科学研究院 | 一种同步相量测量装置的动态相量测量系统及其测量方法 |
| CN107966672A (zh) * | 2017-12-05 | 2018-04-27 | 国网重庆市电力公司电力科学研究院 | 计及同步相量测量的发电机功率变送器测试方法及系统 |
| CN209102804U (zh) * | 2018-08-21 | 2019-07-12 | 国网辽宁省电力有限公司铁岭供电公司 | 一种应用于配电网的同步相量测量装置 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7372875B2 (en) * | 2002-09-30 | 2008-05-13 | Lucent Technologies Inc. | Systems and methods for synchronization in asynchronous transport networks |
| CN103575980B (zh) * | 2012-07-26 | 2016-12-21 | 施耐德电器工业公司 | 系统频率测量方法、同步相量测量方法及设备 |
-
2018
- 2018-08-21 CN CN201810954532.7A patent/CN109254202B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5222189A (en) * | 1989-01-27 | 1993-06-22 | Dolby Laboratories Licensing Corporation | Low time-delay transform coder, decoder, and encoder/decoder for high-quality audio |
| CN101510687A (zh) * | 2009-03-18 | 2009-08-19 | 天津大学 | 电网中利用窗函数实现多采样率信号的频率转换方法 |
| CN105182073A (zh) * | 2015-09-06 | 2015-12-23 | 中国电力科学研究院 | 一种同步相量测量装置的动态相量测量系统及其测量方法 |
| CN107966672A (zh) * | 2017-12-05 | 2018-04-27 | 国网重庆市电力公司电力科学研究院 | 计及同步相量测量的发电机功率变送器测试方法及系统 |
| CN209102804U (zh) * | 2018-08-21 | 2019-07-12 | 国网辽宁省电力有限公司铁岭供电公司 | 一种应用于配电网的同步相量测量装置 |
Non-Patent Citations (9)
| Title |
|---|
| A Hilbert-Huang based Approach for On-line Extraction of Modal Behavior from PMU data;F.Lezama Zarraga 等;41st North Amercian Power Symposium;1-6 * |
| Design and FPGA Implementation of Linear FIR Low-pass Filter Based on Kaiser Window Function;GAO Jinding 等;2011 Fourth International Conference on Intelligent Computation Technology and Automation;496-498 * |
| Dynamic Phasor and Frequency Estimates Through Maximally Flat Differentiators;Miguel A. Platas-Garza等;IEEE Transactions on Instrumentation and Measurement;20101231;第59卷;1803-1811 * |
| Dynamic Phasor Estimates through Maximally Flat Differentiators;Miguel Angel Platas Garza等;2008 IEEE Power and Energy Society General Meeting - Conversion and Delivery of Electrical Energy in the 21st Century;20081231;1-8 * |
| 同步相量测量装置的以太网通信接口的设计;刘燕;王峰;陈金鹏;;计算机测量与控制(第08期);1110-1112 * |
| 基于GPS的同步采样及在保护与控制中的应用;高厚磊，贺家李，江世芳;电网技术;第19卷(第07期);30-32 * |
| 基于改进泰勒加权最小二乘法的相量测量算法;曹磊 等;电力系统保护与控制;第46卷(第7期);24-30 * |
| 应用数字信号处理技术测量电力系统相量;齐俊, 赵玲, 吴文军;华东电力(第05期);13-14 * |
| 面向主动配电网的同步相量测量算法综述;陈毅阳;金涛;;电气开关(第01期);4-7 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109254202A (zh) | 2019-01-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zhan et al. | A Clarke transformation-based DFT phasor and frequency algorithm for wide frequency range | |
| CN103454497B (zh) | 基于改进加窗离散傅立叶变换的相位差测量方法 | |
| Wen et al. | Triangular self-convolution window with desirable sidelobe behaviors for harmonic analysis of power system | |
| CN105137185B (zh) | 一种基于离散傅里叶变换的频域插值电力谐波分析方法 | |
| Tomic et al. | A new power system digital harmonic analyzer | |
| CN109521275B (zh) | 一种同步相量确定方法、系统、装置及可读存储介质 | |
| CN105044456B (zh) | 一种基于正交子带的电网瞬时频率测量与跟踪方法 | |
| WO2016138764A1 (zh) | 一种提升pmu同步相量测量精度的方法 | |
| Wen et al. | Hanning self-convolution window and its application to harmonic analysis | |
| CN105548739B (zh) | 一种避雷器运行状态信号处理方法 | |
| CN107543962B (zh) | 主导间谐波频谱分布的计算方法 | |
| CN102323477A (zh) | 一种电网间谐波测量方法及测量仪 | |
| CN109254202B (zh) | 一种应用于配电网的同步相量测量装置 | |
| CN111693775A (zh) | 一种输电网的谐波检测方法、装置和介质 | |
| CN209102804U (zh) | 一种应用于配电网的同步相量测量装置 | |
| Tarasiuk | Comparative study of various methods of DFT calculation in the wake of IEC Standard 61000-4-7 | |
| CN111830348B (zh) | 一种低频振荡下的动态同步相量测量方法和装置 | |
| CN103543331A (zh) | 一种计算电信号谐波和间谐波的方法 | |
| CN111505375A (zh) | 一种频率测量方法及装置 | |
| CN108334822B (zh) | 基于电动汽车充电非线性负荷特征的卡尔曼和修正小波变换滤波方法 | |
| Kušljević | A simultaneous estimation of frequency, magnitude, and active and reactive power by using decoupled modules | |
| Radil et al. | Methods for estimation of voltage harmonic components | |
| Ferrero et al. | Employment of interpolated DFT-based PMU algorithms in three-phase systems | |
| Grando et al. | A method for synchronized harmonic phasor measurement based on hardware‐enhanced Fast Fourier Transform | |
| Chen et al. | A phase difference measurement method based on the extended kalman filter for Coriolis mass flowmeters |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |