CN104391207A - 一种采用基频单向量s变换的电压暂降检测方法 - Google Patents

一种采用基频单向量s变换的电压暂降检测方法 Download PDF

Info

Publication number
CN104391207A
CN104391207A CN201410740988.5A CN201410740988A CN104391207A CN 104391207 A CN104391207 A CN 104391207A CN 201410740988 A CN201410740988 A CN 201410740988A CN 104391207 A CN104391207 A CN 104391207A
Authority
CN
China
Prior art keywords
fundamental frequency
transformation
module
sequence
calculation
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
Application number
CN201410740988.5A
Other languages
English (en)
Other versions
CN104391207B (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.)
Foshan jintexin Electronic Technology Co.,Ltd.
Original Assignee
Hunan University of Technology
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 Hunan University of Technology filed Critical Hunan University of Technology
Priority to CN201410740988.5A priority Critical patent/CN104391207B/zh
Publication of CN104391207A publication Critical patent/CN104391207A/zh
Application granted granted Critical
Publication of CN104391207B publication Critical patent/CN104391207B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measurement Of Current Or Voltage (AREA)
  • Measuring Frequencies, Analyzing Spectra (AREA)

Abstract

本发明公开了一种采用基频单向量S变换的电压暂降检测方法，该方法包括数据采集及计数模块(1)、计算序列形成模块(2)、基频单向量S变换模块(3)、电压暂降识别及参数计算模块(4)等4个功能模块；数据采集及计数模块(1)根据设定的采样数据量预设值将暂存于数据缓冲区的数据提取至计算序列形成模块(2)，计算序列形成模块(2)根据算法处理完成标志组成计算序列，基频单向量S变换模块(3)对计算序列做基频单向量S变换并求其模向量，电压暂降识别及参数计算模块(4)根据基频单向量S变换模向量幅值判断是否有暂降发生，并计算和保存暂降幅值和起止时刻；本发明为电压暂降的检测提供了一种实时精确的手段。

Description

