基于MNL模型的車車碰撞事故嚴重程度影響因素辨識方法

王健宇，陸化普，孫智源，王天實

(1. 清華大學 土木工程學院，北京 100084；2. 北京工業大學 城市交通學院,北京 100124)

0 引言

交通安全問題是全世界共同面對的嚴峻公共安全挑戰。根據世界衛生組織(World Health Organization，WHO)最新全球交通事故統計報告：2016年世界范圍內因道路交通事故死亡的總人數超過135萬[1]，道路交通事故位列全球十大致死原因的第8名，是全球十大致死原因中唯一的非疾病類影響因素[2]。日益增長的道路交通量導致了交通事故的不斷增加[3]，因此對道路交通安全的研究一直是交通領域緊迫的基礎研究課題[4]。

道路交通事故研究的主要目的是降低事故的發生風險和減少事故的危害性[5]，事故的頻次分析旨在降低事故的發生風險[6]，事故的嚴重程度分析旨在減少事故的危害性[7]。本研究聚焦交通事故嚴重程度的研究。

事故嚴重程度的分析方法主要包含離散選擇模型和數據挖掘方法。離散選擇模型主要用于找尋事故嚴重程度的顯著影響因素，是穩定且使用久遠的統計分析模型。離散選擇模型主要包含logistic模型和probit模型[8]。Al-Ghamdi采用二項logistic回歸模型對9項相關因素(事故位置、事故形態、碰撞方位、事故時間、事故原因、駕駛員國籍、駕照狀態、車輛形態、駕駛員年齡)與事故是否致死之間的關系進行分析，最終發現事故發生地點和事故原因是具有最顯著相關性的[9]。Celik和Oktay采用多項logistic回歸模型對土耳其埃爾祖魯姆省和卡爾斯省發生的11 771例道路交通事故的影響因素進行分析，發現駕駛員年齡達到65歲及以上、駕齡低、高速路、人行橫道是可能導致死亡事故高發的因素[10]。Azimi等(2020)采用有序logistic模型對大型貨車傾覆事故的影響因素進行分析，發現粗心駕駛行為和異常駕駛情況是可能導致嚴重事故高發的因素[11]。參考已有研究，發現logistic模型能夠有效地找尋事故嚴重程度的影響因素。因此，本研究采用Multinomial Logistic Regression(MNL)模型，即多項logistic回歸模型，找尋不同嚴重程度下的顯著影響因素。

交通事故嚴重程度的研究對象主要為交通事故嚴重程度的影響因素[12]。事故嚴重程度的影響因素包含人員屬性[13]、車輛屬性[14]、道路及環境屬性[15]、時間屬性[16]等。人員屬性的影響因素包括：駕駛員年齡[17]、駕齡[18]、性別[19]，是否涉及行人[20]等。車輛屬性的影響因素包括：碰撞類型[21]、事故責任[22]等。道路及環境屬性的影響因素包括：城市分區[23]、道路物理條件[24]、天氣[25]等。時間屬性的影響因素包括：季節[26]、高峰[27]等。參考已有研究中包含的影響因素，本研究整理形成4大類的影響因素進行研究，分別為過失方人員及車輛因素、受害方人員及車輛因素、時間因素、環境因素。

綜上所述，本研究采用MNL模型，對車輛間的交通事故進行分析，辨識影響事故嚴重程度的主要因素。

1 影響因素的定義與數據描述統計

根據沈陽市2015年至2017年的交通事故數據，對兩車之間的不同類型事故進行研究，提取對象為碰撞雙方具有差異化的責任，即同時包含過失方(全責或主責)和受害方(無責或次責)，共提取事故593例。

對593例事故進行標定和結構化處理，將事故嚴重程度定義為因變量Y，將事故嚴重程度的影響因素定義為自變量Xk。影響因素Xk共包含19項，可分為4大類：過失方人員及車輛因素、受害方人員及車輛因素、時間因素、環境因素。其中每大類包含的具體因素如下所示：

(1)過失方人員及車輛因素：車輛類型(X1)、性別(X2)、年齡(X3)、駕齡(X4)、逃逸行為(X5)、碰撞類型(X6)；

