Research Article  Open Access
Ozgur Baskan, Huseyin Ceylan, Cenk Ozan, "A Simultaneous Solution for Reserve Capacity Maximization and Delay Minimization Problems in Signalized Road Networks", Journal of Advanced Transportation, vol. 2019, Article ID 6203137, 18 pages, 2019. https://doi.org/10.1155/2019/6203137
A Simultaneous Solution for Reserve Capacity Maximization and Delay Minimization Problems in Signalized Road Networks
Abstract
In this study, we present a bilevel programming model in which upper level is defined as a biobjective problem and the lower level is considered as a stochastic user equilibrium assignment problem. It is clear that the biobjective problem has two objectives: the first maximizes the reserve capacity whereas the second minimizes performance index of a road network. We use a weightedsum method to determine the Pareto optimal solutions of the biobjective problem by applying normalization approach for making the objective functions dimensionless. Following, a differential evolution based heuristic solution algorithm is introduced to overcome the problem presented by use of biobjective bilevel programming model. The first numerical test is conducted on twojunction network in order to represent the effect of the weighting on the solution of combined reserve capacity maximization and delay minimization problem. Allsop & Charlesworth’s network, which is a widely preferred road network in the literature, is selected for the second numerical application in order to present the applicability of the proposed model on a mediumsized signalized road network. Results support authorities who should usually make a choice between two conflicting issues, namely, reserve capacity maximization and delay minimization.
1. Introduction
As it is well known, road users may be delayed at signalized intersections on urban roads because of implementing inappropriate signal timings even if traffic flow is less than capacity. In order to minimize delay, signal timing optimization concept has widely been used on intersection, arterial or network levels. Implementation of proper signal timing plans on a road network most likely leads to a decrease on flow/capacity ratios and reveals the spare capacity of road segments. This spare capacity can be utilized until the increasing travel demand reaches a certain level. The problem of how to determine this certain level of travel demand, which leads to an flow/capacity ratio equal to one for at least one road segment on the network, has been known as reserve capacity maximization problem (RCMP). The reserve capacity can also be defined as the largest applicable demand multiplier to an existing OriginDestination (OD) demand matrix without violating the capacity of any link subject to signal timing parameters and other related constraints. By means of maximizing the reserve capacity of a road network, decision makers reveal the spare capacity and the network becomes more reliable when taking sudden spikes in travel demand into account. Local authorities or agencies that are responsible for traffic operations can benefit from this improvement. On the other hand, road users may suffer from growing travel demand and correspondingly increasing delay. As it is well known, delay function is mainly used as an objective to optimize signal timings for signalized road networks. Although optimal signal timings can be easily found for isolated intersections, signal timing optimization in coordinated road networks becomes highly challenging issue due to the vast search space considering the offset and common cycle time components. Relevant literature provides a large number of studies concerning both reserve capacity maximization and delay minimization on road networks. From the point of RCMP, Webster and Cobbe [1] firstly studied on this problem and they determined the reserve capacity of a signalized intersection by developing an explicit formulation. In order to generalize this approach, Allsop [2] proposed a new method for junctions with complicated signal plans. Afterward, Yagar [3, 4] improved the linear program proposed by Allsop [2] considering different saturation flows between successive stages. Wong and Yang [5] applied the reserve capacity approach to a signalized network by considering user equilibrium assignment. They aimed to find the maximum OD demand multiplier subject to link flows not exceeding their capacities. Yang et al. [6] proposed a bilevel programming model (BLPM) by using combined assignment and trip distribution model to find the network reserve capacity considering equilibrium constraints. Yang and Wang [7] investigated the results of applying two different objectives, namely, reserve capacity maximization and travel cost minimization. It was found that the relationship between these objectives depended on the level of congestion. At the same year, Ziyou and Yifan [8] dealt with the RCMP in a different viewpoint; that is, they considered different demand multipliers between OD pairs to create a condition closer to the reality. Ge et al. [9] investigated whether the reserve capacity was affected by the road users informed about the network or not. Results showed that the reserve capacity did not increase in a similar way with the information level of road users. Ceylan and Bell [10] solved the RCMP by optimizing traffic signal timings including offset term in different stages. They developed a twostage algorithm in the stochastic user equilibrium (SUE) manner. Chen et al. [11] proposed an approach considering both reliability and capacity to solve the RCMP. For this purpose, they developed a new index to determine the demand multiplier by using a BLPM. Chiou [12] solved the RCMP based on a projected gradient method. The proposed method was compared with traditional methods and results showed that it outperformed other methods by means of reducing delay.
As good examples in terms of biobjective problems, Chiou [13, 14] presented novel algorithms for the solution of the RCMP with toll settings. Encouraging results were attained from two studies through the provided numerical examples. Similarly, Chiou [15] introduced a hybrid approach to maximize the reserve capacity and to minimize delay based on link capacity expansions. Another type of biobjective problems has been tackled with minimizing delay and maximizing reserve capacity with timeinvariant flows by Chiou [16]. From a different viewpoint, Miandoabchi and Farahani [17] solved the RCMP in the context of discrete network design problem using a bilevel solution approach. In a similar manner, Miandoabchi et al. [18] handled this combined problem related to road network design to optimize network reserve capacity and some related performance measures. Chiou [19] defined the reserve capacity maximization and delay minimization problems as a minmax problem subject to the equilibrium constraints. Proposed model was applied to two signal controlled benchmark networks, which have widely been used in studies concerning network design problem, and results showed that the model was capable of solving the RCMP by taking delay minimization problem into account. Wang et al. [20] solved the RCMP with SUE link flows and subsequently extended it by considering maximum link capacity expansions with limited budget. Xiao et al. [21] investigated the RCMP from a different viewpoint, and they developed a zonebased model to reveal the effects of using maximum capacity on the landuse development. Recently, Han and Cheng [22] studied on the effectiveness of the tradable credit scheme, which is a particular type of congestion pricing applications, by considering the RCMP with SUE link flows. To this end, they developed a bilevel model which was applied to two road networks to illustrate the positive effect of the tradable credit scheme. Baskan and Ozan [23] aimed to maximize the reserve capacity by considering equity issue with deterministic user equilibrium link flows. They presented a bilevel heuristic solution algorithm based on Harmony Search (HS) technique. Results indicated that the value of OD demand multiplier was highly sensitive to the equity parameter. Recently, Li and Sun [24] proposed a novel optimization method for signal control in urban road networks by drawing attention to the lack of studies on multiobjective optimization. Road networks with different scales were used in order to reveal the applicability of the proposed method by considering network throughput, delay, and spillbacks. Results showed that the proposed multiobjective optimization method was quite promising.
From the point of delay minimization problem, TRANSYT7F traffic software, which consists of a traffic model and a signal timing optimizer, is nearly most preferred tool for signalized networks (see for details McTrans Center [25]) in order to optimize signal timings. It simulates traffic in a signalized road network to produce flow profiles which are then used to determine a performance index (PI) by considering specific stage plans and signal timing parameters. The PI, which is used to define the cost of traffic congestion in a road network, is defined as the sum of a weighted combination of delay and number of stops. In this context, several studies have been presented in the literature. As known, Webster [26] presented one of the first studies using mathematical framework in order to minimize delay on an isolated junction. This study was further extended by Allsop [27] to develop a more general approach in which total delay was minimized by formulating signal setting problem as a convex mathematical programming. After these first attempts to optimize traffic signal timings, Wong [28] carried perspective on signal timing optimization a step further and developed approximate mathematical expressions for the derivatives of PI by considering phasebased control variables. It was concluded from this study that the use of approximate expressions provided benefit in terms of computational effort. A year later, Heydecker [29] applied decomposition approach to optimize traffic signal timings for road networks. Application of this method to a small example network showed significant benefits in terms of network performance. Wong [30] used groupbased control variables to optimize signal timings with fixed link flows. For this purpose, the PI was considered as a function of groupbased control variables. To overcome some drawbacks in the study by Wong [30] in terms of computational effort, parallel computing technique was presented to reduce computational time by Wong [31]. In Wong et al. [32], TRANSYT traffic model was converted to a timedependent model to determine PI. Girianna and Benekohal [33] developed two different heuristic algorithms based on Genetic Algorithm (GA) for oversaturated signalized networks. Results indicated that microGA method was able to reach quickly to the nearglobal signal timings in comparison with the traditional GA. Ceylan and Bell [34] introduced a GAbased model to optimize signal timings with SUE link flows. They employed TRANSYT software to calculate performance index of a road network whereas SUE link flows were calculated using path flow estimator. Similarly, Ceylan and Bell [35], a year later, used the same GAbased model to reveal its capability in finding optimal signal timings under congestion with different demand conditions. Ceylan [36] proposed a GAbased model using TRANSYT HillClimbing (HC) optimization tool for area traffic control. The proposed algorithm outperformed TRANSYT in terms of performance index although it had some deficiencies regarding computational burden due to the use of HC optimization tool. Chen and Xu [37] used Particle Swarm Optimization (PSO) algorithm in order to find optimal signal timings. Their results revealed that the proposed algorithm was applicable especially in different demand conditions. Dan and Xiaohong [38] investigated the use of improved version of the GA to solve signal timing problem by taking the offset term into account and they reached valuable results in terms of minimizing delay and improving the network capacity. Li [39] developed an optimization model considering an arterial road to find optimal signal settings under oversaturated conditions. Moreover, developed model was able to take complex flow interaction between lanes into account with the cell transmission concept. Liu and Chang [40] presented a new signal optimization model considering mutual interaction between queues occurred on different lanes on arterial links. Results showed that their model outperformed TRANSYT7F for planning arterial signals. Ceylan and Ceylan [41] developed a hybrid algorithm integrating HS and HC optimization methods. In order to find optimal signal settings, the performance index was determined by TRANSYT traffic software. Results revealed that the proposed algorithm outperformed conventional HS and GA methods in terms of PI. Dell’Orco et al. [42] presented a HSbased method using BLPM for signalized road networks by considering SUE link flows. Similarly, Dell’Orco et al. [43] attempted to solve area traffic control problem by using Artificial Bee Colony (ABC) algorithm with TRANSYT7F. Results revealed that the ABC algorithm outperformed the GA and HC methods in terms of PI. Zhu et al. [44] developed Reinforcement Learning (RL) based Junction Tree Algorithm for solving signal control problem in which agents present traffic signals. Results proved that the proposed algorithm was able to produce better results than other RLbased methods. Ozan et al. [45] developed a modified RL algorithm to optimize traffic signal timings for road networks. Results showed that the proposed algorithm was capable of finding optimal signal parameters even in the case of increasing demand conditions. Christofa et al. [46] presented a realtime signal control system aiming to minimize person delay by optimizing signal timings based on mixed integer linear programming approach. Field tests proved that the person based minimization approach outperforms fixedtime signal timings. Zhang et al. [47] developed two novel models in order to solve signal coordination problem especially for long arterials and grid networks. Proposed models were in the form of a mixedinteger linear program. Numerical experiments revealed that both models were able to produce consistent coordination plans in terms of travel speed, stop, and delay. Srivastava and Sahana [48] proposed a hybrid evolutionary approach to optimize signal timings using BLPM. Comparative experiments showed that the hybrid model outperforms Ant Colony Optimization (ACO) and GA techniques. Recently, Memoli et al. [49] focused on signal setting design at networklevel by introducing a stagebased method. In proposed method, stage sequences, stage lengths, and offsets were optimized simultaneously. In order to show the applicability of the proposed model, two example test networks have been analyzed. Results indicated that the developed method was an effective way to design network signal settings.
In the relevant literature, reserve capacity maximization or delay minimization has been mostly performed by optimizing traffic signal timings. However, to the best of our knowledge, there is no study in which they are simultaneously considered. Although Ceylan and Bell [10] presented a twostage approach to the solution of the RCMP by taking delay minimization into account, their study had a limitation, that is, the signal timings were optimized in one stage by minimizing PI. The largest OD matrix multiplier was found in the second stage subject to the best signal timings found in the first stage. This means that the biobjective BLPM was separately solved in their study. Therefore, we define a new approach for investigating the advantage of the simultaneous solution of the RCMP and delay minimization problem which are both considered as a biobjective problem. To evaluate the users’ reaction to this arrangement performed at the upper level, the lower level is presented as a SUE assignment problem. In this context, it is assumed that the local authority tries to simultaneously optimize traffic signal timings by maximizing the reserve capacity and minimizing PI by taking into account the users’ reactions in SUE manner. This mutual interaction between users and local authority is presented by using biobjective BLPM. In summary, the major contributions to the current literature made in this paper are as follows. (i) A biobjective BLPM is proposed. The upper level maximizes the reserve capacity of a road network and minimizes PI. The lower level solves SUE problem by using the path flow estimator (PFE) which follows logit route choice model. (ii) A weightedsum method is applied to find Pareto optimal solutions of the biobjective problem by normalizing the objective functions. (iii) A biobjective BLPM is applied to a mediumsized signalized road network.
The rest of this paper is organized as follows. The next section is about the problem formulation. Section 3 presents detailed explanations of the proposed bobjective BLPM. The numerical experiments are given in Section 4. The conclusions and future directions are given in the last section.
2. Problem Formulation
Considering that the local authorities or agencies are responsible for traffic planning, data collection, sign and signal operations, etc., they desire robust road networks in terms of congestion and delay even in case of sudden spikes in travel demand. On the other hand, traffic flow is distributed on a road network based on route choice behaviour of users that aim to complete their travels within the shortest possible travel time. Therefore, maximizing reserve capacity and minimizing delay provide reasonable benefits for local authorities and road users, respectively. However, maximizing reserve capacity may reduce network performance by increasing delay because of growing travel demand. In this study, this tradeoff is investigated by defining a biobjective BLPM. At the upper level of the proposed model, combined reserve capacity maximization and delay minimization problem is solved, while SUE assignment problem is performed using PFE at the lower level. Considering a road network with a set of OD pairs, , a set of directed links, , a set of paths, , and a set of nodes, , the proposed biobjective BLPM is presented as follows:subject towhere (1) represents the problem of reserve capacity maximization and (2) is related to delay minimization problem. In addition, is the demand matrix multiplier, is delay on link , , is the number of stops on link , and are weighting factors for and for on link , respectively, is stop penalty factor, is vector of equilibrium link flows, is common cycle time, and are possible bounds for common cycle time, is stage green time, is the minimum green time, is intergreen time, is offset, is vector of signal timings, is feasible region for signal timings, and is the number of stages. In addition to signal timings related constraints, (4) represents link capacity constraint in which is the saturation flow, is the capacity, and is equilibrium flow on link . Link flows can be obtained by solving SUE problem using (5)(6) at the lower level [50]. subject towhere is the vector of travel demand, and represent vectors of link and path travel times for the given vector of link flows , respectively, is the vector of path flows, represents the link/path incidence matrix where if link is on path , and otherwise , and is the OD/path incidence matrix . The SUE problem given in (5) can be solved by the PFE [51]. The PFE is a logit based traffic assignment tool and its most important advantage is that it does not require path enumeration. Therefore it is most preferable option for researches for the solution of SUE assignment [34, 35, 41, 50].
Since there are two objectives in the upper level in (1)(4), which are mutually interacted, it is not easy to find optimal solution between two objectives. Thus, a suitable solution has to be found among Pareto optimal solutions. Although there are several efficient methods to find the set of Pareto optimal solutions of biobjective problems, the weightedsum method can be used as an efficient tool to convert a biobjective problem into singleobjective by using convex linear combination of objectives [52]. On the other hand, since each singleobjective problem (delay minimization and reserve capacity maximization) has a different unit, they cannot be directly added. To overcome this issue, they can be normalized before applying the weightedsum method as given in where can be stated as an upper limit of reserve capacity maximization problem which is the solution of the singleobjective problem given in (1) subject to (3)(6). On the other hand, is the lower limit of delay minimization problem and can be found by solving (2) subject to (3) and (5)(6). It should be emphasized that the use of converts the reserve capacity maximization problem into a minimization problem. Thus, objective functions and can be stated as minimization problems. After normalizing the objective functions, they are considered as scalar, and thus there is no further obstacle for applying the weightedsum method to convert the biobjective problem into a singleobjective problem as given in (9) subject to (3)(6).For applying weightedsum method we have used a weighting factor in order to present relative importance of one objective to other. In other words, the value of is specified to find the balance between solutions of network reserve capacity maximization and delay minimization problems. It should be noted that a larger value of implies that the local authority should be more attentive to the improvement of network reserve capacity. On the other side, the authority prioritizes to minimize network overall delay in comparison with reserve capacity maximization when a smaller value of is used. Choosing a suitable value of weighting factor can be performed by generating a uniform set of weights as , where is a positive integer number [53].
3. Model Development
In this study, the proposed biobjective BLPM is developed based on differential evolution (DE) optimization framework. As known, DE is referred as a robust and simple metaheuristic algorithm proposed by [54]. On the other hand, several studies in the relevant literature support this claim by introducing the advantages of DE in solving different complex optimization problems. Nolle et al. [55] have made a comparison of selforganizing migration algorithm (SOMA) with simulated annealing (SA) and DE algorithms for automated waveform tuning. Although a slight difference has been observed between performances of DE and SOMA algorithms, they produced similar computational performance in this application outperforming SA. In the study by Vincenzi et al. [56] the performance of DE algorithm has been compared with Coupled Local Minimizers (CLM) for solving dynamic damage detection problem. Results revealed that DE outperforms CLM method especially in case the number of decision variables is increased. Similarly, a comparison of DE algorithm with superior metaheuristic algorithms has been performed on different benchmark functions by Civicioglu and Besdok [57]. According to the results, DE produces more robust results compared with other metaheuristic algorithms. Deb et al. [58] have compared the performances of DE, PSO, and GA in solving design problem of circularly polarized microstrip antennas. DE based algorithm outperformed GA and PSO according to the numerical applications. Recently, Dezelak et al. [59] tested the performances of DE and PSO algorithms for power plants modelling. Results have shown that DE produces better results in reasonable computation time than PSO. In the light of short review of DE algorithm, it has been chosen for solving the biobjective BLPM by considering its successful applications.
3.1. Reserve Capacity Maximization Model
As previously mentioned, values of the maximum OD multiplier, , and the minimum network performance index, , are required to set up the singleobjective model given in (9) for obtaining the Pareto optimal solutions of the biobjective optimization problem. In this context, both problems are handled with a BLPM based on the DE solution framework to obtain and , in which the network OD multiplier and PI are optimized at the upper level. On the other hand, SUE assignment problem is solved to determine equilibrium link flows at the lower levels of both models. In this context, maximization of the OD multiplier is formulated as a minimization problem subject to (3)(6) as follows:where is a constant penalty weighting factor. Fivestep solution framework is illustrated in Figure 1.
Step 1. The objective function representing the maximization of network reserve capacity given in (10) and the possible bounds for decision variables, which are common cycle time, stage green times, and OD multiplier, are primarily initialized. Subsequently, network parameters, namely, saturation flows, freeflow travel times, and travel demand, are represented. Initialization of three DE algorithm parameters is also carried out at this step. The first one is population size, , and the next one is mutation factor, F, which is used to create mutant vector. Crossover rate () is the third parameter, which is employed in creating the trial vector [54]. At the end of this step, generation number (gen) is set to 1, iteration number () is set to 0, and maximum number of generations (maxgen) is identified as the stopping criterion.
Step 2. At this step, initial solution vectors including OD multipliers and signal timings are generated considering preset bounds of decision variables. Subsequently, objective function values are calculated and represented as given in where is the stage green time of intersection and is the number of stages at intersection ( and ). In this study, initial values of the OD multipliers and signal timing parameters are obtained in the following way.
(i) OD multiplier is generated between 0 and possible upper bound as given in
(ii) Cycle time is randomly generated between the possible bounds as given in
(iii) Stage green times at intersections are generated between and as shown in Subsequently, generated green timing values are revised since all green plus intergreen times must be equal to the cycle time for a given intersection. This procedure is carried out by As can be seen in (11) the equilibrium link flows, , are required to calculate the objective function values. In this context, the SUE assignment problem, which was given in (5), is solved by the PFE traffic assignment tool, the pseudo code of which is given in Algorithm 1.
Define dispersion parameter  
Repeat  
Update link costs  
For each path  
Calculate new path costs  
Next  
For each path  
Calculate new path flows  
Next  
For each link  
Calculate new link flows  
Next  
Until no new path and link flows converged 
In Algorithm 1, represents dispersion parameter. Note that if is zero then the path costs are not fully taken into account in path choice which means even highest cost paths can be chosen by road users. On the contrary, drivers become sensitive to cost in path choice while tends to infinity.
Step 3. The mutation is carried out at this step by using three different, randomly chosen solution vectors which are different from target vector (i.e., solution vector) and a mutant vector, , is created as given in where , , and are indices of randomly chosen solution vectors.
Step 4. At this step, crossover is applied by choosing each member of the trial vector, , from the target or the mutant vectors with the probabilities of or 1, respectively, as given in where the condition of provides that target and trial vectors are definitely different from each other. After determining the members of the trial vector, their related SUE link flows and the objective function value are calculated.
Step 5. At the last step, the target vector, , for the next generation is selected by comparing the objective function values of trial and target vectors as given in
The DE process is repeated until the maximum number of generations, maxgen, is reached.
3.2. Delay Minimization Model
Following the representation of the reserve capacity maximization model, flowchart of the delay minimization model for the solution of (2) is illustrated in Figure 2.
Step 1. At this step, the objective function given in (2) subject to (3) and (5)(6) representing the minimization of network, , and the possible bounds for signal timing variables are primarily initialized. Note that, unlike the reserve capacity maximization model, OD multiplier is set to 1, which means that the base travel demand is loaded into the network. The remaining of this step is the same as the first step of the reserve capacity maximization model.
Step 2. Initial signal timings are generated between preset bounds based on (13)(15). Moreover, the offset variables are randomly generated as shown in
As can be seen in (2) the equilibrium link flows are required to calculate network value. Thus, the SUE assignment problem, which was given in (5), is solved by the PFE for each solution vector in the population pool. Subsequently, TRANSYT7F is employed to calculate the related values and the target vectors, namely, the initial solutions, are stored in the memory as given in
Step 3. The mutation is carried out to create a mutant vector, , which is given in
Step 4 . At this step, crossover is applied to determine the members of the trial vector, , as given in
Step 5. Selection is carried out in terms of the objective function values of trial vector, , and its target vector, , in the following way.
The DE process is repeated until the maximum number of generations, maxgen, is reached.
3.3. Combined Model
Once the values of and are obtained, the combined reserve capacity maximization and delay minimization problem is solved using singleobjective BLPM in (9) subject to (3)(6) by taking a particular weighting factor into account. The corresponding stepwise procedure is given below.
Step 1. Generate the initial signal timing vectors, , and OD multipliers, , with given possible bounds and calculate objective function values, which is given in (9), as many as population size.
Step 2. Set and .
Step 3. Set . Then, generate a new signal timing vector, , and a new OD multiplier, , based on mutation and crossover operations of DE.
Step 4. Calculate based on the OD multiplier, , which has been generated in the previous step, and as given in (7).
Step 5. Run PFE to carry out a SUE assignment for estimating link flows based on the OD multiplier, , and the signal timing vector, , which have been generated in tep 3. Subsequently, run TRANSYT7F to obtain the network value and then calculate based on and as given in (8).
Step 6. Calculate objective function as given in (9).
Step 7. Compare the solution vector in the population pool and newly generated solution vector in terms of their objective function values. Subsequently, keep the one with a better objective function value and discard the other.
Step 8. Go to tep 3 if .
Step 9. Set and and then go to tep 3 if .
Step 10. Terminate the DE process and output the optimized signal timing vector and OD multiplier.
4. Numerical Applications
In this section, two signalized benchmark road networks are chosen to investigate the performance of the proposed biobjective BLPM. Firstly, a twojunction network is used to show the effect of weighting on the balance of solution of reserve capacity maximization and delay minimization problems. In other words, it is proposed to find the set of Pareto optimal solutions of the biobjective problem. With this application, it is also aimed at investigating how the network performance is affected by applying the biobjective BLPM depending upon the weighting factor . Secondly, the effectiveness of the proposed model is investigated on a mediumsized road network presented by Allsop and Charlesworth [60].
4.1. TwoJunction Network
In order to test the capability of the proposed model in solving the biobjective problem, it has been firstly applied to the twojunction network taken from Ceylan [61]. Layout of the network and the input data are given in Figure 3 and Table 1, respectively. Travel demand between the single OD pair is taken as 1500 veh/hr.
 
