电线积冰及路面温度研究的新进展

牛生杰 王天舒 吕晶晶 周悦 王元

摘要 2009年起，研究组开展了电线积冰野外观测试验及道路结冰规律的观测研究，将电线积冰观测研究从传统的积冰气象条件和积冰厚度观测提升到了积冰气象条件、积冰厚度加积冰天气云降水微物理综合观测研究的新高度，揭示了积冰发生的微物理机制，研究积冰增长率及其影响因子，建立的积冰增长模型较好地模拟了积冰增长过程。对沥青、水泥、土壤三种下垫面温度进行了观测研究，观测高速公路和桥梁不同路基深度的温度变化，并对桥面比路面更易结冰的现象从能量平衡方面做了理论解释。本文以作者团队取得的成果为主线，不求大而全，学习梳理相关的代表性研究成果，主要包括积冰发生频次的时空分布、积冰天气微物理特征、积冰增长率及其影响因子、积冰气象条件、积冰增长模型构建、道路结冰及积冰数值预报等几个方面，并对电线积冰未来的研究提出了建议。

关键词 电线积冰; 微物理特征; 积冰气象条件; 积冰增长模型; 路面温度; 研究进展

2008年1月10日至2月2日，我国南方地区连续发生四次罕见低温雨雪冰冻天气，灾害范围广，强度大，持续时间长，造成了极其严重的经济损失和人员伤亡。冻雨、雨夹雪、过冷雾等造成严重的电线覆冰并发生次生灾害，造成输电线塔、通信塔等倒伏，造成供电、通信中断;造成道路结冰，使得公路交通中断。电线积冰为雾凇、雨凇在导线上凝附或湿雪在导线上冻结的现象（中央气象局，2003），通常发生在过冷雾、冻雨、湿雪三种天气条件下（Kringlebotn et al.，2013），其严重时会给电力输送造成极大损失，例如电缆舞动、断线、杆塔倒塌和绝缘子闪络等现象，严重危害国民经济（陈柏堃等，2009;Niu et al.，2012）。20世纪70年代末以来，我国各地尤其是南方地区积冰事故时有发生。根据国家电网运行管理处的统计，2006年1月至2007年6 月，电线覆冰导致13次500 kV电线跳闸，在总跳闸次数中占比8.8%;导致500 kV电线4次非计划停运，在总停运次数中占比11.1%（李庆峰等，2008）。2008年冬季由于适宜的大气环流形势、天气系统及温湿层结条件，1月10日到2月2日我国南方地区发生了4次大面积冰冻雨雪灾害（杨贵名等，2008），13个省遭遇电线覆冰事故，贵州、湖南两省大面积停电，4 216条重要输电线路损毁，公路因结冰导致交通中断，造成超过1 100亿元的直接经济损失（杜骦等，2019）。

提高对积冰现象的认知，需要深入探索积冰期间气象要素变化规律以及雨雾滴微物理特征，揭示冰冻灾害形成机理，从而开发积冰模拟预报方法，提高预报预警水平。本文主要包括积冰发生频次的时空分布、积冰天气微物理特征、积冰增长率及其影响因子、积冰气象条件、积冰增长模型构建、道路结冰及积冰数值预报等几个方面，并对未来电线积冰的研究方向提出了建议。

1 电线积冰发生频次的时空分布

冰冻天气在中国的分布具有较强的区域性特征。王遵娅（2011）对中国1954—2009年雨凇和雾凇的空间分布及不同时间尺度的变率进行全面细致的分析后指出：中国的冰冻主要出现在江南、新疆、西北东部、内蒙古东部、东北、华北、淮河流域等地，年均冰冻日数约有1～5 d，5 d以上的重冰区主要分布在云南东北部、贵州、新疆北部、陕西南部、东北中部、华北东部、秦岭等地（图1）。中国冰冻、雨凇和雾凇日数在1954—2009 年间出现了显著减少的变化趋势（王遵娅，2011）。2008 年1月突如其来的低温雨雪冰冻灾害席卷中国南方多省，此次冰冻灾害不是一个局地或地区性现象，它是同期发生的亚洲大范围冰雪灾害链中的一环，在影响范围和灾害程度上是最严重的一环（丁一汇等，2008）。随后，冰冻雨雪天气研究引起广泛关注，成为研究热点。

由2011—2013 年中国民航机场的例行观测、特殊观测资料得出的冰粒降水天气分布显示：冰粒降水天气集中分布在中纬度的黄河下游、长江中下游地区。冰粒降水发生的次数与其持续时间相关性不强;冰粒降水天气在1—2月比较常见，单次冰粒降水天气过程的持续时间比较短，一般不超过1 h（李杰等，2016）。

王颖等（2012）在綜合考虑1960—2009年各种强度等级发生频率及其产生的危害强弱基础上，建立了一个积冰危险程度的量化指标。分析结果显示，中高危险区主要集中在南方地区，主要分布在长江中下游的山地丘陵地区，包括四川中南部、重庆、云南东北部、贵州、湖北西南部、湖南东南部、江西北部、安徽南部、浙江等地。

由此可见，积冰发生频次的时空分布研究，需要使用更长时间序列的资料使其系统化，给出更高分辨率的积冰发生频次时空分布，揭示其形成的天气气候规律，评估不同地区不同温度湿度条件下电线积冰可能造成的灾害等级，要特别关注雪花中的含水量数值，雪花含水量影响雪花在电线表面的黏附性，进而影响积冰增长速度，这是目前积冰成灾研究中的盲点之一。

2 电线积冰天气微物理特征

研究积冰形成机制，离不开对过冷云雾降水宏微观特征的研究，过冷液滴、湿雪是否能够在电线表面碰冻是积冰发展的关键，因此探讨积冰发生、发展和维持阶段的微物理特征显得十分重要。

2.1 过冷雾

