基于粒子群优化的弱暂态零序电流故障特征提取方法技术

本发明专利技术涉及一种基于粒子群优化的弱暂态零序电流故障特征提取方法。首先，设置势函数参数为优化对象、随机初始化粒子组成一个种群；然后，利用变尺度双稳态系统对弱暂态零序电流进行特征提取，计算初始电流与该电流之间的互相关系数；接着，按照互相关系数的最大化作为暂态零序电流检测效果的度量准则和参数优化选取的依据，根据适应度值更新粒子位置和飞行速度；最后，当进化次数达到最大阈值，停止更新并输出最优化参数，在该条件下，对暂态零序电流进行特征提取，若互相关系数大于0.85，则故障特征提取完成；反之，则再进行优化，直至大于0.85，故障特征提取完成。

【技术实现步骤摘要】

本专利技术涉及基于粒子群优化的弱暂态零序电流故障特征提取方法，属于电力系统配电网故障选线领域。
技术介绍
配电网故障选线具有以下难点和问题：1)信号的故障特征不明显：单相接地故障后，稳态电流一般小于30A甚至只有几A，此外，配电网络结构复杂导致有时候故障特征不明显，虽然故障暂态零序电流信号比稳态零序电流信号大，但是持续时间短，有时难以检测；2)我国配电网运行方式多变，各配电线路的长短不一、数量也会经常发生变化，其线路的谐波电流和分布电容电流也随之发生变。另外，外界噪声的强度、负荷的影响、母线电压的波动和故障点接地电阻的不确定等因素均会影响故障零序电流的变化。综上所述，如何在强随机噪声背景下提取暂态零序电流的故障特征是解决配电网故障选线的关键技术。随机共振技术在上世纪八十年代由意大利学者Benzi等人在研究地球古气象冰川问题时提出的。它指一个非线性双稳态系统，当仅在噪声或仅在小周期信号作用下都不足以使系统输出在两个稳态之间跃迁，而在噪声和小周期信号的共同作用下，系统输出的功率谱中，在信号的频率处出现一峰值，当噪声强度达到某一合适值时，输出功率谱的峰值达到最大。随机共振利用噪声增强微弱信号传输的优点，使其与其他的微弱信号检测方法相比具有独特的优势。然而，现有文献只是根据经验选取势函数参数，有可能导致随机共振提取的特征信号不够准确，大大影响了随机共振的应用效果。因此，如何选取合理的势函数参数，最有效地利用随机共振来增强强噪声背景下弱暂态零序电流信号的检测，是配电网故障诊断领域需要解决的问题之一。
技术实现思路
本专利技术的目的在于利用最优化参数条件下的变尺度双稳态系统提取强噪声背景下弱暂态零序电流的故障特征，为了达到上述目的，本专利技术提出的技术方案为：基于粒子群优化的弱暂态零序电流故障特征提取方法，包括如下步骤：步骤1设定双稳态系统中的势函数参数a、b为优化对象，其中，双稳态系统表达式为：dx/dt＝-dV(x)/dx+s(t)V(x)为势函数，V(x)＝-ax2/2+bx4/4，a和b为势函数参数，s(t)代表暂态零序电流iz(t)，其中，iz(t)不含噪声信号；步骤2设定粒子的维数为2维，其中，a代表第1维，b代表第2维；步骤3设定粒子种群规模及进化次数，其中，粒子个数为25个，进化次数为100次；步骤4初始化粒子种群中各粒子的位置Xj(k)及飞行速度Vj(k)，其中，Xj(k)和Vj(k)分别表示k时刻种群中第j个粒子的位置和速度，j为粒子编号；步骤5将强噪声背景下的弱暂态零序电流izg(t)替换双稳态系统表达式中的s(t)，并对双稳态系统进行四阶龙格-库塔算法求解，求解所得的解为初始电流is(t)，其中，izg(t)为对暂态零序电流iz(t)添加强噪声所得的含有噪声的暂态零序电流；步骤6设定个体适应度评价指标，其中，计算is(t)与iz(t)之间的互相关系数ρsz，并将ρsz设置为个体适应度评价指标；步骤7更新所有粒子的飞行速度和位置，其中，粒子飞行速度更新计算式为：Vj(k+1)＝Vj(k)+c1×R1×[Pj(k)-Xj(k)]+c2×R2×[G(k)-Xj(k)]其中，R1、R2为介于0～1之间相互独立均匀分布的随机数；c1、c2为学习因子，且c1＝c2＝1.49445；Xj(k)和Vj(k)分别表示k时刻种群中第j个粒子的位置和速度，其中，j为粒子编号；Pj(k)表示k时刻粒子本身迄今为止搜索到的最优解；G(k)表示k时刻整个种群迄今为止搜索到的最优解，其中，k表示当前进化次数；Vj(k+1)表示k+1时刻种群中第j个粒子的速度；粒子位置更新计算式为：Xj(k+1)＝Xj(k)+Vj(k+1)其中，Xj(k+1)、Vj(k+1)分别表示调整为k+1时刻后粒子的位置和粒子的飞行速度；步骤8终止条件判定，其中，当进化次数达到设定的最大阈值Tmax时，迭代停止，Tmax为100次；步骤9输出最优化势函数参数ay、by，将ay替换双稳态系统表达式中势函数的参数a，将by替换双稳态系统表达式中势函数的参数b；将强噪声背景下弱暂态零序电流izg(t)替换双稳态系统表达式中的s(t)，并对双稳态系统进行四阶龙格-库塔算法求解，求解所得的解为特征电流ic(t)；步骤10计算ic(t)与iz(t)之间的互相关系数ρcz的数值，若ρcz＞0.85，则暂态零序电流故障特征提取完成；若ρcz≤0.85，则转入步骤4，直至ρcz＞0.85，故障特征提取结束。本专利技术工作原理1暂态零序电流基于暂态量选线法所依据的暂态特征比稳态值大几倍甚至几十倍，且不受消弧线圈的影响，无需添加额外设备，因此具有更高的可靠性及应用价值。因此，利用单相接地零序暂态等值电路进行暂态零序电流分析，如图1所示。其中：C0为线路零序电容；L0为线路零序等值电感；Rg为接地点的过渡电阻；Rp和Lp分别为消弧线圈的等效电阻和电感；e(t)为零序电压。在补偿电网发生故障的瞬间，由图1可得流过故障点的暂态零序电流iz(t)为：其中：i0L.t为暂态零序电流中的电感电流分量；i0C.t为暂态零序电流的电容电流分量；ILm和ICm分别为电感电流和电容电流的初值(ICm＝UphmωC，ILm＝Uphm/ωL)；Uphm为相电压的幅值；ω为工频角频率；ωf和δ分别为暂态零序电流容性分量的振荡角频率和衰减系数；τL为电感电流的衰减时间常数；为接地时故障线路相电压的初始相位。由式(1)可知，当小电流接地系统发生单相接地故障时，暂态电容电流具有周期性的衰减振荡特性。此外，一般架空线路的自由振荡频率为300～1500Hz，电缆线路的电感远小于架空线路，而对地电容却较后者大许多倍，故电容电流暂态过程的振荡频率很高，持续时间很短，其自由振荡频率一般为1500Hz～3000Hz。2.双稳态系统用于研究随机共振的双稳态系统：dx/dt＝-dU(x)/dx+s(t)+Γ(t)(2)其中，U(x)＝-ax2/2+bx4/4，s(t)代表输入信号，本专利技术中s(t)代表暂态零序电流iz(t)，Γ(t)代表噪声，本专利技术将a＝b＝1称为经验参数。尽管随机共振现象与人的直觉有差异，但它的基本原理还是比较简单的。首先，可通过一个简单的模型来解释随机共振的基本原理，示意图如图2所示。单个布朗粒子在图2所示的对称双势阱中运动，当本文档来自技高网...

