# An Optimization Method of Time Window Based on Travel Time and Reliability

Authors: Fengjie Fu; Dongfang Ma; Dianhai Wang; Wei Qian
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

The dynamic change of urban road travel time was analyzed using video image detector data, and it showed cyclic variation, so the signal cycle length at the upstream intersection was conducted as the basic unit of time window; there was some evidence of bimodality in the actual travel time distributions; therefore, the fitting parameters of the travel time bimodal distribution were estimated using the EM algorithm. Then the weighted average value of the two means was indicated as the travel time estimation value, and the Modified Buffer Time Index (MBIT) was expressed as travel time variability; based on the characteristics of travel time change and MBIT along with different time windows, the time window was optimized dynamically for minimum MBIT, requiring that the travel time change be lower than the threshold value and traffic incidents can be detected real time; finally, travel times on Shandong Road in Qingdao were estimated every 10 s, 120 s, optimal time windows, and 480 s and the comparisons demonstrated that travel time estimation in optimal time windows can exactly and steadily reflect the real-time traffic. It verifies the effectiveness of the optimization method.


---

## Body

## 1. Introduction

The time window determines the period of time that should be considered when estimating the current traffic information [1]. The selection of a time window meant that the traffic information, such as travel time, would be reestimated every time window. Most travel time estimation or prediction method adopted a constant time window subjectively: the time window was set as 2 min in the TransGuide algorithm [2]; TranStar algorithm estimated travel times every 30 s; Transmit algorithm used a time window of 15 min [3]; Ma and Koutsopoulos demonstrated that a larger time window gave smoother travel time with longer delay and a smaller one reduced delay with volatile travel time [4]. Sun and Zhang divided the travel time prediction model into two components based on the corresponding signal state: a green phase model with three different time windows and a red phase model with the time window equal to red time. Nevertheless, the values of the four time windows merely depended on the cycle length and were not the optimal [5]. Jiang et al. presented an optimization method of the time window of floating car data for the minimum relative errors and delayed rates of ILTT (Individual Link Travel Time) [6]. However, the method determined the optimal value simply based on the relationship between the quality of the traffic information and the time window. Therefore, it is necessary to propose optimizing indexes and an optimization algorithm to get the optimal time window for travel time estimation and prediction.Jackson and Jucker found that there was a trade-off between the result of the travel time estimation or prediction and travel time reliability [7]. Many studies defined travel time reliability on the basis of the mean, standard deviation, and percentiles [8–10]. For example, the Federal Highway Administration (FHWA) presented the five measures of travel time reliability: (1) the 90th or the 95th percentile travel time; (2) the travel time index (TTI), defined as the ratio of the mean to the free-flow travel time; (3) the Buffer Time Index (BTI), namely, the difference between the 95th percentile and the mean divided by the mean; it indicates the amount of extra time needed to be on time for 95% of trips; (4) the planning time index (PTI), namely, the ratio of the 95th percentile travel time to the free-flow travel time; (5) congestion frequency, defined as the ratio of congestion period to the total days [11]. But Van Lint and Van Zuylen indicated that indexes based on the mean and standard deviation may not reflect travel time reliability well since the travel time distribution showed significant skewness. So they proposed the so-called skew-width method composed of the skewness of the travel time and the width of travel time [12]. And Emam and Al-Deek pointed out that log-normal distribution emerged as a useful model representing travel time by examining the use of different statistical distributions including log-normal, exponential, gamma, and Weibull distributions for empirical travel time data [13]. However, many research found that there was some evidence of bimodality in the actual travel time distributions and developed many models for estimating the bimodal distribution [14–16]. Yang et al. demonstrated that a constant bimodal distribution could not always provide good fit for travel times of interrupted flow on different urban roads or during different periods because of the traffic state on the road, the delay at the downstream intersection, peak hours, and so on. They modeled six types of bimodal distributions which applied normal and lognormal distribution as typical distribution to fit RFID (Ratio Frequency Identification) data in Nanjing. The distribution could distinguish the travel time samples of the rapid flow and slow flow and provided the foundation for redefinition of travel time reliability [17]. However there exists inconsistency in the mixture model of normal and lognormal distribution because of the different units. Then the use of different dimensions in the data in the model can lead to different results [18].The objective of this research is to develop a method to optimize the value of time window on the basis of travel time and reliability. The remainder of the paper is organized as follows: the second section outlines the traffic data source and travel time calculation; the third section proposes a new index reflecting travel time reliability in accord with the bimodal distribution; in the fourth section, an optimization method is presented to obtain the optimal time window; finally, the fifth section gives numerical experiments to demonstrate effectiveness of the optimization method and the optimal time window.


