CN103901273B - 电网谐波检测方法和电网谐波检测装置 - Google Patents

电网谐波检测方法和电网谐波检测装置 Download PDF

Info

Publication number
CN103901273B
CN103901273B CN201210586346.5A CN201210586346A CN103901273B CN 103901273 B CN103901273 B CN 103901273B CN 201210586346 A CN201210586346 A CN 201210586346A CN 103901273 B CN103901273 B CN 103901273B
Authority
CN
China
Prior art keywords
harmonic
vectorial
function component
mode function
intrinsic mode
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.)
Expired - Fee Related
Application number
CN201210586346.5A
Other languages
English (en)
Other versions
CN103901273A (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.)
Beijing Zhizhong Energy Technology Development Co ltd
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CN201210586346.5A priority Critical patent/CN103901273B/zh
Publication of CN103901273A publication Critical patent/CN103901273A/zh
Application granted granted Critical
Publication of CN103901273B publication Critical patent/CN103901273B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明提供了一种电网谐波检测方法和电网谐波检测装置。该电网谐波检测方法,包括:步骤S11,获取待测电网中的交流电量模拟信号,并将模拟信号进行离散化处理,得到离散数据x(t);步骤S13,在离散数据x(t)的基础上叠加多组不同的高斯白噪声信号ωi(t),得到多组原始信号Xi(t),并对多组原始信号Xi(t)分别进行经验模态分解,得到各阶固有模态函数分量;步骤S15,将多组原始信号Xi(t)分解得到的固有模态函数分量进行均值运算,得到第n次固有模态函数分量cn(t);步骤S17,将第n次固有模态函数分量cn(t)代入一组原始信号X(t),得到该组X(t)的谐波分解结果。本发明计算方法简单快捷,谐波检测精度高。

Description

