CN115078909B - 一种基于能量拐点的非理想脉冲前沿检测方法及装置 - Google Patents

一种基于能量拐点的非理想脉冲前沿检测方法及装置 Download PDF

Info

Publication number
CN115078909B
CN115078909B CN202210694580.3A CN202210694580A CN115078909B CN 115078909 B CN115078909 B CN 115078909B CN 202210694580 A CN202210694580 A CN 202210694580A CN 115078909 B CN115078909 B CN 115078909B
Authority
CN
China
Prior art keywords
inflection point
window
ideal pulse
leading edge
energy
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.)
Active
Application number
CN202210694580.3A
Other languages
English (en)
Other versions
CN115078909A (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.)
State Grid Jiangsu Electric Power Co ltd Suqian Power Supply Branch
Suqian Electric Power Design Institute Co ltd
Original Assignee
Suqian Electric Power Design Institute Co ltd
State Grid Jiangsu Electric Power Co ltd Suqian Power Supply Branch
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 Suqian Electric Power Design Institute Co ltd, State Grid Jiangsu Electric Power Co ltd Suqian Power Supply Branch filed Critical Suqian Electric Power Design Institute Co ltd
Priority to CN202210694580.3A priority Critical patent/CN115078909B/zh
Publication of CN115078909A publication Critical patent/CN115078909A/zh
Application granted granted Critical
Publication of CN115078909B publication Critical patent/CN115078909B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/08Locating faults in cables, transmission lines, or networks
    • G01R31/081Locating faults in cables, transmission lines, or networks according to type of conductors
    • G01R31/083Locating faults in cables, transmission lines, or networks according to type of conductors in cables, e.g. underground
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/08Locating faults in cables, transmission lines, or networks
    • G01R31/088Aspects of digital computing
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/08Locating faults in cables, transmission lines, or networks
    • G01R31/11Locating faults in cables, transmission lines, or networks using pulse reflection methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/15Correlation function computation including computation of convolution operations

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
  • Locating Faults (AREA)

Abstract

本发明公开一种基于能量拐点的非理想脉冲前沿检测方法及装置。对于一个具有非理想脉冲前沿的波形,该方法将波形累积能量的拐点视为其脉冲前沿的开始时刻,通过对累积能量拐点的检测达到非理想脉冲前沿检测的目的。本发明实现了基于能量拐点的非理想脉冲前沿检测方法并设计了配套装置,该方法具有抗噪声能力强、准确度高的优点,有效地减小了非理想脉冲前沿检测的误差。

Description

