Real-Time Arterial Coordination Control Based on Dynamic Intersection Turning Fractions Estimation Using Genetic Algorithm

Real-time arterial coordination control is crucial for urban transportation systems and is partially dependent on dynamic turning flows at intersections. Few existing researches employ such information due to the restrictions of traffic surveillance systems. This paper presents a model framework for real-time arterial coordination control based on dynamic intersection turning fraction estimation, including three submodels: (1) a parameter optimizationmodel to estimate dynamic intersection turning fractions using detected link counts at entering and exiting approaches; (2) a nonlinear model using minimum delay as an objective to optimize the time-varying public cycle for the arterial road based on the estimated turning flows; and (3) a revised optimization model to achieve real-time offset and split for the arterial road using the novel uninterrupted ratio as objective function. Two revised genetic algorithms are developed to solve the first and third submodels, respectively, and an ordinary optimization algorithm is designed for the second submodel. Time-varying public cycle, offset, and split constitute the real-time arterial coordination control scheme together. The general model framework removes most of the assumptions of conventional arterial control models and provides a time-varying timing plan. Simulation experiments using actual data indicate that the proposed model yields much better results than the existing methods.


Introduction
Real-time traffic signal control, especially arterial coordination control, is important for intelligent transportation systems (ITS), and the dynamic turning movement flows, that is, the dynamic origin-destination (O-D) flows at intersections, are valuable input data for signal control.However, under current traffic surveillance systems, the dynamic turning movement flows are impossible to be collected directly.Therefore, the dynamic turning fraction estimation from detected real-time link counts has been studied extensively during the past two decades.
Most researches about dynamic fractions estimation have constructed time-varying interrelations between intersection turning flows to be estimated and detected link counts at entering and exiting approaches and have proposed a series of optimization models, for example, Nihan and Davis (1987) [1], Nihan and Davis (1989) [2], Bell (1991) [3], and Jiao et al. (2005) [4].All these models were presented as parameter optimization formulations and were solved using traditional optimization or heuristic approaches.
To improve the estimation efficiency for on-line applications, some other works about dynamic O-D flows estimation have fallen within the scope of state-space methods and have formulated several efficient estimation models using Kalman filtering, for example, Ashok and Ben-Akiva (2002) [5], Bierlaire and Crittin (2004) [6], Lin and Chang (2007) [7] and Lou and Yin (2010) [8].All these models focused on dynamic O-D flows estimation for freeway corridors or general road networks with rather high efficiency.Of course, they can also be transformed to estimate dynamic turning fractions at intersections; for example, Jiao et al. (2014) [9] proposed a Bayesian combined model to estimate intersection fractions, integrating Kalman filtering and back propagation neural model running simultaneously.
With respect to the arterial coordination traffic signal control, mainly two groups of models have been developed, including maximum green wave band (MGWB, or 2 Mathematical Problems in Engineering MAXBAND) method [10,11] and minimum delay method [12].Most existing arterial coordination control models have been further developed based on the above two methods, including some real-world traffic signal control systems, such as Sydney Coordinated Adaptive Traffic System (SCATS) [13] and Split, Cycle and Offset Optimization Technique (SCOOT) [14].All these models and systems have contributed a great deal to urban traffic management systems.However, very few of them employed the valuable information of dynamic O-D flows or dynamic intersection turning fractions due to the restrictions of current traffic surveillance systems.Moreover, existing arterial control models have some limitations as follows: the distance between two adjacent intersections should be approximately equal, and not longer than 800 meters; arriving vehicles must follow some given distribution; the timing plan remains unchanged during some intervals; and so forth.Therefore, the revised model which can eliminate these limitations is quite necessary.This paper will develop a real-time arterial coordination control model with the estimated dynamic turning fractions as input data, and both the arterial control model and the turning fraction estimation model lend themselves to formulations as rather complex optimization models.Existing researches have proved that heuristic approaches or swarm intelligence algorithms are rather suitable for complicated optimization models, for example, genetic algorithm (GA) [15][16][17] and ant colony optimization algorithm (ACO) [18][19][20].This paper will also design two revised genetic algorithms to solve the dynamic turning fraction estimation model and real-time arterial control model, respectively.
The rest of this paper is organized in the following sections.Section 2 describes the basic problem, as well as the general model architecture.Section 3 presents the model framework for real-time arterial coordination control, including dynamic intersection turning fraction estimation model, real-time public signal cycle optimization model, and real-time offset and split optimization model.Section 4 designs three algorithms to solve the above three models, respectively, including two revised genetic algorithms and a classical optimization algorithm.Section 5 illustrates the results of a case study through simulation experiments using practical data.Section 6 concludes this paper and suggests some future research directions.