过冷雾滴谱及其演变直接影响着积冰增长率及积冰属性的变化规律，且过冷雾滴谱微物理特征对云中积冰过程的影响更为复杂。目前，我国学者针对这一科学问题开展了一些观测试验，对造成云中积冰过程的山区云雾宏微观特征获得了一定的认识。邓雪娇等（2002）指出南岭山区的雾过程均以小滴为主，大部分液滴的粒径小于12 μm，平流因素和地形作用对微物理量的起伏变化有重要影响。山顶的雾过程主要是冷锋过境带来的云系接地形成，受平流的作用，相对于山坡更易出现浓厚稳定的雾过程（史月琴等，2006;吴兑等，2007;邓雪娇等，2007）。Niu et al.（2012）对湖北恩施山区的过冷雾微物理特征进行了探讨（表1），指出雾滴谱的微物理特征量变化范围均较小，与平流雾类似，液水含量与积冰厚度的增长率呈显著正相关关系。贾然等（2010）对湖北恩施雷达站处的一次积冰过程中雨天、雾天和雨夹雪天气条件下，雾滴谱和雨滴谱的微物理特征进行了初步的分析，发现在雨夹雪天气时积冰的增长速度要明显大于雨天和雾天。贵州山区雾过程的数浓度、液水含量和平均直径分别为 234 cm-3、0.2 g·m-3和7.5 μm，且以小滴为主，峰值直径为4 μm，峰值浓度占总浓度的 63%，粒径大于14 μm的雾滴对积冰增长的贡献明显大于小雾滴（罗宁等，2008）。当存在充足水汽供应时，山区雾的液水含量会显著增加（Hodges and Pu，2016），Drage and Hauge（2008）对沿海山区的过冷雾液水含量进行了观测和模拟，在积冰快速增长期的液水含量达到了0.4 g·m-3左右。

2.2 冻雨

冻雨、冻毛毛雨滴谱的演变直接影响着降水积冰过程中积冰属性的变化规律，大雨滴受风速和电线直径（简称线径）影响较小，更易直接碰冻在导线表面;而小雨滴在风速较大、线径较小时，则易绕过导线无法碰冻，但同时较大的风速又可以吹来更多的雨滴，增加可能碰冻的雨滴个数（Jones and Eylander，2017）。2008年初安徽潜山的一次冻雨过程雨滴谱与层状云降水类似，其雨强、雨水含量和雨滴粒径都较小，这次冻雨过程是由于小冰晶融化导致的（Chen et al.，2011）。而山区过冷降水过程宏微观物理特征量的数值则明显小于平原地区的冻雨过程，其形态以冻毛毛雨为主，持续时间为 1～3 h，雨强基本小于 0.5 mm·h-1，平均粒径则主要在 0.4～0.6 mm（Zhou et al.，2016）。黄钦等（2018）通过分析 2014年2月庐山地区的雨滴谱观测资料，发现存在低落速冻雨滴向冰粒和干雪转化的过程，Gamma分布对该种类型的降水过程拟合效果较好。

2.3 湿雪

雪花下落过程中经过温度稍高于0 ℃的层结时表面会融化，形成湿雪，湿雪引起的积冰密度较大，对输电线路影响较严重。李遥等（2019）给出了南京地区雪花谱的峰值下落末速度为1.1 m/s，且存在小雪花碰并形成大雪花及雪花结霜增长等两种大雪花形成机制，并归纳得出南京地区稳定的层状云强降雪过程的Z-I关系为Z=1 708I1.51。葛森等（2016）指出运用欧拉二维正定平流输送法可以合理地模拟出群粒子在凝华增长过程中纵横比的演变。徐戈等（2016）总结了产生冻滴的两个主要过程：半径大于80 μm 的液滴冻结形成冻滴;霰与半径大于100 μm 的液滴碰撞转化为冻滴。并发现霰的生成机制受到抑制，浓度减小;冻滴产生、繁生产生的冰晶增多;冰晶聚并形成雪，雪的浓度增加;冻滴消耗过冷水增长，导致液态水浓度减小。

利用先进的高分辨率气溶胶云降水粒子观测设备，在积冰高发区观测了积冰过程中的微物理结构，研究了冰粒子产生及增长的规律，为研究影响积冰增长的气象和微物理因子，构建积冰增长模型奠定了基础。今后的积冰微物理研究，应当与宏观天气过程结合起来，提高积冰模型对积冰物理过程描述的准确性。

3 电线积冰与气象条件的关系

3.1 地面气象要素

积冰过程的必要气象条件是：1）具有液滴可冻结的气温，即0 ℃以下;2）具有较湿润的环境空气，即相对湿度一般在85%以上;3）具有可使空气中液滴运动的动力条件，即大于1 m/s 的风速（蒋兴良和易辉，2001）。20世纪70年代起，我国电线积冰气象条件研究取得了许多有意义的结果。滕中林（1978）讨论了风速、地形、导线直径、电场等对积冰的影响，并给出不同高度积冰的变化。谭冠日（1982）进一步探讨了积冰冰厚随距地高度的变化，并分析了风与电线交角、电线直径和冰厚之间的关系。江祖凡（1983）利用1978—1981年冬季庐山云雾试验资料对电线积冰增长率与气象条件之间的关系进行了较为定量的研究。刘雪静和牛生杰（2016）指出气温和风速是影响积冰过程的关键气象要素。高压输电线出现积冰时的环境气温约为-2 ℃，而脱落时气温为-2～-1 ℃（图2），明显低于积冰架积冰发生和脱落时的气温（Zhou et al.，2012）。

由于风速计在冰冻天气下极易冻结，风速作为影响积冰过程的重要气象要素，在早期的研究中相对较少。一方面，在积冰的发生和发展阶段，风将大量过冷液滴不断地输送向电线，与电线碰冻，加速积冰的增长;另一方面，过大的风速（尤其是瞬时风速）会破坏冰层的结构，加速积冰的脱落，这一现象主要发生在积冰的维持和脱落阶段，对于雾凇积冰过程最为显著。Jones（1998）建立的冻雨积冰物理模型则直接给出了风速对积冰增长的贡献项，其对冻雨积冰的模拟效果良好，并得到了Degaetano et al.（2008）的验证和应用。Savadjiev and Farzaneh（2004）则认为风速对云中积冰过程中冰厚的增长作用强于降水积冰过程中的数值。

相对湿度对积冰过程的直接影响较弱，主要表征了积冰环境中的水汽特征。通常情况下云中积冰过程的出现和发展阶段过冷云雾长时间笼罩在观测点，相对湿度维持在100%左右，饱和的水汽条件使得小液滴能够不断生成，为积冰过程提供足够的过冷液水，而在维持和脱落阶段云雾通常开始消散，相对湿度也会逐渐减小。降水积冰过程则与相对湿度的变化关系较小。当相对湿度较小且气温低于0℃时，积冰的升华作用很可能是导致积冰脱落消失的重要原因之一，且这一现象多出现在北美的外场观测中（Druez et al.，1995;Ryerson and Kenyon，1998）。

