基于回歸模型的青藏鐵路水害氣象風險評估

張 靜，保廣裕*，周 丹，馬守存，劉 瑋

（青海省氣象服務中心，青海 西寧 810001）

鐵路是國民經濟的大動脈，其運輸安全對國民經濟發展具有重大意義。青藏鐵路是中國21世紀以來最重要的鐵路建設工程，是內陸地區聯系西藏的唯一鐵路通道[1-3]。青藏鐵路從青海省民和縣開始，途經湟水河谷地、青海湖流域、柴達木盆地等高原內陸盆地，翻越藏北高原以及昆侖山、唐古拉山、念青唐古拉山等高山，最終到達西藏自治區首府拉薩，全程約2052 km，是世界上海拔最高、線路最長、穿越多年凍土里程最長的高原鐵路[4-5]。其中青海境內青藏鐵路長度為1507 km，占青藏鐵路總長度的73%以上。近年，隨著中國高速鐵路網的快速布局，新建的蘭新高鐵在青海境內的里程也達到268 km。此外，青海境內哈木支線、敦格鐵路、茶卡支線、寧大支線的青海境內鐵路總里程達到2441 km。因此，如何面對復雜自然地理環境和安全有效的生產運營，成為人們關注的熱點。

青海省位于青藏高原東北部邊緣地帶，地質環境復雜，是現今地球表層構造性活動最強的地理單元[6-7]。同時，青海境內的大部分鐵路通過的都是高山峽谷，兩側地勢陡峭，巖石破碎，路段土質松散，具備誘發滑坡和泥石流發育的條件[8]。加之高原天氣變化多端，氣候差異明顯，尤其是夏季高原中、小尺度天氣系統活動頻繁，降水發生突然、歷時短且強度大，極易造成滑坡和泥石流等水害，對鐵路路基、橋梁、隧道造成很大的威脅。我國鐵路每年都遭受不同程度的水害侵襲，國內學者做了很多關于降水特征與水害的關系、降水致災的發生頻數與斷道時間之間的關系、地質災害時空分布特征與氣象條件的關系及鐵路地質災害與氣象預警模型[10-15]等方面的研究，但在具備誘發滑坡和泥石流發育條件的青海境內，鐵路水害氣象風險評估方面的研究較少。因此開展青藏鐵路青海境內降水引發的水害氣象風險評估研究，有助于提升鐵路部門應該對氣象災害防御的能力，滿足對鐵路部門的服務需求，提高鐵路專業氣象的服務能力。

1 資料數據和分析方法

1.1 資料數據

1.1.1 水害資料

鐵路水害是指發生在鐵路沿線的由于高強度或持續性降水（包括因降水誘發的滑坡、泥石流、崩坍等地質災害）造成鐵路受損或鐵路運輸中斷的現象[9-10]。如果一天內因水害造成鐵路同一路段范圍內多次受損，統計時按一次計算。選取2000—2014年青藏鐵路的西寧工務段、德令哈工務段和格爾木工務段提供的鐵路沿線水害統計資料，具有水害出現的時間和地點資料總計176次，其中因降水導致的鐵路水害為158次，占整個鐵路水害的89.8%。

1.1.2 氣象數據

選取2000—2014年青藏鐵路沿線青海境內民和、樂都、平安、西寧、湟源、海晏、剛察、天峻（天棚）、茶卡、烏蘭、德令哈、格爾木、五道梁、沱沱河14個氣象站點逐日降水資料，選取2012—2014年青藏鐵路沿線青海境內28個加密氣象站降水資料進行參考分析。

1.2 分析方法

1.2.1 相關分析及其T檢驗

用來描述兩個時間序列之間相互關系的方法，主要用相關系數Rkl來表示，Rkl的絕對值越大表示兩者之間關系越密切，一般用T檢驗來判斷其相關的可信程度。假設兩個時間序列Xkl、Xli，其樣本長度均為n，其中必須有一個是一維時間序列，那么相關系數表達式為：

