陳 凱，徐國慶，熊 燾，侯紀(jì)勇，王曦晨，戚國輝，劉 青

（1.國家電網(wǎng)公司交流建設(shè)分公司，北京 100052；2.紫光軟件系統(tǒng)有限公司，北京 100084）

0 引言

我國北方冬季積雪和融水易在路面形成結(jié)冰，造成路面抗滑能力顯著降低，車輛制動(dòng)距離延長，影響行車安全[1]。尤其我國東北部地區(qū)，冬春季節(jié)氣溫低，常出現(xiàn)降雪、凍雨[2]，道路結(jié)冰現(xiàn)象突出。研究道路結(jié)冰預(yù)報(bào)方法，對(duì)道路結(jié)冰情況進(jìn)行提前預(yù)報(bào)，能夠有效減少因道路結(jié)冰引起的交通事故。

美國及歐洲一些國家向來重視道路結(jié)冰的監(jiān)測和預(yù)報(bào)工作，并較早地在公路交通氣象學(xué)領(lǐng)域開展研究，研發(fā)出了一些路面溫度測試和結(jié)冰預(yù)報(bào)技術(shù)[3-4]。如美國等國家建立的道路氣象信息系統(tǒng)（Road Meteorological Information System），可根據(jù)路域氣象資料預(yù)報(bào)路面狀況。Sass[5]基于地面熱傳導(dǎo)方程和地表能量平衡方程，進(jìn)行了敏感性實(shí)驗(yàn)，獲取了詳細(xì)的溫度和濕度分析結(jié)果，并利用丹麥道路站的數(shù)據(jù)對(duì)道路溫度和結(jié)冰預(yù)測進(jìn)行了數(shù)值模擬，建立了道路結(jié)冰預(yù)報(bào)方法，取得了較好的效果。Shao等[6]通過因特網(wǎng)收集公路養(yǎng)護(hù)信息而建立的道路結(jié)冰預(yù)報(bào)模型能夠自動(dòng)調(diào)整所需物理參數(shù)，準(zhǔn)確及時(shí)地預(yù)測路面結(jié)冰時(shí)間，為歐洲許多國家冬季公路養(yǎng)護(hù)提供了技術(shù)支持。

目前，我國在道路結(jié)冰預(yù)報(bào)研究方面仍處于初期階段[7-8]。楊武杰等[9]利用1990—2007年內(nèi)蒙古高等級(jí)公路氣象站監(jiān)測的結(jié)冰和積雪天氣數(shù)據(jù)及1990—2000年降水天氣數(shù)據(jù)，分析了引起道路結(jié)冰的氣象和環(huán)流因素，構(gòu)建了道路結(jié)冰統(tǒng)計(jì)模型，對(duì)包頭至興和高等級(jí)公路道路結(jié)冰情況進(jìn)行了預(yù)報(bào)。孫翠梅等[10]利用2003年12月至2010年3月鎮(zhèn)江市公路地面氣溫?cái)?shù)據(jù)，分析了鎮(zhèn)江市道路結(jié)冰的氣候特征，結(jié)合降水性質(zhì)、是否有積雪、日最低氣溫和日平均氣溫，定義了道路結(jié)冰影響指數(shù)，建立了統(tǒng)計(jì)模型進(jìn)行分級(jí)預(yù)報(bào)，對(duì)當(dāng)?shù)氐缆方Y(jié)冰情況具備一定的預(yù)報(bào)能力。馬吉暉[11]從天氣學(xué)角度分析總結(jié)出3種天氣概念模型，得出由弱降水所致的道路結(jié)冰預(yù)報(bào)方法。雷建軍等[12]將遺傳微粒群算法與支持向量機(jī)（Support Vector Machines，簡稱SVM）相結(jié)合，利用遺傳微粒群算法的全局搜索能力從多個(gè)氣象因子中自適應(yīng)地選擇合適的因子組合及其對(duì)應(yīng)的最優(yōu)SVM參數(shù)，分別建立了武漢市和十堰市道路結(jié)冰預(yù)報(bào)模型，取得了較好的預(yù)報(bào)效果。