3.2 大气层结

大气层结特征也与冻雨天气的形成密切相关。按照大气层结特点，冻雨形成途径分为“冰相融化”和“过冷雨滴”（Roberts and Stewart，2008;陶玥等，2012，2013）。“冰相融化” 过程中，冰粒子从高空下落到温度在0 ℃以上的大气层中，融化为雨滴，随后下落到近地面再冻结层成为过冷雨滴，碰撞到物体上即发生冻结。 “过冷雨滴” 过程为地面至高空气温均低于0 ℃，冻雨由液滴碰并形成。

发现两种冻雨形成机制后，冻雨层结特征成为冻雨研究的重要方面。Zerr（1997）分析了美国和加拿大南部34次冻雨、冻雨及冰丸、冰丸个例的雷达和无线电探空仪资料，地面降水类型对融化特性更为敏感，并用两种热传递模型模拟融化和再冻结，模拟结果较符合地面降水观测事实。Thériault et al.（2006）用一维云模式调整大气温湿廓线模拟冬季不同类型降水，增加地面温度、融化层温度和厚度，同时减少再冻结层厚度，导致冻雨累积降水量增加，冰丸量减少。2008年初我国南方冰灾期间倾斜的锋区造成我国从北向南逆温层和融化层高度、厚度、强度的差异，依次形成了雪、冻雨和雨（曾明剑等，2008），可见逆温层高度、厚度和强度对降水类型的决定性（Bourgouin，2000;Kovacik et al.，2011）。凍雨形成中逆温层对聚集水汽起到重要作用，探空曲线表明逆温层以上相对湿度随高度快速下降，湿度大值区位于逆温层以下（李登文等，2009）。

3.3 积冰的梯度观测

Wang et al.（2019）对庐山一次寒潮积冰过程（2016年1月19—25日）进行了积冰梯度观测，研究了不同距地高度（1.5、10 m）微环境差异及其造成的积冰差别（图3）。1.5 m和10 m高度的气温T1.5、T10全过程均值分别为-6.2、-7.0 ℃。全过程1.5 m高度的风速V1.5总体低于10 m高度的风速V10，均值分别为 1.5、2.8 m·s-1。

对比低处、高处积冰的起步和演进特征。太阳辐射加强时，低处积冰受直接太阳辐射、显著提高的T1.5、地面长波辐射和反射太阳辐射的多重影响，冰厚显著下降，夜晚冰厚回升，全过程冰厚最高达1.2 mm。高处积冰生长优势明显：积冰起步方面，高处适宜的低温使积冰起步早于低处 6 h，为高处后续冰累积奠定了良好的基础，冰筒有了更大的表面得以粘附更多降水粒子，逐渐加大了其与低处的冰厚差异;冰增长方面，高处经过两个增长阶段冰厚达到 18.0 mm，两个增长阶段都对应着两高度风速差大值区，可见风对高處冰增长起到了关键作用;并且高处积冰在强太阳辐射时的损失较少，有效保存了已有的积冰。对比低处、高处积冰的结束方式，停止增长后，低处积冰直接脱落;高处冰厚维持14 h 后在强太阳辐射和强风影响下逐渐融化。

4 电线积冰模型及预报

基于野外观测资料的研究，揭示了积冰增长的基本物理过程。基于对成冰基本原理的掌握，众多学者开始尝试建立数学物理模型实现电线积冰模拟，逐步完善冰冻雨雪灾害预报预警方法。依据积冰形成机制，可将积冰分为冻雨型积冰、湿雪型积冰和综合型积冰，根据三类积冰形成机制，分别建立了三种不同类型的致冰模型。

4.1 冻雨积冰模型

冻雨积冰出现条件为湿球温度低于0 ℃、有液态降水发生（Makkonen，2000）。Jones（1998）将冻雨积冰的物理过程简化，以降水率和风速为参数提出了冻雨型积冰增长模型，并讨论了与其他模型的差异及差异的原因。该模型计算简便，在冻雨积冰预报时使用较广泛。温华洋等（2011）通过人工神经网络和逐步多元线性回归两种统计方法对电线积冰资料进行了分析，资料为安徽有电线覆冰观测的15个气象站自建站至2008年的观测资料，为安徽省电线积冰灾害建立了三个区域的标准冰厚估计模式，相比人工神经网络方法，逐步多元线性回归方法对建立安徽省电线积冰标准冰厚估算模式更为合适，对非高山站具有较好的模拟效果。Abdelaal et al.（2019）选取不同电缆倾斜角度进行风洞结冰试验，指出电缆倾斜度影响冰厚及冰钉长度，考虑电缆的倾斜度、积冰增长方向和形状，开发了冻雨积冰模型，通过美国托莱多的两次积冰事件对模型进行了评估，模拟与观测结果较为一致。

4.2 湿雪积冰模型

湿雪积冰出现条件为有降雪发生、湿球温度在0 ℃以上（Makkonen，2000）。Poots and Skelton （1995）对雪粒下落轨迹与电线夹角以及积冰形状进行了研究，用余弦定理计算雪质量增长率，对湿雪积冰过程进行了模拟。Makkonen（1989）以湿球温度作为判断湿雪出现的指标，基于湿雪积冰增长的基本物理过程，以能见度等作为输入量，计算湿雪质量增长率，提出了湿雪型积冰增长模型。Makkonen and Wichura （2010）对Makkonen（1989）模型进行了改进，在计算雪浓度时分白天和夜间两种情况来讨论，并将判断湿雪出现的标准由湿球温度大于0 ℃变为湿球温度大于-0.2 ℃。用改进后的模型对2005年11月德国西北部湿雪积冰灾害进行了模拟，取得了较好的效果。Kringlebotn et al.（2013）利用冰岛南部50 a湿雪积冰观测资料和WRF输出量对Makkonen（1989）湿雪模型进行了改进和讨论，用WRF输出量中液态降水和固态降水的比例来判断湿雪，并对比了三种捕获率的计算方法，模拟取得了较好的效果。

4.3 综合积冰模型

