CN100516897C - 一种谐波电能的计量方法 - Google Patents

一种谐波电能的计量方法 Download PDF

Info

Publication number
CN100516897C
CN100516897C CNB2007101199154A CN200710119915A CN100516897C CN 100516897 C CN100516897 C CN 100516897C CN B2007101199154 A CNB2007101199154 A CN B2007101199154A CN 200710119915 A CN200710119915 A CN 200710119915A CN 100516897 C CN100516897 C CN 100516897C
Authority
CN
China
Prior art keywords
harmonic
voltage
current
sub
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.)
Active
Application number
CNB2007101199154A
Other languages
English (en)
Other versions
CN101144841A (zh
Inventor
门长有
王荣华
杨昆
谭年熊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hangzhou hi tech Limited by Share Ltd
Original Assignee
Wangao (hangzhou) Technology Co Ltd
Wangong Sci & Tech Co Ltd Beijing
Hangzhou Vango Science & Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wangao (hangzhou) Technology Co Ltd, Wangong Sci & Tech Co Ltd Beijing, Hangzhou Vango Science & Technology Co Ltd filed Critical Wangao (hangzhou) Technology Co Ltd
Priority to CNB2007101199154A priority Critical patent/CN100516897C/zh
Publication of CN101144841A publication Critical patent/CN101144841A/zh
Application granted granted Critical
Publication of CN100516897C publication Critical patent/CN100516897C/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明提出的谐波电能的计量方法,属于电能计量技术领域,首先对来自电网的电压、电流信号中通过低通滤波对基波进行分离,采用高通滤波对谐波信号进行分离,得到基波电流和基波电压以及总谐波电流和总谐波电压;使用无限长脉冲响应带通滤波方法,从上述总谐波电流和总谐波电压中分离得到需计量的子谐波电流和子谐波电压;根据不同的需计量的谐波频率得到不同次谐波的子谐波电流和子谐波电压;根据基波电流信号和基波电压信号,计算出电网的各种数据。本发明方法的开销小于传统的快速傅立叶变换方法。

Description

一种谐波电能的计量方法
技术领域
本发明涉及一种谐波电能的计量方法,属于电能计量技术领域。
背景技术
随着电网上电力机车、炼钢中的电弧等大功率非线性负载越来越多,电力系统谐波污染日益严重,对电网的安全、经济运行造成了极大的危害。谐波问题涉及谐波源分析、畸变波形分析方法、谐波补偿与抑制、谐波限制标准和谐波下的电能计量等。电能计量是发电企业、输变电企业和用电客户三方结算的依据,在谐波条件下电能计量不能用传统的电子数字式计量所采用的简单功率代数和来计算,须区分用户是谐波用户还是非谐波用户,对谐波用户进行限制和处罚,从而实现用电的科学管理。
目前被学术界和商业界普遍采用的基波、谐波计量和电能质量分析方法是快速傅立叶变换(以下简称FFT)。然而当电网频率不准确等原因造成非同步采样时,FFT存在的频谱泄漏的问题会使得基波成分泄漏到谐波中去,导致基波和谐波计量的不准确。改进的做法是对FFT加窗和插值来解决频谱泄漏的问题,但这两种做法会使计算的复杂性和存储量大大增加。此外,对电网频率漂移造成的频谱泄漏问题无法直接通过上述方法修正,必须借由锁相环(以下简称PLL)技术通过改变采样频率来解决,但由于PLL的设计需要大量模拟电路,会使设计成本显著增加。
FFT的方法应用于基波、谐波计量时总是把所有基波整倍数频率的幅度和相位一并求解再从中提取有用的信息,实际上电网90%以上的能量集中在基波,完全不需要也不应该对基波和谐波的计量付出同样的代价。对于电力企业或用户来说,他们所关心的是基波以及低频的几个奇次谐波的情况。
发明内容
本发明的目的是针对现有技术的不足,提出了一种谐波电能的计量方法,将基波和电网中含量相对较高的低频各奇次谐波分离出来单独计量,以较小的代价计量谐波下基波和各次谐波的有功功率,电流、电压有效值及谐波失真度等数据。
本发明提出的谐波电能的计量方法,包括以下步骤:
(1)对来自电网的电压、电流信号中通过低通滤波对基波进行分离,采用高通滤波对谐波信号进行分离,得到基波电流和基波电压以及总谐波电流和总谐波电压;
(2)使用无限长脉冲响应带通滤波方法,从上述总谐波电流和总谐波电压中分离得到需计量的子谐波电流和子谐波电压,所述的带通滤波的传递函数为:
H BP ( z ) = 1 - a 2 1 - z - 2 1 + b z - 1 + az - 2 ,
式中,a和b为根据需计量的谐波成分设定的系数,0<a<1,b=-cos(2πfe)·(1+a)=-cos(2πfh/fs)·(1+a),fe=fh/fs,π为圆周率,fh为需计量的谐波频率,fe为需计量的谐波归-化频率,fs为对来自电网的电压、电流信号进行采样的频率,z-1,z-2为延迟单元;
(3)根据不同的需计量的谐波频率fh1,fh2,fh3,……,重复步骤(2),得到不同次谐波的子谐波电流和子谐波电压;
(4)当电网的电流、电压的基波频率发生变化时,测量得到电网电流、电压的基波频率变化量Δfh,根据变化量对上述系数b进行修正,得到修正后的系数:
B=b+Δb=-cos(2π(fh+Δfh)/fs)·(1+a),Δb为b的修正量,Δfh为谐波频率的变化量,并根据修正后的系数B,重复上述步骤(2)和(3),计算不同次谐波的子谐波电流和子谐波电压;
(5)根据上述基波电流信号和基波电压信号,计算出电网的基波的功率、基波电流有效值、基波电压有效值、基波能量;
根据上述总谐波电流信号和总谐波电压信号,计算出电网的谐波的总功率、谐波电流总有效值、谐波电压总有效值、谐波总能量;
根据上述各次谐波的子谐波电流信号和各次谐波的子谐波电压信号,分别计算出电网的各次谐波的功率、各次谐波电流有效值、各次谐波电压有效值、各次谐波能量。
本发明提出的谐波电能的计量方法,相对于传统的快速傅立叶变换方法的优点是:本发明克服了传统的快速傅立叶变换方法频谱泄漏的问题;在计算出所有次谐波的情况下,本发明计算的开销小于传统的快速傅立叶变换方法;本发明可以选择仅仅计算关心的谐波成分,近一步降低计算的开销,传统的快速傅立叶变换方法无法实现。
附图说明
图1为本发明方法的原理框图。
图2为本发明方法中低通滤波过程的幅度-频率响应曲线。
图3为本发明方法中高通滤波过程的幅度-频率响应曲线。
图4为本发明方法中带通滤波过程的幅度-频率响应曲线。
具体实施方式
本发明提供了一种谐波电能的计量方法,其原理框图如图1所示,包括:
(1)对来自电网的电压、电流信号中通过低通滤波对基波进行分离,采用高通滤波对谐波信号进行分离,得到基波电流和基波电压以及总谐波电流和总谐波电压;基波电流和基波电压提取的过程中,采用的低通滤波方法在信号频率50赫兹附近的区域都具有平的增益特性,当基波电流和基波电压的频率在50赫兹附近波动时,低通滤波器不会带来增益的误差。
(2)使用无限长脉冲响应带通滤波方法,从上述总谐波电流和总谐波电压中分离得到需计量的子谐波电流和子谐波电压,所述的带通滤波的传递函数为:
H BP ( z ) = 1 - a 2 1 - z 2 1 + bz - 1 + az - 2 ,
式中,a和b为根据需计量的谐波成分设定的系数,0<a<1,b=-cos(2πfe)·(1+a)=-cos(2πfh/fs)·(1+a),fe=fh/fs,π为圆周率,fh为需计量的谐波频率,fe为需计量的谐波归一化频率,fs为对来自电网的电压、电流信号进行采样的频率,z-1,z-2为延迟单元;
(3)根据不同的需计量的谐波频率fh1,fh2,fh3,……,重复步骤(2),得到不同次谐波的子谐波电流和子谐波电压;
(4)为保证电网频率在波动时谐波计量的准确性,防止频谱泄漏,就必需保证带通滤波的中心频率与谐波频率一致。当电网的电流、电压的基波频率发生变化时,测量得到电网电流、电压的基波频率变化量Δfh,根据变化量对上述系数b进行修正,改变滤波器的中心频率,得到修正后的系数:
B=b+Δb=-cos(2π(fh+Δfh)/fs)·(1+a),Δb为b的修正量,Δfh为谐波频率的变化量,并根据修正后的系数B,重复上述步骤(2)和(3),计算不同次谐波的子谐波电流和子谐波电压,采用测量频率修正B系数后,可以消除因为电网频率变化引起的频谱泄漏的问题;
(5)根据上述基波电流信号和基波电压信号,计算出电网的基波的功率、基波电流有效值、基波电压有效值、基波能量;
根据上述总谐波电流信号和总谐波电压信号,计算出电网的谐波的总功率、谐波电流总有效值、谐波电压总有效值、谐波总能量;
根据上述各次谐波的子谐波电流信号和各次谐波的子谐波电压信号,分别计算出电网的各次谐波的功率、各次谐波电流有效值、各次谐波电压有效值、各次谐波能量。下面结合附图详细介绍本发明方法的实施步骤:
图1中输入的电流和电压信号,包含了基波和谐波信号,本发明方法对该信号进行了如下处理:
用低通和高通滤波的方法将电压、电流信号中的基波和谐波信号分离:
(1)低通滤波的目的是滤除各次谐波成分,并且保证47.5Hz-52.5Hz之间有平坦的幅频响应。它可以采用有限长脉冲响应(以下简称FIR)的方法实现,也可以采用IIR的方法实现。
(2)高通滤波的目的是滤除基波成分,使其衰减30分贝(以下简称dB)以上,并且保证100Hz以上有平坦的幅频响应。它可以FIR的方法实现,也可以采用IIR的方法实现。
图2是低通滤波过程的幅度-频率响应曲线,低通滤波是采用的是切比雪夫II型IIR滤波,在50赫兹附近保持平特性。
图3是高通滤波过程的幅度-频率响应曲线,高通滤波是采用的是切比雪夫II型IIR滤波,在100赫兹和100赫兹以上的频带保持平特性。
用IIR带通滤波的方法选择性分离各次谐波成分:
(1)带通滤波的原理:为了节省面积,在实现时可采用2阶椭圆带通滤波结构,其原型传递函数是:
H BP ( z ) = 1 - a 2 1 - z - 2 1 + bz - 1 + az - 2 , 式中,a和b是根据滤波特性而设定的系数,z-1是离散信号的标准延迟单元。
上述的通带转移函数主要特性如下:
对直流OHz的信号衰减无穷大;即M(0)=HBP(1)=0,即z=1,f=0;
对奈奎斯特频率信号的衰减无穷大;即M(0.5)=HBP(-1)=0,即z=-1,f=0.5;
通带的中心频率为: Me = max 0 ≤ f ≤ 0.5 ( M ( f ) ) = | H BP ( e j 2 πfe ) | = 1 , 即z=ej2xfe f = fe = 1 2 π arccos - b 1 + a ; 式中fe为数字归一化频率;
带通滤波的品质因数Q为: Qp = ( 1 + a ) 2 - b 2 2 ( 1 - a ) ; 且Q值越高,带宽越窄,谐波的选择性也就越好;
此外a还需要满足0<a<1;a与b的关系式为:b=-cos(2πfe)□(1+a)。
图4是带通滤波过程的幅度-频率响应曲线,带通滤波过程设计了150赫兹的通带,可以提取三次谐波的能量进行计量。
(2)为保证带通滤波的特性不随电网频率的变化而改变,需要根据频率改变滤波的系数,利用频率测量得到的实际电网频率与理想电网频率的差值作为判据去矫正谐波带通滤波系数。这种机理的原理如下:
当电网的电流、电压的基波频率发生变化时,测量得到电网电流、电压的基波频率变化量Δfh,根据变化量对上述系数b进行修正,得到修正后的系数:
B=b+Δb=-cos(2π(fh+Δfh)/fs)·(1+a),Δb为b的修正量,Δfh为谐波频率的变化量,并根据修正后的系数B校正带同滤波环节,保持滤波器中心频率跟随信号频率的变化而变化。
利用上述得到的基波电压和电流信号以及上述得到的根据频率变化修正带通滤波系数后得到的各次谐波的电压和电流信号计算基波和各次谐波的功率,电压、电流有效值,基波、谐波能量以及脉冲/频率转换输出等谐波下电能计量所需要的数据。

Claims (1)

1、一种谐波电能的计量方法,其特征在于该方法包括以下步骤:
(1)对来自电网的电压、电流信号中,通过低通滤波对基波进行分离,采用高通滤波对谐波信号进行分离,得到基波电流和基波电压以及总谐波电流和总谐波电压;
(2)使用无限长脉冲响应带通滤波方法,从上述总谐波电流和总谐波电压中分离得到需计量的子谐波电流和子谐波电压,所述的带通滤波的传递函数为:
H BP ( z ) = 1 - a 2 1 - z - 2 1 + bz - 1 + az - 2 ,
式中,a和b为根据需计量的谐波成分设定的系数,0<a<1,b=-cos(2πfe)·(1+a)=-cos(2πfh/fs)·(1+a),fe=fh/fs,π为圆周率,fh为需计量的谐波频率,fe为需计量的谐波归一化频率,fs为对来自电网的电压、电流信号进行采样的频率,z-1,z-2为延迟单元;
(3)根据不同的需计量的谐波频率fh1,fh2,fh3,……,重复步骤(2),得到不同次谐波的子谐波电流和子谐波电压;
(4)当电网的电流、电压的基波频率发生变化时,测量得到电网电流、电压的基波频率变化量Δfh,根据变化量对上述系数b进行修正,得到修正后的系数:
B=b+Δb=-cos(2π(fh+Δfh)/fs)·(1+a),Δb为b的修正量,Δfh为谐波频率的变化量,并根据修正后的系数B,重复上述步骤(2)和(3),计算不同次谐波的子谐波电流和子谐波电压;
(5)根据基波电流信号和基波电压信号,计算出电网的基波的功率、基波电流有效值、基波电压有效值、基波能量;
根据总谐波电流信号和总谐波电压信号,计算出电网的谐波的总功率、谐波电流总有效值、谐波电压总有效值、谐波总能量;
根据各次谐波的子谐波电流信号和各次谐波的子谐波电压信号,分别计算出电网的各次谐波的功率、各次谐波电流有效值、各次谐波电压有效值、各次谐波能量。
CNB2007101199154A 2007-08-03 2007-08-03 一种谐波电能的计量方法 Active CN100516897C (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2007101199154A CN100516897C (zh) 2007-08-03 2007-08-03 一种谐波电能的计量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2007101199154A CN100516897C (zh) 2007-08-03 2007-08-03 一种谐波电能的计量方法

Publications (2)

Publication Number Publication Date
CN101144841A CN101144841A (zh) 2008-03-19
CN100516897C true CN100516897C (zh) 2009-07-22

Family

ID=39207470

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2007101199154A Active CN100516897C (zh) 2007-08-03 2007-08-03 一种谐波电能的计量方法

Country Status (1)

Country Link
CN (1) CN100516897C (zh)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101852826B (zh) * 2009-03-30 2012-09-05 西门子公司 一种电力系统的谐波分析方法及其装置
CN102435845A (zh) * 2011-10-19 2012-05-02 广西电网公司电力科学研究院 一种Nuttall自卷积窗加权傅里叶变换的谐波电能计量方法
CN102608423B (zh) * 2012-03-19 2014-05-07 上海电力通信有限公司 一种谐波电能的检测方法
CN102662106B (zh) * 2012-05-30 2014-09-17 重庆市电力公司长寿供电局 谐波电网电能计量方法
CN103105529A (zh) * 2013-01-22 2013-05-15 广西电网公司电力科学研究院 一种基于参数分析的谐波电能计量系统及其控制方法
CN103293378B (zh) * 2013-05-03 2015-08-19 西南交通大学 一种牵引供电系统谐波谐振频率的辨识方法
CN103974240B (zh) * 2014-05-09 2017-09-29 武汉软件工程职业学院 一种基于谐波能量分布的通信设备个体识别方法
CN105911346A (zh) * 2016-07-05 2016-08-31 国网重庆市电力公司市区供电分公司 一种分别计量基波和谐波的双向智能电表及其计量方法
CN106990286B (zh) * 2017-03-14 2019-08-02 合肥工业大学 基于欧拉方法的四象限谐波电能计量装置及方法
CN110514899A (zh) * 2019-08-28 2019-11-29 哈尔滨理工大学 一种畸变信号条件下电能计量方法
CN111257635B (zh) * 2020-01-17 2021-04-30 电子科技大学 一种多端口非线性器件的谐波检测装置
CN114646805B (zh) * 2022-05-23 2022-09-16 四川大学 超高次谐波测量装置及其校正方法
CN115356534A (zh) * 2022-10-21 2022-11-18 成都特隆美储能技术有限公司 一种减少电网电流fft谐波提取时频谱泄露的方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2164053Y (zh) * 1993-03-19 1994-05-04 湖北省电力试验研究所 工频基波电子电能表
WO2001077696A1 (fr) * 2000-04-12 2001-10-18 Central Research Institute Of Electric Power Industry Systeme et procede permettant d'estimer la consommation en energie d'un appareil electrique, et dispositif d'alarme faisant appel auxdits systeme et procede
CN2833607Y (zh) * 2005-08-15 2006-11-01 刘永强 多功能电能计量仪

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2164053Y (zh) * 1993-03-19 1994-05-04 湖北省电力试验研究所 工频基波电子电能表
WO2001077696A1 (fr) * 2000-04-12 2001-10-18 Central Research Institute Of Electric Power Industry Systeme et procede permettant d'estimer la consommation en energie d'un appareil electrique, et dispositif d'alarme faisant appel auxdits systeme et procede
CN2833607Y (zh) * 2005-08-15 2006-11-01 刘永强 多功能电能计量仪

Also Published As

Publication number Publication date
CN101144841A (zh) 2008-03-19

Similar Documents

Publication Publication Date Title
CN100516897C (zh) 一种谐波电能的计量方法
CN102662106B (zh) 谐波电网电能计量方法
CN101915874B (zh) 一种基于傅立叶变换的谐波检测方法
CN103257271B (zh) 一种基于stm32f107vct6的微电网谐波与间谐波检测装置及检测方法
CN103869162B (zh) 一种基于时域准同步的动态信号相量测量方法
CN101261292A (zh) 基于5项Rife-Vincent(I)窗双谱线插值FFT的基波与谐波检测方法
CN103575984A (zh) 基于凯塞窗双谱线插值fft的谐波分析方法
CN104635094A (zh) 一种提升pmu同步相量测量精度的方法
CN203287435U (zh) 一种基于stm32f107vct6的微电网谐波与间谐波检测装置
CN109831184A (zh) 一种带通滤波方法、系统、存储介质及终端设备
CN102338827A (zh) 用于电网波形畸变分析和电力谐波参数自动监测方法
CN109900959A (zh) 一种动态正弦畸变信号中谐波成分的提取方法
CN103399204A (zh) 一种基于Rife-Vincent(II)窗插值FFT的谐波与间谐波检测方法
CN101788604A (zh) 一种基于频域分析和序分量分析的电功率及电量测量方法
CN102998528B (zh) 动态正弦畸变信号中频率分量的提取方法及其实现装置
CN104090159B (zh) 电能计量方法及装置
CN112630527A (zh) 基于经验小波变换的畸变信号电量计量方法
CN108776263A (zh) 基于高阶汉宁自卷积窗及改进插值算法的谐波检测方法
CN101063695B (zh) 无功功率计算电路和方法
CN103412209B (zh) 一种调谐滤波器的失谐度检测方法
CN101464478B (zh) 基于数字同步载波相干解调技术的电压波动闪变测量方法
CN101308175A (zh) 相位谱分析仪
CN105334381A (zh) 一种交流有功功率的测量方法和装置
CN104483524A (zh) 电能表中无功功率的计量系统
CN105242242A (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
ASS Succession or assignment of patent right

Owner name: WANGONG SCI. + TECH. CO., LTD., BEIJING

Free format text: FORMER OWNER: WANGONG SCI. + TECH. CO., LTD., BEIJING WANGAO (HANGZHOU) TECHNOLOGY CO., LTD.

C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20110603

Address after: Hangzhou City, Zhejiang province 310012 Xihu District Huaxing Road No. 99 Neusoft building A408 room

Co-patentee after: Wangong Sci. & Tech. Co., Ltd., Beijing

Patentee after: Hangzhou Vango Science & Technology Co., Ltd.

Address before: Hangzhou City, Zhejiang province 310012 Xihu District Huaxing Road No. 99 Neusoft building A408 room

Co-patentee before: Wangong Sci. & Tech. Co., Ltd., Beijing

Patentee before: Hangzhou Vango Science & Technology Co., Ltd.

Co-patentee before: Wangao (Hangzhou) Technology Co., Ltd.

ASS Succession or assignment of patent right

Free format text: FORMER OWNER: WANGONG SCI. + TECH. CO., LTD., BEIJING

Effective date: 20130410

Owner name: WANGAO (HANGZHOU) TECHNOLOGY CO., LTD.

Free format text: FORMER OWNER: HANGZHOU VANGO SCIENCE + TECHNOLOGY CO., LTD.

Effective date: 20130410

C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20130410

Address after: Hangzhou City, Zhejiang province 310012 Binjiang District Dangfeng Road No. 788, building 901, 902 sea

Patentee after: Wangao (Hangzhou) Technology Co., Ltd.

Address before: Hangzhou City, Zhejiang province 310012 Xihu District Huaxing Road No. 99 Neusoft building A408 room

Patentee before: Hangzhou Vango Science & Technology Co., Ltd.

Patentee before: Wangong Sci. & Tech. Co., Ltd., Beijing

C56 Change in the name or address of the patentee
CP03 Change of name, title or address

Address after: Hangzhou City, Binjiang District province six and 310053 Zhejiang Road No. 368 building B4004 North Sea Base

Patentee after: Hangzhou hi tech Limited by Share Ltd

Address before: Hangzhou City, Zhejiang province 310012 Binjiang District Dangfeng Road No. 788, building 901, 902 sea

Patentee before: Wangao (Hangzhou) Technology Co.Ltd