總體來看，目前有關(guān)道路結(jié)冰預(yù)報(bào)的研究方法主要有三種：統(tǒng)計(jì)模型[13-14]、理論分析模型[15]和數(shù)據(jù)挖掘模型[16-17]。統(tǒng)計(jì)模型在進(jìn)行道路結(jié)冰預(yù)報(bào)時(shí)，數(shù)據(jù)量大，雖然具有一定的預(yù)報(bào)能力，但空報(bào)率較高。理論分析模型是基于熱力學(xué)和天氣學(xué)基本原理，利用數(shù)值方法建立溫度、降水的預(yù)報(bào)模型，但是該方法參數(shù)復(fù)雜，獲取困難，不適合大型主設(shè)備運(yùn)輸期的道路結(jié)冰災(zāi)害預(yù)報(bào)。數(shù)據(jù)挖掘法可以處理大規(guī)模數(shù)據(jù)，對(duì)應(yīng)的模型在屬性選擇合理的情況下準(zhǔn)確率相對(duì)較高，并且在分析非線性氣象問題上有著廣泛的應(yīng)用。決策樹算法是數(shù)據(jù)挖掘方法的一種，能夠克服數(shù)據(jù)要求高、參數(shù)復(fù)雜、計(jì)算繁瑣等困難[18-19]。

本文將基于C4.5決策樹算法，通過對(duì)道路結(jié)冰樣本進(jìn)行規(guī)則學(xué)習(xí)，訓(xùn)練決策樹模型，將道路是否結(jié)冰抽象為一個(gè)二元分類問題，構(gòu)建道路結(jié)冰預(yù)報(bào)模型，最后利用特高壓輸變電主設(shè)備運(yùn)輸期的道路結(jié)冰預(yù)報(bào)對(duì)模型的準(zhǔn)確率加以驗(yàn)證。

1 C4.5決策樹算法原理

決策樹算法屬于數(shù)據(jù)挖掘中一種重要的分類方法[20-21]，目前在遙感解譯領(lǐng)域得到了廣泛應(yīng)用。本文結(jié)合遙感解譯思想，引入C4.5決策樹算法，通過對(duì)日最低氣溫、日積雪深度、日降水量與對(duì)應(yīng)當(dāng)天道路結(jié)冰關(guān)系進(jìn)行學(xué)習(xí)訓(xùn)練，構(gòu)建道路結(jié)冰二元分類模型，為道路結(jié)冰預(yù)報(bào)提供理論算法基礎(chǔ)。

在決策樹算法中，令S為包括s類樣本的數(shù)據(jù)訓(xùn)練集，m為訓(xùn)練集中的類別數(shù)量，si為s中第i類樣本的個(gè)數(shù)，i=1,2…,m，即：

此時(shí)訓(xùn)練集S的信息（嫡）H(S)為：

根據(jù)屬性A，令S={S1,S2,…,Sv}，其中v為屬性A包含不同值的數(shù)目。于是，訓(xùn)練集在屬性A上的信息熵為：

對(duì)應(yīng)的信息增益為：

2 模型構(gòu)建

2.1 道路結(jié)冰因素分析

道路結(jié)冰是指雨、雪、凍雨或霧滴碰到溫度低于0℃的地面而出現(xiàn)的積雪或結(jié)冰現(xiàn)象，包括凍結(jié)的殘雪、凸凹的冰轍、雪融水或其他原因的道路積水在寒冷冬季形成的結(jié)冰[22]。我國東北部地區(qū)冬季的主要?dú)夂蛱攸c(diǎn)是風(fēng)大、干旱、寒冷、冰期長。顯然，冬季的低溫環(huán)境給道路結(jié)冰災(zāi)害的形成創(chuàng)造了條件。本文以通遼地區(qū)為例，統(tǒng)計(jì)了2015年和2016年各自2月份和3月份的氣溫、降水、降雪與道路結(jié)冰情況，并分析得出這三因素與道路結(jié)冰之間的關(guān)系，在時(shí)間上與主設(shè)備運(yùn)輸時(shí)間保持一致，氣象條件類似，有利于道路結(jié)冰預(yù)報(bào)研究。圖1～圖4中縱軸分別表示最低氣溫、降雨量、積雪深度和是否結(jié)冰（1表示結(jié)冰，0表示不結(jié)冰），橫軸表示日期。

