Research Article  Open Access
Robust Evaluation for Transportation Network Capacity under Demand Uncertainty
Abstract
As more and more cities in worldwide are facing the problems of traffic jam, governments have been concerned about how to design transportation networks with adequate capacity to accommodate travel demands. To evaluate the capacity of a transportation system, the prescribed origin and destination (OD) matrix for existing travel demand has been noticed to have a significant effect on the results of network capacity models. However, the exact data of the existing OD demand are usually hard to be obtained in practice. Considering the fluctuation of the real travel demand in transportation networks, the existing travel demand is represented as uncertain parameters which are defined within a bounded set. Thus, a robust reserve network capacity (RRNC) model using min–max optimization is formulated based on the demand uncertainty. An effective heuristic approach utilizing cutting plane method and sensitivity analysis is proposed for the solution of the RRNC problem. Computational experiments and simulations are implemented to demonstrate the validity and performance of the proposed robust model. According to simulation experiments, it is showed that the link flow pattern from the robust solutions to network capacity problems can reveal the probability of high congestion for each link.
1. Introduction
The capacity of transportation network reflects the supply ability of its infrastructure and service to the travel demand which is generated from the zones covered by the transportation system in a specific period. For many years, transportation planners and managers wanted to understand how many trips can be accommodated at the most by the current or designed network in a certain period of time. This need is more necessary in those developing regions which are confronted with rapid growth of private vehicles and increased urban congestion. Meanwhile, the researchers made a longterm effort to model and estimate the maximum throughput of transportation networks. The achievements include maxflow mincut theorem [1], incremental assignment approach [2], and later bilevel programming models [3–5].
For the network capacity model, the most popular formulation in passenger transportation system is the bilevel model, which maximizes the traffic flows under the equilibrium constraints. Wong and Yang [3] first incorporated the reserve capacity concept into a traffic signal control network. The reserve capacity is defined as the largest multiplier applied to a given OD demand matrix without violating capacity constraints, so the solution is significantly affected by the predetermined OD matrix. Ziyou and Yifan [6] extended the reserve capacity model by considering OD specific demand multipliers, and all demand multipliers should be ensured not lower than a predetermined minimum value. In order to avoid assuming that all OD flows increase in a same rate, another concept of ultimate capacity was proposed [5]. But it assumes that the OD distribution is totally variable, which may produce unrealistic results that cause the trip productions at some origins below their current levels. Furthermore, Yang et al. [4] suggested that the new increased OD demand pattern should be variable in both level and distribution, while the current travel demand is fixed. Later, Yang’s model was also referred to as the practical capacity by Kasikitwiwat and Chen [5]. In summary, although unrealistic, the reserve capacity model is more easytouse and has been adopted widely in many researches [7–10]. The ultimate capacity and practical capacity model are more practical but have more parameters to be calibrated when applied, and the formulated models are still difficult to solve [11].
While the deterministic network capacity problem has been explored extensively, few studies have investigated the issue of uncertainties in demand data associated with this problem. The ultimate capacity and practical capacity model are only concerned with the uncertainties related to the new increased travel demand by using combined models [5], while the uncertainties in the current (or existing) demand are not considered. In reality, travel demands in transportation system are always fluctuant day by day, even hour by hour. Besides, errors of survey data also affect the accuracy of the existing OD matrix. As a consequence, the existing travel demands are usually difficult to be obtained in actual transportation projects and then are not easy to be represented using fixed values. As the existing OD matrix is usually used as the reference matrix in reserve capacity or practical capacity model and its pattern significantly influences the result of the models, we first consider it as an uncertain variable in this study. And thus the network capacity model is extended to be an optimization with parameter uncertainty.
Researches on other areas of transportation network optimization typically adopted two methods to address the uncertain OD demand [12]: (i) stochastic optimization aims at maximizing the expected profit by assuming that the demand follows a known probability distribution; (ii) robust optimization aims at maximizing the profit with the worstcase scenario of the demand pattern. Considering the exact probability distribution of the OD demand is still hard to be obtained, the robust optimization is more effective in dealing with this problem. If a limited number of discrete scenarios of OD demand patterns are detected, the scenariobased robust optimization [13] is conducted, which is a practical approach usually implemented in transportation projects. It is more general to assume the possibility of the travel demand to be a continuous variable within a bounded set, and the setbased robust optimization can be used for decisionmaking [14]. The uncertainty set is constructed to include most of possible values of the travel demand. The decisionmakers’ attitudes to risk should be considered as well when deciding the shapes and size of the uncertainty sets. It is important to make a tradeoff between the system performance and the level of robustness achieved [13].
In this study, we propose a robust optimization model for the network capacity problem by using the existing OD travel demands as uncertain parameters. The existing demand between each OD pair is assumed to be variable between its upper and lower limits. Besides, three typical uncertainty regions are introduced to provide a bounded set for the uncertain demand. A heuristic solution is developed for the solution to the robust network capacity model. In the next section, the concept of network spare capacity is revisited based on the reserve capacity model. Then, the robust model for network capacity estimation is presented, and the three typical uncertainty sets of existing travel demand are defined. After that, the solution algorithm is described. Computational experiments show the validation and justification of the robust model. Conclusions and perspectives for further research are provided in the last section.
2. Network Spare Capacity and Its Flexibility
The reserve capacity was proposed as the largest multiplier applied to a given existing OD demand matrix that can be allocated to a transportation network without violating any individual link capacity [3]. The product of the largest multiplier and the existing OD demand (represented by vector q) gives the maximum travel demand which can be loaded to the network. For clarity sake we refer to the maximum travel demand as the value of network capacity in rest of this paper. For passenger network, it is well known that multiple OD pairs exist and demands between different OD pairs are not exchangeable or substitutable. Thus, the travel demand pattern or matrix reflects both its quantity and spatial distribution. The method of reserve capacity assumes that the existing OD demand is scaled with a uniform OD growth. The largest value of indicates whether the current network has spare capacity or not. So the network spare capacity is generally explained as follows: if > 1, then the network can be loaded more travel demand and the additional demand can be accommodated by the network which is ; otherwise, that is, < 1, the network is overloaded and the existing OD demand should decrease by to satisfy the capacity constraints [10]. In some researches, the demand multiplier is regarded as the uncertainties in the future OD demand [9, 15].
The classical model of reserve network capacity (RNC) is defined as follows:where is obtained by solving the following user equilibrium problem:where is the OD demand multiplier to all OD demands; is the set of all routes in the network; is the origin index, , and is the set of all origin nodes; is the destination index, , and is the set of all destination nodes; is the capacity of link ; is the flow on link , ; is the vector of all link flows; is the existing trip demand between OD pair ; is the vector of all OD demand; is the flow on route , , between OD pair associated with ; is the vector of flows of all route in ; is the linkroute incidence indicator: 1 if link is on route between OD pair and 0 otherwise; is the travel cost function for link .
In the above model, the upperlevel model maximizes the OD matrix multiplier without violating the capacity constraints (2) for every individual link. The parameter gives the prescribed OD travel demand in the network, which can be obtained according to the current trip demand or a predicted demand pattern accordingly. Route choice behavior and congestion effect are considered in the user equilibrium (UE) model as the lowerlevel model in (3)–(6). Generally, other traffic assignment methods, such as stochastic user equilibrium (SUE) model, can be used in place of the above deterministic UE model as required [10].
The result of the reserve capacity model which is considered may underestimate the capacity of the passenger network, because only the existing OD demand pattern that is more congruous with the network topology would achieve a higher value of network capacity [16]. Basically, the reserve capacity depends on the initial OD demand patterns and route choice behavior of the users. Given the lowerlevel traffic assignment method, the existing OD demand should be the only determinant to the result of the above model. It means that if the given OD matrix is not consistent with the network, the reserve capacity model will produce a result having a low level of maximum demand. Otherwise, if the OD pattern is determined according to the network spatial structure, the travel demand can grow to a very high amount.
Directly applying the result of the reserve capacity may have the following problems. (i) It is hard to decide an exact existing (or predetermined) OD matrix, because the real travel demand pattern is changing at different hours every day and different days every week. Also, it is still very difficult to obtain the full data of the OD demands covering many different hours. (ii) In realworld applications, decisionmakers tend to be risk averse and may be more concerned with the worst cases. Using only a few situations of the OD demand pattern may not provide a robust answer to the network capacity estimation. Conversely, as long as the system performance reaches an acceptable level, it does not matter how much it changes above that level. Thus, it may be more desirable to have an optimization result that performs better in the worst case.
When estimating the capacity of transportation systems, decisionmakers are not only concerned with the extreme results that the total trips can be allocated to a transportation network but also need to evaluate the unknown situations resulted from the fluctuation of the travel demand. Thus, to measure the ability of transportation networks that can deal with the variation of travel demand, Chen and Kasikitwiwat [16] discussed the concept of the network capacity flexibility using three typical network capacity models. The network capacity flexibility is defined as the ability of a transport system to accommodate changes in traffic demand while maintaining a satisfactory level of performance [16, 17]. In this study, integrated with the uncertainties from the existing demand in transportation networks, the network capacity flexibility is further illustrated in Figure 1. On the basis of this, the robust estimation of network capacity is defined as the maximum travel demand can be allocated to a transportation network when satisfying all the possibilities of the uncertain changes in the quantitative and spatial demand pattern. The robust value of the network capacity is also illustrated in Figure 1.
In this study, we extended the reserve capacity model by considering the existing OD demands as uncertain parameters within a certain bounded region. Robust solutions to the network capacity can be conducted using the robust optimization. We utilize the classical reserve capacity model to conduct the robust network capacity for two reasons: (i) the reserve capacity is easy to solve, and the OD travel demand is allowed either increasing or decreasing by applying an OD matrix multiplier greater than one or less than one; (ii) as the existing OD matrix is extended to be an uncertain parameter in the reserve capacity model, the OD distribution is no longer fixed but a variable pattern within some range given by the uncertainty set.
3. Robust Network Capacity Estimation under Demand Uncertainty
In this section, we assume that the prescribed OD trip demand is unknown but bounded within an uncertainty set . Mathematically, the uncertainty set should be closed and convex. In practice, the set of uncertain demand should be derived based on the transportation planners’ knowledge on the uncertainty associated with both the current and future OD travel demand. Nevertheless, it is very difficult to obtain the exact probability distribution of the trip demand between all OD pairs. If the random demand follows a continuous distribution defined from zero to infinity, for example, normal distribution [12], this would require an infinity capacity to meet all possible demand realizations. Thus, the uncertainty set of travel demand is defined as a bounded region, typically utilizing the highest travel demand, , and the lowest demand, , for each OD pair. Using the bounded uncertainty region, the events with low probabilities can be excluded, and then the robust optimization would not provide overly conservative results.
In this study, three typical uncertainty sets were constructed for the existing travel demands.
(1) Interval Constraint [13]. The travel demand between each OD pair which is assumed varies independently within a given interval of . The interval could be the confidence interval of an estimated demand obtained from a survey or by using an OD estimation model. Without additional restraints, the whole uncertainty set will be a box centered at the average travel demand. In this case, it is simple to set for all OD pairs and solve the resulting reserve capacity problem. The capacity value would be estimated in the worst case, which is too conservative. In reality, it is never possible that the demand for all OD pairs reaches their estimated upper bound at the same time. The travel demand pattern is always fluctuant among the OD pairs in different directions. Therefore, it is more reasonable to set an upper limit to the summation of the existing travel demand in network, that is, . The upper limits could be estimated by using the maximum trip volume investigated for the existing travel demand.
(2) Ellipsoid [18]. An ellipsoidal set is generally defined as follows: where , the average OD demand; is a parameter that reflects decisionmakers’ attitudes to risk; and larger indicates that it is more adverse to risk. The value of is from zero to , where denotes the number of OD pairs. When = 1, the uncertainty region is the largest ellipsoid contained in the box region .
(3) Polyhedron. The polyhedron is a set of a finite number of linear equalities and inequalities that restrains the travel demand. It is a generalized form of the box uncertainty set. For example, Sun et al. [19] constructed the following polyhedron region for uncertain OD demand: where is the nominal value of the travel demand for OD pair . It may choose the mean value of the interval constraints or an observed result from a survey. This uncertainty set allows the demand pattern varying entirely around its nominal value and involves the implicit possible interactions among OD demands. Therefore, the overly conservative results may be avoided. The last two sets of constraints require that the uncertain travel demands meet the conservation condition with the nominal demand matrix at the zonal production and attraction.
Note that the shape of uncertainty set affects the efficiency and robustness of network capacity value. BenTal and Nemirovski [14] suggested applying the min–max optimization model. Once the uncertainty set of the travel demand, , is determined, the min–max model will find a robust solution that tolerates changes in travel demand up to the given bound. Using any type of the uncertainty sets, , the robust reserve network capacity (RRNC) problem can be formulated as follows:where solves
The above model is referred to as the robust counterpart of the original reserve network capacity problem. The solution of the robust counterpart results in a maximum total travel demand scheme under the corresponding worstcase demand pattern.
4. Solution Algorithm
A heuristic algorithm is proposed to solve the above robust optimization model. It takes a similar framework as the procedure presented in [18], which is referred to as the cutting plane algorithm to robust optimization. The algorithm involves an iterative procedure to solve two inner optimization problems alternately until the convergence criterion is satisfied. The algorithm is presented as follows.
Step 0 (initialization). Give the initial values of the OD demand (usually , ) and solve the following reserve network capacity (RNC) problem to produce an initial demand multiplier :Set the iteration counter .
Step 1 (direction finding).
Step 1.1. Solve the following inner (worstcase scenario (WCS)) problem with the determined to obtain the worstcase demand scenario:Step 1.2. Formulate a RNC problem with the scheme of existing demand produced from the WCS problem in Step 1.1. Solve this inner problem to find a search direction of the maximum demand multiplier .
Step 2 (move). Compute , where is the step length. In this study, the step length is chosen as , which is used in the method of successive averages.
Step 3 (convergence check). If the objective value of the WCS problem or reaches the maximum iterations, then stop, where there is a predetermined convergence criterion. Otherwise, denote the solution of the WCS problem as , and go to Step 1; set .
Remark 1. In the above steps, the WCS problem is formulated to find a solution of , in which case the traffic flows on all links which exceed their capacities the most. If is an optimal solution to the RRNC, the corresponding optimal objective value of the WCS problem must be zero. Otherwise, an improved solution may be obtained by solving the RNC problems in Step 1.2 which is a relaxation of the RRNC problems with a specific demand pattern. In the process of the algorithm, each of these relaxed RRNC problems can approximate the original RRNC better than its predecessors. Although it is still difficult in practice to find a global optimum of the relaxed RRNC and WSC, Yin et al. showed that the cutting plane algorithm is effective in providing a good solution to the robust optimization problem [18]. The relaxed RRNC problems and the WSC problems are solved by the sensitivity analysis based (SAB) algorithm [20].
Remark 2. The second inner problem is a standard RNC model when the existing OD demand is determined. The RNC can be solved efficiently by applying the SAB algorithm [3]. The SAB algorithm locally approximates the original bilevel problem as a singlelevel optimization by using firstorder Taylor expansion. The derivatives of lowerlevel decision variables with respect to upperlevel ones are utilized for the linear approximation. The derivatives can be conducted from the sensitivity analysis of the lowerlevel model.
In this study, we used the restriction approach for the sensitivity analysis of the lowerlevel UE model. The restriction approach was proposed by Tobin and Friesz [21] and then corrected by Yang and Bell [22] for its flaws on selecting the nondegenerate extreme point. One can also refer to Du et al. [23] for the details of this approach. In this section, some necessary results are present without proof.
For the reserve capacity model, the link flows in upperlevel, , are represented as an implicit function of the OD matrix multiplier as constraint (2) shows. Using the firstorder Taylor expansion, it can be approximated aswhere is the given solution of the OD demand multiplier at the current iteration of SAB algorithm.
From the results in Tobin and Friesz [21], the derivatives of the route flows, , with respect to the OD demand multiplier , are derived as follows:where the superscript “” denotes that the variables or matrices are only associated with the restricting subproblem derived by the restriction approach (applying the correction in Yang and Bell [22]) and the superscript “” represents the transposed matrix. Other notations are defined as follows: is the linkroute incidence matrix; is the ODroute incidence matrix, where equals 1 if OD pair is connected by route , and 0 otherwise; is the Lagrangian multiplier associated with constraint (4); is the vector of the travel cost function of all links with the equilibrium link flow for the perturbation parameters at 0.Thus, the derivatives of the link flows to the multiplier are obtained by . Based on the above derivations, the SAB algorithm can be used for solving the RNC problem.
Remark 3. The WCS problem is also formulated as a bilevel programming using equilibrium constraints, so the SAB method can also be modified for its solution. The implicit relationship is firstorder approximated aswhere is the solution of the existing OD demand at the current iteration of SAB algorithm.
The derivatives of the route flows, , with respect to the existing OD demand , are derived from the following equation:Because in this inner problem the value of is fixed, the derivative will be a diagonal matrix with the value of on its diagonal. Then, the derivatives of to are calculated by . Thus, at each iteration of the SAB method, the approximating WCS problem isNote that the above localized approximation problem is a nonlinear problem and a number of optimization tools in commercial software packages could be used for its solution. In this study the approximate WCS problem is solved using MATLAB builtin functions which were converted to a .NET component and used in our solution program in C# language. Similar to [18], we randomly generate 100 vectors of from the uncertainty region . For each , the user equilibrium problem with is solved. Let be the random with the maximum objective value of the WCS problem, and denote the objective value as . If , then set . Otherwise, use the MATLAB functions to solve the approximate WCS problem with as an initial solution. If it gives a solution with an objective larger than , the solution is used as for the next iteration; otherwise, the algorithm is terminated with an optimal solution .
5. Computational Experiments
5.1. Experiment 1: NguyenDupuis Network
Computational experiments are presented in this section to illustrate the results of the robust network capacity model. The example is based on a road network which is adopted from Nguyen and Dupuis [24] as Figure 2 shows. It consists of 13 nodes, 19 links, and 4 OD pairs. The nominal value of the existing travel demand is given by the OD matrix in Figure 2 (denoted by ). The characteristics of the links are listed in Table 1. The Bureau of Public Road link performance function was used in the experiments:where is the freeflow travel time for link .

