陳 堅(jiān),邱智宣,彭 濤,劉柯良,傅志妍,庹永恒

(1. 重慶交通大學(xué) 交通運(yùn)輸學(xué)院,重慶 400074;2. 重慶第二師范學(xué)院 經(jīng)濟(jì)與工商管理學(xué)院,重慶 400067; 3. 重慶市公安局 大渡口區(qū)分局交巡警支隊(duì),重慶 400084)

隨著城鎮(zhèn)化的快速推進(jìn),城市機(jī)動(dòng)車保有量迅猛增長(zhǎng),交通安全問(wèn)題也日益突出。近年來(lái),我國(guó)道路交通死亡人數(shù)長(zhǎng)期居于世界首位,每年有6萬(wàn)人以上死于道路交通事故,直接財(cái)產(chǎn)損失達(dá)13億元,交通安全形勢(shì)依然嚴(yán)峻。因此,從多維視角解析交通事故十分必要。

國(guó)內(nèi)外學(xué)者圍繞交通安全進(jìn)行了系統(tǒng)研究。已有研究從事故構(gòu)成要素入手,將人、車、路及自然環(huán)境等方面與事故的關(guān)系進(jìn)行分析。董傲然等[1]從人、車、路和自然環(huán)境這4個(gè)方面出發(fā),構(gòu)建了機(jī)動(dòng)車與行人事故嚴(yán)重程度的分析模型;S.KUMAR等[2]考慮了受傷者的性別與年齡、照明條件、道路情況等11個(gè)因素,對(duì)印度車輛交通事故嚴(yán)重程度進(jìn)行了對(duì)比分析。也有學(xué)者從區(qū)域視角出發(fā),探索區(qū)域人口經(jīng)濟(jì)、土地利用特征及區(qū)域路網(wǎng)特征等因素與交通安全之間的關(guān)系。王雪松等[3]選擇了區(qū)域人口特征、道路及交通特性、路網(wǎng)形態(tài)特征等指標(biāo)對(duì)交通安全進(jìn)行分析,結(jié)果發(fā)現(xiàn)路網(wǎng)形態(tài)和交通安全顯著相關(guān);P.NAJAF等[4]基于年齡、人口、人均收入等測(cè)度,分析了城市形態(tài)特征與交通安全的關(guān)系;潘義勇等[5]考慮了公園、醫(yī)院等建成環(huán)境因素,探索了老年人因事故造成傷害嚴(yán)重程度的影響因素。

現(xiàn)有研究?jī)H考慮了部分建成環(huán)境因素對(duì)事故的影響,如道路特征[3]、人口密度[4]、土地利用情況[6]等。尚缺少對(duì)建成環(huán)境因素系統(tǒng)的考慮及與其它因素的對(duì)比,導(dǎo)致難以從城市規(guī)劃本源出發(fā),來(lái)規(guī)避交通事故的發(fā)生。筆者基于重慶市某區(qū)交通事故信息,運(yùn)用機(jī)器學(xué)習(xí)方法,探究建成環(huán)境對(duì)交通事故嚴(yán)重程度的非線性關(guān)系,得到了變量的重要度。

1 建成環(huán)境

1.1 研究單元?jiǎng)澐?/h3>

空間數(shù)據(jù)分析結(jié)果會(huì)隨著面積單元定義的不同而發(fā)生變化,其影響主要體現(xiàn)在尺度效應(yīng)與劃區(qū)效應(yīng)兩個(gè)方面,因此適宜的空間尺度對(duì)空間數(shù)據(jù)分析及結(jié)果獲取十分重要。楊文越等[7]對(duì)比分析了多尺度建成環(huán)境對(duì)居民通勤出行的狀況,發(fā)現(xiàn)不同尺度建成環(huán)境的顯著性為:社區(qū)尺度>社區(qū)邊界1 km緩沖范圍尺度>街道尺度。由于文中研究的對(duì)象為事故點(diǎn),社區(qū)尺度相對(duì)狹窄,故筆者選取事故點(diǎn)1 km緩沖區(qū)作為建成環(huán)境的研究單元。

1.2 建成環(huán)境指標(biāo)選取

