基于TFSR-CEEMDAN的齿轮箱早期故障预警方法.pdf

收稿日期:2 0 2 3-0 6-2 7作者简介:马 健(1 9 9 7),男,硕士生,主要研究方向为风机齿轮箱故障诊断,E-m a i l:1 0 6 1 4 7 4 2 7 1 q q.c o m 指导教师:徐 旭(1 9 7 9),男,副教授,博士,主要研究方向为物流系统设计与优化、供应链韧性与风险管理、精益生产系统可视化与数字孪生,E-m a i l:x u x s d j u.e d u.c n文章编号 2 0 9 5-0 0 2 0(2 0 2 3)0 5-0 2 6 2-0 7基于T F S R-C E EMD AN的齿轮箱早期故障预警方法马 健,徐 旭,何晨光,苏晨雅(上海电机学院 商学院,上海 2 0 1 3 0 6)摘 要 针对强噪声环境下,齿轮箱早期故障预警困难的问题,提出一种基于时滞反馈随机共振(T F S R)自适应噪声完备集合经验模态分解(C E EMD A N)的齿轮箱早期故障预警方法。首先通过T F S R的长记忆系统,根据历史信息增强反馈过程中的信号周期性,产生低通滤波下增强的故障特征信号;再结合C E EMD A N将增强后的故障特征信号自适应分解;并以信噪比为指标筛选I M F分量,得到的重构信号可以有效增强故障特征频率及其谐频,达到早期故障预警的目的。实验结果表明:通过信号发生器仿真分析,并以真实风机齿轮箱的数据,对比正常和磨损状态下的齿轮,检测出强噪声背景下齿轮箱的微弱局部磨损故障,有效提高故障特征频率处的信噪比,验证了齿轮箱早期故障预警的灵敏性。关键词 齿轮箱;时滞反馈随机共振;完备集合经验模态分解;信噪比;早期故障预警中图分类号 TH 1 6 5+.3文献标志码 AR e s e a r c h o n e a r l y f a i l u r e w a r n i n g m e t h o d o f g e a r b o x b a s e d o n T F S R-C E E M D A NMA J i a n,X U X u,HE C h e n g u a n g,S U C h e n g y a(S c h o o l o f B u s i n e s s,S h a n g h a i D i a n j i U n i v e r s i t y,S h a n g h a i 2 0 1 3 0 6,C h i n a)A b s t r a c t T o s o l v e t h e p r o b l e m o f d i f f i c u l t y i n e a r l y f a i l u r e w a r n i n g o f g e a r b o x e s i n a s t r o n g n o i s e e n v i r o n m e n t,a n e a r l y f a u l t w a r n i n g m e t h o d f o r g e a r b o x e s b a s e d o n t i m e-d e l a y f e e d b a c k s t o c h a s t i c r e s o n a n c e(T F S R)-a d a p t i v e n o i s e c o m p l e t e e n s e m b l e e m p i r i c a l m o d e d e c o m p o s i t i o n(C E EMD A N)i s p r o p o s e d.F i r s t,t h e s i g n a l p e r i o d i c i t y i n t h e f e e d b a c k p r o c e s s i s e n h a n c e d t h r o u g h t h e l o n g m e m o r y s y s t e m o f T F S R b a s e d o n h i s t o r i c a l i n f o r m a t i o n,a n d t h e e n h a n c e d f a u l t c h a r a c t e r i s t i c s i g n a l u n d e r l o w-p a s s f i l t e r i n g i s g e n e r a t e d.T h e n c o m b i n e d w i t h C E EMD A N,t h e e n h a n c e d f a u l t c h a r a c t e r i s t i c s i g n a l i s a d a p t i v e l y d e c o m p o s e d.T h e I M F c o m p o n e n t i s s c r e e n e d w i t h t h e s i g n a l-t o-n o i s e r a t i o a s t h e i n d e x,a n d t h e r e c o n s t r u c t e d s i g n a l o b t a i n e d c a n e f f e c t i v e l y e n h a n c e t h e c h a r a c t e r i s t i c f r e q u e n c y o f t h e f a u l t a n d i t s h a r m o n i c f r e q u e n c y t o a c h i e v e t h e p u r p o s e o f e a r l y f a u l t w a r n i n g.T h e e x p e r i m e n t a l r e s u l t s s h o w t h a t t h r o u g h s i m u l a t i o n a n a l y s i s o f t h e s i g n a l 第2 6卷 第5期2 0 2 3年上 海 电 机 学 院 学 报J O U R N A L O F S HA N G HA I D I A N J I U N I V E R S I T YV o l.2 6 N o.5 2 0 2 3 g e n e r a t o r a n d t h e d a t a f r o m r e a l w i n d t u r b i n e g e a r b o x e s,w e a k l o c a l w e a r f a u l t s o f t h e g e a r b o x u n d e r a s t r o n g n o i s e b a c k g r o u n d c a n b e d e t e c t e d b y c o m p a r i n g t h e g e a r s u n d e r n o r m a l a n d w e a r c o n d i t i o n s,t h e s i g n a l-t o-n o i s e r a t i o a t f a u l t c h a r a c t e r i s t i c f r e q u e n c y f a u l t i s e f f e c t i v e l y i m p r o v e d a n d t h e s e n s i t i v i t y o f t h e e a r l y f a i l u r e w a r n i n g o f t h e g e a r b o x i s v e r i f i e d.K e y w o r d s g e a r b o x;t i m e-d e l a y f e e d b a c k s t o c h a s t i c r e s o n a n c e(T F S R);c o m p l e t e e n s e m b l e e m p i r i c a l m o d e d e c o m p o s i t i o n(C E EMD A N);s i g n a l-t o-n o i s e r a t i o;e a r l y f a i l u r e w a r n i n g 由于齿轮传动比大、承载能力强等优点1,齿轮箱被广泛应用在汽车、飞机和风电机组中,但是齿轮箱的关键部件齿轮和轴承,容易受到疲劳裂纹和点蚀等损坏的影响而产生故障2。齿轮箱的故障信号具有非平稳性、频域宽和强度弱等特点,故障信息获取困难和采样环境噪声过大等问题都增加了故障早期预警的难度。L u等3提出一种基于时滞反馈随机共振(T i m e-D e l a y e d F e e d b a c k S t o c h a s t i c R e s o n a n c e,T F S R)模型的非平稳弱信号检测策略,考虑到故障振动成分的周期性和微弱性,将历史信息适当地叠加到随机共振(S t o c h a s t i c R e s o n a n c e,S R)输出中,获得低噪声干扰的规则输出波形,验证了该算法适用于检测非平稳特性的信号或受到严重多尺度噪声干扰的信号。H u等4提出了一种基于多时滞反馈随机共振模型(M u l t i-t i m e-d e l a y e d F e e d b a c k S t o c h a s t i c R e s o n a n c e,MT F S R)的微弱信号检测方法,此方法可利用多个延迟反馈项叠加形成的历史信息来增强信号的周期性,适用于齿轮箱故障诊断。韩雪飞等5利用多点最优最小熵解卷积对原始信号进行降噪,并采用完全自适应噪声集合经验模态分解(C o m p l e t e E EMD w i t h A d a p t i v e N o i s e,C E EMD A N)处理信号,证明在滚动轴承故障诊断方面的有效性。蒋玲莉等6提出一种基于C E EMD A N排列熵与支持向量机(S u p p o r t V e c t o r M a c h i n e,S VM)的齿轮故障辨别方法,相较于集合经验模态分解排列熵具有更高的齿轮故障诊断准确度。本文基于T F S R和C E EMD A N算法,提出一种基于T F S R-C E EMD A N的齿轮箱早期故障预警的方法。首先,通过选择适当的时延、反馈强度和计算步长,实现高频噪声能量向低频区域的转移以获取最优的微弱故障特征增强效果,并结 合C E EMD A N分解降噪信号,以信噪比为指标来重构I M F分量。其次,通过计算故障频率及其倍频处的信噪比,放大微弱的故障特征变化,实现故障预警诊断。最后,通过仿真和实验分析验证了该方法在齿轮箱早期故障预警诊断的灵敏性。1 算法原理1.1 时滞反馈随机共振齿轮箱的故障经常以周期性冲击信号存在,并且故障频率及其谐频成分蕴含着故障信息,考虑到早期微弱故障的故障特征频率对应幅值非常微弱,因此有必要进行微弱故障信号增强7。非线性系统输出表达式为X0=a x-b x3+s(t)+n(t)(1)式中:a、b为系统参数;x为影响系统的物理学变量;s(t)为待检测的故障信息信号;n(t)为噪声。待检测的故障信息信号为s(t)=As i n(2 fot)(2)式中:A为信号幅值;fo为信号频率;t为时间。噪声为n(t)=2D g(t)(3)式中:D为噪声强度;g(t)为自噪声8。经典S R的数学模型9在非线性系统中的表达,原理如图1所示。图中U(x)=-a2x2+b4x4为系统模型,以=a/b为系统参数,表述质点在势阱中的过阻尼布朗运动。通过在势阱中的切换,使信号、噪声和系统产生同步,令输出得到增强。图1 S R模型时滞反馈随机共振在经典S R的基础上,通过362 2 0 2 3年第5期马 健,等:基于T F S R-C E EMD A N的齿轮箱早期故障预警方法 加入一个反馈项,提高增强弱信号的能力,原理如图2所示。T F S R系统可用下式表示:ddX0(t)dt=-U(x)x+x(t-)+s(t)+n(t)(4)式中:为反馈强度;为滞后时间。图2 T F S R模型这些结构参数影响粒子在势场中的振荡效率,进而影响从背景噪声中提取故障信号的性能。可保留输出信号的信息,并以调节来增强输出。经典S R的响应本质上是系统方程的解。在积分过程中,不同输入的剧烈抖动将被平滑。而输出是从输入累积的,信号的高频部分在积分处理之后被擦除1 0,因此,T F S R可以看作是一种特殊的低通滤波器,而引入的时延反馈项对S R的输出以及信号的滤波性能都有影响。1.2 C E E MD A N算法T o r r e s1 1等在集合经验模态分解(E n s e m b l e E m p i r i a l M o d e D e c o m p o s i t i o n,E E EMD)基础上提出了C E E M D A N,本质是向E E M D过程的每个阶段逐一加入自适应的白噪声,并且以较低的计算代价实现了本征模式函数更好的谱分离。C E E M D A N因其高分解效率和高重建准确性被广泛应用1 2,但如果直接采用C E E M D A N分解信号,I M F分量中的高频噪声会导致故障特征提取效果欠佳或根本难以提取。因此,先将信号通过T F S R系统,对含强噪声的信号进行消噪处理,以突出其原有的故障冲击成分,再对此时的信号进行C E E M D A N自适应分解,能保持信号的完整,也可解决模态混叠的现象。2 算法实现原始信号经过T F S R增强后信号的幅值已经衰减,需要一个描述故障特征频率处的能量变化的指标,信噪比则适用于此时的输出。信噪比1 3是信号功率谱中特征频率处的幅值与噪声之比,为使得系统输出信号具有最佳信噪比,需要调节T F S R的反馈强度和滞后时间,信噪比S为S=f(,),(-,+),(0,T)(5)式中:T为输入信号周期。假设反馈强度(-2 0 0 0 0,2 0 0 0 0),滞后时间(0,0.0 0 0 3),S、与的关系如图3所示。图3 信噪比、反馈强度和滞后时间的函数图像由于恶劣的运行环境下齿轮箱早期故障信息微弱,而啮合频率处的信噪比在强噪声的影响下难以展现真实的故障特征信息变化。因此,以信噪比为目标函数,遍历选择最佳反馈强度、滞后时间,使T F S R系统中的微弱信号、噪声和电势得到最大程度的匹配,以得到最佳的故障信息增强效果。再对滤波信号进行C E EMD A N自适应分解,重构大于平均信噪比的分量。最后计算重构后的信号啮合频率处及其倍频的信噪比,与正常运行状态下齿轮的指标分析对比,达到早期故障预警的效果。方法的流程如图4所示。图4 早期故障预警流程462上 海 电 机 学 院 学 报 2 0 2 3年第5期3 仿真验证为进一步验证基于T F S R-C E EMD A N的早期故障预警指标诊断方法在真实环境下的有效性,本研究使用型号D G 1 0 2 2 Z的信号发生器与工控机模拟平台模拟齿轮箱早期磨损故障1 4,通道分别产生啮合频率2 0 0H z和边频2 0H z的齿轮磨损故障信号,故障特征频率为8 7H z的轴承外圈故障信号,方差为5的白噪声信号,并通过A R T动态采集卡采集,采样频率1 2k H z,采样点数为8 1 9 2,得到的信号如图5所示。图5 采集信号由图5可见,信号的时域被噪声干扰严重,并且由于轴承故障冲击信号的影响,啮合频率周围的干扰较多。通过对工控机采集信号进行T F S R-C E EMD A N算法处理,得到的信号如图6所示。由图6可见,由于轴承冲击信号和滤波的影响,时域信号同样具有部分失真现象,但是降噪的效果明显,并且在频谱图上能观察到在基频2 0 0H z附近,有1 8 0.1 7 6H z和2 1 9.7 2 7H z的谱峰。这和设置的2 0H z的调制频率相吻合。但是由于干扰信号和滤波的影响,边频带的幅值并不对称,并且在边频带产生轴承二倍频1 7 4.3 1 6H z的较高谱峰,这是由于轴承的二倍频和待增强的微弱故障信号处于同一频带。这时只能初步判 图6 重构信号断齿轮发生早期磨损故障,再计算此时基频处的信噪比为-9.6 8 7 7,对比初始值-1 9.1 0 19相比,同样得到了提高。在现场监测环境中,结合频谱图和信噪比的变化便可以对齿轮的早期故障进行预警诊断。在工程现场监测环境下,要求预警诊断不仅具有有效性,二期需要具备可行性,算法的运行时间就必须在合理的范围内。因此,在这里固定采样频率为1 2k H z,实验不同采样点数下的执行时间,得到的结果如表1所示,根据执行结果,建议在采用本算法时采样点数设置在8 1 9 2以内。表1 不同采样点数下的执行时间采样点数1 0 2 42 0 4 84 0 9 68 1 9 21 6 3 8 4执行时间/s2.2 6 8 13.0 5 6 95.1 3 8 29.7 9 1 42 1.2 5 2 44 实验验证本文实验数据来源于蒙东某风场,风机齿轮箱类型为二级行星二级平行齿轮箱(即主轴连接一级行星轮系,一级行星轮系连接二级行星轮系,行星轮系连接中速轴,中间轴连接高速轴),传感器安装在二级内齿圈上径向采集振动加速度信号,传感器类型为中高传感器,数据单位均为g,测得一组齿轮早期磨损数据和一组正常运行条件下的健康数据。其中正常运行时数据采样频率为5.1 2k H z,采样点数为3 2 7 6 9;早期磨损数562 2 0 2 3年第5期马 健,等:基于T F S R-C E EMD A N的齿轮箱早期故障预警方法 据采样采样频率为5.1 2k H z,采样点数为1 6 3 9 4。齿轮的特征频率,即齿轮的旋转频率和齿轮间的啮合频率1 5为f2=f1i1=N1N2f1(6)f3=f2=N1N2f1(7)fm1=N1f1(8)式中:N1、N2为齿轮的齿数;f1、f2、f3为齿轮的旋转频率;i1为啮合副的传动比;fm1为齿轮的啮合频率。计算出磨损故障下的啮合频率为1 5 7.2 3H z,正常运行下的啮合频率为1 5 3.9 1H z。其中磨损故障的时域和频域如图7所示,可以看出,因为机械传动系统受到风扰动的影响,风电机组的转动速度不断变化,采集到的齿轮箱信号具有非平稳性和复杂性的特点,需要对其进行滤波处理。图7 磨损故障信号经过T F S R的处理,结果如图8所示,可以看出干扰噪声得到抑制和削弱,信号呈现周期性,并且高速轴的转频及其二倍频和二级行星齿轮的啮合频率较为突出,低频带的信息得到极大的增强,中高频带的信息能量大幅衰减,但依旧存在时域信号周期性不够明显、信号中含有少量冲击等问题,接下来使用C E EMD A N进一步的筛选处理。滤波增强后的信号经过C E EMD A N处理,平均实现次数为5 0 0,筛选迭代次数为50 0 0,得到I M F分量。图8 T F S R增强后的信号滤波增强被分解成1 7个I M F分量,其中第1 7个分量为余项,从I M F 1到I M F 1 6波形的混叠现象逐渐减弱。将所有I M F分量啮合频率处的平均信噪比计算出来,筛选出第1 8一共8个I M F分量作为重构信号的基础,进行叠加,得到的信号如图9所示。图9 重构信号由图9可见,发生磨损故障时,信号的时域波形呈现周期性伴随明显的调幅调频特性,并且经过此算法的处理,已基本消除噪声的干扰,可以呈现有效信息。通过频域图可以看出,齿轮的转频及其倍频较为明显,这是因为经过T F S R后,低频带的突出频率都得到了增强。如图1 0所示,为正常运行条件下齿轮的时域和频域图,将正常运行下的频域和发生磨损时的频域对比可以看出,发生磨损故障时,啮合频率的右侧边频带明显突出,此时可以初步判断齿轮发生了662上 海 电 机 学 院 学 报 2 0 2 3年第5期 图1 0 正常运行信号磨损故障。对比正常条件下和磨损时的信噪比指标分析,如表2所示。表2 行星齿轮两种状态信噪比指标对比状态信噪比(啮合频率处)原始信号一倍频算法处理信号一倍频算法处理信号倍频之和正常运行-3 7.8 8 8 5-1 9.0 7 4 0-9 5.6 0 1 6早期磨损-4 7.1 8 9 5-2 7.7 6 6 2-1 3 9.0 1 9 9 由表2可知,经过T F S R-C E EMD A N算法的处理,啮合频率处的信噪比得到很大的提高,能突出有效故障信息。通过正常运行和早期磨损啮合频率处一倍频的信噪比对比,可以发现,发生磨损时,信噪比有所降低,这是因为发生故障时啮合频率周围的边频带开始凸显,涌现调幅调频信号,这使得信噪比开始降低。但是由于时早期磨损故障,在信噪比指标上对比存在差值过小,无法准确判断是否存在故障。在齿轮发生故障时,不仅是啮合频率周围会涌现边频带,在啮合频率的倍频周围同样会出现故障信息。因此引入啮合频率和其二倍频率、三倍频率之和作为健康指标来放大故障信息的变化。通过正常运行和早期磨损的数据对比可知,此时已经发生早期故障磨损。5 结 论本文提出了一种基于T F S R-C E EMD A N的齿轮箱早期故障预警指标诊断方法,通过仿真和实验证明能够有效地检测出强噪声背景下齿轮箱的微弱局部磨损故障,并得出以下结论:(1)针对齿轮箱干扰噪声过大,采用经典的随机共振系统响应信噪比较低的问题,该方法充分利用了时滞反馈随机共振的增强故障信息的作用,而C E EMD A N利用了T F S R的非线性低通滤波特性,可以逐步滤除干扰,实现高频噪声能量向低频区域的进一步转移,从而削弱系统响应中的高频抖动,筛选出信噪比高的重构信号。(2)针对早期故障不明显的问题,通过对啮合频率及其倍频的信噪比考虑,放大微弱故障信息的变化,监测系统的输出信噪比来判断齿轮箱状态,是一种新的齿轮箱早期故障预警的方法。本文以故障频率处的信噪比作为指标,可有效地突出故障特征频率处的能量变化。因此,本研究后续可编写一个相关的L A B V I EW软件包,应用于实际的工程诊断,具有一定的工程价值。参 考 文 献1 L E I Y,L I N J,Z U O M J,e t a l.C o n d i t i o n m o n i t o r i n g a n d f a u l t d i a g n o s i s o f p l a n e t a r y g e a r b o x e s:a r e v i e wJ.M e a s u r e m e n t,2 0 1 4,4 8(1):2 9 2-3 0 5.2 C HA A R F,F A KH F A KH T,HA D D A R M.D y n a m-i c a n a l y s i s o f a p l a n e t a r y g e a r f a i l u r e c a u s e d b y t o o t h p i t t i n g a n d c r a c k i n gJ.J o u r n a l o f F a i l u r e A n a l y s i s&P r e v e n t i o n,2 0 0 6,6(2):7 3-7 8.3 L U S L,H E Q B,Z HA N G H B,e t a l.E n h a n c e d r o-t a t i n g m a c h i n e f a u l t d i a g n o s i s b a s e d o n t i m e-d e l a y e d f e e d b a c k s t o c h a s t i c r e s o n a n c eJ.J o u r n a l o f V i b r a-t i o n a n d A c o u s t i c s:T r a n s a c t i o n s o f t h e A S M E,2 0 1 5,1 3 7(5):5 1 0 0 8.4 HU B B,L I B.F a u l t d i a g n o s i s o f g e a r b o x u s i n g m u l t i-t i m e-d e l a y e d f e e d b a c k s t o c h a s t i c r e s o n a n c eJ.A R C H I V E P r o c e e d i n g s o f t h e I n s t i t u t i o n o f M e c h a n i c a l E n g i n e e r s P a r t C J o u r n a l o f M e c h a n i c a l E n g i n e e r i n g S c i e n c e,2 0 1 7,2 3 1(1 9):3 5 4 1-3 5 5 2.5 韩雪飞,施展,华云松.基于参数优化MOM E D A与C E E MD A N的滚动轴承微弱故障特征提取研究J.机械强度,2 0 2 1,4 3(5):1 0 4 1-1 0 4 9.6 蒋玲莉,谭鸿创,李学军,等.基于C E E MD A N排列熵与S VM的螺旋锥齿轮故障识别J.振动.测试与诊断,2 0 2 1,4 1(1):3 3-4 0,1 9 8-1 9 9.762 2 0 2 3年第5期马 健,等:基于T F S R-C E EMD A N的齿轮箱早期故障预警方法 7 谯自健,陈帅,马莉,等.多稳态匹配随机共振在机械早期故障特征提取中的应用J.振动与冲击,2 0 2 3,4 2(1 1):8 7-9 5.8 Y A N G D X,HU N Q.N u m e r i c a l s i m u l a t i o n o f s t o-c h a s t i c r e s o n a n c e i n b i s t a b l e s y s t e m f o r d e t e c t i n g w e a k s i g n a lJ.J o u r n a l o f N a t i o n a l U n i v e r s i t y o f D e f e n s e T e c h n o l o g y,2 0 0 3(6):1 0 0 1-2 4 8 6.9 曾小玲,任芮彬,邓科.含相关色噪声和周期方波信号的双稳系统的随机共振J.四川大学学报(自然科学版),2 0 2 3,6 0(4):2 3-3 0.1 0 Y A N G D X,HU N Q.D e t e c t i n g w e a k s i g n a l b a s e d o n a n a l o g s i m u l a t i n g o f s t o c h a s t i c r e s o n a n c eJ.J o u r n a l o f C i r c u i t s a n d S y s t e m s,2 0 0 4,9(6):1 3 5-1 3 8.1 1 T O R R E S M E,C O L OM I N A S M A,S C H L O T-T HA U E R G,e t a l.A c o m p l e t e e n s e m b l e e m p i r i c a l m o d e d e c o m p o s i t i o n w i t h a d a p t i v e n o i s eC I E E E I n t e r n a t i o n a l C o n f e r e n c e o n A c o u s t i c s,S p e e c h a n d S i g n a l P r o c e s s i n g.P r a g u e:I E E E,2 0 1 1:4 1 4 4-4 1 4 7.1 2 安冬,魏 亚 静,张 倩,等.轴 系 轮 廓 重 构 信 号 的C E E MD A N小波阈值去噪方法J.工具技术,2 0 2 3,5 7(8):1 5 3-1 5 9.1 3 钟建军,宋健,由长喜,等.基于信噪比评价的阈值优选小波去噪法J.清华大学学报(自然科学版),2 0 1 4,5 4(2):2 5 9-2 6 3.1 4 李肖,冯志鹏.行星齿轮箱齿轮组合故障振动频谱特征J.振动与冲击,2 0 2 0,3 9(1):1 5-2 3.1 5 李成冬,吕红明,杨旭,等.变速器齿轮非线性振动机理研究及稳定性分析J.机械设计与制造,2 0 2 3(1):2 3 7-2 4 1.(上接第2 6 1页)参 考 文 献1 何瑞文,龙隆,张宝仁,等.电力信息物理系统中信息系统物理化的建模及分析方法J/O L.(2 0 2 2-0 9-0 9)2 0 2 3-0 8-0 4.h t t p:k n s.c n k i.n e t/k c m s/d e t a i l/1 1.2 1 0 7.TM.2 0 2 2 0 9 0 9.1 5 0 0.0 0 2.h t m l.2 王琦,李梦雅,汤奕,等.电力信息物理系统网络攻击与防御研究综述(一)建模与评估J.电力系统自动化,2 0 1 9,4 3(9):9-2 1.3 娄文静,王蕾,焦提操,等.基于电网拓扑优化的电压稳定预防控制方法J.南方电网技术,2 0 2 2,1 6(1):7 5-8 2.4 肖飞,杨国健,邓祥力,等.基于电网故障拓扑分析及多数据综合的复杂故障诊断方法J.水电能源科学,2 0 2 0,3 8(2):1 8 9-1 9 2.5 刘迪,张强,吕干云,等.基于支路有功功率的配电网拓扑辨识方法J.电力工程技术,2 0 2 1,4 0(3):9 2-9 8.6 A S HO K A,G O V I N R A S U M.C y b e r a t t a c k s o n p o w e r s y s t e m s t a t e e s t i m a t i o n t h r o u g h t o p o l o g y e r-r o r sC2 0 1 2 I E E E P o w e r a n d E n e r g y S o c i e t y G e n e r a l M e e t i n g.S a n D i e g o,C A:I E E E,2 0 1 2:1-8.7 K I M J,T O N G L.O n t o p o l o g y a t t a c k o f a s m a r t g r i d:u n d e t e c t a b l e a t t a c k s a n d c o u n t e r m e a s u r e sJ.I E E E J o u r n a l o n S e l e c t e d A r e a s i n C o mm u n i c a t i o n s,2 0 1 3,3 1(7):1 2 9 4-1 3 0 5.8 D E N G R L,Z HU A N G P,L I A N G H.C C P A:c o o r d i-n a t e d c y b e r-p h y s i c a l a t t a c k s a n d c o u n t e r m e a s u r e s i n s m a r t g r i dJ.I E E E T r a n s a c t i o n s o n S m a r t G r i d,2 0 1 7,8(5):2 4 2 0-2 4 3 0.9 L I A N G G Q,WE L L E R S R,Z HA O J H,e t a l.A f r a m e w o r k f o r c y b e r-t o p o l o g y a t t a c k s:l i n e-s w i t c h i n g a n d n e w a t t a c k s c e n a r i o sJ.I E E E T r a n s a c t i o n s o n S m a r t G r i d,2 0 1 9,1 0(2):1 7 0 4-1 7 1 2.1 0 王胜锋,丁洲,吴劲松,等.基于拓扑篡改的电力市场虚假数据注入攻击方案J.电力自动化设备,2 0 2 1,4 1(1 1):1 4 7-1 5 2.1 1 何宗伦,高仕斌,韦晓广,等.支路与保护协同篡改的虚假拓扑攻击攻防博弈模型研究J.电网技术,2 0 2 2,4 6(1 1):4 3 4 6-4 3 5 5.1 2 蔡晔,刘放,曹一家,等.电力信息物理系统低代价多阶段高危攻击策略研究J.电力系统自动化,2 0 2 1,4 5(2 0):1-8.1 3 凌霞.多支路故障潮流转移危险线路识别及过载消除控制研究D.北京:华北电力大学,2 0 1 8.1 4 张婧.计及拓扑结构与运行状态的电网关键环节辨识D.长沙:湖南大学,2 0 2 0.1 5 阳育德,蓝水岚,覃智君,等.电力信息物理融合系统的网络 物理协同攻击J.电力自动化设备,2 0 2 0,4 0(2):9 7-1 0 3.1 6 余杰.基于帝国竞争算法的冷热电联供优化调度研究D.南昌:南昌大学,2 0 2 1.862上 海 电 机 学 院 学 报 2 0 2 3年第5期