一种基于能量拐点的非理想脉冲前沿检测方法及装置
技术领域
本发明涉及边缘检测算法领域,具体的是一种基于能量拐点的非理想脉冲前沿检测方法及装置。
背景技术
电缆广泛应用于电网配电系统和铁路、航空、制造等行业的供电系统中。为了避免雷击、风害等外力的影响,电缆通常埋设在系统内部或地下,所以电缆的故障无法像架空线路一样通过巡检进行识别和定位,而是需要特殊的方法。
行波法是一种得到广泛应用的电缆故障定位方法,它利用电缆不连续点处的反射行波来对故障进行检测和定位,需要对反射波的到达时间,即反射波的前沿进行精确的测量。然而,由于电缆本身的色散作用、接受设备的带宽限制和传播路径上的信号干扰以及发送等影响,反射波在传播过程中会发生畸变,脉冲前沿不再满足理想阶跃信号的特征,造成基于行波到达时间的电缆故障定位算法存在偏差。因此,亟需一种具有强抗噪声能力、高准确度的非理想脉冲前沿检测新方法,减小非理想脉冲前沿检测的误差。
发明内容
为解决上述背景技术中提到的不足,本发明的目的在于提供一种基于能量拐点的非理想脉冲前沿检测方法及装置,该方法具有抗噪声能力强、准确度高的优点,有效地减小了非理想脉冲前沿检测的误差。
本发明是采取以下技术方案实现的:
一种基于能量拐点的非理想脉冲前沿检测方法,是对于一个具有非理想脉冲前沿的波形,将波形累积能量的拐点视为其脉冲前沿的开始时刻,通过对累积能量拐点的检测达到非理想脉冲前沿检测的目的。
具体包括以下步骤:
S1、采样得到原始波形;
S2、检测波形累积能量的拐点;
S3、依据累积能量拐点得到非理想脉冲前沿开始时刻。
所述步骤S2中,累积能量拐点的检测,具体为采用滑动窗口算法对累积能量进行拐点检测,具体的是:
S21、设置两个相邻的窗口沿着波形的时间方向滑动;
S22、计算第一窗口与第二窗口的差异;
S23、依据差异检测拐点。
所述第一窗口与第二窗口的差异,具体定义为:
d(fa..t,ft..b)=|Ea..t-Et..b| (1)
式中,d表示第一窗口与第二窗口的差异,下标“a..t”表示“从a时刻到t时刻”且a和t分别表示一个时间窗口的左右边沿,f表示信号本身,E表示累积能量。
优选地,累积能量E的计算方法,具体为:
在实际应用中,可根据效果与计算量更换被积分项,包括且不限于f(τ)2、|f(τ)|等。
在步骤S23中,所述拐点的定义具体为:将第一窗口与第二窗口的差异最大的位置视为拐点。
在步骤S23中,依据差异检测拐点的步骤,具体为:
给定滑动时间窗口的宽度为w,定义滑动窗口的拐点检测函数z(t)为:
z(t)=d(ft-w..t,ft..t+w) (3)
则拐点所在时刻,即非理想脉冲前沿的开始时刻为:
式中,tdect为检测到的非理想脉冲前沿开始时刻。
本发明的有益效果:
本发明基于能量拐点,实现了非理想脉冲前沿的准确检测,提升了非理想脉冲前沿检测的抗噪声能力,有效地减小了非理想脉冲前沿检测的误差。
附图说明
下面结合附图对本发明作进一步的说明。
图1是本发明基于能量拐点的非理想脉冲前沿检测方法及装置的总体流程图;
图2是本发明基于能量拐点的非理想脉冲前沿检测方法及装置的硬件装置框图;
图3是本发明基于能量拐点的非理想脉冲前沿检测方法及装置的具体实施效果图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
如图1所示,一种基于能量拐点的非理想脉冲前沿检测方法,包括以下步骤:
S1、采样得到原始波形;
S2、设置第一、第二滑动时间窗口,计算差异;
S3、求解差异最大位置,得到累积能量拐点;
S4、得到非理想脉冲前沿开始时刻。
如图2所示,一种基于能量拐点的非理想脉冲前沿检测装置,包括分压器、信号采集设备、信号传输设备和信号处理设备。
在本实施例中,为实现高压信号的衰减,满足信号采集设备的耐压范围,使用高压宽带阻容分压器对原始信号分压;使用高速数据采集卡采集信号;所述信号传输设备使用无线传输模块传输信号,实现数据的远程监控和上传;并在STM32板卡上进行信号处理,对非理想脉冲前沿进行检测,从而实现针对ns级别高压脉冲的非理想前沿检测。
所述高速数据采集卡可采用市售的采样频率为百兆Hz级别的高速数据采集卡,并使用LabView或MATLAB驱动板卡,以在ns级别有效检测非理想脉冲前沿。
信号处理设备采用MCU(Micro controller unit),以降低计算资源消耗。
一种基于能量拐点的非理想脉冲前沿检测方法,具体为:对于一个具有非理想脉冲前沿的波形,将波形累积能量的拐点视为其脉冲前沿的开始时刻,通过对累积能量拐点的检测达到非理想脉冲前沿检测的目的。
在本实施例中,累积能量拐点的检测,具体为采用滑动窗口算法对累积能量进行拐点检测:
设置两个相邻的窗口沿着波形的时间方向滑动,计算第一窗口与第二窗口间的差异,并依据差异检测拐点。
第一窗口与第二窗口的差异,具体定义为:
d(fa..t,ft..b)=|Ea..t-Et..b| (1)
式中,d表示第一窗口与第二窗口的差异,下标“a..t”表示“从a时刻到t时刻”且a和t分别表示一个时间窗口的左右边沿,f表示信号本身,E表示累积能量。
在本实施例中,累积能量E的计算方法,具体为:
在实际应用中,可根据效果与计算量更换被积分项,包括且不限于f(τ)2、|f(τ)|等。
拐点的定义具体为:将第一窗口与第二窗口间差异最大的位置视为拐点。
依据差异检测拐点,具体为:
给定滑动时间窗口的宽度为w,在本实施例中,滑动时间窗口的宽度w选择为100Δt,其中,Δt为采样间隔,定义滑动窗口的拐点检测函数z(t)为:
z(t)=d(ft-w..t,ft..t+w) (3)
则拐点所在时刻,即非理想脉冲前沿的开始时刻tdect为:
使用上述方法对一个具有非理想脉冲前沿的波形进行检测,该非理想脉冲前沿上升速率为25ns,波形信噪比为20dB,并与归一化阈值比较方法、赤池信息准则方法进行比较,其中,归一化阈值选择为25%,得到图3所示的实施效果图。
以上实例仅为说明本发明方法的具体流程,不能以此限定本发明的权利保护范围,凡是按照本发明提出的技术思想,在本发明权利要求基础上所做的任何改动,均落入本发明的保护范围之内。