與自然環(huán)境不同,建成環(huán)境是土地利用、城市設(shè)計(jì)及交通系統(tǒng)等一系列相關(guān)要素的集合。建成環(huán)境的定量測(cè)度一直是城市規(guī)劃學(xué)的研究熱點(diǎn),從20世紀(jì)80年代開始建成環(huán)境的內(nèi)涵不斷被豐富。R.EWING等[8]將建成環(huán)境歸結(jié)為“3D”維度,包括密度、多樣性和設(shè)計(jì),而后在“3D”維度基礎(chǔ)上增加了目的地可達(dá)性和到站點(diǎn)距離,形成“5D”維度[9]。

筆者所指的建成環(huán)境測(cè)度指標(biāo)是在“5D”維度基礎(chǔ)上,對(duì)建成環(huán)境要素進(jìn)行重構(gòu),甄選人口密度、土地利用混合度、路網(wǎng)密度、距市中心距離、公交站點(diǎn)密度等11個(gè)變量用于度量建成環(huán)境的密度、多樣性、交通設(shè)計(jì)、可達(dá)性及公共交通鄰近度5個(gè)維度。建成環(huán)境變量描述見表1,其中密度為點(diǎn)密度,通過(guò)計(jì)算POI點(diǎn)個(gè)數(shù)與研究區(qū)域面積之比得到;多樣性通過(guò)土地利用混合度表征[10],計(jì)算如式(1);交通設(shè)計(jì)涵蓋交叉口密度與路網(wǎng)密度;可達(dá)性通過(guò)計(jì)算各研究單元距離城市中心的歐幾里得距離得到;公共交通鄰近度用研究單元內(nèi)的公交站點(diǎn)密度及地鐵站點(diǎn)密度進(jìn)行度量。

(1)

式中:Li為第i個(gè)研究單元的土地利用熵指數(shù);Pij為第i個(gè)研究單元內(nèi)第j種POI數(shù)量占總數(shù)的比例;N為POI興趣點(diǎn)個(gè)數(shù)。

表1 建成環(huán)境變量

1.3 建成環(huán)境數(shù)據(jù)獲取與處理

建成環(huán)境數(shù)據(jù)處理主要借助Arc GIS平臺(tái),涉及數(shù)據(jù)包括:街道行政區(qū)數(shù)據(jù)、道路網(wǎng)絡(luò)數(shù)據(jù)和POI興趣點(diǎn)數(shù)據(jù)。街道行政區(qū)數(shù)據(jù)以面為單元,用于研究空間范圍的構(gòu)建;道路網(wǎng)絡(luò)數(shù)據(jù)包括不同等級(jí)的道路線網(wǎng),以測(cè)度研究區(qū)域的路網(wǎng)密度;POI興趣點(diǎn)數(shù)據(jù)通過(guò)Python中的爬蟲工具在Open Street Map中獲取,用于計(jì)算不同研究區(qū)域點(diǎn)密度及可達(dá)性。建成環(huán)境數(shù)據(jù)處理流程見圖1。

圖1 建成環(huán)境數(shù)據(jù)的處理過(guò)程

2 城市道路事故數(shù)據(jù)處理與分析

2.1 事故分類與統(tǒng)計(jì)

交通事故按事故嚴(yán)重程度分為輕微事故、一般事故、重大事故和特大事故4類,不同等級(jí)事故劃分標(biāo)準(zhǔn)見表2。

表2 事故等級(jí)分類

2.2 事故數(shù)據(jù)地理編碼

由于事故原始數(shù)據(jù)中只有對(duì)事故發(fā)生點(diǎn)的地址描述,沒(méi)有經(jīng)緯度數(shù)據(jù),因此需要對(duì)事故原始數(shù)據(jù)進(jìn)行地理編碼。將表示地址的文字在地圖中進(jìn)行相似性匹配,并在這些表示地址的文字中匹配到相似度最高的文字,最后返回匹配結(jié)果和相應(yīng)的經(jīng)緯度信息作為最終結(jié)果。田沁等[11]比較了百度、高德、搜狗和騰訊等地圖廠商的地理編碼應(yīng)用程序編程接口的服務(wù)品質(zhì),認(rèn)為騰訊的地理編碼API整體服務(wù)品質(zhì)較好,具有較高的數(shù)據(jù)質(zhì)量和較完備的地址數(shù)據(jù)。故筆者采用騰訊的地理編碼API進(jìn)行事故發(fā)生點(diǎn)地址到經(jīng)緯度坐標(biāo)轉(zhuǎn)換。