圖1 2015年2月份最低氣溫、降雨量、積雪深度與道路結(jié)冰關(guān)系

圖2 2015年3月份最低氣溫、降雨量、積雪深度與道路結(jié)冰關(guān)系

圖3 2016年2月份最低氣溫、降雨量、積雪深度與道路結(jié)冰關(guān)系

圖4 2016年3月份最低氣溫、降雨量、積雪深度與道路結(jié)冰關(guān)系

由圖1～圖4可知，通遼地區(qū)在2015年2月、3月以及2016年2月的最低氣溫和3月前半月的最低氣溫均低于0℃，3月后半月最低氣溫波動(dòng)較大，有時(shí)高于0℃；對(duì)應(yīng)月份的降雨量較少，每月不超過2次，且降雨當(dāng)天氣溫低于0℃，出現(xiàn)道路結(jié)冰；降雪雖然頻次不高，由于氣溫較低，雪融較慢，單次降雪會(huì)造成連續(xù)幾天積雪，長期積雪則會(huì)導(dǎo)致道路長期連續(xù)結(jié)冰。根據(jù)數(shù)據(jù)統(tǒng)計(jì)，道路結(jié)冰出現(xiàn)在有積雪或者降雨并且氣溫低于0℃的情況下，由于內(nèi)蒙和遼寧冬季和初春降水較少，故引起道路結(jié)冰的主要因素是低溫和降雪。

2.2 基于C4.5決策樹算法的道路結(jié)冰預(yù)報(bào)模型構(gòu)建

研究表明，地表溫度是道路結(jié)冰的主要影響因素[23]。由于地表溫度相對(duì)最低氣溫較難測得，不過日最低氣溫和日最低地表溫度間存在較好的對(duì)應(yīng)關(guān)系[24]，因此可用日最低氣溫替代日最低地表溫度來研究道路結(jié)冰問題。由前文可知，道路結(jié)冰天氣出現(xiàn)在有積雪或降雨的條件下，因此，可利用滿足日最低氣溫與日積雪深度或日最低氣溫與降水的關(guān)系來判斷道路是否結(jié)冰，即將通遼地區(qū)2015—2016年2月和3月日最低氣溫、日積雪深度、日降水量作為決策樹模型的輸入變量，道路是否結(jié)冰作為目標(biāo)變量（結(jié)冰用1表示，不結(jié)冰用0表示）來訓(xùn)練并構(gòu)建道路結(jié)冰預(yù)報(bào)模型。

基于C4.5決策樹算法的道路結(jié)冰預(yù)報(bào)模型的構(gòu)建具體可分以下三個(gè)步驟，對(duì)應(yīng)的模型運(yùn)行原理如圖5和圖6所示。

（1）收集日最低氣溫、日降水量、日積雪深度以及道路是否結(jié)冰等數(shù)據(jù)；

（2）對(duì)最低氣溫、日降水量、日積雪深度以及對(duì)應(yīng)的道路結(jié)冰情況進(jìn)行分類預(yù)處理；

（3）將日最低氣溫、日積雪深度、日降水量作為決策樹模型的輸入變量，將道路是否結(jié)冰作為目標(biāo)變量，構(gòu)建道路結(jié)冰預(yù)報(bào)模型。

圖5 日最低氣溫與日積雪深度道路結(jié)冰預(yù)報(bào)模型示意圖

圖6 日最低氣溫與日降水量道路結(jié)冰預(yù)報(bào)模型示意圖

3 基于C4.5決策樹算法的道路結(jié)冰預(yù)報(bào)模型應(yīng)用及精度分析

為了驗(yàn)證基于C4.5決策樹算法的道路結(jié)冰預(yù)報(bào)模型的準(zhǔn)確率，現(xiàn)分別對(duì)點(diǎn)狀區(qū)域和線狀區(qū)域兩種類型的道路結(jié)冰情況進(jìn)行預(yù)報(bào)，并與實(shí)際情況作一比較。