## 2. Study Data

### 2.1. Data Collection System

The Multifunction Automatic Detecting and Recording System is installed at intersections with video pattern and capture pattern, so it can not only detect and capture illegal driving but also monitor and record all passing vehicles. The front-end devices transfer videos and images to the back-end computer; and then traffic data, including plate number, intersection ID, departure time, departure lane ID, departure entrance ID, departure speed, and position obtained by using ANNR, is available. The collected traffic information are stored as database, pictures, and videos. The data collection system becomes wide and the installation is as shown in Figure1.Figure 1
Installations of the Multifunction Automatic Detecting and Recording System.


### 2.2. Data Precision

Travel time information data used in the paper is collected from cameras at intersections on the Shandong Road in Qingdao, from December 1, 2014, to December 31, 2014. Firstly, traffic data was collected at Shandong-Jiangxi Intersection, including database and videos, for analyzing the practical detection precision of the system. Figure2(a) shows time-of-day fluctuations of the number of identified vehicles (IV), passing traffic volumes (PV), and their average percentage during 7:00–19:00 on December 22, 2014, at the southern approach. The number of identified vehicles (i.e., vehicles captured and recorded and license plate positioned, including ones with license plate numbers incorrectly recognized) is calculated from the database and passing traffic volume is observed from videos by manual counting. A typical time-of-day change with two peaks can be found in the passing traffic volumes and the number of identified vehicles. The great fluctuations between time windows are due to the signal control. The average percentage of IV with PV without considering the errors of the identified plate number is about 93.28%. The accuracy of plate number can be larger than 90% once identified and the software makers have verified it. Moreover, it went dark after 5:30 and the percentage is still high with the flash or everlasting lamp working as shown in Figure 2(b). Therefore, the author can demonstrate that the real average sample rate can be up to 85.8%, and it can be guaranteed at night with the flash or everlasting lamp working.Figure 2
(a) The detection precision result at Shandong-Jiangxi Intersection during 7:00–19:00 on December 22, 2014 and (b) pictures captured by the system at 18:36.
(a)
(b)


### 2.3. Travel Time Acquisition

The travel time of a vehicle on a link, namely, the time difference from the upstream (A) to downstream (B) intersection, is computed as shown in (1) by matching the plate number from a pair of cameras A and B:(1)Ti=tiB-tiA,where tiA and tiB are the time that vehicle i passes the stop line at the upstream and downstream intersection, respectively and Ti is the travel time of vehicle i on the link.


## 2.1. Data Collection System

The Multifunction Automatic Detecting and Recording System is installed at intersections with video pattern and capture pattern, so it can not only detect and capture illegal driving but also monitor and record all passing vehicles. The front-end devices transfer videos and images to the back-end computer; and then traffic data, including plate number, intersection ID, departure time, departure lane ID, departure entrance ID, departure speed, and position obtained by using ANNR, is available. The collected traffic information are stored as database, pictures, and videos. The data collection system becomes wide and the installation is as shown in Figure1.Figure 1
Installations of the Multifunction Automatic Detecting and Recording System.


## 2.2. Data Precision