對于氣候變量不同時刻的線性相關和兩氣候變量的線性相關是否顯著，即相關系數達到多少算是存在顯著相關關系，必須進行統計檢驗。由相關系數r，依據下式反算t值：

顯然 t-t（N-1）分布，給出信度 α，查表可以得到臨界值 tα，比較t和 tα：當t大于tα，則說明通過信度為α的顯著性檢驗，反之則未通過顯著性檢驗。

1.2.2 邏輯回歸方法

對于鐵路水害預測模型而言，所要回歸的隨機變量為水害發生的概率Pi，Pi取值范圍（0，1），水害發生與否不是一個連續的變量（假設水害發生時Pi=1，不發生時Pi=0），故普通的多元線性回歸方法將不適用于推導此類自變量和因變量之間的關系。而邏輯回歸[16]可以解決此類問題，只需引進連接函數，將Pi取值范圍映射到（-∞，+∞），則邏輯回歸函數可表示為：

式中xi，j（j=1，…，m）為m個預報因子對應于第i個觀測樣本的災害發生概率，Pj（j=0，…，m）為需要擬合的模型參數，模型觀測樣本為觸發或未觸發水害的雨量數據。

由（3）式可得：

利用（4）式便可預測鐵路水害發生概率，其難點是引入邏輯模型中雨量因子的篩選。

1.2.3 逐步回歸方法

在建立多元回歸方程的過程中，按偏相關系數的大小次序將自變量逐個引入方程，對引入方程中的每個自變量偏相關系數進行統計檢驗，效應顯著的自變量留在回歸方程內，循此繼續遴選下一個自變量。如果效應不顯著，停止引入新自變量。由于新自變量的引入，原已引入方程中的自變量由于變量之間的相互作用其效應有可能變得不顯著者，經統計檢驗確證后要隨時從方程中剔除，只保留效應顯著的自變量。直至不再引入和剔除自變量為止，從而得到最優的回歸方程[17]。

2 青海境內鐵路水害特征分析

2.1 鐵路水害時間分布特征

青藏鐵路沿線2000—2014年水害發生次數年際變化曲線（圖1），由圖可以看出，2004、2008、2012年水害發生次數較多，其中在2004年水害發生次數最多，發生了31次，而對應5—9月平均降水量來看，并非完全對應，但也存在一定關系，原因是水害發生的頻次不僅與降水量級有關，還與降水的性質和持續時間有關。據統計1980—2000年中國鐵路氣象因子導致的水害共發生1376次，而其余因子導致的發生水害339次[10]，由此可見氣象因子是鐵路水害發生的最主要的致災因子。

圖1 青藏鐵路沿線2000—2014年水害發生次數與5—9月降水量年際變化曲線

圖2 是2000—2014年青藏鐵路沿線青海境內鐵路水害各月分布圖。由圖可以看出，鐵路水害發生在月分布上存在著顯著差異，1—3月和10—12月沒有鐵路水害發生，4月開始發生鐵路水害，5—6月呈顯著上升趨勢，7—8月鐵路水害發生次數達到最高值，9月開始急劇下降。鐵路沿線水害主要集中在7—8月，其中7月和8月各發生70次，這兩月水害發生次數占總發生次數的79.5%。其次是5、6月和9月，水害發生次數分別為12、15次和8次，水害發生次數占總發生次數的19.9%，4月只發生過1次鐵路水害，出現在2013年。可見，鐵路沿線水害主要發生在夏季的7—8月，這與青海省境內降水量的季節變化相一致。

圖2 青藏鐵路沿線2000—2014年水害發生次數月分布

2.2 鐵路水害空間分布特征