3.1 點(diǎn)狀區(qū)域道路結(jié)冰預(yù)報(bào)模型應(yīng)用及精度分析

以赤峰地區(qū)為例，將其2016年2月和3月的道路實(shí)際結(jié)冰數(shù)據(jù)（每日最低氣溫、每日早上8時(shí)積雪深度、每日累積1d降水量）輸入預(yù)報(bào)模型進(jìn)行結(jié)冰情況預(yù)報(bào)。預(yù)報(bào)結(jié)果和實(shí)際結(jié)果如表1和表2所示。

表1 2016年2月赤峰道路模型預(yù)報(bào)結(jié)冰與實(shí)際結(jié)冰情況

表2 2016年3月赤峰道路模型預(yù)報(bào)結(jié)冰與實(shí)際結(jié)冰情況

由表1和表2可知，2月共29d，準(zhǔn)確預(yù)報(bào)有28d，預(yù)報(bào)準(zhǔn)確率為96.6%，3月共31d，準(zhǔn)確預(yù)報(bào)有28d，預(yù)報(bào)準(zhǔn)確率為90.3%，2月和3月共60d，準(zhǔn)確預(yù)報(bào)有56d，總體預(yù)報(bào)準(zhǔn)確率為93.3%。分析可知，2月9日與3月13日模型計(jì)算結(jié)果與實(shí)際結(jié)果不符，這都是因?yàn)楫?dāng)日積雪量較少，氣溫較高，積雪融化成水，路面未能結(jié)冰；3月17日有少量降雨，最低氣溫為0℃，而道路結(jié)冰模型訓(xùn)練數(shù)據(jù)在降雨條件下的最低氣溫為-2℃，故判斷為不結(jié)冰，而道路結(jié)冰是由綜合因素引起的，導(dǎo)致預(yù)報(bào)錯(cuò)誤。3月27日預(yù)報(bào)無結(jié)冰，是因?yàn)?6日中午氣溫較高，積雪融化成水，27日無降水，輸入的日累計(jì)降雨量為0，實(shí)際路面有積水且當(dāng)天最低氣溫較低，引起道路結(jié)冰，造成預(yù)報(bào)結(jié)果與實(shí)際不符?？傮w來看，盡管模型對(duì)部分日期道路結(jié)冰判斷錯(cuò)誤，但模型預(yù)報(bào)的整體準(zhǔn)確率仍較高。

3.2 線狀區(qū)域道路結(jié)冰預(yù)報(bào)模型應(yīng)用及精度分析

現(xiàn)以2017年冬季特高壓交流輸變電主設(shè)備運(yùn)輸?shù)缆方Y(jié)冰預(yù)報(bào)為例，對(duì)所建模在線狀區(qū)域預(yù)報(bào)中的精度進(jìn)行驗(yàn)證。

特高壓交流輸變電主設(shè)備在2017年的運(yùn)輸時(shí)間分別為2月12日—22日和3月20日—30日，運(yùn)輸線路從沈陽到錫林浩特，沿線經(jīng)過阜新、通遼和赤峰。為了保障運(yùn)輸安全，會(huì)提前1d預(yù)報(bào)次日沿線運(yùn)輸?shù)缆返慕Y(jié)冰情況。將本次運(yùn)輸過程中日最低氣溫、日積雪深度和日降水量輸入道路結(jié)冰預(yù)報(bào)模型，對(duì)道路是否結(jié)冰進(jìn)行預(yù)報(bào)，結(jié)果見表3和表4。

表3 2017年2月12日—22日設(shè)備運(yùn)輸沿線道路預(yù)報(bào)結(jié)冰與實(shí)際結(jié)冰情況

表4 2017年3月20日—30日設(shè)備運(yùn)輸沿線道路預(yù)報(bào)結(jié)冰與實(shí)際結(jié)冰情況