General Model Architecture
The general model architecture of the arterial coordination control is shown in Figure 1.
Figure 1 illustrates a typical arterial corridor, the layout scheme of detectors, and the flow of the arterial coordination control.Detectors are placed at both entering and exiting approaches of each isolated intersection to collect the entrance and exit link counts, which are fundamental input data for the whole model framework.The general model consists of the following three submodels: (1) dynamic turning fraction estimation model: to estimate dynamic turning fractions or turning flows based on detected link counts; (2) nonlinear signal control model for intersection: to achieve the optimized cycle length of each intersection based on the estimated intersection turning flows, with the maximum one as the public cycle length of the arterial road; (3) arterial coordination control model: to optimize the offset and split of the arterial road based on the achieved public cycle length.
All three models will be illustrated in detail in the following sections.

Real-Time Arterial Coordination
Control Model Without considering travel time to cross the intersection, similar estimation problems have been studied in several papers.As described in our previous work [4], to accommodate the possible outliers in the detected input data, a least absolute deviation (LAD) formulation is much more robust than traditional least square (LSQ) form [3] in the objective function.Therefore, the LAD model to estimate dynamic turning fractions at intersection level is formulated as follows: In ( 2), the turning fractions are assumed to be constant within continuous  intervals for each step to make the model overdetermined, along with a dynamic updating mechanism to incorporate the time-varying turning fractions.the delay at intersections, a nonlinear optimization model is formulated to optimize the signal cycle length: Based on the estimated dynamic turning flows in Section 3.2, VOC  can be achieved easily.Together with other detected parameters, (3) can be solved, and the optimized cycle length of each single intersection can be determined.Furthermore, since the turning flows keep changing in each interval, the optimized cycle length updates momentarily; that is, the optimized cycle length for each single intersection is real-time in nature.

Optimization of
The maximum cycle length is employed as the public cycle length of the arterial road for further researches.

Optimization of Real-Time Offset and Split for Arterial
Coordination Control.Based on the public cycle length, we further formulate an optimization model for arterial coordination control to optimize the real-time split and offset.
In drivers' common sense, if they do not or seldom need to stop at intersections, they will feel rather comfortable, even if the speed is not very high and there exist some delays.Therefore, different from existing models, this paper defines a new index, uninterrupted ratio, to describe the effects of the model.It is equal to the percentage of vehicles which cross the intersection without interruptions.Consider where Num   is the number of vehicles crossing the intersection  without interruption, Num  is the total number of vehicles crossing the intersection , and   is the uninterrupted ratio at intersection .
According to Figure 1, to formulate the arterial coordination control model, we must know the vehicle dispersion time first.Taking  as the key intersection, there are two situations when vehicles arrive at it along the arterial direction.Situation 1. Travel time from intersection  − 1 to  is shorter than dispersion time of vehicles accumulated at intersection  during last cycle, and vehicles from intersection  − 1 have to wait for some time to cross the intersection .Consider where    is the travel time for vehicles to travel from intersection −1 to ;   is the offset between intersection −1 and ;  is the public cycle length obtained from Section 3.2;   is the split of intersection  along the arterial direction;  is an integer number;   is the volume adjustment coefficient;   () is the arrival flow rate function at intersection , which can be assumed to follow any random distribution, and the Poisson distribution is used in this paper;   is the capacity of intersection  along the arterial direction.
Then, the dispersion time of vehicles at intersection  is formulated as follows: In the numerator of ( 6), there are a total of two items.The first item is the undispersed queue vehicles during last cycle, and the second item is the vehicles arriving at intersection  from the upstream intersection  − 1.
Situation 2. Travel time from intersection  − 1 to  is longer than dispersion time of vehicles accumulated during last cycle, and vehicles from intersection  − 1 can cross intersection  directly without waiting.Consider In this situation, According to the definition of the uninterrupted ratio, we can formulate Furthermore, to maximize the sum of uninterrupted ratio at all intersections, the real-time arterial coordination control model is formulated as where  ,max denotes the maximum split along the arterial direction (0.8 in this paper) of intersection .
Based on the public cycle length  and the vehicle dispersion time   , we can obtain the real-time offset   and the split   of each intersection along the arterial by solving the nonlinear optimization problem in (10).
From the model formulation, we can find out that there is no limitation about geographic scheme of the arterial road or vehicle arrival distribution, and the optimized offset and split remain updated along with the time interval.

