2. 重庆大学 输配电装备及系统安全与新技术国家重点实验室, 重庆 400044
2. State Key Laboratory of Power Transmission Equipment & System Security and New Technology, Chongqing University, Chongqing 400044, P. R. China
舞动是指导线在风力及(或)不对称覆冰作用下引起的低频率、大幅度的自激振动,是一种空气动力失稳现象[1, 2]。舞动多发生在冬季,其能量很大,且持续时间长,易对输电线路造成机械损坏和电气故障,轻则引起相间闪络、损伤导线、地线及金具等,重则导致断股、断线、甚至倒塔等恶劣事故,严重威胁着输电线路的安全稳定运行[3, 4]。运行观测和统计资料表明,中国是舞动灾害最严重的国家之一[5, 6]。随着中国电网规模的发展和恶劣天气的频繁出现,输电线舞动事故的发生频率和危害程度均有明显的增加,且舞动区域已不仅局限在少数范围内,也遍及到国家电网的大部分地区[7]。因此,对输电线路舞动及其预防措施的研究具有重要的理论意义和工程实用价值。
近年来,国内外学者对导线舞动激发机理、计算仿真及导线防舞动措施等进行了多方面的研究,取得了很多重要成果并应用于实际工程。遗憾的是,由于输电线和气流相互作用所造成的耦合以及输电线的大幅运动所造成的几何非线性等,使舞动问题变得十分复杂,迄今为止尚无统一、普适的输电线舞动激发理论。目前具有代表性的主要有4种[8],分别为Hartog[9]提出的横向起舞机理、Nigol等[10]提出的扭转振动起舞机理、Yu等[11]提出的偏心惯性耦合起舞机理及中国学者尤传永[12]提出的动力失稳起舞机理。
现行的防舞措施概括起来可分为3大类[13]:一是从气象条件考虑,避开易于形成舞动的覆冰区域与线路走向;二是从机械与电气的角度,设计提高线路系统的抗舞动能力;三是采取各种防舞装置,抑制舞动的发生。但值得注意的是:其一,考虑到节省线路走廊的经济性设计要求以及施工便宜性等因素,使得部分输电线路是无法完全避开舞动区的;其二,在实际应用中,加强输电线路质量的技术及相关的抗舞设计尚不够详实和规范,经济性以及可操作性较差,同时也缺乏实践经验;其三,防舞器是基于不同舞动机理开发得到的,导致目前应用较多的几种防舞装置都具有其明显的设计特点和应用局限性,防舞效果也存在很大差异[7]。显然,要想完全制服舞动还需很长一段时间的努力,而目前迫切需要一种主动性更强、应用范围更广的输电线防舞动辅助决策方法,来减轻输电系统遭受的舞动灾害。
一方面,气象条件是影响舞动激发的最重要因素,只要满足舞动条件,各区域、各电压等级的输电线均可能发生舞动[14];另一方面,近年来输电线舞动在线监测装置和方法得到了蓬勃发展[15, 16],而舞动在线监测系统的成功运行为研究导线舞动积累了宝贵的气象资料。因此,为了弥补现有防舞措施的不足,加强输电线路的防舞能力,笔者在分析影响舞动的环境气象因素基础上,提出一种基于Adaboost学习算法的输电线舞动气象预警模型。预警结果可为电网运维人员提供决策支撑,提前做好有针对性的抗风、除融冰等防舞措施,避免防范不足引起事故及过度防范浪费资源,保障输电线路的安全稳定运行。
1 线路舞动的影响因素分析导线舞动是一种复杂的流固耦合振动,其形成因素很多,经过多年来国内外学者对舞动的研究,归纳得出了引起舞动的3个重要因素[2, 5, 7]。
1)导线覆冰:引起导线舞动的决定性因素。导线上要形成覆冰,必须具备3个条件:一是空气湿度较大,一般在85%以上;二是合适的温度,一般为0~-5 ℃;三是可使空气中水滴运动的合适风速,一般为1~15 m/s。
2)风的激励:导线舞动的直接原因。一段线路的舞动状态及其强弱,除了与风速大小有关外,还取决于风向对导线轴线(线路走向)的夹角。
3)线路结构与参数:引起舞动的内因。不合理的线路结构参数组合易引起舞动,这些因素主要包括:导线的类型、导线截面积、线路档距等。
综上可见,导线舞动的形成与外界气象环境密不可分,这些气象条件包括:风速、风向对导线轴线夹角、温度和相对湿度。引起导线舞动的内因是线路结构与参数,而当内因相对不变时,事物的变化由外因决定。因此,笔者结合文献[17]给出的舞动导线参数的相关资料,对引起导线舞动的内因(导线结构、导线截面和档距)按表 1所示进行划分,然后据此将线路历史舞动资料进行归类,得到共计18类线路(如:单导线、小截面、小档距)的舞动情况统计数据,在此基础上建立输电线舞动的气象预警模型。
![]() |
表1 导线舞动内因分类统计表 Table 1 Classification of galloping internal cause |
由于高度的变化对风速有影响,而气象站观测得到的风速以及气象预报风速一般是默认为离地10 m高的风速,因此,应按式(1)统一折算至导线高度处的风速vl[18]。
vl=vq⋅(H10)μ | (1) |
由于输电线舞动的物理模型不够精确,且模型中的部分参量在实际线路上难以通过测量获取,导致利用物理模型进行舞动预警的实用性和准确性较低。但是,由于对总体起作用的环境和条件的变化必然影响着每一个个体,因此,受动力学规律支配的单个客体的行为只能在一定范围内偏离总的方向,且它们在总体上仍表现出统计的必然性[19]。因此,统计学习理论是面对数据而又缺乏理论模型时最基本的分析手段。
机器学习基于以往的观察来获取较准确的预测,它提供了一种从观测数据出发得到目前尚不能通过原理分析得到的规律,进而利用这些规律预测未来数据的方法,其核心任务是把样例分类到各可能的离散值对应的类别中,并对未知样本做出尽可能精确的预测[20]。
综上所述,输电线舞动预警问题可归结为有监督学习下的分类预测问题:对照实际输电线路情况和表 1,选择18种线路类型中与之对应的一类,以该类线路的历史舞动情况统计数据(风速、风向对导线轴线夹角、温度和相对湿度)构成训练样本集,以下述的Adaboost集成学习算法形成分类学习器,并根据气象预报数据进行导线舞动预警。
2.2 Adaboost集成学习算法对一个预测分类问题,统计学习的基本目标是根据观测数据建立具有较强泛化能力的学习器。但在很多情况下,由于学习器的精确性受领域知识和训练数据及其分布(尤其是对那些还未完全了解其本质的预测问题,如导线覆冰舞动)的影响很大,使得很难直接构造具有高精度的学习器。然而,产生数个只比随机猜测略好的弱学习器却很容易,因此,寻找一般的提高已有学习器精确性的方法是很有价值的,这在直接构造强学习器非常困难的情况下,为学习算法的设计提供了一种有效的新思路和新方法,集成学习正是在这一背景下应运而生的[21]。
自集成学习的概念被提出以来,其在很多领域得到了快速的发展和广泛的应用。目前常用的集成学习算法有Bagging算法(bootstrap aggregation自助聚集法)和Boosting算法(增强法),Boosting算法中最流行的一种就是Adaboost(adaptive boosting自适应增强)算法。
2.2.1 Adaboost算法原理Adaboost算法建立在Schapire的关于“弱分类器”的性能分析及Freund的关于“学习理论”的早期研究基础上于1996年被提出。由于该算法不要求事先知道弱学习算法预测精度的下限,只要求基本分类器的正确识别率略大于“随机猜想”即可(也被称为弱分类器),因此,能够更好地适用于实际问题[22]。它具有泛化错误率低、易编码等优点,被评为数据挖掘十大算法之一。
Adaboost算法的基本思想是利用大量的分类能力一般的弱分类器通过一定方法叠加起来,构成一个分类能力很强的强分类器[23]。
Adaboost算法描述如下:
输入:弱分类器设计方法C;训练次数(弱分类器个数)T;样本集X={(x1,y1),(x2,y2),…,(xN,yN)}。其中,xi为第i个样本的气象特征向量;yi={-1,1}表示第i个样本的类别标号:-1表示未发生舞动,1表示有发生舞动;N为样本数。
初始化:样本权值分布w1(i)=1/N,i=1,2,…,N。
当t=1,2,…,T时:
1)根据wt(i)从X中进行有放回的抽样生成新的样本集合Xt。
2)在Xt上训练弱分类器Ct(X),并用Ct(X)对原始训练样本集X进行分类。
3)计算Ct(X)的分类错误率。
εt=N∑i=1ωt(i)I(Ct(xi)≠yi) | (2) |
4)计算Ct(X)的系数。
at=12ln(1−εtεt) | (3) |
5)更新权值分布。
ωt+1(i)=ωt(i)Zt×{e−at,Ct(xi)=yieat,Ct(xi)≠yi=ωt(i)Zt⋅e−atyiCt(xi),i=1,2,⋯,N | (4) |
式中:Zt=N∑i=1ωt(i)⋅e−atyiCt(xi)是归一化因子,可使得N∑i=1ωt+1(i)=1。
6)最终分类器。
y=C(X)=sgn[T∑t=1atCt(X)] | (5) |
作为弱分类器,简单分类器的效果更好,因此,笔者选择最常用的单层决策树作为弱分类器。该决策树仅基于单个输入特征并采用阈值划分方法来做决策,即只有一个节点,且由于这棵树只有一次分裂过程,与树桩形似,因此,它也被称作决策桩(decision stump)[24]。
决策桩构造的最关键问题是如何判断阈值划分结果的好坏,以便选择最佳分割点。目前,一般采用基于信息量或错误率(如Gini指标)的衡量准则,而文献[25]指出:Gini不纯度指标比信息量指标性能更好,且计算方便,其最大的特点是计算时只需考虑类值在被划分时每一部分的分布情况。因此,笔者采用Gini指标来评估分割规则的优劣程度,对于包含c个类别的数据集S,其定义[26, 27]如下:
gini(S)=1−c∑j=1p2j | (6) |
gini(S,rule)=n1n⋅gini(S1)+n2n⋅gini(S2) | (7) |
对于一个数值型属性,基于Gini指数的决策桩分割思想是:在遍历所有可能的分割方法后,选择使评估值Gini(S,rule)达到最小的作为此节点处的最优划分规则[28]。其流程描述为
1) 对数值型属性的样本值进行排序,假设排序后的结果是(x1,y1),(x2,y2),…,(xn,yn)。
2) 由于分割只发生在2个数据点间,因此,通常取中点(xi+xi+1)/2作为分割点,然后从小到大依次取不同的分割点,并计算各分割规则的Gini值。
3) 取使Gini值最小的点作为最佳分割点。
2.3 基于Adaboost算法的输电线舞动预警笔者提出的基于Adaboost方法的输电线舞动预警的流程如图 1所示。
![]() |
图1 输电线舞动预警流程图 Fig. 1 Flow chart of early warning for transmission line galloping |
对由气象预报输入向量x所得的输电线舞动预测结果y,可用式(8)计算其置信度:
margin(x,y)=∑tatCt(x)∑t|at| | (8) |
最后,根据舞动预测结果的margin值,暂按表 2所示设定输电线舞动的预警等级。
![]() |
表2 输电线舞动预警等级表 Table 2 Early warning levels of transmission line galloping |
以河南省电力公司及其“输电线路舞动防治技术”国网重点实验室统计的近3年来110 kV及以上电压等级输电线路覆冰舞动调查情况为例,根据其中类型为220 kV电压等级、双分裂、大截面(S>300)导线(由于记录数据缺乏线路的具体档距,因此,这里暂不分档距大小)的舞动情况,利用Matlab 7.6.0 工具对笔者所述方法的分类结果进行编程验证。具体为选取导线出现覆冰情况下的其中“25(舞动)+10(未舞动,指疑似发生舞动但巡线并未发现任何损害)”组数据构成训练集,“7(舞动)+5(未舞动)”组数据构成测试集。验证结果如图 2所示,方法与其他常用分类算法的效果对比见表 3。
![]() |
图2 文中方法对输电线舞动预测分类效果图 Fig. 2 Effect of classification of transmission line galloping using Adaboost method |
![]() |
表3 不同分类算法的效果对比 Table 3 Comparison of different classification methods |
由图 2和表 3可以看出,所述方法在河南省220 kV大截面分裂导线舞动样本训练集和测试集上都表现出较好的稳定性和较高的正确率,分别达到100.00%和83.3%,与其他常用分类方法相比都有所提高,适用于输电线路舞动预警。
笔者还采用Adaboost学习算法的置信度来近似舞动发生的概率,并据此进行预警。由表 4可以看出:①测试集中分类错误的第9、12两组数据的置信度均较小(0.191 1和0.376 2),但是某些分类正确的数据组(如第3、10、11组)的置信度有可能更小;②就现有结果来看,置信度大于40%的数据均分类正确(这也是本文按表 2所示设定预警等级的重要依据)。因此,该值可作为划分处理预警目标优先度的参考标准。
![]() |
表4 测试集上的分类置信度 Table 4 Classification confidence of test set |
笔者提出将输电线舞动预警问题归结为有监督学习下的分类预测问题的总体思路和提出以引发线路舞动的气象环境因素的预报值为输入,以Adaboost集成学习法作为分类器算法进行导线舞动预警的技术方法都是科学和有效的。对置信度大于40%的Ⅰ、Ⅱ级舞动预警线路应重点防控、加强巡视监测、做好调度应急预案;同时,对置信度较小的Ⅲ级舞动预警也不可掉以轻心,应做好充分准备,尽力将输电线路舞动的损害降至最低。值得注意的是,输电线路舞动的发生是线路系统能量不断积累的过程,因此,除了用具体时刻的温度、风速等气象参数描述线路舞动时的外部环境信息,还应加入刻画积累过程的物理量,如持续时间、覆冰厚度等,这也是下一步研究的方向。
[1] | 郭应龙,李国兴,尤传永.输电线路舞动[M].北京:中国电力出版社,2003. GOU Yinglong,LI Guoxing,LONG Chuanyong.Transmission line dancing[M].Beijing:China Electric Power Press,2003.(in Chinese)(![]() |
[2] | 严波,李文蕴,张宏雁,等.一种模拟覆冰双分裂导线舞动的数值方法[J].重庆大学学报,2009,32(7):787-792. YAN Bo,LI Wenyun,ZHANG Hongyan,et al.Numerical method for galloping of icing twin bundle conductors [J].Journal of Chongqing University,2009,32(7):787-792.(in Chinese)(![]() |
[3] | 刘小会.覆冰导线舞动非线性数值模拟方法及风洞模型试验[D].重庆:重庆大学,2011. LIU Xiaohui.Ice conductor wave nonlinear numerical simulation and wind tunnel model test[D].Chongqing:Chongqing University,2011.(in Chinese)(![]() |
[4] | 庄为栋,马向辉,胡晓军.输电线路舞动技术探讨[J].湖州师范学院学报,2013,35(3):45-49. ZHUANG Weidong,MA Xianghui,HU Xiaojun.A discussion on the transmission line galloping technology[J].Journal of Huzhou Teachers College,2013,35(3):45-49.(in Chinese)(![]() |
[5] | 王少华.架空输电线路导线舞动及防治技术分析[J].高压电器,2010,46(12):63-67. WANG Shaohua.Analysis of conductor galloping phenomena on overhead transmission lines and anti-galloping techniques [J].High Voltage Apparatus,2010,46(12):63-67.(in Chinese)(![]() |
[6] | 朱宽军,刘彬,刘超群,等.特高压输电线路防舞动研究[J].中国电机工程学报,2008,28(34):12-20. ZHU Kuanjun,LIU Bin,LIU Chaoqun,et al.Research on anti-galloping for UHV transmission line[J].Proceedings of the CSEE,2008,28(34):12-20.(in Chinese)(![]() |
[7] | 赵汝祥,胥婷,王德洲.输电线路舞动概况及故障分析[J].价值工程,2012,31(9):18-19. ZHAO Ruxiang,XU Ting,WANG Dezhou.The situation and failure analysis of conductor galloping on transmission line [J].Value Engineering,2012,31(9):18-19.(in Chinese)(![]() |
[8] | 胡景,严波,张宏雁,等.覆冰导线舞动数值仿真分析[J].重庆大学学报,2010,33(3):76-81. HU Jing,YAN Bo,ZHANG Hongyan,et al.Numerical simulation on galloping of iced conductors[J].Journal of Chongqing University,2010,33(3):76-81.(in Chinese)(![]() |
[9] | Hartog J P D.Transmission line vibration due to sleet[J].Transactions of the American Institute of Electrical Engineers,1932,51(4):1074-1076.(![]() |
[10] | Nigol O,Clarke G J.Conductor galloping and control based on torsional mechanism[C]//IEEE transactions on power apparatus and systems.345 E 47TH ST,New York,NY 10017-2394:IEEE-INST electrical electronics engineers inc,1974(6):1729-1729.(![]() |
[11] | Yu P,Desai Y M,Shah A H,et al.Three-degree-of-freedom model for galloping (Part I:formulation and Part II:solution) [J].Journal of Engineering Mechanics,1993,119(12):2404-2448.(![]() |
[12] | 尤传永.导线舞动稳定性机理及其在输电线路上的应用[J].电力设备,2004,5(6):13-17. YOU Chuanyong.Stability mechanism of conductor galloping and its application on transmission line[J].Electrical Equipment,2004,5(6):13-17.(in Chinese)(![]() |
[13] | 石吉汉,吴继云.导线舞动的防治[J].电力建设,2005,26(12):39-42. SHI Jihan,WU Jiyun.Prevention of conductor galloping[J].Electric Power Construction,2005,26(12):39-42.(in Chinese)(![]() |
[14] | 李新民,朱宽军,李军辉.输电线路舞动分析及防治方法研究进展[J].高电压技术,2011,37(2):484-490. LI Xinmin,ZHU Kuanjun,LI Junhui.Review on analysis and prevention measures of galloping for transmission line[J].High Voltage Engineering,2011,37(2):484-490.(in Chinese)(![]() |
[15] | LIU Jianjun,ZHU Kuanjun,WANG Jingchao,et al.Application and Research on Conductor Galloping Online-monitoring in Power Transmission Lines [C]//2010 International Conference on Remote Sensing (ICRS),[S.L.]:[s.n.],2010:625-628.(![]() |
[16] | 黄新波,陶保震,冯玲.基于光流场的输电导线舞动监测方法[J].电力自动化设备,2012,32(7):99-103. HUANG Xinbo,TAO Baozhen,FENG Ling.Transmission line galloping monitoring based on optical flow field[J].Electric Power Automation Equipment,2012,32(7):99-103.(in Chinese)(![]() |
[17] | 黄经亚.架空送电线路导线舞动的分析研究[J].中国电力,1995(2):21-26. Huang Jingya.Analysis of overhead transmission line analysis galloping[J].Electric Power,1995(2):21-26.(in Chinese)(![]() |
[18] | 中华人民共和国住房和城乡建设部.GB 50009-2012 建筑结构荷载规范[S].北京:中国建筑工业出版社,2012. Ministry of Housing and Urban-Rural Development of the People's Republic of China.GB 50009-2012 load code[S].Beijing:China Building Industry Press,2012.(in Chinese)(![]() |
[19] | 钱时惕.规律及其三种主要形式——科学与人文漫话之十[J].物理通报,2010(7):85-87. Qian Shixi.Law and its three primary forms[J].Physics Bulletin,2010(7):85-87.(in Chinese)(![]() |
[20] | Tom M.Mitchell.机器学习(计算机科学丛书)[M].曾华军,张银奎,等,译.北京:机械工业出版社,2003. Tom M.Mitchell.Machine Learning (Computer science books) [M].ZENG Huajun,ZHANG Yingkui,et al.Beijing:China Machine Press,2013.(in Chinese)(![]() |
[21] | 张春霞.集成学习中有关算法的研究[D].西安:西安交通大学,2010. ZHANG Chunxia.The integrated learning about the research of the algorithm[D].Xi'an:Xi'an Jiaotong University,2010.(in Chinese)(![]() |
[22] | 孙即祥.现代模式识别[M].2版.北京:高等教育出版社,2008. SUN Jixiang.Modern pattern recognition [M].2th ed.Beijing:Higher Education Press,2008.(in Chinese)(![]() |
[23] | 何光辉,张太平.局部匹配的人脸识别方法[J].重庆大学学报,2012,35(12):133-138. HE Guanghui,ZHANG Taiping.Face recognition method with local matching[J].Journal of Chongqing University,2012,35(12):133-138.(in Chinese)(![]() |
[24] | 哈林顿.机器学习实战[M].李锐,李鹏,等,译.北京:人民邮电出版社,2013. Harrington P.Machine Learning in Action[M].LI Rui,LI Peng,et al.Beijing:Posts and Telecom Press,2013.(in Chinese)(![]() |
[25] | 万源.一个基于SLIQ算法的模型及应用决策桩选择[J].信息技术,2005(12):60-62. WAN Yuan.A model and its application based on SLIQ algorithm [J].Information Technology,2005(12):60-62.(in Chinese)(![]() |
[26] | 谢金梅,王燕妮.决策树算法综述[J].软件导刊,2008,7(11):83-85. XIE Jinmei,WANG Yanni.The decision tree algorithm [J].Software Guide,2008,7(11):83-85.(in Chinese)(![]() |
[27] | 王康,孙宏斌,张伯明,等.基于二维组合属性决策树的暂态稳定评估[J].中国电机工程学报,2009,29(Sup 1):17-24. WANG Kang,SUN Hongbin,ZHANG Boming,et al.Transient stability assessment based on 2D combined attribute decision tree[J].Proceedings of the CSEE,2009,29(Sup 1):17-24.(in Chinese)(![]() |
[28] | 陈云樱,吴积钦,徐可佳.决策树中基于基尼指数的属性分裂方法[J].微机发展,2004,14(5):66-68. CHEN Yunying,WU Jiqin,XU Kejia.Using gini-index for attribute selection in decision trees[J].Microcomputer Development,2004,14(5):66-68.(in Chinese)(![]() |