由表3和表4所示的預(yù)報(bào)結(jié)冰和實(shí)際結(jié)冰結(jié)果對(duì)比可知，特高壓交流輸變電主設(shè)備運(yùn)輸時(shí)間共計(jì)22d，準(zhǔn)確預(yù)報(bào)19d，錯(cuò)誤預(yù)報(bào)3d，預(yù)報(bào)準(zhǔn)確率為86.4%?？梢?，該道路結(jié)冰預(yù)報(bào)模型具有較高的準(zhǔn)確率，當(dāng)然日最低氣溫、日積雪深度、日降水量3項(xiàng)輸入數(shù)據(jù)較高的準(zhǔn)確度也為道路結(jié)冰預(yù)報(bào)的準(zhǔn)確性提供了保障。錯(cuò)誤預(yù)報(bào)原因分析如下：模型預(yù)報(bào)2月16日道路有結(jié)冰，而實(shí)際當(dāng)天無結(jié)冰，是因?yàn)樘鞖忸A(yù)報(bào)當(dāng)天有降雪，而真正的降雪出現(xiàn)在晚上（設(shè)備運(yùn)輸暫停），因此預(yù)報(bào)結(jié)果與實(shí)際情況不符；3月24日預(yù)報(bào)道路無結(jié)冰而實(shí)際道路有結(jié)冰，是因?yàn)樵O(shè)備運(yùn)輸經(jīng)過國道時(shí)該地區(qū)23日有少量降水，24日氣溫較低且路面有積水引起了道路結(jié)冰；3月29日預(yù)報(bào)道路結(jié)冰而實(shí)際當(dāng)天無結(jié)冰，是因?yàn)楫?dāng)日氣溫較高，且伴有大風(fēng)，在設(shè)備運(yùn)輸開始前路面降水已經(jīng)消失，沒有引起結(jié)冰，導(dǎo)致預(yù)報(bào)結(jié)果與實(shí)際不符。綜上可見，本文所建道路結(jié)冰預(yù)報(bào)模型在線狀區(qū)域道路結(jié)冰預(yù)報(bào)中可靠性高，能滿足道路結(jié)冰預(yù)報(bào)需求。

4 結(jié)語

本文采用C4.5決策樹算法，利用日最低氣溫、日降水量、日積雪深度與對(duì)應(yīng)的道路結(jié)冰之間的關(guān)系，構(gòu)建了基于C4.5決策樹算法的道路結(jié)冰預(yù)報(bào)模型。最后，通過對(duì)赤峰地區(qū)和沈陽到錫林浩特特高壓交流輸變電主設(shè)備運(yùn)輸期間道路結(jié)冰情況的預(yù)報(bào)對(duì)所建模型進(jìn)行實(shí)用性和精度驗(yàn)證，結(jié)果表明該道路結(jié)冰預(yù)報(bào)模型具有較高的準(zhǔn)確率，可為我國東北地區(qū)冬季道路運(yùn)輸提供有效、安全的技術(shù)保障。利用數(shù)據(jù)挖掘方法，通過分析歷史數(shù)據(jù)可發(fā)現(xiàn)潛在的信息和規(guī)律，在局部區(qū)域可以實(shí)現(xiàn)對(duì)道路結(jié)冰的有效預(yù)報(bào)。今后若結(jié)合與氣象相關(guān)的理論知識(shí)，進(jìn)行更加深入的信息挖掘，探索本質(zhì)信息，將會(huì)對(duì)道路結(jié)冰有更加科學(xué)理性的認(rèn)識(shí)和預(yù)報(bào)。

參考文獻(xiàn)

[1]王學(xué)強(qiáng)，王澄海，孟雪峰，等.內(nèi)蒙古錫林郭勒盟2012年冬季暴雪過程天氣學(xué)特征研究[J].冰川凍土，2013，35（6）：1446-1453.

[2]孫婷，張昌利，孟穎，等.面向公路主動(dòng)式防冰處置的行業(yè)信息物理系統(tǒng)[J].交通運(yùn)輸研究，2017，3（3）：49-54.

[3]NORMAN J.Slipperiness on Roads-An Expert Systerm Classification[J].Meteorological Applications,2000,7(1):27-36.

[4]SHAO J,LISTER P J.An Automated Nowcasting Model of Road Surface Temperature and State for Winter Road Main?tenance[J].Journal of Applied Meteorology,1996,35(8):1352-1361.

[5]SASS B H.A Numarical Model for Prediction of Road Tem?prature and Ice[J].Journal of Applied Meteorology,2007,91(4):475-499.

