基于IC卡數據的公交OD推導與應用

張鐵巖 高洪振 禚保玲 胡 倩

（青島市城市規劃設計研究院，青島 266071）

引言

公交客流出行調查難度較大，隨著公交IC卡收費系統和車輛運行調度管理系統的發展普及，產生了大量公交客流和車輛運行數據，基于IC卡數據的公交非集計分析手段在國內外各城市中得到廣泛研究和應用[1-3]，通過采用公交IC卡數據融合技術，推算整個公交系統乘客出行起訖點，可以更加準確掌握公交客流的時空分布特征及公交運力供需關系，同時也將節省大量人工調查時間和成本。

由于各城市在公交刷卡收費模式和車輛運行調度管理系統完善程度上存在差異，造成了IC卡數據分析方法差別較大。對于分段計價線路，乘客上下車均需要刷卡，通過關聯車輛GPS數據和簡單的聚類排序則可以得到完整公交出行鏈信息。而目前國內公交運行系統主要采用的是單一票制，乘客僅在上車時刷卡，且刷卡記錄本身不包含上、下車站點位置信息。

在公交系統的IC卡分析技術研究中，CHUKKA[1]提出前期數據處理質量對融合結果有重要影響；PAEZ[2]、TREPANIER[3]從公交底層數據存儲、IC卡數據深度應用角度對公交數據提出規范化要求；Chu[4]、LEE[5]、Barry[6]等人分布利用GPS數據和IC卡數據，提出融合分析方法；吳子嘯[7]以鄭州為例詳細闡述了公交數據與交通模型之間關系和應用；Takahiko[8]、馬曉磊[9]在公交IC卡分析中提出貝葉斯算法，以彌補信息缺失問題；陳紹輝等[10]利用禁忌搜索算法對系統偏差數據進行二次匹配，提升了匹配精度；劉文芳[11]、胡郁蔥[12]、吳祥國[13]、楊萬波[14]等人分別針對成都、廣州、重慶、佛山數據特點進行了針對性的分析研究；陳學武及其合作者[15-17]針對不同的數據場景，分別闡述了公交IC卡數據融合方法。

青島公交收費系統采用的是單一票制，同時擁有較為完善的GPS數據和運營調度數據，本文結合青島公交系統數據特點，上車站點識別方法采用IC卡數據與GPS數據關聯推算，下車站點識別方法采用出行鏈算法。并利用交通宏觀需求分析軟件EMME，對模型參數進行標定和交通分配，將分析結果與手機信令數據和建筑數據進行相關性分析，證明了結果有效性。

1 數據特征及分析流程

青島市公交IC數據主要字段包括線路名稱、卡號、刷卡時間、POS機編號、卡類型等，GPS數據主要字段包括線路名稱、車輛編號、到站時間、離站時間、站點編號。結合站點、車次信息表，將公交IC卡數據與GPS數據進行時間匹配，從而確定每位乘客上車刷卡線路和站點。本次采集公交數據為青島市2015年9月14日-20日一個完整周的IC卡和GPS數據，其中IC卡數據共計1361萬條記錄，日均記錄194萬條，平均每秒約22.5次刷卡記錄，共包含334條公交線路刷卡數據。GPS數據共計1396萬條記錄，日均記錄199萬條。公交IC卡刷卡數據樣表、公交GPS數據樣表分別見表1、表2。

為方便管理各類數據表，基礎數據存儲采用SQL Server、算法處理采用Python語言編程、交通分配采用EMME4中5.11模塊（Prepare for standard traffic or transit assignment）和5.31模塊（Standard transit assignment）。數據處理總流程如圖1所示。

上車站點推算流程如圖2所示，首先對GPS數據進行預處理，主要包括剔除重復報站、統一線路編號、規范時間格式及空值記錄清洗等工作，將處理好的公交車輛到離站位置信息表同IC卡數據按照車次和時間進行匹配，并記錄各線路原始刷卡量，進而對上車站點初始匹配數據再次進行數據清洗，包括剔除一次刷卡記錄數據，將數據按照卡號、時間排序，形成預處理后上車站點數據。