2000—2014年青藏鐵路沿線青海境內水害發生次數空間分布（圖3），可以看出，青藏鐵路沿線青海境內水害的空間分布各地差異較大。鐵路水害發生頻率最高的路段是樂都—平安、湟源—海晏，15 a中總計發生鐵路水害分別為46、47次；其次是天棚—烏蘭、烏蘭—德令哈、格爾木—玉珠峰，15 a中總計發生鐵路水害分別為31、32次和31次。民和—樂都、平安—西寧、西寧—湟源、海晏—剛察、天棚—茶卡、德令哈—飲馬峽、玉珠峰—五道梁、五道梁—沱沱河段15 a間鐵路水害發生次數均在10～30次之間，西寧—大通、飲馬峽—格爾木段15a間鐵路水害發生次數均不足10次。由此可見，青藏鐵路沿線青海境內鐵路水害主要集中在西寧工務段的樂都—平安、湟源—海晏，德令哈工務段的天棚—烏蘭、烏蘭—德令哈，格爾木工務段的格爾木—玉珠峰，這與青海省強降水導致洪澇災害及其所引發的滑坡、泥石流等地質災害的空間分布[18-19]規律是基本相符的。

圖3 青藏鐵路沿線2000—2014年水害發生次數空間分布

2.3 鐵路水害類型及影響

據調查鐵路沿線水害的主要類型為路基工程水害、涵洞工程水害、橋梁工程水害和水漫線路水害。具體包括山洪暴發、泥石流淤埋線路、沖空道床、路塹坍塌、路基沉陷沖毀、水漫線路、泥石流沖擊橋墩、道床懸空、淤積物掩埋線路、山體滑坡路塹邊坡滑塌、水毀便道等。由此可知，鐵路水害的主要誘發因子就是降水及其次生災害。

3 降水與鐵路水害的關系

3.1 降水類型與鐵路水害的關系

根據統計分析表明，導致鐵路沿線水害的主要降水類型為區域性強降水和連陰雨。在158次鐵路水害中，由區域性強降水引起的鐵路水害130次，由連陰雨造成的鐵路水害28次。區域性強降水造成的水害占總次數的82.2%，其原因可能是由于過程雨量強度大、時間短，易引發次生災害。

根據降水量等級標準，統計2000—2014年青藏鐵路沿線逐日降水量，沿線各站各量級降水日數與水害發生次數（表1），可見，在降水引發的158次鐵路水害中，136次鐵路水害當日降水量達到了中到大雨及以上級別（≥17.0 mm），占鐵路水害總次數的86.1%，其中中到大雨級別（17.0～24.9 mm）導致鐵路水害發生48次，大雨級別（25.0～49.9 mm）導致鐵路水害發生84次，暴雨及其以上級別（≥50.0mm）導致鐵路水害4次。通過比較不同量級降水與鐵路水害關系表明：中到大雨、大雨、大到暴雨與鐵路水害的相關性顯著，其中大雨與鐵路水害的相關系數達到0.69。由于暴雨及其以上降水在青海境內是小概率事件，樣本數少，因而導致水害的次數也較少。

表1 青藏鐵路沿線各站各量級降水日數與水害發生次數

3.2 降水持續時間與鐵路水害的關系

通過對青藏鐵路沿線強降水過程分析發現，水害的發生發展與降水天氣過程的持續時間關系密切。經統計，青海境內的一次強降水天氣過程的維持時間一般小于5 d。因此，根據鐵路沿線各地水害易發程度，選取1～5 d相應影響系數、水害發生等級條件，以24 h為間隔5 d逐日降水量，計算鐵路水害誘發指數[20]，由此來表征強降水持續時間對水害的影響。由于某日的降雨量對成災的影響隨時間的延長而逐漸減弱，因此本文取當日影響系數為1.0，以距離災害發生日天數為分母逐日遞減，建立青海境內鐵路水害誘發指數方程[15]如下：

