CN102135567A - 一种电力***交流采样实时频率跟踪及谐波测量方法 - Google Patents
一种电力***交流采样实时频率跟踪及谐波测量方法 Download PDFInfo
- Publication number
- CN102135567A CN102135567A CN2011100464591A CN201110046459A CN102135567A CN 102135567 A CN102135567 A CN 102135567A CN 2011100464591 A CN2011100464591 A CN 2011100464591A CN 201110046459 A CN201110046459 A CN 201110046459A CN 102135567 A CN102135567 A CN 102135567A
- Authority
- CN
- China
- Prior art keywords
- digital
- sampling
- fast fourier
- frequency
- signal
- 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.)
- Pending
Links
Images
Landscapes
- Measurement Of Current Or Voltage (AREA)
Abstract
本发明公开了一种电力***交流采样实时频率跟踪及谐波测量方法,其特征为:取交流采样输入进入固定采样率模数转换器,固定采样率模数转换器分别输出第一段N点长度的数字信号、第二段N点长度数字信号进入第一快速傅里叶模块;第一快速傅里叶模块输出第一段N点长度的数字信号、第二段N点长度数字信号的最高频率分量相角先后送入鉴相器,鉴相器输出校正频率送入数字信号重构滤波器,数字信号重构滤波器根据校正后频率对固定采样率模数转换器的第一段N点长度的数字输出信号进行数字重构采样,将得到的N点重构采样信号送入第二快速傅里叶模块;第二快速傅里叶模块的最高频率分量相角输出送入鉴相器,模数转换器的第三段N点长度的数字信号进入第一快速傅里叶模块,鉴相器根据精度要求判断后输出校正频率送入数字信号重构滤波器。上述过程达到预定迭代次数M后,第二快速傅里叶模块输出即为交流采样信号的各次谐波分量值,鉴相器输出频率即为交流采样信号的频率。采用本发明方法可以实时跟踪测量电网频率波动,实现对电力***交流工频信号的高精度测量。
Description
技术领域
本发明属于交流采样测量应用技术领域,更具体地涉及电力***工频交流采样测量实时频率跟踪及谐波测量技术。
背景技术
随着电力***的快速发展,电力网容量不断增大,结构日趋复杂,电力***中实时监控、调度的自动化就显得十分重要,而数字采集又是实现自动化的重要环节,尤其是如何准确、快速地采集***中的各个模拟电量,一直是电力工作者关注的热点。根据采样信号的不同,可分为直流采样和交流采样2大类。直流采样是把交流电压、电流信号转化为0~5V的直流电压,这种方法的主要优点是算法简单,便于滤波,但投资较大,维护复杂,无法实现实时信号采集,因而在电力***中的应用受到限制。交流采样是把交流量转化为±5V(或0~5V)的交流电压进行采集,主要优点是实时性好,相位失真小,投资少、便于维护;其缺点是算法复杂,精度难以提高,对A/D转换速度要求较高。随着微机技术的发展,交流采样以其优异的性能价格比,有逐步取代直流采样的趋势。
发明内容
本发明的目的是为了克服电力工频信号交流电量的数据采集对A/D转换速度要求较高,精度难以提高等传统交流采样技术的缺点,提供一种快速实时的频率跟踪和高精度交流电量测量方法。
本发明的一种电力***交流采样实时频率跟踪及谐波测量方法是通过如下技术方案实现的,所述的实时同步交流采样方法,包括步骤如下:
1、取交流采样输入1进入固定采样率模数转换器2。
2、固定采样率模数转换器2分别输出第一段N点长度的数字信号、第二段N点长度数字信号进入第一快速傅里叶模块(FFT)3,第一快速傅里叶模块(FFT)3输出第一段N点长度的数字信号、第二段N点长度数字信号的最高频率分量相角8先后送入鉴相器4,鉴相器4输出校正频率10送入数字信号重构滤波器5。
3、数字信号重构滤波器5根据校正后频率10对固定采样率模数转换器2的第一段N点长度的数字输出信号12进行数字重构采样,将得到的N点重构采样信号送入第二快速傅里叶模块6。
4、第二快速傅里叶模块6的最高频率分量相角输出9送入鉴相器4,模数转换器2第三段N点长度的数字信号进入第一快速傅里叶模块3,鉴相器4根据精度要求判断后输出校正频率10送入数字信号重构滤波器5。
5、重复步骤2、3和4,达到预定迭代次数M后,第二快速傅里叶模块6输出7即为交流采样信号的各次谐波分量值,鉴相器4输出频率10即为交流采样信号的频率。
附图说明
图1是本发明结构示意图。
图2是本发明的一个含谐波噪声交流采样输入信号的测量实例。
具体实施方式
图1是本发明方法结构示意图。稳态电网的交流采样输入1表示为x(t),其工频频率表示为f0。由于谐波干扰和其他因素干扰所造成的频率偏差使得电网信号的频率出现波动,模数转换器1以固定采样率f0的N倍速率对信号x(t)的一个周期信号进行采样,数字长度是N点。采样信号表示为xi(nTs),其中n=0,1,2,L N-1为整数,Ts为采样时间间隔。
对模拟电网信号x(t)进行采样,在时间轴上分别取两段长度都为N点的数字得到xm(n)和xn(n),n=0,1,2,L N-1。然后将xm(n)和xn(n)送入第一快速傅里叶变换(FFT)模块3,分别输出为Xm(k)和Xn(k),k=0,1,2,L,N-1
m时刻最高频率分量表示为Xm(k),k表示为输入信号的频率分量。从快速傅里叶变换模块(FFT)输出中得到m时刻N点输入信号的高频率分量的初始相位为:
式中,Im和Re分别表示复数的实部和虚部。在m时刻之后的n时刻,N点数字的快速傅里叶变换模块(FFT分量)最高频率分量表示为Xn(k),其初始相角为:
利用这两个相位角可得到跟踪频率偏差
由(1-3)式得到的最小频率偏差即频率分辨率远高于传统快速傅里叶变换FFT得到的频率分辨率精度,同时具有快速的收敛速度。固定采样频率f0与精细频率fΔ相加得到电网信号的新的频率:
fnew=f0+fΔ (1-4)
这样就可以得到一个相对逼近电网工频信号频率的估算频率值,对电网信号x(t)随后两个周期信号按新采样率fs采样,其中
fs=fnew·N (1-5)
从重构滤波器5的输出信号13中取两段连续N点数字信号,送入模块6重新做快速傅里叶变换FFT,利用(1-1)至(1-5)过程得到校正后的新的采样速率,重复迭代上述过程就可以得到无限逼近真实电网信号的频率,从而实现交流采样输入信号1的高精度实时频率跟踪和测量。
鉴相器2估算得到新的采样率fnew后,对进入模数转换器2之前的交流采样输入1按新采样率fnew进行重新采样,该过程是通过重构滤波器5实现的。重构滤波器5是一种数字信号重构滤波器,对模数转换器2的输出采样信号12进行数字恢复重构。
图2是采用本发明方法的一个含谐波噪声交流采样输入信号的测量实例。稳态电网的交流采样输入1为:
其中f0取52,N取128,按照本发明方法实施如下步骤:
1、取交流采样输入1进入固定采样率模数转换器2。
2、固定采样率模数转换器2分别输出第一段N点长度的数字信号、第二段N点长度数字信号进入第一快速傅里叶模块(FFT)3,第一快速傅里叶模块(FFT)3输出第一段N点长度的数字信号、第二段N点长度数字信号的最高频率分量相角8先后送入鉴相器4,鉴相器4输出校正频率10送入数字信号重构滤波器5。
3、数字信号重构滤波器5根据校正后频率10对固定采样率模数转换器2的第一段N点长度的数字输出信号12进行数字重构采样,将得到的N点重构采样信号送入第二快速傅里叶模块6。
4、第二快速傅里叶模块6的最高频率分量相角输出9送入鉴相器4,模数转换器2第三段N点长度的数字信号进入第一快速傅里叶模块3,鉴相器4根据精度要求判断后输出校正频率10送入数字信号重构滤波器5。
5、重复步骤2、3和4,经过3次迭代后频率跟踪和测量误差达到2.7e-12,各次谐波测量最大误差1e-15。
Claims (4)
1.一种电力***交流采样实时频率跟踪及谐波测量方法,其特征在于:该方法包括步骤如下:
(1)取交流采样输入1进入固定采样率模数转换器2。
(2)固定采样率模数转换器2分别输出第一段N点长度的数字信号、第二段N点长度数字信号进入第一快速傅里叶模块(FFT)3,第一快速傅里叶模块(FFT)3输出第一段N点长度的数字信号、第二段N点长度数字信号的最高频率分量相角8先后送入鉴相器4,鉴相器4输出校正频率10送入数字信号重构滤波器5。
(3)数字信号重构滤波器5根据校正后频率10对固定采样率模数转换器2的第一段N点长度的数字输出信号12进行数字重构采样,将得到的N点重构采样信号送入第二快速傅里叶模块6。
(4)第二快速傅里叶模块6的最高频率分量相角输出9送入鉴相器4,模数转换器2第三段N点长度的数字信号进入第一快速傅里叶模块3,鉴相器4根据精度要求判断后输出校正频率10送入数字信号重构滤波器5。
(5)重复步骤2、3和4,达到预定迭代次数M后,第二快速傅里叶模块6输出7即为交流采样信号的各次谐波分量值,鉴相器4输出频率10即为交流采样信号的频率。
2.根据权利要求1所述的方法,其特征在于:预定迭代次数M为正整数,一般取M值为3至6比较合适。
3.根据权利要求1所述的方法,其特征在于:第一快速傅里叶模块(FFT)3和第二快速傅里叶模块6的最高频率分量相角送入鉴相器4,利用二者的相角之差获得跟踪频率偏差。
4.根据权利要求1所述的方法,其特征在于:数字信号重构滤波器5以校正频率10为采样速率,对固定采样率模数转换器2输出的数字信号进行变速率采样。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011100464591A CN102135567A (zh) | 2011-02-25 | 2011-02-25 | 一种电力***交流采样实时频率跟踪及谐波测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011100464591A CN102135567A (zh) | 2011-02-25 | 2011-02-25 | 一种电力***交流采样实时频率跟踪及谐波测量方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102135567A true CN102135567A (zh) | 2011-07-27 |
Family
ID=44295409
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011100464591A Pending CN102135567A (zh) | 2011-02-25 | 2011-02-25 | 一种电力***交流采样实时频率跟踪及谐波测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102135567A (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103063913A (zh) * | 2012-12-07 | 2013-04-24 | 深圳市金宏威技术股份有限公司 | 用于傅里叶变换的频率跟踪方法 |
CN103728523A (zh) * | 2014-01-21 | 2014-04-16 | 国家电网公司 | 一种城市电网合格率的测量方法 |
CN103941086A (zh) * | 2014-03-25 | 2014-07-23 | 中国科学院长春光学精密机械与物理研究所 | 一种超高精度频率测量仪及其测量方法 |
CN105004920A (zh) * | 2015-07-10 | 2015-10-28 | 国网天津市电力公司 | 傅里叶修正系数频率测量方法 |
CN105004913A (zh) * | 2015-07-21 | 2015-10-28 | 黎文安 | 一种电力***交流采样频率实时跟踪方法 |
CN105044459A (zh) * | 2015-07-21 | 2015-11-11 | 青岛艾诺智能仪器有限公司 | 一种谐波分析方法 |
CN106597095A (zh) * | 2016-12-22 | 2017-04-26 | 南京因泰莱电器股份有限公司 | 一种双cpu等间隔采样的频率跟踪实现方法 |
EP3287795A1 (fr) * | 2016-08-24 | 2018-02-28 | Schneider Electric Industries SAS | Determination de la frequence d'un signal alternatif |
CN108037357A (zh) * | 2017-11-29 | 2018-05-15 | 西电通用电气自动化有限公司 | 一种电力***自动化装置的频率跟踪方法 |
CN109030957A (zh) * | 2015-05-19 | 2018-12-18 | 江苏理工学院 | 介质损耗测量方法 |
CN109030942A (zh) * | 2015-05-19 | 2018-12-18 | 江苏理工学院 | 谐相角分析方法 |
CN110187176A (zh) * | 2019-05-31 | 2019-08-30 | 中国电力科学研究院有限公司 | 一种电网谐波测量装置及方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101221201A (zh) * | 2008-01-25 | 2008-07-16 | 杭州万工科技有限公司 | 一种基波频率的测量方法 |
US20100069030A1 (en) * | 2006-02-10 | 2010-03-18 | Thierry Briand | Frequency measuring broadband digital receiver |
CN101852826A (zh) * | 2009-03-30 | 2010-10-06 | 西门子公司 | 一种电力***的谐波分析方法及其装置 |
-
2011
- 2011-02-25 CN CN2011100464591A patent/CN102135567A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100069030A1 (en) * | 2006-02-10 | 2010-03-18 | Thierry Briand | Frequency measuring broadband digital receiver |
CN101221201A (zh) * | 2008-01-25 | 2008-07-16 | 杭州万工科技有限公司 | 一种基波频率的测量方法 |
CN101852826A (zh) * | 2009-03-30 | 2010-10-06 | 西门子公司 | 一种电力***的谐波分析方法及其装置 |
Non-Patent Citations (3)
Title |
---|
李振然等: "利用富里叶变换实现电力***频率的测量", 《广西大学学报》 * |
牟龙华等: "基于傅里叶变换的精确频率测量算法", 《电力***自动化》 * |
王柏林等: "电网谐波分析的频率估计与跟踪", 《电测与仪表》 * |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103063913A (zh) * | 2012-12-07 | 2013-04-24 | 深圳市金宏威技术股份有限公司 | 用于傅里叶变换的频率跟踪方法 |
CN103063913B (zh) * | 2012-12-07 | 2016-01-20 | 深圳市金宏威技术有限责任公司 | 用于傅里叶变换的频率跟踪方法 |
CN103728523A (zh) * | 2014-01-21 | 2014-04-16 | 国家电网公司 | 一种城市电网合格率的测量方法 |
CN103941086A (zh) * | 2014-03-25 | 2014-07-23 | 中国科学院长春光学精密机械与物理研究所 | 一种超高精度频率测量仪及其测量方法 |
CN109030957A (zh) * | 2015-05-19 | 2018-12-18 | 江苏理工学院 | 介质损耗测量方法 |
CN109030957B (zh) * | 2015-05-19 | 2020-12-25 | 江苏理工学院 | 介质损耗测量方法 |
CN109030942B (zh) * | 2015-05-19 | 2020-12-04 | 江苏理工学院 | 谐相角分析方法 |
CN109030942A (zh) * | 2015-05-19 | 2018-12-18 | 江苏理工学院 | 谐相角分析方法 |
CN105004920A (zh) * | 2015-07-10 | 2015-10-28 | 国网天津市电力公司 | 傅里叶修正系数频率测量方法 |
CN105004920B (zh) * | 2015-07-10 | 2017-11-17 | 国网天津市电力公司 | 傅里叶修正系数频率测量方法 |
CN105044459A (zh) * | 2015-07-21 | 2015-11-11 | 青岛艾诺智能仪器有限公司 | 一种谐波分析方法 |
CN105044459B (zh) * | 2015-07-21 | 2017-09-29 | 青岛艾诺智能仪器有限公司 | 一种谐波分析方法 |
CN105004913A (zh) * | 2015-07-21 | 2015-10-28 | 黎文安 | 一种电力***交流采样频率实时跟踪方法 |
FR3055417A1 (fr) * | 2016-08-24 | 2018-03-02 | Schneider Electric Industries Sas | Determination de la frequence d'un signal alternatif |
CN107782966A (zh) * | 2016-08-24 | 2018-03-09 | 施耐德电器工业公司 | 确定交流信号的频率 |
EP3287795A1 (fr) * | 2016-08-24 | 2018-02-28 | Schneider Electric Industries SAS | Determination de la frequence d'un signal alternatif |
US10613127B2 (en) | 2016-08-24 | 2020-04-07 | Schneider Electric Industries Sas | Determining the frequency of an alternating signal |
CN106597095B (zh) * | 2016-12-22 | 2019-09-10 | 南京因泰莱电器股份有限公司 | 一种双cpu等间隔采样的频率跟踪实现方法 |
CN106597095A (zh) * | 2016-12-22 | 2017-04-26 | 南京因泰莱电器股份有限公司 | 一种双cpu等间隔采样的频率跟踪实现方法 |
CN108037357A (zh) * | 2017-11-29 | 2018-05-15 | 西电通用电气自动化有限公司 | 一种电力***自动化装置的频率跟踪方法 |
CN110187176A (zh) * | 2019-05-31 | 2019-08-30 | 中国电力科学研究院有限公司 | 一种电网谐波测量装置及方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102135567A (zh) | 一种电力***交流采样实时频率跟踪及谐波测量方法 | |
de la O Serna et al. | Instantaneous Oscillating Phasor Estimates With Taylor $^ K $-Kalman Filters | |
CN103869162B (zh) | 一种基于时域准同步的动态信号相量测量方法 | |
CN103257271A (zh) | 一种基于stm32f107vct6的微电网谐波与间谐波检测装置及检测方法 | |
CN203287435U (zh) | 一种基于stm32f107vct6的微电网谐波与间谐波检测装置 | |
CN103543335B (zh) | 一种同步相量测量方法 | |
CN101915874A (zh) | 一种基于傅立叶变换的谐波检测方法 | |
CN103941088A (zh) | 一种基于三相信号的电力***频率快速测量方法 | |
CN203909138U (zh) | 一种通用三相交流信号同步采样装置 | |
CN107478896A (zh) | 一种基于级联广义积分器的频率自适应谐波电流检测方法 | |
CN109142863B (zh) | 一种电力***测频方法及*** | |
CN102236048A (zh) | 一种电力***相量频率测量方法 | |
CN102035554A (zh) | 电能计量和分析***、方法和模数转换电路 | |
CN103018555A (zh) | 一种高精度的电力参数软件同步采样方法 | |
CN102243272A (zh) | 一种高精度采样数据同步的谐波分析方法 | |
CN105629060A (zh) | 基于最优基带滤波的电网频率测量方法和装置 | |
CN104330623A (zh) | 电力***中正弦波信号的参数测量方法及*** | |
CN104155517A (zh) | 一种数字化电能表非整周期采样误差补偿方法及*** | |
CN102495250A (zh) | 一种基于Hilbert算法的准同步宽频无功电能表及其采样方法 | |
CN102809687A (zh) | 一种交流电频率的数字化测量方法 | |
CN103969508A (zh) | 一种实时高精密的电力谐波分析方法及装置 | |
CN114966194A (zh) | 一种基于adc的低频计量三相电能表 | |
CN105044459A (zh) | 一种谐波分析方法 | |
CN105004913A (zh) | 一种电力***交流采样频率实时跟踪方法 | |
CN105182072A (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 | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20110727 |