一种电力系统的次/超同步谐波检测方法技术方案

本发明专利技术公开了一种电力系统的次/超同步谐波检测方法，包括：提取并分离工频分量；根据模态混叠条件进行分组；提取含有最高频率谐波分量的分组信号并判断是否发生模态混叠；发生模态混叠时进行频移调制，经验模态分解得到该组所含的各次/超同步谐波对应的固有模态函数；进行幅值相位和频率补偿；希尔伯特变换得到该组次/超同步谐波的瞬时频率和幅值；分离该组次/超同步谐波；重复上述步骤直至分离所有分组。本发明专利技术具有如下优点：基于信号调制的频移方法，通过增大两个次/超同步谐波分量的频率比，扩展经验模态分解的信号分解范围，适用于电力系统次/超同步谐波的实时精确检测。

Sub / Ultra synchronous harmonic detection method for power system

The invention discloses a power system / super synchronous harmonic detection method, including: the extraction and separation of power frequency component; according to the mode mixing conditions of the packet contains the highest frequency harmonic component extraction; packet signal and to determine whether the occurrence of modal aliasing; frequency shift modulation mode mixing, EMD obtained the intrinsic mode function corresponding to the harmonic / super synchronous group contained; amplitude phase and frequency compensation; Hilbert transform to get the instantaneous frequency and amplitude of the second harmonic / super synchronous; the second harmonic separation / supersynchronous; repeat the above steps until the separation of all groups. The invention has the following advantages: signal modulation method based on frequency shift, by increasing two times / super synchronous harmonic frequency ratio, extended signal EMD decomposition range, suitable for power system harmonic / super synchronous real-time accurate detection.

