許倫輝，郭雅婷

（華南理工大學(xué) 土木與交通學(xué)院，廣州510641）

網(wǎng)約車服務(wù)，指以互聯(lián)網(wǎng)技術(shù)構(gòu)建服務(wù)平臺，整合供需，通過移動設(shè)備與應(yīng)用連接出行乘客與車輛和駕駛員的服務(wù)[1]。2018年，我國主要城市的網(wǎng)約車應(yīng)用總滲透率高達(dá)11%[2]。海量的網(wǎng)約車數(shù)據(jù)，為出租車供需研究提供了新思路，也為新形勢下提高差異化出行的供需匹配效率提供了研究基礎(chǔ)[3]。

1 網(wǎng)約車需求分析

1.1 網(wǎng)約車供需研究現(xiàn)狀

從供需平衡角度，出租車供需研究分供給預(yù)測、需求預(yù)測與供需缺口預(yù)測三類。

1）供給預(yù)測 常采用經(jīng)典的供需函數(shù)關(guān)系開展宏觀研究。

2）供需缺口預(yù)測 適用于區(qū)域內(nèi)供不應(yīng)求的情況，如文獻(xiàn)[4]基于時空數(shù)據(jù)挖掘的深度學(xué)習(xí)網(wǎng)約車需求預(yù)測模型。

3）需求預(yù)測 與交通流預(yù)測類似，即網(wǎng)約車出行需求的時空狀態(tài)分布預(yù)測，分為模型驅(qū)動、數(shù)據(jù)驅(qū)動兩類。其中，模型驅(qū)動主要通過建立交通各項參數(shù)間的數(shù)學(xué)模型以實現(xiàn)交通流預(yù)測[5]；數(shù)據(jù)驅(qū)動即機器學(xué)習(xí)，通過對大數(shù)據(jù)進(jìn)行智能計算以挖掘數(shù)據(jù)的隱含信息。文獻(xiàn)[6]提出單個、組合及提升等多種決策樹算法進(jìn)行網(wǎng)約車數(shù)據(jù)的短時需求預(yù)測；文獻(xiàn)[7]用LASSO（least absolute shrinkage and selection operator）進(jìn)行數(shù)據(jù)特征排序，建立基于隨機數(shù)森林算法的預(yù)測模型；文獻(xiàn)[8]以紐約市出租車數(shù)據(jù)為基礎(chǔ)，提出長短期記憶神經(jīng)網(wǎng)絡(luò)LSTM 的出租車需求模型。然而，傳統(tǒng)LSTM 在網(wǎng)絡(luò)參數(shù)尋優(yōu)中采用反向傳播BPTT（back propagation through time）算法，復(fù)雜度高且存在易收斂于局部最優(yōu)解。故在此，針對上述缺陷，提出了采用收斂速度快、全局搜索能力強的灰狼優(yōu)化算法對長短期記憶網(wǎng)絡(luò)進(jìn)行參數(shù)優(yōu)化，建立GWO-LSTM 網(wǎng)約車出行需求短時預(yù)測模型進(jìn)行改進(jìn)。

為把握乘客出行需求特性，以滴滴出行2016年1月份M 市的訂單數(shù)據(jù)作為研究對象，進(jìn)行數(shù)據(jù)分析。該數(shù)據(jù)含30 d日均40 萬條的訂單量、66 個出行區(qū)域、各區(qū)域內(nèi)信息點PoI（point of information）類別及數(shù)量、道路擁堵程度與天氣狀況等。通過對超出定義時空范圍內(nèi)數(shù)據(jù)進(jìn)行刪除，短時間內(nèi)重復(fù)訂單數(shù)據(jù)識別與剔除、噪聲平均化填充等操作進(jìn)行數(shù)據(jù)清洗。對出行數(shù)據(jù)進(jìn)行區(qū)域匹配、15 min 時間分片的短時需求量統(tǒng)計。

1.2 網(wǎng)約車需求時變規(guī)律