Claims (4)

1.一种基于能量拐点的非理想脉冲前沿检测方法,其特征在于:对于一个具有非理想脉冲前沿的波形,将波形累积能量的拐点视为其脉冲前沿的开始时刻,通过对累积能量拐点的检测达到非理想脉冲前沿检测的目的;
所述检测方法具体包括以下步骤:
S1、采样得到原始波形;
S2、检测波形累积能量的拐点;
S3、依据累积能量拐点得到非理想脉冲前沿开始时刻;
所述步骤S2中,累积能量拐点的检测,具体为采用滑动窗口算法对累积能量进行拐点检测,具体的是:
S21、设置两个相邻的窗口沿着波形的时间方向滑动;
S22、计算第一窗口与第二窗口的差异;
S23、依据差异检测拐点;
所述第一窗口与第二窗口的差异,具体定义为:
式中,d表示第一窗口与第二窗口的差异,下标“a..t”表示“从a时刻到t时刻”且at分别表示一个时间窗口的左右边沿,f表示信号本身,E表示累积能量;
在步骤S23中,所述拐点的定义具体为:将第一窗口与第二窗口的差异最大的位置视为拐点;
在步骤S23中,依据差异检测拐点的步骤,具体为:
给定滑动时间窗口的宽度为w,定义滑动窗口的拐点检测函数z(t)为:
则拐点所在时刻,即非理想脉冲前沿的开始时刻为:
式中,t dect为检测到的非理想脉冲前沿开始时刻。
2.根据权利要求1所述的基于能量拐点的非理想脉冲前沿检测方法,其特征在于,所述累积能量E的计算方法,具体为:
3.一种用于权利要求1所述的基于能量拐点的非理想脉冲前沿检测方法的检测装置,其特征在于,包括分压器、信号采集设备、信号传输设备和信号处理设备;
所述分压器用于对原始信号分压,实现高压信号的衰减,满足信号采集设备的耐压范围;
所述信号采集设备用于采集信号,采用高速数据采集卡;
所述信号传输设备用于实现数据的远程监控和上传;
所述信号处理设备用于对非理想脉冲前沿进行检测,从而实现针对ns级别高压脉冲的非理想前沿检测。
4.根据权利要求3所述的基于能量拐点的非理想脉冲前沿检测装置,其特征在于,所述高速数据采集卡采用采样频率为百兆Hz级别的高速数据采集卡。
CN202210694580.3A 2022-06-16 2022-06-16 一种基于能量拐点的非理想脉冲前沿检测方法及装置 Active CN115078909B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210694580.3A CN115078909B (zh) 2022-06-16 2022-06-16 一种基于能量拐点的非理想脉冲前沿检测方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210694580.3A CN115078909B (zh) 2022-06-16 2022-06-16 一种基于能量拐点的非理想脉冲前沿检测方法及装置