式中：Y表示水害誘發因素大小，R0～R4表示水害潛伏因素大小，單位為mm，R0為水害發生前24 h雨量，R1為前25～48 h雨量，R2為前49～72 h雨量，R3為前73～96 h小時雨量，R4為前97～120 h雨量。1.0、1/2、1/3、1/4、1/5為經驗系數，表示該降雨日有效降水量權重大小，R0、R1、R2、R3、R4為以北京時間20：00為日界的1～5 d降水量。

3.3 前期降水量與鐵路水害的關系

鐵路水害發生時，當日降水量和前期的降水量都對水害發生有一定影響，但如何將這些因子在聚類分析中聚為一類，在此引入有效降水量[21]，將水害發生當日和前4 d逐日降水量結合起來，作為5 d有效降水量因子，具體計算公式如下[21]：

式中，R為5 d有效降水量，R0為水害發生前24 h雨量，R1為前25～48 h雨量，R2為前49～72 h雨量，R3為前73～96 h雨量，R4為前97～120 h雨量。此經驗公式說明，距水害發生時間愈近，日降水量對水害發生的作用就愈大，并將前期降水量和當天降水量對鐵路水害的影響有效統一，因此選用5 d的有效降水量作為鐵路水害的致災的因子之一。

4 水害氣象風險評估模型建立

首先應用邏輯回歸模型判斷鐵路沿線各站水害是否有發生的可能，若無發生可能，則定義鐵路水害等級為4級，發生鐵路水害概率低，不需要采取防范措施；若水害有發生的可能，則應用逐步回歸模型計算鐵路沿線各站水害等級（1～3級，其中1級為發生鐵路水害概率很高、封鎖警戒，2級為發生鐵路水害概率較高、限速警戒，3級為發生鐵路水害概率一般、出巡警戒），并綜合考慮水害致災臨界值研判指標判斷是否出警）。等級劃分參考《鐵路運輸安全保護條例》、《地質災害防治條例》、《鐵路實施〈中華人民共和國防汛條例〉細則》、《鐵路技術管理規程》，以及《鐵路防洪標準化管理辦法》。

4.1 鐵路水害評估邏輯回歸模型

邏輯回歸評估模型是對水害發生與否的判定，是利用各降雨因子與災害發生的關系，選擇相關性較大的雨量因子，納入概率預測模型[22]，對鐵路水害發生與否進行預測。在研究中假定同一區段災點的地質條件一致，分析降雨因子對水害的影響[22]。考慮到短時強降水和持續性降雨造成水害的特點不同，前者是即發型，后者具有一定的滯后效應，此時誘發鐵路水害的雨量需要一定的時間累積，故模型選取與當日降水（R0）、前1天降水（R1）、前2天降水（R2）、前3天降水（R3）、前4天降水（R4）、5 d的降水合計（R5）、降水誘發指數（R6）、有效降水（R7）共計8個因子，運用相關性分析和邏輯回歸方法篩選對災害發生貢獻率較大的降水量因子，分區段建立鐵路沿線水害概率預評估模型。運用SPSS軟件進行邏輯回歸分析，得到各區段相關統計量。

湟源—海晏段鐵路水害發生頻率是整個青藏鐵路沿線最高的，因此本文以湟源—海晏段為例計算其回歸方程相關統計量。由表2可以看出，對于湟源-海晏段水害的發生與當日降水（R0）、前1天降水（R1）、有效降水（R7）關系密切。方差分析結果表明，上述路段的因子通過0.05的顯著性檢驗。

表2 湟源—海晏段邏輯回歸方程相關統計量

把自變量雨量因子代入模型，對邏輯回歸評估模型進行預測效果評估，預測評估效果見表3（0表示無水害，1表示有水害）。當降雨因子組合為R0、R1和R7時，湟源—海晏段發生鐵路水害的判對率是89.6%，模型總判對率為90.4%，評估模型的概率預測結果相對較好。

表3 湟源—海晏段邏輯回歸評估模型預測效果