Travel time information data used in the paper is collected from cameras at intersections on the Shandong Road in Qingdao, from December 1, 2014, to December 31, 2014. Firstly, traffic data was collected at Shandong-Jiangxi Intersection, including database and videos, for analyzing the practical detection precision of the system. Figure2(a) shows time-of-day fluctuations of the number of identified vehicles (IV), passing traffic volumes (PV), and their average percentage during 7:00–19:00 on December 22, 2014, at the southern approach. The number of identified vehicles (i.e., vehicles captured and recorded and license plate positioned, including ones with license plate numbers incorrectly recognized) is calculated from the database and passing traffic volume is observed from videos by manual counting. A typical time-of-day change with two peaks can be found in the passing traffic volumes and the number of identified vehicles. The great fluctuations between time windows are due to the signal control. The average percentage of IV with PV without considering the errors of the identified plate number is about 93.28%. The accuracy of plate number can be larger than 90% once identified and the software makers have verified it. Moreover, it went dark after 5:30 and the percentage is still high with the flash or everlasting lamp working as shown in Figure 2(b). Therefore, the author can demonstrate that the real average sample rate can be up to 85.8%, and it can be guaranteed at night with the flash or everlasting lamp working.Figure 2
(a) The detection precision result at Shandong-Jiangxi Intersection during 7:00–19:00 on December 22, 2014 and (b) pictures captured by the system at 18:36.
(a)
(b)


## 2.3. Travel Time Acquisition

The travel time of a vehicle on a link, namely, the time difference from the upstream (A) to downstream (B) intersection, is computed as shown in (1) by matching the plate number from a pair of cameras A and B:(1)Ti=tiB-tiA,where tiA and tiB are the time that vehicle i passes the stop line at the upstream and downstream intersection, respectively and Ti is the travel time of vehicle i on the link.


## 3. Travel Time Reliability

### 3.1. The Basic Unit of Time Window

Points in Figures3(a) and 3(b) indicate travel times and departure times of all vehicles during the four different periods (22:00–23:00, 7:30–8:30, 13:00–14:00, and 17:30–18:30) of day on the Shandong Road from Shandong-Jiangxi Intersection to Shandong-Minjiang Intersection and from Shandong-Minjiang Intersection to Shandong-Jiangxi Intersection separately. It was obvious that traffic flow was interrupted and travel times showed cyclic change due to traffic state, signal timing, and signal cycle length. Therefore, the signal cycle length at the upstream intersection should be set as the basic unit of the time window for travel time estimation and prediction.Figure 3
(a) Individual travel time along with time from Shandong-Jiangxi Intersection to Shandong-Minjiang Intersection and (b) individual travel time along with time from Shandong-Minjiang Intersection to Shandong-Jiangxi Intersection.
(a)
(b)


### 3.2. Travel Time Distribution

