Discrete Dynamics in Nature and Society

Volume 2012 (2012), Article ID 858731, 14 pages

doi:10.1155/2012/858731

Research Article

Convergence Study of Minimizing the Nonconvex Total Delay Using the Lane-Based Optimization Method for Signal-Controlled Junctions

C. K. Wong and Y. Y. Lee

Department of Civil and Architectural Engineering, City University of Hong Kong, Tat Chee Road, Kowloon, Hong Kong

Received 12 January 2012; Revised 8 March 2012; Accepted 20 March 2012

Academic Editor: Beatrice Paternoster

Copyright © 2012 C. K. Wong and Y. Y. Lee. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This paper presents a 2D convergence density criterion for minimizing the total junction delay at isolated junctions in the lane-based optimization framework. The lane-based method integrates the design of lane markings and signal settings for traffic movements in a unified framework. The problem of delay minimization is formulated as a Binary Mix Integer Non Linear Program (BMINLP). A cutting plane algorithm can be applied to solve this difficult BMINLP problem by adding hyperplanes sequentially until sufficient numbers of planes are created in the form of solution constraints to replicate the original nonlinear surface in the solution space. A set of constraints is set up to ensure the feasibility and safety of the resultant optimized lane markings and signal settings. The main difficulty to solve this high-dimension nonlinear nonconvex delay minimization problem using cutting plane algorithm is the requirement of substantial computational efforts to reach a good-quality solution while approximating the nonlinear solution space. A new stopping criterion is proposed by monitoring a 2D convergence density to obtain a converged solution. A numerical example is given to demonstrate the effectiveness of the proposed methodology. The cutting-plane algorithm producing an effective signal design will become more computationally attractive with adopting the proposed stopping criterion.

1. Introduction

Over the past decades, various traffic modeling techniques have been developed and studied by many researchers (e.g., [1–8]). In the context of signal settings optimization at isolated junctions, it is always assumed that traffic pattern is arrived dynamically according to Poisson distribution. Pioneering works can be found in [9–11] for the stage-based signal timing design approach and incorporated into a commercial package called OSCADY. Improta and Cantarella [12] proposed a better phase-based approach to the problem of calculating signal settings. The cycle structure was specified by a set of binary variables relating to incompatible signal groups. A branch-and-bound technique was employed to solve the Binary Mixed Integer Linear Program (BMILP). Gallivan and Heydecker [13] developed a similar approach, in which the cycle structure was specified by one of the stage sequences. Heydecker [14] further reduced this problem difficulty by introducing a procedure to group all possibilities into a much smaller number of equivalence classes. The cycle structure of each equivalence class was represented by a “successor function” that makes the method computationally faster.

Wong and Heydecker [15] and Wong and Wong [16] proposed a lane-based optimization method for determining the capacity-maximizing and cycle-length-minimizing signal settings for an isolated junction. The present study extends the lane-based optimization method for determining a combined set of lane markings and signal settings to minimize the traffic delay at an isolated signal-controlled junction. A cutting plane algorithm with a faster convergence criterion is proposed to solve the delay minimization problem. A numerical example is given to demonstrate the effectiveness of the proposed methodology. On the other hand, the concept of the convergence density in this paper originates from the posterior probability densities of the parameter values of a model in the probabilistic system identification method (e.g., [17–20]). Instead of pinpointing particular values assigned to the parameters of a model, the probabilistic identification method focuses on those parameter values with the highest posterior probability densities. This is similar to that the stopping criterion in this paper focuses on the subdomain formed by a set of gridlines with the highest density and the corresponding parameter values.

2. The Lane-Based Optimization Method

Saturation flow of a traffic lane prescribes the maximum number of vehicles that can leave a junction if full-green signal is given, and this discharge rate is defined by the following expression, which is generalized from that of Kimber et al. [21], that gives the maximum discharge flow rate of traffic flow on a traffic lane of different lane type, width, turning traffic proportion, and physical turning radius:

where is the saturation flow of lane in arm i, and are the radius of the turning trajectory (= ∞ for straight-ahead movement) and the proportion of flow from arm to arm via lane k, and is the saturation flow of the lane if it is for straight-ahead movement only. In the lane-based optimization method, the control variables can be specified as follows. Let be the set of control variables, where and are the subsets of binary and continuous variables, respectively.

The subset consists of the following binary variables.

Permitted movements: .Successor functions: .

The subset consists of the following continuous variables.

Designed lane flows: .Common flow multiplier: μ.Cycle length: ζ.Starts of green for movements: and .Durations of green for movements: and .Starts of green for traffic lanes: .Durations of green for traffic lanes: ,