Genetic Algorithm for Dynamic Turning Fraction
Estimation at Intersections.Since the objective function of ( 2) is a LAD formulation and there are equality constraints in the model, it is very difficult to be solved using traditional optimization method.Therefore, we develop a revised genetic algorithm (GA) for solution.The turning fractions   () are assumed to be 0; that is, the U-turn phenomena are neglected here.
We borrow the GA from our previous work [4] directly, and two important issues are described as following: 4.1.1.Encoding and Decoding.We use the binary encoding method in this paper.According to the equality constraint in (2), as well as the assumption that there is no U-turn at intersections, there are −2 independent turning fractions for each entering approach.Totally, there are (−2) independent turning fractions for an intersection during one time interval.These turning fractions to be estimated can be written as a matrix with  rows and  − 2 columns.
The process is further represented as where ch  is a decimal real number transformed from the binary code of   and  GA is the length of each chromosome.As stated before, there is one turning fraction which is dependent on other  − 1 results, and we assume it to be as the one corresponding to the exiting approach with the index .
Through the decoding method in (11), the estimated turning fractions will satisfy both equality and inequality constraints in (2); therefore, we do not need to perform the truncation and normalization processes [1], which are usually employed in turning proportions estimation.

Fitness Function.
Since GA always tries to search for the maximum fitness of chromosome in the evolution process, we introduce a positive number to transform the objective function to the following fitness function: where  is a positive constant and Fit is the fitness function.
The overall steps and other issues are similar to our previous work and some other existing GAs.We code this GA using M language of MATLAB software and then obtain the estimated dynamic turning fractions for each single intersection.

Solution of Real-Time Public Cycle Length Optimization at
Intersections.Since ( 3) is an ordinary nonlinear optimization problem, it can be solved directly using existing mathematical methods.We code this model using Lingo software, and then obtain the optimized cycle length of each single intersection.The maximum one is taken as the public cycle length for the arterial road.

Genetic Algorithm for Real-Time Offset and Split Optimization at Arterial Corridors.
Due to the integral formulation in the objective function of (10), its solution is very difficult.Here, we also design a GA for solving [15][16][17].The overall steps are described as follows.
Step 1 (encoding).A real-coded scheme is adopted here to represent the feasible solutions during the GA evolution process.
Step 2 (fitness function).Since the objective function of ( 10) is a maximum formulation, we take it as the fitness function directly, as shown in the following equation: Step 3 (selection).The classical Roulette wheel selection approach is employed for the selection operation in this paper.
Step 4 (crossover).Crossover operation is employed to pass the excellent genes of parent chromosomes to children chromosomes and to generate new chromosomes.In this paper, we use the following method to implement the crossover operation: where  is a random number between [0, 1], chr +1  and chr +1  are genes from children chromosomes, chr   and chr   are genes from parent chromosomes,  is the index of current generation, and  is the index of gene in the chromosome.
Step 5 (mutation).Mutation operation is implemented to ensure the local random search ability of GA, as well as to produce new children chromosomes.The following method is designed for the mutation operation: where  and  are random numbers between [0, 1], chr ,min is the lower bound of chr  , chr ,max is the upper bound of chr  , and Δ max is the maximum generations.
Some important controlling parameters of the GA are set as follows: the population size is 60, the maximum number of generations is 120, the crossover rate is 0.7, and the mutation rate is 0.01.
We also code this GA using M language of MATLAB software and then obtain the optimized offset and split for the arterial road.