It can be seen that points surrounded by a red outline in Figure3 were composed of one cluster or multiclusters: (1) vehicles can pass through the downstream intersection without signal delay or wait for at most a short duration of the red phase during 22:00–23:00 and 13:00–14:00 because of unsaturated traffic condition and short cycles. Therefore, travel time data fell into no more than two clusters and each cluster was centered. (2) Some vehicles may wait for two or more red lights during 7:30–8:30 and 17:30–18:30 because of the great traffic demand and oversaturated traffic state. Consequently, travel time points can be divided into two or more groups and each group was discrete. A bimodal distribution is fitted against the travel time data in each cycle. The model of the bimodal distribution is as presented in(2)fTT=ω112πσ1exp⁡TT-μ122σ12+ω212πσ2exp⁡TT-μ222σ22,where f(TT) is the probability density function of travel time; TT indicates the estimation value of travel time; μ expresses the mean of the travel time distribution and μ1 is the larger one; σ2 is the variance of travel time with σ1 and σ2 corresponding to μ1 and μ2, respectively; likewise, ω1 and ω2, namely, the ratio of each subsample to the whole sample, correspond to μ1 and μ2 separately with the sum of them equal to 1.Parameters in (2) are estimated using the Expectation Maximization (EM) algorithm, and the process is as follows:(1)
Initialize the parametersθ(μ,σ,ω) of the Gaussian mixture distribution model.(2)
E-step: given the current estimate of the parameters in step(1), the conditional distributions of the two distributions p(i,1) and p(i,2) are determined by Bayes theorem to using their own probability density functions g(TTi∣μ1,σ1) and g(TTi∣μ2,σ2) weighted by ω1 and ω2 (see (3)).(3)
The M-step computes parameters maximizing the expected log-likelihood found on the E-step. The estimation results of the first distribution withm samples are presented as shown in (4):(3)pi,1=ω1gTTi∣μ1,σ1ω1gTTi∣μ1,σ1+ω2gTTi∣μ2,σ2,pi,2=ω2gTTi∣μ2,σ2ω1gTTi∣μ1,σ1+ω2gTTi∣μ2,σ2,(4)μ1=1ω1∑i=1mpi,1TTi,σ12=1ω1∑i=1mpi,1TTi-μ1TTi-μ1T,ω1=∑i=1mpi,1.(4)
Repeat from step(1) again until the value of θN-θN-1 meets the threshold.Consequently, the bimodal distributions of travel times in the certain signal cycles surrounded by a red outline in Figure3 were got using the EM algorithm. Figure 4 gives their frequency histogram and probability density curve. The results are tested by K-S (Kolmogorov-Smirnov) test. The null hypothesis is that the samples have a standard normal distribution. The alternative hypothesis is that the samples do not have that distribution. The test accepts the null hypothesis if the significance probabilityProb is larger than the significance level and rejects it otherwise. As Table 1 shows, the values ofProb in the K-S test were all larger than 0.05. It indicates that the samples do have a standard normal distribution at 0.05 significance level and the bimodal distribution results fit well.Table 1
The results of K-S test.
Data
Figure4(a)-1
Figure4(a)-2
Figure4(a)-3
Figure4(a)-4
Figure4(b)-1
Figure4(b)-2
Figure4(b)-3
Figure4(b)-4
N
28
154
77
48
58
100
74
225
Prob
0.11
0.39
0.17
0.21
0.23
0.18
0.30
0.12Figure 4
(a) Fitting results in a single cycle marked in Figure3(a) (from Shandong-Jiangxi Intersection to Shandong-Minjiang Intersection) and (b) fitting results in a single cycle marked in Figure 3(b) (from Shandong-Minjiang Intersection to Shandong-Jiangxi Intersection).
(a)
(b)


### 3.3. Travel Time Reliability Index

A new travel time reliability index, the Modified Buffer Time Index, is proposed on the basis of the characteristic of the bimodal distribution and its common ground with the unimodal distribution as follows:(1)
When travel times are bimodal-distributed, the travel time estimation andMBIT are computed as presented in(5)TT=ω1μ1+ω2μ2,MBIT=μ2+σ2-μ1-σ1TTp1=p2,whereMBIT is the value of the Modified Buffer Time Index and TTp1=p2 indicates the travel time estimation at the crossing point of the two distributions, namely, the estimation value when p1=p2 (see Figure 5).(2)
When travel times are unimodal-distributed, the computations of travel time estimation andMBIT are as follows:(6)TT=μ,MBIT=2σμ.Figure 5
Parameters of travel time distribution.
(a)
(b)


## 3.1. The Basic Unit of Time Window

Points in Figures3(a) and 3(b) indicate travel times and departure times of all vehicles during the four different periods (22:00–23:00, 7:30–8:30, 13:00–14:00, and 17:30–18:30) of day on the Shandong Road from Shandong-Jiangxi Intersection to Shandong-Minjiang Intersection and from Shandong-Minjiang Intersection to Shandong-Jiangxi Intersection separately. It was obvious that traffic flow was interrupted and travel times showed cyclic change due to traffic state, signal timing, and signal cycle length. Therefore, the signal cycle length at the upstream intersection should be set as the basic unit of the time window for travel time estimation and prediction.Figure 3
(a) Individual travel time along with time from Shandong-Jiangxi Intersection to Shandong-Minjiang Intersection and (b) individual travel time along with time from Shandong-Minjiang Intersection to Shandong-Jiangxi Intersection.
(a)
(b)


## 3.2. Travel Time Distribution