where are subscripts denoting arms connecting to a junction, is subscript denoting a traffic lane, is total number of traffic arms in a junction, is number of traffic lane from arm , is assigned lane flow on traffic lane for the traffic turning from arm to arm , is given demand of traffic flow turning from arm to arm , is the common flow multiplier, is the permitted movement on a traffic lane for the traffic turning from arm to arm in the form of a lane marking (= 1 exists/permitted or = 0 otherwise), is a user-defined mathematical function to find out the exit direction (i.e., which arm the approach traffic is leaving) from the jth destination arm counting from arm i, is maximum number of exit lane for leaving traffic toward arm , is the operating cycle length (reciprocal of cycle time), is the maximum allowable limit of cycle time, is the minimum allowable limit of cycle time, is an arbitrary large integral number, is the start of green for a signal phase controlling the traffic turn from arm to arm , is the start of green for a signal phase controlling the traffic turn from arm on lane , is the duration of green for a signal phase controlling the traffic turn from arm to arm , is the duration of green for a signal phase controlling the traffic turn from arm on lane , is the minimum duration of green for a signal phase for the traffic turn from arm to arm , is a successor function governing the sequence of signal display within a signal cycle (= 1 if start of green of signal group follows that of signal group or = 0 if the opposite is true), is set of incompatible signal groups, Ψ is set of incompatible traffic movements, is saturation flow (maximum discharge flow rate) for straight-ahead movement on a traffic lane from arm , is turning radius for traffic turning from arm to arm one lane k, is difference between effective green time and display green time durations (usually 1 second), and is maximum allowable degree of saturation (demand flow to saturation flow ratio) on a traffic lane from arm (usually take 90%).

To formulate the present total delay optimization problem in the form of a Binary Mixed Integer Linear Program, linear constraint sets below are required to outline the feasible solution region ensuring safe operation of traffic signal plans and logical sequence of signal timings.

Flow Conservation Balancing the Traffic Demand and Optimal Lane Flow Pattern

Consider

Equation (2.2) is an equality constraint set mainly to ensure the sum of all optimized lane flows to be assigned on each traffic lane that will eventually be equal to the input turning demand flows.

Minimum Permitted Movement on a Lane

Consider

This set of constraint is developed to ensure each traffic lane will be utilized to let through at least one traffic movement.

Maximum Permitted Movements at the Exit

Consider

Equation (2.4) is a constraint set to make sure that there is sufficient number of traffic lanes at downstream exit to manage the vehicles leaving from upstream approach lanes for various turning directions.

Permitted Movements Across Adjacent Lanes

Consider

To prevent internal conflicts, for example, a right turn on left-hand lane and a left turn on a right-hand lane in which vehicles will be clashed during green signal phases, constraint set in (2.5) is necessary to restrict the formation of lane marking on adjacent traffic lane from the same approach.

Cycle Length

Consider

Practically, all signal timings will be repeated in cycles in which a signal cycle should be selected within maximum and minimum bounds.Identical lane signal settings for common turns on lanes:

In the present formulation, green times will be defined on both lane basis or turning basis. Once a specific traffic turn is assigned on a particular traffic lane, their signal settings should be identical and constraint sets in (2.7) should be included.

Start of Green

Consider

Having said that all signal timings will be repeated in signal cycles and thus the start of green time should always be chosen within a signal cycle, (2.8) is given to maintain this design criteria.

Duration of Green

Consider

Similarly, the assigned green time duration for a turning movement should always be less than the duration of a signal cycle. For safety concern, traffic signals should not be changed too frequently so that the assigned green time should always be greater than a user-defined minimum green time. Equation (2.9) is developed for this purpose.

Order of Signal Displays

Consider

Since conflicting traffic movements from different approaches cannot receive green times simultaneously, a successor function and constraint set in (2.10) should be defined to order the signal display in a safe sequence.

Clearance Time

Consider

With the introduction of a successor function, the order of signal displays can be designed. For those conflicting movements, a clearance time should be given so that the traffic in the common area of a junction can be cleared before other conflicting traffic movement can enter the junction for discharging.

Prohibited Movement

Consider

To design optimal lane marking pattern, some vehicle turnings will be banned on different traffic lanes. If the lane marking for a particular vehicle turn is prohibited, the corresponding assigned lane flow must be forced to be zero. Hence, (2.12) is formulated as one of the design constraint sets.

Flow Factor

Consider

If an identical lane marking is assigned on two adjacent lanes, the flow factors for the two traffic lanes must be identical to ensure the queue lengths are equal. Equation (2.13) should be included to prevent unequal vehicle queue development.

