Stochastic scenario decomposition for multistage stochastic programs

Abstract

Over the past decade, several stochastic approaches have been proposed for two-stage stochastic programs. Many of these algorithms have attractive computational as well as conceptual properties (e.g. convergence with probability one). This paper expands the realm of such approaches to multistage convex stochastic programming problems. We present a stochastic scenario decomposition (SSD) algorithm which is a statistically motivated cutting plane algorithm for the solution of multistage stochastic programs. The method is based on solving a dual problem in which the variables correspond to multipliers associated with the non-anticipativity constraints of the primal problem. Our analytical results verify conditions under which SSD identifies an optimal solution asymptotically. We also overcome some computational hurdles resulting from increases in the column dimension of the SSD master problem. We propose a variable aggregation scheme that allows us to solve much smaller master programs without sacrificing solution quality. Our computational results demonstrate the effectiveness of this aggregation scheme in solving the SSD master program.