[6]SHAO J,LISTER P.Automated Road Ice Nowcasts and In?ternet[C]//Proceedingofthe11thInternationalRoadWeather Conference.Sapporo:Standing International Road Weather Commission,,2002.

[7]李蕊，牛生杰，汪玲玲，等.三種下墊面溫度對(duì)比觀測及結(jié)冰氣象條件分析[J].氣象，2011，37（3）：325-333.

[8]張朝林，張利娜，程叢蘭，等.高速公路氣象預(yù)報(bào)系統(tǒng)研究現(xiàn)狀與未來趨勢[J].熱帶氣象學(xué)報(bào)，2007，23（6）：652-658.

[9]楊武杰，司瑤冰，胡英華.內(nèi)蒙古高等級(jí)公路道路結(jié)冰預(yù)報(bào)方法初探[J].內(nèi)蒙古氣象，2010（5）：45-47.

[10]孫翠梅，馬俊峰，孔啟亮，等.鎮(zhèn)江市道路結(jié)冰影響指數(shù)分級(jí)預(yù)報(bào)方法[J].氣象與環(huán)境學(xué)報(bào)，2013，29（4）：84-88.

[11]馬吉暉.弱降水造成通化市道路結(jié)冰的天氣學(xué)預(yù)報(bào)方法[J].吉林氣象，2013（3）：17-20.

[12]雷建軍，韋惠紅，李劍.一種遺傳微粒群算法-支持向量機(jī)集成道路結(jié)冰預(yù)測系統(tǒng)[J].華中師范大學(xué)學(xué)報(bào)（自然科學(xué)版），2010，44（3）：392-396.

[13]黃立葵，賈璐，萬劍平，等.瀝青路面溫度狀況的統(tǒng)計(jì)分析[J].中南公路工程，2005，30（3）：8-9.

[14]宋麗萍，韓基良，劉宇飛，等.哈爾濱市冬季地面溫度變化特征及預(yù)報(bào)模型[J].黑龍江氣象，2014（4）：28-30.

[15]曹廣真，侯鵬，毛顯強(qiáng).北京市城市化對(duì)地表溫度時(shí)空特征的影響[J].氣象，2010，36（3）：19-26.

[16]劉梅，尹東屏，王清樓，等.南京地區(qū)冬季路面結(jié)冰天氣標(biāo)準(zhǔn)及其預(yù)測[J].氣象科學(xué)，2007（6）：685-690.

[17]熊竹，秦杰君.基于模擬實(shí)驗(yàn)的路面結(jié)冰預(yù)測研究[J].公路工程，2014，39（1）：216-220.

[18]ZHANG W,GAO S,CHEN B,et al.The Application of De?cision Tree to Intensity Change Classification of Tropical Cyclones in Western North Pacific[J].Geophsical Research Letters,2013,40(9):1883-1887.

[19]劉忠陽，杜子璇，陳懷亮，等.基于決策樹的鄭州市土地利用與覆蓋分類研究[J].氣象與環(huán)境科學(xué)，2008，31（3）：48-53.

[20]QUINLAN J R.C4.5 Programs for Machine Learning[M].San Francisco:Morgan Kaufmann Publishers Inc.,1993.

[21]史達(dá)偉，耿煥同，吉辰，等.基于C4.5決策樹算法的道路結(jié)冰預(yù)報(bào)模型構(gòu)建及應(yīng)用[J].氣象科學(xué)，2015，35（2）：204-209.

[22]保廣裕，王文英，張景華，等.青海省道路結(jié)冰變化時(shí)空分布特征及其影響等級(jí)劃分[J].氣象科技，2016，44（1）：104-110.

[23]蔣建瑩，史歷，倪允琪.一次“高影響天氣”的弱降雪過程的數(shù)值研究[J].應(yīng)用氣象學(xué)報(bào)，2005，16（2）：231-237.

[24]楊曉玲，蘭曉波，丁文魁，等.河西走廊東部地面溫度的變化特征及影響因子分析[J].干旱地區(qū)農(nóng)業(yè)研究，2014，32（1）：223-229.