技术领域
[0001] 本
发明涉及输电线路故障设备检测技术领域，具体是一种基于波权重的输电线路故障电流行波波形的重现方法。
背景技术
[0002] 输电线路是电
力系统中最易发生故障的设备，对于输电线路的故障测距和故障
定位，学者进行了诸多研究。对于输电线路的故障而言，只可被动增加其防御能力以减少输电线路发生故障的概率。而在实际线路局部改造工程应用时，由于缺乏故障数据作为
支撑，使得改造工作无法做到有的放矢，从而浪费大量的人力物力投资。若能重演故障的发展过程，实现输电线路故障整个过程的状态情景可视，对输电线路故障分析、诊断、后期的线路防御方案的制定及故障事故演习等方面都有重大的现实意义。
[0003] 对于单相线路，若设该段输电线路均匀、统一，其单位长度的
电阻、电感、电容和电导分别为R、L、G、C，输电线路上取一段长度为dx的线路，该线路在
频率域上的传播方程为：
[0004]
[0005] 式(1)的解最终可以写成如下形式:
[0006]
[0007] 其中， 为线路传播系数，x为传播距离，Zc为波阻抗。A1、A2为由边界条件确定的积分常数。
[0008] 对于单个故障行波，若不需考虑沿x反方向传播的反行波，式(2)可以写为：
[0009]
[0010] 由式(3)可得，对于线路上相距为x的两点1和2，他们之间的电流波和
电压波有如下关系
[0011]
[0012] 由式(4)可以看出，波在单相均匀输电线路的传播与其频率密切相关，且随着传播距离x的增加，电压、电流逐渐衰减。记H＝e-λx为输电线路频率域沿线传播的传递函数，则该传播模型称为输电线路频依函数模型。
[0013] 而对于三相输电线路来说，三相间存在耦合关系，需通过相模变换变成三个独立的分量，以便于分析。
[0014] 对于线路上相距为d的两点1和2，由(4)式经过相模变换后最终得到如下电流行波波形关系
[0015]
[0016] 其中，上标i(i＝0,1,2)表示i模分量。
[0017] 若想知道输电线路的传输规律，λ中的四个分布参数必须已知，而从式(5)可以看到，在传递函数中λ有四个变量且相互耦合影响，传统的数学方法对于这类方程的求解已经变得力不从心，且很难得到全局最优解。故需要借用全局优化智能
算法——模拟
退火算法。
[0018] 模拟退火算法(Simulated Annealing，SA)最早的思想是由N.Metropolis等人于1953年提出。1983年,S.Kirkpatrick等成功地将退火思想引入到组合优化领域。它是基于Monte-Carlo
迭代求解策略的一种随机寻优算法，其出发点是基于物理中固体物质的退火过程与一般组合优化问题之间的相似性。模拟退火算法从某一较高初温出发，伴随
温度参数的不断下降,结合概率突跳特性在解空间中随机寻找目标函数的全局最优解，即在局部最优解能概率性地跳出并最终趋于全局最优。
[0019] 在故障重演方面，国内外学者更多的关注系统层面的故障重现，如华北电力大学张东英依靠
电网一次模型、电网实时运行、故障录波信息、保护动作信息及故障过程中
开关动作等信息先构建故障区域确定可疑设备，然后根据可疑设备主动搜集与可疑设备相关的信息，综合利用故障录波结果作为中间结论，最后利用证据理论、正反混合推理等方法最终确定故障设备和对故障过程进行初步判断。但是对于输电线路故障行波波形信息的更多利用和
数据挖掘并未有充分挖掘。
发明内容
[0020] 为了解决故障点信息重演这
块技术的缺失，本发明旨在提供一种基于波权重的输电线路故障电流行波波形的重现方法，以故障沿线检测点检测到的行
波数据为依据，通过智能算法模拟退火算法反演出输电线路上的分布参数，再结合输电线路频依模型反演出线路故障点的行波波形。
[0021] 本发明的技术解决方案如下：
[0022] 一种基于波权重的输电线路故障电流行波波形的重现方法，该方法包括以下步骤：
[0023] 步骤S1：在输电线路上依次设定多个故障检测点，用于故障电流行波波形的采集；
[0024] 步骤S2：对输电线路上距离故障实际发生处同侧除最近一个检测点的另两个故障检测点的波形进行卡伦布尔变换三相解耦；
[0025] 步骤S3：对经三相解耦后的行波各模量进行小波包变换后进行快速傅里叶变换，分频率段代入线路频依函数模型，具体步骤：
[0026] 步骤S3.1：以N层小波包变换为例，对于N层小波包变换后将两个检测点各分成的2N个频段中的每一个频段波形，进行快速傅里叶变换生成各自2N段频域数据。
[0027] 步骤S3.2：将靠近故障点的那处检测点的各模量的频域数据乘以线路依频函数H，并与另一检测点对应模量的对应频段的数据做差得Δd(i)，i表示第i个频段；
[0028] 其中H＝e-λx，x是两个检测点之间的距离， 对于每个频段，f为该频段的中心频率，R是电阻，L是电感，G是电导，C是电容。
[0029] 步骤S4：对各个频率段通过靠近故障点的检测波形的波形
能量占比作为系数进行加权组合生成最终的目标寻优函数；
[0030] 步骤S4.1：以N层小波包变换为例，对靠近故障点出的检测点的解耦波形进行卡伦布尔变换解耦。
[0031] 步骤S4.2：对解耦后的各模量波形进行N层小波包变换分成2N个频段。
[0032] 步骤S4.3：各个频段的分量所拥有的波形能量占总波形能量的百分比就是加权组合的权重a(i)，i表示第i个频段；
[0033] 步骤S4.4：最终加权后的目标函数如下： *表示共轭转置。
[0034] 步骤S5：通过模拟退火算法对目标寻优函数进行全局求解，得到该段输电线路的分布参数；
[0035] 步骤S6：将得到的分布参数代入线路频依函数模型，根据已知点波形结合线路频依参数模型反演出未知点波形，具体步骤：
[0036] S6.1：在S5中得到线路分布参数之后，确定线路频依函数模型H，其中H＝e-λx，x是行波传输的距离，将用于反演的已知点波形进行卡伦布尔变换解耦成为各模量波形；
[0037] S6.2：对三相解耦后的行波线模波形进行快速傅里叶变换，根据未知点与已知点的
位置关系与线路频依函数模型做运算：
[0038] 若未知点在已知点的下游，则已知点各模量波形乘以线路频依函数模型H；
[0039] 若未知点在已知点与线路故障点中间，则已知点各模量波形除以线路频依函数模型H；
[0040] S6.3：根据运算后未知点的线模波形，结合故障类型边界条件，推得未知点的零模波形；
[0041] S6.4：根据未知点的各模量波形进行卡伦布尔反变换得到未知点的三相故障电流行波波形。
[0042] 分布参数反演原理是基于使得行波传递模型理论推得的行波和观测测量得到的实际波形的残差达到最小来进行目标函数寻优，本发明所用为最小二乘目标函数：
[0043] Q＝Δd·Δd* (6)
[0044] 式中Δd＝dcal-dobs，dcal是当前传递模型对应的正演数据，dobs是实际测量到的数据，*表示共轭转置。
[0045] 而本发明分析的电流行波波形是个高频暂态
信号，在对于其细节提取过程中我们用到了小波包变换对信号进行分频段处理，从而最后目标函数是由各个频率段的分量所组合形成的，形式如下：
[0046]
[0047] 式子中a(i)表示第n-1层小波包变换后第i段的权重系数，由波形各个频段的能量占总能量的百分比确定；d为经过卡伦布尔三相解耦并快速傅里叶变换之后的波形数据。
[0048] 通过对(7)式中的目标函数用智能算法模拟退火算法进行全局寻优，得到线路上两个检测点间线路的各模量分布参数。
[0049] 求出线路上各模量的分布参数之后，传递函数H就已经确定，从而可以根据已知检测点的波形结合频依模型H计算出线路上其他点的波形并与该点的真实检测波形相对比，来验证反演参数的正确性和波形反演方法的可行性。
[0050] 本发明利用行波在输电线路中的频依传播特性构建故障电流波形的反演模型，通过对线路分布参数的求取和故障电流行波波形的反演，技术效果如下：
[0051] 1)通过构建故障电流行波波形的反演模型，结合小波包变换和快速傅里叶变换设立最小二乘目标函数，并通过智能搜索算法模拟退火算法全局寻优，能有效并准确地根据两个检测点波形求出输电线路上的分布参数。
[0052] 2)先求线路分布参数，再结合行波传递模型能根据已知点的波形精确反演出未知点的波形。
[0053] 3)能精确得到线路上各点的行波，为后期全面利用行波故障信息，进行故障分析做技术支撑。
附图说明
[0054] 图1是输电线路pscad模型示意图
[0055] 图2是检测点1、2、3电流行波波形
[0056] 图3是用于参数反演的各模量波形图
[0057] 图4是检测点1的电流行波波形对比图
具体实施方式
[0058] 下面结合附图，给出本发明的较佳
实施例，并予以详细描述。
[0059] 在PSCAD中建立如图1的输电线路模型。线路采用频依模型，F为
雷击故障点，在故障点沿线20km，40km和70km依次设立三个检测点，用于故障电流行波波形的采集。输电线路的杆塔为实际杆塔ZB1的模型。
[0060] 线路在0.2秒时在点F处遭受雷击，在pscad中用一个高频的脉冲信号进行模拟，产生沿线传播的暂态电流行波，在下游检测点1,2,3分别检测到衰减和畸变后的电流行波波形如图2所示，行波波形的
采样频率取1MHz.
[0061] 由图2可以看出，F点发生故障之后，检测点1,2,3的第一个电流行波波形波头和波尾大部都相对完好，但是由于线路行波存在折反射，故无法得到单一得到首个波的完整波形。为了使得反演更加精确，本文拟截去首个波头和第二个波头重叠部分，并对其进行卡伦布尔变换三相解耦，解耦后的各模量波形如图3所示
[0062] 对经三相解耦后的行波各模量进行5层小波包变换分成32个频段，并进行进行快速傅里叶变换，分频率段代入线路频依函数模型H。
[0063] 对各个频率段通过靠近故障点的检测波形(此处为检测点2)的波形能量占比作为系数进行加权组合生成最终的目标寻优函数如下：
[0064]
[0065] 通过模拟退火算法对目标寻优函数进行全局求解，得到该段输电线路的各模量分布参数如表1所示：
[0066] 表1检测点2,3之间的分布参数
[0067]
[0068] 针对波形反演的效果，本发明以以下5方面来进行评价。
[0069] 1)波头起始时间ts。
[0070] 2)波头上升时间 即在波的幅值从最大幅值的0.1倍上升到0.9倍所用的时间。
[0071] 3)故障行波峰值Im。
[0072] 4)峰值的位置tm。
[0073] 5)半
波长度thm＝th-ts。其中th为幅值增加到最高再下降到一半时的时刻。
[0074] 通过以上5方面的对比评价，能全方位的反映波形的各方面特性。
[0075] 由于零模分量的反演受到各方面的制约，并不能准确反演，故本发明放弃了对零模的反演，转而通过不同故障类型时零模和线模的边界条件，结合线模的反演结果得到零模分量。本例故障类型为A相接地，通过边界条件及相模关系可以求得线路反演零模分量。再对各模量波形进行卡伦布尔反变换得到
三相电流波形对比图如图4所示。
[0076] 对图4中第一个波头的数据，本发明的5个评价体系评价如下：
[0077] 1)波头起始时间ts。三相反演波形和检测波形绝对误差分别为A相0.5μs，B相1μs，C相2μs。
[0078] 2)波头上升时间 A、B、C三相的相对误差分别为1.698％，3.425％，4.875％。
[0079] 3)故障行波峰值Im。A、B、C三相的相对误差分别为1.774％，1.852％，4.691％。
[0080] 4)峰值的位置tm。A、B、C三相的绝对误差分别2μs，1μs，1μs。
[0081] 5)半波长度thm＝th-ts。A、B、C三相的相对误差分别为0.723％，2.820％，4.669％。
[0082] 经过大量仿真实验证明，该故障电流行波波形反演方法对五项评价指标绝对误差不超过3μs，相对误差不超过6％。
[0083] 要注意的是，以上列举的仅为本发明的具体实施例，显然本发明不限于以上实施例，随之有着许多的类似变化。本领域的技术人员如果从本发明公开的内容直接导出或联想到的所有
变形，均应属于本发明的保护范围。