在工作日與非工作日，居民對網(wǎng)約車需求的時變規(guī)律如圖1所示，圖中將一天24 h 分為0～95 共96 個時間段，每個時段15 min。在工作日，居民對于網(wǎng)約車的出行需求趨勢與城市交通早晚高峰一致。7：30—9：30 為早高峰，需求量最大；17：00—19：00為晚高峰。其中周一和周五較其他工作日需求量要大。周末白天無明顯早高峰，總體變化平穩(wěn)；人們傾向于下午至夜間出行，總需求量與工作日相差不大。

圖1 工作日與非工作日網(wǎng)約車需求時變圖Fig.1 Time-varying graph of demand for network booking taxi on weekdays and non weekdays

1.3 網(wǎng)約車需求區(qū)域分布規(guī)律

網(wǎng)約車需求的區(qū)域分布與土地利用類型相關(guān)，土地利用情況主要反映在POI 數(shù)量及類型上。通過對比分析各區(qū)域數(shù)據(jù)，主要變化類型如圖2所示。

圖2 典型區(qū)域網(wǎng)約車需求時變圖Fig.2 Time-varying graph of demand for network booking taxi in typical regions

圖中，類型1 為商住混用型，該類區(qū)域早高峰與晚高峰對于網(wǎng)約車需求量較大且數(shù)量相當(dāng)；類型2 為以居住用地為主的典型區(qū)域，網(wǎng)約車需求以通勤出行為主，需求數(shù)劇增；類型3 主要為休閑娛樂類服務(wù)場所用地，以夜間需求為主；類型4 主要為商務(wù)用地，相較于其它區(qū)域，晚高峰的網(wǎng)約車需求量較大。

1.4 影響因素相關(guān)性分析

網(wǎng)約車需求量受時間、空間、交通擁堵狀況、天氣、溫度與污染指數(shù)等外部因素的影響。選用spearman 相關(guān)系數(shù)進(jìn)行分析，即

相關(guān)性分析結(jié)果見表1。

表1 網(wǎng)約車需求的影響因素相關(guān)性分析Tab.1 Correlation analysis of influencing factors of demand for network booking taxi

2 網(wǎng)約車需求短時預(yù)測模型

2.1 長短期記憶神經(jīng)網(wǎng)絡(luò)

循環(huán)神經(jīng)網(wǎng)絡(luò)RNN（recurrent neural networ）是以序列數(shù)據(jù)為輸入，在序列的演進(jìn)方向進(jìn)行遞歸且所有循環(huán)單元按鏈?zhǔn)竭B接的遞歸神經(jīng)網(wǎng)絡(luò)[9]。基于RNN 進(jìn)行改進(jìn)，通過引入靈活可控的自循環(huán)設(shè)計，LSTM 產(chǎn)生了讓梯度能夠可持續(xù)長期流動的途徑，可從輸入的歷史數(shù)據(jù)中學(xué)習(xí)長時間的計算記憶信息。

LSTM 內(nèi)部結(jié)構(gòu)如圖3所示，遺忘門ft，輸入門it，記憶單元ct和輸出門ot組成了LSTM 的內(nèi)部結(jié)構(gòu)，其中σ 為激活函數(shù)sigmoid。

圖3 LSTM 內(nèi)部結(jié)構(gòu)Fig.3 LSTM internal structure

具體計算如下：

遺忘門ft決定保留或丟棄歷史信息，輸出區(qū)間為［0，1］，其中1 表示保留，0 表示全部丟棄。ft為

輸入門it添加當(dāng)前時刻有效輸入，決定了xt添加到信息流中的程度，sigmoid 函數(shù)部分決定需輸入的值。tanh 函數(shù)部分創(chuàng)建臨時狀態(tài)ct′為

記憶單元ct更新網(wǎng)絡(luò)細(xì)胞狀態(tài)：

輸出門ot控制可用于下一層網(wǎng)絡(luò)更新而輸出的信息：