一种采用基频单向量S变换的电压暂降检测方法
技术领域
本发明涉及短时电能质量扰动检测领域，特别涉及一种电压暂降检测方法。
背景技术
电压暂降是最常见和危害很大的电能质量问题，快速检测电压暂降并准确估计其参数，是治理该类电能质量的前提，对保障敏感设备的安全可靠运行具有重要意义。
在分析短时电能质量扰动时常采用时频分析工具，近年S变换由于其优良的时频特性，已成为电能质量扰动分析的热点工具；在利用S变换分析电能质量扰动过程中，通常利用S变换求模时频矩阵，并从中提取基频幅值向量、时间幅值包络、频率幅值包络等特征向量，然后采用决策树等分类器根据特征向量的特征值对扰动进行识别或计算扰动参数。
S变换在求解扰动数据过程中采用固定窗宽的高斯函数，其在低频部分(包括基频)求解后的时间分辨率较低，难以精确检测持续时间仅有几个基频周波的电压暂降扰动；通过增加高斯函数窗宽调整系数可以提高S变换低频部分的时间分辨率，使其精确检测短时电压暂降幅值和起止时刻成为可能。
S变换的求解需要遍历输入数据FFT频谱的所有频率点，而在检测和计算电压暂降参数时只需用到基频向量，所以在此类应用中，如果在求解数据的S变换时仅计算基频向量，可以大大缩短运算时间，保障算法的实时应用。
发明内容
本发明提供一种电压暂降检测和参数计算方法，通过只求解改进S变换的基频向量，提高电压暂降检测和参数计算的准确性和实时性。
本发明所述的采用基频单向量S变换的电压暂降检测方法，包括数据采集及计数模块(1)、计算序列形成模块(2)、基频单向量S变换模块(3)、电压暂降识别及参数计算模块(4)，所述4个功能模块依次循环联系实现电压暂降的实时检测，数据采集及计数模块(1)根据设定的采样数据量计数器预设值将暂存于数据缓冲区的数据截取至计算序列形成模块(2)，计算序列形成模块(2)根据算法处理完成标志P组成计算序列，基频单向量S变换模块(3)对计算序列做基频单向量S变换并求其模向量，电压暂降识别及参数计算模块(4)根据基频单向量S变换模向量幅值判断是否有暂降发生，若无暂降则以20个基频周波为处理周期循环检测是否有暂降发生，若有暂降且不存在暂降结束时刻，以0.5个基频周波为处理周期循环检测计数序列直到有暂降结束时刻为止；若有暂降且存在暂降结束时刻，则计算并保存暂降幅值和起止时刻，并保存计算序列数据。
所述采用基频单向量S变换的电压暂降检测方法，其初始化过程中采样数据量计数器预设值为20个基频周波，初始化过程中处理完成标志P的预设值为1。
所述采用基频单向量S变换的电压暂降检测方法，其数据采集及计数模块(1)通过数据过零点检测采样数据量的周波数，并触发截取数据事件，将数据从缓冲区截取到计算区。
所述采用基频单向量S变换电压暂降检测方法，其计算序列形成模块(2)根据处理完成标志P组成计算序列h(k)，若P=1，则计算序列h(k)只包括新截取的数据；若P=0，则由新截取的数据续接上次计算序列组成新的计算序列h(k)。
所述采用基频单向量S变换的电压暂降检测方法，其基频单向量S变换模块(3)的计算步骤如下：
S1：检测计算序列h(k)的数据长度l，进入步骤S2；
S2：判断l是否等于20个基频周波的数据长度，若相等则进入步骤S4，若不相等则进入步骤S3；
S3：将20个基频周波长度的高斯函数序列g(m)两端补零使其长度为l，进入步骤S4；
S4：对h(k)求FFT频谱序列Ds(m)，进入S5；
S5：将高斯函数序列g(m)中心循环移位到对齐Ds(m)基频频率点，进入步骤S6；
S6：将高斯函数序列g(m)与Ds(m)做数组乘（加窗处理），得序列Dg(m)，进入步骤S7；
S7：对序列Dg(m)求FFT反变换得基频单向量S变换结果，其为复数向量Vc(k)，进入步骤S8；
S8：对Vc(k)求模得模向量Vm(k)，结果送入电压暂降及参数计算模块(4)。
所述的基频单向量S变换模块(3)中高斯序列g(m)表达式为g(m)=exp(-wm 2)，其中w为窗宽调整系数，其取值为0.0037到0.0075之间。
所述采用基频单向量S变换电压暂降检测方法，其电压暂降识别及参数计算模块(4)通过检测模向量Vm(k)的最小值来判断电压暂降的发生，即Vm(k)的最小值小于0.45，则有电压暂降发生，否则未发生电压暂降扰动。
所述电压暂降识别及参数计算模块(4)通过查找模向量Vm(k)的差分向量Vd(k)的有效极大值点来判断电压暂降的结束时刻，所述有效极大值点是差分向量Vd(k)幅值大于0.0005的极大值对应的时间采样点。
附图说明
图1为本发明所述电压暂降检测方法的功能模块及流程图；
图2为本发明所述基频单向量S变换求解流程；
图3为本发明的实施举例说明。
具体实施方式
为解决电压暂降的精确检测和实时应用问题，本发明提供一种采用基频单向量S变换的电压暂降检测方法，以下结合附图对本发明的优选实施实例进行阐述，需要说明的是，优选实施实例是为了进一步说明本发明，而不是为了限制本发明的保护范围。
如图1为本发明所述一种采用基频单向量S变换的电压暂降检测方法的功能模块及流程图。
所述采用基频单向量S变换的电压暂降检测方法，包括数据采集及计数模块(1)、计算序列形成模块(2)、基频单向量S变换模块(3)、电压暂降识别及参数计算模块(4)，所述4个功能模块依次循环联系实现电压暂降的实时检测，数据采集及计数模块(1)根据设定的采样数据量计数器预设值将暂存于数据缓冲区的数据截取至计算序列形成模块(2)，计算序列形成模块(2)根据算法处理完成标志P组成计算序列，基频单向量S变换模块(3)对计算序列做基频单向量S变换并求其模向量，电压暂降识别及参数计算模块(4)根据基频单向量S变换模向量幅值判断是否有暂降发生，若无暂降则以20个基频周波为处理周期循环检测是否有暂降发生，若有暂降且不存在暂降结束时刻，以0.5个基频周波为处理周期循环检测计数序列直到有暂降结束时刻为止；若有暂降且存在暂降结束时刻，则计算并保存暂降幅值和起止时刻，并保存计算序列数据。
所述采用基频单向量S变换的电压暂降检测方法，其初始化过程中采样数据量计数器预设值为20个基频周波，初始化过程中处理完成标志P的预设值为1。
所述采用基频单向量S变换的电压暂降检测方法，其数据采集及计数模块(1)通过检测数据过零点确定采样数据量的周波数，并据此触发截取数据事件，将数据从缓冲区截取到计算区。
所述采用基频单向量S变换的电压暂降检测方法，其计算序列形成模块(2)根据处理完成标志P组成计算序列h(k)，若P=1，则计算序列h(k)只包括新截取的数据；若P=0，则由新截取的数据续接上次计算序列组成新的计算序列h(k)。
如图2为所述基频单向量S变换模块(3)求解流程，其计算步骤如下：
S1：检测计算序列h(k)的数据长度l，进入步骤S2；
S2：判断l是否等于20个基频周波的数据长度，若相等则进入步骤S4，若不相等则进入步骤S3；
S3：将20个基频周波长度的高斯函数序列g(m)两端补零使其长度为l，进入步骤S4；
S4：对h(k)求FFT频谱序列Ds(m)，进入S5；
S5：将高斯函数序列g(m)中心循环移位到对齐Ds(m)基频频率点，进入步骤S6；
S6：将高斯函数序列g(m)与Ds(m)做数组乘（加窗处理），得序列Dg(m)，进入步骤S7；
S7：对序列Dg(m)求FFT反变换得基频单向量S变换结果，其为复数向量Vc(k)，进入步骤S8；
S8：对Vc(k)求模得模向量Vm(k)，结果送入电压暂降及参数计算模块(4)。
所述的基频单向量S变换模块(3)中高斯序列g(m)表达式为g(m)=exp(-wm 2)，其中w为窗宽调整系数，其取值为0.0037到0.0075之间。
所述采用基频单向量S变换电压暂降检测方法，其电压暂降识别及参数计算模块(4)通过检测模向量Vm(k)的最小值来判断电压暂降的发生，即Vm(k)的最小值小于0.45，则有电压暂降发生，否则未发生电压暂降扰动。
所述电压暂降识别及参数计算模块(4)通过查找模向量Vm(k)的差分向量Vd(k)的有效极大值点来判断电压暂降的结束时刻，所述有效极大值点是差分向量Vd(k)幅值大于0.0005的极大值对应的时间采样点。
如图3为本发明的实施举例说明，其中，图3(a)为电压跌落原始信号，电压暂降幅值由图3(b)所示跌落幅值计算区间的模向量Vm(k)幅值的均值确定，电压暂降的开始时刻由差分向量Vd(k)的有效极小值确定，如图3(c)所示，所述有效极小值是差分向量Vd(k)幅值小于-0.0005的极小值对应的时间采样点。 
以上所述的本发明的实施方式，并非成为本发明保护范围的限定，倘若对本发明实施方式进行各种变形或修改，但尚在本发明的精神和原则之内，均应包含在本发明的权利要求保护范围之内。