电网谐波检测方法和电网谐波检测装置
技术领域
本发明涉及电力领域,具体而言,涉及一种电网谐波检测方法和电网谐波检测装置。
背景技术
近年来,由于非线性负载(如整流器、变频器等)的大量使用,用户侧的电能质量问题日益严重,其中以谐波污染最为严重。随着智能电网建设的推进,各种分布式电源(如小型光伏发电、光热发电、小型风力发电)、储能设备(如锂电池组、铅酸蓄电池组等)和电动汽车充电站、充电桩开始出现,这些设备并网运行,也给电网带来一定的谐波影响。
谐波污染给电网带来的各种影响具体表现在:降低供配电系统的效率,增加各级电网的损耗;给各种用电设备带来附加损耗,缩短寿命;影响各种表计的计量精度。另外,谐波还可能使各级继电保护装置误动和拒动,给电网带来严重后果。
为提高电能质量,必须首先准确检测出各种情况下的谐波分量,才有可能补偿这些谐波分量。因此,如何快速方便的检测出电网中谐波分量,一直是个重要的研究问题。但是由于电网谐波受到随机性、分布性和非平稳性等因素的影响,难以对其进行同步精确测量。
针对现有技术中的电网谐波检测方法中检测精度较差的问题,目前尚未提出有效的解决方案。
发明内容
本发明旨在提供一种电网谐波检测方法和电网谐波检测装置,以解决电网谐波检测方法中检测精度较差的问题。
为了实现上述目的,根据本发明的一个方面,提供了一种电网谐波检测方法。该电网谐波检测方法,包括:
步骤S11,获取待测电网中的交流电量模拟信号,并将模拟信号进行离散化处理,得到离散数据x(t);
步骤S13,在离散数据x(t)的基础上叠加多组不同的高斯白噪声信号ωi(t),得到多组原始信号Xi(t),Xi(t)=x(t)+ωi(t),并对多组原始信号Xi(t)分别进行经验模态分解,得到各阶固有模态函数分量
步骤S15,将多组原始信号Xi(t)分解得到的固有模态函数分量进行均值运算,得到第n次固有模态函数分量cn(t),表示为:
步骤S17,将第n次固有模态函数分量cn(t)代入一组原始信号X(t),则该组X(t)的谐波分解结果为:其中,X(t)=x(t)+ω(t),rm(t)为高斯白噪声信号ω(t)的各阶固有模态函数分量。
进一步地,步骤S11包括:初始化维数为2N的向量I2N;获取连续的2N个交流电量的采样值,这2N个采样值记为ik-(2N-1),L,ik-1,ik,将ik-(N-1),L,ik-1;将ik-(N-1),L,ik-1赋值给向量I2N的前N-1个元素,将ik-(2N-1),L,ik-N赋值给向量I2N的后N个元素,向量I2N表示为:[ik-(N-1) Lik-1 ik ik-2N+1 L ik-N];将赋值后的向量I2N作为离散数据x(t)。
进一步地,将赋值后的向量I2N作为离散数据x(t)之后还包括:对赋值后的向量I2N依次送入经验模态分解模块、均值运算模块、结果输出模块进行计算,得到谐波分解结果;从得到的谐波分解结果中查询到工频50Hz的固有模态函数分量cf,cf中第N个元素即为k时刻的基波分量ikf;利用公式ikh=ik-ikf,计算得到k时刻的总谐波分量ikh
进一步地,得到谐波分解结果之后还可以包括:从得到的谐波分解结果中查询到特定次谐波的固有模态函数分量cfn,cfn中第N个元素即为k时刻的基波分量ikfn
进一步地,获取连续的2N个交流电量的采样值包括:实时采集交流电量,并记录交流电量的采样值;将当前采集到的采样值记录为ik,从获取当前采样时刻之前连续2N-1个采样时刻的采样数据。
进一步地,获取待测电网中的交流电量模拟信号包括:获取一下任意一种电网电路中的交流电量模拟信号:单相电路、三相三线制电路和三相四线制电路中,其中交流电量包括电流信号和/或电压信号。
基于本发明的另一个方面,提供了一种电网谐波检测装置。该电网谐波检测装置包括:数据离散模块,用于获取待测电网中的交流电量模拟信号,并将模拟信号进行离散化处理,得到离散数据x(t);经验模态分解模块,用于在离散数据x(t)的基础上叠加多组不同的高斯白噪声信号ωi(t),得到多组原始信号Xi(t),Xi(t)=x(t)+ωi(t),并对多组原始信号Xi(t)分别进行经验模态分解,得到各阶固有模态函数分量均值运算模块,用于将多组原始信号Xi(t)分解得到的固有模态函数分量进行均值运算,得到第n次固有模态函数分量cn(t),其公式表示为:结果输出模块,用于将第n次固有模态函数分量cn(t)代入一组原始信号X(t),则该组X(t)的谐波分解结果为:其中,X(t)=x(t)+ω(t),rm(t)为高斯白噪声信号ω(t)的各阶固有模态函数分量。
进一步地,数据离散模块包括:向量初始子模块,用于初始化维数为2N的向量I2N;数据采样子模块,用于获取连续的2N个交流电量的采样值,这2N个采样值记为ik-(2N-1),L,ik-1,ik,将ik-(N-1),L,ik-1;向量构建子模块,用于将ik-(N-1),L,ik-1赋值给向量I2N的前N-1个元素,将ik-(2N-1),L,ik-N赋值给向量I2N的后N个元素,向量I2N表示为:[ik-(N-1) L ik-1 ik ik-2N+1L ik-N];将赋值后的向量I2N作为离散数据x(t)。
进一步地,本发明提供的电网谐波检测装置还包括:总谐波分量计算模块,用于对赋值后的向量I2N执行步骤S13至步骤S17,得到谐波分解结果;从得到的谐波分解结果中查询到工频50Hz的固有模态函数分量cf,cf中第N个元素即为k时刻的基波分量ikf;利用公式ikh=ik-ikf,计算得到k时刻的总谐波分量ikh
进一步地,本发明提供的电网谐波检测装置还包括:特定次谐波分析模块,用于从得到的谐波分解结果中查询到特定次谐波的固有模态函数分量cfn,cfn中第N个元素即为k时刻的基波分量ikfn
应用本发明的技术方案,基于经验模态分解(Empirical Mode Decomposition,EMD)方法的不足,提出了利用高斯白噪声辅助进行数据分析的方法。利用了高斯白噪声频谱的零均值原理,通过多组信号的叠加均值计算消除高斯白噪声作为时域分布参考结构带来的影响,对原始电流或电压信号的采样要求低,不需要电网同步角度,可用于单相、三相三线制或三相四线制系统中检测谐波分量。计算方法简单快捷,谐波检测精度高。可以应用于各种谐波检测仪器、谐波检测软件程序、有源电力滤波器等不同的场合。
附图说明
构成本申请的一部分的说明书附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1是根据本发明实施例的电网谐波检测装置的示意图;
图2是根据本发明实施例的电网谐波检测方法的示意图;
图3A是根据本发明实施例的电网谐波检测方法的待分析波形的波形图;
图3B是根据本发明实施例的电网谐波检测方法提取的谐波的波形图;
图3C是根据本发明实施例的电网谐波检测方法提取的基波的波形图。
具体实施方式
需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
图1是根据本发明实施例的电网谐波检测装置的示意图,如图所示,本发明实施例的电网谐波检测装置包括:数据离散模块11,用于获取待测电网中的交流电量模拟信号,并将模拟信号进行离散化处理,得到离散数据x(t);经验模态分解模块13,用于在离散数据x(t)的基础上叠加多组不同的高斯白噪声信号ωi(t),得到多组原始信号Xi(t),Xi(t)=x(t)+ωi(t),并对多组原始信号Xi(t)分别进行经验模态分解,得到各阶固有模态函数分量均值运算模块15,用于将多组原始信号Xi(t)分解得到的固有模态函数分量进行均值运算,得到第n次固有模态函数分量cn(t),其公式表示为:结果输出模块17,用于将第n次固有模态函数分量cn(t)代入一组原始信号X(t),则该组X(t)的谐波分解结果为:其中,X(t)=x(t)+ω(t),rm(t)为高斯白噪声信号ω(t)的各阶固有模态函数分量。
本发明实施例的电网谐波检测装置通过数据离散模块11可以将各种方式测量任意电压、电流信号,并进行离散化,均值运算模块15利用高斯白噪声频谱的零均值原理,消除高斯白噪声作为时域分布参考结构带来的影响,其原理为:当附加的白噪声均匀分布在整个时频空间时,该时频空间就由滤波器组分割成的不同尺度成分组成。当信号加上均匀分布的白噪声背景时,不同尺度的信号区域将自动映射到与背景白噪声相关的适当尺度上去。当然,每个独立的测试都可能会产生非常嘈杂的结果,这是因为每个附加噪声的成分都包括了信号和附加的白噪声。既然在每个独立的测试中噪声是不同的,当使用足够测试的全体均值时,噪声将会被消除。全体的均值最后将会被认为是真正的结果,唯一持久稳固的部分是信号本身,所加入的多次测试是为了消除附加的噪声。而高斯白噪声具有零均值特性,在使用足够测试的高斯白噪声均值时,累加的噪声分量为零,从而得到待分析信号自身的性质。
优选地,数据离散模块11可以具体包括:向量初始子模块,用于初始化维数为2N的向量I2N;数据采样子模块,用于获取连续的2N个交流电量的采样值,这2N个采样值记为ik-(2N-1),L,ik-1,ik,将ik-(N-1),L,ik-1;向量构建子模块,用于将ik-(N-1),L,ik-1赋值给向量I2N的前N-1个元素,将ik-(2N-1),L,ik-N赋值给向量I2N的后N个元素,向量I2N表示为:[ik-(N-1) Lik-1 ik ik-2N+1 L ik-N];将赋值后的向量I2N作为离散数据x(t)。这样组成的向量可以保证使当前时刻采样值ik始终处于向量I2N的中间位置。
在实时采集交流电量时,数据离散模块11可以并记录所述交流电量的采样值;将当前采集到的采样值记录为ik,从获取当前采样时刻之前连续2N-1个采样时刻的采样数据。从而可以实现在线检测。
本发明提供的电网谐波检测装置还可以包括:总谐波分量计算模块,用于对赋值后的向量I2N依次送入经验模态分解模块13、均值运算模块15、结果输出模块17进行计算,得到谐波分解结果;从得到的谐波分解结果中查询到工频50Hz的固有模态函数分量cf,cf中第N个元素即为k时刻的基波分量ikf;利用公式ikh=ik-ikf,计算得到k时刻的总谐波分量ikh。从而简便的计算得出k时刻的总谐波分量ikh
此外,本发明提供的电网谐波检测装置还可以包括:特定次谐波分析模块,用于从得到的谐波分解结果中查询到特定次谐波的固有模态函数分量cfn,cfn中第N个元素即为k时刻的基波分量ikfn
上述电网谐波检测装置中使用的经验模态分解(Empirical ModeDecomposition,EMD)是美籍华人N.E.Huang等人于1998年提出的,适合于分析非线性、非平稳信号序列,具有很高的信噪比。该方法的关键是经验模式分解,它能使复杂信号分解为有限个固有模态函数(Intrinsic Mode Function,IMF),所分解出来的各固有模态函数分量包含了原信号的不同时间尺度的局部特征信号。
本发明实施例提供的电网谐波检测装置,在经验模态分解基础上利用高斯白噪声辅助进行数据分析的方法。利用了高斯白噪声频谱的零均值原理,通过多组信号的叠加均值计算消除高斯白噪声作为时域分布参考结构带来的影响,对原始电流或电压信号的采样要求低,不需要电网同步角度,可用于单相、三相三线制或三相四线制系统中检测谐波分量。计算方法简单快捷,谐波检测精度高。
本发明实施例还提供了一种电网谐波检测方法,电网谐波检测方法可以通过本发明上述实施例所提供的任一种电网谐波检测装置来执行,图2是根据本发明实施例的电网谐波检测方法的示意图,如图所示:
本发明实施例的电网谐波检测方法包括:
步骤S11,获取待测电网中的交流电量模拟信号,并将模拟信号进行离散化处理,得到离散数据x(t);
步骤S13,在离散数据x(t)的基础上叠加多组不同的高斯白噪声信号ωi(t),得到多组原始信号Xi(t),Xi(t)=x(t)+ωi(t),并对多组原始信号Xi(t)分别进行经验模态分解,得到各阶固有模态函数分量
步骤S15,将多组原始信号Xi(t)分解得到的固有模态函数分量进行均值运算,得到第n次固有模态函数分量cn(t),表示为:
步骤S17,将第n次固有模态函数分量cn(t)代入一组原始信号X(t),则该组X(t)的谐波分解结果为:其中,X(t)=x(t)+ω(t),rm(t)为高斯白噪声信号ω(t)的各阶固有模态函数分量。
通过以上步骤,本发明实施例的电网谐波检测方法,将任意电网中的交流电量信号做离散化处理,利用高斯白噪声具有零均值特性,在使用足够测试的高斯白噪声均值时,累加的噪声分量为零,经过一系列运算得到待分析信号自身的性质。
其中,步骤S11可以具体包括:初始化维数为2N的向量I2N;获取连续的2N个交流电量的采样值,这2N个采样值记为ik-(2N-1),L,ik-1,ik,将ik-(N-1),L,ik-1;将ik-(N-1),L,ik-1赋值给向量I2N的前N-1个元素,将ik-(2N-1),L,ik-N赋值给向量I2N的后N个元素,向量I2N表示为:[ik-(N-1) L ik-1 ik ik-2N+1 L ik-N];将赋值后的向量I2N作为离散数据x(t)。这样组成的向量可以保证使当前时刻采样值ik始终处于向量I2N的中间位置。
如果第k个采样值ik为当前时刻采集到的采样值,那么利用最近2N个交流电量的采样值进行计算就可以实现在线检测,此时获取连续的2N个交流电量的采样值的步骤包括:实时采集交流电量,并记录交流电量的采样值;将当前采集到的采样值记录为ik,从获取当前采样时刻之前连续2N-1个采样时刻的采样数据。
上述将赋值后的向量I2N作为离散数据x(t)的步骤之后还包括:对赋值后的向量I2N依次送入经验模态分解模块、均值运算模块、结果输出模块进行计算,得到谐波分解结果;从得到的谐波分解结果中查询到工频50Hz的固有模态函数分量cf,cf中第N个元素即为k时刻的基波分量ikf;利用公式ikh=ik-ikf,计算得到k时刻的总谐波分量ikh。从而可以简单地得到总谐波分量。
上述得到谐波分解结果的步骤之后还可以包括:从得到的谐波分解结果中查询到特定次谐波的固有模态函数分量cfn,cfn中第N个元素即为k时刻的基波分量ikfn。从而实现了检测特定次谐波的效果。
上述获取待测电网中的交流电量模拟信号的步骤具体包括:获取一下任意一种电网电路中的交流电量模拟信号:单相电路、三相三线制电路和三相四线制电路中,其中交流电量包括电流信号和/或电压信号。因此,本实施例的谐波分析方法并不限于任何一种电网形式,适用广泛。
本实施例的电网谐波检测方法,可以实现谐波的离线检测、谐波的在线检测、以及特定频率的检测。
其中,谐波的离线检测方法的具体步骤为:
步骤S21,通过各种方式测量任意电压、电流信号,并进行离散化,形成离散数据x(t)。
步骤S22,通过给原始信号x(t)叠加一组高斯白噪声信号ω(t)获得一个总体信号X(t)。
X(t)=x(t)+ω(t)
步骤S23,对X(t)进行EMD分解,得到各阶IMF分量
步骤S24,给原始信号加入不同的白噪声ωi(t),重复以上步骤S22和步骤S23;
步骤S25,利用高斯白噪声频谱的零均值原理,消除高斯白噪声作为时域分布参考结构带来的影响,原始信号对应的IMF分量cn(t)可表示为
步骤S26,最后,原始信号x(t)可分解为:
谐波的在线检测方法包括以下步骤:
步骤S31,初始化维数为2N的向量I2N
步骤S32,当前时刻的采样值为ik,前面2N-1个时刻的数据分别为ik-(2N-1),L,ik-1,将ik-(N-1),L,ik-1赋值给向量I2N的前N-1个元素,ik-(2N-1),L,ik-N赋值给向量I2N的后N个元素,使当前时刻采样值ik始终处于向量I2N的中间位置。这样,2N维向量I2N可表示为:
[ik-(N-1) L ik-1 ik ik-2N+1 L ik-N]
步骤S33,将步骤S32中赋值后的向量I2N作为原始信号,按照离线方法进行计算,分解出各阶IMF,并计算各阶IMF的频率,找到工频50Hz附近的IMF分量cf,其中第N个元素即为k时刻的基波分量ikf
ikf=cf(N)
步骤S34,利用当前时刻的采样值ik减掉步骤S33中计算出的基波分量ikf,即可得到信号中的总谐波分量ikh:ikh=ik-ikf
谐波的特定频率检测方法为:修改在线检测方法步骤S33中IMF的频率即可求出对应的谐波分量
下面以典型三相可控整流负载的任一相电流为例对使用本发明实施例的技术效果进行说明,图3A是根据本发明实施例的电网谐波检测方法的待分析波形的波形图,图3B是根据本发明实施例的电网谐波检测方法提取的谐波的波形图,图3C是根据本发明实施例的电网谐波检测方法提取的基波的波形图。由以上附图可以看出,使用本实施例的电网谐波检测方法对电网模拟信号处理后,得到的谐波检测结果精确。
应用本发明的技术方案,基于经验模态分解方法的不足,提出了利用高斯白噪声辅助进行数据分析的方法。利用了高斯白噪声频谱的零均值原理,通过多组信号的叠加均值计算消除高斯白噪声作为时域分布参考结构带来的影响,对原始电流或电压信号的采样要求低,不需要电网同步角度,可用于单相、三相三线制或三相四线制系统中检测谐波分量。计算方法简单快捷,谐波检测精度高。可以应用于各种谐波检测仪器、谐波检测软件程序、有源电力滤波器等不同的场合。
显然,本领域的技术人员应该明白,上述的本发明的各模块或各步骤可以用通用的计算装置来实现,它们可以集中在单个的计算装置上,或者分布在多个计算装置所组成的网络上,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样,本发明不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (10)