综合积冰模型可模拟多种类型的积冰增长。Makkonen （2000）在假设覆冰形状为圆柱体的基础上，考虑了水成物粒子与导线之间碰并的微物理过程与热平衡过程，建立了一套适用于多种覆冰情况的覆冰预报模型，被大部分组织和学者推荐使用。Niu et al.（2012）计算积冰增长中的碰撞率、捕获率和冻结率等参量，在计算捕获率时区别干增长和湿增长，模拟了湖北恩施积冰厚度，模拟结果能较好地反映实际冰厚变化。邓芳萍等（2017）将雨凇、雾凇增长模型与融化、升华模型相结合，建立了一个基于常规气象资料的小时标准冰厚模型，用其模拟浙江2008和2013年两次电线积冰事件，并用事故线路调查资料、电线积冰观测站和模拟导线拉力监测点的观测资料进行验证分析，模拟值与事故线路最大标准冰厚观测值相关性通过0.01的显著性水平检验，模型较准确地抓住了天气过程对电线积冰的影响。Szilder （2019）研发了新分析模型来预测电线积冰，对垂直、倾斜和水平电缆上冻结的积雪和冻雨的分析表明，两个最有影响力的参数是风速与粒子下落末速度之比以及冰筒圆柱体倾斜角。

4.4 基于中尺度数值模式的积冰预报

随着对积冰形成的物理机制认识的不断深化，可以利用WRF等数值天气诊断预报模式进行冰冻雨雪天气的模拟和积冰预报。Musilek et al.（2009）、Pytlak et al.（2010）、Hosek et al.（2011）将WRF模式与积冰增长模型相结合，用WRF模式输出量建立了冻雨判别指标，发展了积冰增长预报系统（IAFS）。杜骦等（2019）、刘善峰等（2019）将IAFS积冰预报系统分别应用到我国河南省和贵州省积冰灾害模拟中，得到了初步应用。Davis et al.（2014）将积冰增长、融化模型与WRF模式相结合，并用其预测风力发电场风力机积冰，发现该方法对积冰短时预报效果良好，微物理与边界层方案的选择对冰重模拟有较大影响。Lamraoui et al.（2014）将积冰增长模型与电力损耗模型相结合，对电力损耗进行了估计。结果表明，冻结比为0.88时对应电力损失的峰值。Barszcz et al.（2018） 研究了加拿大气象中心数值天气预报系统严重低估的一次冰冻事件，该模式直接根据Milbrandt-Yau微物理方案对降水类型进行预报，冻雨量的低估主要是由于融化的雪和霰的过度再冻结，通过降低雨-霰收集效率和设置-5 ℃的温度阈值，微物理方案中过度的雨-霰聚集和霰形成得到控制。Kmrinen et al.（2018）进行了七区域气候模型模拟，评估了欧洲中等（RCP4.5）和强（RCP8.5）未来排放背景下气候变化对冻雨和冰冻事件的影响。在RCP8.5排放背景下，到本世纪末，西欧、中欧和东南欧的冻雨频率下降了20%～55%，而该大陆的北部和东北部则增加了0～50%，在RCP4.5排放背景下变化相似但量级偏小。

已有的积冰研究按照积冰形成的三种机制即冻雨积冰、湿雪积冰和过冷雾积冰分别建立了积冰模型，在积冰监测预警中发挥了一定的作用。影响积冰增长的因子很复杂，首先要揭示有利于形成积冰的大尺度环流形势天气系统，对冰冻雨雪天气成因的研究要把宏观和微观结合起来，还要考虑地理条件，输电线路的物理属性等因素，在具有代表性的积冰高发地区进行多站点、长时间、高分辨率的智能化积冰增长综合观测，建立具有中国自主知识产权的积冰增长模型，并进行数值预报试验，以提高积冰预警水平;研究气候变化条件下积冰灾害发生趋势，为防灾减灾服务。

5 道路结冰及路桥温度观测试验

道路结冰和电线积冰为伴生灾害，道路结冰亦可参考电线积冰的研究方法和形成原理。道路结冰严重影响交通运输安全（舒斯等，2019），准确预报路面温度和路面状况对降低冬季道路维护成本和预防减少交通事故发生具有重要作用，这就需要研究影响路面温度變化的因素，揭示下垫面特性对道路结冰气象条件的影响机制（李蕊等，2011;牛生杰等，2011;吕晶晶等，2013）。牛生杰等（2011）以水泥、土壤和沥青三种下垫面温度作为研究对象，在结冰预报模型中充分考虑三种不同下垫面特性，研究结果显示下垫面温度受不同材料性质引起的误差较小。吕晶晶等（2013）提出由于桥梁材质混凝土的比热容低于路面材质沥青，考虑到桥梁与空气的接触面大于路面，因此夜间桥梁以感热形式与空气进行热交换的能量应高于相同气象条件的路面，解释了桥面比路面更易结冰的现象，提出云量是影响路、桥面温度变化的主要因素之一，并根据总云量定义晴空、多云和阴天三种天空状况。

路面温度与路面状况预报技术目前以统计分析法（牛生杰等，2011）和基于地表能量辐射平衡理论的物理模式预报方法（李蕊等，2011;吕晶晶等，2013）两种为主。李蕊等（2011）建立不同下垫面温度与气象要素的多元回归模型，针对路面温度取得了较好的预报效果，并指出路面0 ℃低温与路面积水对结冰的产生和维持影响密切。田华等（2019）指出统计方法虽然运算简便，预测结果比较精准，但受路面实测数据观测数量和时间长度的限制，存在地区适用性问题。

基于地表能量辐射平衡理论的数值模式预报方法及其改进，本质就是以能量平衡方程为基础，考虑太阳短波辐射、大气和地面的长波辐射、潜热、感热传输以及下垫面的热传导等能量之间的平衡，建立数值预报模型，从而探究路面温度和路面状况变化的物理规律。牛生杰等（2011）在基本能量平衡方程的基础上，考虑水汽、气溶胶、浮尘以及云等对太阳短波辐射的吸收和散射，建立了一种适用于实际预报冬季高速公路路面结冰状况的下垫面温度预报模型。吕晶晶等（2013）利用武英高速凤凰关水库自动气象站逐时路面、桥面温度和常规气象资料，分析了气温、路面温度和桥面温度的变化规律。在冬季典型天空状况（晴空、阴天） 和天气过程（雾、降雨、降雪） 中，将实测值与下垫面能量平衡模型计算值对比，验证模型计算的准确性和可靠性。应用该模型分析路面、桥面辐射能量收支变化，取得了较好的效果。