We applied the proposed approach for robust network spare capacity estimation with the three typical uncertainty sets which are described in this paper. Assume that the intervals for the OD travel demands are , , , and for OD pairs 12, 13, 42, and 43, respectively.
(1) Interval Region. Theoretically, under the interval constrains for each OD pair, the total of the existing OD demand can vary from 1400 to 2975, which covers a wide range. However, according to our practical experience, the travel demands between the OD pairs may not reach their maximum simultaneously. Thus, an additional constraint was introduced to give an upper bound of the total existing demand. To be consistent with the interval constraints, was set to change from 1400 to 3000. Figure 3 illustrates how uncertainty of the existing travel demand affects the reserve capacity results of the example network. The reserve capacity value is calculated by , where is the robust solution of the RRNC problem and is the corresponding travel demand pattern. Note that the reserve capacity value and the decrease synchronously along with the growth of . Because when is increased, the vector of travel demand pattern, , will have more space to change its spatial distribution, and this makes the existing demand easier to archive a demand pattern whose corresponding multiplier can reach a smaller value. Besides, when exceeds 2700, the value of cannot get worse. At this value , the corresponding existing travel demand is , which produces a smallest reserve capacity value 363.42 (referred to as the conservative solution in the rest part of this section). Furthermore, for the robust solutions when > 2700, the change of will have no effect on the value of . Here, we use the lower bound of to conduct the worstcase results of the reserve capacity. According to the properties of the robust optimization, the solutions of the existing OD demand may not be unique, so no specific solution of the existing demand is presented in this paper.
(2) Ellipsoidal Region. For the ellipsoidal uncertainty set, the parameter is set to be 0, 0.2, 0.5, 1.0, and 2.0. The center of the ellipsoid is decided by the boundary for each OD demand (not the nominal value for this example). The computational results in Table 2 show that the robust solution to changes between 0.1795 and 0.1286 with a descending trend, as well as the corresponding network capacity value.