1.一种电网谐波检测方法,其特征在于,包括
步骤S11,获取待测电网中的交流电量模拟信号,并将所述模拟信号进行离散化处理,得到离散数据x(t);
步骤S13,在所述离散数据x(t)的基础上叠加多组不同的离散高斯白噪声信号ωi(t),得到多组原始信号Xi(t),所述Xi(t)=x(t)+ωi(t),并对所述多组原始信号Xi(t)分别进行经验模态分解,得到各阶固有模态函数分量
步骤S15,将所述多组原始信号Xi(t)分解得到的固有模态函数分量进行均值运算,得到第n次固有模态函数分量cn(t),表示为:
步骤S17,将所述第n次固有模态函数分量cn(t)代入一组原始信号X(t),则该组X(t)的谐波分解结果为:其中,X(t)=x(t)+ω(t),rm(t)为模拟高斯白噪声信号ω(t)的各阶固有模态函数分量。
2.根据权利要求1所述的电网谐波检测方法,其特征在于,所述步骤S11包括:
初始化维数为2N的向量I2N
获取连续的2N个所述交流电量的采样值,这2N个采样值记为ik-(2N-1),…,ik-1,ik,将ik-(N-1),…,ik-1
将ik-(N-1),…,ik-1赋值给所述向量I2N的前N-1个元素,将ik-(2N-1),…,ik-N赋值给向量I2N的后N个元素,所述向量I2N表示为:[ik-(N-1) … ik-1 ik ik-2N+1 … ik-N];
将赋值后的所述向量I2N作为所述离散数据x(t)。
3.根据权利要求2所述的电网谐波检测方法,其特征在于,将赋值后的所述向量I2N作为所述离散数据x(t)之后还包括:
对赋值后的所述向量I2N执行所述步骤S13至所述步骤S17,得到谐波分解结果;
从得到的谐波分解结果中查询到工频50Hz的固有模态函数分量cf,cf中第N个元素即为k时刻的基波分量ikf
利用公式ikh=ik-ikf,计算得到k时刻的总谐波分量ikh
4.根据权利要求3所述的电网谐波检测方法,其特征在于,得到谐波分解结果之后还包括:
从得到的谐波分解结果中查询到特定次谐波的固有模态函数分量cfn,cfn中第N个元素即为k时刻的基波分量ikfn
5.根据权利要求2至4中任一项所述的电网谐波检测方法,其特征在于,获取连续的2N个所述交流电量的采样值包括:
实时采集交流电量,并记录所述交流电量的采样值;
将当前采集到的采样值记录为ik,从获取当前采样时刻之前连续2N-1个采样时刻的采样数据。
6.根据权利要求1至4中任一项所述的电网谐波检测方法,其特征在于,获取待测电网中的交流电量模拟信号包括:
获取一下任意一种电网电路中的交流电量模拟信号:单相电路、三相三线制电路和三相四线制电路中,其中所述交流电量包括电流信号和/或电压信号。
7.一种电网谐波检测装置,其特征在于,包括:
数据离散模块,用于获取待测电网中的交流电量模拟信号,并将所述模拟信号进行离散化处理,得到离散数据x(t);
经验模态分解模块,用于在所述离散数据x(t)的基础上叠加多组不同的高斯白噪声信号ωi(t),得到多组原始信号Xi(t),所述Xi(t)=x(t)+ωi(t),并对所述多组原始信号Xi(t)分别进行经验模态分解,得到各阶固有模态函数分量
均值运算模块,用于将所述多组原始信号Xi(t)分解得到的固有模态函数分量进行均值运算,得到第n次固有模态函数分量cn(t),其公式表示为:
结果输出模块,用于将所述第n次固有模态函数分量cn(t)代入一组原始信号X(t),则该组X(t)的谐波分解结果为:其中,X(t)=x(t)+ω(t),rm(t)为高斯白噪声信号ω(t)的各阶固有模态函数分量。
8.根据权利要求7所述的电网谐波检测装置,其特征在于,所述数据离散模块包括:
向量初始子模块,用于初始化维数为2N的向量I2N
数据采样子模块,用于获取连续的2N个所述交流电量的采样值,这2N个采样值记为ik-(2N-1),…,ik-1,ik,将ik-(N-1),…,ik-1
向量构建子模块,用于将ik-(N-1),…,ik-1赋值给所述向量I2N的前N-1个元素,将ik-(2N-1),…,ik-N赋值给向量I2N的后N个元素,所述向量I2N表示为:[ik-(N-1) … ik-1 ik ik-2N+1… ik-N];
将赋值后的所述向量I2N作为所述离散数据x(t)。
9.根据权利要求8所述的电网谐波检测装置,其特征在于,还包括:
总谐波分量计算模块,用于将赋值后的所述向量I2N依次送入所述经验模态分解模块、所述均值运算模块、所述结果输出模块进行计算,得到谐波分解结果;从得到的谐波分解结果中查询到工频50Hz的固有模态函数分量cf,cf中第N个元素即为k时刻的基波分量ikf;利用公式ikh=ik-ikf,计算得到k时刻的总谐波分量ikh
10.根据权利要求9所述的电网谐波检测装置,其特征在于,还包括:
特定次谐波分析模块,用于从得到的谐波分解结果中查询到特定次谐波的固有模态函数分量cfn,cfn中第N个元素即为k时刻的基波分量ikfn
CN201210586346.5A 2012-12-28 2012-12-28 电网谐波检测方法和电网谐波检测装置 Expired - Fee Related CN103901273B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210586346.5A CN103901273B (zh) 2012-12-28 2012-12-28 电网谐波检测方法和电网谐波检测装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210586346.5A CN103901273B (zh) 2012-12-28 2012-12-28 电网谐波检测方法和电网谐波检测装置