(2)受害方人員及車輛因素：車輛類型(X7)、性別(X8)、年齡(X9)、駕齡(X10)；

(3)時間因素：工作日(X11)、季節(X12)、時段(X13)、高峰時刻(X14)；

(4)環境因素：城市分區(X15)、道路橫斷面位置(X16)、道路物理隔離(X17)、路面狀況(X18)、天氣(X19)。

表1所示的結果為車輛間事故嚴重程度和影響因素不同取值的定義。其中，根據貨運、客運車輛的差異性將車輛類型分為小客車、大客車、小貨車、大貨車和其他等5類，對年齡和駕齡通過合理分組進行分析。

表1 車輛間事故嚴重程度和影響因素不同取值的定義Tab.1 Definitions of different values of accident severity and influencing factors between vehicles

2 多項logistic回歸模型

本研究以沈陽交通事故為數據來源進行分析，令Y為事故嚴重程度。依照財產損失、受傷人數和死亡人數的參考值，事故的嚴重程度被設置為3分類變量。令Y的取值為1，2，3，分別對應僅財產損失、有受傷但無死亡、有死亡這3種事故嚴重程度的情況，數據中僅財產損失(Y=1)占比12.14%、有受傷但無死亡(Y=2)占比47.55%、有死亡(Y=3)占比40.30%。Xk為K個影響因素的第k項，bjk為影響因素Xk的影響系數，則對j=1，2，3的無分類變量，則可構建2個分類變量的模型，其模型表達式為：

(1)

(2)

式中，a1,a2為常數項；α為模型參數，即Logistic回歸系數。

對全部J分類中的第j項概率計算表達式為：

(3)

通過計算，可獲得bjk的值，即為各項影響因素的系數，其指數函數Exp(bjk)即為各影響因素的OR值，也稱優勢比，表示相對于被參考項，該影響因素在此取值情況下可能導致Y發生變化的相對概率。通常顯著水平P的臨界值設為0.05來篩選顯著的影響因素。對建立的模型需進行似然比檢驗，當P<0.05時說明模型具有意義。

3 虛擬變量的建立

對影響因素Xk需對其中的多分類變量進行處理，引入虛擬變量，使多分類變量轉化多個為0-1分類的變量，代入進行計算。通常引入虛擬變量的規則為：對某含有n分類的變量因素Xk，引入n-1個虛擬變量，使得該變量避免多重共線性。以過失方車輛類型(X1)為例建立虛擬變量，共5種取值，分別為1=小客車(低速汽車、微型客車、小型客車)；2=小貨車(輕型貨車)；3=大客車(中型客車、大型客車)；4=大貨車(汽車列車、中型貨車、重型貨車)；5=其他。其中，以1=小客車為參考，引入4個虛擬變量，如表2所示。

表2 過失方車輛類型(X1)引入虛擬變量后的結果Tab.2 result of introducing dummy variable into vehicle type of fault party (X1)

依照上述方法對過失方年齡(X3)、過失方駕齡(X4)、碰撞類型(X6)、受害方車輛類型(X7)、受害方年齡(X9)、受害方駕齡(X10)、季節(X12)、時段(X13)、城市分區(X15)、道路橫斷面位置(X16)、道路物理隔離(X17)進行數據處理，引入虛擬變量代入計算。

4 MNL模型結果分析

MNL模型在加入二元自變量參量和經虛擬變量處理后的自變量參量后，擬合優化，最終模型擬合似然比檢驗P<0.05，說明構建模型具有統計學意義，模型分析結果具有可靠性。

令P=0.05為顯著性界限值，在全部19項影響因素中，經模型計算可得過失方車輛類型(X1)、過失方年齡(X3)、碰撞類型(X6)、受害方車輛類型(X7)、城市分區(X15)這5項影響因素為顯著的影響因素，其對應P值分別為0.004，0.021，0.003，0.000，0.000，Y=1與Y=2情況下的計算結果如表3和表4所示。其中，B表示影響因素的對應參數；P為顯著值；Wald，即參數B與其標準誤平方值的比值，表征對參數B的檢驗值；OR值為參數B的指數函數，即優勢比，表征相對的發生概率倍數。本研究對Y=1和Y=2情況下的OR值進行分析，以Y=3為參照，表征不同場景發生的相對概率。