where is free flow travel time in seconds on link . 
The intergreen time, , between stages was selected as 5 sec. Constraints for signal timing variables are set as given below:Setting of mutation factor , crossover rate , and population size is crucial for the performance of DE algorithm. Storn and Price [54] recommend that the ranges of and may be chosen for and , respectively. In this context, an analysis is performed for 18 parameter combinations where and are varied from 0.5 to 1.0 and from 0.8 to 1.0 in steps of 0.1, respectively. 30 independent runs are performed for every parameter combination for solving delay minimization problem and the corresponding results are given in Table 2. is set as 20 based on the recommendation given in Rönkkönen et al. [62] and maxgen is set as 100. Note that the average computational time for complete run of the proposed model is about 28 minutes. It means that each generation takes about 16.8 seconds of CPU. The proposed model has been executed in MATLAB programming and performed on PC with Intel Core i7 2.10 GHz, RAM 8 GB.

It can be seen in Table 2 that value is obtained as 7.02 with 7 different parameter combinations. In addition to that, the setting resulted in the lowest mean and standard deviation values which are 7.06 and 0.04, respectively. Therefore, this parameter combination is selected for further applications. Corresponding results for at twojunction network are presented in Tables 3 and 4. As can be seen in Table 3, the demand was almost equally shared by links 1 and 4. Other links in the network do not carry flow according to the solution for .