式（2）～式（7）中：ωhf，ωxf，bf，ωhi，ωxi，bi，ωhc，ωxc，bc，ωho，ωxo，bo分別為門限單元中相應(yīng)下標(biāo)的權(quán)重矩陣和偏置量。標(biāo)準(zhǔn)LSTM 采用式（2）～式（7）進(jìn)行計算得到模型輸出數(shù)據(jù)，并利用基于時間的BPTT 算法訓(xùn)練以求解上述參數(shù)。但在試驗過程中，BPTT 算法常收斂于局部最優(yōu)解，且計算復(fù)雜度較高。

2.2 灰狼優(yōu)化算法

灰狼優(yōu)化算法根據(jù)灰狼群體協(xié)作機制，模擬灰狼捕食獵物以實現(xiàn)目標(biāo)，具有較強的收斂性、參數(shù)少、易實現(xiàn)，且能收斂于全局最優(yōu)解[10]。在一定數(shù)量的群體中，按照職責(zé)分工，灰狼被分為α，β，δ 及ω共4 個等級。α 是唯一的領(lǐng)導(dǎo)者，負(fù)責(zé)群體決策；β級狼輔助領(lǐng)導(dǎo)者管理狼群，也是領(lǐng)導(dǎo)者替補；δ 服從α 與β 的命令，負(fù)責(zé)偵察、捕獵及看護(hù)等；最底層ω服從上級安排并維系種族平衡。

GWO 優(yōu)化的算法步驟包括尋找、跟蹤、包圍、狩獵以及攻擊獵物等。圍獵時，群體中每一只狼即為潛在解，按照等級，α，β，δ，ω 優(yōu)先級依次遞減。其他個體根據(jù)前三只狼的位置，計算與其距離，并不斷靠近更新位置以完成尋優(yōu)。在Q 維空間中，各個體Wi=（wi，1，wi，2，…，wi，q）組成狼群Wi=｛W1，W2，…，WN｝。有

其中

式中：Dm為灰狼Wi在第t 次迭代時分別距離α，β，δ狼的位置；Ck，Ak為系數(shù)向量；r1，r2均為［0，1］的隨機取值；Wik′為灰狼向α，β，δ 靠近時移動的矢量；Wi（t+1）為向3 只高等級灰狼移動的矢量的平均值，得到灰狼Wi下一次迭代的移動值；a 為收斂因子，隨迭代次數(shù)t 的增加線性減小，取a 值范圍［0，2］。

2.3 GWO-LSTM 網(wǎng)約車需求預(yù)測模型

基于GWO-LSTM 的網(wǎng)約車需求預(yù)測流程如圖4所示，該預(yù)測模型由數(shù)據(jù)處理、參數(shù)優(yōu)化、網(wǎng)絡(luò)訓(xùn)練與數(shù)據(jù)預(yù)測四大部分構(gòu)成。參數(shù)優(yōu)化部分采用GWO 算法優(yōu)化LSTM 網(wǎng)絡(luò)中的參數(shù)，使模型得出的預(yù)測值最大限度逼近真實值。為得到最小目標(biāo)函數(shù)，適應(yīng)度函數(shù)設(shè)為

式中：y′為實測數(shù)據(jù)；y 為模型輸出；S 為訓(xùn)練長度。

圖4 基于GWO-LSTM 的網(wǎng)約車需求預(yù)測流程Fig.4 Flow chart of forecasting demand for network booking taxi based on GWO-LSTM