4.1 僅財產損失事故的計算結果

Y=1表示事故嚴重程度為“僅財產損失”，其顯著值P小于0.05處對應自變量的系數有顯著統計意義，即對因變量有顯著性的影響。OR值表示Y=1“僅財產損失”場景在該參考項下對Y=3“有死亡”場景發生概率相對的可能性倍數。根據表3所示的結果，正面碰撞(X6=1)、受害方車輛類型為小客車(X7=1)、主城區(X15=1)、近郊區(X15=2)為具有顯著意義的參考項，過失方車輛類型(X1)、過失方年齡(X3)的各類取值為無顯著意義的參考項。

表3 MNL模型下Y=1的參數估計值Tab.3 Estimated parameters in MNL model when Y=1

(1)碰撞類型(X6)

對于碰撞類型(X6)，其取值為1時為具有顯著影響的參考項取值(P=0.033)，表明“僅財產損失”事故與“有死亡”事故相比，當碰撞事故類型為正面碰撞(X6=1)時，更偏向于出現“僅財產損失”的事故，這種可能性為其他碰撞類型(X6=4)情況的4.854倍。因此對于不同碰撞類型而言，正面碰撞更多出現的是輕微事故。

(2)受害方車輛類型(X7)

對于受害方車輛類型(X7)，其取值為1時為具有顯著影響的參考項取值(P=0.001)，表明“僅財產損失”事故與“有死亡”事故相比，當受害方車輛類型為小客車(X7=1)時，更偏向出現“僅有財產損失”的事故，這種可能性為其他車輛類型(X7=4)情況的9.986倍。因此，對于不同類型的車輛而言，小客車發生事故類型多為僅有財產損失的輕微事故，相比于小客車，其他類型的車輛更應該關注其產生嚴重交通事故的可能性。

(3)城市分區(X15)

對于城市分區(X15)，其取值為1和2都具有顯著影響的參考項取值(P=0.003和P=0.012)，表明“僅財產損失”事故與“有死亡”事故相比，當城市分區分別為主城區(X15=1)和近郊區(X15=2)時，更偏向于出現“僅財產損失”事故，這種可能性分別為市域高速(X15=4)情況的35.292倍、18.549倍。因此，相比于在市域高速出現的交通事故，出現在主城區和近郊區的交通事故都主要為僅有財產損失的輕微事故。王建軍等通過統計分析方法發現高速公路發生的事故具有重大、特大事故占比高的特點，部分原因可能為追尾碰撞、氣候情況、駕駛員特征及車輛本身特征，需相應找尋預防事故和提高保護措施的方法[28]。對于沈陽市，主城區和近郊區的經濟發展水平高于遠郊區，基礎設施防護水平相對較高，因而更傾向于出現輕微事故。

(4)對整體顯著但對僅財產損失事故不顯著

過失方車輛類型(X1)和過失方年齡(X3)這兩項影響因素對模型整體為顯著的影響因素，但在Y=1情況與Y=3情況的比對分析中未展現其顯著性。因而對“僅財產損失”事故和“有死亡”事故相比較的場景，過失方車輛類型(X1)和過失方年齡(X3)并未展現其顯著性。

4.2 有受傷但無死亡事故的計算結果

Y=2表示事故嚴重程度為“有受傷但無死亡”，其顯著值P小于0.05處對應自變量的系數有顯著統計意義，即對因變量有顯著性的影響。OR值表示Y=2“有受傷但無死亡”場景在該參考項下對Y=3“有死亡”場景發生概率相對的可能性倍數。根據表4的結果，過失方年齡為19至45歲(X3=2)、側面碰撞(X6=2)為具有顯著意義的參考項。過失方車輛類型(X1)、受害方車輛類型(X7)、城市分區(X15)的各類取值為無顯著意義的參考項。