【技术实现步骤摘要】
一种电力系统的次/超同步谐波检测方法
本专利技术涉及电力系统动态监测与电能质量
，具体涉及一种电力系统的次/超同步谐波检测方法。
技术介绍
近年来，风电、光伏等新型可再生能源发电迅速扩展，其渗透率不断提高，由风电机组控制器与交流电网之间的相互作用引起的次同步振荡问题日益凸显，影响电网和设备的安全运行，甚至造成严重的稳定性事故或设备损伤，危害极大。次同步振荡在电网中主要表现为幅值较大的持续增长(发散)或恒幅的次同步和/或超同步电流、电压和功率谐波，准确检测次/超同步谐波的幅值、频率等特征参数，有助于深入分析电力系统次同步振荡特性，从而有效抑制振荡，保障机组及电力系统稳定运行。电力系统次/超同步谐波具有非线性、非平稳性特性，而传统的傅里叶方法不能准确得到振荡模式的时变特性，因此难以有效检测次/超同步谐波。1998年，美国宇航局NordenE.Huang等人提出了针对非线性、非平稳信号的自适应分析方法——希尔伯特·黄变换(Hilbert-HuangTransform，HHT)，HHT分为两步，第一步为经验模态分解(EmpiricalModeDecomposition,EMD)，第二步为希尔伯特变换。其中EMD是HHT的核心部分，将信号中多个组成分量分解成一系列固有模态函数(IntrinsicModeFunction,IMF)，使得每一瞬间每一个IMF只有一个频率，然后经过希尔伯特变换得到有实际物理意义的每一个IMF的瞬时幅值和瞬时频率。但EMD方法应用中存在端点效应、模态混叠等问题，尤其是当信号中两个组成分量的频率在2倍频内时，EMD无法将两者分解开，产生模态混叠现象，降低了频率分辨率。由于次/超同步谐波与基波或者次/超同步谐波之间的频率接近2倍频或在2倍频内，因此这个缺点制约了HHT在次/超同步谐波检测中的应用。
技术实现思路
本专利技术旨在至少解决上述技术问题之一。为此，本专利技术的目的在于提出一种电力系统的次/超同步谐波检测方法，解决模态混叠问题。为了实现上述目的，本专利技术的实施例公开了一种电力系统的次/超同步谐波检测方法，包括以下步骤：S1：假设实测信号包含工频分量和多个次/超同步谐波分量，如果工频分量幅值超过频率与其相近的次/超同步谐波幅值预设倍数时，将所述工频分量从原信号中提取并分离；S2：通过对所述实测信号进行快速傅里叶变换(FastFourierTransform,FFT)，得到所述实测信号所含的各次/超同步谐波的大致频率，将频率比小于频率比阈值的两个谐波分量分为一组，直接进行EMD会发生混叠，将剩余谐波分量根据各自频率所在区间进行分组；从而将所述实测信号所含的次/超同步谐波分量全部进行分组；S3：提取含有所述实测信号未分离部分最高频率谐波分量的分组信号，根据预设混叠设定条件判断所述分组信号是否使用信号频移方法进行调制；S4：如果所述分组信号不需要调制，EMD得到所述分组所含的各次/超同步谐波对应的IMF，如果所述分组信号需要调制，则对经过调制的信号进行EMD，得到所述分组所含的各次/超同步谐波经过调制后对应的IMF；S5：对所有IMF的幅值、相位和频率根据预设补偿条件进行补偿；S6：对补偿后的函数进行希尔伯特变换得到所述分组所含的各次/超同步谐波的瞬时频率和瞬时幅值；S7：将所述分组信号从所述实测信号中分离；S8：重复S3-S7，直至将所述实测信号中全部分组进行分离。根据本专利技术实施例的电力系统的次/超同步谐波检测方法，通过增大两个次/超同步谐波分量的频率比，扩展EMD的信号分解范围，适用于电力系统次/超同步谐波的实时精确检测另外，根据本专利技术上述实施例的电力系统的次/超同步谐波检测方法，还可以具有如下附加的技术特征：进一步地，在步骤S2中，通过以下步骤对所述实测信号所含的次/超同步谐波分量进行分组：将所述实测信号所含的次/超同步谐波分量按照频率的高低进行降序排列，计算两个相邻谐波分量的频率比；将频率比小于预设频率比阈值的两个谐波分量分为一组，对得到的多个分组按照组成分量频率的高低进行降序排列；将频率位于两个相邻分组所含谐波频率之间的剩余谐波分量分为一组。进一步地，在步骤S3中，使用高通滤波器得到含有最高频率谐波分量的分组信号。进一步地，在步骤S4中，在得到所述多个IMF之后还包括：筛除包含混叠频率的IMF。进一步地，所述预设倍数为5倍。本专利技术的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本专利技术的实践了解到。附图说明本专利技术的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：图1是本专利技术一个实施例的电力系统的次/超同步谐波检测方法的流程图。具体实施方式下面详细描述本专利技术的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，仅用于解释本专利技术，而不能理解为对本专利技术的限制。参照下面的描述和附图，将清楚本专利技术的实施例的这些和其他方面。在这些描述和附图中，具体公开了本专利技术的实施例中的一些特定实施方式，来表示实施本专利技术的实施例的原理的一些方式，但是应当理解，本专利技术的实施例的范围不受此限制。相反，本专利技术的实施例包括落入所附加权利要求书的精神和内涵范围内的所有变化、修改和等同物。以下结合附图描述本专利技术。图1是本专利技术一个实施例的电力系统的次/超同步谐波检测方法的流程图。如图1所示，一种电力系统的次/超同步谐波检测方法，包括以下步骤：S1：假设实测信号包含工频分量和多个次/超同步谐波分量，如果工频分量幅值超过频率与其相近的次/超同步谐波幅值预设倍数时，将所述工频分量从原信号中提取并分离。具体地，实测信号s(t)中，如果工频分量幅值是频率与其相近的次/超同步谐波幅值的5倍以上，为了避免工频分量对该次/超同步谐波的特征参数提取产生干扰，将工频分量从原信号中提取并分离，具体过程如下：首先使用带通滤波器获得所述工频分量，经过EMD得到相应的IMF，记为g(t)，进行幅值相位补偿后，通过希尔伯特变换得到所述工频分量的瞬时幅值和瞬时频率，再将所述工频分量从实测信号中分离。S2：通过对所述分离出工频分量的实测信号进行FFT，得到所述实测信号所含的各次/超同步谐波的大致频率，将频率比小于频率比阈值的两个谐波分量分为一组，直接进行EMD会发生混叠，将剩余谐波分量根据各自频率所在区间进行分组；从而将所述实测信号所含的次/超同步谐波分量全部进行分组。具体地，假设s(t)中含有N个次/超同步谐波分量，根据预设的频率分辨率，选择相应的数据窗长对s(t)进行FFT计算，得到N个次/超同步谐波分量的大致频率值fn(n＝1,2,…,N),其中f1>f2>f3>…>fN，将N个次/超同步谐波分量按照频率的高低进行降序排列，计算两个相邻谐波分量的频率比mn,mn＝fn/fn+1(n＝1,2,…,N-1)；根据混叠现象的产生条件，将所述实测信号所含的次/超同步谐波分量进行分组，具体过程如下：将mn(1≤n≤N-1)依次与2进行比较，若mn≤2，则将频率为fn和fn+1的组成分量分为一组，若此时mn+1≤2同时成立，由于每个次/超同步谐波分量都唯一存在于某一分组中本文档来自技高网...

