CN105260777A - 一种多参数优化的弱暂态零序电流故障特征提取方法 - Google Patents
一种多参数优化的弱暂态零序电流故障特征提取方法 Download PDFInfo
- Publication number
- CN105260777A CN105260777A CN201510615562.1A CN201510615562A CN105260777A CN 105260777 A CN105260777 A CN 105260777A CN 201510615562 A CN201510615562 A CN 201510615562A CN 105260777 A CN105260777 A CN 105260777A
- Authority
- CN
- China
- Prior art keywords
- bistable system
- sequence current
- transient zero
- optimization
- parameter
- 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.)
- Granted
Links
Landscapes
- Complex Calculations (AREA)
Abstract
本发明涉及一种多参数优化的弱暂态零序电流故障特征提取方法,包括:步骤S1:设定双稳态系统中的势函数参数a、b和四阶龙格-库塔算法求解双稳态系统时的数值计算步长h为优化对象;步骤S2:采用遗传算法对参数a、b和h进行优化,输出最优化参数aJ、bJ和hJ;步骤S3:在最优化参数条件下,利用变尺度方式下的双稳态系统对噪声背景下的暂态零序电流进行特征提取,得到变尺度方式下的双稳态系统的输出信号ic(t)。与现有技术相比,本发明利用遗传算法的全局最优搜索性能以及多参数同时优化的思想,寻找最优化势函数参数和数值计算步长,以便双稳态系统在最优化参数条件下提取强噪声背景下微弱暂态零序电流的故障特征,有利于提高配电网故障选线的准确度。
Description
技术领域
本发明涉及电力系统配电网故障选线领域,尤其是涉及一种多参数优化的弱暂态零序电流故障特征提取方法。
背景技术
配电网故障选线具有以下难点和问题:1)信号的故障特征不明显:单相接地故障后,稳态电流一般小于30A甚至只有几A,此外,配电网络结构复杂导致有时候故障特征不明显,虽然故障暂态零序电流信号比稳态零序电流信号大,但是持续时间短,有时难以检测;2)我国配电网运行方式多变,各配电线路的长短不一、数量也会经常发生变化,其线路的谐波电流和分布电容电流也随之发生变。另外,外界噪声的强度、负荷的影响、母线电压的波动和故障点接地电阻的不确定等因素均会影响故障零序电流的变化。综上所述,如何在强随机噪声背景下提取微弱暂态零序电流(TransientZero-SequenceCurrent,TZSC)的故障特征是解决配电网故障选线的关键技术。
所谓随机共振指一个非线性双稳态系统,当仅在噪声或仅在小周期信号作用下都不足以使系统输出在两个稳态之间跃迁,而在噪声和小周期信号的共同作用下,系统输出的功率谱中,在信号的频率处出现一峰值,当噪声强度达到某一合适值时,输出功率谱的峰值达到最大。随机共振利用噪声增强微弱信号传输的优点,使其与其他的微弱信号检测方法相比具有独特的优势。然而,现有文献只是根据经验选取势函数参数,有可能导致随机共振提取的特征信号不够准确,大大影响了随机共振的应用效果。此外,在求解双稳态系统时,四阶龙格—库塔算法的数值计算步长也会影响随机共振的效果。因此,如何选取合理的势函数参数及数值计算步长,最有效地利用随机共振来增强噪声背景下微弱暂态零序电流信号的检测,是配电网故障诊断领域需要解决的问题之一。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种多参数优化的弱暂态零序电流故障特征提取方法,利用遗传算法的全局最优搜索性能以及多参数同时优化的思想,寻找最优化势函数参数和数值计算步长,以便双稳态系统在最优化参数条件下提取强噪声背景下微弱暂态零序电流的故障特征,有利于提高配电网故障选线的准确度。
本发明的目的可以通过以下技术方案来实现:
一种多参数优化的弱暂态零序电流故障特征提取方法包括:
步骤S1:设定双稳态系统中的势函数参数a、b和四阶龙格-库塔算法求解双稳态系统时的数值计算步长h为优化对象,其中,双稳态系统X的表达式为:
X=dx/dt=-dV(x)/dx+s(t)+Γ(t)
式中:dx/dt为双稳态系统的输出信号,V(x)为势函数,V(x)=-ax2/2+bx4/4;a和b为势函数参数;s(t)代表无噪声的暂态零序电流,Γ(t)代表噪声信号,t为时间,x为布朗粒子运动的速度;
步骤S2:采用遗传算法对双稳态系统X的表达式中势函数参数a、b和四阶龙格-库塔算法求解双稳态系统时的数值计算步长h进行优化,输出最优化参数aJ、bJ和hJ;
步骤S3:在最优化参数条件下,利用变尺度方式下的双稳态系统对噪声背景下的暂态零序电流进行特征提取,得到变尺度方式下的双稳态系统的输出信号ic(t),ic(t)即为微弱暂态零序电流故障特征提取的特征电流。
所述步骤S2包括以下步骤:
201:设定遗传算法的参数,所述遗传算法的参数包括种群规模N、交叉概率PJ、变异概率PB、优化对象a的搜索范围、优化对象b的搜索范围、优化对象h的搜索范围和最大进化次数Gmax,同时设定理想暂态零序电流iz(t),iz(t)中加入强噪声得到带有强噪声的暂态零序电流izg(t);
202:初始化种群,随机产生N个个体,其中,每个个体j对应一组优化对象组[aj,bj,hj],每组优化对象组对应一个双稳态系统Xj;
203:将izg(t)替换每个双稳态系统Xj的表达式中的s(t)+Γ(t),按个体编号顺序,依次对所有的Xj进行四阶龙格-库塔算法求解,求解双稳态系统时的计算步长记为hj,求解得到每个双稳态系统Xj的输出信号isj(t),isj(t)作为初始电流;
204:获取isj(t)与iz(t)之间的互相关系数ρsj,并将ρsj的值作为个体适应度;
205:根据个体适应度、交叉概率PJ和变异概率PB依次对N个个体进行选择操作、交叉操作和变异操作,完成一次的进化;
206:判断进化次数是否达到设定的最大进化次数,若是,则进化停止,执行步骤207,若否,则返回步骤205;
207:将ρsj数值最大的个体作为最优化个体J,J所对应的双稳态系统为最优化双稳态系统XJ,提取XJ的势函数参数aJ、bJ和求解XJ时的数值计算步长hJ,aJ、bJ和hJ即为最优化参数。
所述步骤S2还包括步骤207之后执行的步骤208,所述步骤208为:根据aJ、bJ和hJ,利用变尺度方式下的双稳态系统对izg(t)进行特征提取,得到变尺度方式下的双稳态系统的输出信号iz′(t),获取iz′(t)与iz(t)之间的互相关系数ρz′z,判断ρz′z是否大于设定阈值,若是,则执行步骤S3,若否,则跳转步骤201。
所述设定阈值为0.92。
所述优化对象a的搜索范围为[-10,10],所述优化对象b的搜索范围为[0,10],所述优化对象h的搜索范围为[0.0001,0.0219]。
所述最大进化次数Gmax的取值范围为80~100次。
所述步骤201中的强噪声的噪声强度取值为30db~500db。
与现有技术相比,本发明具有以下优点:
1)基于遗传算法优化的变尺度下双稳态系统能够根据信号和噪声的特点自适应地选取合适的势函数参数和四阶龙格-库塔算法的数值计算步长,利用最优化参数条件下的双稳态系统提取强噪声背景下微弱暂态零序电流的故障特征,获得的特征电流曲线可以很好的反应微弱暂态零序电流的状况,相比采用经验参数的势函数,本发明方法利用遗传算法的全局最优搜索性能和多参数优化思想,在选取合理的势函数参数基础上进一步选取合适的数值计算步长,可以最有效地利用随机共振来增强噪声背景下微弱暂态零序电流信号的检测,增强变尺度下双稳态系统提取弱暂态零序电流故障特征的能力,更精准的特征电流有利于提高配电网故障选线的准确度。
2)基于遗传算法的变尺度下双稳态系统能有效地检测出强噪声背景下的微弱暂态零序电流,并设定相应的设定阈值,使得最优化参数条件下的互相关系数都维持在0.9以上,通过实验,得到的与无噪声背景下的暂态零序电流之间的互相关系数高达0.9593,说明本发明方法特征提取的特征电流具有良好的抗噪性能。
附图说明
图1为单相接地零序暂态等值电路;
图2为对称双势阱中随机共振现象的示意图;
图3为遗传算法交叉操作示意图;
图4为遗传算法变异操作示意图;
图5为本发明基于多参数优化的微弱暂态零序电流故障特征提取方法流程图;
图6为无噪声的暂态零序电流示意图;
图7为强噪声背景下的暂态零序电流示意图;
图8为在经验参数下经变尺度双稳态系统所提取的特征电流;
图9为遗传算法适应度曲线;
图10为在最优化参数下经变尺度双稳态系统所提取的特征电流;
图11为在不同噪声强度下由经验参数所提取的特征电流;
其中,(11a)为噪声强度为50db时的特征电流;(11b)为噪声强度为100db时的特征电流;(11c)为噪声强度为200db的特征电流;(11d)为噪声强度为300db的特征电流;
图12为在不同噪声强度下由最优化参数所提取的特征电流;
其中,(12a)为噪声强度为50db时的特征电流;(12b)为噪声强度为100db时的特征电流;(12c)为噪声强度为200db的特征电流;(12d)为噪声强度为300db的特征电流。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
本发明从暂态零序电流的基本原理出发,涉及以下几个方面:双稳态系统的基本原理、龙格-库塔方法的基本原理、变尺度下双稳态系统的基本原理及评价指标选取和遗传算法的基本原理。
1)暂态零序电流
基于暂态量选线法所依据的暂态特征比稳态值大几倍甚至几十倍,且不受消弧线圈的影响,无需添加额外设备,因此具有更高的可靠性及应用价值。因此,利用单相接地零序暂态等值电路进行暂态零序电流分析,如图1所示。图中:C0为线路零序电容;L0为线路零序等值电感;Rg为接地点的过渡电阻;Rp和Lp分别为消弧线圈的等效电阻和电感;e(t)为零序电压。
在补偿电网发生故障的瞬间,由图1可得流过故障点的暂态零序电流iz(t)为:
式中:i0L.t为暂态零序电流中的电感电流分量;i0C.t为暂态零序电流的电容电流分量;ILm和ICm分别为电感电流和电容电流的初值(ICm=UphmωC0,ILm=Uphm/ωL0);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)代表噪声,x为布朗粒子运动的速度。本发明将a=b=1称为经验参数。
尽管随机共振现象与人的直觉有差异,但它的基本原理还是比较简单的。首先,可通过一个简单的模型来解释随机共振的基本原理,示意图如图2所示。
单个布朗粒子在图2所示的对称双势阱中运动,当它不受任何外力作用时,粒子将最终停留于其中的一个势阱内,而位于哪个势阱将由初始位置决定。但当存在随机扰动时,粒子在随机力的作用下会有一定的机率在两势阱间跳跃。当粒子仅受周期外力作用时,如果周期外力的强度很小,那么布朗粒子将在某个势阱内做小范围的振动,而不会有跨势阱的大范围运动。而当周期外力和噪声同时作用时,上述情况将会发生改变:随机力诱导的势阱间的跃迁和周期外力发生同步,粒子将以外驱动力频率在两个势阱间做大范围运动,弱的输入周期信号得以放大,于是便发生了随机共振。
3)龙格-库塔方法
为了避免计算高阶导数,龙格-库塔方法利用F(x,y)在某些点处的值的线性组合,构造一类计算公式,使其按泰勒级数展开后,与初值问题的解的泰勒展开式比较,存在尽可能多的项完全相同,从而保证算式有较高的精度,x,y表示横纵坐标。这种方法间接利用了泰勒展开的思想,避免了计算高阶导数的困难。
一般的龙格-库塔方法的形式为:
式中:αn、μn,i、cn均为待定参数,i、n为自然数;h为数值计算步长。选取这些参数的原则,是要求式(3)中第1式右端在(xn,yn)处作泰勒展开式,并按h的幂次从低到高的排列式与微分方程解的泰勒展开式 有尽可能多的项重合,也就是要求符合下面关系式:
χ1=Fn,χ2=F′n,χ3=F″n,…(4)
式中:Fn、F′n、F″n…分别为表示为y'(xn)=F(xn,yn)、y”(xn)、y”'(xn)…。
通常把式(3)称为n级龙格-库塔方法,简记为n级龙格-库塔方法。更高阶的龙格-库塔方法由于计算量较大,一般不采用。本发明采用四阶龙格-库塔方法,其计算式为:
4)变尺度下双稳态系统及评价指标
变尺度思想的本质:将大频率转换为低频率,使得大参数信号频率接近或符合随机共振所要求的小参数条件,也即将频率压缩到双稳态系统所能检测的频带范围内。变尺度的具体运算过程是:根据信号的频率和采样频率fs确定一个暂态频率压缩比CR(ContractionRatio,CR),然后根据CR再确定压缩采样频率fsc=fs/CR,由压缩采样频率fsc进一步得到压缩数值计算步长hsc=1/fsc,最后数值求解双稳系统的响应输出。
初始电流is(t)定义:在遗传算法优化过程中,对带有强噪声背景下的暂态零序电流的变尺度下双稳态系统进行四阶龙格-库塔算法求解,求解所获得输出信号即为初始电流。其中,对应第j个方案(势函数)获得的初始电流标记为isj(t)。
特征电流ic(t)定义:在最优化参数条件下,对带有未知噪声的暂态零序电流的变尺度下双稳态系统进行四阶龙格-库塔算法求解,求解所获得的解即为特征电流。
传统的信噪比测度适用于输入/输出信号的频谱中具有较清晰谱线的情况,而暂态零序电流一般都是宽带信号,信号频率范围不是集中在一个或几个可数的频率上,而是分布在很宽一段频带内。因此,传统的信噪比测度将难以适用,需要发展其它的测度指标。非线性朗之万方程虽然不能准确预测布朗粒子的运动,但是能很好地预言粒子轨道的统计性质,于是本发明利用互相关系数作为测度来描述变尺度下双稳态系统对于微弱非周期输入的响应。初始电流isj(t)和无噪声暂态零序电流iz(t)的协方差Cov(iz(t),isj(t))和互相关系数ρsj可表示为:
其中,E为数学期望,D为方差。
5)遗传算法基本原理
微弱暂态零序电流的随机共振不同于周期信号,它对系统参数的合理选取更为敏感。此外,在变尺度思想条件下,四阶龙格-库塔算法的数值计算步长h的选择也十分关键。据已有文献可知,遗传算法本质上是一种启发式的随机搜索算法,由遗传算法得出的结果每次都不尽相同,再加上自变量在给定的约束条件下进行了无缝编码,所以从理论上讲,遗传算法总有很多机会得到全局最优结果。因此,本发明采用遗传算法对双稳态系统表达式中势函数参数a、b以及四阶龙格-库塔算法的数值计算步长h进行优化,使得变尺度下双稳态系统能够更好地检测强噪声背景下的弱暂态零序电流。
遗传算法的概念是由Michigan大学的J.Holland教授在1962年提出的。遗传算法是模拟达尔文生物进化论的自然选择和孟德尔遗传学机理的生物进化过程的计算模型,是一种通过模拟自然进化过程搜索最优解的方法。它把自然界“优胜劣汰,适者生存”的生物进化原理引入优化参数形成的编码串联群体中,按照所选择的适应度函数并通过遗传中的选择、交叉和变异对个体进行筛选,使适应度值好的个体被保留,适应度差的个体被淘汰,新的群体既继承了上一代的信息,又优于上一代。这样反复循环,直至满足条件。遗传算法基本的操作分为:
选择操作:选择操作是指从旧群体中以一定概率选择个体到新群体中,个体被选中的概率跟适应度值有关,个体适应度值越好,被选中的概率越大。遗传算法选择操作有轮盘赌法、锦标赛法等多种方法,本发明选择轮盘赌法,即基于适应度比例的选择策略,每个个体j的选择概率Pj为:
oj=k1Sj(7)
式中:Sj为个体j的适应度值,本发明中令互相关系数ρsj为Sj;k1为系数;N为种群个体数目;j为个体编号。
交叉操作:交叉操作是指从个体中选择两个个体,通过两个染色体的交换组合,来产生新的优秀个体,本发明将交叉操作的概率记为PJ。交叉过程为从群体中任选两个染色体,随机选择一点或多点染色体位置进行交换。交叉操作如图3所示。由于个体采用实数编码,所以交叉操作方法采用实数交叉法,第j个染色体rj和第e个染色体re在m位的交叉操作方法为:
rjm=rjm(1-k2)+rimk2
(9)
rem=rem(1-k2)+remk2
式中:e、j为个体编号;m为染色体编号;k2为[0,1]间的随机数。
变异操作:变异操作是指从群体中任选一个个体,选择染色体中的一点进行变异以产生更优秀的个体,本发明将变异操作的概率记为PB。变异操作如图4所示。选择第j个个体的第m个基因rjm进行变异,变异操作方法为:
式中:rmax为基因rjm的上界;rmin为基因rjm的下界;f(g)=k4(1-g/Gmax)2;k4为一个随机数;g为当前迭代次数;Gmax为最大迭代次数;k3为[0,1]间的随机数。
根据以上基本原理,本发明提出一种多参数优化的弱暂态零序电流故障特征提取方法,如图5所示,包括:
步骤S1:设定双稳态系统中的势函数参数a、b和四阶龙格-库塔算法求解双稳态系统时的数值计算步长h为优化对象,其中,双稳态系统X的表达式为:
X=dx/dt=-dV(x)/dx+s(t)+Γ(t)
式中:dx/dt为双稳态系统的输出信号,V(x)为势函数,V(x)=-ax2/2+bx4/4;a和b为势函数参数;s(t)代表无噪声的暂态零序电流,Γ(t)代表噪声信号,t为时间,x为布朗粒子运动的速度。
步骤S2:采用遗传算法对双稳态系统X的表达式中势函数参数a、b和四阶龙格-库塔算法求解双稳态系统时的数值计算步长h进行优化,输出最优化参数aJ、bJ和hJ,具体包括以下步骤:
201:设定遗传算法的参数,遗传算法的参数包括种群规模N、交叉概率PJ、变异概率PB、优化对象a的搜索范围、优化对象b的搜索范围、优化对象h的搜索范围和最大进化次数Gmax,同时设定理想暂态零序电流iz(t),iz(t)中加入强噪声得到带有强噪声的暂态零序电流izg(t),加入的强噪声的噪声强度取值为30db~500db;
202:初始化种群,随机产生N个个体,其中,每个个体j对应一组优化对象组[aj,bj,hj],每组优化对象组对应一个双稳态系统Xj,j为个体编号;
203:将izg(t)替换每个双稳态系统Xj的表达式中的s(t)+Γ(t),按个体编号顺序,依次对所有的Xj进行四阶龙格-库塔算法求解,求解双稳态系统时的计算步长记为hj,求解得到每个双稳态系统Xj的输出信号isj(t),isj(t)作为初始电流;
204:获取isj(t)与iz(t)之间的互相关系数ρsj,并将ρsj的值作为个体适应度,计算种群适应度值,从中找出初始最优个体J1及初始最优适应度值这里的种群适应度值和最优个体是为了步骤205中进行选择、交叉、变异提供比较值的;
205:根据步骤204中的初始最优适应度值交叉概率PJ和变异概率PB,依次对N个个体进行选择操作、交叉操作和变异操作,完成一次的进化,进化所得的适应度值记为过程最优适应度值 所对应的个体为过程最优个体J2,若大于则用J2替换J1后进行下一次进化,反之,直接进行下一次进化;
206:判断进化次数是否达到设定的最大进化次数,最大进化次数Gmax的取值范围为80~100次,若是,则进化停止,执行步骤207,若否,则返回步骤205;
207:将ρsj数值最大的个体作为最优化个体J,J所对应的双稳态系统为最优化双稳态系统XJ,提取XJ的势函数参数aJ、bJ和求解XJ时的数值计算步长hJ,aJ、bJ和hJ即为最优化参数。
208:步骤207得到aJ、bJ和hJ之后,利用变尺度方式下的双稳态系统对izg(t)进行特征提取,即将aJ替换双稳态系统表达式中势函数的参数a,将bJ替换双稳态系统表达式中势函数的参数b,将带有强噪声的暂态零序电流izg(t)替换双稳态系统表达式中的s(t)+Γ(t),采用四阶龙格-库塔算法的数值计算步长取值为hJ,求解得到变尺度方式下的双稳态系统的输出信号i′z(t),获取i′z(t)与iz(t)之间的互相关系数ρz′z,判断ρz′z是否大于设定阈值,设定阈值选取为0.92,若是,则执行步骤S3,若否,则跳转步骤201,通过遗传算法再次解算。此步骤208是为了对得到的最优化参数进行评价,评价标准为互相关系数的数值,若评价结果不能达到预定效果,则重新解算。
步骤S3:在最优化参数条件下,利用变尺度方式下的双稳态系统对噪声背景下的暂态零序电流进行特征提取,得到变尺度方式下的双稳态系统的输出信号ic(t),ic(t)即为微弱暂态零序电流故障特征提取的特征电流。
为了验证本发明方法可以获得更好的特征提取,分别对采用无优化的方法和本发明优化的方法进行暂态零序电流检测的试验:
1)无优化的暂态零序电流检测
定义理想暂态零序电流iz(t)为:
iz(t)=x1(t)+x2(t)+x3(t)+x4(t)+Γ(t)(11)
x1(t)=5.6cos(2π×50t+60°)(12)
x2(t)=40e-56tcos(2π×250t+30°)(13)
x3(t)=72e-102tcos(2π×315t)(14)
x4(t)=10e-5.5t(15)
式中:x1(t)为幅值较小的工频信号;x2(t)为幅值较大的5次谐波;x3(t)为幅值较大的非整次谐波;x4(t)为衰减直流分量;Γ(t)为噪声。值得指出的是,该暂态零序电流理想信号是按照小电流接地系统单相接地故障时的特点定义的。可见,信号iz(t)由5个信号叠加而成的,具有多频率和衰减特性,是一种非周期信号,如图6所示。在iz(t)中加入强噪声,得到带有强噪声的暂态零序电流izg(t),如图7所示,将izg(t)输入到式(2)中,其对应的参数分别为a=b=1,D=200db,采样频率fs=100000Hz,选择压缩频率CR为2000并进行数值仿真,如图8所示。
计算初始电流is(t)与无噪声背景下的暂态零序电流iz(t)之间的互相关系数,得互相关系数为0.7338。从图8可知,is(t)还是带有一定的噪声,且幅值相对于原来信号来说降幅较大。
2)经遗传算法优化的暂态零序电流检测
采用本发明所述的基于遗传算法优化的微弱暂态零序电流故障特征提取方法,暂态零序电流为式(11)所示,暂态零序电流的噪声强度为200,设置压缩频率CR为2000,优化对象a的搜索范围为[-10,10],优化对象b的搜索范围为[0,10],优化对象h的搜索范围为[0.0001,0.0219],种群规模N为55,交叉概率PJ为0.39,变异概率PB为0.18,最大进化次数Gmax为90,经程序运行,可得图9。由图9可知,经过11次迭代,互相关系数已经达到最大,也即算法收敛,此时的参数最优值为a=-7.8796,b=0.0036,h=0.0136,这也表明最大进化次数设置为90能够满足强噪声背景下的微弱暂态零序电流的提取。将最优化的参数代入变尺度下双稳态系统方程,进行四阶龙格-库塔方程求解,得出结果如图10所示。计算最优化参数下,变尺度下双稳态系统处理所得的特征电流与原始无噪声下的暂态零序电流之间的互相关系数,得互相关系数为0.9593。
对比图10与图8,从波形方面来分析,图10的波形比图8的波形更加清晰,也即图10波形所含的噪声比图8少;从幅值方面来分析,图10波形的幅值要远远大于图8;从互相关系数方面分析,在最优化参数条件下,互相关系数值为0.9593,在经验参数条件下,互相关系数值为0.7338,可见,经遗传算法优化后,互相关系数有了很大程度的提高。因此可得以下结论:1)双稳态系统具有提取强噪声背景下弱暂态零序电流的能力;2)合适的双稳态系统参数及数值计算步长可增强变尺度双稳态系统提取强噪声背景下弱暂态零序电流特征的能力。
分别设置噪声强度为50db、100db、200db、300db,对经验参数条件下和本发明所述特征提取方法进行对比分析,分别计算不同参数条件下的特征电流与无噪声条件下的暂态零序电流之间的互相关系数,其中,经验参数条件下的特征电流如图11所示,最优化参数条件下的特征电流如图12所示,互相关系数如表1所示:
表1不同参数条件下的互相关系数
对比图11和图12可知,图12各噪声强度下的特征电流所含的噪声均比图11少,图12中各特征电流的幅值均比图11大。此外,随着噪声强度的增加,图11和图12中的特征电流所包含的噪声也有所增加,说明变尺度双稳态系统提取强噪声背景下的弱暂态零序电流是有一定限度的。
从表1可知,随着噪声强度的增加,各条件下的互相关系数均有所下降;最优化参数条件下的互相关系数都维持在0.94以上,均属于强相关,互相关系数远远大于经验参数条件下的互相关系数,因此,经遗传算法优化后,变尺度下双稳态系统提取弱暂态零序电流故障特征的能力增强,因此,选取合理的势函数参数及数值计算步长是非常有必要的。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种多参数优化的弱暂态零序电流故障特征提取方法,其特征在于,包括:
步骤S1:设定双稳态系统中的势函数参数a、b和四阶龙格-库塔算法求解双稳态系统时的数值计算步长h为优化对象,其中,双稳态系统X的表达式为:
X=dx/dt=-dV(x)/dx+s(t)+Γ(t)
式中:dx/dt为双稳态系统的输出信号,V(x)为势函数,V(x)=-ax2/2+bx4/4;a和b为势函数参数;s(t)代表无噪声的暂态零序电流,Γ(t)代表噪声信号,t为时间,x为布朗粒子运动的速度;
步骤S2:采用遗传算法对双稳态系统X的表达式中势函数参数a、b和四阶龙格-库塔算法求解双稳态系统时的数值计算步长h进行优化,输出最优化参数aJ、bJ和hJ;
步骤S3:在最优化参数条件下,利用变尺度方式下的双稳态系统对噪声背景下的暂态零序电流进行特征提取,得到变尺度方式下的双稳态系统的输出信号ic(t),ic(t)即为微弱暂态零序电流故障特征提取的特征电流。
2.根据权利要求1所述的一种多参数优化的弱暂态零序电流故障特征提取方法,其特征在于,所述步骤S2包括以下步骤:
201:设定遗传算法的参数,所述遗传算法的参数包括种群规模N、交叉概率PJ、变异概率PB、优化对象a的搜索范围、优化对象b的搜索范围、优化对象h的搜索范围和最大进化次数Gmax,同时设定理想暂态零序电流iz(t),iz(t)中加入强噪声得到带有强噪声的暂态零序电流izg(t);
202:初始化种群,随机产生N个个体,其中,每个个体j对应一组优化对象组[aj,bj,hj],每组优化对象组对应一个双稳态系统Xj;
203:将izg(t)替换每个双稳态系统Xj的表达式中的s(t)+Γ(t),按个体编号顺序,依次对所有的Xj进行四阶龙格-库塔算法求解,求解双稳态系统时的计算步长记为hj,求解得到每个双稳态系统Xj的输出信号isj(t),isj(t)作为初始电流;
204:获取isj(t)与iz(t)之间的互相关系数ρsj,并将ρsj的值作为个体适应度;
205:根据个体适应度、交叉概率PJ和变异概率PB依次对N个个体进行选择操作、交叉操作和变异操作,完成一次的进化;
206:判断进化次数是否达到设定的最大进化次数,若是,则进化停止,执行步骤207,若否,则返回步骤205;
207:将ρsj数值最大的个体作为最优化个体J,J所对应的双稳态系统为最优化双稳态系统XJ,提取XJ的势函数参数aJ、bJ和求解XJ时的数值计算步长hJ,aJ、bJ和hJ即为最优化参数。
3.根据权利要求2所述的一种多参数优化的弱暂态零序电流故障特征提取方法,其特征在于,所述步骤S2还包括步骤207之后执行的步骤208,所述步骤208为:根据aJ、bJ和hJ,利用变尺度方式下的双稳态系统对izg(t)进行特征提取,得到变尺度方式下的双稳态系统的输出信号i′z(t),获取i′z(t)与iz(t)之间的互相关系数ρz′z,判断ρz′z是否大于设定阈值,若是,则执行步骤S3,若否,则跳转步骤201。
4.根据权利要求3所述的一种多参数优化的弱暂态零序电流故障特征提取方法,其特征在于,所述设定阈值为0.92。
5.根据权利要求2所述的一种多参数优化的弱暂态零序电流故障特征提取方法,其特征在于,所述优化对象a的搜索范围为[-10,10],所述优化对象b的搜索范围为[0,10],所述优化对象h的搜索范围为[0.0001,0.0219]。
6.根据权利要求2所述的一种多参数优化的弱暂态零序电流故障特征提取方法,其特征在于,所述最大进化次数Gmax的取值范围为80~100次。
7.根据权利要求2所述的一种多参数优化的弱暂态零序电流故障特征提取方法,其特征在于,所述步骤201中的强噪声的噪声强度取值为30db~500db。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510615562.1A CN105260777B (zh) | 2015-09-24 | 2015-09-24 | 一种多参数优化的弱暂态零序电流故障特征提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510615562.1A CN105260777B (zh) | 2015-09-24 | 2015-09-24 | 一种多参数优化的弱暂态零序电流故障特征提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105260777A true CN105260777A (zh) | 2016-01-20 |
CN105260777B CN105260777B (zh) | 2018-02-02 |
Family
ID=55100457
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510615562.1A Active CN105260777B (zh) | 2015-09-24 | 2015-09-24 | 一种多参数优化的弱暂态零序电流故障特征提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105260777B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105825197A (zh) * | 2016-03-29 | 2016-08-03 | 西安交通大学 | 一种线性化势阱壁的路径扩展随机共振微弱特征提取方法 |
CN105913042A (zh) * | 2016-04-28 | 2016-08-31 | 上海电力学院 | 一种基于优化级联双稳态的弱暂态零序电流特征提取方法 |
CN106019082A (zh) * | 2016-05-26 | 2016-10-12 | 上海电力学院 | 一种基于暂态零序电流的含dg配电网故障选线方法 |
CN106385272A (zh) * | 2016-09-14 | 2017-02-08 | 厦门大学 | 基于随机共振和时间反转镜的水下信号增强方法 |
CN106959684A (zh) * | 2017-03-29 | 2017-07-18 | 清华大学 | 基于压缩感知和联合范数优化的故障诊断方法和装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104007307A (zh) * | 2014-05-06 | 2014-08-27 | 西安理工大学 | 基于随机共振的小电流单相接地故障电流检测电路及方法 |
CN104462769A (zh) * | 2014-11-05 | 2015-03-25 | 中国南方电网有限责任公司电网技术研究中心 | 基于遗传算法的变电站电流量测回路故障检测方法 |
-
2015
- 2015-09-24 CN CN201510615562.1A patent/CN105260777B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104007307A (zh) * | 2014-05-06 | 2014-08-27 | 西安理工大学 | 基于随机共振的小电流单相接地故障电流检测电路及方法 |
CN104462769A (zh) * | 2014-11-05 | 2015-03-25 | 中国南方电网有限责任公司电网技术研究中心 | 基于遗传算法的变电站电流量测回路故障检测方法 |
Non-Patent Citations (4)
Title |
---|
JIE YANG: "Weak Signal Detection Based On Stochastic Resonance Combining With Genetic Algorithm", 《IEEE SINGAPORE INTERNATIONAL CONFERENCE ON COMMUNICATION SYSTEMS》 * |
任立通等: "基于随机共振的振动故障特征提取及可行性分析", 《空军工程大学学报(自然科学版)》 * |
吴利平等: "一种基于经典遗传算法的自适应随机共振系统", 《计算机科学》 * |
王晶等: "采用遗传算法的自适应随机共振系统弱信号检测方法研究", 《西安交通大学学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105825197A (zh) * | 2016-03-29 | 2016-08-03 | 西安交通大学 | 一种线性化势阱壁的路径扩展随机共振微弱特征提取方法 |
CN105825197B (zh) * | 2016-03-29 | 2019-01-08 | 西安交通大学 | 一种线性化势阱壁的路径扩展随机共振微弱特征提取方法 |
CN105913042A (zh) * | 2016-04-28 | 2016-08-31 | 上海电力学院 | 一种基于优化级联双稳态的弱暂态零序电流特征提取方法 |
CN105913042B (zh) * | 2016-04-28 | 2018-12-21 | 上海电力学院 | 一种基于优化级联双稳态的弱暂态零序电流特征提取方法 |
CN106019082A (zh) * | 2016-05-26 | 2016-10-12 | 上海电力学院 | 一种基于暂态零序电流的含dg配电网故障选线方法 |
CN106385272A (zh) * | 2016-09-14 | 2017-02-08 | 厦门大学 | 基于随机共振和时间反转镜的水下信号增强方法 |
CN106959684A (zh) * | 2017-03-29 | 2017-07-18 | 清华大学 | 基于压缩感知和联合范数优化的故障诊断方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN105260777B (zh) | 2018-02-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105260777A (zh) | 一种多参数优化的弱暂态零序电流故障特征提取方法 | |
CN105160405A (zh) | 基于遗传算法优化的微弱暂态零序电流故障特征提取方法 | |
CN109271975B (zh) | 一种基于大数据多特征提取协同分类的电能质量扰动识别方法 | |
CN103177188B (zh) | 一种基于特征映射的电力系统负荷动特性分类方法 | |
CN103245861B (zh) | 一种基于贝叶斯网络的变压器故障诊断方法 | |
CN102809718B (zh) | 一种gis超高频局部放电信号识别方法 | |
CN105300692B (zh) | 一种基于扩展卡尔曼滤波算法的轴承故障诊断及预测方法 | |
CN105259471A (zh) | 一种基于随机共振和暂态电流信号的三维故障选线方法 | |
CN106093724B (zh) | 一种基于混合粒子群的局部放电小波去噪方法 | |
CN103323755A (zh) | Gis特高频局部放电信号识别方法及系统 | |
CN102087337B (zh) | 非线性模拟电路诊断激励的退火遗传优化方法 | |
CN103728535A (zh) | 一种基于小波变换暂态能量谱的特高压直流输电线路故障测距方法 | |
CN105388390B (zh) | 基于粒子群优化的弱暂态零序电流故障特征提取方法 | |
CN104616061B (zh) | 基于小波包对数能量熵及遗传算法优化的孤岛检测方法 | |
CN102520342A (zh) | 基于动态反馈神经网络建模的模拟电路测试节点选择方法 | |
CN114692505A (zh) | 一种电容式电压互感器误差预测方法及装置 | |
CN103678886B (zh) | 一种基于地面测试数据的卫星贝叶斯网络健康确定方法 | |
CN112994017A (zh) | 基于配电网概率潮流计算的分布式光伏电源选址优化方法 | |
CN105447243A (zh) | 基于自适应分数阶随机共振系统的微弱信号检测方法 | |
CN106203382A (zh) | 一种基于核函数极限学习机的励磁涌流和故障电流识别方法 | |
CN104268629A (zh) | 一种基于先验信息和网络固有信息的复杂网络社区检测方法 | |
CN104503420A (zh) | 一种基于新型fde-elm和时延efsm的非线性过程工业故障预测方法 | |
CN108108672A (zh) | 一种基于线性搜索策略的随机共振电流弱信息识别方法 | |
Alshareef et al. | Application of wavelet-based ensemble tree classifier for non-intrusive load monitoring | |
CN107121617B (zh) | 一种采用滤波器支路电流的直流输电线路k近邻测距方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |