CN102841251B - 一种微网中电动汽车充电站谐波检测法 - Google Patents
一种微网中电动汽车充电站谐波检测法 Download PDFInfo
- Publication number
- CN102841251B CN102841251B CN201210289227.3A CN201210289227A CN102841251B CN 102841251 B CN102841251 B CN 102841251B CN 201210289227 A CN201210289227 A CN 201210289227A CN 102841251 B CN102841251 B CN 102841251B
- Authority
- CN
- China
- Prior art keywords
- signal
- continuation
- microgrid
- frequency
- harmonic
- 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
Links
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
本发明公开了电气工程微网技术领域中的一种微网中电动汽车充电站谐波检测法。本发明首先对微网中谐波电流进行预处理,将其分解成窄带频率信号族S1;然后运用自适应分步延拓法对其进行延拓处理,得到频率信号族S;进一步对频率信号族S进行经验模态分解,得到各个固有模态分量;对固有模态分量进行Hilbert变换,最后得到各分量的瞬时频率信号和瞬时幅值信号。本发明不仅适用于分析微网电能质量中的非线性、非平稳信号,而且对于微网电能质量中的线性、平稳信号的分析比其他时频分析方法更好地反映了信号的频率特征以及幅值特征,具有更明确的物理意义。
Description
技术领域
本发明属于电气工程微网技术领域,尤其涉及一种微网中电动汽车充电站谐波检测法。
背景技术
随着石油资源的日益枯竭以及人们对城市空气污染的倍加关注,电动汽车开始受到全世界的青睐,各国政府和工业界均在加大力度支持电动汽车的发展。电动汽车以电代油,能够实现零排放的优点使其成为解决能源和环境问题的重要手段。我国电动汽车起步较晚,但是发展很快。为了满足电动汽车的充电需求,国内正在大力推进电动汽车充电站的建设。2010年4月,国家有关部门审查并通过了《电动汽车传导式充电接口》、《电动汽车充电站通用要求》、《电动汽车电池管理系统与非车载充电机之间的通信协议》、《轻型混合动力电动汽车能量消耗量试验方法》。本次公布的4项国家标准,对电动汽车充电接口、充电站标准和功能等情况进行了细化和规范。这意味着日后国内生产的电动汽车充电接口将统一化,同时将大大促进目前充电站的有序建设。
微网作为分布式电源并网发电规模化应用的有效技术途径,是一种由微电源、储能装置及负荷共同组成的有机系统。其通过有效的协调控制,使主要基于新能源和可再生能源的分布式电源并网所产生的负面问题都在微网内得到解决,减少了分布式电源并网对主电网产生的各种扰动,为电动汽车并网的实现提供了良好的平台。然而,电动汽车充电机是一种典型的非线性设备,其产生 的大量谐波注入微网,不仅会对微网造成污染,同时也将导致充电机功率因数的降低,而谐波和功率因数过低都会对微网和主电网造成危害。为此,有必要对微网中电动汽车充电站产生的谐波电流进行检测,分析影响其变化的主要因素。
传统的谐波检测方法主要有快速傅里叶变换(FFT)以及近些年来被广泛运用于电能质量检测的小波变换(WT)等等。对于整数次谐波,FFT能进行精确的检测,但是对于非整数次谐波,运用FFT方法可能会产生频谱泄露和栅栏效应。虽然采用加窗和插值等方法可以有效进行改善,但是仍有许多不足之处,因为这些算法都是以降低频率分辨率为代价的。与傅里叶变换相比,小波分析方法则是一种时频局部化分析方法。但是在小波变换中,经过高通和低通滤波器分解后的各个频带之间存在严重的混叠现象,因此难以实现各个频带的严格划分。不仅如此,在分析谐波信号时,有很多小波基函数可供选择,小波基的选择对于信号的分析结果有很大的影响,但是最优基的选择没有确切的规律可循,在很大程度上取决于设计者的经验。
HHT(Hilbert-Huang Transform)是Norden E.Huang等人于1998年提出的一种完全自适应的时频分析方法,HHT方法由经验模态分解(Empirical Mode Decomposition,EMD)及Hilbert变换(HT)两部分组成。其核心部分在于对信号进行EMD分解,得到信号的固有模态分量(Intrinsic Mode Function,IMF),然后对IMF分量进行Hilbert变换,从而得到信号的瞬时频率和瞬时幅值。HHT非常适合于处理非线性、非平稳信号,但是在运用其对信号进行处理的同时也存在不少问题,其中端点效应和频谱混叠是重要而且难以解决的问题。在对信号进行EMD分解时,当信号的端点处不是极值点会导致信号包络线拟合的过程中产生误差,而这个误差就是产生端点效应问题的原因。不仅如此,在对信号进行Hilbert变换时,信号的端点处也会产生端点效应问题。另外,EMD分解是一个由高频到低频的过程,高频信号总是出现在第一个IMF分量中。因此,分 解出来的第一个IMF分量往往不是单分量信号,包含了较宽的频率范围,这就产生了频谱混叠现象,使得Hilbert变换所得到的瞬时频率和瞬时幅值失去了原来的物理意义。
发明内容
针对上述背景技术中提到的现有时频分析方法不能分析非线性、非平稳信号等不足,本发明提出了一种微网中电动汽车充电站谐波检测法。
本发明的技术方案是,一种微网中电动汽车充电站谐波检测法,其特征是该方法包括以下步骤:
步骤1:运用小波变换的方法对微网中电动汽车充电机的谐波电流进行预处理,将其分解成不同频率段的信号族S1,S1={s1t},t=1,2,…n;
步骤2:用自适应分步延拓法对步骤1的信号族S1进行延拓处理,得到信号族S,S={s0t};
步骤3:对频率信号族S的每个信号进行经验模态分解,得到每个信号的各个固有模态分量;
步骤4:对所有的固有模态分量进行Hilbert变换,在得到的瞬时频率信号和瞬时幅值信号中去掉因延拓而增加的信号段,得到最终的信号。
所述自适应分步延拓法为:
a:对S1中的每个频率信号进行延拓,得到信号族S2,S2={s2t};
b:对S2中的每个信号进行延拓,包括,
情况一:当信号s2t的左端或右端最后一个极值点为极大值点时,以该极大值点为对称轴进行向左或向右延拓;
情况二:当信号s2t的左端或右端最后一个极值点为极小值点时,以该极小值点为对称轴进行向左或向右延拓;
情况三:当信号s2t的左端或右端处的函数值大于离该端点最近的第一个极大值或者小于离该端点最近的第一个极小值时,把端点当作延拓的第一个极值点或者作为信号s2t的最后一个极值点,然后以该指定端点为对称轴进行向左或向右延拓。
所述步骤3具体为:
步骤3.1:通过三次样条插值法求得频率信号s0t的上包络线和下包络线,求得其上包络线和下包络线的均值;
步骤3.2:求得频率信号s0t与上包络线和下包络线的均值的差值,若所述差值满足指定条件,则所述差值为固有模态分量;否则,用所述差值替换频率信号s0t,返回步骤3.1;
指定条件为:若所述差值与前一个差值代入判据SDk后所得的值小于一个指定门限值时,则该差值为固有模态分量;
步骤3.3:将频率信号s0t和步骤3.2中的固有模态分量的差值记为残差,若该残差满足设定条件,则过程结束;否则,将该残差替换频率信号s0t,返回步骤3.1;
设定条件为:残差为单调函数或残差的模值为一常数。
本发明所要得到的效果是:对电动汽车充电站并网点所产生的谐波电流进行实时的检测。与其他时频分析方法相比,本专利所提出的检测方法不仅适用于分析微网电能质量中的非线性、非平稳信号,而且对于电能质量中的线性、平稳信号的分析比其他时频分析方法更好地反映了信号的频率特征和幅值特 征,具有更明确的物理意义。对于谐波而言,这种时频分析方法能根据信号自身的特性进行自适应分解,且不存在基函数的选择问题,从而可以实现微网中谐波的自动提取。与原HHT方法相比,本发明所提出的检测方法可以有效地抑制原HHT变换过程中所产生的端点效应问题和频谱混叠问题,从而可以更为有效地对谐波电流的频率和幅值进行实时的监测,有利于分析影响其变化的主要因素。为今后电动汽车充电站的谐波治理提供理论依据和技术支撑,使得微网中的谐波含量符合国家电能质量相关标准。
本发明所提出的谐波检测方法简捷易于实现,对微网运行工况适应性强,便于实际工程应用。
附图说明
图1为本发明所提出的检测方法的基本流程图;
图2为EMD的分解流程图;
图3为短时间谐波序列波形图;
图4为信号的包络线示意图;
图5为原信号的EMD分解结果;图5(a)为IMF1分量;图5(b)为IMF2分量;图5(c)为IMF3分量;
图6为延拓信号的EMD分解结果;图6(a)为IMF1分量;图6(b)为IMF2分量;图6(c)为IMF3分量;
图7为延拓前后信号瞬时频率对比图;图7(a)为延拓前信号的瞬时频率;图7(b)为延拓后信号的瞬时频率;
图8为整数次谐波信号波形图;
图9为谐波信号的EMD分解结果;
图10为IMF1的FFT分析结果;
图11为经过小波变换预处理后得到的EMD分解结果;图11(a)为IMF1分 量;图11(b)为IMF2分量;
图12为运用本发明所提出的检测方法所得到的瞬时频率和瞬时幅值;图12(a)为谐波信号的瞬时频率;图12(b)为谐波信号的瞬时幅值;
图13为时变谐波信号波形图;
图14为时变谐波信号的瞬时幅值和瞬时频率;图14(a)为谐波信号的瞬时幅值;图14(b)为谐波信号的瞬时频率;
图15为微网中电动汽车充电站的并网示意图;
图16为高频充电机的原理图;
图17为高频充电机的电路图;
图18为充电机的仿真原理图;
图19为电动汽车充电站的仿真模型;
图20为电动汽车充电站并网点处流过的A相谐波电流波形;
图21为运用本发明所提出的检测方法所得到的谐波电流的瞬时频率和瞬时幅值;图21(a)为谐波电流的瞬时频率;图21(b)为谐波电流的瞬时幅值。
具体实施方式
下面结合附图,对优选实施例作详细说明。应该强调的是,下述说明仅仅是示例性的,而不是为了限制本发明的范围及其应用。
对于微网中电动汽车充电站的谐波电流,根据本发明所提出的检测方法,可以按照图1的流程进行处理,其具体过程如下所述:
步骤1:运用小波变换的方法对电动汽车充电机的谐波电流进行预处理,将其分解成不同频率段的窄带信号族S1,S1={s1t},t=1,2,…n;(n与小波变换的分解层数有关),对原信号的频率进行细化,从而避免频谱混叠现象的产生;
步骤2:为了改善HHT中的端点效应问题,在此采用自适应分步延拓法对 步骤1得到的窄带信号族S1进行延拓,最后得到信号族S,S={s0t}。
谐波电流经过小波变换预处理后成为窄带信号族S1,S1={S1t}。由于S1中的每一个窄带信号的延拓方法都相同,此处设其中一个谐波电流窄带信号为s1t,其数据序列为s1t(1),s1t(2)…,s1t(N),其中,N为采样点数。运用自适应分步延拓法对其进行处理的具体过程如下所示:
首先构造一个如式(1)所示的回归模型,其中,ai为非负的拉格朗日乘子;b为偏置项;k(x,xi)为核函数;l为训练样本数。
在所构造的回归模型的基础上对所要延拓的数据序列进行向左和向右延拓。向左、向右延拓方法相同,下面以向右延拓为例进行分析:首先确定要延拓的点数m和训练样本数l,然后按照一定的规则形成一个训练样本集T={(x1,y1),(x2,y2),…,(xl,yl)},其中,xi=[s1t(i),s1t(i+1),…,s1t(N-l+i-1)]T,yi=s1t(N-l+i),1≤i≤l。利用如式(1)所示的回归模型就可以得到数据序列边界外的第一个预测值s1t(N+1)。
式中,x={x1,x2…,xl}T。再将s1t(N+1)作为原始谐波信号序列新的边界点,就可得到下一个预测值s1t(N+2),以此类推,根据所要延拓的数据个数就可以得到全部的延拓序列s1t(N+1),…,s1t(N+m)。
令s2t={s1t(1),…,s1t(N+m)}。设s2t对应的时间序列为:{t1,t2…,tn},则s2t可 以表示为:{s2t(t1),s2t(t2),…,s2t(tn)}。又设此时该信号有M个极大值和N个极小值,极大值点所对应的时间为:TM(i),函数值为:s2tM(i),1≤i≤M;极小值所对应的时间为:TN(i),函数值为:s2tN(i),1≤i≤N。
对s2t继续进行延拓,由于其端点处的波形和函数值大小情况各不相同,延拓时必须分情况讨论。下面仍以向右延拓为例进行分析:
情况一:当信号s2t右端最后一个极值点为极大值点时,即TM(M)>TN(N),以此处极大值点为对称轴向右延拓:
TM(M+j)=2TM(M)-TM(M-j),s2tM(M+j)=s2tM(M-j);
TN(N+j)=2TM(M)-TN(N-j+1),s2tN(N+j)=s2tN(N-j+1)。
j为延拓的第j个极值点,1≤j≤k,k为所要延拓的极值点数。
情况二:当信号s2t右端最后一个极值点为极小值点时,即TN(N)>TM(M),以此处极小值点为对称轴向右延拓:
TM(M+j)=2TN(N)-TM(M-j+1),s2tM(M+j)=s2tM(M-j+1);
TN(N+j)=2TN(N)-TN(N-j),s2tN(N+j)=s2tN(N-j)。
则右端延拓的极大值点所对应的时间为TM(M+1),TM(M+2),…,TM(M+k);极大值点所对应的函数值为s2tM(M+1),s2tM(M+2),…,s2tM(M+k)。右端延拓的极小值点所对应的时间为TN(N+1),TN(N+2),…,TN(N+k);极小值点所对应的函数值为s2tN(N+1),s2tN(N+2),…,s2tN(N+k)。
情况三:当信号s2t右端点处的函数值大于离端点最近的第一个极大值或者小于离端点最近的第一个极小值时,为了使端点落在包络线的包络范围内,就必须对端点进行特殊处理,这时应当把端点当作延拓的第一个极值点或者作为 信号的最后一个极值点,然后以其为对称轴进行其他极值点的延拓,极值点所对应的时间和函数值的推导与前两种情况类似。
向左延拓与向右延拓情况类似,此处不再赘述。在下面叙述中,用s0t表示s1t经自适应分步延拓法延拓后所得到的信号。
步骤3:对延拓后的信号族S={s0t}进行EMD分解,得到一系列IMF分量。由于每个信号EMD分解过程相同,在此以第t个信号s0t的EMD分解过程为例进行说明。
提取s0t谐波信号IMF分量的计算过程如下所示,其算法流程图如图2所示。
(1)找到信号s0t(t)的所有极大值点和极小值点,分别用三次样条插值的方法拟合上包络线v1(t)和下包络线v2(t),并求出其平均值m1(t)=[v1(t)+v2(t)]/2。
(2)求取信号s0t(t)与m1的差值,即h1(t)=s0t(t)-m1(t)。对于不同的信号,h1(t)可能是一个IMF分量,也可能不是。对于分解得到的IMF分量,必须满足以下两个条件:1)其过零点和极值点的数目必须相等或者最多相差一个;2)连接局部极大值所形成的上包络线和连接局部极小值所形成的下包络线的平均值在任一点处都为零。当不满足IMF所需条件时,把h1(t)当作原信号,重复上述步骤,即得:
h11(t)=h1(t)-m11(t)
式中,m11(t)是h1(t)的上、下包络线的平均值。若h11(t)不是IMF分量,则继续筛选,重复上述方法k次,得到第k次筛选的数据h1k(t):
h1k(t)=h1(k-1)(t)-m1(k-1)(t)
以SDk作为判据,把两个连续处理结果代入SDk表达式,当SDk小于某一指 定门限值时,筛选过程即可结束。
经验表明,当门限值取0.2~0.3时,即可保证IMF的线性和稳定性,又可以使IMF具有相应的物理意义。
(3)当h1k(t)满足筛选终止准则的要求,则h1k(t)视为第一个IMF信号,记为c1(t)。从s0t(t)中减去c1(t)得到剩余信号,即残差r1(t):
r1(t)=s0t(t)-c1(t)
将r1(t)视为一组新信号,重复上述过程,经过多次运算可以得到全部的残差ri(t):
ri(t)=ri-1(t)-ci(t) i=1,2,3,…,n
直到rn为单调函数或者|rn|保持不变时,分解即可停止。最终信号可以表示为:
步骤4:对各个IMF分量作Hilbert变换,得到谐波信号的瞬时频率和瞬时幅值,然后从瞬时频率和瞬时幅值信号中减去因延拓而增加的信号段,得到最终的信号。此时信号的长度仍和原始信号s1t的长度保持一致。
现以任意一个IMF分量ci为例进行分析,对其作Hilbert变换,则有:
其反变换为:
得到的解析信号x(t)为:
x(t)=ci(t)+jH[ci(t)] (5)
上式还可以表示为x(t)=a(t)ejθ(t)。其中,a(t)为瞬时幅值;θ(t)为相位。
对于IMF分量的瞬时频率,可以按式(7)进行计算,这样谐波信号的瞬时频率就得到了。
在此设原始谐波信号s1t的采样点数为N,经过自适应分步延拓法延拓后的数据序列s0t的采样点数为N+2m1。则对于所得到的谐波信号的瞬时频率序列f(i)和瞬时幅值序列a(i),1≤i≤N+2m1,可以分别表示为{f(1),f(2),…,f(N+2m1)}和{a(1),a(2),…,a(N+2m1)}。舍去经自适应分步延拓法延拓出来的那部分信号就得到最终结果。瞬时频率可以表示为{f(m1+1),f(m1+2),…,f(m1+N)},瞬时幅值可以表示为{a(m1+1),a(m1+2),…,a(m1+N)}。
上述即为微网中电动汽车充电站谐波电流检测法的全过程,下面结合附图和实例对本发明进一步说明。
为了充分说明本发明所提出的谐波检测方法在谐波问题分析方面的有效性,现分别对以下谐波信号进行分析。
算例1: 其中A1=2A,A2=1.5A,A3=1.5A;f1=4Hz,f2=2Hz,f3=0.5Hz。该信号为一短时间谐波电流序列,采样频率为300,采样点数为301,其波形图如图3所示。
对该短时间谐波信号进行三次样条插值,得到信号的包络线如图4所示。从图中可以看出,原信号只存在3个极大值点和4个极小值点,由于不能确定端点处是否为极值点,所以端点处信号包络的趋势无从得知,得到的包络线不能包含所有的数据。上下包络线在数据两端均背离了原信号,产生了端点效应。同时由于数据序列很短,产生的端点效应将向内“污染”整个数据序列。图5为信号EMD分解结果与原信号IMF分量的对比图。通过对比可以发现,无论是幅值还是频率,分解结果都出现了严重的偏差。所以在采用EMD对信号进行分解时必须先抑制其端点效应的产生,否则分解毫无意义。
为了改善EMD分解过程中产生的端点效应问题,本发明提出采用自适应分步延拓法对数据序列进行延拓。在此,利用网格搜索和交叉验证算法确定其参数,取训练样本数为100,训练样本集为T={(x1,y1),(x2,y2),…,(x100,y100)}。对延拓后的信号进行EMD分解,得到的分解结果如图6所示。从图中可以清楚地看到,对延拓后的信号进行EMD分解所得到的IMF分量与原信号的IMF分量基本重合,与直接进行EMD分解所得的结果相比,IMF波形得到了很大的改善。特别是对于IMF3,直接EMD分解得到的结果只是一条直线,当采用自适应分步延拓法处理后,虽然还存在一定的误差,但是考到其特征尺度最大且只有半个波形,处理效果已经相当令人满意。
对延拓前后的信号分别进行Hilbert变换,得到的瞬时频率对比图如图7所示,延拓前信号的瞬时频率出现了严重的失真现象,当采用自适应分步延拓法对其进行延拓后,瞬时频率具有了物理意义并与原给定信号相一致,可见该方 法在处理端点效应问题上的有效性。
算例2: 其中,A1=10V,A2=1V;f1=50HZ,f2=150HZ。取采样频率fs=3200HZ,采样点数为1600。其波形如图8所示。
对算例2的波形其直接进行EMD分解,得到的结果如图9所示。接着对分解得到的IMF1进行FFT分析,结果如图10所示,由于出现了频谱混叠现象,IMF1的FFT分析结果中含有两个频率分量,这将会导致Hilbert变换得到的瞬时频率和瞬时幅值失去原来的物理意义。
为了解决这一问题,本发明提出采用小波变换的方法对谐波信号进行预处理,提高信号特征提取的准确性。在此,选择‘db5’小波作为小波基函数,对信号进行5层分解,对于重构后的高低频系数,采用自适应分步延拓法对其进行延拓处理,得到的EMD分解结果如图11所示。接着对IMF1和IMF2进行Hilbert变换,得到的瞬时频率和瞬时幅值如图12所示。从得到的结果可知利用小波变换方法对谐波信号进行预处理使得EMD分解中的频谱混叠问题得到了很好的解决。
算例3:算例1与算例2均为平稳谐波信号,为了验证本文所提出的谐波检测方法在分析非平稳谐波信号上的有效性,在此以一时变间谐波信号为例进行分析,其表达式为 式中,fN=50Hz,f1=110Hz,f2=325Hz。取采样频率fs=3200Hz。该谐波信号的波形如图13所示。
为了分析该谐波信号的幅值特征和频率特征,在此采用本发明所提出的检测方法对其进行检测。所得到的瞬时频率和瞬时幅值如图14所示。从瞬时幅值图中可以清楚地看到,IMF1分量的瞬时幅值在0.2s和0.3s处发生了突变;从瞬 时频率图中可以看到,其瞬时频率在0.3s处发生了突变,与原给定信号相一致。由此可见,本发明所提出的检测方法不仅能精确地检测出扰动发生和终止的时间,而且可以准确地检测出各次谐波的幅值和频率信息。
以上三个仿真算例所得的结果均表明,本发明所提出的微网中谐波检测方法能够准确提取各个频率分量的参数并且有助于得到具有物理意义的瞬时频率和瞬时幅值,不仅适用于微网中线性、平稳谐波信号的分析,而且适用于非线性、非平稳信号的分析,从真正意义上实现了整数次谐波和非整数次谐波的检测。因此,本发明将运用该方法对微网中电动汽车充电站的谐波进行检测。
参见图15,该图为电动汽车充电站并入微网的结构示意图。该微网系统由微电源(微型燃气轮机、燃料电池、光伏发电及风力发电机等)、储能装置(电动汽车充电站)及负荷构成。
参见图16和图17,其分别为高频充电机的原理图和电路图。其中,uA、uB、uC为来自微网的三相交流电。为了简化计算,用一个非线性电阻RC来近似的模拟高频功率变换电路的等效输入阻抗。所用参数如图17所示,非线性电阻RC可以近似的表示为:
式中,U0为充电机的输出电压,I0为充电机的输出电流,II为高频功率变换器的输入电流,P0为充电机的输出功率,η为充电机效率。
本发明采用Matlab/Simulink软件作为研究工具,以图18为原理图搭建仿真模型如图19所示。对上述仿真模型进行试验,由于充电机电路工作在开关状态,其转换效率高,整个工作期间效率都在90%左右,仿真过程取值为90%。供电电源通过10KV母线给充电站供电,经过10/0.4KV变压器降压后给充电机供电, 充电机将电能经过一系列变换后给蓄电池充电。每台充电机的功率都是9KW。通过运行充电机仿真模型,得到的谐波电流如图20所示。运用本发明所提出的检测方法对A相谐波电流进行检测,得到的谐波电流瞬时频率和瞬时幅值如图21所示。从图21中可知:(1)电动汽车充电机产生的谐波电流次数主要为6k±1次,k=1,2,3,…,即5、7、11、13…;(2)谐波电流的幅值与谐波次数成反比,谐波次数越高,谐波电流幅值越小。
对微网中电动汽车充电站并网点产生的谐波电流进行实时、准确地检测不仅是对其进行有效治理的重要前提,而且是微网稳定可靠运行的重要前提之一。为此,本发明提出了一种微网中电动汽车充电站的谐波检测新方法,即首先运用小波变换的方法对电动汽车充电机的谐波电流进行预处理,初步细化其频率避免了EMD过程中出现的频谱混叠现象;接着运用自适应分步延拓法对谐波信号进行延拓避免了端点效应问题的产生。新方法的提出对于原HHT存在的频谱混叠问题和端点效应问题分别进行了有效合理的改进,取得了优良的效果。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (1)
1.一种微网中电动汽车充电站谐波检测法,其特征是:
运用小波变换的方法对微网中电动汽车充电机的谐波电流进行预处理,将其分解成不同频率段的信号族S1;
用自适应分步延拓法对频率信号族S1内的每一个频率信号进行延拓处理,得到信号族S;
对频率信号族S的每个信号进行HHT分析,在得到的瞬时频率信号和瞬时幅值信号中去掉因两次延拓而增加的信号段,得到最终的信号;
所述自适应分步延拓法步骤包括:
步骤1:首先构造回归模型,在所构造的回归模型的基础上对所要延拓的数据序列进行向左和向右延拓得到信号族S2;
步骤2:在步骤1的基础上对得到的频率信号族根据端点处的函数值大小分情况进行第二次延拓得到信号族S;
所述自适应分步延拓法具体过程包括:
首先构造一个如式(1)所示回归模型,其中,ai为非负的拉格朗日乘子;b为偏置项;k(x,xi)为核函数;l为训练样本数;
然后确定要延拓的点数m和训练样本数l,然后按照一定的规则形成一个训练样本集T={(x1,y1),(x2,y2),…,(xl,yl)},其中xi=[s1t(i),s1t(i+1),…,s1t(N-l+i-1)]T,yi=s1t(N-l+i),1≤i≤l,x={x1,x2…,xl}T;谐波电流经过小波变换预处理后成为窄带信号族S1={s1t},其中一个谐波电流窄带信号为s1t,其数据序列为s1t(1)、s1t(2)…s1t(N),其中,N为采样点数;利用如式(1)所示的回归模型得到数据序列边界外的第一个预测值s1t(N+1);
再将s1t(N+1)作为原始谐波信号序列新的边界点,就得到下一个预测值s1t(N+2),以此类推,根据所要延拓的数据个数得到全部的延拓序列s1t(N+1),…,s1t(N+m);
所述自适应分步延拓法对数据序列进行延拓,利用网格搜索和交叉验证算法确定其参数,对延拓后的信号进行EMD分解;
所述检测法还具体包括:选择‘db5’小波作为小波基函数,对于重构后的高低频系数,采用自适应分步延拓法对其进行延拓处理,得到的EMD分解结果,接着对得到的固有模态分量进行Hilbert变换,得到的瞬时频率和瞬时幅值,使得EMD分解中的频谱混叠问题得到了有效解决。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210289227.3A CN102841251B (zh) | 2012-08-14 | 2012-08-14 | 一种微网中电动汽车充电站谐波检测法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210289227.3A CN102841251B (zh) | 2012-08-14 | 2012-08-14 | 一种微网中电动汽车充电站谐波检测法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102841251A CN102841251A (zh) | 2012-12-26 |
CN102841251B true CN102841251B (zh) | 2015-01-07 |
Family
ID=47368795
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210289227.3A Active CN102841251B (zh) | 2012-08-14 | 2012-08-14 | 一种微网中电动汽车充电站谐波检测法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102841251B (zh) |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103901273B (zh) * | 2012-12-28 | 2017-12-22 | 白晓民 | 电网谐波检测方法和电网谐波检测装置 |
CN103091554B (zh) * | 2013-03-01 | 2015-04-22 | 山东电力集团公司电力科学研究院 | 线路末端同步测量信号多母线微电网稳态谐波分析方法 |
CN103412205B (zh) * | 2013-07-10 | 2016-08-10 | 华北电力大学(保定) | 一种电动汽车充电设备的测试方法 |
CN103941091A (zh) * | 2014-04-25 | 2014-07-23 | 福州大学 | 基于改进emd端点效应的电力系统hht谐波检测方法 |
CN104360153B (zh) * | 2014-10-22 | 2017-02-08 | 江苏省电力公司淮安供电公司 | 电网谐波在线检测与分析方法 |
CN104569651B (zh) * | 2014-11-24 | 2017-04-05 | 中车青岛四方机车车辆股份有限公司 | 动车组高次谐波发射测试系统 |
CN105510711A (zh) * | 2015-12-24 | 2016-04-20 | 合肥工业大学 | 一种改进的经验模态分解的谐波分析法 |
CN106707021A (zh) * | 2016-08-14 | 2017-05-24 | 华北电力大学(保定) | 一种分布式光伏发电集群并网系统的谐振检测方法 |
CN108061821A (zh) * | 2016-11-05 | 2018-05-22 | 南京理工大学 | 一种改进的双馈风力发电系统谐波检测方法 |
CN106771595A (zh) * | 2016-12-09 | 2017-05-31 | 国网上海市电力公司 | 一种快充桩谐波检测装置 |
CN107016217B (zh) * | 2017-04-28 | 2020-08-04 | 茂硕电源科技股份有限公司 | 一种采样滤波电流值的计算方法及其系统 |
CN107122553B (zh) * | 2017-05-02 | 2021-11-05 | 国网山东省电力公司电力科学研究院 | 考虑随机特性的充电机谐波特性分析仿真方法 |
CN109212310A (zh) * | 2017-06-30 | 2019-01-15 | 南京理工大学 | 一种微电网谐波检测方法 |
CN109116109B (zh) * | 2018-08-16 | 2021-08-06 | 中国电力科学研究院有限公司 | 一种间谐波在线监测方法和装置 |
CN110514921B (zh) * | 2019-07-22 | 2022-07-26 | 华南理工大学 | 一种电力电子变换器非平稳信号中非线性现象的识别方法 |
CN110824282A (zh) * | 2019-11-20 | 2020-02-21 | 湖南铁路科技职业技术学院 | 一种高压充电桩电能质量监测系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102169476A (zh) * | 2011-04-14 | 2011-08-31 | 哈尔滨工业大学 | 一种基于灰色理论的希尔伯特-黄变换端点效应抑制方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6990436B1 (en) * | 2003-11-28 | 2006-01-24 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Computing frequency by using generalized zero-crossing applied to intrinsic mode functions |
CN100385438C (zh) * | 2006-05-19 | 2008-04-30 | 长安大学 | 基于相似波形的数字信号端点数据延拓方法 |
US7908103B2 (en) * | 2007-05-21 | 2011-03-15 | Nilanjan Senroy | System and methods for determining masking signals for applying empirical mode decomposition (EMD) and for demodulating intrinsic mode functions obtained from application of EMD |
TWI439874B (zh) * | 2009-01-10 | 2014-06-01 | Ind Tech Res Inst | 結構化之經驗模態分解之訊號處理方法及其裝置 |
-
2012
- 2012-08-14 CN CN201210289227.3A patent/CN102841251B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102169476A (zh) * | 2011-04-14 | 2011-08-31 | 哈尔滨工业大学 | 一种基于灰色理论的希尔伯特-黄变换端点效应抑制方法 |
Non-Patent Citations (3)
Title |
---|
一种自适应的EMD端点延拓方法;邵晨曦;《电子学报》;20071031;第35卷(第10期);正文第1944页第1段,第1945页第3段,第6段,第1946页第4段,第1947页第1-2段 * |
一致经验模分解边缘效应的极值点对称延拓法;舒忠平;《西北工业大学学报》;20061031;第24卷(第5期);全文 * |
徐世燕.基于延拓技术的经验模态分解改进算法.《吉林大学学报(信息科学版)》.2009,第27卷(第4期), * |
Also Published As
Publication number | Publication date |
---|---|
CN102841251A (zh) | 2012-12-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102841251B (zh) | 一种微网中电动汽车充电站谐波检测法 | |
CN103050942B (zh) | 一种基于hht的ct饱和检测方法 | |
CN103760464B (zh) | 基于解析图求解与svm的小电流接地系统故障选线方法 | |
CN103116090B (zh) | 一种基于小波包分析和支持向量机的三相pwm整流器故障诊断方法 | |
CN103532205B (zh) | 一种电动汽车三相充电机谐波模型的建模方法 | |
CN107329044B (zh) | 一种基于电弧暂态分量的配电网单相接地故障选线方法 | |
CN103336265B (zh) | 一种谐波背景下电能表计量误差定量化分析方法 | |
CN106645929A (zh) | 一种改进的电动汽车充电非线性负荷谐波电能计量方法 | |
CN102323480A (zh) | 一种基于Hilbert-Huang变换的电能质量分析方法 | |
CN111308260B (zh) | 一种基于小波神经网络的电能质量监测和电器故障分析系统及其工作方法 | |
CN105759117A (zh) | 一种基于谐波分析综合等值电路的谐波责任量化方法 | |
CN104898019A (zh) | 一种应用于主动配电网的谐波源定位和谐波责任划分方法 | |
CN103323747A (zh) | 一种小电流接地系统的单相接地故障选线方法 | |
CN102401870A (zh) | 基于Hilbert-huang变换和遗传算法的配电网单相接地故障测距方法 | |
CN105868490A (zh) | 模块化多电平变换器多目标特定谐波抑制脉宽调制方法 | |
CN106154040A (zh) | 一种用于计算电网单点等值时等效谐波阻抗方法 | |
CN103746382A (zh) | 有源电力滤波系统 | |
CN104638671A (zh) | 基于引入粒子群算法及相位扰动的孤岛检测方法 | |
CN101726690B (zh) | 电力系统短路全电流各分量分解方法 | |
CN106645948B (zh) | 一种基于小波变换的自适应电网信号基波检测方法 | |
CN110334476B (zh) | 一种电磁暂态仿真方法及系统 | |
CN106021805A (zh) | 一种基于hht分析的电能表动态测试信号拟合方法 | |
Damala et al. | Change detection filter technique-based fault analysis of HVDC transmission line | |
Li et al. | Model establishment and harmonic analysis of electric vehicle charger | |
CN109164294A (zh) | 一种用于dvr中的自适应电压跌落检测方法 |
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 |