A modified cell transmission model (CTM) is proposed to depict the temporal-spatial evolution of traffic congestion on urban freeways. Specifically, drivers’ adaptive behaviors and the corresponding influence on traffic flows are emphasized. Two piecewise linear regression models are proposed to describe the relationship of flow and density (occupancy). Several types of cellular connections are designed to depict urban rapid roads with on/off-ramps and junctions. Based on the data collected on freeway of Queen Elizabeth, Ontario, Canada, we show that the new model provides a relatively higher accuracy of temporal-spatial evolution of traffic congestions.
1. Introduction
It is well known that traffic flow dynamics estimation is important to transport engineering. To get an accurate estimation, most existing models are based on macroscopic continuum traffic flow modeling. The study of macroscopic continuum traffic flow model began with the Lighthill-Whitham-Richards (LWR) model independently proposed by Lighthill and Whitham [1] and Richards [2]. The LWR model formulated a hyperbolic partial differential equation (PDE) to describe the temporal-spatial evolution of traffic flow on a homogeneous highway.
Cell transmission model (CTM) is one of convergent numerical approximations to LWR model and was developed by adopting a simplified trapezoidal relationship between traffic flow and traffic density to solve the scalar kinematic wave model [3, 4]. It assumes a constant free flow speed, vf, when density is lower than the critical traffic density and a constant backward shockwave speed, w, at a higher density.
Although the linear structure of CTM has the advantage of simplifying traffic control analysis, design, and data-estimation, the assumption of a simplified trapezoidal fundamental diagram seems too strong for general applications. Some more complex models (e.g., Switching Mode Model (SMM)) derived from CTM were presented to reflect the temporal-spatial dynamic features of traffic congestion [5]. Alternatively, the investigation of oscillating congested traffic on freeways based on CTM, Lo et al. [6], leads to a dynamic traffic control formulation designed as dynamic intersection signal control optimization.
It was gradually found that the simple first-order continuum model does not have the ability to explain the amplification of small disturbances on heavy traffic, because no stability condition can be derived from the model [7]. In addition to the simple continuum model, Payne [8] introduced a higher-order model that is derived from car-following model to overcome the shortcomings of LWR model. Based on Payne’s model, the deterministic and stochastic finite difference equations of validated high-order model were employed to describe freeway traffic flow dynamics, respectively [9–12]. Another higher-order continuum model that includes a dynamics equation developed from car-following theory in addition to the conservation equation was introduced by Jiang et al. in [7]. This model overcomes the backward travel problem that exists in many other higher-order continuum models and can explain nonlinear dynamical spatial-temporal phenomena on freeways like shock waves, rarefaction waves, stop-and-go waves, and local cluster effects.
The uncertainty of traffic flow also attracts increasing concerns. Alecsandru [13] developed a stochastic compositional model of the evolution of traffic flows on freeways to extend the Daganzo cell transmission model by defining sending and receiving functions explicitly as random variables and by also specifying the dynamics of the average speed in each cell. Differently, an entropy solution was proposed for the LWR model with a flow density relationship of being piecewise quadratic, continuous, and concave, but not differentiable at the junction points where two quadratic polynomials meet [14]. Several alternative models were also developed in [15, 16].
In this paper, a modified cell transmission model (CTM) is proposed to depict the temporal-spatial evolution of traffic congestion on urban freeways. Specifically, drivers’ adaptive behaviors and the corresponding influence on traffic flows are emphasized. Two piecewise linear regression models are used to describe the relationship of flow and density (occupancy). Several types of cellular connections are designed to depict urban rapid roads with on/off-ramps and junctions.
The rest of this paper is organized as follows. In Section 2, a new formation of CTM is presented. Section 3 discusses two critical issues when linear initialization of density and local perturbation analysis are conducted to reflect the dynamic and nonequilibrium properties of traffic flow. Section 4 addresses the two applications of spatial-temporal evolution analysis based on the proposed method with real detective data on unban freeways. Finally, main conclusions are summarized in Section 5.
2. An Improved Spatial-Temporal Finite Difference Model
The LWR model and the selected higher-order continuum model can be written in two PDEs as follows:(1a)∂ρ∂t+∂ρv∂x=gx,t,(1b)∂v∂t+v∂v∂x=veρ-vT-cρ,v∂pρ,v∂x,where t and x denote time and distance from the origin, ρ(x,t) is the traffic density, v(x,t) is the space mean speed, g(x,t) is net-merging traffic flow, T is the relaxation time, ve(ρ) is the steady-state speed-density relationship implied by car-following model, and p(ρ,v) denotes the traffic pressure from upstream, and c(ρ,v) means the coefficient of traffic pressure.
2.1. Finite Difference of the First-Order Equation
For the first-order model, CTM is an efficient finite difference solution scheme for realistic freeways and urban road networks. Figure 1 demonstrates three typical cell combinations of freeways, considering on/off-ramps.
Three typical cell combinations of freeways: (a) with no ramps; (b) with one on-ramp; (c) with one off-ramp.
For the basic segments without on/off-ramps demonstrated in Figure 1(a), the conservation equation can be derived from LWR model as follows:(2)ρit+1=ρit+Tsλiliqi-1,it-qi,i+1t=ρit+Tsλiliqi,int-qi,outt,where i denotes the cell number, Ts is the time interval, li is the length of cell i, λi is the lane number of cell i, ρit means the vehicle density of cell i at the time of t. And qi,int and qi,outt, respectively, denote the traffic flow entering and leaving cell i during the tth time interval, Ts[t,t+1), including flows along the mainline and the on- and off-ramps.
In order to refine the transmission rules, it is assumed that two adjacent cells can be divided into one of five connection modes: (1) “free flow-free flow” (FF), (2) “congestion-congestion” (CC), (3) “congestion free flow” (CF), (4) “free flow congestion 1” (FC1), in which the wave front is moving downstream, and (5) “free flow congestion 2” (FC2), in which the wave front is moving downstream [5]. The SMM is composed of several linear models; straightforward linear techniques for model analysis and control design can be applied to the individual linear subsystems. The wave propagation direction is usually difficult to detect exactly. So, four distinguishable states are formulated in this paper as follows:
According to the conservation law at the boundary of cell i-1 and cell i, we get(4)qi,int=qi-1,outt+Rit,where Rit is the on-ramp flow and ρJ,it is the jam density of cell i.
In this paper, we extend the above simple state switch transmission rule to depict more complex cloverleaf junctions. As shown in Figure 2, each ramp is considered as a cell and the other types of junctions can be formulated in the simplified expression. The traffic flow transmission matrix can be denoted as the structure of (5). Here, we mainly consider the properties of cell ik+2(k=1,2,3,4) and ramp-cells by the following rules:(5)Cellikik+1ik+1ik+1ik+1ik0qik,ik+1ik+10qik,ik+1ik+20qik,ik+1ik+30qik,ik+1ik+40,k∈1,2,3,4,Cellik+1+2ik-1+4RrikRlikik00qik,Rrik0ik+2000qik+2,RlikRrik0qRrik,ik-1+400RlikqRlik,ik+1+2000,k∈1,2,3,4,i0=i4,i5=i1,where(6)qik+1,int=qik,ik+1t,qik+2,int=qik+1,ik+2t+qRlik-1,ik+2t,qik+3,int=qik+2,ik+3t,qik+4,int=qik+3,ik+4t+qRrik+1,ik+4t,qRrik,int=qik,Rrikt,qRlik,int=qik+2,Rlikt.Here, Rrik and Rlikrepresent the right/left turning ramp-cell, respectively.
Cell representation of a cloverleaf junction: (a) the layout of junction; (b) the corresponding division of cells; (c) the corresponding CTM cells.
2.2. Finite Difference of the Higher-Order Equation
In this paper, we applied the PDE model proposed in [7]. This PDE model replaces the density gradient term with the speed gradient term to consider drivers’ anticipation as follows:(7)dvx,tdt=∂v∂t+v∂v∂x=veρ-vT+w∂v∂x,where w represents the propagation speed of the disturbance. In this model, the characteristic speed is equal to or smaller than the macroscopic flow velocity. In reality, traffic density ρ(x,t) is hard to investigate or estimated directly (Figure 8); the most popular method is to monitor the time occupancy of vehicles passing detection stations. A number of literatures indicate that the time occupancy at x is proportional to traffic density around x.
Assuming there exist no merging or diverging in the section, say g(x,t)=0, and applying the finite difference method to discretize (1a) and (7), we obtain the following difference equations based on point detectors like inductive loops, microwave detectors, and radar sensors. The fundamental equation of traffic flow is(8)qit=λioitvitgit+ζiqt,where oit denotes the time occupancy at the upstream boundary of cell i, git represents the effective vehicle length which depends on both the actual vehicle length, which varies by lane and over the course of a day, and the loop’s electrical circuit, and ζiqt is the zero-mean white noise acting on the approximation to reflect the measurement inaccuracies and fluctuations in the traffic flow.
During one time interval, the drivers in cell i will adapt their speed to the local traffic density and are thus impacted by the traffic state in front of them. We suppose the adaptive behavior is influenced by the composition of traffic density of these two adjacent cells; it is a weighted average of the density in cell i and in cell i+1: (9)oi¯t=foαit,oit,oi+1t.
To simplify the combinational equations, we assume that this adaptation is linear as(10)oi¯t=αioit+1-αioi+1t,where the parameter 0<αi<1 weights how far ahead the drivers look and how well they anticipate the traffic density downstream. The drivers adapt the local traffic density more habitually; αi is closer to 1.
Since each drivers’ aggressiveness is random, the anticipant average space speed of cell i is relative to the speed of cell i and cell i+1; the weighted factor 0<βi<1 expresses the relative impact on the calculation of the speed of convective and inertial behavior:(11)vi¯t=fvβit,vit,vi+1t.
Similarly, we assume that(12)vi¯t=βivit+1-βivi+1t.
Because most of drivers tend to adapt their behavior whose function reflects an deterministic statistical relationship between the traffic flow and density or occupancy; the anticipant traffic flow yi,out(t) that passes through the boundary of cell i and cell i+1 is approximately equal to oi(t)vi¯(t), when the stationary equivalent traffic flow is taken into consideration.
Meanwhile, the transmission flow is also constrained by the current sending capacity of cell i, that is, oi(t)li/Ts, and by the receiving capacity of cell i+1, that is, [om-oi(t)]li+1/Ts, where om is the maximum occupancy; then it yieldsAyi,outt=minoitliTs,om-oi+1tli+1Ts,γioitvi¯t+1-γiQeoi¯t+ζiyt,where Qe[oi¯(t)] is the stationary equivalent traffic flow and time occupancy relationship function, Qe[oi¯(t)]=Qe1[oi¯(t)], if oi¯(t)<oc,i(t); Qe[oi¯(t)]=Qe1[oi¯(t)], otherwise. oc,i(t) is the critical time occupancy of cell i, and the coefficient 0<γi<1 expresses the relative impact on the calculation of the equivalent traffic flow of convective and macroscopic stable behavior and ζiyt is a zero-mean white noise that describes the uncertainty of traffic transmission between two adjacent cells.
Thus, we can formulate the evolution of time occupancy under the assumption that gi(t) is approximately equal to gi+1(t) by the following finite difference equation:(13)oit+1=oit+Tsλiliyi,int-yi,outt,where yi,int is equal to yi,outt when there exist no on/off-ramps. The formulas of conditions with merging or diverging areas can be obtained similarly as(14)vit+1=vit+Tsliwi-vitvi+1t-vit+TsTVeoit-vit+ζivt,ifvit<wivit+Tsliwi-vitvit-vi-1t+TsTVeoit-vit+ζivt,ifvit≥wi,where wi represents the propagation speed of the disturbance that is the backward speed of traffic shockwave and ζivt is a zero-mean white noise reflecting the fluctuations in the drivers’ speed.
3. Simulation Results
In the following experiments, we apply the initial stationary equivalent traffic speed and occupancy relationship function as (15) and the corresponding traffic flow and time occupancy relationship function as (16) given in [10]:(15)Veoi¯t=vfexp-1αoi¯toc,itα,(16)Qeoi¯t=oi¯tVeoi¯t=vfoi¯texp-1αmoi¯toc,itαm.
In the simulation of traffic density and speed evolution of a segment of urban freeways with periodical boundary conditions, the following constants that keep the same in all cells are applied. Geometric conditions and temporal-spatial discretization are (17)L=10km,ΔL=0.1km,Ts=1s,T=10s,gx,t=6.67m.Coefficients of fundamental diagram are(18)vf=80km/h,w=30km/h,oc,it=20%,Qm=1800veh/h/lane.Parameters of drivers’ adaptive behavior and standard deviations of Gaussian noise are(19)αit=0.5,βit=0.5,γit=0.5,σζivt=σζiqt=0.
Generally, the adaptive time can be chosen as 0 to 1. For numerical tests shown in this section, we just choose adaptive coefficients as 0.5 to show the influence of drivers’ adaptive behavior.
3.1. Experiment 1: Linear Initialization of Density
In Experiment 1, we test the influence of linear initialization of density (Figure 5). During the whole process of simulation, lane number and speed limit (equal to free-flow speed vf) are not changed. Four typical experiments are conducted with the same critical occupancy of 20%; the boundary occupancy conditions are [15%, 25%], [10%, 30%], [5%, 35%], and [0%, 40%], respectively. Therefore, the initial differences between critical occupancy are 5%, 10%, 15%, and 20%.
The results show that a larger initial swing leads to a more obvious backward moving jam that is observed in Figure 3. Waves moving in the opposite directions meet at around 3.5 min and evolve into the backward moving jam with a larger density and a smaller fluctuation that propagates downstream especially obtained in the last experiment. Similarly, the evolution pattern of traffic speed shows the wave fluctuation phenomena of a moving jam in Figure 4.
The initial and finial spatial distributions of occupancy for Experiment 1.
Occupancy evolution under linear initialization for Experiment 1.
Speed evolution under linear initialization for Experiment 1.
When the initial density swing is small, the two reverse waves meet and the gap arises but no moving jam is achieved because of the lack of enough density difference between adjacent cells where the fluctuations cross. The scatter plot of fundamental diagram is illustrated in Figure 6; the widest line represents the equivalent functions of traffic speed/occupancy and traffic flow/occupancy with exponential expression showed in (16). Points of four different colors are four experimental results, respectively. It indicates that a linear initialization with a larger swing is more nonstable and the results scatter in a larger area in the fundamental diagram. The numerical results indicate that a larger initial difference of density generally triggers a larger fluctuation of traffic flow.
Fundamental diagram obtained for Experiment 1.
3.2. Experiment 2: Local Cluster Effect
The local cluster effect corresponds to the stop-and-go phenomena observed in the field due to a small disturbance [7]. In this subsection, we simulate the local cluster effect with respect to a localized perturbation in an initial homogeneous condition.
The following initial variation of the average occupancy o0 is used as in [17]:(20)ox,0=o0+Δo0cosh-2160Lx-3L8-14cosh-240Lx-13L32,where v(x,0)=Ve[o(x,0)]; the other parameters are the same as the simulation of linear initialization and the second term acts as a local perturbation to the unstable constant distribution; then the combinatorial difference methods are used to show the evolution pattern when o(x,0) increases.
Figure 7 shows the initial local perturbations with Δo0=5% eventually evolve into local cluster effects when initial distribution of average time occupancy is limited in a specific extent and periodic boundary conditions are applied. The observed local cluster effects are also consistent with the diverse nonlinear dynamics phenomena observed in realistic traffic flow.
The initial and finial spatial distributions of occupancy for Experiment 2.
Traffic density evolution under small local disturbance for Experiment 2.
4. Model Validation and Application
The real data testing was conducted with traffic measurements collected from a 4.7 km, 3-lane stretch of Queen Elizabeth Way, Canada. As shown in Figure 9, it involves two on-ramps and an off-ramp at stations 50 and 53; another five loop detector stations were installed along this freeway and three at the on/off-ramps. The test stretch is subdivided into 8 segments, each with an approximate length of 650 m. Furthermore, a more refined cell division (where the size is only about 125 m) is conducted to explore the detailed spatial-temporal evolution of traffic jam.
The layout and the corresponding cell representation of Queen Elizabeth Way, Canada.
The individual vehicle data were recorded by the detectors from 6 am to 10 am on December 15, 1998. All data were converted into aggregated traffic measurements of flow and space mean speed with a temporal-resolution-level 20 seconds. Only traffic flow and occupancy data are used for validation.
For the sake of simple and convenient computation, we use piecewise linear assumption of fundamental diagram in this section: A-1Qeox,t=vfx,tox,tifox,t<ocx,t;A-2Qeox,t=wx,toJx,t-ox,t,otherwise.
Figure 10 shows the fundamental diagram defined in this paper, where points A and B are the critical milestone of the formation and dissipation of traffic congestion; the corresponding occupancies are kq and kd where “q” means queuing and “d” means dissipating. Regressive results of the real traffic data show the occupancy gap is accepted that midpoint of line AB is the inflexion of traffic states. As a result, the assumption of critical occupancy is approximately to kc. So, we determine that we meet a congested state if and only if the o(x,t)>oc(orkc) and q(x,t)/o(x,t)<[qq(x,t)-qd(x,t)]/[oq(x,t)-od(x,t)].
The fundamental diagram defined to depict asymmetric formation and dissipation of traffic congestion.
It is supposed that the traffic data collected by station 49~52 is unknown, and only date from stations 48 and 53 and the on/off-ramps is used to validate the finite different method discussed in this paper. The comparative results of measured and simulated occupancy data are shown in Figure 11. The agreement between two kinds of data suggests that this finite difference method is well applied to reflect the traffic congestion and the influences of on/off-ramp.
Spatial-temporal evolution of empirical and simulated occupancy.
Tests also show that the simulated results will significantly deviate from the empirical measurements, if no drivers’ anticipation term is introduced in the traffic flow model.
5. Conclusive Remarks
An improved spatial-temporal finite difference method based on LWR and a higher-order car-following model is proposed. A compositional model combining cell transmission model (CTM) and drivers’ adaptive car-following behaviors is proposed for the spatial-temporal evolution of traffic flows on freeways.
Numerical tests verify that the model is able to simulate complex traffic phenomena observed in the field such as shock waves, stop-and-go waves, and local cluster effects. Empirical traffic data are also used to validate the finite difference method proposed in this paper. The results are consistent with the nonlinear dynamic properties measured in practice.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
The supports of the National Natural Science Foundation of China (71303269), the Fundamental Research Funds for the Central Universities, Central University of Finance and Economics Major Research Task of Fostering Project (14ZZD006) are gratefully acknowledged.
LighthillM. J.WhithamG. B.On kinematic waves. II. A theory of traffic flow on long crowded roads19552291178317345RichardsP. I.Shock waves on the highway19564425110.1287/opre.4.1.42MR0075522DaganzoC. F.The cell transmission model: a dynamic representation of highway traffic consistent with the hydrodynamic theory199428426928710.1016/0191-2615(94)90002-72-s2.0-0028481056DaganzoC. F.A finite difference approximation of the kinematic wave model of traffic flow199529426127610.1016/0191-2615(95)00004-W2-s2.0-0029416696HorowitzR.SunX. T.MuñozL.SkabardonisA.VaraiyaP.ZhangM.MaJ. T.Design, field implementation and evaluation of adaptive ramp metering algorithms: final report2006UCB-ITS-PRR-2006-21LoH. K.ChangE.ChanY. C.Dynamic network traffic control200135872174410.1016/s0965-8564(00)00014-82-s2.0-0035450136JiangR.WuQ.-S.ZhuZ.-J.A new continuum model for traffic flow and numerical tests200236540541910.1016/S0191-2615(01)00010-82-s2.0-0036605010PayneH. J.BekeyG. A.Models of freeway traffic and control197115161Simulation Council Proceedings SeriesWillskyA. S.ChowE. Y.GershwinS. B.GreeneC. S.KurkjianA. L.HouptP. K.Dynamic model-based techniques for the detection of incidents on freeways19802533473602-s2.0-0019031682PapageorgiouM.BlossevilleJ.-M.Hadj-SalemH.Modelling and real-time control of traffic flow on the southern part of Boulevard Périphérique in Paris: Part I: modelling19902453453592-s2.0-0025570776KotsialosA.PapageorgiouM.DiakakiC.PavlisY.MiddelhamF.Traffic flow modeling of large-scale motorway networks using the macroscopic modeling tool METANET20023428229210.1109/tits.2002.8068042-s2.0-0348002067WangY.PapageorgiouM.MessmerA.Real-time freeway traffic state estimation based on extended Kalman filter: adaptive capabilities and real data testing200842101340135810.1016/j.tra.2008.06.0012-s2.0-54049147670AlecsandruC. D.2006Baton Rouge, La, USALouisiana State UniversityLuY. D.WongS. C.ZhangM. P.ShuC.-W.ChenW. Q.Explicit construction of entropy solutions for the Lighthill-Whitham-Richards traffic flow model with a piecewise quadratic flow–density relationship200842435537210.1016/j.trb.2007.08.0042-s2.0-39749105993PanT. L.SumaleeA.ZhongR. X.Indra-PayoongN.Short-term traffic state prediction based on temporal-spatial correlation20131431242125410.1109/tits.2013.22589162-s2.0-84883802662ZhongR. X.SumaleeA.PanT. L.LamW. H. K.Stochastic cell transmission model for traffic network with demand and supply uncertainties20139756760210.1080/18128602.2011.6345562-s2.0-84883498272HerrmannM.KernerB. S.Local cluster effect in different traffic flow models19982551-2163188