CN106771871A - 一种基于vmd与柔性形态学的输电线路故障行波测距方法 - Google Patents
一种基于vmd与柔性形态学的输电线路故障行波测距方法 Download PDFInfo
- Publication number
- CN106771871A CN106771871A CN201611230923.1A CN201611230923A CN106771871A CN 106771871 A CN106771871 A CN 106771871A CN 201611230923 A CN201611230923 A CN 201611230923A CN 106771871 A CN106771871 A CN 106771871A
- Authority
- CN
- China
- Prior art keywords
- flexible
- vmd
- edge detection
- mode
- travelling wave
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/081—Locating faults in cables, transmission lines, or networks according to type of conductors
- G01R31/085—Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution lines, e.g. overhead
Abstract
本发明公开了一种基于VMD与柔性形态学的输电线路故障行波测距方法。本方法为:首先采集故障发生后的两端的电压数据,对原始故障行波信号进行截断,取故障后3ms的数据;然后对故障行波进行相模变换;再然后对其进行VMD(变分模态分解)分解,分解过程中自动滤去白噪声,其分解得到的IMF1分量集中包含了故障信息,再采用柔性形态边缘检测方法对IMF1分量进行处理,结合两者的优势可以很清晰的得到行波波头达到时刻。与以前的测距方法相比行波波头的检测效果更好,精确度更高,误差在工程允许的范围内。
Description
技术领域
本发明属于输电线路故障定位技术领域,特别涉及一种VMD与柔性形态学相结合的输电线路故障定位方法。
背景技术
对于电力系统来说,输电线路故障测距是必不可少的一部分。当发生线路故障后,传统人工巡线方式比较浪费时间和精力,高效、快速的故障定位可以及时发现故障位置,从而可排除故障,及时恢复供电以减少社会经济损失。
故障测距有故障分析法和行波法,故障分析法的精度不高,然而行波的定位精度高,所以现在行波检测成为了热点。目前行波的检测方法主要有小波变换(WT)和希尔伯特-黄变换(HHT)两种。小波模极大值法虽然对故障信号奇异性检测有良好的效果,但是小波变换需要选取适合特定信号的小波基和分解尺度,不能达到需要的满意效果。HHT方法使用EMD进行分解,EMD的分解效果受包络方式影响,存在端点效应和模态混叠问题,难以准确的确定行波波头的到达时刻。变分模态分解(VMD)是这几年新出来的一种自适应信号分解方法,分解模态稳定性好,能够更好地反应信号的奇异性特征。柔性形态学的膨胀和腐蚀操作可以看作是对原信号的扩张和收缩,将膨胀与腐蚀的结果做差值可以有效放大信号的突变点,精确度高,实现故障点的准确定位。
发明内容
本发明的目的在于克服现有技术的缺点与不足,提供了一种基于VMD与柔性形态学的输电线路故障行波测距方法,
本发明所采用的技术方案是:
一种基于VMD与柔性形态学的输电线路故障行波测距方法,具体步骤如下:
步骤1:对原始故障行波信号进行截断,取故障后3ms的数据;
步骤2:对截断的数据进行相模变换;
步骤3:对线模分量进行VMD分解,同时去除了白噪声;
步骤4:寻找固有模态函数IMF1分量进行柔性形态边缘检测;
步骤5:寻找柔性形态边缘检测结果图中第一个幅值最大脉冲确定故障行波波头初始到达两端的时刻tM、tN,根据测距公式计算故障距离,式中:L为线路的水平长度,dMf为故障点到一端的距离,tM、tN为故障点行波到线路两端的时间。
进一步,步骤3中对线模分量进行VMD分解的具体步骤为:
(1)初始化模态各模态对应的中心频率拉格朗日乘子和迭代次数n;模态个数K取3;
(2)计算更新模态i表示第i个,1≤i≤k,1≤k≤K,α表示数据保真约束的平衡参数;
(3)计算中心频率
(4)计算更新拉格朗日乘子τ表示双上升时间步长;
(5)如果则迭代停止,否则返回执行步骤(2)到(4)。
进一步,步骤4中对IMF1分量进行柔性形态边缘检测,具体如下
(1)对IMF1采用适应性较强的膨胀腐蚀型算子进行柔性形态边缘检测:
其中表示柔性膨胀:
其中fΘ[B,A,r]表示柔性腐蚀:
其中[B,A,r]为结构元素;r◇f(a)为重复集;r为重复度;
(2)设定基于噪声强度的软阈值方法
噪声强度定义为:
g(n)表示柔性形态边缘检测结果;N表示采样点个数。
边缘检测算子的输出阈值:
形态边缘检测定位结果:
f(n)表示被检测信号。
本发明的有益效果:
VMD为完全非递归分解模型且运算效率高,具有良好的抗噪性。形态学膨胀、腐蚀运算具有扩张性和收缩性,能够有效放大信号的奇异点。
附图说明
图1是本发明的220KV输电线路仿真模型图;
图2是本发明的两端测得的线模分量图;
其中:图(a)为一端检测到的电压线模分量图;图(b)为另一端检测到的电压线模分量图;
图3是本发明的两端VMD分解后的IMF1分量图;
其中:图(a)为一端VMD分解后的IMF1分量图;图(b)为另一端VMD分解后的IMF1分量图;
图4是本发明的柔性形态边缘检测结果图;
其中图(a)为一端的柔性形态边缘检测结果图;图(b)为另一端的柔性形态边缘检测结果图。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步详细的说明。
本发明具体按照以下步骤实施:
步骤1:采集故障发生后的两端的电压数据,对原始故障行波信号进行截断,取故障后3ms的数据;
步骤2:对截断的数据进行相模变换;
步骤3:对线模分量进行VMD分解,同时去除了白噪声;
VMD为完全非递归分解模型且运算效率高,具有良好的抗噪性,克服了EMD存在的模态混叠等不足。基于干扰信号性质与数量,预设分解模态数,通过提取信号的模态进而将信号从含噪信号中剥离出来。
线模分量进行VMD分解具体步骤为:
(1)初始化模态各模态对应的中心频率拉格朗日乘子和迭代次数n;模态个数K取3;
(2)计算更新模态i表示第i个,1≤i≤k,1≤k≤K,α表示数据保真约束的平衡参数;
(3)计算中心频率
(4)计算更新拉格朗日乘子τ表示双上升时间步长;
(5)如果则迭代停止,否则返回执行步骤(2)到(4)。
步骤4:寻找IMF1分量并进行柔性形态边缘检测;
对IMF1分量进行柔性形态边缘检测,具体如下:
(1)对IMF1采用适应性较强的膨胀腐蚀型算子进行柔性形态边缘检测:
其中表示柔性膨胀:
其中fΘ[B,A,r]表示柔性腐蚀:
其中[B,A,r]为结构元素;r◇f(a)为重复集;r为重复度;
(2)设定基于噪声强度的软阈值方法
噪声强度定义为:
边缘检测算子的输出阈值:
形态边缘检测定位结果:
步骤5:寻找柔性形态边缘检测结果图中第一个幅值最大脉冲确定故障行波波头初始到达两端的时刻tM、tN,根据测距公式计算故障距离,式中:L为线路的水平长度,dMf为故障点到一端的距离,tM、tN为故障点行波到线路两端的时间。
实施例
选取单极接地故障的一个典型算例:220KV输电线路仿真模型,如图1所示,设置距离线路一端30kM处发生单极接地故障,过渡电阻为50Ω,从两端获得行波数据,截取故障发生后的3ms数据,将故障行波数据进行相模变换得到线模分量,如图2所示,对线模分量进行VMD分解获得IMF1分量,如图3所示,对IMF1分量进行柔性形态边缘检测,检测结果,如图4所示。通过图4可知两端到达点对应为102和571,将两个数值代入测距公式得到故障距离线路一端的距离。
本发明对不同故障距离、不同故障类型进行了仿真验证,结果如下表1所示。
表1不同故障下的测距结果
上文所列出的一系列的详细说明仅仅是针对本发明的可行性实施方式的具体说明,它们并非用以限制本发明的保护范围,凡未脱离本发明技艺精神所作的等效实施方式或变更均应包含在本发明的保护范围之内。
Claims (3)
1.一种基于VMD与柔性形态学的输电线路故障行波测距方法,其特征在于,包括如下步骤:
步骤1:采集故障发生后的两端的电压数据,对原始故障行波信号进行截断,取故障后3ms的数据;
步骤2:对截断的数据进行相模变换;
步骤3:对线模分量进行VMD分解,同时去除了白噪声;基于干扰信号性质与数量,预设分解模态数,通过提取信号的模态进而将信号从含噪信号中剥离出来;
步骤4:寻找IMF1分量并进行柔性形态边缘检测;
步骤5:寻找柔性形态边缘检测结果图中第一个幅值最大脉冲确定故障行波波头初始到达两端的时刻tM、tN,根据测距公式计算故障距离,式中:L为线路的水平长度,dMf为故障点到一端的距离,tM、tN为故障点行波到线路两端的时间。
2.根据权利要求1所述的一种基于VMD与柔性形态学的输电线路故障行波测距方法,其特征在于,所述步骤3中对线模分量进行VMD分解的具体步骤为:
步骤3.1:初始化模态各模态对应的中心频率拉格朗日乘子和迭代次数n;模态个数K取3;
步骤3.2:计算更新模态i表示第i个,1≤i≤k,1≤k≤K,α表示数据保真约束的平衡参数;
步骤3.3:计算中心频率
步骤3.4:计算更新拉格朗日乘子τ表示双上升时间步长;
步骤3.5:如果则迭代停止,否则返回执行步骤步骤3.2到步骤3.4。
3.根据权利要求1所述的一种基于VMD与柔性形态学的输电线路故障行波测距方法,其特征在于,所述步骤4中对IMF1分量进行柔性形态边缘检测的具体步骤如下:
步骤4.1:对IMF1采用适应性较强的膨胀腐蚀型算子进行柔性形态边缘检测:
其中表示柔性膨胀:
其中fΘ[B,A,r]表示柔性腐蚀:
f(n)Θ[B,A,r]=min{[r◇f(a),a∈An]∪[f(b):b∈(B/A)n]}
其中[B,A,r]为结构元素;r◇f(a)为重复集;r为重复度;
步骤4.2:设定基于噪声强度的软阈值方法
噪声强度定义为:
边缘检测算子的输出阈值:
形态边缘检测定位结果:
。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611230923.1A CN106771871A (zh) | 2016-12-28 | 2016-12-28 | 一种基于vmd与柔性形态学的输电线路故障行波测距方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611230923.1A CN106771871A (zh) | 2016-12-28 | 2016-12-28 | 一种基于vmd与柔性形态学的输电线路故障行波测距方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106771871A true CN106771871A (zh) | 2017-05-31 |
Family
ID=58922477
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611230923.1A Pending CN106771871A (zh) | 2016-12-28 | 2016-12-28 | 一种基于vmd与柔性形态学的输电线路故障行波测距方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106771871A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109001594A (zh) * | 2018-07-26 | 2018-12-14 | 国网湖南省电力有限公司 | 一种故障行波定位方法 |
CN109239528A (zh) * | 2018-09-18 | 2019-01-18 | 四川大学 | 基于形态学的架空输电线路高阻抗故障分析方法及系统 |
CN109633271A (zh) * | 2019-01-17 | 2019-04-16 | 长沙理工大学 | 基于变分模态分解和维格纳威尔分布的行波时频分析方法 |
CN111308272A (zh) * | 2020-03-09 | 2020-06-19 | 西南交通大学 | 一种小电流接地故障区段定位方法 |
CN112305378A (zh) * | 2020-10-21 | 2021-02-02 | 中国科学院声学研究所南海研究站 | 一种光电复合海缆故障测距系统及故障测距方法 |
CN113075503A (zh) * | 2021-03-30 | 2021-07-06 | 新疆大学 | 一种直流输电线路双端行波测距方法及系统 |
CN113376478A (zh) * | 2021-06-22 | 2021-09-10 | 清华大学 | 一种基于边缘检测的输电线路雷击或短路故障定位方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000200874A (ja) * | 1999-01-06 | 2000-07-18 | Hitachi Ltd | 半導体集積回路装置 |
CN105160077A (zh) * | 2015-08-07 | 2015-12-16 | 中国电力科学研究院 | 一种高压直流线路电晕电流变化趋势提取方法 |
CN105717422A (zh) * | 2015-12-04 | 2016-06-29 | 国家电网公司 | 一种高压电力设备局部放电特征提取方法和装置 |
-
2016
- 2016-12-28 CN CN201611230923.1A patent/CN106771871A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000200874A (ja) * | 1999-01-06 | 2000-07-18 | Hitachi Ltd | 半導体集積回路装置 |
CN105160077A (zh) * | 2015-08-07 | 2015-12-16 | 中国电力科学研究院 | 一种高压直流线路电晕电流变化趋势提取方法 |
CN105717422A (zh) * | 2015-12-04 | 2016-06-29 | 国家电网公司 | 一种高压电力设备局部放电特征提取方法和装置 |
Non-Patent Citations (3)
Title |
---|
包广清 等: "基于经验模态分解和形态学的风电并网电压故障检测", 《农业工程学报》 * |
宋平岗 等: "基于柔性形态滤波与信息熵的电能质量暂态扰动定位分析", 《科学技术与工程》 * |
高艳丰 等: "基于VMD和TEO的高压输电线路雷击故障测距研究", 《基于VMD和TEO的高压输电线路雷击故障测距研究》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109001594A (zh) * | 2018-07-26 | 2018-12-14 | 国网湖南省电力有限公司 | 一种故障行波定位方法 |
CN109001594B (zh) * | 2018-07-26 | 2021-03-16 | 国网湖南省电力有限公司 | 一种故障行波定位方法 |
CN109239528A (zh) * | 2018-09-18 | 2019-01-18 | 四川大学 | 基于形态学的架空输电线路高阻抗故障分析方法及系统 |
CN109633271A (zh) * | 2019-01-17 | 2019-04-16 | 长沙理工大学 | 基于变分模态分解和维格纳威尔分布的行波时频分析方法 |
CN111308272A (zh) * | 2020-03-09 | 2020-06-19 | 西南交通大学 | 一种小电流接地故障区段定位方法 |
CN112305378A (zh) * | 2020-10-21 | 2021-02-02 | 中国科学院声学研究所南海研究站 | 一种光电复合海缆故障测距系统及故障测距方法 |
CN113075503A (zh) * | 2021-03-30 | 2021-07-06 | 新疆大学 | 一种直流输电线路双端行波测距方法及系统 |
CN113376478A (zh) * | 2021-06-22 | 2021-09-10 | 清华大学 | 一种基于边缘检测的输电线路雷击或短路故障定位方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106771871A (zh) | 一种基于vmd与柔性形态学的输电线路故障行波测距方法 | |
CN101509949B (zh) | 直流输电线路双端非同步且参数自适应的故障测距时域法 | |
CN108627741B (zh) | 一种基于故障指示器的行波—阻抗法双端带支路配电网故障定位方法 | |
CN110969194B (zh) | 基于改进型卷积神经网络的电缆早期故障定位方法 | |
CN109524972B (zh) | 基于gso和svm算法的低频振荡参数估计方法 | |
CN105445624A (zh) | 结合小波变换及曲线拟合的电缆故障定位方法 | |
CN109557429A (zh) | 基于改进小波阈值去噪的gis局部放电故障检测方法 | |
CN109387744B (zh) | 基于奇异值分解的配网线路故障点定位方法及装置 | |
CN106093705B (zh) | 一种配电网单相接地故障行波波头的计算方法 | |
CN108957403B (zh) | 一种基于广义互相关的高斯拟合包络时延估计方法及系统 | |
CN104503432A (zh) | 一种基于小波能量的自主式水下机器人故障辨识方法 | |
CN115980514B (zh) | 一种基于多端行波频率矩阵的复杂配电网故障定位方法 | |
CN110837791A (zh) | 一种基于过完备字典的声速剖面反演方法 | |
CN105223467B (zh) | 基于分维计算和mallat分解的配电网故障选线方法 | |
CN104456089A (zh) | 一种流体管道多泄漏点定位方法 | |
CN105929380A (zh) | 一种卫星激光高度计全波形激光雷达数据去噪方法 | |
CN103605856B (zh) | 基于分数阶线路模型的输电线路参数估计方法 | |
CN107179473A (zh) | 一种输电线路故障定位方法 | |
CN112540260B (zh) | 基于行波能量变化特性的高压输电网混联线路故障测距方法、装置及系统 | |
CN109324267B (zh) | 基于双采样率的配网线路故障点定位方法及装置 | |
CN111289991A (zh) | 一种基于多场景的激光测距方法及装置 | |
CN113625101B (zh) | 基于果蝇算法与随机共振的行波信号处理方法 | |
CN104808074A (zh) | 高空电磁脉冲在海水中传播特性的计算方法 | |
CN113098006B (zh) | 新能源发电设备的阻抗特性曲线及传递函数拟合方法 | |
CN108427063A (zh) | 一种基于交叉验证法、经验模态和Teager峭度的首波定位方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for 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: 20170531 |