基于差值能量函数的局部放电特高频信号起始时刻确定方法
【技术领域】：
[0001] 本发明涉及电力设备局部放电检测领域，尤其涉及一种基于差值能量函数的局部 放电特高频信号起始时刻确定方法。
【背景技术】：
[0002] 局部放电是导致电力设备发生重大事故的主要原因之一，快速准确定位局部放电 源对电力设备的维护和保证电力系统安全运行至关重要。局部放电监测对于电力设备的安 全运行有着极其重要的意义。由于特高频法具有检测灵敏度高、信噪比高及可实现局部放 电源定位等优点，其在电力设备局部放电检测中获得广泛应用，其中对于局部放电源的定 位更是研究的热点之一。
[0003] 由于气体绝缘组合电器、变压器等大型设备结构复杂，同时现场干扰强烈，以及传 播过程会使信号产生畸变，这些因素都严重影响对放电源的精确定位。利用特高频法进行 局部放电定位的前提是准确确定局部放电产生的特高频信号的起始时刻。因此，确定局部 放电特高频信号的起始时刻是局部放电定位技术的基础，也是决定定位准确度的关键。目 前的时差确定方法主要有阈值法、累积能量法和相关估计法等。阈值法即首先根据背景噪 声水平合理选定信号的阈值，以特高频传感器接收的信号达到阈值的时刻作为信号起始时 亥IJ。累积能量法即寻找接收信号的累积能量曲线的拐点作为信号起始时刻，将多个不同位 置的特高频传感器接收到信号起始时刻的差值作为信号时差的估计值。发明人发现，以上 几种方法均存在不足之处。其不足之处体现在：（1)要求特高频信号具有较高的信噪比。当 特高频信号的噪声较大时，确定的信号时差有较大误差；(2)当现场存在干扰使特高频信 号波形发生畸变时，应用以上方法确定特高频信号的时差也会产生较大误差。
[0004] 因此，亟需一种方法来确定局部放电特高频信号的起始时刻，并且受噪声的干扰 较小，计算结果精确度较高。
【发明内容】
：
[0005] 本发明的目的在于提供一种基于差值能量函数法的局部放电特高频信号起始时 刻确定方法。
[0006] 为达到上述目的，本发明采用如下技术方案来实现的：
[0007] 基于差值能量函数的局部放电特高频信号起始时刻确定方法，包括如下步骤：
[0008] 1)在被测试设备安装特高频传感器，采用高采样率数字检测设备采集并记录局部 放电产生的特高频信号波形；
[0009] 2)计算局部放电产生的特高频信号波形的累积能量函数；
[0010] 3)优化选取数据窗口的长度WL，从第一个数据窗开始滑动计算特高频信号波形 累积能量函数的差值能量函数；
[0011] 4)确定差值能量函数的首个显著峰值对应的时刻，将其作为特高频信号的起始时 刻。
[0012] 本发明进一步的改进在于，步骤1)中所述特高频传感器，为超宽带特高频信号传 感器，数量为四个或者以上，能够建立多个时差定位方程组，通过求解时差方程组确定局部 放电源的位置。
[0013] 本发明进一步的改进在于，步骤1)中所述采用高采样率数字检测设备，采用的采 样率应不低于5GS/s，且随着采样率的增加起始时刻精度会逐渐提高。
[0014] 本发明进一步的改进在于，步骤2)中所述累积能量函数，具体计算方法为：
[0015] 假定采集的局部放电特高频信号波形为VaiKi= 1，2，一N，N为单次脉冲的采样 点数，h为第i个采样点对应的时刻，则累积能量函数U(tk)计算公式如式⑴：
[0017] 本发明进一步的改进在于，步骤3)中所述的差值能量函数，为一定数据窗口长度 内累积能量最大值与最小值之差，具体计算公式为：
[0018] du(i) =max{Ui+m}-min{Ui+m}m=-WL/2,…，WL/2 (2)
[0019] 式中，du(i)为第i个数据窗内的差值能量值，U1为第i个采样点对应的累积能量 数值，WL为滑动数据窗口的长度。
[0020] 本发明进一步的改进在于，步骤3)中所述的优化选取数据窗口的长度，包括以下 步骤：
[0021] (a)根据步骤1)检测并记录20组特高频信号波形；
[0022] (b)针对测到的每组特高频信号，利用蒙特卡罗统计法分析时差受噪声影响程度： 在测量到的局部放电特高频信号波形的数据中加入指定信噪比的白噪声，利用差值能量函 数法确定加入白噪声后各传感器之间的起始时刻，进而计算天线信号到达时刻的时差；利 用蒙特卡罗法重复以上加入白噪声-确定传感器信号时差过程N次，最后计算N次时差统 计的平均值及标准差；
[0023] (c)改变窗口长度，重复步骤（b)，以计算得到的时差平均值与真实时间差误差最 小、标准差最小为目标，选取最优的窗口长度；
[0024] ⑷对采集的20组特高频信号，重复步骤（b)~（c)，得到各组信号优选的数据窗 口长度；综合比较20组信号，以评估的时差平均值与真实时间差误差最小、标准差最小为 目标，选取最优的窗口长度。
[0025] 本发明进一步的改进在于，步骤4)所述的确定差值能量函数的首个显著峰值对 应的时刻，包括以下步骤：
[0026] (a)对特高频信号差值能量函数进行以下计算：计算其最大值M1，求取差值能量函 数为0.SM1对应的时刻tpt2=ti-lOns时刻之前一般仅包含噪声，提取该部分噪声信号， 计算噪声部分数学形态学梯度的最大值Vn;
[0027] (b)以1.5Vn为阈值，提取差值能量第一个大于该阈值的峰值对应的时刻，将其作 为局部放电特高频信号的起始时刻。
[0028] 相对于现有技术，本发明的有益效果体现在：
[0029] 本发明提供了一种基于差值能量函数的局部放电特高频信号起始时刻确定方法。 差值能量函数可强化累积能量的快速增长位置，减弱累积能量的缓慢增长部分。特高频信 号中的白噪声与低频（〈1GHz)窄带噪声对应累积能量的缓慢增长位置，因此差值能量函数 可有效减低白噪声与窄带白噪声水平。
[0030] 进一步，由于差值能量函数可有效降低噪声水平，使用本发明确定局部放电特高 频信号的起始时刻，相对于传统方法具有较高的精确度，并且受外部噪声的干扰较小，特别 在现场局部放电检测时可有效减低噪声的影响。
[0031] 进一步，本发明优化选取数据窗口的长度步骤中提出了分析噪声对时差确定算法 影响的方法，该方法还可用于对不同时差确定算法受噪声影响程度的评估。
[0032] 进一步，本发明可以应用在电力设备检测领域，对提高局部放电源的定位精度有 较大的实用价值。
【附图说明】：
[0033] 图1为本发明基于差值能量函数的局部放电特高频信号起始时刻确定方法的流 程图。
[0034] 图2为本发明实验中采用的特高频传感器布置方式图。
[0035] 图3为本发明四个特高频传感器测得的局部放电特高频信号波形图。
[0036] 图4为本发明传感器#1信号的累积能量函数与差值能量函数波形图。
[0037] 图5为本发明不同数据窗口长度下传感器#1差值能量函数波形图。
[0038] 图6为本发明采用蒙特卡罗法计算时差平均值与标准差的流程图。
[0039] 图7为本发明计算得到的不同信噪比下传感器#1与#2信号时差平均值与标准 差；其中，图7 (a)、7 (b)分别为时差平均值与标准差变化趋势图。
[0040] 图8为本发明利用差值能量函数确定起始时刻实例图。
【具体实施方式】：
[0041] 下面结合附图和实施例对本发明的【具体实施方式】进行说明。
[0042] 如图1所示，本发明基于差值能量函数的局部放电特高频信号起始时刻确定方 法，包括如下步骤：
[0043] (1)在被测试设备安装特高频传感器，采用高采样率数字检测设备采集并记录局 部放电产生的特高频信号波形；
[0044] 作为一种实施例，图2为本发明搭建的一种基于特高频传感器阵列的变电站局部 放电检测实验系统，该系统由4个全向传感器与数字采集设备组成，4个传感器位于矩形的 四个顶点；利用便携式特高频信号发生器向空间中辐射特高频信号，传感器为盘锥形全向 传感器，其带宽为300MHz~3GHz。利用高速采集数字示波器以10GS/s的采样率检测采集 特高频信号，检测得到4个传感器的特高频信号波形如图3所示。
[0045] (2)计算局部放电产生的特高频信号波形的累积能量函数；
[0046]本步骤中所述的累积能量函数，具体计算方法为：假定采集的局部放电特高频信 号波形为VU1)，i= 1，2,…N，N为单次脉冲的采样点数，h为第i个采样点对应的时刻， 则累积能量函数U(tk)计算公式为：
[0048] 作为一种实施例，计算了传感器#1的累积能量函数，如图4所示。
[0049] (3)优化选取数据窗口的长度WL，从第一个数据窗开始滑动计算特高频信号波形 累积能量函数的差值能量函数；
[0050] 本步骤中所述的差值能量函数，为一定数据窗口长度内累积能量最大值与最小值 之差，具体计算公式为：
[0051] du(i) =max{Ui+m}-min{Ui+m}m=-WL/2,…，WL/2 (2)
[0052] 式中，du(i)为第i个数据窗内的差值能量值，U1为第i个采样点对应的累积能量 数值，WL为滑动数据窗口的长度。
[0053] 作为一种实施例，利用公式（2)计算了窗口长度为2、8及14个采样点时传感器 #1特高频信号的差值能量函数，