CN105976266A - 基于典型电压暂降波形匹配的电压暂降原因识别方法 - Google Patents
基于典型电压暂降波形匹配的电压暂降原因识别方法 Download PDFInfo
- Publication number
- CN105976266A CN105976266A CN201610324282.XA CN201610324282A CN105976266A CN 105976266 A CN105976266 A CN 105976266A CN 201610324282 A CN201610324282 A CN 201610324282A CN 105976266 A CN105976266 A CN 105976266A
- Authority
- CN
- China
- Prior art keywords
- waveform
- temporarily
- voltage
- exemplary voltages
- voltage dip
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000012544 monitoring process Methods 0.000 claims abstract description 20
- 239000000284 extract Substances 0.000 claims description 28
- 238000012360 testing method Methods 0.000 claims description 26
- 239000011159 matrix material Substances 0.000 claims description 15
- 238000005452 bending Methods 0.000 claims description 13
- 230000002085 persistent effect Effects 0.000 claims description 6
- 230000006698 induction Effects 0.000 claims description 4
- 239000000203 mixture Substances 0.000 claims description 4
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000008878 coupling Effects 0.000 claims description 3
- 238000010168 coupling process Methods 0.000 claims description 3
- 238000005859 coupling reaction Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 230000001932 seasonal effect Effects 0.000 claims description 3
- 230000002265 prevention Effects 0.000 abstract description 5
- 230000005611 electricity Effects 0.000 description 6
- 230000006872 improvement Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 2
- 230000007935 neutral effect Effects 0.000 description 2
- 238000009412 basement excavation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- RFHAOTPXVQNOHP-UHFFFAOYSA-N fluconazole Chemical compound C1=NC=NN1CC(C=1C(=CC(F)=CC=1)F)(O)CN1C=NC=N1 RFHAOTPXVQNOHP-UHFFFAOYSA-N 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001537 neural effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/22—Matching criteria, e.g. proximity measures
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Water Supply & Treatment (AREA)
- Public Health (AREA)
- Primary Health Care (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种基于典型电压暂降波形匹配的电压暂降原因识别方法,包括以下步骤,步骤(A),采集典型电压暂降历史事件的历史监测波形;步骤(B),建立形成典型电压暂降事件波形库;步骤(C),提取待识别的电压暂降波形;步骤(D),对提取待识别的电压暂降波形的波段与典型电压暂降波形库进行波形匹配分析;步骤(E),识别当前电压暂降事件波形产生的原因。本发明在波形的层面上对电压暂降进行匹配分析,更为基础也更为准确,能够快速准确的识别电压暂降产生的原因,对于电压暂降的预防和治理具有十分重要的参考价值,具有良好的应用前景。
Description
技术领域
本发明涉及一种基于典型电压暂降波形匹配的电压暂降原因识别方法,属于电能质量分析技术领域。
背景技术
电压暂降,又称为电压骤降、电压凹陷或电压跌落,电压暂降是供电电压方均根值在短时间内突然下降后又迅速恢复到正常状态的事件,一般典型时间为0.5~30个周波的时间。
随着经济的快速发展,敏感电力电子设备在工业中的广泛应用,电压暂降造成的损失也显著增多。根据国外某电力公司的调查,关于电压暂降的投诉占电能质量投诉的80%,电压暂降已经成为最重要的电能质量问题。
引起电压暂降的主要原因分为三大类,分别为:电力系统短路故障、感应电机启动、变压器投运,其中,电力系统短路故障又根据实际情况分为四种,分别为单相接地短路故障、相间短路故障、两相接地短路故障以及三相接地短路故障,其中不同原因引起的电压暂降有不同的特征,在电压暂降幅值、持续时间、电压变化趋势、相位跳变以及三相暂降不平衡度等有着很大的差异。如何根据不同暂降原因引起的电压暂降特征的不同识别电压暂降产生原因,对于电压暂降的预防和治理措施的制定具有十分重要的意义,同时也是协调用户和电力企业之间纠纷的重要依据。
目前,对电压暂降原因的识别已经有了一定的研究,常见的电压暂降原因识别方法主要有以下两步:首先提取电压暂降的主要特征量如电压暂降幅值、持续时间、相位跳变、三相不平衡度、下降斜率以及上升斜率等,然后通过将这些主要特征量与典型常见类型的电压暂降特征量进行相似度的分析来进行电压暂降原因的识别。目前,各个方法的区别主要存在于电压暂降特征量的提取以及与典型事件波形的相似度分析上,其中,电压暂降特征量的提取方法主要有傅立叶变换、S变换、小波变换等方法,匹配识别的方法主要包括神经网络、模糊推理、支持向量机、ID3算法等。然而,上述方法中,对于电压暂降特征量的提取,一般提取的特征量不超过10个,并且一些在实际情况中较为重要的特征量不能被很好的发掘出来,从而影响识别结果的准确性;而采用模糊推理等方法需要设定一定的阈值,在实际情况中,电力系统的短路容量、变压器的接线方式以及负荷类型等都会在一定程度上影响电压暂降的特征,因此,很难确定固定的阈值使得识别结果十分准确;而且,神经网络需要大量的训练样本,而在实际中实测数据不够丰富,从而影响神经网络方法的应用。
发明内容
本发明的目的是为了克服现有技术中对电压暂降原因识别方法的不足。本发明提供的基于典型电压暂降波形匹配的电压暂降原因识别方法,建立典型电压暂降事件波形库,在波形的层面上对电压暂降进行匹配分析,更为基础也更为准确,能够快速准确的识别电压暂降产生的原因,对于电压暂降的预防和治理具有十分重要的参考价值,具有良好的应用前景。
为了达到上述目的,本发明所采用的技术方案是:
基于典型电压暂降波形匹配的电压暂降原因识别方法,其特征在于:包括以下步骤,
步骤(A),采集典型电压暂降历史事件的历史监测波形
通过暂降录波仪采集典型电压暂降历史事件的历史监测波形,并提取历史监测波形所对应的电压暂降产生原因;
步骤(B),建立形成典型电压暂降事件波形库
将步骤(A)中采集到的历史监测波形和对应的电压暂降产生原因,保存到典型电压暂降事件波形库中,建立电压暂降波形特征量和产生原因的对应关系表,形成典型电压暂降事件波形库;
步骤(C),提取待识别的电压暂降波形
将采集检测到的当前电压暂降事件波形作为测试波形;
步骤(D),对提取待识别的电压暂降波形的波段与典型电压暂降波形库进行波形匹配分析
基于动态弯曲距离相似度算法,对步骤(C)中提取待识别的电压暂降波形与典型电压暂降波形库中的典型电压暂降事件波形进行逐一匹配分析,得到一组匹配度系数;
步骤(E),识别当前电压暂降事件波形产生的原因
对步骤(D)中得到的匹配度系数按大小进行排序,选取匹配度最大的典型电压暂降事件波形作为最佳的匹配波形,提取典型电压暂降事件波形对应的电压暂降所产生的原因,为当前电压暂降事件波形产生原因的识别结果。
前述的基于典型电压暂降波形匹配的电压暂降原因识别方法,其特征在于:步骤(A)提取历史监测波形所对应的电压暂降产生原因,包括以下6种原因:三相短路故障、单相接地故障、两相接地故障、相间短路故障、感应电机启动、变压器投运。
前述的基于典型电压暂降波形匹配的电压暂降原因识别方法,其特征在于:步骤(A)为保证所采集到的典型电压暂降历史事件的历史监测波形具有典型性,每种典型电压暂降历史事件至少采集10个波形,并选取其中最为典型的波形作为历史监测波形,选取过程,包括以下步骤,
(A1)计算某一电压暂降产生的原因所有电压暂降波形典型特征量X={x1,x2,x3,x4,x5,x6,x7},其中,x1~x7为列向量,每个列向量表示电压暂降产生的原因中所有电压暂降波形7个典型特征量之一,包括电压暂降平均幅值、最小幅值、持续时间、电压下降斜率、电压恢复斜率、电压相角跳变、三相幅值不平衡度,每个向量的维数表示该产生原因对应的电压暂降中所采集到的波形个数;
(A2)求取所有波形的每个特征值平均值得到向量
(A3)逐一计算所有波形的典型特征量向量与平均值向量的欧氏距离;
(A4)取所有波形中的欧氏距离最小值,并提取欧氏距离最小值对应的波形为最为典型的波形。
前述的基于典型电压暂降波形匹配的电压暂降原因识别方法,其特征在于:步骤(C),提取待识别的电压暂降波形,包括以下步骤,
(C1)在检测到原因未明的当前电压暂降事件并且计划对该电压暂降进行电压暂降原因识别时,记录该电压暂降事件发生时的三相波形;
(C2)提取三相波形其中发生电压暂降的部分波段,作为待识别的电压暂降波形。
前述的基于典型电压暂降波形匹配的电压暂降原因识别方法,其特征在于:步骤(D),基于动态弯曲距离相似度算法,对步骤(C)中提取待识别的电压暂降波形与典型电压暂降波形库中的典型电压暂降事件波形进行逐一匹配分析,得到一组匹配度系数,包括以下步骤,
(D1)设置典型电压暂降事件波形的时间窗长度为Di1个周波,提取各典型电压暂降事件波形每个时间窗的特征量,形成每个典型电压暂降事件波形的时间序列形特征量集合称为样本波形时间序列特征量集合,其中,i=1,2,…,6,表示典型电压暂降事件波形库中共有6种典型波形,mi为第i种典型电压暂降事件波形总的时间窗个数,rik为第i个典型电压暂降事件波形第k个时间窗的特征量,k=1,2,…,mi,特征量rik为以下6个量:最大值、最小值、平均值、方差、峰度和偏度;
(D2)设置待识别的电压暂降波形的时间窗口长度为D2个周波,提取待识别的电压暂降波形的每个时间窗的特征量,形成每个测试电压暂降波形时间序列特征量集合T={t1,t2,…,tn},称为测试波形时间序列特征量集合,其中,n为总的时间窗个数,tn为测试电压暂降波形第N个时间窗的特征量,N=1,2,…,n,特征量tn为以下6个量:最大值、最小值、平均值、方差、峰度和偏度;
(D3)基于动态弯曲距离相似度算法,分别计算每个典型电压暂降事件的样本波形时间序列特征量集合Ri与测试波形时间序列特征量集合T的相似度,得到一组匹配度系数。
前述的基于典型电压暂降波形匹配的电压暂降原因识别方法,其特征在于:(D3)基于动态弯曲距离相似度算法,分别计算每个典型电压暂降事件的样本波形时间序列特征量集合Ri与测试波形时间序列特征量集合T的相似度,得到一组匹配度系数,包括以下步骤,
(1)取两数据维数分别为n和mi的时间序列T’和Ri’,即T’=(t1,t2,…,tn)T,Ri’=(r1,r2,…,rmi)T,其中T’定义为测试时间序列,Ri’定义为样本时间序列;
(2)将时间序列T’和Ri’按照相应的时间序列的位置排序,构造时间序列距离矩阵An×mi,如公式(12)所示,
其中,An×mi中元素表示时间序列点tp和rq之间的欧氏距离;
(3)对于时间序列距离矩阵An×mi,将每一组相邻元素组成的集合成为弯曲路径W=(w1,w2,…,wS),弯曲路径W中的第s’个元素ws’=(apq)s’,且此路径满足如下条件,a)有界性:max{mi,n}≤S≤mi+n-1;b)临界条件:w1=a11,wS=anmi;c)连续性与单调性:对于ws’=(apq)s’,其邻接元素ws’-1=(ap’q’)s’‐1下标满足0≤p‐p′≤1,0≤q‐q′≤1;p、q分别为弯曲路径W的ws’在距离矩阵An×mi中对应元素的行下标和列下标;p’、q’分别为弯曲路径W的ws’‐1在距离矩阵An×mi中对应元素的行下标和列下标;
(4)计算时间序列T’和Ri’的最小弯曲路径θDTW,如公式(13)所示,
其中,p=1,2,3,…,n;q=1,2,3,…,mi;A中弯曲路径最小累加值为D(n,mi);D(p,q)表示距离矩阵元素apq与其前段的弯曲路径部分长度最小累计值之和;<>表示该位置的p或q为取范围内的任意值;
(5)根据(1)-(4)分别计算每个典型电压暂降事件的样本波形时间序列特征量集合Ri与测试波形时间序列特征量集合的最小弯曲距离Di(n,mi);
(6)将最小弯曲距离Di(n,mi)根据公式(14)进行变换,s(T,Ri)表示时间序列T’和Ri’之间的变换相似度,取值范围为[0,100%],其取值越大说明最小弯曲距离Di(n,mi)越小,时间序列T’和Ri’之间的相似度越大,即s(T,Ri)为得到的一组匹配度系数;
本发明的有益效果是:本发明的基于典型电压暂降波形匹配的电压暂降原因识别方法,建立典型电压暂降事件波形库,在波形的层面上对电压暂降进行匹配分析,更为基础也更为准确,能够快速准确的识别电压暂降产生的原因,对于电压暂降的预防和治理具有十分重要的参考价值,并具有以下特点,
(1)有利于提高电压暂降产生原因识别的准确性,从波形层面上对电压暂降进行匹配分析,避免了电压暂降典型特征量提取遗漏等问题;
(2)简化匹配分析的不确定程度,相比于其它人工智能算法,本方法无需使用大量的样本,无需为各个特征量进行阈值的选取等,降低了人为因素造成的判断失误;
(3)有利于提高电压暂降产生原因的识别速度,本方法中的计算主要产生在波形时间窗特征量以及动态弯曲距离的计算过程中,相比于人工智能算法,大大提升了计算速度,有利于电压暂降产生原因识别速度的提高。
附图说明
图1是本发明的基于典型电压暂降波形匹配的电压暂降原因识别方法的流程图。
图2为本发明进行动态弯曲距离计算的路径示意图。
具体实施方式
下面将结合说明书附图,对本发明做进一步说明。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
如图1所示,本发明的基于典型电压暂降波形匹配的电压暂降原因识别方法,包括以下步骤,
步骤(A),采集典型电压暂降历史事件的历史监测波形
通过暂降录波仪采集典型电压暂降历史事件的历史监测波形,并提取历史监测波形所对应的电压暂降产生原因,电压暂降产生原因,包括以下6种原因:三相短路故障、单相接地故障、两相接地故障、相间短路故障、感应电机启动、变压器投运。为保证所采集到的典型电压暂降历史事件的历史监测波形具有典型性,每种典型电压暂降历史事件至少采集10个波形,,选取过程,包括以下步骤,
(A1)计算某一电压暂降产生的原因所有电压暂降波形典型特征量X={x1,x2,x3,x4,x5,x6,x7},其中,x1~x7为列向量,每个列向量表示电压暂降产生的原因中所有电压暂降波形7个典型特征量之一,包括电压暂降平均幅值、最小幅值、持续时间、电压下降斜率、电压恢复斜率、电压相角跳变、三相幅值不平衡度,每个向量的维数表示该产生原因对应的电压暂降中所采集到的波形个数;
(A2)求取所有波形的每个特征值平均值得到向量
(A3)逐一计算所有波形的典型特征量向量与平均值向量的欧氏距离;
(A4)取所有波形中的欧氏距离最小值,并提取欧氏距离最小值对应的波形为最为典型的波形;
详细的具体过程为,
(1)以第P个波形为例,各个特征量的计算方法分别为,
电压暂降平均幅值:设电压暂降从发生到结束共有M个周波,UPQ表示第P个暂降波形中的第Q个周波幅值,在电压暂降发生时刻Q=1,在电压暂降结束时刻Q=M,电压暂降平均幅值UPave表示为:
电压暂降最小幅值:电压暂降M个周波中,最小的周波幅值UPmin表示为:
UPmin=min(UPQ) (2)
其中,Q=1,2,…,M;
电压暂降持续时间:设电压暂降开始时间为tPs,暂降结束时间为tPe,电压暂降持续时间tPdur可表示为:
tPdur=tPe-tPs (3)
电压下降斜率:电压从开始出现下降到电压降到最低值的平均斜率kPdown可以表示为:
kPdown=(1-UPmin)/|tPmin-tPs| (4)
其中,tPmin表示电压暂降幅值达到最小的时间;
电压上升斜率:电压从电压达到最小值恢复到正常电压水平的的平均斜率kPup可以表示为:
kPup=(1-UPmin)/|tPmin-tPe| (5)
电压相角跳变:电压暂降前电压过零点相位角为电压暂降过程中电压过零点相位角为则电压相角跳变为:
三相幅值不平衡度:系统发生电压暂降之后,三相电压暂降情况在很大一部分情况下不尽相同,三相幅值不平衡度UPε可以表示为:
其中,UPA、UPB、UPC分别为tmin时刻A、B、C三相基波电压幅值,max(UPA、UPB、UPC)为UPA、UPB、UPC中的最大值,min(UPA、UPB、UPC)为UPA、UPB、UPC中的最小值,UPN为电压标称值;
(2)求取所有波形某个特征值的平均值可以表示为:
其中,j=1,2,…,7,表示电压暂降波形的7个典型特征量,Q表示该类电压暂降共采集到Q个电压波形,J=1,2,…,Q;
分别计算每个特征量的平均值,得到向量
(3)逐一计算Q个波形的典型特征量向量xJ=(xJ1,xJ2,xJ3,xJ4,xJ5,xJ6,xJ7)与平均值向量的欧氏距离,
其中,J=1,2,…,Q,表示该类电压暂降共采集到Q个电压波形;
(4)取Q个波形中的欧氏距离最小值,并提取相应的波形作为最为典型的电压暂降事件波形;
步骤(B),建立形成典型电压暂降事件波形库
将步骤(A)中采集到的历史监测波形和对应的电压暂降产生原因,保存到典型电压暂降事件波形库中,建立电压暂降波形特征量和产生原因的对应关系表,形成典型电压暂降事件波形库;
步骤(C),提取待识别的电压暂降波形
将采集检测到的当前电压暂降事件波形作为测试波形,包括以下步骤,
(C1)在检测到原因未明的当前电压暂降事件并且计划对该电压暂降进行电压暂降原因识别时,记录该电压暂降事件发生时的三相波形;
(C2)提取三相波形其中发生电压暂降的部分波段,作为待识别的电压暂降波形;
步骤(D),对提取待识别的电压暂降波形的波段与典型电压暂降波形库进行波形匹配分析
基于动态弯曲距离相似度算法,对步骤(C)中提取待识别的电压暂降波形与典型电压暂降波形库中的典型电压暂降事件波形进行逐一匹配分析,得到一组匹配度系数,包括以下步骤,
(D1)设置典型电压暂降事件波形的时间窗长度为Di1个周波,提取各典型电压暂降事件波形每个时间窗的特征量,形成每个典型电压暂降事件波形的时间序列形特征量集合,称为样本波形时间序列特征量集合,其中,i=1,2,…,6,表示典型电压暂降事件波形库中共有6种典型波形,mi为第i种典型电压暂降事件波形总的时间窗个数,rik为第i个典型电压暂降事件波形第k个时间窗的特征量,k=1,2,…,mi,特征量rik为以下6个量:最大值、最小值、平均值、方差、峰度和偏度;
(D2)设置待识别的电压暂降波形的时间窗口长度为D2个周波,提取待识别的电压暂降波形的每个时间窗的特征量,形成每个测试电压暂降波形时间序列特征量集合T={t1,t2,…,tn},称为测试波形时间序列特征量集合,其中,n为总的时间窗个数,tn为测试电压暂降波形第N个时间窗的特征量,N=1,2,…,n,特征量tn为以下6个量:最大值、最小值、平均值、方差、峰度和偏度,
(D1)和(D2)中时间序列特征量集合的具体计算方法如下:
(1)确定滑动窗口的长度L;
(2)滑动窗口从该天第一个周波开始,每滑动一次,窗口就向后移动一个周波,例如对于第v个窗口wv={x(k):k=v-L+1,v-L+2,...,v},分别计算该窗口的最大值、最小值、平均值、方差、偏度和峰度。其中偏度和峰度的具体计算方法如下:
偏度有多种定义,正态分布以及所有对称分布的偏度都为0,并且广泛采用的计算公式为:
其中μ为变量X的平均值,μi为第i阶中心矩,E(·)为期望。
峰度是描述所有取值分布形态陡缓程度的统计量,正态分布的峰值为3,如果某分布比正态分布更为陡峭,则峰值大于3,峰度的计算公式为:
(3)重复(2)步,滑动到下一个窗口,计算该窗口的特征量集合,这样就得到所有窗口的特征量集合;
(D3)基于动态弯曲距离相似度算法,分别计算每个典型电压暂降事件的样本波形时间序列特征量集合Ri与测试波形时间序列特征量集合T的相似度,如图2所示,得到一组匹配度系数,包括以下步骤,
(1)取两数据维数分别为n和mi的时间序列T’和Ri’,即T’=(t1,t2,…,tn)T,Ri’=(r1,r2,…,rmi)T,其中T’定义为测试时间序列,Ri’定义为样本时间序列;
(2)将时间序列T’和Ri’按照相应的时间序列的位置排序,构造时间序列距离矩阵An×mi,如公式(12)所示,
其中,An×mi中元素表示时间序列点tp和rq之间的欧氏距离;
(3)对于时间序列距离矩阵An×mi,将每一组相邻元素组成的集合成为弯曲路径W=(w1,w2,…,wS),弯曲路径W中的第s’个元素ws’=(apq)s’,且此路径满足如下条件,a)有界性:max{mi,n}≤S≤mi+n-1;b)临界条件:w1=a11,wS=anmi;c)连续性与单调性:对于ws’=(apq)s’,其邻接元素ws’-1=(ap’q’)s’‐1下标满足0≤p‐p′≤1,0≤q‐q′≤1,p、q分别为弯曲路径W的ws’在距离矩阵An×mi中对应元素的行下标和列下标;p’、q’分别为弯曲路径W的ws’‐1在距离矩阵An×mi中对应元素的行下标和列下标。
(4)计算时间序列T’和Ri’的最小弯曲路径θDTW,如公式(13)所示,
其中,p=1,2,3,…,n;q=1,2,3,…,mi;A中弯曲路径最小累加值为D(n,mi);D(p,q)表示距离矩阵元素apq与其前段的弯曲路径部分长度最小累计值之和;<>表示该位置的p或q为取范围内的任意值;
(5)根据(1)-(4)分别计算每个典型电压暂降事件的样本波形时间序列特征量集合Ri与测试波形时间序列特征量集合的最小弯曲距离Di(n,mi);
(6)将最小弯曲距离Di(n,mi)根据公式(14)进行变换,s(T,Ri)表示时间序列T’和Ri’之间的变换相似度,取值范围为[0,100%],其取值越大说明最小弯曲距离Di(n,mi)越小,时间序列T’和Ri’之间的相似度越大,即s(T,Ri)为得到的一组匹配度系数;
步骤(E),识别当前电压暂降事件波形产生的原因
对步骤(D)中得到的匹配度系数按大小进行排序,选取匹配度最大的典型电压暂降事件波形作为最佳的匹配波形,提取典型电压暂降事件波形对应的电压暂降所产生的原因,为当前电压暂降事件波形产生原因的识别结果。
综上所述,本发明的基于典型电压暂降波形匹配的电压暂降原因识别方法,建立典型电压暂降事件波形库,在波形的层面上对电压暂降进行匹配分析,更为基础也更为准确,能够快速准确的识别电压暂降产生的原因,对于电压暂降的预防和治理具有十分重要的参考价值,并具有以下特点,
(1)有利于提高电压暂降产生原因识别的准确性,从波形层面上对电压暂降进行匹配分析,避免了电压暂降典型特征量提取遗漏等问题;
(2)简化匹配分析的不确定程度,相比于其它人工智能算法,本方法无需使用大量的样本,无需为各个特征量进行阈值的选取等,降低了人为因素造成的判断失误;
(3)有利于提高电压暂降产生原因的识别速度,本方法中的计算主要产生在波形时间窗特征量以及动态弯曲距离的计算过程中,相比于人工智能算法,大大提升了计算速度,有利于电压暂降产生原因识别速度的提高。
以上显示和描述了本发明的基本原理、主要特征及优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (6)
1.基于典型电压暂降波形匹配的电压暂降原因识别方法,其特征在于:包括以下步骤,
步骤(A),采集典型电压暂降历史事件的历史监测波形
通过暂降录波仪采集典型电压暂降历史事件的历史监测波形,并提取历史监测波形所对应的电压暂降产生原因;
步骤(B),建立形成典型电压暂降事件波形库
将步骤(A)中采集到的历史监测波形和对应的电压暂降产生原因,保存到典型电压暂降事件波形库中,建立电压暂降波形特征量和产生原因的对应关系表,形成典型电压暂降事件波形库;
步骤(C),提取待识别的电压暂降波形
将采集检测到的当前电压暂降事件波形作为测试波形;
步骤(D),对提取待识别的电压暂降波形的波段与典型电压暂降波形库进行波形匹配分析
基于动态弯曲距离相似度算法,对步骤(C)中提取待识别的电压暂降波形与典型电压暂降波形库中的典型电压暂降事件波形进行逐一匹配分析,得到一组匹配度系数;
步骤(E),识别当前电压暂降事件波形产生的原因
对步骤(D)中得到的匹配度系数按大小进行排序,选取匹配度最大的典型电压暂降事件波形作为最佳的匹配波形,提取典型电压暂降事件波形对应的电压暂降所产生的原因,为当前电压暂降事件波形产生原因的识别结果。
2.根据权利要求1所述的基于典型电压暂降波形匹配的电压暂降原因识别方法,其特征在于:步骤(A)提取历史监测波形所对应的电压暂降产生原因,包括以下6种原因:三相短路故障、单相接地故障、两相接地故障、相间短路故障、感应电机启动、变压器投运。
3.根据权利要求1所述的基于典型电压暂降波形匹配的电压暂降原因识别方法,其特征在于:步骤(A)为保证所采集到的典型电压暂降历史事件的历史监测波形具有典型性,每种典型电压暂降历史事件至少采集10个波形,并选取其中最为典型的波形作为历史监测波形,选取过程,包括以下步骤,
(A1)计算某一电压暂降产生的原因所有电压暂降波形典型特征量X={x1,x2,x3,x4,x5,x6,x7},其中,x1~x7为列向量,每个列向量表示电压暂降产生的原因中所有电压暂降波形7个典型特征量之一,包括电压暂降平均幅值、最小幅值、持续时间、电压下降斜率、电压恢复斜率、电压相角跳变、三相幅值不平衡度,每个向量的维数表示该产生原因对应的电压暂降中所采集到的波形个数;
(A2)求取所有波形的每个特征值平均值得到向量
(A3)逐一计算所有波形的典型特征量向量与平均值向量的欧氏距离;
(A4)取所有波形中的欧氏距离最小值,并提取欧氏距离最小值对应的波形为最为典型的波形。
4.根据权利要求1所述的基于典型电压暂降波形匹配的电压暂降原因识别方法,其特征在于:步骤(C),提取待识别的电压暂降波形,包括以下步骤,
(C1)在检测到原因未明的当前电压暂降事件并且计划对该电压暂降进行电压暂降原因识别时,记录该电压暂降事件发生时的三相波形;
(C2)提取三相波形其中发生电压暂降的部分波段,作为待识别的电压暂降波形。
5.根据权利要求1所述的基于典型电压暂降波形匹配的电压暂降原因识别方法,其特征在于:步骤(D),基于动态弯曲距离相似度算法,对步骤(C)中提取待识别的电压暂降波形与典型电压暂降波形库中的典型电压暂降事件波形进行逐一匹配分析,得到一组匹配度系数,包括以下步骤,
(D1)设置典型电压暂降事件波形的时间窗长度为Di1个周波,提取各典型电压暂降事件波形每个时间窗的特征量,形成每个典型电压暂降事件波形的时间序列形特征量集合称为样本波形时间序列特征量集合,其中,i=1,2,…,6,表示典型电压暂降事件波形库中共有6种典型波形,mi为第i种典型电压暂降事件波形总的时间窗个数,rik为第i个典型电压暂降事件波形第k个时间窗的特征量,k=1,2,…,mi,特征量rik为以下6个量:最大值、最小值、平均值、方差、峰度和偏度;
(D2)设置待识别的电压暂降波形的时间窗口长度为D2个周波,提取待识别的电压暂降波形的每个时间窗的特征量,形成每个测试电压暂降波形时间序列特征量集合T={t1,t2,…,tn},称为测试波形时间序列特征量集合,其中,n为总的时间窗个数,tn为测试电压暂降波形第N个时间窗的特征量,N=1,2,…,n,特征量tn为以下6个量:最大值、最小值、平均值、方差、峰度和偏度;
(D3)基于动态弯曲距离相似度算法,分别计算每个典型电压暂降事件的样本波形时间序列特征量集合Ri与测试波形时间序列特征量集合T的相似度,得到一组匹配度系数。
6.根据权利要求5所述的基于典型电压暂降波形匹配的电压暂降原因识别方法,其特征在于:(D3)基于动态弯曲距离相似度算法,分别计算每个典型电压暂降事件的样本波形时间序列特征量集合Ri与测试波形时间序列特征量集合T的相似度,得到一组匹配度系数,包括以下步骤,
(1)取两数据维数分别为n和mi的时间序列T’和R’i,即T’=(t1,t2,…,tn)T,R’i=(r1,r2,…,rmi)T,其中T’定义为测试时间序列,R’i定义为样本时间序列;
(2)将时间序列T’和R’i按照相应的时间序列的位置排序,构造时间序列距离矩阵An×mi,如公式(12)所示,
其中,An×mi中元素表示时间序列点tp和rq之间的欧氏距离;
(3)对于时间序列距离矩阵An×mi,将每一组相邻元素组成的集合成为弯曲路径W=(w1,w2,…,wS),弯曲路径W中的第s’个元素ws’=(apq)s’,且此路径满足如下条件,a)有界性:max{mi,n}≤S≤mi+n-1;b)临界条件:w1=a11,wS=anmi;c)连续性与单调性:对于ws’=(apq)s’,其邻接元素ws’-1=(ap’q’)s’‐1下标满足0≤p‐p′≤1,0≤q‐q′≤1;p、q分别为弯曲路径W的ws’在距离矩阵An×mi中对应元素的行下标和列下标;p’、q’分别为弯曲路径W的ws’‐1在距离矩阵An×mi中对应元素的行下标和列下标;
(4)计算时间序列T’和R’i的最小弯曲路径θDTW,如公式(13)所示,
其中,p=1,2,3,…,n;q=1,2,3,…,mi;A中弯曲路径最小累加值为D(n,mi);D(p,q)表示距离矩阵元素apq与其前段的弯曲路径部分长度最小累计值之和;<>表示该位置的p或q为取范围内的任意值;
(5)根据(1)-(4)分别计算每个典型电压暂降事件的样本波形时间序列特征量集合Ri与测试波形时间序列特征量集合的最小弯曲距离Di(n,mi);
(6)将最小弯曲距离Di(n,mi)根据公式(14)进行变换,s(T,Ri)表示时间序列T’和Ri’之间的变换相似度,取值范围为[0,100%],其取值越大说明最小弯曲距离Di(n,mi)越小,时间序列T’和Ri’之间的相似度越大,即s(T,Ri)为得到的一组匹配度系数;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610324282.XA CN105976266A (zh) | 2016-05-16 | 2016-05-16 | 基于典型电压暂降波形匹配的电压暂降原因识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610324282.XA CN105976266A (zh) | 2016-05-16 | 2016-05-16 | 基于典型电压暂降波形匹配的电压暂降原因识别方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105976266A true CN105976266A (zh) | 2016-09-28 |
Family
ID=56955835
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610324282.XA Pending CN105976266A (zh) | 2016-05-16 | 2016-05-16 | 基于典型电压暂降波形匹配的电压暂降原因识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105976266A (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107769197A (zh) * | 2017-11-14 | 2018-03-06 | 国网江苏省电力公司电力科学研究院 | 一种基于电网设备拓扑的电压暂降源定位系统 |
CN108051664A (zh) * | 2017-11-15 | 2018-05-18 | 全球能源互联网研究院有限公司 | 一种电能质量监测装置 |
CN108921409A (zh) * | 2018-06-21 | 2018-11-30 | 国网江苏省电力有限公司电力科学研究院 | 一种电压暂降严重度等级判定方法 |
CN108923398A (zh) * | 2018-07-23 | 2018-11-30 | 国网浙江省电力有限公司电力科学研究院 | 一种基于电压特征行波相似测度的直流配电网保护方法 |
CN109146319A (zh) * | 2018-09-12 | 2019-01-04 | 四川大学 | 一种基于欧氏距离法计算电压暂降类型的方法 |
CN109599874A (zh) * | 2018-10-29 | 2019-04-09 | 国网江苏省电力有限公司 | 基于电压暂降事件关联分析的治理设备选择方法 |
CN109782086A (zh) * | 2018-12-25 | 2019-05-21 | 武汉中原电子信息有限公司 | 一种基于多维度信号分析的非侵入负荷识别方法 |
CN110247420A (zh) * | 2019-07-17 | 2019-09-17 | 四川轻化工大学 | 一种hvdc输电线路故障智能识别方法 |
CN110726898A (zh) * | 2018-07-16 | 2020-01-24 | 北京映翰通网络技术股份有限公司 | 一种配电网故障类型识别方法 |
CN111341085A (zh) * | 2017-08-09 | 2020-06-26 | 绿铜能源科技有限公司 | 用于从电力监控器通过无线系统提供波形的系统和方法 |
CN112116013A (zh) * | 2020-09-24 | 2020-12-22 | 四川大学 | 一种基于波形特征的电压暂降事件归一化方法 |
CN112131956A (zh) * | 2020-08-27 | 2020-12-25 | 国网湖北省电力有限公司电力科学研究院 | 一种基于差异哈希算法的电压暂降源分类方法 |
CN112415258A (zh) * | 2020-11-05 | 2021-02-26 | 山东山大电力技术股份有限公司 | 电压暂态事件的特征幅值和持续时间的判别方法及系统 |
CN114118404A (zh) * | 2021-10-20 | 2022-03-01 | 广西电网有限责任公司 | 基于神经网络的电能表相位不平衡状态症结研判方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101572411A (zh) * | 2009-06-10 | 2009-11-04 | 天津市电力公司 | 基于Mamdani模糊推理的电压暂降源识别方法 |
CN103578050A (zh) * | 2013-11-14 | 2014-02-12 | 国家电网公司 | 一种电压暂降原因的识别方法 |
CN105277852A (zh) * | 2015-11-09 | 2016-01-27 | 广州航海学院 | 配电网线路工况分类识别方法 |
-
2016
- 2016-05-16 CN CN201610324282.XA patent/CN105976266A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101572411A (zh) * | 2009-06-10 | 2009-11-04 | 天津市电力公司 | 基于Mamdani模糊推理的电压暂降源识别方法 |
CN103578050A (zh) * | 2013-11-14 | 2014-02-12 | 国家电网公司 | 一种电压暂降原因的识别方法 |
CN105277852A (zh) * | 2015-11-09 | 2016-01-27 | 广州航海学院 | 配电网线路工况分类识别方法 |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11431601B2 (en) | 2017-08-09 | 2022-08-30 | Verdigris Technologies, Inc. | System and methods for providing waveforms over wireless systems from power monitors |
CN111341085B (zh) * | 2017-08-09 | 2022-05-17 | 绿铜能源科技有限公司 | 用于从电力监控器通过无线系统提供波形的系统和方法 |
CN111341085A (zh) * | 2017-08-09 | 2020-06-26 | 绿铜能源科技有限公司 | 用于从电力监控器通过无线系统提供波形的系统和方法 |
CN107769197A (zh) * | 2017-11-14 | 2018-03-06 | 国网江苏省电力公司电力科学研究院 | 一种基于电网设备拓扑的电压暂降源定位系统 |
CN108051664A (zh) * | 2017-11-15 | 2018-05-18 | 全球能源互联网研究院有限公司 | 一种电能质量监测装置 |
CN108921409B (zh) * | 2018-06-21 | 2021-10-15 | 国网江苏省电力有限公司电力科学研究院 | 一种电压暂降严重度等级判定方法 |
CN108921409A (zh) * | 2018-06-21 | 2018-11-30 | 国网江苏省电力有限公司电力科学研究院 | 一种电压暂降严重度等级判定方法 |
CN110726898A (zh) * | 2018-07-16 | 2020-01-24 | 北京映翰通网络技术股份有限公司 | 一种配电网故障类型识别方法 |
CN110726898B (zh) * | 2018-07-16 | 2022-02-22 | 北京映翰通网络技术股份有限公司 | 一种配电网故障类型识别方法 |
CN108923398A (zh) * | 2018-07-23 | 2018-11-30 | 国网浙江省电力有限公司电力科学研究院 | 一种基于电压特征行波相似测度的直流配电网保护方法 |
CN109146319A (zh) * | 2018-09-12 | 2019-01-04 | 四川大学 | 一种基于欧氏距离法计算电压暂降类型的方法 |
CN109599874A (zh) * | 2018-10-29 | 2019-04-09 | 国网江苏省电力有限公司 | 基于电压暂降事件关联分析的治理设备选择方法 |
CN109599874B (zh) * | 2018-10-29 | 2022-07-12 | 国网江苏省电力有限公司 | 基于电压暂降事件关联分析的治理设备选择方法 |
CN109782086B (zh) * | 2018-12-25 | 2020-12-18 | 武汉中原电子信息有限公司 | 一种基于多维度信号分析的非侵入负荷识别方法 |
CN109782086A (zh) * | 2018-12-25 | 2019-05-21 | 武汉中原电子信息有限公司 | 一种基于多维度信号分析的非侵入负荷识别方法 |
CN110247420A (zh) * | 2019-07-17 | 2019-09-17 | 四川轻化工大学 | 一种hvdc输电线路故障智能识别方法 |
CN112131956A (zh) * | 2020-08-27 | 2020-12-25 | 国网湖北省电力有限公司电力科学研究院 | 一种基于差异哈希算法的电压暂降源分类方法 |
CN112131956B (zh) * | 2020-08-27 | 2022-08-26 | 国网湖北省电力有限公司电力科学研究院 | 一种基于差异哈希算法的电压暂降源分类方法 |
CN112116013A (zh) * | 2020-09-24 | 2020-12-22 | 四川大学 | 一种基于波形特征的电压暂降事件归一化方法 |
CN112415258A (zh) * | 2020-11-05 | 2021-02-26 | 山东山大电力技术股份有限公司 | 电压暂态事件的特征幅值和持续时间的判别方法及系统 |
CN112415258B (zh) * | 2020-11-05 | 2021-08-24 | 山东山大电力技术股份有限公司 | 电压暂态事件的特征幅值和持续时间的判别方法及系统 |
CN114118404A (zh) * | 2021-10-20 | 2022-03-01 | 广西电网有限责任公司 | 基于神经网络的电能表相位不平衡状态症结研判方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105976266A (zh) | 基于典型电压暂降波形匹配的电压暂降原因识别方法 | |
CN104865499B (zh) | 一种特高压直流输电线路区内外故障识别方法 | |
CN104184144B (zh) | 一种用于多电压等级电网模型的抗差状态估计方法 | |
CN101227092B (zh) | 一种配电网络中的谐波源识别方法 | |
CN104050366B (zh) | 一种动态无功补偿装置响应时间测试方法 | |
CN106250660A (zh) | 一种海上风电场谐振分析方法 | |
CN103578050A (zh) | 一种电压暂降原因的识别方法 | |
CN109840691B (zh) | 基于深度神经网络的非侵入式分项电量估计方法 | |
CN103824129A (zh) | 一种基于动态阈值的高铁电能质量异常状态预警方法 | |
CN106443244A (zh) | 用电设备类型辨识方法和系统 | |
CN101572411A (zh) | 基于Mamdani模糊推理的电压暂降源识别方法 | |
CN110244099A (zh) | 基于用户电压的窃电检测方法 | |
CN103116097A (zh) | 基于多断面混合量测信息的设备参数在线辨识方法 | |
CN104714075A (zh) | 一种电网电压闪变包络参数提取方法 | |
CN105141253A (zh) | 基于光伏出力曲线斜率的光伏窃电识别方法 | |
CN106443253A (zh) | 一种基于pmu数据的输电线路参数辨识方法 | |
CN108614154B (zh) | 一种基于频谱分析的电力用户相位识别方法 | |
CN109444505A (zh) | 一种基于变分贝叶斯参数学习方法的电动汽车充电站谐波电流检测算法 | |
CN105048473B (zh) | 变电站自动电压控制中低压无功设备电压灵敏度统计方法 | |
CN107167689B (zh) | 考虑联络变压器变比的概率短路灵敏度计算方法 | |
Kazaki et al. | Cluster analysis of university campus smart meter data | |
CN113702767B (zh) | 基于小波滑窗能量的孤岛直流微电网故障诊断方法 | |
CN107482778A (zh) | 一种改进的电力系统健康状态监测的方法及系统 | |
CN106291239A (zh) | 一种采用滤波器支路电流和主成分分析方法的直流输电线路故障识别方法 | |
CN103913654B (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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160928 |