It can be seen that points surrounded by a red outline in Figure3 were composed of one cluster or multiclusters: (1) vehicles can pass through the downstream intersection without signal delay or wait for at most a short duration of the red phase during 22:00–23:00 and 13:00–14:00 because of unsaturated traffic condition and short cycles. Therefore, travel time data fell into no more than two clusters and each cluster was centered. (2) Some vehicles may wait for two or more red lights during 7:30–8:30 and 17:30–18:30 because of the great traffic demand and oversaturated traffic state. Consequently, travel time points can be divided into two or more groups and each group was discrete. A bimodal distribution is fitted against the travel time data in each cycle. The model of the bimodal distribution is as presented in(2)fTT=ω112πσ1exp⁡TT-μ122σ12+ω212πσ2exp⁡TT-μ222σ22,where f(TT) is the probability density function of travel time; TT indicates the estimation value of travel time; μ expresses the mean of the travel time distribution and μ1 is the larger one; σ2 is the variance of travel time with σ1 and σ2 corresponding to μ1 and μ2, respectively; likewise, ω1 and ω2, namely, the ratio of each subsample to the whole sample, correspond to μ1 and μ2 separately with the sum of them equal to 1.Parameters in (2) are estimated using the Expectation Maximization (EM) algorithm, and the process is as follows:(1)
Initialize the parametersθ(μ,σ,ω) of the Gaussian mixture distribution model.(2)
E-step: given the current estimate of the parameters in step(1), the conditional distributions of the two distributions p(i,1) and p(i,2) are determined by Bayes theorem to using their own probability density functions g(TTi∣μ1,σ1) and g(TTi∣μ2,σ2) weighted by ω1 and ω2 (see (3)).(3)
The M-step computes parameters maximizing the expected log-likelihood found on the E-step. The estimation results of the first distribution withm samples are presented as shown in (4):(3)pi,1=ω1gTTi∣μ1,σ1ω1gTTi∣μ1,σ1+ω2gTTi∣μ2,σ2,pi,2=ω2gTTi∣μ2,σ2ω1gTTi∣μ1,σ1+ω2gTTi∣μ2,σ2,(4)μ1=1ω1∑i=1mpi,1TTi,σ12=1ω1∑i=1mpi,1TTi-μ1TTi-μ1T,ω1=∑i=1mpi,1.(4)
Repeat from step(1) again until the value of θN-θN-1 meets the threshold.Consequently, the bimodal distributions of travel times in the certain signal cycles surrounded by a red outline in Figure3 were got using the EM algorithm. Figure 4 gives their frequency histogram and probability density curve. The results are tested by K-S (Kolmogorov-Smirnov) test. The null hypothesis is that the samples have a standard normal distribution. The alternative hypothesis is that the samples do not have that distribution. The test accepts the null hypothesis if the significance probabilityProb is larger than the significance level and rejects it otherwise. As Table 1 shows, the values ofProb in the K-S test were all larger than 0.05. It indicates that the samples do have a standard normal distribution at 0.05 significance level and the bimodal distribution results fit well.Table 1
The results of K-S test.
Data
Figure4(a)-1
Figure4(a)-2
Figure4(a)-3
Figure4(a)-4
Figure4(b)-1
Figure4(b)-2
Figure4(b)-3
Figure4(b)-4
N
28
154
77
48
58
100
74
225
Prob
0.11
0.39
0.17
0.21
0.23
0.18
0.30
0.12Figure 4
(a) Fitting results in a single cycle marked in Figure3(a) (from Shandong-Jiangxi Intersection to Shandong-Minjiang Intersection) and (b) fitting results in a single cycle marked in Figure 3(b) (from Shandong-Minjiang Intersection to Shandong-Jiangxi Intersection).
(a)
(b)


## 3.3. Travel Time Reliability Index