Publications (2)

Publication Number Publication Date
CN103901273A CN103901273A (zh) 2014-07-02
CN103901273B true CN103901273B (zh) 2017-12-22

Family

ID=50992727

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210586346.5A Expired - Fee Related CN103901273B (zh) 2012-12-28 2012-12-28 电网谐波检测方法和电网谐波检测装置

Country Status (1)

Country Link
CN (1) CN103901273B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104360153B (zh) * 2014-10-22 2017-02-08 江苏省电力公司淮安供电公司 电网谐波在线检测与分析方法
CN105067881A (zh) * 2015-07-24 2015-11-18 淄博格瑞恩电气科技有限公司 一种基于谐波检测的用电安全监控方法
CN107046445B (zh) * 2016-02-06 2020-07-28 富士通株式会社 谐波失真分离方法、非线性特性确定方法、装置和系统
CN105866571B (zh) * 2016-03-25 2018-09-21 浙江工业大学 一种基于高频谐波补偿迭代emd的暂态电能质量信号分析方法
CN106253275A (zh) * 2016-08-02 2016-12-21 上海电机学院 一种基于瞬时功率模型的谐波检测方法以及有源电力滤波器
CN106546818B (zh) * 2016-10-20 2019-04-09 南京航空航天大学 一种基于差分非线性模式分解的谐波信号检测方法
CN107064633B (zh) * 2017-03-29 2019-10-18 广西电网有限责任公司电力科学研究院 城市轨道交通负荷谐波电流迭加系数确定方法
CN113030568B (zh) * 2021-02-25 2022-01-07 南方电网科学研究院有限责任公司 一种直流输电系统谐波风险评估方法及装置

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005015759A1 (en) * 2003-07-16 2005-02-17 Ericsson Technology Licensing Ab Transceiver architecture with reduced vco-pulling sensitivity
CN1609630A (zh) * 2004-10-21 2005-04-27 上海交通大学 提取混沌干扰下的谐波信号的方法
CN101408567A (zh) * 2008-11-28 2009-04-15 北京四方继保自动化股份有限公司 基于经验模态分解的大规模电网低频振荡分频段检测法
CN102170126A (zh) * 2011-04-06 2011-08-31 华中电网有限公司 一种电力系统低频振荡在线分析及预警方法
CN102831433A (zh) * 2012-06-06 2012-12-19 西南交通大学 一种基于时频域多特征量的电能质量混合扰动分类方法
CN102841251A (zh) * 2012-08-14 2012-12-26 华北电力大学 一种微网中电动汽车充电站谐波检测法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6498736B1 (en) * 2001-03-27 2002-12-24 Baldor Electric Company Harmonic filter with low cost magnetics

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005015759A1 (en) * 2003-07-16 2005-02-17 Ericsson Technology Licensing Ab Transceiver architecture with reduced vco-pulling sensitivity
CN1609630A (zh) * 2004-10-21 2005-04-27 上海交通大学 提取混沌干扰下的谐波信号的方法
CN101408567A (zh) * 2008-11-28 2009-04-15 北京四方继保自动化股份有限公司 基于经验模态分解的大规模电网低频振荡分频段检测法
CN102170126A (zh) * 2011-04-06 2011-08-31 华中电网有限公司 一种电力系统低频振荡在线分析及预警方法
CN102831433A (zh) * 2012-06-06 2012-12-19 西南交通大学 一种基于时频域多特征量的电能质量混合扰动分类方法
CN102841251A (zh) * 2012-08-14 2012-12-26 华北电力大学 一种微网中电动汽车充电站谐波检测法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
EEMD在电能质量扰动检测中的应用;张杨 等;《电力自动化设备》;20111231;第31卷(第12期);86-91 *