訓(xùn)練數(shù)據(jù)時，先確定LSTM 中輸入層的歷史數(shù)據(jù)步長lookback、隱藏層單元個數(shù)hidden 及預(yù)測步長predstep；確定GWO 的迭代次數(shù)、維度以及灰狼個數(shù)。通過隨機分布函數(shù)，對灰狼位置進(jìn)行種群初始化，并分解為LSTM 網(wǎng)絡(luò)所求的權(quán)重矩陣與偏置值等未知參數(shù)。將參數(shù)代入網(wǎng)絡(luò)中計算相應(yīng)的輸出值y 并根據(jù)式（13）計算各個體的適應(yīng)度值。對比適應(yīng)度值，標(biāo)記其中適應(yīng)度最佳的前3 個體位置為Wα，Wβ，Wδ。當(dāng)?shù)螖?shù)不斷迭加，由計算式（8）～式（12）更新每一個體的位置，并多次計算適應(yīng)度值，直至迭代結(jié)束。返回全局最佳的Wα即為LSTM 記憶網(wǎng)絡(luò)參數(shù)最優(yōu)解。最終，將測試數(shù)據(jù)輸入到最優(yōu)參數(shù)解的訓(xùn)練網(wǎng)絡(luò)中，得到預(yù)測結(jié)果，將均方根誤差RMSE（root mean square error）與平均絕對誤差MAE（mean absolute error）作為模型評估指標(biāo)。其計算如下：

式中：Yi′為測試集實測數(shù)據(jù)；Yi為預(yù)測數(shù)據(jù)。

3 算例分析

3.1 網(wǎng)約車需求量數(shù)據(jù)處理

試驗仍以滴滴平臺M 市2016年1月份的訂單數(shù)據(jù)為例，70%為訓(xùn)練數(shù)據(jù)，30%用作預(yù)測。通過初步數(shù)據(jù)分析，區(qū)域劃分、時間分片得網(wǎng)約車需求量時間序列。采用spearman 相關(guān)系數(shù)法對網(wǎng)約車需求量的影響因素進(jìn)行相關(guān)性分析，選定PoI 數(shù)量與5個歷史時間段數(shù)據(jù)為輸入指標(biāo)，并采用min-max 標(biāo)準(zhǔn)化方法對數(shù)據(jù)進(jìn)行歸一化處理，即

3.2 試驗對比模型及平臺搭建