【技术保护点】
一种电力系统的次/超同步谐波检测方法，其特征在于，包括以下步骤：S1：假设实测信号包含工频分量和多个次/超同步谐波分量，如果工频分量幅值超过频率与其相近的次/超同步谐波幅值预设倍数时，将所述工频分量从原信号中提取并分离；S2：通过对所述实测信号进行快速傅里叶变换，得到所述实测信号所含的各次/超同步谐波的大致频率，将频率比小于频率比阈值的两个谐波分量分为一组，直接进行经验模态分解会发生混叠，将剩余谐波分量根据各自频率所在区间进行分组；从而将所述实测信号所含的次/超同步谐波分量全部进行分组；S3：提取含有所述实测信号未分离部分最高频率谐波分量的分组信号，根据预设混叠设定条件判断所述分组信号是否使用信号频移方法进行调制；S4：如果所述分组信号不需要调制，经验模态分解得到所述分组所含的各次/超同步谐波对应的IMF，如果所述分组信号需要调制，则对经过调制的信号进行经验模态分解，得到所述分组所含的各次/超同步谐波经过调制后对应的固有模态函数；S5：对所有固有模态函数的幅值、相位和频率根据预设补偿条件进行补偿；S6：对补偿后的函数进行希尔伯特变换得到所述分组所含的各次/超同步谐波的瞬时频率和瞬时幅值；S7：将所述分组信号从所述实测信号中分离；S8：重复S3‑S7，直至将所述实测信号中全部分组进行分离。...

【技术特征摘要】
1.一种电力系统的次/超同步谐波检测方法，其特征在于，包括以下步骤：S1：假设实测信号包含工频分量和多个次/超同步谐波分量，如果工频分量幅值超过频率与其相近的次/超同步谐波幅值预设倍数时，将所述工频分量从原信号中提取并分离；S2：通过对所述实测信号进行快速傅里叶变换，得到所述实测信号所含的各次/超同步谐波的大致频率，将频率比小于频率比阈值的两个谐波分量分为一组，直接进行经验模态分解会发生混叠，将剩余谐波分量根据各自频率所在区间进行分组；从而将所述实测信号所含的次/超同步谐波分量全部进行分组；S3：提取含有所述实测信号未分离部分最高频率谐波分量的分组信号，根据预设混叠设定条件判断所述分组信号是否使用信号频移方法进行调制；S4：如果所述分组信号不需要调制，经验模态分解得到所述分组所含的各次/超同步谐波对应的IMF，如果所述分组信号需要调制，则对经过调制的信号进行经验模态分解，得到所述分组所含的各次/超同步谐波经过调制后对应的固有模态函数；S5：对所有固有模态函数的幅值、相位和频率根据预设补偿条件进行补偿；S6：对补偿后的函数进行希尔伯...

【专利技术属性】
技术研发人员：谢小荣，贺静波，宋瑞华，张锋，郭小龙，陈军，占颖，
申请(专利权)人：清华大学，国家电网公司，中国电力科学研究院，国网新疆电力公司，
类型：发明
国别省市：北京,11