所以，通過對8個降雨因子進行邏輯回歸，得出湟源—海晏段的水害預評估模型為：

類似上述方法，分別對其他區段進行邏輯回歸分析，得到各區段的預測模型的總判對率（表4）所示，其中格爾木—五道梁段預測模型總判對率最高為100%，五道梁—沱沱河、天棚—茶卡、樂都—平安的預測模型總判對率均不到90%。

表4 青藏鐵路沿線各段邏輯回歸評估模型預測效果

4.2 鐵路水害評估逐步回歸模型

逐步回歸評估模型是在邏輯回歸判定有水害發生的基礎上建立的，模型計算鐵路沿線各站水害等級（1～3級），其中1級為封鎖警戒，2級為限速警戒，3級為出巡警戒。

4.2.1 水害與氣象因子的相關性分析

根據青藏鐵路沿線各站點的資料統計，利用公式（1）對沿線各站當日降水量、前4 d逐日的降水量、過程累計降水量、水害誘發指數、有效降水量8個因子與鐵路水害發生時相應的等級進行相關性分析，得出：鐵路沿線水害等級與當日降水量（R0）、過程累計降水量（R5）、水害誘發指數（R6）和有效降水量（R7）的相關性較高，相關系數均通過了0.001（閾值是0.380）的顯著性檢驗。相關系數見表5。

4.2.2 逐步回歸評估模型建立

目前回歸分析方法按照涉及自變量的多少，分為一般回歸和逐步回歸分析，本研究涉及自變量較多，因此采用逐步回歸分析方法建立青藏鐵路沿線各段水害等級評估模型，模型方程如下：

式中：Y為鐵路水害等級，P、P1、P2、P3、P4分別代表相關因子通過逐步回歸方法計算的各因子系數，R0表示當日降水量，R5表示過程累計降水量，R6表示水害誘發指數，R7表示有效降水量。青藏鐵路沿線各區段相關因子逐步回歸系數如表6所示。

表5 鐵路沿線各氣象因子與水害等級相關分析結果

表6 青藏鐵路沿線各段逐步回歸評估模型參數值

4.2.3 逐步回歸評估模型檢驗

利用2015年6—9月實況水害資料對鐵路水害逐步回歸評估模型預報為1級的情況進行檢驗（表7）。總計預報1級封鎖警戒9次，實況出現了8次，預報準確率為88.9%，另外一次預報為1級，而實況為2級限速警戒。由此可見，逐步回歸評估模型在預報和評估青藏鐵路青海境內鐵路水害工作中具有較高的預報準確率和應用價值。

表7 2015年6—9月青藏鐵路沿線水害預報1級實況檢驗

5 結論

本文以青藏鐵路青海段民和到沱沱河為研究對象，基于回歸模型開展青藏鐵路水害氣象風險評估研究，主要得到以下結論：

（1）鐵路水害的發生與降水因子密不可分，并且由降水引發的次生災害是鐵路水害的主要誘發因子。水害發生的頻次不僅與降水量級有關，還與降水的性質和持續時間有關。

（2）青藏鐵路沿線青海境內水害主要集中在7—8月，這兩月水害發生次數占總發生次數的79.5%，這與青海省境內降水量的季節變化相一致；鐵路沿線水害的空間分布各地差異較大，發生頻率最高的路段是樂都—平安、湟源—海晏。

（3）青藏鐵路青海境內水害的主要降水類型為區域性強降水和連陰雨，中到大雨、大雨、大到暴雨與鐵路水害的相關性顯著。

（4）根據降水對鐵路水害的各種影響關系，結合沿線各地水害易發程度，參考鐵路部門四級警戒雨量標準，設定降雨誘發鐵路水害氣象預警等級分別為4、3、2、1級，對應發生鐵路水害概率低、一般、較高和很高。

（5）采用邏輯回歸模型分析得到的青藏鐵路青海境內各區段是否發生水害，預測模型總判對率均超過86.2%，對水害發生頻次最高路段（湟源—海晏段）的判對率是89.6%，評估模型的概率預測效果較好。

