Article ID: | iaor20171670 |
Volume: | 51 |
Issue: | 2 |
Start Page Number: | 518 |
End Page Number: | 535 |
Publication Date: | May 2017 |
Journal: | Transportation Science |
Authors: | Homem-de-Mello Tito, Zhang Leilei |
Keywords: | combinatorial optimization, stochastic processes, risk, simulation |
In this paper we discuss stochastic optimal path problems where the goal is to find a path that has minimal expected cost and at the same time is less risky (in terms of travel time) than a given benchmark path. The model is suitable for a risk‐averse traveler who prefers a path with a more guaranteed travel time to another path that could be faster but could also be slower. Such risk attitude is incorporated using the concept of second order stochastic dominance constraints. A recently developed theory for optimization problems with stochastic dominance constraints ensures that the resulting problem can be written as a large linear integer program with binary variables; for networks of realistic size, however, such a direct approach is not practical because of the size of the resulting optimization problem. Moreover, the solution returned by the model may contain cycles, which is clearly undesirable from a practical perspective. A number of strategies are explored to solve the problem. First, we prove that cycles can be prevented by a simple modification of the model if the arc travel times are mutually independent. We then propose a sample average approximation approach to the problem using samples from the distribution of travel times. Because of the randomness resulting from sampling, it is important that statistical guarantees for the solution returned by the algorithm be given, and we provide heuristic procedures to deal with stochastic constraints. We also incorporate a branch‐and‐cut approach that exploits the structure of the problem to deal with the integrality constraints more efficiently. We present some numerical experiments for a 1,522‐arc system that corresponds to a large portion of the Chicago area network. The results show that our approach can solve the problem very effectively, producing solutions with statistical guarantees of optimality within a reasonable computational time.