Claims (4)

1.一种采用基频单向量S变换的电压暂降检测方法，包括数据采集及计数模块(1)、计算序列形成模块(2)、基频单向量S变换模块(3)、电压暂降识别及参数计算模块(4)，所述4个功能模块依次循环联系实现电压暂降的实时检测，其特征在于，数据采集及计数模块(1)根据设定的采样数据量计数器预设值将暂存于数据缓冲区的数据截取至计算序列形成模块(2)，计算序列形成模块(2)根据算法处理完成标志P组成计算序列，基频单向量S变换模块(3)对计算序列做基频单向量S变换并求其模向量，电压暂降识别及参数计算模块(4)根据基频单向量S变换模向量幅值判断是否有暂降发生，若无暂降则以20个基频周波为处理周期循环检测是否有暂降发生，若有暂降且不存在暂降结束时刻，以0.5个基频周波为处理周期循环检测计数序列直到有暂降结束时刻为止；若有暂降且存在暂降结束时刻，则计算并保存暂降幅值和起止时刻，并保存计算序列数据。
2.根据权利要求1所述的一种采用基频单向量S变换电压暂降检测方法，其特征在于，所述方法的初始化过程中采样数据量计数器预设值为20个基频周波，初始化过程中处理完成标志P的预设值为1。
3.根据权利要求1所述的一种采用基频单向量S变换电压暂降检测方法，其特征在于，所述基频单向量S变换模块(3)的计算步骤如下：
S1：检测计算序列h(k)的数据长度l，进入步骤S2；
S2：判断l是否等于20个基频周波的数据长度，若相等则进入步骤S4，若不相等则进入步骤S3；
S3：将20个基频周波长度的高斯函数序列g(m)两端补零使其长度为l，进入步骤S4；
S4：对h(k)求FFT频谱序列Ds(m)，进入S5；
S5：将高斯函数序列g(m)中心循环移位到对齐Ds(m)基频频率点，进入步骤S6；
S6：将高斯函数序列g(m)与Ds(m)做数组乘（加窗处理），得序列Dg(m)，进入步骤S7；
S7：对序列Dg(m)求FFT反变换得基频单向量S变换结果，其为复数向量Vc(k)，进入步骤S8；
S8：对Vc(k)求模得模向量Vm(k)，结果送入电压暂降及参数计算模块(4)。
4.根据权利要求2所述的基频单向量S变换模块(3)，所述高斯序列g(m)表达式为g(m)=exp(-wm 2)，其中w为窗宽调整系数，其特征在于，w窗宽调整系数取值为0.0037到0.0075之间。
CN201410740988.5A 2014-12-09 2014-12-09 一种采用基频单向量s变换的电压暂降检测方法 Active CN104391207B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410740988.5A CN104391207B (zh) 2014-12-09 2014-12-09 一种采用基频单向量s变换的电压暂降检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410740988.5A CN104391207B (zh) 2014-12-09 2014-12-09 一种采用基频单向量s变换的电压暂降检测方法

