CN104391207A - 一种采用基频单向量s变换的电压暂降检测方法 - Google Patents
一种采用基频单向量s变换的电压暂降检测方法 Download PDFInfo
- 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
Links
- 230000009466 transformation Effects 0.000 title claims abstract description 30
- 238000001514 detection method Methods 0.000 title claims abstract description 23
- 238000004364 calculation method Methods 0.000 claims abstract description 55
- 238000000034 method Methods 0.000 claims abstract description 22
- 238000001228 spectrum Methods 0.000 claims description 4
- 238000011897 real-time detection Methods 0.000 claims description 3
- 230000000295 complement effect Effects 0.000 claims 1
- 238000012545 processing Methods 0.000 abstract description 13
- 238000005070 sampling Methods 0.000 abstract description 10
- 230000015572 biosynthetic process Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 3
- 125000004122 cyclic group Chemical group 0.000 description 2
- 238000013480 data collection Methods 0.000 description 1
- 238000003066 decision tree Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
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变换的求解需要遍历输入数据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之间。
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)
| 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)
| 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 | 西南交通大学 | 一种复合电能质量扰动信号的识别方法 |
-
2014
- 2014-12-09 CN CN201410740988.5A patent/CN104391207B/zh active Active
Patent Citations (5)
| 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)
| Title |
|---|
| 全惠敏 等: "电压暂降的改进S变换分析与分类", 《湖南大学学报（自然科学版）》》 * |
| 徐方维 等: "基于改进S变换的电压暂降识别", 《电力系统保护与控制》 * |
| 易吉良 等: "基于广义S变换的短时电能质量扰动信号分类", 《电网技术》 * |
| 易吉良 等: "基于广义S变换的短时电能质量扰动的模糊模式识别", 《电力系统保护与控制》 * |
| 易吉良 等: "采用不完全S变换的电能质量扰动检测方法", 《高电压技术》 * |
| 李立 等: "采用改进不完全S变换估计电能质量扰动参数", 《电工技术学报》 * |
| 汪洋 等: "基于改进S变换的电压暂降起止时刻检测", 《电力系统及其自动化学报》 * |
Cited By (8)
| 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 |
|---|---|---|
| CN104391207B (zh) | 一种采用基频单向量s变换的电压暂降检测方法 | |
| CN103245832B (zh) | 基于快速s变换的谐波时频特性参数估计方法及分析仪 | |
| CN104459560B (zh) | 基于小波降噪和相关向量机的锂电池剩余寿命预测方法 | |
| CN108319981A (zh) | 一种基于密度的时序数据异常检测方法及装置 | |
| CN103941196A (zh) | 锂离子电池荷电状态估计方法 | |
| CN103018555B (zh) | 一种高精度的电力参数软件同步采样方法 | |
| CN108416141A (zh) | 一种线性时变结构模态振型辨识方法 | |
| CN105929201B (zh) | 一种基于细化谱分析的加速度计动态模型结构参数辨识的方法 | |
| CN105137180B (zh) | 基于六项余弦窗四谱线插值的高精度谐波分析方法 | |
| CN104979829B (zh) | 一种基于改进的FastICA算法的谐波电流估计方法 | |
| CN104808035A (zh) | 一种基于广义双曲s变换的电压暂降检测方法 | |
| Akçay | Subspace-based spectrum estimation in frequency-domain by regularized nuclear norm minimization | |
| CN104977515A (zh) | 一种三相共筒gis局部放电监测中放电类型的识别方法 | |
| CN106154037A (zh) | 一种基于校验的同步相量自适应计算方法 | |
| CN104320144B (zh) | 稀疏度自适应信号重构方法 | |
| CN109359822B (zh) | 电子式电压互感器测量状态评估方法及系统 | |
| CN103543026B (zh) | 基于振动传递率函数和支持向量机的结构损伤识别方法 | |
| CN110705041A (zh) | 一种基于easi的线性结构工作模态参数识别方法 | |
| CN104089791A (zh) | 基于振动的改进损伤定位和损伤程度识别方法 | |
| CN108845182B (zh) | 一种应用于夹具参数提取的相位计算方法 | |
| CN106970265B (zh) | 一种采用多时间尺度不完全s变换估计谐波参数的方法 | |
| CN106446405A (zh) | 一种集成电路器件神经网络建模样本选择方法及装置 | |
| CN104820074B (zh) | 一种基于移动窗函数的电子鼻特征提取方法 | |
| CN104062498A (zh) | 对称频谱信号的中心频率的质心估计方法 | |
| CN109270455B (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 | ||
| 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. |