Study
To investigate the performance of the proposed model framework, this paper implements a case study using practical data and the effects of the arterial control are evaluated through traffic simulation.
An arterial corridor around Tianshuiyuan block in Chaoyang district, in Beijing city, is taken as the case study area.Figure 2 presents a graphical illustration of the area.
Tianshuiyuan Street is the arterial road to be studied, and there are three intersections along it.Intersection (I-2) is employed as the key intersection, which is a crossroad with four approaches.The arterial direction is from south to north, that is, from I-1 to I-3.
To collect the field data, we implemented a survey around the case area and collected a great deal of data, including the time-varying entering and exiting link flows at all three intersections, the geographical information of the arterial, the existing signal timing plan, and queue length, delay time, and number of stops at all three intersections.All information needed in the case study can be extracted from the above survey data.
We further implemented the proposed model, as well as the existing MAXBAND method.Table 1 reports cycle length, arterial split, and offset of each intersection from current scheme, MAXBAND model, and the proposed model.
In the current scheme, all three intersections have fixed timing plans, and the arterial road is not coordinated at all.For the proposed model, we take 5 cycles, as examples here, and the public cycle length, arterial split, and offset are all illustrated in Table 1.We can also find out that all signal timings of the proposed model are updated in each cycle; therefore, they are real-time in nature.
To further evaluate the effects of the proposed model, we select three indices as evaluation criteria, including queue length, delay time, and number of stops.All these indices of current scheme were obtained through the field survey.Since it is impossible to adjust the signal timings in real world just for this case study, we simulated both MAXBAND method and the proposed model using Vissim software [21].
Figure 3 shows the simulation environment.The road network was created based on an original AutoCAD file, and the surveyed entering link counts at three intersections were taken as the input data, together with the current turning proportions and signal timings.Using the above three indices as evaluation criteria, the simulation model was calibrated by adjusting some parameters, such as speed, acceleration, and priority rules.The timings from the proposed methodology were then input to the simulation model, and all three evaluation indices under proposed timings were extracted through some detectors laid in the road network in Vissim.
We further compared the three simulated evaluation indices with those from current traffic control scheme.Tables 2, 3, and 4 report queue length, delay time, and number of stops, respectively.
The detailed graphical comparisons of the three evaluation indices are further illustrated in Figures 4, 5, and 6.
From Tables 2, 3, and 4 and Figures 4, 5, and 6, we can reach the following results.(1) Generally, the proposed arterial coordination control model yields the best results during all cycles in terms of all three indices, including queue length, delay time, and number of stops.
(2) The MAXBAND method also outperforms current scheme during most cycles in all three indices, because the current timing plan is not coordinated along the arterial.
(3) For delay time, the proposed model gets a little worse result than MAXBAND in the first cycle, which is due to some unknown reasons, but it is still better than the current scheme and does not influence the general outstanding performances of the proposed model.

Conclusions
This paper presents a general model framework for realtime arterial coordination control.To provide important input data for the arterial control, this study first proposes a dynamic intersection turning fraction estimation model using the least absolute deviation formulation and designs a genetic algorithm for solution, integrating revised encoding and decoding methods.Based on the estimated timevarying turning proportions, this paper develops a nonlinear optimization model based on minimum delay to optimize the real-time cycle of each single intersection along the arterial road, and the maximum cycle is taken as the public arterial cycle.Furthermore, this paper puts forward a novel optimization model based on minimum uninterrupted ratio and develops a genetic algorithm to optimize both offset and split of the arterial road, which are also time-varying.Time-dependent public cycle, offset, and split constitute the real-time timing plan of arterial coordination control scheme together.The proposed model removes most of the restrictions of conventional arterial control models and is a real-time control method in nature.The simulation experiments based on field data have confirmed the outstanding performances of the proposed model framework compared with both current scheme and MAXBAND method.
This paper can be enhanced in following directions.The first is to further consider the opposite direction of the arterial corridor and construct a bidirectional arterial coordination control model.The second is to extend the model to a wide area and develop a regional coordination control model.And the third is to integrate all three submodels in the framework together and formulate an accurate and efficient combined model.

Figure 1 :
Figure 1: Architecture of the arterial coordination control.

Figure 2 :
Figure 2: Graphical illustration of the case study arterial road.

Figure 3 :Figure 4 :
Figure 3: Graphical illustration of the simulation environment.

Figure 5 :
Figure 5: Graphical illustration of delay time.

Figure 6 :
Figure 6: Graphical illustration of number of stops.

Table 2 :
Comparison of queue length.

Table 3 :
Comparison of delay time.

Table
Comparison of number of stops.