2.3 事故空間分布特征

筆者采用2016—2018年重慶市某區(qū)102 km2范圍內(nèi)的道路交通事故數(shù)據(jù),共計(jì)154起,經(jīng)剔除缺失值、地理編碼后,最終確定了126起事故數(shù)據(jù)作為研究對(duì)象?；趐ython對(duì)事故傷亡情況進(jìn)行識(shí)別并做出統(tǒng)計(jì),不同等級(jí)事故統(tǒng)計(jì)情況見表3。

表3 事故數(shù)據(jù)統(tǒng)計(jì)

將含經(jīng)緯度的事故數(shù)據(jù)導(dǎo)入GIS中,運(yùn)用核密度分析工具描述事故的空間分布特征。不同等級(jí)事故分布情況見圖2,圖2中:顏色越深的地方表示事故數(shù)量越多;從分布情況來(lái)看,西部的事故數(shù)量明顯高于東部,特別是靠近北側(cè)的事故數(shù)量最多,同時(shí)這3種類型事故均在不同區(qū)域出現(xiàn)了集聚現(xiàn)象。

圖2 事故分布

3 事故非線性分析模型

3.1 模型選取

傳統(tǒng)統(tǒng)計(jì)學(xué)模型由于變量的度量方式(如單位、數(shù)量級(jí)等)的不同,難以體現(xiàn)出各變量對(duì)因變量的重要程度。而梯度提升決策樹(gradient boosting decision tree, GBDT)模型具有以下優(yōu)勢(shì):① 沒(méi)有預(yù)先確定兩個(gè)變量之間的結(jié)構(gòu)關(guān)系假設(shè),可以實(shí)現(xiàn)數(shù)據(jù)的自我學(xué)習(xí);② 可以提供各變量在結(jié)果中相對(duì)重要性估計(jì);③ 與回歸模型不同,該模型對(duì)建成環(huán)境與交通問(wèn)題分析中經(jīng)常存在的多重共線性問(wèn)題不敏感。因此,筆者選擇GBDT模型對(duì)建成環(huán)境因素與事故嚴(yán)重程度的非線性關(guān)系進(jìn)行分析。

3.2 變量描述

城市道路交通事故是人、車、路、環(huán)境等諸多因素耦合作用的結(jié)果,因此除了考慮建成環(huán)境因素外,同時(shí)將個(gè)體行為、道路情況、事故類型和自然環(huán)境因素納入到模型中,其它變量見表4。

表4 其它變量描述與賦值

3.3 模型構(gòu)建與求解

GBDT是基于Boosting算法的集成決策樹模型,Boosting算法依據(jù)上一次訓(xùn)練的殘差生成基學(xué)習(xí)器。GBDT是在Boosting基礎(chǔ)上,沿著殘差減小的梯度方向上建立新的決策樹。GBDT方法的核心目標(biāo)是不斷減小訓(xùn)練過(guò)程產(chǎn)生的殘差,模型可由式(2)表示[12]。

(2)

(3)

式中:γj、c分別為待估計(jì)參數(shù);hj(x)為弱分類器j的估計(jì)結(jié)果;J為弱分類器數(shù)量。

GBDT模型求解過(guò)程如下[13]:

2)計(jì)算第m次迭代的殘差,如式(4):

(4)

3)利用弱分類器hm(x)擬合步驟二中得到的殘差rim;

4)計(jì)算估計(jì)乘子,如式(5):

(5)

5)根據(jù)得到的γm對(duì)Fm(x)進(jìn)行更新,Fm(x)=Fm-1(x)+γmhm(x);

6)判斷是否達(dá)到預(yù)先設(shè)置的迭代次數(shù)與精度要求,若滿足,則得到最終估計(jì)結(jié)果;否則,返回2)。

GBDT模型通過(guò)計(jì)算所有加法樹的均值,從而得到每個(gè)變量對(duì)因變量的相對(duì)重要程度,如式(6):

(6)

4 實(shí)證分析

4.1 模型訓(xùn)練