A new travel time reliability index, the Modified Buffer Time Index, is proposed on the basis of the characteristic of the bimodal distribution and its common ground with the unimodal distribution as follows:(1)
When travel times are bimodal-distributed, the travel time estimation andMBIT are computed as presented in(5)TT=ω1μ1+ω2μ2,MBIT=μ2+σ2-μ1-σ1TTp1=p2,whereMBIT is the value of the Modified Buffer Time Index and TTp1=p2 indicates the travel time estimation at the crossing point of the two distributions, namely, the estimation value when p1=p2 (see Figure 5).(2)
When travel times are unimodal-distributed, the computations of travel time estimation andMBIT are as follows:(6)TT=μ,MBIT=2σμ.Figure 5
Parameters of travel time distribution.
(a)
(b)


## 4. The Optimization of the Time Window

### 4.1. The Optimization Index

As (7) presents, the value of kth time window, TWk, is the sum of several cycles:(7)TWk=∑j=1mkcj,where mk expresses the number of signal cycles in kth time window and cj expresses the length of jth signal cycle with j smaller than mk; when the signal timing plan is fixed and the cycle length is c, TWk is equal to the product of mk and c.Travel time information data during 17:30–18:30 on the link from Shandong-Jiangxi Intersection to Shandong-Minjiang Intersection were analyzed to determine the optimization index. As Figure6 shows, the values of MBIT surrounded by the red line reduced with the number of signal cycles in the first time window (m1), and then traffic demand changed greatly resulting in a great change in MBIT and TT (surrounded by the blue line).Figure 6
The variation ofMBIT and TT with m1.Overall, when traffic demand changes greatly, traffic flow runs unstable and travel times change greatly, so a larger time window reduces the travel time reliability instead and cannot reflect the real change of the traffic condition. Although traffic flow runs stable with stable travel times when the traffic demand change little, a larger time window cannot ensure timely response to traffic events. However, the two conditions have one thing in common: the travel time reliability in a small time window can be too low. Therefore, the optimization of the time window should be based on the travel time and its reliability. And it is realized by using the relative variation ratios of the travel time estimation and Modified Buffer Time IndexRTTmk and RMBITmk as presented in(8)RTTmk=TTmk+1-TTk1TTk1·100%,RMBITmk=MBITmk+1-MBITk1MBITk1·100%,where TTk1, TTmk+1 are travel time estimations of the first signal cycle and (mk+1)th signal cycle in kth time window separately and MBITk1 and MBITmk+1, respectively, indicate the Modified Buffer Time Index corresponding to TTk1 and TTmk+1.Traffic events can be detected using the relative variation ratios of the travel time estimations in two successive time windows at least. So as (9) demonstrates, a timely response to traffic events requires the sum of the two successive time windows to be less than the minimum response time to traffic events:(9)Δtk=2∑j=1m1cj-Tk=1TWk-1+∑j=1mkcj-Tk>1,where Δtk is the time difference between the sum of the two successive time windows and the minimum response time to traffic events; T is the minimum response time to traffic events, and under urban roads it is taken between 15 min and 30 min [19]; TWk-1 indicates the value of (k-1)th time window.


### 4.2. The Optimization Process

The values ofRTTm1 and RMBITm1 when the number of signal cycles of the first time window (m1) ranges from 1 to 22 are as shown in Figure 7. On the one hand, RMBITm1 was less than 0 when m1≤9 and kept reducing when m1≤6, so m1 should not be larger than 6 to ensure a higher travel time reliability; on the other hand, RTTm1 did not exceed 5% (for instance) until m1>7 and Δt1 did not exceed 0 until m1>3, so m1 should not be larger than 3. In general, the optimization process on the basis of the variation of RTTmk and RMBITmk and the value of Δtk is as follows:(1)
The internalization:k=1, m1=1, RTTm1=0, RMBITm1=0.(2)
Calculate the values ofRTTmk, RMBITmk, and Δtk.(3)
The threshold judgment:
ifRTTmk≤ε&RMBITmk≤0&RMBITmk-RMBITmk-1≤0, mk=mk+1, and repeat from Step (2);ifRTTmk>ε∣RMBITmk>0∣Δtk>0∣RMBITmk-RMBITmk-1>0, mk=mk-1, k=k+1, and repeat from Step (1).Figure 7
The variation ofRTTm1, RMBITm1, and Δt1 with m1.The value ofε can be set according to traffic control requirements. It is taken as 5% in the paper and T=15 min. Then the optimal number of signal cycles in the first time window during 17:30–18:30 on the link from Shandong-Jiangxi Intersection to Shandong-Minjiang Intersection is equal to 3, so the optimal time window is 450 s.