Publications (2)

Publication Number Publication Date
CN104391207A true CN104391207A (zh) 2015-03-04
CN104391207B CN104391207B (zh) 2018-03-02

Family

ID=52609132

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410740988.5A Active CN104391207B (zh) 2014-12-09 2014-12-09 一种采用基频单向量s变换的电压暂降检测方法

Country Status (1)

Country Link
CN (1) CN104391207B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104749432A (zh) * 2015-03-12 2015-07-01 西安电子科技大学 基于聚焦s变换的多分量非平稳信号瞬时频率估计方法
CN108983046A (zh) * 2018-08-16 2018-12-11 国网山东省电力公司泰安供电公司 一种基于奇异值分解法的电压暂降状态估计方法及系统
CN109188081A (zh) * 2018-09-12 2019-01-11 四川大学 一种基于空间向量法的电压暂降波形点检测方法
CN111257619A (zh) * 2020-02-17 2020-06-09 南京工程学院 一种基于多属性决策和改进s变换方法的电压暂降检测方法
CN112540220A (zh) * 2020-11-05 2021-03-23 广东电网有限责任公司广州供电局 电压暂降检测电路和装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02161366A (ja) * 1988-12-14 1990-06-21 Kansai Electric Power Co Inc:The 変電所の事故点判別方法
CN102608441A (zh) * 2011-01-25 2012-07-25 华北电力科学研究院有限责任公司 基于s变换的sdpqd信号识别方法、设备及系统
CN102831433A (zh) * 2012-06-06 2012-12-19 西南交通大学 一种基于时频域多特征量的电能质量混合扰动分类方法
CN103995178A (zh) * 2014-05-20 2014-08-20 江苏大学 一种基于时频聚集特性准则s变换的电压暂降检测方法
CN104034974A (zh) * 2014-05-09 2014-09-10 西南交通大学 一种复合电能质量扰动信号的识别方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02161366A (ja) * 1988-12-14 1990-06-21 Kansai Electric Power Co Inc:The 変電所の事故点判別方法
CN102608441A (zh) * 2011-01-25 2012-07-25 华北电力科学研究院有限责任公司 基于s变换的sdpqd信号识别方法、设备及系统
CN102831433A (zh) * 2012-06-06 2012-12-19 西南交通大学 一种基于时频域多特征量的电能质量混合扰动分类方法
CN104034974A (zh) * 2014-05-09 2014-09-10 西南交通大学 一种复合电能质量扰动信号的识别方法
CN103995178A (zh) * 2014-05-20 2014-08-20 江苏大学 一种基于时频聚集特性准则s变换的电压暂降检测方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
全惠敏 等: "电压暂降的改进S变换分析与分类", 《湖南大学学报（自然科学版）》》 *
徐方维 等: "基于改进S变换的电压暂降识别", 《电力系统保护与控制》 *
易吉良 等: "基于广义S变换的短时电能质量扰动信号分类", 《电网技术》 *
易吉良 等: "基于广义S变换的短时电能质量扰动的模糊模式识别", 《电力系统保护与控制》 *
易吉良 等: "采用不完全S变换的电能质量扰动检测方法", 《高电压技术》 *
李立 等: "采用改进不完全S变换估计电能质量扰动参数", 《电工技术学报》 *
汪洋 等: "基于改进S变换的电压暂降起止时刻检测", 《电力系统及其自动化学报》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104749432A (zh) * 2015-03-12 2015-07-01 西安电子科技大学 基于聚焦s变换的多分量非平稳信号瞬时频率估计方法
CN104749432B (zh) * 2015-03-12 2017-06-16 西安电子科技大学 基于聚焦s变换的多分量非平稳信号瞬时频率估计方法
CN108983046A (zh) * 2018-08-16 2018-12-11 国网山东省电力公司泰安供电公司 一种基于奇异值分解法的电压暂降状态估计方法及系统
CN109188081A (zh) * 2018-09-12 2019-01-11 四川大学 一种基于空间向量法的电压暂降波形点检测方法
CN109188081B (zh) * 2018-09-12 2019-08-20 四川大学 一种基于空间向量法的电压暂降波形点检测方法
CN111257619A (zh) * 2020-02-17 2020-06-09 南京工程学院 一种基于多属性决策和改进s变换方法的电压暂降检测方法
CN111257619B (zh) * 2020-02-17 2022-06-14 南京工程学院 一种基于多属性决策和改进s变换方法的电压暂降检测方法
CN112540220A (zh) * 2020-11-05 2021-03-23 广东电网有限责任公司广州供电局 电压暂降检测电路和装置