由上可知，基于地表能量辐射平衡理论的数值模式预报方法具有较好的普适性和较高的时空分辨率，但是建立的模型形式复杂，输入参数多且不易获得（秦健和孙立军，2005）;而统计分析法克服了参数多、计算繁琐等缺点，但存在地区适用性问题。因此，路面温度与路面状态预报技术研究以两者相互补充为主（田华等，2019）。

6 小结与展望

本文从积冰发生频次的时空分布、积冰形成的温度湿度层结特征、积冰气象条件、微物理过程、积冰模型构建及数值预报试验、道路结冰规律等对已有成果进行了归纳总结，这些成果揭示了水凝物碰冻在电线表面并持续累积的机制，并建立了道路温度预报方法，在气象服务中发挥了较好的作用。

目前的积冰观测研究，主要是在积冰多发地区的气象站开展，实际上属于积冰模拟观测，主要目的是揭示积冰形成的宏微观气象条件和积冰演变规律。要进一步提高积冰研究针对性，提高积冰预警准确率，需要与电力公司密切合作，在高压输电线路上安装自动气象站、拉力器和摄像头，在真实电力输送环境下，针对积冰形成的不同气候背景、天气系统、微物理过程、地理环境等因子，进行多站点、长时间、高分辨率的智能化积冰增长观测研究，积累更多具有多方面代表性的资料，构建具有中国知识产权的积冰增长模型，建设符合我国天气气候条件、地理环境等因素的积冰数值预报系统，提高积冰灾害监测预警水平，为冰冻雨雪灾害防灾减灾服务，将科研成果运用在保障电力输送网安全稳定高效运行的服务系统中。

参考文献（References）

Abdelaal A，Nims D，Jones K，et al.，2019.Prediction of ice accumulation on bridge cables during freezing rain：a theoretical modeling and experimental study[J].Cold Reg Sci Technol，164：102782.doi：10.1016/j.coldregions.2019.102782.

Barszcz A，Milbrandt J A，Thériault J M，2018.Improving the explicit prediction of freezing rain in a kilometer-scale numerical weather prediction model[J].Wea Forecasting，33（3）：767-782.doi：10.1175/waf-d-17-0136.1.

Bourgouin P，2000.A method to determine precipitation types[J].Wea Forecasting，15（5）：583-592.doi：10.1175/1520-0434（2000）015<0583：amtdpt>2.0.co;2.

Chen B J，Hu W，Pu J P，2011.Characteristics of the raindrop size distribution for freezing precipitation observed in southern China[J].J Geophys Res：Atmos，116（D6）：D06201.doi：10.1029/2010JD015305.

陈柏堃，郜庆林，吴明江，2009.金华近56年电线积冰气候特征及灾害防御[J].气象，35（8）：85-90. Chen B K，Gao Q L，Wu M J，2009.Climatic characteristics analysis and the disaster defense of wire icing in Jinhua City in the recent 56 years[J].Meteor Mon，35（8）：85-90.（in Chinese）.

Davis N，Hahmann A N，Clausen N E，et al.，2014.Forecast of icing events at a wind farm in Sweden[J].J Appl Meteorol Climatol，53（2）：262-281.doi：10.1175/jamc-d-13-09.1.

DeGaetano A T，Belcher B N，Spier P L，2008.Short-term ice accretion forecasts for electric utilities using the weather research and forecasting model and a modified precipitation-type algorithm[J].Wea Forecasting，23（5）：838-853.doi：10.1175/2008waf2006106.1.

鄧芳萍，康丽莉，姜瑜君，等，2017.基于常规气象资料的小时标准冰厚模型及验证[J].应用气象学报，28（2）：142-156. Deng F P，Kang L L，Jiang Y J，et al.，2017.An hourly standard ice thickness model using conventional meteorological data with its validation[J].J Appl Meteor Sci，28（2）：142-156.doi：10.11898/1001-7313.20170202.（in Chinese）.

邓雪娇，吴兑，叶燕翔，2002.南岭山地浓雾的物理特征[J].热带气象学报，18（3）：227-236. Deng X J，Wu D，Ye Y X，2002.Physical characteristics of dense fog at nanling mountain region[J].J Trop Meteor，18（3）：227-236.（in Chinese）.

邓雪娇，吴兑，唐浩华，等，2007.南岭山地一次锋面浓雾过程的边界层结构分析[J].高原气象，26（4）：881-889. Deng X J，Wu D，Tang H H，et al.，2007.Analyses on boundary layer structure of a frontal heavy fog process in nanling mountain area[J].Plateau Meteor，26（4）：881-889.（in Chinese）.

丁一汇，王遵娅，宋亚芳，等，2008.中国南方2008年1月罕见低温雨雪冰冻灾害发生的原因及其与气候变暖的关系[J].气象学报，66（5）：808-825. Ding Y H，Wang Z Y，Song Y F，et al.，2008.Causes of the unprecedented freezing disaster in January 2008 and its possible association with the global warming[J].Acta Meteorol Sin，66（5）：808-825.（in Chinese）.

Druez J，Louchez S，McComber P，1995.Ice shedding from cables[J].Cold Reg Sci Technol，23（4）：377-388.doi：10.1016/0165-232X（94）00024-R.

Drage M A，Hauge G，2008.Atmospheric icing in a coastal mountainous terrain.Measurements and numerical simulations，a case study[J].Cold Reg Sci Technol，53（2）：150-161.doi：10.1016/j.coldregions.2007.12.003.

杜骦，周宁，韩永翔，等，2019.河南省一次冻雨过程中电线积冰厚度模拟[J].气象，45（5）：641-650. Du S，Zhou N，Han Y X，et al.，2019.Simulation of wire ice thickness during a freezing rain process in Henan Province[J].Meteor Mon，45（5）：641-650.doi：10.7519/j.issn.1000-0526.2019.05.005.（in Chinese）.

葛森，孙继明，牛生杰，2016.典型非球形冰晶粒子的凝华增长数值模拟试验[J].大气科学，40（3）：617-629. Ge S，Sun J M，Niu S J，2016.Numerical simulations of ice crystal growth in the water vapor deposition process for typical ice habits[J].Chin J Atmos Sci，40（3）：617-629.doi：10.3878/j.issn.1006-9895.1508.15151.（in Chinese）.