## 4.1. The Optimization Index

As (7) presents, the value of kth time window, TWk, is the sum of several cycles:(7)TWk=∑j=1mkcj,where mk expresses the number of signal cycles in kth time window and cj expresses the length of jth signal cycle with j smaller than mk; when the signal timing plan is fixed and the cycle length is c, TWk is equal to the product of mk and c.Travel time information data during 17:30–18:30 on the link from Shandong-Jiangxi Intersection to Shandong-Minjiang Intersection were analyzed to determine the optimization index. As Figure6 shows, the values of MBIT surrounded by the red line reduced with the number of signal cycles in the first time window (m1), and then traffic demand changed greatly resulting in a great change in MBIT and TT (surrounded by the blue line).Figure 6
The variation ofMBIT and TT with m1.Overall, when traffic demand changes greatly, traffic flow runs unstable and travel times change greatly, so a larger time window reduces the travel time reliability instead and cannot reflect the real change of the traffic condition. Although traffic flow runs stable with stable travel times when the traffic demand change little, a larger time window cannot ensure timely response to traffic events. However, the two conditions have one thing in common: the travel time reliability in a small time window can be too low. Therefore, the optimization of the time window should be based on the travel time and its reliability. And it is realized by using the relative variation ratios of the travel time estimation and Modified Buffer Time IndexRTTmk and RMBITmk as presented in(8)RTTmk=TTmk+1-TTk1TTk1·100%,RMBITmk=MBITmk+1-MBITk1MBITk1·100%,where TTk1, TTmk+1 are travel time estimations of the first signal cycle and (mk+1)th signal cycle in kth time window separately and MBITk1 and MBITmk+1, respectively, indicate the Modified Buffer Time Index corresponding to TTk1 and TTmk+1.Traffic events can be detected using the relative variation ratios of the travel time estimations in two successive time windows at least. So as (9) demonstrates, a timely response to traffic events requires the sum of the two successive time windows to be less than the minimum response time to traffic events:(9)Δtk=2∑j=1m1cj-Tk=1TWk-1+∑j=1mkcj-Tk>1,where Δtk is the time difference between the sum of the two successive time windows and the minimum response time to traffic events; T is the minimum response time to traffic events, and under urban roads it is taken between 15 min and 30 min [19]; TWk-1 indicates the value of (k-1)th time window.


## 4.2. The Optimization Process

The values ofRTTm1 and RMBITm1 when the number of signal cycles of the first time window (m1) ranges from 1 to 22 are as shown in Figure 7. On the one hand, RMBITm1 was less than 0 when m1≤9 and kept reducing when m1≤6, so m1 should not be larger than 6 to ensure a higher travel time reliability; on the other hand, RTTm1 did not exceed 5% (for instance) until m1>7 and Δt1 did not exceed 0 until m1>3, so m1 should not be larger than 3. In general, the optimization process on the basis of the variation of RTTmk and RMBITmk and the value of Δtk is as follows:(1)
The internalization:k=1, m1=1, RTTm1=0, RMBITm1=0.(2)
Calculate the values ofRTTmk, RMBITmk, and Δtk.(3)
The threshold judgment:
ifRTTmk≤ε&RMBITmk≤0&RMBITmk-RMBITmk-1≤0, mk=mk+1, and repeat from Step (2);ifRTTmk>ε∣RMBITmk>0∣Δtk>0∣RMBITmk-RMBITmk-1>0, mk=mk-1, k=k+1, and repeat from Step (1).Figure 7
The variation ofRTTm1, RMBITm1, and Δt1 with m1.The value ofε can be set according to traffic control requirements. It is taken as 5% in the paper and T=15 min. Then the optimal number of signal cycles in the first time window during 17:30–18:30 on the link from Shandong-Jiangxi Intersection to Shandong-Minjiang Intersection is equal to 3, so the optimal time window is 450 s.