Maximum Acceptable Degree of Saturation

Consider

In the present design framework, it is necessary to provide sufficient green time for different traffic movements such that their demand flow to capacity ratios are always within a user acceptable level. Equation (2.14) is developed to ensure the assignment of enough green time to handle the demand traffic.For isolated signal control junctions, there are three common criteria for optimizing the signal settings: capacity maximization, cycle length minimization, and delay minimization. Alternative signal settings can be generated to fulfill specific operating requirements if different optimizing objectives are chosen. In general, the delay function that is used as the objective for optimization is nonlinear. The problem has to be formulated as the BMINLP that is shown as follows:

subject to the linear constraints in (2.2)–(2.14) and , where defines the total delay of the junction.

3. Delay Optimization and 2D Convergence Density

Due to the complicated high dimension solution space and the nonconvexity of the objective delay function, the piecewise linearization approach is not applicable to the present problem. A classical cutting plane algorithm with a 2D convergence density is proposed and discussed in the following sections to solve the BMINLP problem. Now, consider the auxiliary continuous variable as the objective function which is the approximated delay found in the cutting plane algorithm (it is approximated as the solution is from the set of linearized hyperplane but not from the original nonlinear delay function), with the following additional nonlinear constraint:

Then, the stopping criterion is considered and the 2D convergence density function is defined by

where and ; is the number of iterations; Count is a function which counts the times of and that match the specified criteria. and are the coordinates of any grid point within the domains of and . and are the numerical values at the solution point. Note that the domains and are divided by a set of rectangular grid lines; the interval between any two grid lines is . The stopping requirement is that there is one grid box of and with the highest convergence density value for the last n iterations. If the interval of is smaller (i.e., the grid box size is smaller), it requires more iterations for the coordinates of and staying the same grid box in the last iterations to achieve the highest convergence density value. It is noted that the nonlinear constraint is linearized by a first-order Taylor’s series expansion to form a linear function , where

in which specifies the th linearization point, , and is a parameter that controls the degree of a cut in the feasible region. The lane-based delay minimization problem becomes the minimization of , subject to constraints in (2.2)–(3.1), , and all of the constraints (3.3) that were generated in previous iterations.

4. Numerical Example

The following steady-state delay formula, based on random arrivals and regular departure patterns, is used [9, 10]:

where is the rate of delay for the traffic lane on arm of a signal junction. The total rate of delay of the junction,, is the sum of the delays for all traffic lanes on all arms . This formula is widely used for the estimation of junction delay.

Consider a four-arm junction with four traffic lanes on each arm, as shown in Figure 1. The maximum cycle length is set at 120 seconds, and the maximum acceptable degree of saturation is 90% in all lanes. The minimum greens are 5 seconds for traffic movements. The traffic demands are given in Table 1. The required clearance time for any two conflicting movements (including both traffic and pedestrian movements) is 6 seconds. We also take the solution of capacity maximization as the starting point for linearization in the cutting plane algorithm. Set all initial ’s = 1.0, and a multiplication factor 10.0 is applied to update the ’s each time there is an invalid underestimator. The cutting plane algorithm is then solved by a commercial package called “The mathematical programming language” in which a CPLEX solver is installed to solving the Binary Mixed Integer Linear Program.

Table 1: The traffic demand for the example junction.

Traffic demand  in veh/hr
To arm

1
2
3
4

From arm
1
—
500
200
100

2
100
—
100
500

3
300
300
—
300

4
100
400
400
—

Figure 1: The layout of the four-arm example junction.

Now, consider the stopping requirements in (3.2). In this example (see Figure 2), and = 2 (the domain of and are 0 to 50). The two-dimensional plot of the convergence density is shown in Figure 3. It can be seen that from 84th to 93rd iterations, the convergence density value of the marked grid box is the highest. The minimum delay from the cutting plane algorithm is 26.683 veh-s/s, and the corresponding optimized cycle length to operate the junction is 71.24 seconds.

(a)

(b)

(c)

(d)

Figure 2: Two-dimensional plot of the convergence density.

Figure 3: Optimized lane markings for the example problem.

Typical results of the proposed solution algorithms, including the lane permitted movements, designed lane flows, and the signal settings, are given in Table 2. The total delays, based on Webster’s delay function, are calculated and put in the last column of the tables. It can also be verified that the degrees of saturation for all traffic lanes are within the specified upper limit of 90%. Corresponding optimized lane marking pattern is given in Figure 3 diagrammatically. A summary of the model statistics for the example calculations is presented in Table 3.

Table 2: Optimization results from the cutting plane algorithm.