Hodges D，Pu Z X，2016.The climatology，frequency，and distribution of cold season fog events in northern Utah[J].Pure Appl Geophys，173（9）：3197-3211.doi：10.1007/s00024-015-1187-6.

Hosek J，Musilek P，Lozowski E，et al.，2011.Forecasting severe ice storms using numerical weather prediction：the March 2010 Newfoundland event[J].Nat Hazards Earth Syst Sci，11（2）：587-595.

黄钦，牛生杰，吕晶晶，等，2018.庐山一次积冰天气过程冻雨滴谱及下落末速度物理特征个例研究[J].大气科学，42（5）：1023-1037. Huang Q，Niu S J，Lü J J，et al.，2018.Physical characteristics of freezing raindrop size distribution and terminal velocity in two ice weather cases in Lushan area[J].Chin J Atmos Sci，42（5）：1023-1037.（in Chinese）.

賈然，牛生杰，李蕊，2010.鄂西电线积冰微物理观测研究[J].气象科学，30（4）：481-486. Jia R，Niu S J，Li R，2010.Observational study on microphysical characteristics of wire icing in west Hubei[J].Sci Meteorol Sin，30（4）：481-486.doi：10.3969/j.issn.1009-0827.2010.04.008.（in Chinese）.

蒋兴良，易辉，2001.输电线路覆冰及防护[M].北京：中国电力出版社. Jiang X L，Yi H，2001.Transmission line icing and protection[M].Beijing：China Electric Power Press.（in Chinese）.

江祖凡，1983.电线积冰增长速度的研究[J].科学通报，28（15）：928-931. Jiang Z F，1983.The growth rate of ice accretion on power line [J].Chin Sci Bull，28（15）：928-931.（in Chinese）.

Jones K F，1998.A simple model for freezing rain ice loads[J].Atmos Res，46（1/2）：87-97.doi：10.1016/S0169-8095（97）00053-7.

Jones K，Eylander J，2017.Vertical variation of ice loads from freezing rain[J].Cold Reg Sci Technol，143：126-136.doi：10.1016/j.coldregions.2017.07.008.

Kmrinen M，Hyvrinen O，Vajda A，et al.，2018.Estimates of present-day and future climatologies of freezing rain in Europe based on CORDEX regional climate models[J].J Geophys Res：Atmos，123（23）：13291-13304.doi：10.1029/2018JD029131.

Kovacik C，Hocker J，Shafer M，2011.An analysis of southern U.S.ice storm frequency from 2000—2009[C]//Sixth Symposium on Policy and Socio-economic Research.American Meteorological Society，Seattle.

Kringlebotn N B E，gústsson H，Somfalvi-Tóth K，2013.Modeling wet snow accretion on power lines：improvements to previous methods using 50 years of observations[J].J Appl Meteor Climatol，52（10）：2189-2203.doi：10.1175/jamc-d-12-0332.1.

Lamraoui F，Fortin G，Benoit R，et al.，2014.Atmospheric icing impact on wind turbine production[J].Cold Reg Sci Technol，100：36-49.doi：10.1016/j.coldregions.2013.12.008.

李登文，喬琪，魏涛，2009.2008年初我国南方冻雨雪天气环流及垂直结构分析[J].高原气象，28（5）：1140-1148. Li D W，Qiao Q，Wei T，2009.Analyses on the freezing rain and snow weather circulation and vertical structure of Southern China in early 2008[J].Plateau Meteor，28（5）：1140-1148.（in Chinese）.

李杰，郭学良，盛日峰，等，2016.我国冰粒降水天气的观测特征统计分析[J].大气科学学报，39（3）：349-360. Li J，Guo X L，Sheng R F，et al.，2016.Statistical analysis of observed properties of ice-pellet precipitation in China[J].Trans Atmos Sci，39（3）：349-360.doi：10.13878/j.cnki.dqkxxb.20140418032.（in Chinese）.

李庆峰，范峥，吴穹，等，2008.全国输电线路覆冰情况调研及事故分析[J].电网技术，32（9）：33-36. Li Q F，Fan Z，Wu Q，et al.，2008.Investigation of ice-covered transmission lines and analysis on transmission line failures caused by ice-coating in China[J].Power Syst Technol，32（9）：33-36.（in Chinese）.

李蕊，牛生杰，汪玲玲，等，2011.三种下垫面温度对比观测及结冰气象条件分析[J].气象，37（3）：325-333. Li R，Niu S J，Wang L L，et al.，2011.Comparison of three kinds of underlying surface temperatures and analysis of the road icing meteorological conditions[J].Meteor Mon，37（3）：325-333.（in Chinese）.

李遥，牛生杰，吕晶晶，等，2019.2018年冬季南京三次暴雪过程微物理特征分析[J].大气科学，43（5）：1095-1108. Li Y，Niu S J，Lü J J，et al.，2019.Analysis on microphysical characteristics of three blizzard processes in Nanjing in the Winter of 2018 [J].Chin J Atmos Sci，43（5）：1095-1108.（in Chinese）.

刘善峰，陆正奇，韩永翔，等，2019.2008年初贵州电线积冰厚度的模拟研究：基于天气研究和预报（WRF） 模式耦合电线积冰预报系统[J].科学技术与工程，19（8）：303-309. Liu S F，Lu Z Q，Han Y X，et al.，2019.Simulation on ice thickness of transmission line in Guizhou Province in early 2008—based on weather research and forecasting （WRF） model couple with ice acceleration forecast system[J].Sci Technol Eng，19（8）：303-309.doi：10.3969/j.issn.1671-1815.2019.08.047.（in Chinese）.

刘雪静，牛生杰，2016.两次高压电线积冰过程气象成因分析[J].气象科学，36（2）：230-235. Liu X J，Niu S J，2016.Analysis of the causes for two high-voltage wire icing cases[J].J Meteor Sci，36（2）：230-235.（in Chinese）.

吕晶晶，牛生杰，周悦，等，2013.冬季高速公路路桥温度变化规律及能量平衡分析[J].大气科学学报，36（5）：546-553. Lü J J，Niu S J，Zhou Y，et al.，2013.Characteristics of bridge and road surface temperature changes in winter and energy budget analysis[J].Trans Atmos Sci，36（5）：546-553.doi：10.13878/j.cnki.dqkxxb.2013.05.001.（in Chinese）.

