Scenario reduction in stochastic programming: An approach using probability metrics

Given a convex stochastic programming problem with a discrete initial probability distribution, the problem of optimal scenario reduction is stated as follows: Determine a scenario subset of prescribed cardinality and a probability measure based on this set that is the closest to the initial distribution in terms of a natural (or canonical) probability metric. Arguments from stability analysis indicate that Fortet-Mourier type probability metrics may serve as such canonical metrics. Efficient algorithms are developed that determine optimal reduced measures approximately. Numerical experience is reported for reductions of electrical load scenario trees for power management under uncertainty. For instance, it turns out that after 50% reduction of the scenario tree the optimal reduced tree still has about 90% relative accuracy.

[1]  J. Dupacová Stability and sensitivity-analysis for stochastic programming , 1991 .

[2]  Werner Römisch,et al.  Stability analysis for stochastic programs , 1991, Ann. Oper. Res..

[3]  Jim Freeman Probability Metrics and the Stability of Stochastic Models , 1991 .

[4]  Roger J.-B. Wets,et al.  Quantitative stability of variational systems: III.ε-approximate solutions , 1993, Math. Program..

[5]  S. Zenios,et al.  Constructing Optimal Samples from a Binomial Lattice , 1993 .

[6]  Alexander Shapiro,et al.  Quantitative stability in stochastic programming , 1994, Math. Program..

[7]  Jitka Dupačová Scenario-Based Stochastic Programs: Strategies for Deleting Scenarios , 1995 .

[8]  William T. Ziemba,et al.  Concepts, Technical Issues, and Uses of the Russell-Yasuda Kasai Financial Planning Model , 1998, Oper. Res..

[9]  Werner Römisch,et al.  Qantitative stability for scenario-based stochastic programs , 1998 .

[10]  S. Rachev,et al.  Mass transportation problems , 1998 .

[11]  Michael A. H. Dempster,et al.  EVPI‐based importance sampling solution proceduresfor multistage stochastic linear programmeson parallel MIMD architectures , 1999, Ann. Oper. Res..

[12]  Jitka Dupacová,et al.  Portfolio optimization via stochastic programming: Methods of output analysis , 1999, Math. Methods Oper. Res..

[13]  Stavros A. Zenios,et al.  Scenario modeling for the management ofinternational bond portfolios , 1999, Ann. Oper. Res..

[14]  Werner Römisch,et al.  Stochastic Lagrangian Relaxation Applied to Power Scheduling in a Hydro-Thermal System under Uncertainty , 2000, Ann. Oper. Res..

[15]  Jitka Dupacová,et al.  Sensitivity of Bond Portfolio's Behavior with Respect to Random Movements in Yield Curve: A Simulation Study , 2000, Ann. Oper. Res..

[16]  Jitka Dupacová,et al.  Scenarios for Multistage Stochastic Programs , 2000, Ann. Oper. Res..

[17]  Rüdiger Schultz,et al.  Some Aspects of Stability in Stochastic Programming , 2000, Ann. Oper. Res..

[18]  Georg Ch. Pflug,et al.  Scenario tree generation for multiperiod financial optimization by optimal discretization , 2001, Math. Program..

[19]  K. Kiwiel,et al.  Power management in a hydro-thermal system under uncertainty by Lagrangian relaxation , 2002 .

[20]  Svetlozar T. Rachev,et al.  Quantitative Stability in Stochastic Programming: The Method of Probability Metrics , 2002, Math. Oper. Res..