(3) Polyhedral Region. For the robust network capacity estimation with the polyhedral set described in previous section, the robust results at = 0, 0.25, 0.5, 0.75, and 1.0 are computed and presented in Table 2. The nominal demand (given in Figure 2) is used in this type of uncertainty. The value of declines from 0.1944 to 0.1753. The same trend has been observed on the robust values of network capacity, which varies from 388.89 to 350.54.
Table 3 reports the solutions of the reserve capacity value from two RNC problems with and separately and two RRNC problems with defined in the ellipsoid region () and the polyhedral region (). We refer to the first two solutions as “nominal” and “conservative” estimation, respectively, and the latter two as “robuste” and “robustp” separately. For the results, the robust solution for the OD demand multiplier indicates that all the values of the existing demand within the uncertainty set can be applied by a multiplier larger than this robust value without violating the capacity constraints. A corresponding maximum travel demand is produced as the robust network capacity solution for the same uncertainty set. Since the same maximum travel demand may produce by different combinations of the multiplier and feasible existing demand, the proposed algorithm only finds the smallest value of the maximum multiplier. Besides, if the shape of the uncertainty set changes, different solutions for the robust multiplier and maximum demand will be produced. Table 3 shows these features.

To evaluate the four estimations of the network capacity, we randomly generate 500 samples of OD demand as the possible realizations of the existing travel demand pattern. The samples are from the normal distribution, , within its interval for all OD pairs. Then, for each sample, one has the following.
Firstly, the user equilibrium assignment associated with each reserve network capacity estimation (i.e., the largest ) is computed. The number of failures is counted whenever the link flow exceeds its capacity at each sample. In the end of the simulation, the successful rates are computed as [1 − total number of failures/(number of samples × number of links)]. The results are referred to as “percentage of meeting capacity constraints” in Table 3. Obviously, the conservative estimation gives the worst value of the network capacity; the nominal estimation produces a medium level of result. By comparison, the robust results from the ellipsoid and polyhedron are more reasonable. Note that all samples generated in the above simulations belong to a box region comprised of the intervals for OD demands. The aforementioned uncertainty sets are used as replacements of the box regions, so as to prevent overly conservative robust resolutions. The results in Table 3 indicate that the ellipsoid approximate the box region better than the polyhedron, because more random samples will not exceed the network capacity after they are scaled by the robust multiplier = 0.1498. On the other hand, although = 0.1843 is the worst case in the polyhedral set, when the boundary is extended to the box region, only 72.4% of the random samples can be covered.
Furthermore, considering that the demand multiplier in essence is a relative value, the reserve capacity results are derived. Therefore, the reserve capacity problem with the every random existing OD demand is solved in our test. For each estimation value of the reserve capacity: the number of successes is counted whenever the reserve capacity value of the sample exceeds the robust capacity estimation. Consequently, the successful rates are computed as [total number of successes/number of samples]. These results are also presented in Table 3 as “percentage of aboverobustestimation.” From this aspect, the polyhedron provides a more robust estimation of the network capacity, which can be met by most realizations of the existing travel demand pattern (98.8% versus 54.6% compared with the ellipsoid). We may suppose that the worst cases of network capacity values exist in the corner area of the boxshaped uncertain region, where it is not easy to be covered by the corresponding ellipsoidal set. Consequently, the choice of the robust results depends on the decisionmakers’ attitude to risk and the usage of the robust results. From the computational perspective, the robust optimizations with polyhedral uncertainty sets have more advantages. More effective solution methods could be developed in future studies.
5.2. Experiment 2: SiouxFalls Network
Experiments are further presented on the SiouxFalls network [25]. The network contains 24 nodes, 76 links, and 528 OD pairs. The characteristics of the links and travel demands are also provided in BarGera [25]. In this experiment, we used the default values of the travel demand in as the nominal value of the existing demand in the network. If ≥ 1000, then the demand for OD pair is uncertain and its upper limit and lower limit are set to and , separately. Hence, there are 104 OD pairs with uncertain existing demands. According to the previous section, four estimation methods are employed to evaluate the network capacity of the SiouxFalls.
Table 4 reports the solutions of the reserve capacity from the four estimation methods. Note that the robuste solution gives a moderate robust result compared to the others. The robustp solution provides the lowest estimation of network capacity. One may use this lowest value as the worstcase performance that the network can serve the travel demand. Besides, note that the lowest estimation of the network capacity was not derived from the conservative solutions, in which the uncertain OD demands are set to its upper limits. Although the conservative solution is corresponding to the lowest demand multiplier, it may not reflect the most unfavorable situation which is possible to be resulted from the changes of the network demands.