罗宁，文继芬，赵彩，等，2008.导线积冰的云雾特征观测研究[J].应用气象学报，19（1）：91-95. Luo N，Wen J F，Zhao C，et al.，2008.Observation study on properties of cloud and fog in ice accretion areas[J].J Appl Meteor Sci，19（1）：91-95.doi：10.3969/j.issn.1001-7313.2008.01.012.（in Chinese）.

Makkonen L，1989.Estimation of wet snow accretion on structures[J].Cold Reg Sci Technol，17（1）：83-88.doi：10.1016/S0165-232X（89）80018-7.

Makkonen L，2000.Models for the growth of rime，glaze，icicles and wet snow on structures[J].Philos Trans Royal Soc Lond Ser A：Math Phys Eng Sci，358（1776）：2913-2939.doi：10.1098/rsta.2000.0690.

Makkonen L，Wichura B，2010.Simulating wet snow loads on power line cables by a simple model[J].Cold Reg Sci Technol，61（2/3）：73-81.doi：10.1016/j.coldregions.2010.01.008.

Musilek P，Arnold D，Lozowski E P，2009.An ice accretion forecasting system （IAFS） for power transmission lines using numerical weather prediction[J].SOLA，5：25-28.doi：10.2151/sola.2009-007.

牛生杰，李蕊，呂晶晶，等，2011.三种下垫面温度及结冰预报模型研究[J].地球物理学报，54（4）：909-917. Niu S J，Li R，Lü J J，et al.，2011.Research on a numerical model for predicting three types of underlying surface temperature and ice[J].Chin J Geophys，54（4）：909-917.doi：10.3969/j.issn.0001-5733.2011.04.004.（in Chinese）.

Niu S J，Zhou Y，Jia R，et al.，2012.The microphysics of ice accretion on wires：observations and simulations[J].Sci China Earth Sci，55（3）：428-437.doi：10.1007/s11430-011-4325-8.

Poots G，Skelton P L I，1995.Simulation of wet-snow accretion by axial growth on a transmission line conductor[J].Appl Math Model，19（9）：514-518.doi：10.1016/0307-904X（95）00012-9.

Pytlak P，Musilek P，Lozowski E，et al.，2010.Evolutionary optimization of an ice accretion forecasting system[J].Mon Wea Rev，138（7）：2913-2929.doi：10.1175/2010mwr3130.1.

秦健，孙立军，2005.国外沥青路面温度预估方法综述[J].中外公路，25（6）：19-23. Qin J，Sun L J，2005.A summary of foreign asphalt pavement temperature prediction methods [J].J China Foreign Highw，25（6）：19-23.doi：10.3969/j.issn.1671-2579.2005.06.006.（in Chinese）.

Roberts E，Stewart R E，2008.On the occurrence of freezing rain and ice pellets over the eastern Canadian Arctic[J].Atmos Res，89（1/2）：93-109.doi：10.1016/j.atmosres.2007.11.032.

Ryerson C C，Kenyon P，1998.Atmospheric ice ablation processes on Mt Equinox，Vermont，USA[J].Atmos Res，46（1/2）：75-86.doi：10.1016/S0169-8095（97）00052-5.

Savadjiev K，Farzaneh M，2004.Modeling of icing and ice shedding on overhead power lines based on statistical analysis of meteorological data[J].IEEE Trans Power Deliv，19（2）：715-721.doi：10.1109/TPWRD.2003.822527.

史月琴，邓雪娇，胡志晋，等，2006.一次山地浓雾的三维数值研究[J].热带气象学报，22（4）：351-359. Shi Y Q，Deng X J，Hu Z J，et al.，2006.Three-demensional numerical study on dense fog over mountain area[J].J Trop Meteor，22（4）：351-359.（in Chinese）.

舒斯，熊守权，陈英英，等，2019.湖北省高速公路道路结冰预警模型[J].气象，45（11）：1589-1599. Shu S，Xiong S Q，Chen Y Y，et al.，2019.Warning model of road surface icing of expressway in Hubei Province[J].Meteor Mon，45（11）：1589-1599.doi：10.7519/j.issn.1000-0526.2019.11.009.（in Chinese）.

Szilder K，2019.Snow accretion prediction on an inclined cable[J].Cold Reg Sci Technol，157：224-234.doi：10.1016/j.coldregions.2018.11.001.

谭冠日，1982.电线积冰若干小气候特征的探讨[J].气象学报，40（1）：13-23. Tan G R，1982.On the microclimatological characteristics of icing on wires[J].Acta Meteorol Sin，40（1）：13-23.（in Chinese）.

陶玥，史月琴，刘卫国，2012.2008年1月南方一次冰凍天气中冻雨区的层结和云物理特征[J].大气科学，36（3）：507-522. Tao Y，Shi Y Q，Liu W G，2012.Characteristics of stratification structure and cloud physics of the freezing rain over Southern China in January 2008[J].Chin J Atmos Sci，36（3）：507-522.doi：10.3878/j.issn.1006-9895.2011.11082.（in Chinese）.

陶玥，李宏宇，刘卫国，2013.南方不同类型冰冻天气的大气层结和云物理特征研究[J].高原气象，32（2）：2501-2518. Tao Y，Li H Y，Liu W G，2013.Characteristics of atmospheric stratification and cloud physics of different types of freezing rain over Southern China[J].Plateau Meteor，32（2）：2501-2518.（in Chinese）.

滕中林，1978.架空线路的结冰及冰厚计算[J].电力技术，11（3）：53-56. Teng Z L，1978.Calculation of ice thickness on overhead lines[J].Electric Power，3：53-56.（in Chinese）.

Thériault J M，Stewart R E，Milbrandt J A，et al.，2006.On the simulation of winter precipitation types[J].J Geophys Res：Atmos，111（D18）：D18202.doi：10.1029/2005JD006665.

田华，张楠，张楷翊，等，2019.国外公路交通气象研究进展[J].气象与环境学报，35（3）：79-86. Tian H，Zhang N，Zhang K Y，et al.，2019.Advance on highway traffic meteorological research in foreign country[J].J Meteor Environ，35（3）：79-86.doi：10.3969/j.issn.1673-503X.2019.03.010.（in Chinese）.