【技术保护点】
一种基于粒子群优化的弱暂态零序电流故障特征提取方法，其特征在于，所述弱暂态零序电流故障特征提取方法包括如下步骤：步骤1设定双稳态系统中的势函数参数a、b为优化对象，其中，双稳态系统表达式为：dx/dt＝‑dV(x)/dx+s(t)V(x)为势函数，V(x)＝‑ax2/2+bx4/4，a和b为势函数参数，s(t)代表暂态零序电流iz(t)，其中，iz(t)不含噪声信号；步骤2设定粒子的维数为2维，其中，a代表第1维，b代表第2维；步骤3设定粒子种群规模及进化次数，其中，粒子个数为25个，进化次数为100次；步骤4初始化粒子种群中各粒子的位置Xj(k)及飞行速度Vj(k)，其中，Xj(k)和Vj(k)分别表示k时刻种群中第j个粒子的位置和速度，j为粒子编号；步骤5将强噪声背景下的弱暂态零序电流izg(t)替换双稳态系统表达式中的s(t)，并对双稳态系统进行四阶龙格‑库塔算法求解，求解所得的解为初始电流is(t)，其中，izg(t)为对暂态零序电流iz(t)添加强噪声所得的含有噪声的暂态零序电流；步骤6设定个体适应度评价指标，其中，计算is(t)与iz(t)之间的互相关系数ρsz，并将ρsz设置为个体适应度评价指标；步骤7更新所有粒子的飞行速度和位置，其中，粒子飞行速度更新计算式为：Vj(k+1)＝Vj(k)+c1×R1×[Pj(k)‑Xj(k)]+c2×R2×[G(k)‑Xj(k)]其中，R1、R2为介于0～1之间相互独立均匀分布的随机数；c1、c2为学习因子，且c1＝c2＝1.49445；Xj(k)和Vj(k)分别表示k时刻种群中第j个粒子的位置和速度，其中，j为粒子编号；Pj(k)表示k时刻粒子本身迄今为止搜索到的最优解；G(k)表示k时刻整个种群迄今为止搜索到的最优解，其中，k表示当前进化次数；Vj(k+1)表示k+1时刻种群中第j个粒子的速度；粒子位置更新计算式为：Xj(k+1)＝Xj(k)+Vj(k+1)其中，Xj(k+1)、Vj(k+1)分别表示调整为k+1时刻后粒子的位置和粒子的飞行速度；步骤8终止条件判定，其中，当进化次数达到设定的最大阈值Tmax时，迭代停止，Tmax为100次；步骤9输出最优化势函数参数ay、by，将ay替换双稳态系统表达式中势函数的参数a，将by替换双稳态系统表达式中势函数的参数b；将强噪声背景下弱暂态零序电流izg(t)替换双稳态系统表达式中的s(t)，并对双稳态系统进行四阶龙格‑库塔算法求解，求解所得的解为特征电流ic(t)；步骤10计算ic(t)与iz(t)之间的互相关系数ρcz的数值，若ρcz＞0.85，则暂态零序电流故障特征提取完成；若ρcz≤0.85，则转入步骤4，直至ρcz＞0.85，故障特征提取结束。...