We selected the conservative, robuste, and robustp solutions to further inspect the link flow patterns at the maximum travel demand situations (i.e., the reserve capacity). The link flow patterns are shown in Figures 4(a), 4(b), and 4(c). The width of the line indicates the traffic volume through the link. The red lines show the links whose V/C (volume/capacity) ratio is greater than 0.9, and the black lines denote the links with zero flow. As a reference, we also randomly generate 500 samples of OD demand which are from the normal distribution, , within its interval . For each sample, the reserve capacity problem is solved. Then, at each link, the highly saturated number is counted whenever its V/C ratio is greater than 0.9. In consequence, the link saturation rates are computed as [total the saturated number/number of samples]. The results are also shown in Figure 4(d). The width of the line indicates the saturation rate of the link. The red lines indicate links with saturation rates greater than 50%. The orange lines denote links with saturation rates lower than 50% but greater than zero. The black lines mean the links have no chance to be blocked. From Figure 4, we conclude that the solution of the robuste will produce a link flow pattern most likely to be realized when the network reaches its capacity (compared to the simulation results). The robustp solution provides an extreme result that the network is congested only because of a very few links. These links could be considered as the most critical links which restrict the capacity of the entire road network. By contrast, the conservative solution seems not to fit the simulation results very well. In conclusion, the robust estimation to the network capacity problem could be much more practical than the results from any specific scenario.
(a) Conservative solution
(b) Robust solution with ellipsoid ()
(c) Robust solution with polyhedron ()
(d) Simulation results
6. Conclusions and Perspectives
In this study, a robust network capacity model with uncertain demand has been proposed. The robust optimization is formulated using the min–max model with a bounded uncertainty set of the existing OD travel demands. With the uncertainty set, the lowprobability realizations of the travel demand pattern are excluded, and thus the robust model can produce a proper estimation of network capacity which can be achieved with a large probability. Then, a heuristic algorithm has been proposed for the proposed robust model. It solves two inner problems iteratively: one is the worstcase scenario problem; and the other is the relaxed robust optimization, namely, the standard reserve capacity model. At each iteration, the cutting plane method has been adopted to generate the worstcase demand scenario, and the sensitivity analysis based approach has been developed for the solution of the worstcase model and the reserve capacity model. The validity and performance of the proposed robust model have been demonstrated in the computational experiments. Different results under three typical uncertainty sets, say interval, ellipsoidal, and polyhedral region, have been conducted and compared. The interval set is simple but easy to produce too conservative results; the ellipsoidal set is a good approximation to the uncertain region and produces results with moderate robustness, but its solution is more complicated due to the nonlinear constraints; the polyhedral set is considered if a high level of robustness is required, and its linear formulation makes the robust model easier to solve. Furthermore, by conducting computational experiments on the SiouxFalls network, robust solutions shown can provide more practical results of the link flow patterns. In applications, these uncertainty sets and their parameters should be selected according to the desired level of robustness.
The robust model based on the reserve capacity model has been proposed and explored in this study. Future researches should focus on more efficient solution approaches for the robust problem with the min–max model. The experiments on largescale networks are also needed. Besides, the demand uncertainties existing in other network capacity models are also expected to be detected and discussed. Alterative traffic assignment model, such as the stochastic user equilibrium, could also be discussed for the network capacity problems. The robust solution of network capacity gives a lower bound to the possible schemes of the maximum demand in a transportation network. These possibilities constitute a range where the robust solution can be most likely to be reached in reality. Therefore, the robust solution to network capacity problems needs to receive more attentions in transportation planning applications.
Disclosure
An initial version of this paper was presented at the Transportation Research Board (TRB) 96th Annual Meeting. It therefore also appears in the proceedings of the TRB 96th Annual Meeting Compendium of Papers.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This research is supported by “the National Natural Science Foundation of China (no. 51508161),” “the Natural Science Foundation of Jiangsu Province (no. BK20150817),” and “the Fundamental Research Funds for the Central Universities (no. 2017B12414).”
References
 R. K. Ahuja, T. L. Magnanti, and J. Orlin, Network Flows: Theory, Algorithms, and Applications, PrenticeHall, New Jersey, NJ, USA, 1993. View at: MathSciNet
 Y. Asakura, “Maximum capacity of road network constrained by user equilibrium conditions,” in Proceedings of the 24th Annual Conference of the UTSG, 1992. View at: 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
 P. Kasikitwiwat and A. Chen, “Analysis of transportation network capacity related to different system capacity concepts,” Journal of the Eastern Asia Society of Transportation Studies, vol. 6, pp. 1439–1454, 2005. View at: Google Scholar
 G. Ziyou and S. Yifan, “A reserve capacity model of optimal signal control with userequilibrium route choice,” Transportation Research B: Methodological, vol. 36, no. 4, pp. 313–323, 2002. 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, vol. 8, no. 2, pp. 87–99, 2004. 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  MathSciNet
 A. Chen, P. Kasikitwiwat, and C. Yang, “Alternate capacity reliability measures for transportation networks,” Journal of Advanced Transportation, vol. 47, no. 1, pp. 79–104, 2013. View at: Publisher Site  Google Scholar
 X. Xu, A. Chen, S. Jansuwan, K. Heaslip, and C. Yang, “Modeling transportation network redundancy,” Transportation Research Procedia, vol. 9, pp. 283–302, 2015. View at: Publisher Site  Google Scholar
 M. Du, L. Cheng, X. Jiang, and Z. Li, “Sensitivity based heuristics for network capacity estimation in transportation,” in Proceedings of the 94th Transportation Research Board Annual Meeting, Washington, D.C., USA, 2015. View at: Google Scholar
 K. An and H. K. Lo, “Twophase stochastic program for transit network design under demand uncertainty,” Transportation Research Part B: Methodological, vol. 84, pp. 157–181, 2016. View at: Publisher Site  Google Scholar
 Y. Yin, S. M. Madanat, and X.Y. Lu, “Robust improvement schemes for road networks under demand uncertainty,” European Journal of Operational Research, vol. 198, no. 2, pp. 470–479, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 A. BenTal and A. Nemirovski, “Robust optimization—methodology and applications,” Mathematical Programming, vol. 92, no. 3, pp. 453–480, 2002. View at: Publisher Site  Google Scholar  MathSciNet
 S.W. Chiou, “Optimization of robust area traffic control with equilibrium flow under demand uncertainty,” Computers and Operations Research, vol. 41, no. 1, pp. 399–411, 2014. View at: Publisher Site  Google Scholar
 A. Chen and P. Kasikitwiwat, “Modeling capacity flexibility of transportation networks,” Transportation Research A: Policy and Practice, vol. 45, no. 2, pp. 105–117, 2011. View at: Publisher Site  Google Scholar
 E. K. Morlok and D. J. Chang, “Measuring capacity flexibility of a transportation system,” Transportation Research Part A: Policy and Practice, vol. 38, no. 6, pp. 405–420, 2004. View at: Publisher Site  Google Scholar
 Y. Yin, S. Lawphongpanich, and Y. Lou, “Estimating investment requirement for maintaining and improving highway systems,” Transportation Research Part C: Emerging Technologies, vol. 16, no. 2, pp. 199–211, 2008. View at: Publisher Site  Google Scholar
 H. Sun, Z. Gao, and J. Long, “The robust model of continuous transportation network design problem with demand uncertainty,” Journal of Transportation Systems Engineering and Information Technology, vol. 11, no. 2, pp. 70–76, 2011. View at: Publisher Site  Google Scholar
 T. L. Friesz, R. L. Tobin, H.J. Cho, and N. J. Mehta, “Sensitivity analysis based heuristic algorithms for mathematical programs with variational inequality constraints,” Mathematical Programming, vol. 48, no. 2, pp. 265–284, 1990. View at: Publisher Site  Google Scholar  MathSciNet
 R. L. Tobin and T. L. Friesz, “Sensitivity analysis for equilibrium network flow,” Transportation Science, vol. 22, no. 4, pp. 242–250, 1988. View at: Publisher Site  Google Scholar  MathSciNet
 H. Yang and M. G. H. Bell, “Sensitivity analysis of network traffic equilibria revisited: the corrected approach,” in Selected Proceedings of the 4th IMA International Conference on Mathematics in Transport, pp. 373–395, 2007. View at: Google Scholar
 M. Du, L. Cheng, and H. Rakha, “Sensitivity analysis of combined distributionassignment model with applications,” Transportation Research Record, no. 2284, pp. 10–20, 2012. View at: Publisher Site  Google Scholar
 S. Nguyen and C. Dupuis, “An efficient method for computing traffic equilibria in networks with asymmetric transportation costs,” Transportation Science, vol. 18, no. 2, pp. 185–202, 1984. View at: Publisher Site  Google Scholar
 H. BarGera, Transportation Network Test Problems, 2001, 2017, http://www.bgu.ac.il/~bargera/tntp.
Copyright
Copyright © 2017 Muqing Du 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.