Optimal cycle length = 71.24 seconds

(2.1)
(2.2)
(2.3)–(2.6)
(2.7)
(8)
(2.8)
(2.9)
(2.10)
(2.11)
(2.12)
(2.13)

From arm
Lane
To arm
Sum of lane flow
Turning proportion
Saturation flow
y
Effective green
Degree of sat.
Green start
Total delay

1
2
3
4

1
1

241.4

241.4
1.00
1746.7
0.1382
27.8
0.3541
0.0
1.0148

1
2

258.6

258.6
1.00
1871.1
0.1382
27.8
0.3541
0.0
1.0808

1
3

200.0

200.0
0.00
2105.0
0.0950
11.0
0.6129
16.8
1.8419

1
4

100.0
100.0
1.00
1871.1
0.0534
6.0
0.6346
35.4
1.2849

2
1

48.3

48.3
1.00
1746.7
0.0276
9.4
0.2086
56.8
0.3575

2
2

51.7

51.7
1.00
1871.1
0.0276
9.4
0.2086
56.8
0.3812

2
3

500.0
500.0
0.00
2105.0
0.2375
19.9
0.8514
46.4
5.2316

2
4
100.0

100.0
1.00
1871.1
0.0534
11.8
0.3238
0.0
0.7256

3
1

300.0
300.0
1.00
1746.7
0.1718
30.4
0.4029
0.0
1.1842

3
2
300.0

300.0
0.00
2105.0
0.1425
13.6
0.7462
16.8
3.0261

3
3

150.0

150.0
1.00
1871.1
0.0802
8.6
0.6669
32.8
1.7250

3
4

150.0

150.0
1.00
1871.1
0.0802
8.6
0.6669
32.8
1.7250

4
1
100.0

100.0
1.00
1746.7
0.0573
30.9
0.1321
35.4
0.3123

4
2

400.0

400.0
0.00
2105.0
0.1900
19.9
0.6811
46.4
2.9411

4
3

200.0

200.0
1.00
1871.1
0.1069
11.8
0.6475
71.2
1.9255

4
4

200.0

200.0
1.00
1871.1
0.1069
11.8
0.6475
71.2
1.9255

26.6830

Table 3: Statistics of the example problem.

Model variables
Number of variables
Constraints
Number of constraints

Binary variables
142
Cutting planes
93

    Permitted movements
40
Starts of green
58

    Successor functions
102
Durations of green
60

Order of signal displays
51

Continuous variables
110
Clearance times
68

    Lane flows
48
Cycle length
1

    Starts of green
29
Maximum acceptable degree of saturation
16

    Durations of green
30
Flow conservations
12

    Cycle length
1
Minimum permitted movements in a lane
8

    Common flow multiplier
1
Maximum permitted movements at exits
3

    Auxiliary variable
1
Prohibited movements
40

Permitted movements across adjacent lanes
72

Flow factors
72

Lane signal settings
176

Others
53

Total number
252

783

Figures 4(a) and 4(b) show the number of all iterations and the minimum delay against the number of iterations in the highest density box, respectively. It can be seen that the smaller the grid size is, the higher the number of all iterations is; the higher the number of iterations required in the highest density box is, the higher the number of all iterations is. Figures 4(c) and 4(d) show the number of all iterations and the minimum delay against the relative error which is the main parameter in the stopping criteria in [22]. Similar to the observation in Figures 4(a) and 4(b), the minimum delay can converge to a certain value when the number of all iterations increases. The main differences between the convergence patterns in the present study and [22] are that the curves representing the numbers of all iterations in Figure 4(a) are more linear; the minimum delay in Figure 4(b) can converge faster to a certain value against the number of iterations required in the highest density box (the curve looks more concave).

(a)

(b)

(c)

(d)

Figure 4: (a) The number of all iterations against the number of iterations in the highest density box. (b) The minimum delay against the number of iterations in the highest density box. (c) The number of all iterations plotted against the relative error (the stop criteria in [22]). (d) The minimum delay against the relative error.

5. Conclusions

Lane-based optimization models have been presented in this paper. While the capacity maximization and cycle length minimization for isolated signal-control junctions are formulated as BMILP problems that can be solved by the standard branch-and-bound technique [16], the delay minimization problem for isolated junctions is a BMINLP because of the nonlinear delay function [22]. The cutting plane algorithms with a 2D convergence density stopping criterion have been proposed to solve the problem. The nonlinear delay function is approximated by successive linearizations. The BMINLP problem is transformed into a series of BMILP problems that are solved by the standard branch-and-bound technique. To achieve a more cost-effective method for the lane-based optimization of delay at isolated signalized junctions, the cutting plane algorithm should be applied together with the new proposed stopping criterion taking the 2D convergence density into considerations.