Also Published As

Publication number Publication date
CN103901273A (zh) 2014-07-02

Similar Documents

Publication Publication Date Title
CN103901273B (zh) 电网谐波检测方法和电网谐波检测装置
Lee et al. Online embedded impedance measurement using high-power battery charger
Mazumdar et al. Recurrent neural networks trained with backpropagation through time algorithm to estimate nonlinear load harmonic currents
Rens et al. On techniques for the localisation of multiple distortion sources in three‐phase networks: Time domain verification
Mohammadi-Ivatloo et al. Online small signal stability analysis of multi-machine systems based on synchronized phasor measurements
CN103401238B (zh) 一种基于总体测辨法的电力负荷建模方法
CN107851995A (zh) 用于检测供给电网中的电压的方法和设备
CN104793082A (zh) 基于谐波相关分析的电力系统谐波源识别装置
CN203133168U (zh) 电力谐波检测仪
Gallo et al. Low cost smart power metering
CN104267258A (zh) 一种利用不完全s变换的谐波瞬时功率计算方法
CN114518488A (zh) 基于电能参数搜索量测高精度物联智能电表分布式光储充放系统
CN109687474A (zh) 基于卡尔曼滤波的指令电流检测算法及系统
Zhang et al. Research on power system harmonic state estimation
CN104850751A (zh) 一种电流质量评估方法
Folly et al. Critical clearing time prediction using recurrent neural networks
An et al. Supraharmonics Measurement Based on Hybrid Online Measurement and Offline Analysis
Wang et al. Total harmonic distortion (THD) estimation technique based on power concept for smart power meters
CN203069669U (zh) 一种基于参数分析的谐波电能计量系统
Kumar et al. A Comparative Analysis of Harmonics Distortion and Transient Monitor in LABVIEW
Pourarab et al. Impact of variations in the utility side on the harmonic contribution of customer installations
CN108008219A (zh) 风光互补仿真实验装置
Yuan et al. Harmonic characteristics of distributed generation and electric vehicle integrating to the grid
CN109831027A (zh) 一种配电网电能质量灵活检测分析系统
CN106291103A (zh) 一种基于eemd的谐波检测方法

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20200907

Address after: 102199 b-702, building 1, No.2 Yuncheng street, Beijing Economic and Technological Development Zone, Daxing District, Beijing

Patentee after: BEIJING ZHIZHONG ENERGY TECHNOLOGY DEVELOPMENT Co.,Ltd.

Address before: 100084 Qingyuan 12-3-201, Tsinghua University, Beijing, Haidian District

Co-patentee before: Dong Weijie

Patentee before: Bai Xiaomin

Co-patentee before: Zhu Ninghui

CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20171222

Termination date: 20211228