CN105334381A - 一种交流有功功率的测量方法和装置 - Google Patents
一种交流有功功率的测量方法和装置 Download PDFInfo
- Publication number
- CN105334381A CN105334381A CN201510424513.XA CN201510424513A CN105334381A CN 105334381 A CN105334381 A CN 105334381A CN 201510424513 A CN201510424513 A CN 201510424513A CN 105334381 A CN105334381 A CN 105334381A
- Authority
- CN
- China
- Prior art keywords
- sequence
- frequency
- digital notch
- signal sequence
- unit
- 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
Landscapes
- Measuring Frequencies, Analyzing Spectra (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
Abstract
本发明公开了一种交流有功功率的测量方法及测量装置,所述交流有功功率的测量方法根据初步频率和乘法序列中的混频干扰频率成分设置数字陷波参数,利用该数字陷波参数构造陷波器,利用该陷波器对乘法序列进行数字陷波,得到数字陷波序列。由于该方法的陷波器是根据输入信号序列的初步频率和乘法序列中的混频干扰频率设计的,针对性强,能够对大量复杂的混频频率成分产生深度的抑制作用,因此通过该方法测量的交流有功功率具有准确度高、实时性好的特点。
Description
技术领域
本发明涉及电力信号测量技术领域,特别是一种交流有功功率的测量方法和装置。
背景技术
在电力工程实践中,交流有功功率测量对电力***具有重要意义。传统测量方法主要基于电压信号与电流信号的乘法运算、对乘法信号进行平滑陷波运算,平滑陷波运算信号的终值即为交流有功功率测量结果。当交流电压信号和交流电流信号中存在直流成分、次谐波成分及分次谐波成分时,在两个信号乘法运算后将产生大量复杂的混频频率成分,需要进行平滑滤波处理。典型的平滑滤波运算包括:直接积分运算、在积分运算前加窗口函数运算,典型的窗口函数如汉宁窗函数、布莱克曼窗函数等,具有良好的低通滤波效果,但这些平滑滤波运算没有针对性,不能够对所述大量复杂的混频频率成分产生深度的抑制作用,因此不能满***流有功功率的高准确度测量要求。
发明内容
本发明所要解决的问题在于提供一种准确度高、实时性好的交流有功功率的测量方法和装置。
为解决上述技术问题,本发明采用如下技术方案。
一种交流有功功率的测量方法,包括以下步骤:
根据预设信号时间长度和预设采样频率,对电力信号进行采样,获得输入信号序列,所述电力信号包括电压信号和电流信号,所述输入信号序列分别包括获得电压信号序列和电流信号序列;
测量所述电压信号序列或电流信号序列的频率,作为所述电力信号的初步频率;
将所述电压信号序列和电流信号序列同步相乘,得到乘法序列;
根据所述初步频率和所述乘法序列中的混频干扰频率成分设置数字陷波参数;
利用所述数字陷波参数构造陷波器,利用该陷波器对所述乘法序列进行数字陷波,得到数字陷波序列;所述数字陷波序列终值为所述交流有功功率测量结果。
进一步地,通过零交法、基于滤波的算法、基于小波变换算法、基于神经网络的算法、基于DFT变换的频率算法或基于相位差的频率算法对所述输入信号序列进行频率初测,获取所述初步频率。
进一步地,所述数字陷波采用算术平均陷波算法,将所述乘法序列的若干个连续离散值相加,然后取其算术平均值作为本次数字陷波的输出值。
进一步地,所述数字陷波参数指所述乘法序列的若干个连续离散值相加的长度。
进一步地,所述输入信号序列包括直流信号、1/3次谐波成分、1/2次谐波成分、基波、2次谐波成分、3次谐波成分、4次谐波成分、5次谐波成分。
进一步地,所述乘法序列中的混频干扰频率以基波频率为计算单位,包括分母为6的分次单位混频频率、分母为3的分次单位混频频率、分母为2的分次单位混频频率、次单位混频频率。
进一步地,所述数字陷波过程为:所述乘法序列经过三种陷波参数的陷波器,分别完成对分母为6的分次单位混频频率、分母为3的分次单位混频频率、分母为2的分次单位混频频率和所有次单位混频频率的陷波抑制。
本发明还公开了一种交流有功功率的测量装置,其包括:
采样模块,用于根据预设信号时间长度和预设采样频率,对电力信号进行采样,获得输入信号序列,所述电力信号包括电压信号和电流信号,所述输入信号序列分别包括获得电压信号序列和电流信号序列;
频率初测模块,用于测量所述电压信号序列或电流信号序列的频率,作为所述电力信号的初步频率;
乘法序列模块,用于将所述电压信号序列和电流信号序列同步相乘,得到乘法序列;
陷波模块,用于根据所述初步频率和所述乘法序列中的混频干扰频率成分设置数字陷波参数;利用所述数字陷波参数构造陷波器,利用该陷波器对所述乘法序列进行数字陷波,得到数字陷波序列;所述数字陷波序列终值为所述交流有功功率测量结果。
进一步地,所述陷波模块进行数字陷波时采用算术平均陷波算法,将所述乘法序列的若干个连续离散值相加,然后取其算术平均值作为本次数字陷波的输出值。
进一步地,所述陷波模块设置的数字陷波参数指若干个连续离散值相加的长度。
进一步地,所述采样模块获得的所述输入信号序列包括直流信号、1/3次谐波成分、1/2次谐波成分、基波、2次谐波成分、3次谐波成分、4次谐波成分、5次谐波成分。
进一步地,所述乘法序列模块得到的所述乘法序列中的混频干扰频率以基波频率为计算单位,包括分母为6的分次单位混频频率、分母为3的分次单位混频频率、分母为2的分次单位混频频率、次单位混频频率。
进一步地,所述陷波模块的数字陷波过程为:所述乘法序列经过三种陷波参数的陷波器,分别完成对分母为6的分次单位混频频率、分母为3的分次单位混频频率、分母为2的分次单位混频频率和所有次单位混频频率的陷波抑制。
本发明与现有技术相比的有益效果是:上述交流有功功率的测量方法根据初步频率和乘法序列中的混频干扰频率成分设置数字陷波参数,利用该数字陷波参数构造陷波器,利用该陷波器对乘法序列进行数字陷波,得到数字陷波序列。由于该方法的陷波器是根据输入信号序列的初步频率和乘法序列中的混频干扰频率设计的,针对性强,能够对大量复杂的混频频率成分产生深度的抑制作用,因此通过该方法测量的交流有功功率具有准确度高、实时性好的特点。
附图说明
图1是本发明交流有功功率的测量方法的流程示意图。
图2是一些实施例中三级数字陷波域幅频特性仿真图。
图3是本发明交流有功功率相对误差绝对值|Xerr(f)|随信号基波频率f变化特性的实验结果图。
图4是本发明交流有功功率的测量装置的结构示意图。
具体实施方式
为了使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步地详细描述。
请参阅图1,图1是一些实施方式中交流有功功率的测量方法的流程示意图。
步骤S101,根据预设信号时间长度和预设采样频率,对电力信号进行采样,获得电力信号序列,所述电力信号包括电压信号和电流信号,所述电力信号序列分别包括获得电压信号序列和电流信号序列;
步骤S102,测量所述电压信号序列或电流信号序列的频率,得到所述电力信号的初步频率;
步骤S103,将所述电压信号序列和电流信号序列同步相乘,得到乘法序列;
步骤S104,根据所述初步频率和所述乘法序列中的混频干扰频率成分设置数字陷波参数;
步骤S105,利用所述数字陷波参数构造陷波器,利用该陷波器对所述乘法序列进行数字陷波,得到数字陷波序列;所述数字陷波序列终值为所述交流有功功率测量结果。
上述交流有功功率的测量方法根据初步频率和乘法序列中的混频干扰频率成分设置数字陷波参数,利用该数字陷波参数构造陷波器,利用该陷波器对乘法序列进行数字陷波,得到数字陷波序列。由于该方法的陷波器是根据输入信号序列的初步频率和乘法序列中的混频干扰频率设计的,针对性强,能够对大量复杂的混频频率成分产生深度的抑制作用,因此通过该方法测量的交流有功功率具有准确度高、实时性好的特点。
其中,对于步骤S101,可通过电网领域的采样设备对所述电压信号和电流信号进行采样,分别获得电压信号序列和电流信号序列。
优选地,可根据在额定频率50Hz,采样频率远大于电力***额定频率的原则设置预设数的采样频率。
进一步地,为了保证一定的频率测量实时性,电力信号对应的输入时间可取0.25s。
更进一步地,电力***额定频率50Hz,为了提高性能,采样频率应远大于50Hz,优选地,设置采样频率等于fn=10KHz,采样间隔表达为式(1):
式(1)中,Tn为采样间隔,单位s;fn为所述预设采样频率,单位Hz。
所述样电力信号序列长度表达为式(2):
N=Tsfn(2);
式(2)中,N为电力信号序列长度,单位无量纲;Ts为电力信号对应的输入时间,单位s;fn为所述预设采样频率,单位Hz。在采样频率fn=10KHz,输入信号对应的输入时间Ts=0.25s时,则输入信号序列长度N=2500。
在一些实施例中,在不考虑直流成分和谐波成分时,所述电力信号序列表达为式(3):
XI(n)=Icos(ωTnn+β)
(3);
n=0,1,2,.....,N-1
式(3)中,XU(n)为电压信号序列;U为电压幅值,单位V;为电压信号初相位,单位rad;XI(n)为电流信号序列;I为电流幅值,单位A;β为电流信号初相位,单位rad;ω为信号基波频率,单位rad/s;Tn为采样间隔,单位s;n为序列离散数,单位无量纲;N为输入信号序列长度。
对于步骤S102,可通过零交法对所述电压信号序列或者电流信号序列进行频率测量,以作为所述初步频率。由于电流信号序列的频率较难测准,优选的,测量电压信号序列的频率。另外,还可通过本领域技术人员惯用的其他频率测量方法如基于滤波的算法、基于小波变换算法、基于神经网络的算法、基于DFT变换的频率算法、基于相位差的频率算法等。对所述输入信号序列进行频率初测,允许频率初测存在±0.25%以内相对误差。
所述初步频率表达式为式(4):
ωo(4);
式(4)中,ωo为初步频率,单位rad/s;
对于步骤S103,将所述电压信号序列和电流信号序列同步相乘,得到乘法序列。
在一些实施例中,在不考虑直流成分和谐波成分时,所述乘法序列表达为式(5):
P(n)=XU(n)XI(n)
(5);
n=0,1,2,.....,N-1
式(5)中,XU(n)为电压信号序列;XI(n)为电流信号序列;P(n)为所述乘法序列。
在一些实施例中,在电力信号序列含有直流成分、谐波成分时,设电力信号序列由直流、1/3次谐波、1/2次谐波、基波、2次谐波、3次谐波、4次谐波、5次谐波成分所构成,电压信号序列表达为式(6),电流信号序列表达为式(7):
式(6)中,XU+(n)为电压信号序列;U0为电压信号中的直流幅值,单位V;U1为电压信号中的基波幅值,单位V;U1/3为电压信号中的1/3分次谐波幅值,单位V;U1/2为电压信号中的1/2分次谐波幅值,单位V;U2为电压信号中的2次谐波幅值,单位V;U3为电压信号中的3次谐波幅值,单位V;U4为电压信号中的4次谐波幅值,单位V;U5为电压信号中的5次谐波幅值,单位V;ω为信号基波频率,单位rad/s;Tn为采样间隔,单位s;n为序列离散数,单位无量纲;N为输入信号序列长度。
式(7)中,XI+(n)为电流信号序列;I0为电流信号中的直流幅值,单位A;I1为电流信号中的基波幅值,单位A;I1/3为电流信号中的1/3分次谐波幅值,单位A;I1/2为电流信号中的1/2分次谐波幅值,单位A;I2为电流信号中的2次谐波幅值,单位A;I3为电流信号中的3次谐波幅值,单位A;I4为电流信号中的4次谐波幅值,单位A;I5为电流信号中的5次谐波幅值,单位A;ω为信号基波频率,单位rad/s;Tn为采样间隔,单位s;n为序列离散数,单位无量纲;N为输入信号序列长度。
由式(6)、式(7)可知,电压信号序列与电流信号序列的乘法序列的成分将非常复杂,为了简化问题,将电压信号序列与电流信号序列的乘法序列简化为式(8):
P[n]=P0+P~[n](8);
式(8)中,P0为直流成分,P~[n]为混频干扰成分。其中,P0为直流成分为式(9):
对于步骤S104,根据所述初步频率设置数字陷波参数。
优选地,数字陷波具体采用算术平均陷波算法,即将所述乘法序列的若干个连续离散值相加,然后取其算术平均值作为本次陷波值输出。所述数字陷波参数指平均值陷波序列长度NT,其需要根据对所述输入信号序列进行频率初测得到的初步频率ωo来设置,具体还需要根据混频频率设置该数字陷波参数。
式(7)中的混频干扰成分即为式(6)序列中各频率成分相互混频得到的混频频率,如表1所示:
表1
表1中,a0、a1/3、a1/2、a1、a2、a3、a4、a5分别代表电压信号中的直流、1/3次谐波、1/2次谐波、基波、2次谐波、3次谐波、4次谐波、5次谐波频率成分;b0、b1/3、b1/2、b1、b2、b3、b4、b5分别代表电流信号中的直流、1/3次谐波、1/2次谐波、基波、2次谐波、3次谐波、4次谐波、5次谐波频率成分。
表1中混频频率分成以基波频率为计算单位,可分为:分次单位混频频率,次单位混频频率,具体地可分为5类:零频率成分、分母为6的分次单位混频频率、分母为3的分次单位混频频率、分母为2的分次单位混频频率、次单位混频频率。其中零频率成分为有用成分,即所述的直流成分,其它的成分均为混频干扰频率成分。
NT为采用算术平均陷波算法时若干个连续离散值相加的长度,现设NT1为六分之一初步频率的单位周期序列长度,该序列用于对分母为6的分次单位混频频率进行陷波抑制,陷波参数计算为式(10):
式(10)中,int代表取整数,fn为所述预设采样频率,单位Hz;ωo为初步频率,单位rad/s。在采样频率fn=10KHz,初步频率ωo=100πrad/s时,则NT1=1200。
同理,设NT2为三分之一初步频率的单位周期序列长度时,该序列用于对分母为3的分次单位混频频率进行陷波抑制,陷波参数计算为式(11):
式(11)中,int代表取整数,fn为所述预设采样频率,单位Hz;ωo为初步频率,单位rad/s。在采样频率fn=10KHz,初步频率ωo=100πrad/s时,则NT2=600。
同理,设NT3为二分之一初步频率的单位周期序列长度时,该序列用于对分母为2的分次单位混频频率和所有次单位混频频率进行陷波抑制。陷波参数计算为式(12):
式(12)中,int代表取整数,fn为所述预设采样频率,单位Hz;ωo为初步频率,单位rad/s。在采样频率fn=10KHz,初步频率ωo=100πrad/s时,NT3=400。
对于步骤S105,利用数字陷波参数构造陷波器,利用该陷波器对乘法序列进行数字陷波,得到数字陷波序列;数字陷波序列终值为所述交流有功功率测量结果。
为了提高混频干扰的抑制性能,数字陷波由三种陷波参数的陷波器所构成,三级数字陷波的表达式为(13):
[Xi+(n)]2n=0,1,2,3,....,N-1(13);
Xo(n)n=0,1,2,3,....,N-NT1-NT2-NT3-1
式(13)中,Xo(n)为输出的数字陷波序列、序列长度为N-TT1-TT2-TT3-1。NT1为数字陷波参数1,NT2为数字陷波参数2,NT3为数字陷波参数3。
在一些实施例中,数字陷波参数NT1取值为六分之一参考频率的单位周期序列长度、NT2取值为三分之一参考频率的单位周期序列长度、NT3取值为二分之一参考频率的单位周期序列长度。在三级数字陷波参数NT1=1200、NT2=600、NT3=400时,得到三级数字陷波域幅频特性仿真图,如图2所示。
在频率初测没有误差、三级数字陷波参数没有误差时,混频干扰频率点正好处于图2给出的频域陷波特性的陷波频率点,因此图2给出的频域陷波特性对混频干扰频率成分具有完全的抑制特性。但由于实际存在误差,包括频率初测误差,三级陷波参数整数化误差。在频率初测相对误差±0.25%以内,三级数字陷波参数整数化误差±1时,图2给出的频域陷波特性对混频干扰频率成分仍然具有良好的抑制特性。
由于三级数字陷波滤除了式(8)中波动频率成分,所述数字陷波序列终值等于式(8)中直流成分P0,式(14)为所述交流有功功率测量结果。
式(14)得到交流有功功率测量结果,交流有功功率为直流功率、基波有功功率、所有分次谐波有功功率、所有次谐波有功功率的总和。
为了验证本发明具有较高的准确度,以式(15)所表示的具体信号为例进行验证:
对所述具体信号,在信号时间长度等于0.25s、采样频率10KHz、离散数据量化位数16bit、信号基波频率f变化45Hz—55Hz、频率初测相对误差≤|±0.25%|、得到的交流有功功率相对误差绝对值|Perr(f)|随信号基波频率f变化特性的实验结果图,图3所示。图3给出的交流真有效值测量准确度在10-5量级。
本发明还公开了一种交流真有效值的测量装置,如图4所示,在一些实施例中该测量装置包括采样模块1010、频率初测模块1020、乘法序列模块1030、陷波模块1040,其中:
采样模块1010,用于根据预设信号时间长度和预设采样频率,对电力信号进行采样,获得输入信号序列,所述电力信号包括电压信号和电流信号,所述输入信号序列分别包括获得电压信号序列和电流信号序列;
频率初测模块1020,用于测量所述电压信号序列或电流信号序列的频率,得到所述电力信号的初步频率;
乘法序列模块1030,用于将所述电压信号序列和电流信号序列同步相乘,得到乘法序列;
陷波模块1040,用于根据所述初步频率和所述乘法序列中的混频干扰频率成分设置数字陷波参数;利用所述数字陷波参数构造陷波器,利用该陷波器对所述乘法序列进行数字陷波,得到数字陷波序列;所述数字陷波序列终值为所述交流有功功率测量结果。
在上述交流有功功率的测量装置中,陷波模块1040根据初步频率和乘法序列中的混频干扰频率成分设置数字陷波参数,利用该数字陷波参数构造陷波器,利用该陷波器对乘法序列进行数字陷波,得到数字陷波序列。由于该陷波器是根据输入信号序列的初步频率和乘法序列中的混频干扰频率设计的,针对性强,能够对大量复杂的混频频率成分产生深度的抑制作用,因此由该测量装置测量的交流有功功率具有准确度高、实时性好的特点。
其中,对于采样模块1010,可通过电网领域的采样设备对所述电力信号进行采样,获得输入信号序列。
优选地,可根据在额定频率50Hz,采样频率远大于电力***额定频率的原则设置预设数的采样频率。
进一步地,为了保证一定的频率测量实时性,输入信号对应的输入时间可取0.25s。
更进一步地,电力***额定频率50Hz,为了提高性能,采样频率应远大于50Hz,优选地,设置采样频率等于fn=10KHz,采样间隔表达为式(1):
式(1)中,Tn为采样间隔,单位s;fn为所述预设采样频率,单位Hz。
所述采样输入信号序列长度表达为式(2):
N=Tsfn(2);
式(2)中,N为输入信号序列长度,单位无量纲;Ts为输入信号对应的输入时间,单位s;fn为所述预设采样频率,单位Hz。在采样频率fn=10KHz,输入信号对应的输入时间Ts=0.25s时,则输入信号序列长度N=2500。
在一些实施例中,在不考虑直流成分和谐波成分时,所述电力信号序列表达为式(3):
XI(n)=Icos(ωTnn+β)
(3);
n=0,1,2,.....,N-1
式(3)中,XU(n)为电压信号序列;U为电压幅值,单位V;为电压信号初相位,单位rad;XI(n)为电流信号序列;I为电流幅值,单位A;β为电流信号初相位,单位rad;ω为信号基波频率,单位rad/s;Tn为采样间隔,单位s;n为序列离散数,单位无量纲;N为输入信号序列长度。
对于频率初测模块1020,可通过零交法对所述电压信号序列或者电流信号序列进行频率测量,以作为所述初步频率。由于电流信号序列的频率较难测准,优选的,测量电压信号序列的频率。另外,还可通过本领域技术人员惯用的其他频率测量方法如基于滤波的算法、基于小波变换算法、基于神经网络的算法、基于DFT变换的频率算法、基于相位差的频率算法等。对所述输入信号序列进行频率初测,允许频率初测存在±0.25%以内相对误差。
所述初步频率表达式为式(4):
ωo(4);
式(4)中,ωo为初步频率,单位rad/s;
对于乘法序列模块1030,将所述电压信号序列和电流信号序列同步相乘,得到乘法序列。
在一些实施例中,在不考虑直流成分和谐波成分时,所述乘法序列表达为式(5):
P(n)=XU(n)XI(n)
(5);
n=0,1,2,.....,N-1
式(5)中,XU(n)为电压信号序列;XI(n)为电流信号序列;P(n)为所述乘法序列。
在一些实施例中,在电力信号序列含有直流成分、谐波成分时,设电力信号序列由直流、1/3次谐波、1/2次谐波、基波、2次谐波、3次谐波、4次谐波、5次谐波成分所构成,电压信号序列表达为式(6),电流信号序列表达为式(7):
式(6)中,XU+(n)为电压信号序列;U0为电压信号中的直流幅值,单位V;U1为电压信号中的基波幅值,单位V;U1/3为电压信号中的1/3分次谐波幅值,单位V;U1/2为电压信号中的1/2分次谐波幅值,单位V;U2为电压信号中的2次谐波幅值,单位V;U3为电压信号中的3次谐波幅值,单位V;U4为电压信号中的4次谐波幅值,单位V;U5为电压信号中的5次谐波幅值,单位V;ω为信号基波频率,单位rad/s;Tn为采样间隔,单位s;n为序列离散数,单位无量纲;N为输入信号序列长度。
式(7)中,XI+(n)为电流信号序列;I0为电流信号中的直流幅值,单位A;I1为电流信号中的基波幅值,单位A;I1/3为电流信号中的1/3分次谐波幅值,单位A;I1/2为电流信号中的1/2分次谐波幅值,单位A;I2为电流信号中的2次谐波幅值,单位A;I3为电流信号中的3次谐波幅值,单位A;I4为电流信号中的4次谐波幅值,单位A;I5为电流信号中的5次谐波幅值,单位A;ω为信号基波频率,单位rad/s;Tn为采样间隔,单位s;n为序列离散数,单位无量纲;N为输入信号序列长度。
由式(6)、式(7)可知,电压信号序列与电流信号序列的乘法序列的成分将非常复杂,为了简化问题,将电压信号序列与电流信号序列的乘法序列简化为式(8):
P[n]=P0+P~[n](8);
式(8)中,P0为直流成分,P~[n]为混频干扰成分。其中,P0为直流成分为式(9):
对于陷波模块1040,根据所述初步频率和所述乘法序列中的混频干扰频率成分设置数字陷波参数;利用所述数字陷波参数构造陷波器,利用该陷波器对所述乘法序列进行数字陷波,得到数字陷波序列;所述数字陷波序列终值为所述交流有功功率测量结果。
优选地,数字陷波具体采用算术平均陷波算法,即将若干个连续离散值相加,然后取其算术平均值作为本次陷波值输出。所述数字陷波参数指平均值陷波序列长度NT,其需要根据对所述输入信号序列进行频率初测得到的初步频率ωo来设置,具体还需要根据混频频率设置该数字陷波参数。
乘法序列模块1030的式(8)中的混频干扰成分即式为(6)、式(7)序列中各频率成分相互混频得到的混频频率,如表2所示:
表2
表2中,a0、a1/3、a1/2、a1、a2、a3、a4、a5分别代表电压信号中的直流、1/3次谐波、1/2次谐波、基波、2次谐波、3次谐波、4次谐波、5次谐波频率成分;b0、b1/3、b1/2、b1、b2、b3、b4、b5分别代表电流信号中的直流、1/3次谐波、1/2次谐波、基波、2次谐波、3次谐波、4次谐波、5次谐波频率成分。
表2中混频频率分成以基波频率为计算单位,可分为:分次单位混频频率,次单位混频频率,具体地可分为5类:零频率成分、分母为6的分次单位混频频率、分母为3的分次单位混频频率、分母为2的分次单位混频频率、次单位混频频率。其中零频率成分为有用成分,即所述的直流成分,其它的成分均为混频干扰频率成分。
NT为采用算术平均陷波算法时若干个连续离散值相加的长度,现设NT1为六分之一初步频率的单位周期序列长度,该序列用于对分母为6的分次单位混频频率进行陷波抑制,陷波参数计算为式(10):
式(10)中,int代表取整数,fn为所述预设采样频率,单位Hz;ωo为初步频率,单位rad/s。在采样频率fn=10KHz,初步频率ωo=100πrad/s时,则NT1=1200。
同理,设NT2为三分之一初步频率的单位周期序列长度时,该序列用于对分母为3的分次单位混频频率进行陷波抑制,陷波参数计算为式(11):
式(11)中,int代表取整数,fn为所述预设采样频率,单位Hz;ωo为初步频率,单位rad/s。在采样频率fn=10KHz,初步频率ωo=100πrad/s时,则NT2=600。
同理,设NT3为二分之一初步频率的单位周期序列长度时,该序列用于对分母为2的分次单位混频频率和所有次单位混频频率进行陷波抑制。陷波参数计算为式(12):
式(12)中,int代表取整数,fn为所述预设采样频率,单位Hz;ωo为初步频率,单位rad/s。在采样频率fn=10KHz,初步频率ωo=100πrad/s时,NT3=400。
为了提高混频干扰的抑制性能,数字陷波由三种陷波参数的陷波器所构成,三级数字陷波的表达式为(13):
[Xi+(n)]2n=0,1,2,3,....,N-1(13);
Xo(n)n=0,1,2,3,....,N-NT1-NT2-NT3-1
式(13)中,Xo(n)为输出的数字陷波序列、序列长度为N-TT1-TT2-TT3-1。NT1为数字陷波参数1,NT2为数字陷波参数2,NT3为数字陷波参数3。
在一些实施例中,数字陷波参数NT1取值为六分之一参考频率的单位周期序列长度、NT2取值为三分之一参考频率的单位周期序列长度、NT3取值为二分之一参考频率的单位周期序列长度。在三级数字陷波参数NT1=1200、NT2=600、NT3=400时,得到三级数字陷波域幅频特性仿真图,如图2所示。
由于三级数字陷波滤除了式(8)中波动频率成分,所述数字陷波序列终值等于式(8)中直流成分P0,式(14)为所述交流有功功率测量结果。
式(14)得到交流有功功率测量结果,交流有功功率为直流功率、基波有功功率、所有分次谐波有功功率、所有次谐波有功功率的总和。
以上陈述仅以实施例来进一步说明本发明的技术内容,以便于读者更容易理解,但不代表本发明的实施方式仅限于此,任何依本发明所做的技术延伸或再创造,均受本发明的保护。
Claims (10)
1.一种交流有功功率的测量方法,其特征在于,包括以下步骤:
根据预设信号时间长度和预设采样频率,对电力信号进行采样,获得输入信号序列,所述电力信号分别包括电压信号和电流信号,所述输入信号序列分别包括获得电压信号序列和电流信号序列;
测量所述电压信号序列或电流信号序列的频率,作为所述电力信号的初步频率;
将所述电压信号序列和电流信号序列同步相乘,得到乘法序列;
根据所述初步频率和所述乘法序列中的混频干扰频率成分设置数字陷波参数;
利用所述数字陷波参数构造陷波器,利用该陷波器对所述乘法序列进行数字陷波,得到数字陷波序列;所述数字陷波序列终值为所述交流有功功率测量结果。
2.如权利要求1所述的交流有功功率的测量方法,其特征在于,所述数字陷波采用算术平均陷波算法,将所述乘法序列的若干个连续离散值相加,然后取其算术平均值作为本次数字陷波的输出值。
3.如权利要求1所述的交流有功功率的测量方法,其特征在于,所述数字陷波参数指所述乘法序列的若干个连续离散值相加的长度。
4.如权利要求1所述的交流有功功率的测量方法,其特征在于,所述乘法序列中的混频干扰频率以基波频率为计算单位,包括分母为6的分次单位混频频率、分母为3的分次单位混频频率、分母为2的分次单位混频频率、次单位混频频率。
5.如权利要求1所述的交流有功功率的测量方法,其特征在于,所述乘法序列经过三种陷波参数的陷波器,分别完成对分母为6的分次单位混频频率、分母为3的分次单位混频频率、分母为2的分次单位混频频率和所有次单位混频频率的陷波抑制。
6.一种交流有功功率的测量装置,其特征在于,所述测量装置包括:
采样模块,用于根据预设信号时间长度和预设采样频率,对电力信号进行采样,获得输入信号序列,所述电力信号分别包括电压信号和电流信号,所述输入信号序列分别包括获得电压信号序列和电流信号序列;
频率初测模块,用于测量所述电压信号序列或电流信号序列的频率,作为所述电力信号的初步频率;
乘法序列模块,用于将所述电压信号序列和电流信号序列同步相乘,得到乘法序列;
陷波模块,用于根据所述初步频率和所述乘法序列中的混频干扰频率成分设置数字陷波参数;利用所述数字陷波参数构造陷波器,利用该陷波器对所述乘法序列进行数字陷波,得到数字陷波序列;所述数字陷波序列终值为所述交流有功功率测量结果。
7.如权利要求6所述的交流有功功率的测量装置,其特征在于,所述陷波模块进行数字陷波时采用算术平均陷波算法,将所述乘法序列的若干个连续离散值相加,然后取其算术平均值作为本次数字陷波的输出值。
8.如权利要求6所述的交流有功功率的测量装置,其特征在于,所述陷波模块设置的数字陷波参数指若干个连续离散值相加的长度。
9.如权利要求6所述的交流有功功率的测量装置,其特征在于,所述乘法序列模块得到的所述乘法序列中的混频干扰频率以基波频率为计算单位,包括分母为6的分次单位混频频率、分母为3的分次单位混频频率、分母为2的分次单位混频频率、次单位混频频率。
10.如权利要求6所述的交流有功功率的测量装置,其特征在于,所述陷波模块的数字陷波过程为:所述乘法序列经过三种陷波参数的陷波器,分别完成对分母为6的分次单位混频频率、分母为3的分次单位混频频率、分母为2的分次单位混频频率和所有次单位混频频率的陷波抑制。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510424513.XA CN105334381B (zh) | 2015-07-17 | 2015-07-17 | 一种交流有功功率的测量方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510424513.XA CN105334381B (zh) | 2015-07-17 | 2015-07-17 | 一种交流有功功率的测量方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105334381A true CN105334381A (zh) | 2016-02-17 |
CN105334381B CN105334381B (zh) | 2018-04-13 |
Family
ID=55285027
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510424513.XA Active CN105334381B (zh) | 2015-07-17 | 2015-07-17 | 一种交流有功功率的测量方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105334381B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106771556A (zh) * | 2016-12-23 | 2017-05-31 | 中国计量科学研究院 | 一种基于量子技术的交流功率差分测量***及方法 |
CN109633258A (zh) * | 2019-01-09 | 2019-04-16 | 彭志辉 | 一种有功功率测量电路 |
CN109633257A (zh) * | 2018-12-25 | 2019-04-16 | 彭志辉 | 一种测量准确的有功功率测量方法 |
CN109752584A (zh) * | 2019-01-14 | 2019-05-14 | 彭志辉 | 一种周期信号有效值测量方法 |
CN117388569A (zh) * | 2023-12-11 | 2024-01-12 | 浙江宏仁电气有限公司 | 一种电网波形畸变下的电能计量方法、电能计量箱及介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE3625011A1 (de) * | 1985-11-08 | 1987-05-21 | Aeg Elotherm Gmbh | Einrichtung zur messung der wirkleistung in einem als schwingkreis ausgebildeten lastkreis eines statischen elektrischen umrichters |
US5508617A (en) * | 1993-03-19 | 1996-04-16 | Hewlett-Packard Company | Electric power measuring apparatus and method |
CN1815247A (zh) * | 2006-03-10 | 2006-08-09 | 中国航天科技集团公司第五研究院第五一四研究所 | 一种电功率测量方法及其测量装置 |
CN101627312A (zh) * | 2007-02-27 | 2010-01-13 | 三菱电机株式会社 | 交流电量测量装置 |
CN102187586A (zh) * | 2008-10-23 | 2011-09-14 | Nxp股份有限公司 | 用于自适应射频滤波的***和方法 |
CN102682781A (zh) * | 2012-04-01 | 2012-09-19 | 杭州联汇数字科技有限公司 | 一种音频信号中时变正弦干扰的自适应抑制方法 |
CN103439691A (zh) * | 2013-08-30 | 2013-12-11 | 西安电子科技大学 | 宽带组网雷达抑制窄带干扰的方法 |
-
2015
- 2015-07-17 CN CN201510424513.XA patent/CN105334381B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE3625011A1 (de) * | 1985-11-08 | 1987-05-21 | Aeg Elotherm Gmbh | Einrichtung zur messung der wirkleistung in einem als schwingkreis ausgebildeten lastkreis eines statischen elektrischen umrichters |
US5508617A (en) * | 1993-03-19 | 1996-04-16 | Hewlett-Packard Company | Electric power measuring apparatus and method |
CN1815247A (zh) * | 2006-03-10 | 2006-08-09 | 中国航天科技集团公司第五研究院第五一四研究所 | 一种电功率测量方法及其测量装置 |
CN101627312A (zh) * | 2007-02-27 | 2010-01-13 | 三菱电机株式会社 | 交流电量测量装置 |
CN102187586A (zh) * | 2008-10-23 | 2011-09-14 | Nxp股份有限公司 | 用于自适应射频滤波的***和方法 |
CN102682781A (zh) * | 2012-04-01 | 2012-09-19 | 杭州联汇数字科技有限公司 | 一种音频信号中时变正弦干扰的自适应抑制方法 |
CN103439691A (zh) * | 2013-08-30 | 2013-12-11 | 西安电子科技大学 | 宽带组网雷达抑制窄带干扰的方法 |
Non-Patent Citations (2)
Title |
---|
徐垦 等: "有功功率测量的异步采样方法", 《电子测量技术》 * |
郭磊: "基于自适应陷波器的微电网信息提取及应用", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106771556A (zh) * | 2016-12-23 | 2017-05-31 | 中国计量科学研究院 | 一种基于量子技术的交流功率差分测量***及方法 |
CN106771556B (zh) * | 2016-12-23 | 2019-02-22 | 中国计量科学研究院 | 一种基于量子技术的交流功率差分测量***及方法 |
CN109633257A (zh) * | 2018-12-25 | 2019-04-16 | 彭志辉 | 一种测量准确的有功功率测量方法 |
CN109633257B (zh) * | 2018-12-25 | 2020-10-09 | 彭志辉 | 一种测量准确的有功功率测量方法 |
CN109633258A (zh) * | 2019-01-09 | 2019-04-16 | 彭志辉 | 一种有功功率测量电路 |
CN109633258B (zh) * | 2019-01-09 | 2020-10-02 | 彭志辉 | 一种有功功率测量电路 |
CN109752584A (zh) * | 2019-01-14 | 2019-05-14 | 彭志辉 | 一种周期信号有效值测量方法 |
CN109752584B (zh) * | 2019-01-14 | 2020-10-09 | 彭志辉 | 一种周期信号有效值测量方法 |
CN117388569A (zh) * | 2023-12-11 | 2024-01-12 | 浙江宏仁电气有限公司 | 一种电网波形畸变下的电能计量方法、电能计量箱及介质 |
CN117388569B (zh) * | 2023-12-11 | 2024-03-01 | 浙江宏仁电气有限公司 | 一种电网波形畸变下的电能计量方法、电能计量箱及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN105334381B (zh) | 2018-04-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104635044B (zh) | 基于幅值调制的电力信号频率检测方法和*** | |
CN103869162B (zh) | 一种基于时域准同步的动态信号相量测量方法 | |
CN104635045B (zh) | 基于相位调制的电力信号频率检测方法和*** | |
CN104502700A (zh) | 电力信号的正弦参数测量方法和*** | |
CN105334381A (zh) | 一种交流有功功率的测量方法和装置 | |
CN104502706B (zh) | 电力信号的谐波幅值测量方法和*** | |
CN104535836B (zh) | 电力信号的基波频率测量方法和*** | |
CN103454497B (zh) | 基于改进加窗离散傅立叶变换的相位差测量方法 | |
CN103257271A (zh) | 一种基于stm32f107vct6的微电网谐波与间谐波检测装置及检测方法 | |
CN104459321B (zh) | 电力信号的基波相位测量方法和*** | |
CN203287435U (zh) | 一种基于stm32f107vct6的微电网谐波与间谐波检测装置 | |
CN109831184B (zh) | 一种带通滤波方法、***、存储介质及终端设备 | |
CN104502698A (zh) | 电力信号的频率测量方法和*** | |
CN104502675B (zh) | 电力信号的基波幅值测量方法和*** | |
CN104502701B (zh) | 基于相位调制检测电力信号频率的方法和*** | |
CN110221116A (zh) | 基于加窗插值和解析模式分解的电压闪变包络检测方法 | |
CN105445541A (zh) | 一种任意频率下自适应功率计算方法 | |
CN104111373A (zh) | 一种用于智能变电站数字化电能的计量方法 | |
CN104502702B (zh) | 检测电力信号的频率的方法和*** | |
CN105137186A (zh) | 一种微机自动准同期装置的同期电压相位差测量法 | |
CN103592514A (zh) | 一种新型谐波高精度检测方法 | |
CN104483547A (zh) | 电力信号的滤波方法及*** | |
CN104459318B (zh) | 电力谐波的测量方法和*** | |
CN105353197A (zh) | 一种交流真有效值的测量方法和装置 | |
CN103308759A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |