CN112230054A - 一种谐波含有量和功率计算方法 - Google Patents
一种谐波含有量和功率计算方法 Download PDFInfo
- Publication number
- CN112230054A CN112230054A CN201911335483.XA CN201911335483A CN112230054A CN 112230054 A CN112230054 A CN 112230054A CN 201911335483 A CN201911335483 A CN 201911335483A CN 112230054 A CN112230054 A CN 112230054A
- Authority
- CN
- China
- Prior art keywords
- harmonic
- voltage
- current
- value
- calculating
- 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.)
- Withdrawn
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R21/00—Arrangements for measuring electric power or power factor
- G01R21/06—Arrangements for measuring electric power or power factor by measuring current and voltage
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R21/00—Arrangements for measuring electric power or power factor
- G01R21/133—Arrangements for measuring electric power or power factor by using digital technique
- G01R21/1331—Measuring real or reactive component, measuring apparent energy
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Mathematical Physics (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及一种谐波含有量和功率计算方法,包括如下步骤:S1.采样被测系统的一个周波的电压和电流值,采样点数为正整数M,将电压、电流信号模数转换成数字量,系统的采样率Fs需满足Fs=5000Hz;S2.对电压、电流采样点进行50n Hz频点的DFT运算,得各次谐波的同相分量幅度和正交分量幅度,各频点DFT所需的三角函数计算通过查找表法进行;S3.计算电压、电流各次谐波的谐波含有量;S4.计算各次谐波含有量的相角的正弦值和余弦值;S5.采用和角公式对相角正弦值和余弦值进行修正;S6.利用和角公式计算各次谐波电压电流之间的夹角;S7.计算各次谐波的有功功率、无功功率、视在功率。
Description
技术领域
本发明属于电力系统电能计量领域,具体涉及一种谐波含有量和功率计算方法。
背景技术
众所周知,谐波污染问题对电力系统的安全、稳定、经济运行构成了一定威胁;所以,基波、谐波电参量准确检测和计量是电网电能计量的一项重要技术指标。近几年,新颁布的技术标准也对这一指标提出了更高的要求,如 GB/T17215.302—2013对有功基波、谐波电能计量精度有明确要求。
电力谐波分析通常采用FFT(快速傅里叶变换)算法,当出现非整周波的截断 (即非同步采样)时会发生栅栏效应和频谱泄露,影响谐波分析精度。针对FFT 分析算法的不足,目前普遍采用的解决办法是加窗频域插值法,即通过加窗函数可以消除频谱泄露引起的误差,通过频域插值算法可以消除栅栏效应引起的误差。已有技术的不足之处在于这些操作运算量都很大,就算得以实现,也会要占用绝大部分处理器资源,实时性会大打折扣。所以设计一种实现更简便,运算量更小,而且能确保谐波含有量、失真度、功率计算满足精度要求的算法尤为重要。
发明内容
本发明的目的在于提供一种实现更加简便、运算量更小且计算准确的谐波含有量和功率的计算方法。
本发明的技术方案为:一种谐波含有量和功率计算方法,包括如下步骤:
S1.采样被测系统的一个周波的电压和电流值,采样点数为正整数M,将电压、电流信号模数转换成数字量,系统的采样率Fs需满足Fs=5000Hz,由此的M=100;
S2.对电压、电流采样点进行50n Hz频点的DFT运算,其中 1<n<HARM_NUM,HARM_NUM为所需计量的谐波最高次数,本系统HARM_NUM最高可为50,可得各次谐波的同相分量幅度Re_I(n)、Re_U(n),正交分量幅度Im_I(n)、 Im_U(n),幅度Am_I(n)、Am_U(n),各频点DFT所需的三角函数计算,可通过查找表法进行,可节省计算时间;
S3.计算电压、电流各次谐波的谐波含有量,计算公式为Countent_I(n)= Am_I(n)/Am_I(1),Countent_U(n)=Am_U(n)/Am_U(1);
S4.由步骤2获得的数据可以计算各次谐波含有量的相角的正弦值和余弦值;
S5.采用和角公式对步骤4获得的相角正弦值和余弦值进行修正,获得各次谐波的准确相角正弦值和余弦值,其中各次谐波正余弦值的修正值可通过给系统加相角为0的各次谐波工况来测得;
S6.由步骤5算得的各次谐波的电压和电流值的正余弦值,利用和角公式计算各次谐波电压电流之间的夹角;
S7.由基波电压、电流的幅值,结合计算出的各次谐波含有量的值和相角正余弦值,可计算各次谐波的有功功率、无功功率、视在功率。
本发明的特征在于计算量小,可在成本低廉的低性能芯片上部署,同时,本发明也有一定的局限性,为了节省计算资源,限定了系统的采样率;本发明的方法适用于工况较为稳定不会短时突变的电力系统谐波计量。
附图说明
图1为本发明方法的流程图。
具体实施方式
以下通过具体事例和附图对本发明进行详细的说明。系统采样率5000Hz。如图1所示为部署的算法的执行过程:
S1.采样被测系统的一个周波的电压和电流值,采样点数为正整数M,将电压、电流信号模数转换成数字量,系统的采样率Fs需满足Fs=5000Hz,由此的M=100;
S2.对电压、电流采样点进行50n Hz频点的DFT运算,其中 1<n<HARM_NUM,HARM_NUM为所需计量的谐波最高次数,本系统HARM_NUM最高可为50,可得各次谐波的同相分量幅度Re_I(n)、Re_U(n),正交分量幅度Im_I(n)、 Im_U(n),幅度Am_I(n)、Am_U(n),各频点DFT所需的三角函数计算,可通过查找表法进行,可节省计算时间;
S3.计算电压、电流各次谐波的谐波含有量,计算公式为Countent_I(n)=Am_I(n)/Am_I(1),Countent_U(n)=Am_U(n)/Am_U(1);
S4.由步骤2获得的数据可以计算各次谐波含有量的相角的正弦值和余弦值;
S5.采用和角公式对步骤4获得的相角正弦值和余弦值进行修正,获得各次谐波的准确相角正弦值和余弦值,其中各次谐波正余弦值的修正值可通过给系统加相角为0的各次谐波工况来测得;
S6.由步骤5算得的各次谐波的电压和电流值的正余弦值,利用和角公式计算各次谐波电压电流之间的夹角;
S7.由基波电压、电流的幅值,结合计算出的各次谐波含有量的值和相角正余弦值,可计算各次谐波的有功功率、无功功率、视在功率。
采用该算法的电能表,谐波计量可满足1级精度要求。
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (1)
1.一种谐波含有量和功率计算方法,其特征在于包括以下步骤:
S1.采样被测系统的一个周波的电压和电流值,采样点数为正整数M,将电压、电流信号模数转换成数字量,系统的采样率Fs需满足Fs=5000Hz,由此的M=100;
S2.对电压、电流采样点进行50n Hz频点的DFT运算,其中1<n<HARM_NUM,HARM_NUM为所需计量的谐波最高次数,本系统HARM_NUM最高可为50,可得各次谐波的同相分量幅度Re_I(n)、Re_U(n),正交分量幅度Im_I(n)、Im_U(n),幅度Am_I(n)、Am_U(n),各频点DFT所需的三角函数计算,可通过查找表法进行,可节省计算时间;
S3.计算电压、电流各次谐波的谐波含有量,计算公式为Countent_I(n)=Am_I(n)/Am_I(1),Countent_U(n)=Am_U(n)/Am_U(1);
S4.由步骤2获得的数据可以计算各次谐波含有量的相角的正弦值和余弦值;
S5.采用和角公式对步骤4获得的相角正弦值和余弦值进行修正,获得各次谐波的准确相角正弦值和余弦值,其中各次谐波正余弦值的修正值可通过给系统加相角为0的各次谐波工况来测得;
S6.由步骤5算得的各次谐波的电压和电流值的正余弦值,利用和角公式计算各次谐波电压电流之间的夹角;
S7.由基波电压、电流的幅值,结合计算出的各次谐波含有量的值和相角正余弦值,可计算各次谐波的有功功率、无功功率、视在功率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911335483.XA CN112230054A (zh) | 2019-12-20 | 2019-12-20 | 一种谐波含有量和功率计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911335483.XA CN112230054A (zh) | 2019-12-20 | 2019-12-20 | 一种谐波含有量和功率计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112230054A true CN112230054A (zh) | 2021-01-15 |
Family
ID=74111680
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911335483.XA Withdrawn CN112230054A (zh) | 2019-12-20 | 2019-12-20 | 一种谐波含有量和功率计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112230054A (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100088049A1 (en) * | 2008-10-03 | 2010-04-08 | Analog Devices, Inc. | Apparatus and method for measuring active and reactive powers |
CN102998523A (zh) * | 2012-12-25 | 2013-03-27 | 上海贝岭股份有限公司 | 一种用于电能计量的谐波功率计算方法 |
JP2015219189A (ja) * | 2014-05-20 | 2015-12-07 | 一般財団法人電力中央研究所 | 高調波推定装置、高調波推定方法及び高調波推定プログラム |
CN107390022A (zh) * | 2017-08-23 | 2017-11-24 | 威胜集团有限公司 | 基于离散频谱校正的电能计量方法 |
CN109444515A (zh) * | 2018-09-26 | 2019-03-08 | 徐文涛 | 一种基于sdft算法的无功、不平衡与谐波检测方法 |
CN110470926A (zh) * | 2019-07-15 | 2019-11-19 | 国网浙江省电力有限公司电力科学研究院 | 一种交流电网背景谐波阻抗的实用估算方法 |
-
2019
- 2019-12-20 CN CN201911335483.XA patent/CN112230054A/zh not_active Withdrawn
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100088049A1 (en) * | 2008-10-03 | 2010-04-08 | Analog Devices, Inc. | Apparatus and method for measuring active and reactive powers |
CN102998523A (zh) * | 2012-12-25 | 2013-03-27 | 上海贝岭股份有限公司 | 一种用于电能计量的谐波功率计算方法 |
JP2015219189A (ja) * | 2014-05-20 | 2015-12-07 | 一般財団法人電力中央研究所 | 高調波推定装置、高調波推定方法及び高調波推定プログラム |
CN107390022A (zh) * | 2017-08-23 | 2017-11-24 | 威胜集团有限公司 | 基于离散频谱校正的电能计量方法 |
CN109444515A (zh) * | 2018-09-26 | 2019-03-08 | 徐文涛 | 一种基于sdft算法的无功、不平衡与谐波检测方法 |
CN110470926A (zh) * | 2019-07-15 | 2019-11-19 | 国网浙江省电力有限公司电力科学研究院 | 一种交流电网背景谐波阻抗的实用估算方法 |
Non-Patent Citations (2)
Title |
---|
王茂海,孙元章: "基于DFT的电力系统相量及功率测量新算法", 《电力系统自动化》 * |
袁石良,董杰,徐志强,宋美杰: "基于频率测量值的相量及电气量的DFT修正算法", 《电力系统保护与控制》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103869162B (zh) | 一种基于时域准同步的动态信号相量测量方法 | |
CN104897960B (zh) | 基于加窗四谱线插值fft的谐波快速分析方法及系统 | |
Wen et al. | Hanning self-convolution window and its application to harmonic analysis | |
CN101113995A (zh) | 基于Nuttall窗双峰插值FFT的基波与谐波检测方法 | |
CN106405229B (zh) | 一种基波和谐波电能计量方法 | |
CN110837001B (zh) | 一种电力系统中谐波和间谐波的分析方法与装置 | |
CN107589299B (zh) | 基于多频率相量模型的电力信号同步相量测量方法 | |
CN108490257B (zh) | 一种基于频谱曲线拟合的短时窗间谐波测量方法 | |
CN108896944B (zh) | 一种同步测量装置实验室校准仪及其同步相量测量方法 | |
CN109061345B (zh) | 适用于电力系统的有效值测量方法与系统 | |
CN104502707A (zh) | 一种基于三次样条插值的电力系统同步相量测量方法 | |
CN102890190A (zh) | 基于多小波的随机非线性负荷有功电能计量方法 | |
CN103941090A (zh) | 基于谱线能量插值的谐波测量方法 | |
CN102495285B (zh) | 对称窗函数功率重心估计电力谐波参数的方法 | |
CN110967658A (zh) | 一种基于数字微差法的模拟量输入合并单元校验仪溯源的方法 | |
CN103543331B (zh) | 一种计算电信号谐波和间谐波的方法 | |
CN110763930A (zh) | 基于布莱克曼双峰插值谐波分析法的避雷器阻性电流在线监测系统 | |
CN112798861B (zh) | 电压闪变参数识别方法、装置、计算机设备和存储介质 | |
CN112230054A (zh) | 一种谐波含有量和功率计算方法 | |
CN102759659B (zh) | 一种电力系统中电气信号谐波瞬时值的提取方法 | |
Shao et al. | Power harmonic detection method based on dual HSMW Window FFT/apFFT comprehensive phase difference | |
Jiao et al. | An approach for electrical harmonic analysis based on interpolation DFT | |
CN103983852B (zh) | 电能质量谐波分析仪的谐波分析方法 | |
CN108037351A (zh) | 一种变频电源功率因素测量方法 | |
CN106771584A (zh) | 一种应用于配电网的谐波检测方法 |
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 | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20210115 |
|
WW01 | Invention patent application withdrawn after publication |