對比預(yù)測模型為標(biāo)準(zhǔn)LSTM 和BPNN。標(biāo)準(zhǔn)LSTM 通過BPTT 算法尋找最優(yōu)參數(shù)解，而本文模型采用灰狼優(yōu)化算法尋找最優(yōu)參數(shù)。BP 神經(jīng)網(wǎng)絡(luò)為多層前饋神經(jīng)網(wǎng)絡(luò)，在網(wǎng)絡(luò)訓(xùn)練過程中以誤差逆?zhèn)鞑サ姆椒ㄗ钚』`差，并調(diào)整優(yōu)化連接權(quán)值以減小預(yù)測誤差。

這3 個對比模型均使用python 語言編程實現(xiàn)。為建立預(yù)測模型而搭建實驗平臺的軟硬件信息如下：系統(tǒng)為Windows 10（64 位），處理器為AMD Ryzen 7 PRO 2700U，內(nèi)存為8.0 GB；程序語言版本為Python 3.6.7，開發(fā)環(huán)境采用Anaconda 包jupyter 1.0.0。深度學(xué)習(xí)框架采用CPU 版tensorflow 2.0 中的keras 2.3.1。

3.3 模型訓(xùn)練參數(shù)優(yōu)化

訓(xùn)練過程如下：

步驟1確定lookback，hidden 與predstep，歷史數(shù)據(jù)步長為6 （區(qū)域PoI 數(shù)量qp及5 個歷史時間網(wǎng)約車需求量，預(yù)測步長為1，即當(dāng)前預(yù)測時間段網(wǎng)約車需求量。由經(jīng)驗法則，隱藏神經(jīng)元數(shù)量為

步驟2確定灰狼算法搜索空間維度、迭代次數(shù)及個體數(shù)：設(shè)灰狼個體數(shù)為50，迭代次數(shù)為800。GWO-LSTM 訓(xùn)練網(wǎng)絡(luò)參數(shù)優(yōu)化誤差曲線如圖5所示。

圖5 GWO-LSTM 訓(xùn)練網(wǎng)絡(luò)參數(shù)優(yōu)化誤差曲線Fig.5 Error curve of training network parameter optimization in GWO-SLTM

當(dāng)?shù)螖?shù)大于600 時，誤差值逐漸變小且趨于穩(wěn)定； 搜索空間維度為LSTM 各權(quán)重矩陣與偏置量的維度之和。如ωxi，維度hidden*lookback 為30，其他參數(shù)依此類推，因此得到灰狼群體中搜索空間維度為

步驟3初始化種群操作，采用隨機分布函數(shù)對群體進(jìn)行初始化，即為Wij～U（uj，lj）。其中，uj，lj分別為分布區(qū)間的上、下界，將得到的初始種群值分解為LSTM 參數(shù)。

步驟4當(dāng)?shù)螖?shù)≤800 時，計算每個灰狼個體參數(shù)相應(yīng)的輸出值y 及適應(yīng)度值，對比各適應(yīng)值并排序，取得最佳的3 個個體，依次標(biāo)記為Wtα，Wtβ，Wtδ，其中t 為當(dāng)前迭代次數(shù)。更新其余灰狼位置，并迭代當(dāng)前步驟，直至迭代結(jié)束，尋找到最優(yōu)解Wα′，進(jìn)而分解為LSTM 網(wǎng)絡(luò)參數(shù)。

步驟5預(yù)測階段，將測試集數(shù)據(jù)輸入到網(wǎng)絡(luò)中，得到預(yù)測值與實際值進(jìn)行誤差計算，即得模型預(yù)測準(zhǔn)確度。

3.4 模型預(yù)測結(jié)果分析

GWO-LSTM，標(biāo)準(zhǔn)LSTM 和BPNN 三種模型對示例區(qū)域全天網(wǎng)約車需求量預(yù)測結(jié)果的對比如圖6所示。3 種模型預(yù)測趨勢與真實數(shù)據(jù)變化趨勢均呈一致性；在早晚高峰等起伏較大的時間段，3 種模型均體現(xiàn)了不同程度的預(yù)測偏差。相對地，所提出的GWO-LSTM 模型相較其他模型，誤差較小，預(yù)測精準(zhǔn)度較高，預(yù)測效果較優(yōu)。其數(shù)值表現(xiàn)見表2。

圖6 預(yù)測結(jié)果對比Fig.6 Comparison of prediction results

表2 預(yù)測結(jié)果評價Tab.2 Evaluation of prediction results

對于模型評估指標(biāo)RMSE 與MAE 而言，數(shù)值越小則表明模型預(yù)測效果越好。由表2可知，GWO算法改進(jìn)的LSTM 預(yù)測值的RMSE 為9.8858，相對于傳統(tǒng)的LSTM 網(wǎng)絡(luò)提升了34.45%，相對于BP 神經(jīng)網(wǎng)絡(luò)提升了48.16%； 在MAE 方面，GWO-LSTM預(yù)測結(jié)果相對于真實值為7.7708，相較LSTM 提升了36.89%，相較于BPNN 提升了52.12%，預(yù)測效果顯著。

4 結(jié)語

基于實際數(shù)據(jù)，分析了網(wǎng)約車需求時空分布特性及影響因素與網(wǎng)約車需求量的相關(guān)性；針對標(biāo)準(zhǔn)LSTM 模型在時序數(shù)據(jù)預(yù)測中存在的問題，提出了GWO 改進(jìn)LSTM 網(wǎng)絡(luò)參數(shù)的網(wǎng)約車短時需求預(yù)測模型，并對模型效果進(jìn)行實例驗證。結(jié)果表明，該模型相較于標(biāo)準(zhǔn)LSTM 和BPNN，預(yù)測效果均有較大提升；尋優(yōu)過程中，灰狼優(yōu)化算法有種群多樣性差與后期收斂速度慢等缺陷，將對灰狼優(yōu)化算法進(jìn)行改進(jìn)以優(yōu)化模型。