Once the network has been calculated, the reserve capacity maximization problem can be solved for twojunction network assuming that the bounds for OD multiplier are (0,5). According to the analysis, was found as 2.1370 which is the solution of the singleobjective problem given in (10). Other corresponding results obtained with the solution are given in Tables 5 and 6, respectively.


As shown in Table 5, the values of degree of saturation on the links indicate that none of the links is oversaturated although the OD demand was multiplied by 2.1370. In other words, flows do not exceed their corresponding link capacities. In the case that OD demand matrix was multiplied by , the cycle time was found as 100 which is its upper bound for this example. This means that the common cycle time was increased, together with the OD multiplier to overcome the increased demand which was more than double of the base demand on the network.
Once the values of and are obtained, there is no further obstacle to solve the combined reserve capacity maximization and delay minimization problem by using (9). Optimized values of the objective functions, , , and , are given for different values of weighting factors in Table 7. It can be seen in Table 7 that, in case is zero, any improvement in the network reserve capacity is not allowed. Thus, the value of becomes equal to since is 1.0 which represents the base demand. On the other hand, the value of equals 1.0 indicating that the optimized value of is equivalent to which is the upper bound of . In brief, the authority is fully concerned with the delay minimization and no additional travel demand is added to the network. On the contrary, when equals 1.0, the authority is not further concerned with the delay minimization and tries to maximize the reserve capacity in the network and consequently the value of equals 1.0 which means that becomes equal to .