Also Published As

Publication number Publication date
CN104391207B (zh) 2018-03-02

Similar Documents

Publication Publication Date Title
CN104391207A (zh) 一种采用基频单向量s变换的电压暂降检测方法
CN103245832B (zh) 基于快速s变换的谐波时频特性参数估计方法及分析仪
US9563530B2 (en) Device state estimation apparatus, device power consumption estimation apparatus, and program
CN109375060B (zh) 一种配电网故障波形相似度计算方法
Reddy et al. Fast tracking of power quality disturbance signals using an optimized unscented filter
CN105137180B (zh) 基于六项余弦窗四谱线插值的高精度谐波分析方法
CN105954695B (zh) 一种基于同步的同质传感器突变参数识别方法与装置
Jin et al. A robust and adaptive detection scheme for interharmonics in active distribution network
CN106137184B (zh) 基于小波变换的心电信号qrs波检测方法
CN104979829B (zh) 一种基于改进的FastICA算法的谐波电流估计方法
CN103809023A (zh) 基于二分搜索的电网同步谐波相量测量方法
CN106487026B (zh) 一种基于广域测量信息的同调机群快速识别方法
CN108281961B (zh) 一种自适应鲁棒扩展卡尔曼的参数辨识方法
Mendonça et al. Least squares optimization of zero crossing technique for frequency estimation of power system grid distorted sinusoidal signals
CN106842099A (zh) 一种电流互感器饱和识别方法及装置
Zhu et al. Injection amplitude guidance for impedance measurement in power systems
CN104330644B (zh) 检测电力系统中正弦波信号幅值的方法和系统
CN106970265B (zh) 一种采用多时间尺度不完全s变换估计谐波参数的方法
CN108845182B (zh) 一种应用于夹具参数提取的相位计算方法
KR101207495B1 (ko) 전기품질 분석을 위한 시간영역에서 왜곡파형의 고조파 요소 측정방법
CN108008262B (zh) 一种局部放电快速检测方法
CN112051479A (zh) 一种配电网运行状态识别方法及系统
CN104808055A (zh) 一种电信号频率的数字化测量方法
de Souza et al. A modified time-frequency method for testing wide-sense stationarity
CN105078442A (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
TR01 Transfer of patent right

Effective date of registration: 20200904

Address after: 233000 Xianghe Home Floor No.9, Xinhuai Road, Dongsheng Street, Longzihu District, Bengbu City, Anhui Province

Patentee after: Bengbu Qibang Science and Technology Information Consulting Co.,Ltd.

Address before: 412007 Zhuzhou Research Institute, Hunan University of Technology, 88 West Taishan Road, Tianyuan District, Hunan

Patentee before: HUNAN University OF TECHNOLOGY

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20220402

Address after: 528300 one of the fourth floor, No. 10, South 1st Road, guanglongwei Industrial Zone, Yongxing community, Chencun Town, Shunde District, Foshan City, Guangdong Province (residence declaration)

Patentee after: Foshan jintexin Electronic Technology Co.,Ltd.

Address before: 233000 No.9, 1st floor, Xianghe Jiayuan, Xinhuai Road, Dongsheng Street, Longzihu District, Bengbu City, Anhui Province

Patentee before: Bengbu Qibang Science and Technology Information Consulting Co.,Ltd.

TR01 Transfer of patent right