Acknowledgment

This research was supported by a research grant (SRG 7008031) from the City University of Hong Kong.

References

A. S. Alfa, “Modelling traffic queues at a signalized intersection with vehicle-actuated control and Markovian arrival processes,” Computers & Mathematics with Applications, vol. 30, no. 8, pp. 105–119, 1995.

S. Blandin, G. Bretti, A. Cutolo, and B. Piccoli, “Numerical simulations of traffic data via fluid dynamic approach,” Applied Mathematics and Computation, vol. 210, no. 2, pp. 441–454, 2009.

S.-W. Chiou, “Reserve capacity of signal-controlled road network,” Applied Mathematics and Computation, vol. 190, no. 2, pp. 1602–1611, 2007.

S.-W. Chiou, “An optimization model for area traffic control with link tolls,” Applied Mathematics and Computation, vol. 192, no. 2, pp. 520–532, 2007.

A. de Palma, C. Lefèvre, and M. Ben-Akiva, “A dynamic model of peak period traffic flows and delays in a corridor,” Computers & Mathematics with Applications, vol. 14, no. 3, pp. 201–223, 1987.

W.-C. Hong, Y. Dong, F. Zheng, and S. Y. Wei, “Hybrid evolutionary algorithms in a SVR traffic flow forecasting model,” Applied Mathematics and Computation, vol. 217, no. 15, pp. 6733–6747, 2011.

D. C. Park, “Prediction of MPEG video traffic over ATM networks using dynamic bilinear recurrent neural network,” Applied Mathematics and Computation, vol. 205, no. 2, pp. 648–657, 2007.

A. H. Shahpar, H. Z. Aashtiani, and A. Babazadeh, “Dynamic penalty function method for the side constrained traffic assignment problem,” Applied Mathematics and Computation, vol. 206, no. 1, pp. 332–345, 2008.

F.V. Webster, “Traffic signal settings,” Road Research Technical Paper number 39, HMSO, London, UK, 1958.

R. E. Allsop, “Delay-minimizing settings for fixed-time traffic signals at a single road junction,” Journal of the Institute of Mathematics and its Applications, vol. 8, pp. 164–185, 1971.

R. E. Allsop, “Estimating the traffic capacity of a signalized road junction,” Transportation Research, vol. 6, pp. 245–255, 1972.

G. Improta and G. E. Cantarella, “Control system design for an individual signalized junction,” Transportation Research, vol. 18, no. 2, pp. 147–167, 1984.

S. Gallivan and B. G. Heydecker, “Optimising the control performance of traffic signals at a single junction,” Transportation Research, vol. 22, pp. 357–370, 1988.

B. G. Heydecker, “Sequencing of traffic signals,” in Mathematics in Transport and Planning and Control, J. D. Griffiths, Ed., pp. 57–67, Clarendon Press, Oxford, UK, 1992.

C. K. Wong and B. G. Heydecker, “Optimal allocation of turns to lanes at an isolated signal-controlled junction,” Transportation Research Part B, vol. 45, no. 4, pp. 667–681, 2011.

C. K. Wong and S. C. Wong, “Lane-based optimization of signal timings for isolated junctions,” Transportation Research Part B, vol. 37, no. 1, pp. 63–84, 2003.

D. Capecchi and F. Vestroni, “Identification of finite element models in structural dynamics,” Engineering Structures, vol. 15, no. 1, pp. 21–30, 1993.

W. T. Kung, Y. Y. Lee, and H. Y. Sun, “Sound leakage identification for an enclosed room using the probabilistic approach and model class selection index: an experiment,” Journal of Sound and Vibration, vol. 310, pp. 776–781, 2008.

V. Peterka, “Bayesian system identification,” Automatica, vol. 17, no. 1, pp. 41–53, 1981.

H. Y. Sun, Y. Y. Lee, A. Y. T. Leung, and C. F. Ng, “Distinguishing between the interior pressures induced from two independent sources within a room using the probabilistic approach,” Applied Acoustics, vol. 69, no. 6, pp. 496–505, 2008.

R. M. Kimber, M. McDonald, and N. Hounsell, “The prediction of saturation flows for road junctions controlled by traffic signals,” TRRL Report RR 67, Transport and Road Research Laboratory, Crowthorne, 1986.

C. K. Wong and S. C. Wong, “A lane-based optimization method for minimizing delay at isolated signal-controlled junctions,” Journal of Mathematical Modelling and Algorithms, vol. 2, no. 4, pp. 379–406, 2003.