In addition, the graphical representation for the values of and regarding different values of is given in Figure 4 which also shows the Pareto optimal solutions for the objective functions, and . It can be emphasized that and are in conflict with each other; that is, shows a decreasing tendency while has a tendency to increase. Thus, regarding the dilemma between reserve capacity maximization and delay minimization problems, the authority can choose one among the Pareto optimal solutions according to its own preferences. It can be seen that there is no significant change in the values of and when the value of increases from 0 to 0.5. Afterwards, the value of tends to increase in case is higher than 0.5. Likewise, it can be clearly observed that the value of dramatically increases when is higher than 0.8. This indicates that the authority pays more attention to the network reserve capacity maximization than delay minimization.
Depending on the variation of values between 0 and 1, values of and , their changes, and cumulative changes are given in Table 8. As can be seen in Table 8 and increase about 8% and 9%, respectively, while the objective weight increases from 0 to 0.4. That is, approximately 2% change is observed for both and for every 0.1 change in . However, when is selected as the solution point, change in and is realized to be about 4% and 6%, respectively. This means that the users’ benefit is badly affected by growing travel demand. Examining the rest of Table 8, it can be seen that there are two significant breaking points in terms of the changes in and . The first one is ; that is, the increases in and are about 4% and 21%, respectively. The second point is ; that is, and increases are about 16% and 118%, respectively. This result can be clearly seen in Figure 5 which represents the change of with the OD multiplier .