【技术特征摘要】
1.一种基于粒子群优化的弱暂态零序电流故障特征提取方法，其特征
在于，所述弱暂态零序电流故障特征提取方法包括如下步骤：
步骤1设定双稳态系统中的势函数参数a、b为优化对象，其中，双稳
态系统表达式为：
dx/dt＝-dV(x)/dx+s(t)
V(x)为势函数，V(x)＝-ax2/2+bx4/4，a和b为势函数参数，s(t)代
表暂态零序电流iz(t)，其中，iz(t)不含噪声信号；
步骤2设定粒子的维数为2维，其中，a代表第1维，b代表第2维；
步骤3设定粒子种群规模及进化次数，其中，粒子个数为25个，进化
次数为100次；
步骤4初始化粒子种群中各粒子的位置Xj(k)及飞行速度Vj(k)，其中，
Xj(k)和Vj(k)分别表示k时刻种群中第j个粒子的位置和速度，j为粒子编
号；
步骤5将强噪声背景下的弱暂态零序电流izg(t)替换双稳态系统表达式
中的s(t)，并对双稳态系统进行四阶龙格-库塔算法求解，求解所得的解为
初始电流is(t)，其中，izg(t)为对暂态零序电流iz(t)添加强噪声所得的含有
噪声的暂态零序电流；
步骤6设定个体适应度评价指标，其中，计算is(t)与iz(t)之间的互相
关系数ρsz，并将ρsz设置为个体适应度评价指标；
步骤7更新所有粒子的飞行速度和位置，其中，粒子飞行速度更新计
算式为：
Vj(k+1)＝Vj(k)+c1×R1×[...

【专利技术属性】
技术研发人员：王晓卫，高杰，魏向向，韦延方，曾志辉，
申请(专利权)人：河南理工大学，
类型：发明
国别省市：河南;41