## 5. Empirical Analysis

The optimal time windows during 17:30–18:30 on the link from Shandong-Jiangxi Intersection to Shandong-Minjiang Intersection can be obtained using the optimization method in Section3 and the results are as presented in Table 2. Each column in the table gives the travel time reliability index values in each time window and the time window value. For example, in the second column, k=2 indicates that it is the optimization of the second time window and then mk should just be m2; m2 increases gradually and when it is equal to 2, Δtk, RTTmk, and RMBITmk-1 do not exceed their thresholds, while RMBITmk is over 10%, therefore m2 should be set as 2, namely, 450 s. Thus all the optimized time windows values can be got and they are shown as in the last two rows.Table 2
The optimal time intervals on Shandong Road during 17:30–18:30.
k
1
2
3
4
5
6
7
8
9
10
11
12
Δ
t
k
0
−300
−300
−300
−600
−600
0
0
−600
−600
−300
−600
R
T
T
m
k (%)
1.6
0.2
8.8
7.7
7.1
2.7
8.1
0.1
0.2
15.0
13
4.5
R
MBIT
m
k
-
1 (%)
−0.2
−3.3
−12.4
−22.1
0
0
−5.6
−34
0
0
−3.6
0
R
MBIT
m
k (%)
−0.5
10.4
−17.0
−22.3
−21
1.4
−5.8
−38
6.6
−18.2
−10.1
9.9
m
k
3
2
2
2
1
1
3
3
1
1
2
1
T
W
k
450
300
300
300
150
150
450
450
150
150
300
150Figure8 shows travel time estimations with the time window equal to 10 s and 120 s, the optimal values of TWk, and 480 s from top to bottom. Red points denote travel time sample data of all vehicles on the road and each black bar shows the travel time estimation value in each time window. Time windows in all the cases except for the third one were set as fixed values. It can be seen from the figure that the widths of black bars in each of these three cases are the same except the one using the optimal values of TWk. When the time window was set as 10 s or 120 s, travel time sample data in a single time window was so few that the travel time estimation was directly taken as the average travel time. Figure 8(a) presents that there existed gaps and great fluctuations among travel time estimations (TW=10 s) since traffic flow was interrupted due to traffic signal control and the time window was too small. The difference of the sample data among each time window was big and correspondingly several travel time estimations (TW=120 s) were too large or too small and cannot demonstrate the real value, since TW was not a multiple of the signal cycle length and the travel times changed periodically with the signal timing, while a time window of 480 s was so large that there was hardly any difference or change among travel time estimations and the result cannot make any sense to traffic state identification or traffic events detection. However, when the value of each time window was taken as the optimal one, the estimation result can reflect the change of the traffic condition with stable values as shown in Figure 8(c).Figure 8
Travel time estimations in different time intervals.
(a)
TW = 10 s
(b)
TW = 120 s
(c)
TW = The optimal values
(d)
TW = 480 s


## 6. Conclusion

The authors provided an optimization model for the determination of the time window: firstly, by exploring the characteristic of the travel time bimodal distribution, the model used the means, standard deviations, and crossing point value of the two distributions to compute the Modified Buffer Time Index indicating the travel time reliability on urban roads; secondly, the model took the signal cycle length as the basic unit of the time window, the relative variation ratios of the travel time estimation, and the Modified Buffer Time Index as the optimizing index for a minimum travel time reliability and timely response to traffic events; finally, the empirical analysis verified the optimization method since travel time estimations using the optimal time window can reflect the variation of the real traffic state better with much stable values.However, much more travel time information data should be collected and analyzed to confirm whether all the travel time distributions can be modeled by the six types of bimodal distributions applying normal and lognormal distribution. Correspondingly, the estimation of travel time and the computation of travel time reliability should be completed. Moreover, study on the relationship between the decision threshold and different urban roads should be added to realize the automatic determination of the decision threshold.



---

*Source: 921480-2015-10-12.xml*