表1 公交IC卡刷卡數據樣表

表2 公交GPS數據樣表

圖1 數據處理流程

圖2 公交IC卡數據上車站點處理流程

下車站點推算采用的出行鏈算法，并且基于以下兩個假設條件：假設1：公交出行者第一次刷卡站點也是其當日出行終點；假設2：上一次刷卡線路中，離下一次刷卡上車站點最近的站即為上一次公交乘車下車點。

另外，對于換乘乘客來說，上一次公交乘車距離或路況不同會導致前后兩次刷卡時間差存在較大差異，單純采用兩次刷卡時間差或換乘距離作為判別換乘的約束參數，勢必會造成對部分長距離出行者誤判。基于以上考慮，本文采用換乘距離和速度兩個約束參數作為換乘點判別參數。即作為換乘點需滿足兩個條件：上一次乘車的下車點和下一次乘車上車點距離應在一定閾值范圍內，用Δl表示；上一次乘車出行距離和前后兩次乘車刷卡時間差比值應在合理范圍內，公式表示為：

l/(T-Δt）≤v臨界

式中，l為上一次乘車距離，T為兩次乘車刷卡時間差，Δt為換乘等待時間，v臨界為速度閾值。

下車站點推算具體步驟如下：

（1）基于以上假設，將上車站點數據按卡號、時間排序后，每次提取同一卡號相鄰兩條刷卡記錄。

（2）根據下一條記錄上車點B，搜索上一條記錄中對應線路的最近站點A2，作為上一次乘車下車點，并計算上一條記錄中上車站點A1和站點B站點間距離lA1B。

（3）若lA1B＞Δl則無法判斷下車點，轉（5）。

（4）若lA1B≤Δl則計算上一次出行距離l，若l/(T-Δt）≤v臨界則標記為一次出行，否則標記為換乘點。

（5）重復上述操作，直至最后一條記錄。

公交IC卡數據下車站點處理流程如圖3所示。

圖3 公交IC卡數據下車站點處理流程

2 多源數據校驗

眾所周知，公交出行量與公交站點和線路布設關系密切，而在便捷公交出行條件下，則與該區域的建筑規模及人口、崗位數量有直接關系，基于以上理解，選取中心城區范圍內有公交站點的1平方公里網格聚合數據，與相同網格范圍內手機信令數據、建筑地塊數據進行相關性分析，對公交IC卡分析結果開展多源數據校驗。

青島市手機信令數據基于2015年9月數據分析結果[18]（與公交IC卡數據時間一致），信令數據字段包括手機號偽碼、手機號前七位、信息采集時間、LAC區號、基站編號、基站經緯度；建筑數據為2015年青島市中心城區建筑地塊集合，共計14.9萬個建筑物，附帶建筑基底面積、層數、用地性質等屬性。

將青島市中心城區按照1km×1km網格劃分，分別統計各網格單元內公交IC卡乘降量、手機信令人口崗位數、建筑面積，并基于ArcGIS，采用自然斷點分級法（Natural Breaks），對比各指標整體空間分布特征，如圖4所示，各指標均總體呈現由南部沿海核心區向外圍遞減趨勢。

采用SPSS對上述網格數據指標進行Pearson 相關性分析，每類指標均包含576個有效樣本。結果表明：手機信令人口崗位數、建筑面積均與公交IC卡乘降量在0.01水平（雙側）上顯著相關，相關系數分別為0.778、0.736，屬于強正相關關系（表3），即手機信令人口崗位數或建筑面積越高，公交乘降量越高；反之，手機信令人口崗位數或建筑面積越低，公交乘降量越低（圖5）。通過對多源數據校驗表明：公交IC卡數據分析結果具有較高可信度，可以清晰反映城市人口分布、用地開發強度與公交客流規模關系，為精細化城市公交規劃、建設、管理提供科學決策依據。

圖4 青島中心城區1km網格公交乘降量、手機信令人口崗位、建筑面積空間分布

圖5 公交乘降量與手機信令人口崗位、建筑面積相關關系散點圖

表3 公交乘降量與手機信令人口崗位、建筑面積相關系數

表4 客流指標匯總

3 結果分析

工作日公交總體客流早高峰集中在7:20-8:20時段，占全日客流的12.3%；晚高峰集中在17:10-18:10，占全日客流的10.2%，晚高峰低于早高峰。學生卡和普通卡有明顯早晚高峰，符合通勤、通學的出行特點。普通卡刷卡量約占全部樣本的66%，決定了全部刷卡人口的總體趨勢。學生公交客運量占比最低，僅占總客運量的5%，受學生上下課時間制約，學生卡早高峰出現在6:30-7:30，明顯早于其他兩種卡類型，晚高峰出行時段在17:20-18:20，與普通卡時間疊加，且學生出行高度集中，早晚高峰兩個小時出行量占學生總出行量的50%，早晚高峰之間幾乎沒有出行。老年卡持有者出行時間聚集度低，且早高峰較滯后，晚高峰提前，說明老年人對出行時間要求較低，體現出老年人退休生活特點。客流指標匯總見表4；工作日不同卡類型客運量變化如圖6所示。

圖6 工作日不同卡類型客運量變化

從公交客流集散量空間分布上看，東岸城區李村商圈、臺東商圈、火車站周邊、香港路沿線、四流路-人民路-威海路沿線、延安路沿線都是客流集散量比較大的區域；北岸城區集散量較大的區域集中在正陽路沿線和汽車北站周邊，但單個站點集散量明顯低于東岸城區。排名前十位公交站點日集散量均在3萬人次以上，其中最大的客流集散點為臺東，日集散量達5.2萬人次；其次是維客廣場，日集散量為5萬人次（圖7），公交客流集散點前十位見表5。

圖7 公交客流集散點空間分布（單位：次/日）

表5 排名前十位公交客流集散點

基于交通宏觀需求分析軟件EMME對公交OD進行交通流量分配，得到公交流量分布圖（圖8），東岸老城區內客流十分密集，由南向北客流量逐步減少，客流總體分布呈偏心輻射式。其中，遼陽路以北客流以南北向運動為主，遼陽路及以南客流東西向運動也很明顯；北岸城區內部公交客流以正陽路為中軸向外輻射；西岸城區受地理因素的影響，公交客流總體呈現側“U”型分布特點。

青島中心城區圍繞膠州灣分別形成東岸城區、北岸城區、西岸城區三大組團，同青島城市空間發展格局一致，常規公交在三大組團中自成體系。結合線網和站點布局發現：北岸城區建成區范圍線網密度為1.9 km/km2，高于西岸城區的1.5 km/km2，同時北岸城區300m半徑站點覆蓋率為50.8%，高于西岸城區的46.8%，可見北岸城區公交覆蓋水平高于西岸城區，且兩城區人口規模接近，但北岸城區公交客運量遠低于西岸城區，表明北岸城區公交出行比例仍處于較低水平（2015年公交方式結構僅為12.08%），公共交通在城市交通中的重要作用尚未得到很好的發揮。隨著機動化水平的快速提升，北岸城區要使城市交通可持續發展，需要進一步加強引導，強化居民公交出行理念，增強公交對居民出行的吸引力。

4 結束語

本文以青島市為例，介紹了公交IC卡數據融合分析方法，并對結果進行了校驗和展示。公交IC卡收費系統已經在國內得到普及和應用，相對于傳統公交調查，基于該系統的數據融合分析技術，具有樣本量大、低成本見效快、長期持續檢測等優勢。在對公交需求深入分析的基礎上，通過整合公共交通系統的路網、線路、站點及各類運營信息，對公共交通系統供需關系進行建模和仿真，尋找公交服務在時空上的薄弱環節，可以為線網規劃和優化等實際工作提供決策支持。

圖8 公交客流分布圖