Publications (2)

Publication Number Publication Date
CN115078909A CN115078909A (zh) 2022-09-20
CN115078909B true CN115078909B (zh) 2024-03-15

Family

ID=83252852

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210694580.3A Active CN115078909B (zh) 2022-06-16 2022-06-16 一种基于能量拐点的非理想脉冲前沿检测方法及装置

Country Status (1)

Country Link
CN (1) CN115078909B (zh)

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6509740B1 (en) * 1999-04-09 2003-01-21 Jovial Test Equipment, Inc. Methods and apparatus for time domain reflectometry
JP2009216587A (ja) * 2008-03-11 2009-09-24 Central Res Inst Of Electric Power Ind 送電線故障箇所標定方法、装置並びにプログラム
CN103323032A (zh) * 2013-05-17 2013-09-25 厦门科华恒盛股份有限公司 带死区单元双滑动窗滤波的阶跃信号检测方法
CN105223481A (zh) * 2015-10-22 2016-01-06 西安交通大学 基于差值能量函数的局部放电特高频信号起始时刻确定方法
CN105319487A (zh) * 2015-10-30 2016-02-10 国家电网公司 变电站局部放电信号检测和定位系统及方法
CN105574284A (zh) * 2015-12-29 2016-05-11 山东鲁能软件技术有限公司 一种基于趋势特征点的电力设备故障诊断方法
CN107102244A (zh) * 2017-05-25 2017-08-29 国网浙江省电力公司电力科学研究院 一种gis特高频局部放电在线监测装置的放电源定位方法
CN107655533A (zh) * 2017-10-31 2018-02-02 合肥工业大学 一种基于回波能量积分的气体超声波流量计信号处理方法和系统
CN109828187A (zh) * 2018-11-19 2019-05-31 国网浙江省电力有限公司电力科学研究院 一种定位变电设备局部放电故障的时间差提取分析方法
CN110658418A (zh) * 2019-09-30 2020-01-07 山东信通电子股份有限公司 一种线缆故障检测方法和装置
CN110940933A (zh) * 2019-12-02 2020-03-31 国网河北省电力有限公司电力科学研究院 一种用于测量陡脉冲上升沿起始时刻的综合计算方法
KR102106307B1 (ko) * 2019-11-11 2020-05-04 한국건설기술연구원 에너지 데이터 처리 방법 및 그 장치
CN113639208A (zh) * 2021-07-16 2021-11-12 汕头大学 一种基于极限逼近的负压波信号拐点定位方法及系统
CN114386324A (zh) * 2021-12-27 2022-04-22 上海电力大学 一种基于转折性时段识别的超短期风电功率分段预测方法

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6509740B1 (en) * 1999-04-09 2003-01-21 Jovial Test Equipment, Inc. Methods and apparatus for time domain reflectometry
JP2009216587A (ja) * 2008-03-11 2009-09-24 Central Res Inst Of Electric Power Ind 送電線故障箇所標定方法、装置並びにプログラム
CN103323032A (zh) * 2013-05-17 2013-09-25 厦门科华恒盛股份有限公司 带死区单元双滑动窗滤波的阶跃信号检测方法
CN105223481A (zh) * 2015-10-22 2016-01-06 西安交通大学 基于差值能量函数的局部放电特高频信号起始时刻确定方法
CN105319487A (zh) * 2015-10-30 2016-02-10 国家电网公司 变电站局部放电信号检测和定位系统及方法
CN105574284A (zh) * 2015-12-29 2016-05-11 山东鲁能软件技术有限公司 一种基于趋势特征点的电力设备故障诊断方法
CN107102244A (zh) * 2017-05-25 2017-08-29 国网浙江省电力公司电力科学研究院 一种gis特高频局部放电在线监测装置的放电源定位方法
CN107655533A (zh) * 2017-10-31 2018-02-02 合肥工业大学 一种基于回波能量积分的气体超声波流量计信号处理方法和系统
CN109828187A (zh) * 2018-11-19 2019-05-31 国网浙江省电力有限公司电力科学研究院 一种定位变电设备局部放电故障的时间差提取分析方法
CN110658418A (zh) * 2019-09-30 2020-01-07 山东信通电子股份有限公司 一种线缆故障检测方法和装置
KR102106307B1 (ko) * 2019-11-11 2020-05-04 한국건설기술연구원 에너지 데이터 처리 방법 및 그 장치
CN110940933A (zh) * 2019-12-02 2020-03-31 国网河北省电力有限公司电力科学研究院 一种用于测量陡脉冲上升沿起始时刻的综合计算方法
CN113639208A (zh) * 2021-07-16 2021-11-12 汕头大学 一种基于极限逼近的负压波信号拐点定位方法及系统
CN114386324A (zh) * 2021-12-27 2022-04-22 上海电力大学 一种基于转折性时段识别的超短期风电功率分段预测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Selective review of offline change point detection methods;Charles Truong;《Signal Processing》;20190919(第167期);第1-20页 *
最小能量法因子对RSO 法检测匝间短路故障的时间差 定位影响分析;李天辉;《高压电器》;20190516;第55卷(第5期);第54-61页 *