（6）在判定水害發生的基礎上，應用逐步回歸評估模型，計算得出鐵路沿線各站水害等級，在實際檢驗鐵路水害9次預報中，準確預報了8次，預報準確率高達88.9%。

參考文獻：

[1] 李衛東，張玥.青藏鐵路運營對西藏經濟發展的影響[J].中國鐵路，2016，55（5）：12-16.

[2] 張瑞英，席建超，姚予龍，等.基于視覺廊道的青藏鐵路沿線旅游動態景觀評價[J].生態學報，2014，34（12）：3320-3330.

[3] 鄭度，張鐿鋰.青藏鐵路沿線生態與環境安全[M].杭州：浙江科技出版社，2009.

[4] 尹國安，牛富俊，林戰舉，等.青藏鐵路沿線多年凍土分布特征及其對環境變化的響應[J].冰川凍土，2014，36（4）：772-781.

[5] 蔡漢成，李勇，楊永鵬，等.青藏鐵路沿線多年凍土區氣溫和多年凍土變化特征[J].巖土力學與工程學報，2016，35（7）：1434-1444.

[6] 李軍成.青藏鐵路及沿線生態地質環境研究[D].北京：中國地質大學，2014.

[7] 劉鴻源.青藏鐵路地質災害主要類型及危險性評價研究[D].成都：西南交通大學，2016.

[8] 劉峰貴，張海峰，陳瓊，等.青藏鐵路沿線自然災害地理組合特征分析[J].地理科學，2010，30（3）：384-390.

[9] 接小峰，徐林榮，曹祿來，等.鐵路水害預警雨量閾值預測方法及其應用[J].自然災害學報，2015，24（5）：177-183.

[10] 魏慶朝，楊永平，張大煒.中國鐵路水害環境致災因子分析[J].自然災害學報，2002，11（1）：123-127.

[11] 王冀，馬寧，申文軍，等.華北地區汛期降水特征及對鐵路水害的影響[J].自然災害學報，2016，25（4）：30-39.

[12] 王勇，李曉霞，李曉蘋.蘭州鐵路防洪指揮氣象預警服務系統[J].干旱氣象，2009，27（4）：415-418.

[13] 朱良峰，吳信才，劉修國.基于GIS的鐵路地質災害信息管理與預警預報系統[J].山地學報，2004，22（2）：230-253.

[14] 陳穎，陳鵬翔，江遠安，等.烏魯木齊河流域致災洪水臨界雨量分析[J].沙漠與綠洲氣象，2017，11（2）：8-13.

[15] 徐羹慧，毛煒嶧，陸幗英.新疆氣象災害近期變化和防災減災工作綜合評述[J].沙漠與綠洲氣象，2008，2（1）：50-54.

[16] 李倩，趙中英.基于邏輯回歸的信息轉發預測模型研究[J].軟件導刊，2017（2）：4-6.

[17] 保廣裕，張景華，錢有海，等.青藏鐵路沿線風沙指數預報方法研究[J].青海環境，2010（1）：44-46.

[18] 費杜秋，劉峰貴，周強，等.青藏鐵路沿線滑坡泥石流災害風險分析[J].干旱區地理，2016，39（2）：345-352.

[19] 巨克英，海顯蓮，韓有才.青海省暴雨洪澇氣象災害特征分析及其預報[J].青海氣象，2011（3）：7-9.

[20] 周秀美，苗蕓，程林.紅河州地質災害臨界雨量及預報方法初探[J].云南地理環境研究，2012，24（3）：37-42.

[21] 王禮先，王志民.山洪及泥石流災害預報[M].北京：中國林業出版社，2001.

[22] 張清，黃朝迎.我國鐵路水害及其評估模型研究[J].應用氣象學報，1999（4）：498-502.