For further explanation of the changes of and , we present Figure 6 in which their cumulative changes and corresponding change rates are considered with the objective weight of . Bars with three different colours in Figure 6, which represent the cumulative change rates , reveal three different zones where each one has specific values of cumulative change rates. Therefore, the local authority can make decision about by considering the first points of those zones. Two breaking points, which are and as explained based on Table 8, are clearly seen in the figure. When equals 0.9, a more significant increase in cumulative value is observed than that for . Therefore, may be considered as the best objective weight that provides an OD multiplier and signal timing configuration since all other Pareto optimal solutions proportionally reduce users’ benefit more than the gain in network reserve capacity.
4.2. Allsop & Charlesworth’s Network
Allsop & Charlesworth’s test network is chosen as second numerical example since it has widely been used in the studies concerning transportation network design [35, 36, 41, 42, 45]. Thus, we aim to give readers a chance to compare some results drawn by some previous studies. The representation of the test road network and its stage plans are given in Figure 7. The example network consists of 23 links and 21 signal setting variables at six junctions. The input data and relevant OD demand matrix are given in Tables 9 and 10, respectively.


Constraints for each signal timing variable are set as given below:The intergreen time, , between stages was selected as 5 sec. In addition, bounds for OD multiplier were assumed to be (0,3). It should be pointed out that the reason for selecting these upper and lower bounds for OD multiplier on this example is to reduce the feasible search space and thus to decrease the computational burden of proposed solution algorithm. Considering the number of decision variables for Allsop & Charlesworth’s network, and maxgen are selected as 40 and 500, respectively. Firstly, the RCMP was solved by using BLPM given in Figure 1, and was found as 1.29 for Allsop & Charlesworth’s network which is the solution of the singleobjective problem given in (10). Subsequently, was determined as 341.10 by solving the delay minimization problem by using the proposed model given in Figure 2. Similarly to the Pareto efficiency analysis performed for twojunction network, the best values of and with their changes and cumulative changes after 10 independent runs are given in Table 11 for values between 0 and 1. The average computational time for complete run of the proposed model is about 4.7 hours. As can be seen in Table 11, and increased on average by 1.0% and 4.8%, respectively, for every 0.1 change in while it increases from 0 to 0.8. On the other hand, while equals 0.9, the changes of and are about 6 and 8 times higher than their values in the case for . This means that the users’ benefit is badly affected by significant increase in network . Therefore, may be considered as the best objective weight that provides a balance between solutions of reserve capacity maximization and delay minimization problems. This result can be seen in Figure 8 which represents the change of with the OD multiplier .