學(xué)習(xí)率、最大決策樹數(shù)量、樹深度這3個(gè)參數(shù)取值決定了GBDT模型的整體效果。GBDT模型中樹的深度不會(huì)超過(guò)log2(n),其中n為樣本數(shù)量,因此文中樹的深度取7。對(duì)學(xué)習(xí)率、最大決策樹數(shù)量取值進(jìn)行調(diào)參,探究不同取值對(duì)模型得分的影響,從而獲得最優(yōu)參數(shù)值,如圖3。由圖3可看出:有多處位置可取到最大值,最終文中模型的學(xué)習(xí)率取值0.15,樹的數(shù)量取值82。

圖3 參數(shù)調(diào)整

4.2 模型結(jié)果

經(jīng)訓(xùn)練后的模型擬合效果穩(wěn)定,各變量對(duì)交通事故嚴(yán)重程度影響的重要度見表5。

表5 變量重要度排序

建成環(huán)境變量中:土地利用混合度、快速路及主干路密度、次干路及支路密度、人口密度及可達(dá)性的累計(jì)重要度達(dá)60.57%。故這5個(gè)變量是解釋建成環(huán)境對(duì)交通事故嚴(yán)重程度影響的主要因素。其中土地利用混合度對(duì)交通事故嚴(yán)重程度的影響最大,為14.29%;其次是快速路及主干路密度、次干路及支路密度、人口密度和可達(dá)性,分別占比12.43%、11.54%、11.35%、10.96%;剩余變量重要度較低,均低于10%,這說(shuō)明這些變量對(duì)交通事故嚴(yán)重程度影響較小。道路情況與文獻(xiàn)[14]研究結(jié)果中的道路物理隔離與車道數(shù)對(duì)事故嚴(yán)重程度不顯著相吻合;事故發(fā)生時(shí)間變量與文獻(xiàn)[15]研究結(jié)果中的時(shí)間段對(duì)事故嚴(yán)重程度的重要度較低的結(jié)果一致。土地利用混合度、快速路及主干路密度、次干路及支路密度、人口密度及可達(dá)性對(duì)交通事故嚴(yán)重程度的影響曲線見圖4。

圖4 事故嚴(yán)重程度影響曲線

由圖4可知:土地利用混合度在1.72處發(fā)生較大突變,土地利用混合度高的地方,車輛的速度會(huì)受到限制,從而導(dǎo)致較低的事故嚴(yán)重程度;快速路及主干路密度在2.2 km/km2處發(fā)生較大突變,這表明快速路及主干路密度在超過(guò)該閾值后,對(duì)事故嚴(yán)重程度將產(chǎn)生較大影響,因此快速路及主干路密度在前期規(guī)劃中應(yīng)考慮安全合理取值;人口密度方面,在約1 600人/km2和19 000人/km2處發(fā)生突變,在該區(qū)間內(nèi)基本對(duì)事故嚴(yán)重程度影響基本不變,因此人口密度只有達(dá)到閾值后對(duì)事故嚴(yán)重程度的影響才能得到凸顯;次干路及支路密度對(duì)事故嚴(yán)重程度的影響與快速路及主干路相反,整體上表現(xiàn)出下降的趨勢(shì),因此從交通安全的角度出發(fā),次干路及支路的建設(shè)應(yīng)受到重視?？蛇_(dá)性對(duì)交通事故的影響呈現(xiàn)出上升的趨勢(shì),最后趨于穩(wěn)定;即隨著距離市中心距離越遠(yuǎn),事故嚴(yán)重程度越高,且在超過(guò)閾值后基本穩(wěn)定。

5 結(jié) 語(yǔ)

筆者基于重慶市某區(qū)的交通事故數(shù)據(jù),運(yùn)用GBDT模型探索了建成環(huán)境與交通事故嚴(yán)重程度的非線性關(guān)系;得到各變量的重要程度排序,同時(shí)給出了重要度較高的變量對(duì)事故嚴(yán)重程度的影響曲線。從研究結(jié)果可看出:① 建成環(huán)境變量在事故嚴(yán)重程度方面具有重要作用,其中土地利用混合度、快速路及主干路密度、次干路及支路密度、人口密度及可達(dá)性的集體貢獻(xiàn)占比60.57%;② GBDT模型能更加精確地捕獲各變量對(duì)事故嚴(yán)重程度影響的非線性關(guān)系及閾值效應(yīng)。