Also Published As

Publication number Publication date
CN115078909A (zh) 2022-09-20

Similar Documents

Publication Publication Date Title
CN101509972B (zh) 基于高分辨目标距离像修正相关矩阵的宽带雷达检测方法
CN103199944B (zh) 广播式自动相关监视信号检测方法及装置
CN107632236B (zh) 一种基于对端母线反射波识别的单出线输电线路单端故障测距方法
CN102263601B (zh) 一种宽带多信号检测方法
CN210923925U (zh) 一种leu的c接口线缆故障检测的装置
CN110514963B (zh) 一种改进双端行波故障测距方法
CN102377452A (zh) 高速取样有限精度量化脉冲超宽带信号到达时间估计方法
CN108344982A (zh) 基于长时间相参积累的小型无人机目标雷达检测方法
CN104406510A (zh) 一种基于经验包络法ee的电缆长度测量方法
CN103743435A (zh) 一种多传感器数据融合方法
CN105223481A (zh) 基于差值能量函数的局部放电特高频信号起始时刻确定方法
CN104122538A (zh) 确定风廓线雷达噪声功率的方法
CN107179473B (zh) 一种输电线路故障定位方法
CN115097253A (zh) 一种mmc-hvdc直流输电线路故障测距方法及系统
CN104457643A (zh) 一种轨道几何检测数据的脉冲噪声滤波方法及装置
CN115078909B (zh) 一种基于能量拐点的非理想脉冲前沿检测方法及装置
CN109901016B (zh) 一种消除波速影响的双端测距方法
CN108445477A (zh) 机场场面异物检测雷达的高精度测距方法
CN108828398B (zh) 一种基于结构相似性算法的单芯电缆局放定位方法
CN106154257A (zh) 基于FFT与apFFT的精密测量雷达二次测频方法
CN103441965A (zh) 一种Sounding信号的信噪比估计方法
CN102033229B (zh) 一种合成孔径雷达信号跟踪方法
CN109696603B (zh) 一种双端行波测距的两侧外部时钟实时监测方法
CN106209278A (zh) 一种无人直升机旋翼缝隙检测方法
CN104316929B (zh) 一种功率谱识别方法及装置

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
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20240208

Address after: 223800 no.2481 Fazhan Avenue, Sucheng District, Suqian City, Jiangsu Province

Applicant after: STATE GRID JIANGSU ELECTRIC POWER Co.,Ltd. SUQIAN POWER SUPPLY BRANCH

Country or region after: China

Applicant after: SUQIAN ELECTRIC POWER DESIGN INSTITUTE Co.,Ltd.

Address before: 223800 no.2481 Fazhan Avenue, Sucheng District, Suqian City, Jiangsu Province

Applicant before: STATE GRID JIANGSU ELECTRIC POWER Co.,Ltd. SUQIAN POWER SUPPLY BRANCH

Country or region before: China

GR01 Patent grant
GR01 Patent grant