城市轨道交通运营客流预测的探讨

城市轨道交通运营客流预测的探讨

客流预测是城市轨道交通规划设计和运营管理的基本依据，是城市轨道交通建设过程的重要组成部分。传统的客流预测方法中应用最广泛的是时间序列方法,该方法最早起源于1927年,直到70年代末到80年代中后期才在我国得以深入研究和应用。近年来自回归移动平均模型(ARIMA模型)是预测领域的主要方法之一,是时间序列方法中许多基本思想的基础,姚亚夫等应用ARIMA模型进行交通流量趋势预测,得到较好的预测结果;池启水等将ARIMA模型应用于煤炭消费预测。ARIMA方法现在已经有一套完整的建模方法,但是具有一定的局限性,它假定分析的时间序列是平稳的,对于非平稳的时间序列预测结果不能达到最好,要将其转换为平稳序列。神经网络是一种新型的预测方法,适合复杂、非线性的条件,摆脱了建立精确数学模型的困扰,具有良好的自组织性、自适应性,有很强的学习能力、抗干扰能力等,目前已经建立了BP神经网络模型、多层感知器神经网络、径向基函数神经网络、动态神经网络等神经网络模型,并且这些模型已经广泛应用于预测领域。例如,Deng等利用BP神经网络模型对公路客运量预测,取得了较好的预测结果;Xiao等利用RBF神经网络预测高速公路短期交通流量,预测结果显示RBF神经网络能够较准确地预测高速交通流量;李晓俊等提出了基于径向基神经网络的铁路短期客流预测模型。近年来,许多学者考虑将不同预测方法组合在一起,建立组合预测模型从而提高预测的精度。例如Tsai等提出多时序单位神经网络模型和并行组合神经网络,对铁路客运需求进行预测;Zheng等将贝叶斯法和神经网络组合对短期的高速公路交通流预测,结果表明组合模型的预测效果比单个神经网络模型更好。城市轨道交通的客流具有一定的时序特征,即客流量随着年、月、周、日的变化而变化的规律,其中最明显的时序特征是客流以周为周期,每个周期中客流具有相似的变化趋势,由于客流易于受到多种因素的影响,每个周期内客流变化大小均不相同,具有一定的非平稳时间序列特征,因此,本文拟采用ARIMA-RBF组合模型来描述和预测同时具有周期性和非平稳时间序列特征的客流数据,主要做了以下几个方面的工作:(1)从客流量折线图中总结客流的时序特征,选出用于预测的历史客流量数据;(2)用ARIMA模型对客流进行预测,得到包含客流的线性特征的客流预测结果和包含非线性特征的ARIMA预测误差;(3)用RBF神经网络模型对ARIMA预测误差进行预测,建立ARIMA-RBF组合预测模型。1后移算子b一般来说原始的城市轨道交通客流量数据是非平稳时间系列,需要通过(d阶)差分计算将其转化为平稳时间序列{yt}:式中,B被称为后移算子。对于平稳化处理后的序列{yt}采用ARMA(p,q)模型就可以得到ARIMA(p,q)模型,设式中,yt为t时刻的观测值,{εt}为零均值白噪声序列,ϕi(i=1,2,3,…,p)和θj(j=1,2,3,…,q)为参数,p为自回归阶数,q为移动平均阶数,ϕ(B)和θ(B)为算子多项式,并且ARIMA(p,d,q)模型可以表示为:2隐蔽层神经元个数RBF神经网络是一种逼近性能良好的前向神经网络,结构简洁,具有学习速度快,还存在局部极小等优点,RBF神经网络具有3层结构,即:输入层、隐含层和输出层。通常,一个具有m个输入值、n个隐含节点及一个输出值的RBF神经网络结构如图1所示。在图中,{x1,x2,…,xn}为输入数据,输入层的神经元个数由输入数据的种类决定;F(x)为单输出,(ω1,ω2,…,ωn)为隐含层与输出层之间的连接权重。输入层节点只是将数据传递到隐含层,所以输入层节点到隐含层节点间的连接权重为1,隐含层节点i的输出为ui=R[∥X-Ci∥σi]ui=R[∥X−Ci∥σi],(8)其中,i=1,2,…,m;R为径向基函数,Ci为第i个神经元的径向基函数中心,σi为第i个神经元的径向基函数的宽度参数;‖X-Ci‖为欧式范数。隐含层径向基函数可以选用不同的函数,通常使用高斯函数作为径向基函数,即,R(x)=exp(-x22)R(x)=exp(−x22)。(9)所以RBF神经网络的输出可以表示为:F(x)=m∑i=0wiexp(-∥X-Ci∥22σ2i)F(x)=∑i=0mwiexp(−∥X−Ci∥22σ2i),(10)式中,F(x)为RBF神经网络的输出。3客流序列模型预测和描述实际的客流变化不仅有线性规律也有非线性规律,ARIMA模型和RBF神经网络模型分别适用于线性和非线性条件下的客流预测,本文将这两种模型组合,所建立的组合模型既能把握客流序列的线性规律,又能考虑到客流序列的非线性规律。设一个客流序列由线性和非线性组成构成,即:式中,Lt和Nt分别表示客流序列的线性和非线性组成。首先用ARIMA模型对客流序列中的线性组成进行预测,得到具有非线性特征的ARIMA模型的预测误差et:et=yt-ˆLtet=yt−Lˆt,(12)式中,ˆLLˆt为t时刻ARIMA模型的预测值。然后将et其作为RBF神经网络模型的输入单元,则RBF神经网络模型可以表示为:式中,f为神经网络拟合的非线性方程,εt为随机误差,设RBF神经网络的预测结果为ˆΝNˆ,则有:ˆyt=ˆLt+ˆΝtyˆt=Lˆt+Nˆt,(14)式中ˆyyˆt为t时刻组合模型的客流量预测值。由以上分析可知,ARIMA-RBF组合模型预测过程为:首先利用ARIMA模型对客流序列的线性组成进行预测,再用RBF神经网络对包含非线性特征的ARIMA预测误差进行预测,将两部分预测结果相加得到最终的客流预测值,预测流程见图2。42年5月28日6月3日内部客流数据预测本文选取北京地铁1号线2011年2月21日～2012年5月27日的客流量数据预测1号线2012年5月28日～6月3日的客流量,数据由北京市地铁运营有限公司提供。下面对1号线2011年2月21日～2012年5月27日客流量数据进行分析。4.1数据分析4.1.1节假日因素:传统线路是由季节性因素组成的,传统是个本文选取北京地铁1号线2011年2月21日～2012年5月27日的日客流量进行分析,客流的变化见图3。由图3可以看出,客流量基本呈现以周为周期变化的趋势。节假日因素对客流量影响较大,节假日期间客流变化与平常日客流明显不同。元旦、清明节、五一劳动节、端午节、国庆节和中秋节期间客流量较大,五一劳动节和国庆节期间客流量能达到全年最大客流量,春节期间的客流量大幅度降低,出现全年最低值。并且客流变化规律也与平常日不同,除了我国法定节假日外,3月8日妇女节、圣诞节也会出现客流较大的变化。节假日因素不仅影响着节假日期间客流的变化,也会导致节前一周的客流量发生较大的变化。4.1.2日常日客流趋势本文是对2012年5月28日～6月3日平常日的客流进行预测,节假日对客流变化的影响较大,所以本文将1号线2011年2月21日～2012年5月27日期间,受到元旦、3月8日妇女节、清明节、五一、端午节、国庆节、中秋节和春节、圣诞节等节假日影响的节假日及节前一周的客流量剔除,选取2011年2月21日～2012年5月27日期间不包含任何节假日的294个平常日,共计42个完整周的日客流量进行分析,客流变化趋势如图4所示。可以看到,1号线在2011年2月21日～2012年5月27日期间平常日客流量以周为周期发生规律性的变化:每一周内周一到周四是工作日,这几天的客流量相差不是很大,周五的客流量明显变多,周六、周日的客流量又呈明显的下降趋势,且远远低于工作日的客流量;另外可以看出平常日客流量整体呈现一定的非平稳特征,下面用数学方法分析平常日客流序列的平稳性和周期性。4.1.3日常日客流单位根检验自相关函数图、偏自相关函数图和单位根检验是考察时间序列平稳性最重要的工具,本文利用SPSS软件绘制2011年2月21日～2012年5月27日期间284个平常日共计42个完整周的客流量的自相关函数图和偏自相关函数图,见图5。从图5中可以看出,平常日客流的自相关函数图呈明显的周期性波动,周期长度为7,以周期长度及其整数倍数为阶数的自相关和偏自相关函数值均显著地不为0,说明此时间序列为具有周期性的非平稳时间序列。利用Eviews软件对平常日客流进行单位根检验结果见表1。从表1中可以看出单位根统计量ADF=-2.142258,大于Eviews给出的显著性水平1%～10%的ADF临界值,所以认为该平常日客流序列是非平稳的。4.1.4时间序列的平稳化检验由以上分析可知平常日客流序列是以一周为周期的非平稳时间序列,需要将客流数据进行平稳化处理,即进行周期差分。定义经过平稳化处理后的数据序列为{yt},{yt}共有277个数据,{yt,t=1,2,…,277}的自相关函数和偏自相关函数图见图6。可以看出,经过周期差分后得到的时间序列{yt,t=1,2,…,277}的自相关函数在9阶后截尾,偏自相关函数拖尾,说明此序列已成为平稳时间序列,可以用来构建ARIMA模型。对平稳化处理后的时间序列{yt}进行单位根检验,结果见表2。从表2中可以看出{yt}的单位根统计量ADF=-4.305616,小于Eviews给出的显著性水平1%～10%的ADF临界值,所以认为该序列是平稳的。通过以上分析可知,节假日因素对城市轨道交通客流的影响是较大的,不同节假日期间客流的变化规律是不同的;平常日客流虽然呈现以周为周期的变化规律,但是属于非平稳时间序列,经过周期差分处理后的序列成为平稳时间序列。4.2示例预测4.2.1时间/时间本文选取2011年2月21日～2012年5月27日期间294个平常日共计42个完整周的客流量数据预测2012年5月28日～6月3日的客流量,预测过程见图7。284个平常日客流分为3部分:2011年2月21日～2011年12月18日期间224个平常日、2012年1月9日～2012年3月25日期间35个平常日、2012年4月9日～2012年5月27日期间35个平常日。首先利用2011年2月21日～2011年12月18日期间平常日的客流量数据构建ARIMA模型,预测2012年1月9日～2012年6月3日平常日的客流量,由此可以得到2012年1月9日～2012年6月3日期间平常日客流量的预测值,通过和实际客流量对比,再将得到的2012年1月9日～2012年3月25日平常日的ARIMA模型的预测误差作为RBF神经网络的输入数据,2012年4月9日～2012年5月27日期间平常日的ARIMA模型的预测误差作为RBF神经网络的期望输出,从而构建RBF神经网络来预测2012年5月28日～2012年6月3日的ARIMA模型预测客流量的误差,将ARIMA模型的预测值与RBF模型预测的误差值相加得到最终的客流预测值。4.2.2平均绝对相对比本文分别用ARIMA模型、RBF神经网络模型和ARIMA-RBF组合预测模型对北京地铁1号线2012年5月28日～2012年6月3日的客流量进行预测,用平均绝对相对百分比误差(MAPE)作为预测结果的评价指标,其计算公式为:ΜAΡE=|(yi-ˆyi)yi|×100,(i=1,2,3,⋯,7)MAPE=∣∣(yi−yˆi)yi∣∣×100,(i=1,2,3,⋯,7),(15)式中,yi为1号线路所预测的一周内第i天客流量的真实值,ˆyyˆi为预测值,三种模型所预测的一周内日客流量的平均绝对相对百分比误差如图8所示。从图8中可以看出,三种模型对北京地铁1号线2012年5月28日～6月3日的日客流量的预测结果中,ARIMA-RBF模型预测结果的预测误差最小,平均每天的预测误差为3%,ARIMA模型的预测误差最大,平均每天的预测误差为4%,RBF神经网络的预测误差介于ARIMA模型和ARIMA-RBF模型