表4 MNL模型下Y=2的參數估計值Tab.4 Estimated parameters in MNL model when Y=2

(1)過失方年齡(X3)

對于過失方年齡(X3)，其取值為2時為具有顯著影響的參考項取值(P=0.008)，表明“有受傷但無死亡”事故與“有死亡”事故相比，當過失方年齡為19～45歲(X3=2)時，更偏向出現“有受傷但無死亡”的事故，這種可能性為60歲及以上(X3=4)情況的6.671倍。因此，對于不同年齡的過失方，19～45歲的群體容易出現受傷事故而非死亡事故，一種可能的解釋為19～45歲群體精力很旺盛，因而可能在駕駛過程中風格依然較為激進，但由于其身體反應條件和注意力都更加集中，因此，當遇到危險時，通常這一年齡段的駕駛員都可以通過使用果斷的措施來避免致命事故的出現。在未來的監管改進中，可對這一年齡段群體進行持續研究，挖掘其出現更容易出現交通事故，但并非出現致死交通事故的根本原因。

(2)碰撞類型(X6)

對于碰撞類型(X6)，其取值為2時為具有顯著影響的參考項取值(P=0.018)，表明“有受傷但無死亡”事故與“有死亡”事故相比，當碰撞類型為側面碰撞(X6=2)時，更偏向于出現“有受傷但無死亡”的事故，這種可能性為其他碰撞類型(X6=4)的2.037倍。因此，說明參與交通組織的交通工具中對于側面碰撞風險的保護做得較好，可有效地降低死亡事故出現的概率，下一步提高方向為車輛提高側面碰撞的承受能力，進一步將“有死亡”事故轉化為“有受傷但無死亡”事故。通過側面碰撞的分析模型，評價側面碰撞對乘客胸部、髖關節等部位的損傷，發現當側面碰撞被限制在一定的低侵入速度時，通過車門自身剛度可顯著減輕乘客所受到的傷害[29]。同時，側面氣囊在側面碰撞時對乘客的保護效果也是至關重要的[30]。以上研究和本節的結論充分說明在汽車工業設計中，通過合理的工藝設計使得車輛在發生側面碰撞時，有效地降低了致死事件發生的可能性。

(3)對整體顯著但對有受傷但無死亡事故不顯著

過失方車輛類型(X1)、受害方車輛類型(X7)和城市分區(X15)這3個影響因素對模型整體為顯著的影響因素，但在Y=2情況與Y=3情況的比對分析中未展現其顯著性。因而對“有受傷但無死亡”事故和“有死亡”事故相比較的場景，過失方車輛類型(X1)、受害方車輛類型(X7)、城市分區(X15)并未展現其顯著性。

過失方車輛類型(X1)在Y=1和Y=2兩種情況與Y=3的比對分析中均未展現出顯著性，但對模型是顯著的影響因素，則可得其為弱影響因素，僅有部分取值貼近影響因素顯著性的臨界值。

5 結論

本研究以沈陽為例，進行了車輛間碰撞事故影響因素的研究。從過失方人員及車輛屬性、受害方人員及車輛屬性、時間屬性、道路及環境屬性等4個方面構建了影響因素集合。建立MNL模型識別道路交通安全的顯著性影響因素，并根據事故嚴重程度的不同等級進行了差異性討論，獲得研究結論如下：

(1)在降低車輛間事故的安全風險時，需要重點考慮過失方車輛類型、過失方年齡、碰撞類型、受害方車輛類型、城市分區等5個因素。

(2)對僅財產損失事故和有受傷但無死亡事故，影響因素的取值存在差異性。其中對僅財產損失事故，正面碰撞、受害方車輛類型為小客車,主城區、近郊區是可能的顯著影響因素。對有受傷但無死亡事故，過失方年齡為19～45歲,側面碰撞是可能的顯著影響因素。

(3)持續提高對駕駛員的教育警示，使用新技術監督駕駛員的行車狀態，提升車輛防御側面碰撞的能力是可能降低事故嚴重程度的手段。