Resulting SUE link flows, degrees of saturations, and signal timings for values of 0.0, 0.8, and 1.0 are given in Tables 12 and 13, respectively.


When the value of is 0, that is, only the delay minimization is taken into account by the authority, the degrees of saturation of SUE links converge to the lowest values as shown in Table 12. On the contrary, in case of , which represents that the authority is totally concerned with the maximizing of reserve capacity, has substantially increased about 3.2 times, namely, 1084 for the maximum OD multiplier of 1.29 in comparison with the case of as shown in Table 13. Moreover, although the OD demand matrix is increased about 30%, none of the links is oversaturated which means that their degrees of saturation are less than 1.0. For , and are found as 1.08 and 495.40, respectively. It means that when equals 0.8, the road network can accommodate about 8% more travel demand with ensuring links do not exceed their capacities while increases about 45%. This is because equals 0.8 and the authority does not pay attention to the delay minimization as much as the maximization of reserve capacity problem. On the other side, the common cycle time for Allsop & Charlesworth’s network was found as 150, 138, and 150 sec for values of 0, 0.8, and 1.0, respectively.
5. Conclusions and Future Studies
This study deals with the simultaneous solution of reserve capacity maximization and delay minimization problems by optimizing traffic signal timings. For this purpose, a biobjective BLPM is developed and solved by implementing a DE based heuristic solution algorithm. We introduce the weightedsum method to solve the proposed biobjective model and present an approach to make the objective functions of the biobjective model dimensionless. In the upper level, it is aimed to maximize reserve capacity and to minimize delay in a road network, whereas at the lower level the user’s reactions are determined by solving stochastic traffic assignment problem.
Two numerical applications are performed in order to show the applicability of the proposed model. Firstly, the proposed model is applied to a twojunction network. In this application, the cumulative change of value is extremely high when is set to 1.0. It means that the authority is totally concerned with the reserve capacity maximization without taking delay minimization into account. Thus, the value spikes up to about 100. For the twojunction network, is selected as a suitable objective weight that provides a balance between conflicting issues, namely, reserve capacity maximization and delay minimization. Secondly, the proposed model is applied to Allsop & Charlesworth’s road network in order to show its effectiveness on a mediumsized signalized road network. According to the Pareto efficiency analysis for this network, it is found that the case where equals 0.8 is also suitable since all other Pareto solutions proportionally reduce users’ benefit more than the gain in network reserve capacity. Moreover, results indicate that there are no oversaturated links in the network even in the case where equals 1.0, which means that the OD multiplier is maximized.
In conclusion, results show that the use of different α values gives an opportunity to responsible authorities to make a suitable choice for the benefit of road users between two conflicting issues. For future studies, dynamic traffic assignment rather than static assignment may be considered in order to reflect the effect of timevarying flows on the solution of the proposed model.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there is no conflict of interest regarding the publication of this paper.
References
 F. V. Webster and B. M. Cobbe, “Traffic signal,” Road Research Technical Paper 56, HMSO, London, UK, 1966. View at: Google Scholar
 R. E. Allsop, “Estimating the traffic capacity of a signalized road junction,” Transportation Research, vol. 6, no. 3, pp. 245–255, 1972. View at: Publisher Site  Google Scholar
 S. Yagar, “Capacity of a signalized road junction: critique and extensions,” Transportation Research, vol. 8, no. 2, pp. 137–147, 1974. View at: Publisher Site  Google Scholar
 S. Yagar, “Addressing errors and omissions in paper on intersection capacity maximization,” Transportation Research Part B: Methodological, vol. 19, no. 1, pp. 81–84, 1985. View at: Publisher Site  Google Scholar
 S. C. Wong and H. Yang, “Reserve capacity of a signalcontrolled road network,” Transportation Research Part B: Methodological, vol. 31, no. 5, pp. 397–402, 1997. View at: Publisher Site  Google Scholar
 H. Yang, M. G. H. Bell, and Q. Meng, “Modeling the capacity and level of service of urban transportation networks,” Transportation Research Part B: Methodological, vol. 34, no. 4, pp. 255–275, 2000. View at: Publisher Site  Google Scholar
 H. Yang and J. Y. Wang, “Travel time minimization versus reserve capacity maximization in the network design problem,” Transportation Research Record, vol. 1783, no. 1, pp. 17–26, 2002. View at: Publisher Site  Google Scholar
 G. Ziyou and S. Yifan, “A reserve capacity model of optimal signal control with userequilibrium route choice,” Transportation Research Part B: Methodological, vol. 36, no. 4, pp. 313–323, 2002. View at: Publisher Site  Google Scholar
 Y. E. Ge, H. M. Zhang, and W. H. Lam, “Network reserve capacity under influence of traveler information,” Journal of Transportation Engineering, vol. 129, no. 3, pp. 262–270, 2003. View at: Publisher Site  Google Scholar
 H. Ceylan and M. G. H. Bell, “Reserve capacity for a road network under optimized fixed time traffic signal control,” Journal of Intelligent Transportation Systems: Technology, Planning, and Operations, vol. 8, no. 2, pp. 87–99, 2004. View at: Publisher Site  Google Scholar
 A. Chen, P. Chootinan, and S. Wong, “New reserve capacity model of signalcontrolled road network,” Transportation Research Record, vol. 1964, pp. 35–41, 2006. View at: Publisher Site  Google Scholar
 S.W. Chiou, “Reserve capacity of signalcontrolled road network,” Applied Mathematics and Computation, vol. 190, no. 2, pp. 1602–1611, 2007. View at: Publisher Site  Google Scholar
 S. Chiou, “An efficient search algorithm for road network optimization,” Applied Mathematics and Computation, vol. 201, no. 12, pp. 128–137, 2008. View at: Publisher Site  Google Scholar
 S. Chiou, “Optimization of limited network capacity with toll settings,” Information Sciences, vol. 179, no. 12, pp. 109–119, 2009. View at: Publisher Site  Google Scholar
 S. Chiou, “A hybrid approach for optimal design of signalized road network,” Applied Mathematical Modelling, vol. 32, no. 2, pp. 195–207, 2008. View at: Publisher Site  Google Scholar
 S. Chiou, “Optimization for signal setting problems using nonsmooth techniques,” Information Sciences, vol. 179, no. 17, pp. 2985–2996, 2009. View at: Publisher Site  Google Scholar
 E. Miandoabchi and R. Z. Farahani, “Optimizing reserve capacity of urban road networks in a discrete Network Design Problem,” Advances in Engineering Software, vol. 42, no. 12, pp. 1041–1050, 2011. View at: Publisher Site  Google Scholar
 E. Miandoabchi, F. Daneshzand, W. Y. Szeto, and R. Zanjirani Farahani, “Multiobjective discrete urban road network design,” Computers & Operations Research, vol. 40, no. 10, pp. 2429–2449, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 S. Chiou, “Optimal signalsetting for road network with maximum capacity,” Information Sciences, vol. 273, pp. 287–303, 2014. View at: Publisher Site  Google Scholar
 J. Wang, W. Deng, and J. Zhao, “Road network reserve capacity with stochastic user equilibrium,” Transport, vol. 30, no. 1, pp. 103–116, 2015. View at: Publisher Site  Google Scholar
 H. Xiao, J. Gao, and Z. Zou, “Reserve capacity model based on variable demand for landuse development control,” Transportation Planning and Technology, vol. 40, no. 2, pp. 199–212, 2016. View at: Publisher Site  Google Scholar
 F. Han and L. Cheng, “Stochastic user equilibrium model with a tradable credit scheme and application in maximizing network reserve capacity,” Engineering Optimization, vol. 49, no. 4, pp. 549–564, 2016. View at: Publisher Site  Google Scholar
 O. Baskan and C. Ozan, “Reserve capacity model for optimizing traffic signal timings with an equity constraint,” in Highway Engineering, IntechOpen, 2017. View at: Publisher Site  Google Scholar
 X. Li and J. Sun, “Multiobjective optimal predictive control of signals in urban traffic network,” Journal of Intelligent Transportation Systems: Technology, Planning, and Operations, pp. 1–19, 2019. View at: Publisher Site  Google Scholar
 McTrans Center, “TRANSYT7F Release 11.3 Users Guide,” University of Florida, Gaineville, Florida, USA, 2008. View at: Google Scholar
 F. V. Webster, “Traffic signal settings,” Road Research Technical Paper 39, HMSO, London, UK, 1958. View at: Google Scholar
 R. E. Allsop, “Delayminimizing settings for fixedtime traffic signals at a single road junction,” Journal of the Institute of Mathematics and Its Applications, vol. 8, pp. 164–185, 1971. View at: Publisher Site  Google Scholar  MathSciNet
 S. Wong, “Derivatives of the performance index for the traffic model from TRANSYT,” Transportation Research Part B: Methodological, vol. 29, no. 5, pp. 303–327, 1995. View at: Publisher Site  Google Scholar
 B. Heydecker, “A decomposition approach for signal optimisation in road networks,” Transportation Research Part B: Methodological, vol. 30, no. 2, pp. 99–114, 1996. View at: Publisher Site  Google Scholar
 S. C. Wong, “Groupbased optimisation of signal timings using the transyt traffic model,” Transportation Research Part B: Methodological, vol. 30, no. 3, pp. 217–244, 1996. View at: Publisher Site  Google Scholar
 S. C. Wong, “Groupbased optimisation of signal timings using parallel computing,” Transportation Research Part C: Emerging Technologies, vol. 5, no. 2, pp. 123–139, 1997. View at: Publisher Site  Google Scholar
 S. C. Wong, W. T. Wong, C. M. Leung, and C. O. Tong, “Groupbased optimization of a timedependent TRANSYT traffic model for area traffic control,” Transportation Research Part B: Methodological, vol. 36, no. 4, pp. 291–312, 2002. View at: Publisher Site  Google Scholar
 M. Girianna and R. F. Benekohal, “Application of genetic algorithms to generate optimum signal coordination for congested networks,” in Proceedings of the 7th International Conference on Applications of Advanced Technology in Transportation, pp. 762–769, 2002. View at: Google Scholar
 H. Ceylan and M. G. H. Bell, “Traffic signal timing optimisation based on genetic algorithm approach, including drivers' routing,” Transportation Research Part B: Methodological, vol. 38, no. 4, pp. 329–342, 2004. View at: Publisher Site  Google Scholar
 H. Ceylan and M. G. Bell, “Genetic algorithm solution for the stochastic equilibrium transportation networks under congestion,” Transportation Research Part B: Methodological, vol. 39, no. 2, pp. 169–185, 2005. View at: Publisher Site  Google Scholar
 H. Ceylan, “Developing combined genetic algorithm hillclimbing optimization method for area traffic control,” Journal of Transportation Engineering, vol. 132, no. 8, pp. 663–671, 2006. View at: Publisher Site  Google Scholar
 J. Chen and L. Xu, “Roadjunction traffic signal timing optimization by an adaptive particle swarm algorithm,” in Proceedings of the 9th International Conference on Control, Automation, Robotics and Vision, vol. 1–5, pp. 1103–1109, 2006. View at: Google Scholar
 C. Dan and G. Xiaohong, “Study on intelligent control of traffic signal of urban area and microscopic simulation,” in Proceedings of the 8th International Conference of Chinese Logistics and Transportation Professionals, Logistics: The Emerging Frontiers of Transportation and Development in China, pp. 4597–4604, 2008. View at: Google Scholar
 Z. Li, “Modeling arterial signal optimization with enhanced cell transmission formulations,” Journal of Transportation Engineering, vol. 137, no. 7, pp. 445–454, 2011. View at: Publisher Site  Google Scholar
 Y. Liu and G.L. Chang, “An arterial signal optimization model for intersections experiencing queue spillback and lane blockage,” Transportation Research Part C: Emerging Technologies, vol. 19, no. 1, pp. 130–144, 2011. View at: Publisher Site  Google Scholar
 H. Ceylan, “A hybrid harmony search and TRANSYT hill climbing algorithm for signalized stochastic equilibrium transportation networks,” Transportation Research Part C: Emerging Technologies, vol. 25, pp. 152–167, 2012. View at: Publisher Site  Google Scholar
 M. Dell'Orco, O. Baskan, and M. Marinelli, “A Harmony Search Algorithm approach for optimizing traffic signal timings,” PROMET  Traffic&Transportation, vol. 25, no. 4, pp. 349–358, 2013. View at: Publisher Site  Google Scholar
 M. Dell’Orco, Ö. Başkan, and M. Marinelli, “Artificial bee colonybased algorithm for optimising traffic signal timings,” in Soft Computing in Industrial Applications, Advances in Intelligent Systems and Computing, V. Snášel, P. Krömer, M. Köppen, and G. Schaefer, Eds., vol. 223, pp. 327–337, Springer, Berlin, Germany, 2014. View at: Google Scholar
 F. Zhu, H. A. Aziz, X. Qian, and S. V. Ukkusuri, “A junctiontree based learning algorithm to optimize network wide traffic control: a coordinated multiagent framework,” Transportation Research Part C: Emerging Technologies, vol. 58, pp. 487–501, 2015. View at: Publisher Site  Google Scholar
 C. Ozan, O. Baskan, S. Haldenbilen, and H. Ceylan, “A modified reinforcement learning algorithm for solving coordinated signalized networks,” Transportation Research Part C: Emerging Technologies, vol. 54, pp. 40–55, 2015. View at: Publisher Site  Google Scholar
 E. Christofa, K. Ampountolas, and A. Skabardonis, “Arterial traffic signal optimization: a personbased approach,” Transportation Research Part C: Emerging Technologies, vol. 66, pp. 27–47, 2016. View at: Publisher Site  Google Scholar
 L. Zhang, Z. Song, X. Tang, and D. Wang, “Signal coordination models for long arterials and grid networks,” Transportation Research Part C: Emerging Technologies, vol. 71, pp. 215–230, 2016. View at: Publisher Site  Google Scholar
 S. Srivastava and S. K. Sahana, “Nested hybrid evolutionary model for traffic signal optimization,” Applied Intelligence, vol. 46, no. 1, pp. 113–123, 2017. View at: Publisher Site  Google Scholar
 S. Memoli, G. E. Cantarella, S. de Luca, and R. D. Pace, “Network signal setting design with stage sequence optimisation,” Transportation Research Part B: Methodological, vol. 100, pp. 20–42, 2017. View at: Publisher Site  Google Scholar
 M. G. H. Bell and Y. Iida, Transportation Network Analysis, John Wiley and Sons, Chichester, UK, 1997.
 M. G. H. Bell and C. M. Shield, “A loglinear model for path flow estimation,” in Proceedings of the 4th International Conference on the Applications of Advanced Technologies in Transportation Engineering, Y. J. Stephanedes and F. Filippi, Eds., pp. 695–699, Capri, Italy, 1995. View at: Google Scholar
 Z.C. Li and X.Y. Ge, “Traffic signal timing problems with environmental and equity considerations,” Journal of Advanced Transportation, vol. 48, no. 8, pp. 1066–1086, 2014. View at: Publisher Site  Google Scholar
 E. Ulungu, J. Teghem, P. Fortemps, and D. Tuyttens, “MOSA method: a tool for solving multiobjective combinatorial optimization problems,” Journal of MultiCriteria Decision Analysis, vol. 8, no. 4, pp. 221–236, 1999. View at: Publisher Site  Google Scholar
 R. Storn and K. Price, “Differential evolution—a simple and efficient heuristic for global optimization over continuous spaces,” Journal of Global Optimization, vol. 11, no. 4, pp. 341–359, 1997. View at: Publisher Site  Google Scholar
 L. Nolle, I. Zelinka, A. Hopgood, and A. Goodyear, “Comparison of an selforganizing migration algorithm with simulated annealing and differential evolution for automated waveform tuning,” Advances in Engineering Software, vol. 36, no. 10, pp. 645–653, 2005. View at: Publisher Site  Google Scholar
 L. Vincenzi, G. De Roeck, and M. Savoia, “Comparison between coupled local minimizers method and differential evolution algorithm in dynamic damage detection problems,” Advances in Engineering Software, vol. 65, pp. 90–100, 2013. View at: Publisher Site  Google Scholar
 P. Civicioglu and E. Besdok, “A conceptual comparison of the Cuckoosearch, particle swarm optimization, differential evolution and artificial bee colony algorithms,” Artificial Intelligence Review, vol. 39, no. 4, pp. 315–346, 2013. View at: Publisher Site  Google Scholar
 A. Deb, J. S. Roy, and B. Gupta, “Performance comparison of differential evolution, particle swarm optimization and genetic algorithm in the design of circularly polarized microstrip antennas,” IEEE Transactions on Antennas and Propagation, vol. 62, no. 8, pp. 3920–3928, 2014. View at: Publisher Site  Google Scholar
 P. Bracinik, K. Dezelak, A. Otcenasova, and M. Höger, “Comparison between the particle swarm optimisation and differential evolution approaches for the optimal proportional–integral controllers design during photovoltaic power plants modelling,” IET Renewable Power Generation, vol. 10, no. 4, pp. 522–530, 2016. View at: Publisher Site  Google Scholar
 R. E. Allsop and J. A. Charlesworth, “Traffic in a signalcontrolled road network: an example of different signal timings inducing different routeings,” Traffic Engineering & Control, vol. 18, no. 5, pp. 262–265, 1977. View at: Google Scholar
 H. Ceylan, A genetic algorithm approach to the equilibrium network design problem [Ph.D. Thesis], University of Newcastle upon Tyne, Newcastle, UK, 2002.
 J. Ronkkonen, V. Kukkonen, and K. V. Price, “Realparameter optimization with differential evolution,” in Proceedings of the IEEE Congress on Evolutionary Computation, pp. 506–513, Edinburgh, Scotland, September 2005. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Ozgur Baskan et al. 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.