Wang T S，Niu S J，Lü J，et al.，2019.Observational study on the supercooled fog droplet spectrum distribution and icing accumulation mechanism in Lushan，southeast China[J].Adv Atmos Sci，36（1）：29-40.doi：10.1007/s00376-018-8017-6.

王颖，王晓云，江志红，等，2012.中国电线积冰灾害危险性评估和区划[J].自然灾害学报，21（3）：110-118. Wang Y，Wang X Y，Jiang Z H，et al.，2012.Assessment and zonation of wire icing disaster risk in China[J].J Nat Disasters，21（3）：110-118.doi：10.13577/j.jnd.2012.0316.（in Chinese）.

王遵娅，2011.中国冰冻日数的气候及变化特征分析[J].大气科学，35（3）：411-421. Wang Z Y，2011.Climatic characters and changes of ice-freezing days in China[J].Chin J Atmos Sci，35（3）：411-421.（in Chinese）.

温华洋，田红，唐为安，等，2011.安徽省电线积冰标准冰厚的气象估算模型[J].应用气象学报，22（6）：747-752. Wen H Y，Tian H，Tang W A，et al.，2011.Establishment of meteorological model for estimating standard ice thickness in Anhui Province[J].J Appl Meteor Sci，22（6）：747-752.（in Chinese）.

吴兑，赵博，邓雪娇，等，2007.南岭山地高速公路雾区恶劣能见度研究[J].高原气象，26（3）：649-654. Wu D，Zhao B，Deng X J，et al.，2007.A study on bad visibility over foggy section of freeway in nanling mountainous region[J].Plateau Meteor，26（3）：649-654.doi：10.3321/j.issn：1000-0534.2007.03.030.（in Chinese）.

徐戈，孙继明，牛生杰，等，2016.冻滴微物理过程的分档数值模拟试验研究[J].大气科学，40（6）：1297-1319. Xu G，Sun J M，Niu S J，et al.，2016.A numerical study for the microphysical processes of ice pellets with a spectral（Bin） cloud model[J].Chin J Atmos Sci，40（6）：1297-1319.（in Chinese）.

杨贵名，孔期，毛冬艳，等，2008.2008年初“低温雨雪冰冻”灾害天气的持续性原因分析[J].气象学报，66（5）：836-849. Yang G M，Kong Q，Mao D Y，et al.，2008.Analysis of the long-lasting cryogenic freezing rain and snow weather in the beginning of 2008[J].Acta Meteorol Sin，66（5）：836-849.doi：10.3321/j.issn：0577-6619.2008.05.016.（in Chinese）.

曾明剑，陆维松，梁信忠，等，2008.2008年初中国南方持续性冰冻雨雪灾害形成的温度场结构分析[J].气象学报，66（6）：1043-1052. Zeng M J，Lu W S，Liang X Z，et al.，2008.Analysis of temperature structure for persistent disasterous freezing rain and snow over Southern China in early 2008[J].Acta Meteorol Sin，66（6）：1043-1052.（in Chinese）.

Zerr R J，1997.Freezing rain：an observational and theoretical study[J].J Appl Meteor，36（12）：1647-1661.doi：10.1175/1520-0450（1997）036<1647：fraoat>2.0.co;2.

中央气象局，2003.地面气象观测规范[M].北京：气象出版社：89-90. China Meteorological Administration，2003.Specifications for Surface Meteorological Observation[M].Beijing：China Meteorological Press：89-90.（in Chinese）.

Zhou Y，Niu S J，Lü J，et al.，2012.Meteorological conditions of ice accretion based on real-time observation of high voltage transmission line[J].Chin Sci Bull，57（7）：812-818.doi：10.1007/s11434-011-4868-2.

Zhou Y，Niu S J，Lü J，et al.，2016.The effect of freezing drizzle，sleet and snow on microphysical characteristics of supercooled fog during the icing process in a mountainous area[J].Atmosphere，7（11）：143.doi：10.3390/atmos7110143.

New advances in research on power line icing and pavement temperature

NIU Shengjie1，2，WANG Tianshu3，L Jingjing2，4，ZHOU Yue5，6，WANG Yuan4

1College of Safety Science and Engineering，Nanjing Tech University，Nanjing 211816，China;

2Key Laboratory for Aerosol-Cloud-Precipitation of China Meteorological Administration，Nanjing 210044，China;

3China Meteorological Administration Training Centre，Beijing 100081，China;

4School of Atmospheric Physics，Nanjing University of Information Science& Technology，Nanjing 210044，China;

5Wuhan Regional Climate Center，Wuhan 430074，China;

6Hubei Key Laboratory for Heavy Rain Monitoring and Warning Research，China Meteorological Administration，Institute of Heavy Rain，Wuhan 430205，China

From January 10th to February 2nd，2008，four rare low-temperature rain，snow and freezing weather events occurred consecutively in southern China.This disaster was characterized by a wide range，high intensity and long duration，thereby resulting in extremely serious economic losses and casualties.In addition，freezing rain，sleet，supercooled fog and other factors caused serious typical icing and secondary disasters，which in turn led to the lodging of electricity transmission and communication line towers，and interruption of power supply and communication and traffic transportation.Therefore，since 2009，the research team of this paper have carried out field experiments on wire icing and road icing，the results of which have progressed the research of wire icing from the traditional observation of icing meteorological conditions to a new highlight of comprehensive observation research regarding icing meteorological conditions and cloud and precipitation micro physics in icing weather events.The present paper reveals the micro physical mechanism of icing accretion，studies the icing accretion growth rate and its influencing factors，and establishes the icing accretion growth model，so as to better simulate the icing accretion growth process.Furthermore，observations on the temperature of the three underlying surfaces of asphalt，cement and soil and the temperature changes of different subgrade depths of expressways and bridges are carried out，and explanations are made from the perspective of energy balance.Based on the achievements of the research team，this paper studies and organizes the relevant representative research results，including the spatial and temporal distribution of ice accretion frequency，microphysical characteristics of icing weather，growth rate of ice and its influencing factors，meteorological conditions of icing，construction of ice growth model，and law of road icing growth and prediction，and proposes suggestions for the future research of wire icing.

power line icing;microphysical characteristics;icing meteorological factor;icing growth model;pavement temperature;research progress

doi：10.13878/j.cnki.dqkxxb.20210309001

（責任编辑：刘菲）