CN101701985A - 定频变点电网谐波检测方法及其测量仪 - Google Patents
定频变点电网谐波检测方法及其测量仪 Download PDFInfo
- Publication number
- CN101701985A CN101701985A CN200910232708A CN200910232708A CN101701985A CN 101701985 A CN101701985 A CN 101701985A CN 200910232708 A CN200910232708 A CN 200910232708A CN 200910232708 A CN200910232708 A CN 200910232708A CN 101701985 A CN101701985 A CN 101701985A
- Authority
- CN
- China
- Prior art keywords
- sampling
- dft
- constant
- frequency
- zero
- 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.)
- Granted
Links
Images
Landscapes
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
本发明公布了定频变点电网谐波检测方法及其测量仪,测量仪包括两大特殊构件即定频整周采样器和变点数DFT计算器。所述检测方法包括如下步骤:(1)定频整周采样:以恒定采样频率fS对电网电压(电流)波形采样,由过零检测器控制每个数据窗口无限逼近S个整周波;(2)变点DFT计算:电网频率波动时每个数据窗口的采样点数是可变的,能够高速计算不同点的DFT,‘高速’的根源主要来自:采用高速DSP、执行三角函数型DFT、用查表法获得不同点的DFT系数、只计算IEC61000-4-30标准要求的500个频点。本发明采用定频变点法实现同步,同步偏差一定小于1/L0(L0=S×fS/55),只要选取足够大的fS和S采样同步偏差就足够小,不存在时延,减少了谐波测量同步误差。
Description
技术领域
本发明涉及电力***谐波测量仪器和谐波测量方法,属于电能质量监测技术领域。
背景技术
电力***谐波测量是电力***电能质量监测的重要内容。电力***谐波测量的常用方法是FFT(快速傅立叶变换)法,其基本原理是:首先对电网波形采样和ADC(模数转换)得到采样序列,然后对采样序列进行FFT得到各次谐波、间谐波的幅值和相位信息。FFT是DFT(离散傅立叶变换)的快速算法,在数据点数很多时,FFT比标准DFT快得多,但FFT要求数据点数必须是2的整次幂。
理论上,FFT方法要求‘同步采样’,即‘等间隔整周期采样’——同一数据窗口数据的采样周期相等且每个数据窗口的时宽刚好等于整数个电网周波。由于电网频率是不断变化的,‘同步采样’很难实现,所以,数据总存在‘采样同步偏差’,‘采样同步偏差’直接造成各次谐波、间谐波的幅值和相位的测量误差——称为‘谐波测量同步误差’。目前减少采样‘采样同步偏差’的方法主要有:(1)硬件同步——主要指加‘PLL(锁相环)’的同步方法;(2)软件同步——利用DSP的定时器实时改变采样频率。
‘硬件同步’和‘软件同步’的采样频率都是变化的,但每一数据窗口内的采样点数相同且都是2的整次幂。
PLL同步是IEC61000-4-30标准推荐的方法,IEC61000-4-30标准要求PLL的采样同步偏差不大于万分之三。电网频率恒定不变时PLL做到万分之三就不容易了,更何况电网频率时刻都在波动——PLL时刻都处在闭环调节的动态过程之中,这种闭环调节还存在时延,所以实际应用中PLL很难保证采样同步偏差不大于万分之三——传统的‘硬件同步’和‘软件同步’在实际应用中都不容易获得高精度。
IEC61000-4-30(2003)标准还要求FFT加矩形窗,Matlab仿真表明:同步偏差为万分之三时,矩形窗FFT对40次以上的谐波、间谐波的计算结果就是超标的。
软件同步可以看作是一个‘软件PLL’,除了上述弊端,‘软件同步’比‘硬件PLL’的时延大,因为后者是‘锁相’、前者是‘锁周’——用前一个周波的宽度确定下一个周波的采样周期——时延至少是一个周波。软件同步’还有个问题就是定时器的分数定时误差。
发明内容
本发明目的是针对现有的谐波测量仪同步误差很难减小的问题,提出一种谐波测量仪结构——定频变点谐波测量仪。
本发明的另一目的是提出一种定频变点谐波测量仪的检测方法。
定频变点谐波测量仪包括两大特殊构件:(1)定频整周采样器:以恒定采样频率fS对电网电压(电流)波形采样,由过零检测器控制每个数据窗口‘无限逼近’S个整周波;(2)变点DFT计算器:电网频率波动时每个数据窗口的采样点数是可变的,设计了一种‘变点DFT计算器’——能够高速计算不同点的DFT,‘高速’的根源主要来自:采用高速DSP、执行三角函数型DFT、用查表法获得不同点的DFT系数、只计算IEC61000-4-30标准要求的500个频点(谐波和间谐波)。
‘变点DFT计算器’的基本特征是:(1)每个数据窗口DFT的点数是变动的且一般不是2的整次幂;(2)DFT的执行用三角函数型DFT(离散富里叶变换),而不是FFT(快速富里叶变换);(3)不同点数的DFT系数用‘查表法’获得;(4)只计算IEC-4-30要求的500个频点(谐波和间谐波)——剔除了大量的不需要计算的频点——大大减少DFT的计算量。
所述的定频变点电网谐波测量仪的检测方法之一,其特征在于包括如下步骤:
(1)采用可编程片上***建立DFT系数表;
(2)电网三相电流分别经电流互感器送到模/数转换器,电网三相电压分别经电压互感器送到模/数转换器和过零检测器;
(3)S个电网周波的开始(电网波形过零)DSP通过模/数转换器以恒定采样频率fS采当前数据窗口的第1点,S个电网周波的结束(电网波形过零)停止本数据窗口的采样,每个所述采样数据窗口的采样点数为Lk、宽度等于S个电网整周波宽度,Lk、S为大于零的自然数;
(4)DSP根据实时采样数据及采样点数,然后查步骤(1)所述的DFT系数表获得当前测量的谐波参数。
定频变点电网谐波测量仪的检测方法之二,其特征在于包括如下步骤:
(a)采用可编程片上***建立DFT系数表;
(b)电网三相电流分别经电流互感器送到模/数转换器,电网三相电压分别经电压互感器送到模/数转换器和过零检测器;
(c)S个电网周波的开始(电网波形过零)DSP通过模/数转换器以恒定采样频率fS采当前数据窗口的第1点,S个电网周波的结束(电网波形过零)停止本数据窗口的采样,每个所述采样数据窗口的采样点数为Lk、宽度等于S个电网整周波宽度,Lk、S为大于零的自然数;
(d)将步骤(c)所述的采样数据经过可编程片上***三角函数型离散富里叶变换DFT,然后查步骤(a)所述的DFT系数表得到谐波参数。
定频变点电网谐波测量的采样同步偏差一定小于其中,(round{}表示取整运算),所以fS越大、S越大、采样同步偏差越小——这就是‘无限逼近’的含义,而且电网频率波动与电网频率相对稳定时的采样同步偏差基本一致。
本发明不用传统的‘硬件同步’和‘软件同步’,而是采用‘定频变点法’实现同步。即固定采样频率,不需要每个数据窗口的采样点数相同,解决了各个数据窗口的数据点数不同且不是2的整次幂的问题,这种方案即使在电网频率不断波动的情况下,也能实现同步偏差小于万分之一、甚至更小,不存在时延。减少了采样同步偏差就直接减少了谐波测量误差。对于不同点数的DFT的系数利用可编程片上***计算出来存在数据存储器里,供DSP实时计算谐波测量值查表,这样大大节省了运算时间,提高了谐波分析的实时性。
附图说明
图1:定频变点电网谐波测量仪的结构框图。
图2:本发明定频变点电网谐波测量仪的实施结构图。
具体实施方式
下面结合附图对发明的技术方案进行详细说明:
如图1所示,定频变点谐波测量仪包括两大特殊构件:(1)定频整周采样器:以恒定采样频率fS对电网电压(电流)波形采样,由过零检测器控制每个数据窗口非常接近S个整周波;(2)变点DFT计算器:电网频率波动时每个数据窗口的采样点数是可变的,设计了一种‘变点DFT计算器’。
如图2所示,定频整周采样器包括模/数转换器、过零检测器和DSP,变点DFT计算器包括可编程片上***和存储器,定频变点谐波测量仪还包括三个电压互感器和三个电流互感器,其中三个电压互感器和三个电流互感器的输入端分别接电网的三相输出端,三个电压互感器的输出端分别接模/数转换器、过零检测器的输入端,三个电流互感器的输出端分别接模/数转换器的输入端,过零检测器的输出端接DSP的输入端,DSP依次串接可编程片上***、存储器,模/数转换器与DSP双向通信。
本发明检测方法如下:采用可编程片上***建立DFT参数表;电网三相电流分别经电流互感器送到模/数转换器,电网三相电压分别经电压互感器送到模/数转换器和过零检测器;当过零检测器检测到电网三相电压过零点,则DSP通过模/数转换器以恒定采样频率fS采样电网三相电压和三相电流得到采样数据,直至采样宽度大于或等于S个电网整周波宽度并且过零检测器再次检测到电网三相电压过零点,每个所述采样数据窗口的采样点数为Lk、Lk、S为大于零的自然数;DSP根据实时采样数据及采样点数,然后查步骤(1)所述的DFT参数表获得当前测量的谐波参数或则将所述的采样数据经过可编程片上***三角函数型离散富里叶变换DFT,然后查所述的DFT参数表得到谐波参数。
确保每个数据窗口宽度十分逼近S个整周波宽度主要由过零检测器实现:S个电网周波的开始(电网波形过零)采当前数据窗口的第1点,S个电网周波的结束(电网波形过零)停止本数据窗口的采样,这就保证了‘S个整周波’的定时误差小于一个采样周期。
变点DFT的基本特征是:(1)每个数据窗口DFT的点数是变动的且一般不是2的整次幂;(2)DFT的执行用三角函数型DFT(离散富里叶变换),而不是FFT(快速富里叶变换);(3)不同点数的DFT系数用‘查表法’获得;(4)不需要计算的频点从DFT中剔除——大大减少DFT的计算量。
本发明工作步骤和原理是:
用固定的采样频率fS对是电网波形采样,允许采样频率很高——越高越有利。
一个数据窗口的时间宽度规定为S个电网整周波,S个电网周波的开始(电网波形过零)采当前数据窗口的第1点,S个电网周波的结束(电网波形过零)停止本数据窗口的采样,采样同时计算采样点数。按IEC61000-4-30标准的要求S=10。
设当前数据窗口的采样点数为Lk,且Lk个采样数据(采样序列)为:
y=[y(1) y(2)… y(Lk)]
当前数据窗口的采样结束后立即计算Lk点实序列三角函数型DFT:
X=FyT
式中,F为实序列三角函数型DFT的系数矩阵,yT表示向量y的转置:
其中,M是要求的频点数。按IEC61000-4-30标准的要求:需要测量谐波、间谐波共500个频点——M=500——一次DFT需完成500×Lk次‘乘加’运算。
以前由于计算机的数据处理能力不高,DFT曾经是难以实施的,如今高速DSP(数据处理器)的出现已经使得直接用DFT作为谐波分析的工具成为现实。DFT的最大优势在于不要求一个数据窗口的采样点数是2的整次幂,我们则进一步发展为:不要求每个数据窗口的采样点数相同——这就是变点DFT。
为了提高DFT的计算速度,我们采取了两个措施:
(1)不是进行Lk×Lk次‘乘加’运算,而是进行500×Lk次‘乘加’运算,因为高于50次的谐波、间谐波不需要测量。
(2)不同点数的DFT系数(系数矩阵F的元素)用‘查表’法得到。把不同点数的DFT系数预存在大容量flash存储器中,设Lk的可能取值是L0,L0+1,…,L0+n,…,L0+N-1,它们分别对应第1,2,…,n,…,N组DFT系数,根据三角函数的循环性,第n组DFT系数只需存储L0+n个三角函数值——不需要存2M×(L0+n)个三角函数值。
采取了这两个措施的DFT比FFT慢不了多少,用现代DSP足以从容地完成。
考虑电网频率f波动的最大范围为45Hz~55Hz,如果采样频率fS和数据窗口的周波数S都固定不变,那么f=55Hz时的采样点数L0最小:
f=45Hz时的采样点数L0+N-1最大:
其中,round{}表示取整运算。
例如,DSP定时器工作频率为150MHz,取fS=60kHz、S=10,计算出L0=10909,‘采样同步偏差’小于万分之一。下表1是不同电网频率下,本发明所得到的谐波幅值的相对误差(Matlab仿真结果):
表1谐波幅值的相对误差matlab仿真结果表
f | E1 | E2 | E3 | E25 | E31 | E48 | E49 | E50 |
50 | -0.001e-13 | -0.001e-13 | -0.001e-13 | 0.01e-13 | 0.05e-13 | 0.26e-13 | 0.22e-13 | 0.31e-3 |
52 | -0.000038 | -0.000034 | -0.000028 | 0.000861 | 0.001450 | 0.006321 | 0.007439 | 0.009600 |
47 | -0.000073 | -0.000067 | -0.000057 | 0.001324 | 0.002272 | 0.010767 | 0.012816 | 0.016821 |
55 | -0.000008 | -0.000007 | -0.000005 | 0.000208 | 0.000347 | 0.001424 | 0.001662 | 0.002116 |
45 | -0.000024 | -0.000021 | -0.000017 | 0.000581 | 0.000973 | 0.004112 | 0.004819 | 0.006171 |
表中,E1、E2、E3、E25、E31、E48、E49、E50分别代表第1、2、3、25、31、48、49、50次谐波幅值的相对误差,e-13代表10的-13次幂。可见,本方案很容易使谐波精度达标。
总结起来,本发明提出的减少采样同步偏差的方法有如下特点:
1)结构简单、易实现——定频整周采样器可以用过零检测器加ADC构成,变点DFT计算器可以用高速DSP(或PSOC)加flash存储器构成,不需要常见的PLL同步环节;
2)很容易做到采样同步偏差小于万分之一、甚至更小——而传统的‘硬件同步’和‘软件同步’方法很难做到采样同步偏差小于万分之一;
3)电网频率波动与电网频率相对稳定时的采样同步偏差基本一致——而传统的‘硬件同步’和‘软件同步’在电网频率波动时的采样同步偏差比电网频率相对稳定时大得多;
4)同步不存在时延——而‘硬件同步’和‘软件同步’同步都存在时延;
5)极小的采样同步偏差保证了极小的谐波测量误差。
Claims (5)
1.一种定频变点电网谐波检测方法,其特征在于包括如下步骤:
(1)采用可编程片上***建立DFT系数表;
(2)电网三相电流分别经电流互感器送到模/数转换器,电网三相电压分别经电压互感器送到模/数转换器和过零检测器;
(3)S个电网周波的开始(电网波形过零)DSP通过模/数转换器以恒定采样频率fS采当前数据窗口的第1点,S个电网周波的结束(电网波形过零)停止本数据窗口的采样,每个所述采样数据窗口的采样点数为Lk、宽度等于S个电网整周波宽度,Lk、S为大于零的自然数;
(4)DSP根据实时采样数据及采样点数,然后查步骤(1)所述的DFT系数表获得当前测量的谐波参数。
2.一种定频变点电网谐波检测方法,其特征在于包括如下步骤:
(a)采用可编程片上***建立DFT系数表;
(b)电网三相电流分别经电流互感器送到模/数转换器,电网三相电压分别经电压互感器送到模/数转换器和过零检测器;
(c)S个电网周波的开始(电网波形过零)DSP通过模/数转换器以恒定采样频率fS采当前数据窗口的第1点,S个电网周波的结束(电网波形过零)停止本数据窗口的采样,每个所述采样数据窗口的采样点数为Lk、宽度等于S个电网整周波宽度,Lk、S为大于零的自然数;
(d)将步骤(c)所述的采样数据经过可编程片上***三角函数型离散富里叶变换DFT,然后查步骤(a)所述的DFT系数表得到谐波参数。
4.根据权利要求3所述的定频变点电网谐波检测方法,其特征在于所述要求的频点数M按IEC61000-4-30标准的要求:M不大于500。
5.一种基于权利要求1所述的定频变点电网谐波检测方法的测量仪,其特征在于包括:(1)定频整周采样器:以恒定采样频率fS对电网电压或电流波形采样,由过零检测器控制每个数据窗口无限逼近S个整周波;(2)变点DFT计算器:电网频率波动时每个数据窗口的采样点数是可变的,采用高速DSP、执行三角函数型DFT、用查表法获得不同点的DFT系数;所述定频整周采样器包括模/数转换器、过零检测器和DSP,变点DFT计算器包括可编程片上***和存储器,定频变点谐波测量仪还包括三个电压互感器和三个电流互感器,其中三个电压互感器和三个电流互感器的输入端分别接电网的三相输出端,三个电压互感器的输出端分别接模/数转换器、过零检测器的输入端,三个电流互感器的输出端分别接模/数转换器的输入端,过零检测器的输出端接DSP的输入端,DSP依次串接可编程片上***、存储器,模/数转换器与DSP双向通信。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009102327089A CN101701985B (zh) | 2009-11-27 | 2009-11-27 | 定频变点电网谐波检测方法及其测量仪 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009102327089A CN101701985B (zh) | 2009-11-27 | 2009-11-27 | 定频变点电网谐波检测方法及其测量仪 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101701985A true CN101701985A (zh) | 2010-05-05 |
CN101701985B CN101701985B (zh) | 2011-04-20 |
Family
ID=42156905
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009102327089A Expired - Fee Related CN101701985B (zh) | 2009-11-27 | 2009-11-27 | 定频变点电网谐波检测方法及其测量仪 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101701985B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102623968A (zh) * | 2012-03-14 | 2012-08-01 | 上海交通大学 | 基于特征谐波的高压直流输电线路保护方法及其*** |
CN102882211A (zh) * | 2012-10-15 | 2013-01-16 | 北京京仪椿树整流器有限责任公司 | 用于有源电力滤波器的自适应重复控制方法 |
CN103941091A (zh) * | 2014-04-25 | 2014-07-23 | 福州大学 | 基于改进emd端点效应的电力***hht谐波检测方法 |
CN106324342A (zh) * | 2016-08-22 | 2017-01-11 | 江苏南自通华电力自动化股份有限公司 | 一种基于查表的谐波检测方法 |
CN107340454A (zh) * | 2016-04-29 | 2017-11-10 | 中国电力科学研究院 | 一种基于RuLSIF变点探测技术的电力***故障定位分析方法 |
CN109061301A (zh) * | 2018-10-25 | 2018-12-21 | 程海明 | 一种新型电力***间谐波检测装置 |
CN110719043A (zh) * | 2019-11-21 | 2020-01-21 | 深圳国能电气有限公司 | 一种逆变器谐波补偿方法及其控制装置 |
CN113447700A (zh) * | 2021-06-25 | 2021-09-28 | 南京丰道电力科技有限公司 | 一种用于不同频率下对电气量测量的方法及装置 |
-
2009
- 2009-11-27 CN CN2009102327089A patent/CN101701985B/zh not_active Expired - Fee Related
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102623968A (zh) * | 2012-03-14 | 2012-08-01 | 上海交通大学 | 基于特征谐波的高压直流输电线路保护方法及其*** |
CN102623968B (zh) * | 2012-03-14 | 2014-07-02 | 上海交通大学 | 基于特征谐波的高压直流输电线路保护方法及其*** |
CN102882211A (zh) * | 2012-10-15 | 2013-01-16 | 北京京仪椿树整流器有限责任公司 | 用于有源电力滤波器的自适应重复控制方法 |
CN102882211B (zh) * | 2012-10-15 | 2015-09-09 | 北京京仪椿树整流器有限责任公司 | 用于有源电力滤波器的自适应重复控制方法 |
CN103941091A (zh) * | 2014-04-25 | 2014-07-23 | 福州大学 | 基于改进emd端点效应的电力***hht谐波检测方法 |
CN107340454A (zh) * | 2016-04-29 | 2017-11-10 | 中国电力科学研究院 | 一种基于RuLSIF变点探测技术的电力***故障定位分析方法 |
CN107340454B (zh) * | 2016-04-29 | 2020-10-13 | 中国电力科学研究院 | 一种基于RuLSIF变点探测技术的电力***故障定位分析方法 |
CN106324342A (zh) * | 2016-08-22 | 2017-01-11 | 江苏南自通华电力自动化股份有限公司 | 一种基于查表的谐波检测方法 |
CN109061301A (zh) * | 2018-10-25 | 2018-12-21 | 程海明 | 一种新型电力***间谐波检测装置 |
CN110719043A (zh) * | 2019-11-21 | 2020-01-21 | 深圳国能电气有限公司 | 一种逆变器谐波补偿方法及其控制装置 |
CN110719043B (zh) * | 2019-11-21 | 2021-02-05 | 深圳国能电气有限公司 | 一种逆变器谐波补偿方法及其控制装置 |
CN113447700A (zh) * | 2021-06-25 | 2021-09-28 | 南京丰道电力科技有限公司 | 一种用于不同频率下对电气量测量的方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN101701985B (zh) | 2011-04-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101701985A (zh) | 定频变点电网谐波检测方法及其测量仪 | |
CN102435844B (zh) | 一种频率无关的正弦信号相量计算方法 | |
CN102331526B (zh) | 汉宁窗函数连续频谱内插电力谐波参数获取方法 | |
CN102338827B (zh) | 用于电网波形畸变分析和电力谐波参数自动监测方法 | |
CN106932642B (zh) | 电力谐波分析方法 | |
CN103257271A (zh) | 一种基于stm32f107vct6的微电网谐波与间谐波检测装置及检测方法 | |
CN102253282B (zh) | 纳托尔窗函数连续频谱内插电力谐波参数获取方法 | |
CN103454497A (zh) | 基于改进加窗离散傅立叶变换的相位差测量方法 | |
CN104897961A (zh) | 基于互乘法窗函数的三谱线插值fft谐波分析方法及*** | |
CN103575984A (zh) | 基于凯塞窗双谱线插值fft的谐波分析方法 | |
CN103941088A (zh) | 一种基于三相信号的电力***频率快速测量方法 | |
CN106841778A (zh) | 基于pmu实现的次同步和超同步谐波参数的处理方法 | |
CN105004939A (zh) | 一种复合电能质量扰动信号定量分析方法 | |
CN101216511A (zh) | 一种数字同步采样方法 | |
CN104237622A (zh) | 基于软件频率跟踪的采样方法和宽频电压/功率校准装置 | |
CN103018555A (zh) | 一种高精度的电力参数软件同步采样方法 | |
CN101807918A (zh) | 基于同步坐标系的单相锁相环及其实现方法 | |
CN115236392A (zh) | 多特征量的电能计量方法、装置、终端及存储介质 | |
CN102495285B (zh) | 对称窗函数功率重心估计电力谐波参数的方法 | |
CN103969508A (zh) | 一种实时高精密的电力谐波分析方法及装置 | |
CN106569543B (zh) | 一种双通道信号发生器及其输出波形同步方法 | |
CN103543331B (zh) | 一种计算电信号谐波和间谐波的方法 | |
CN104111373A (zh) | 一种用于智能变电站数字化电能的计量方法 | |
CN117452053A (zh) | 一种电力***多回路高速切换采样方法与装置 | |
CN103105532A (zh) | 一种改进傅里叶变换的谐波电能计量***及其控制方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110420 Termination date: 20131127 |