Title: A DecisionMaking Framework for More Sustainable Urban Environments
1A Decision-Making Framework for More Sustainable
Urban Environments
Victoria Chen and Zehua Yang Industrial
Manufacturing Systems Engineering University of
Texas at Arlington (UTA)vchen_at_uta.edu
Julia Tsai Purdue University
Michael Chang, Ellis Johnson, Eva Lee, and Terry
Murphy Georgia Institute of Technology
Supported by TSE/EPA Contract R-82820701-0
2Outline
- Decision-Making Framework (DMF)
- Stochastic Dynamic Programming (SDP)
- Wastewater Treatment Application
- Metropolitan Atlanta Ozone Pollution Application
3A Modular Decision-Making Framework
- For time period/level/stage t
- xt state of system
ut decision/control
4Applications
- Inventory Forecasting
- Determine optimal order quantities using
forecasts for product demand. - Airline Yield Management (Sabre)
- Develop an on-line YM policy for a real airline
network. - Water Resources (Italy)
- Determine optimal control of a water reservoir
network. - Wastewater Treatment System (EPA)
- Evaluate current and future technologies
according to economic and sustainability
measures. - Ozone Pollution (EPA)
- Identify key times and locations for emissions
reductions to avoid ground-level ozone
exceedances.
5Stochastic Dynamic Programming (SDP)
- Objective Minimize expected cost over T stages.
- Optimal Value Function Ft(xt) in Stage t Minimum
expected cost to operate the system over stages t
through T.
6SDP Optimization Formulation
Objective Transition Constraints
For stage t xt is the vector of state
variables ut is the vector of decision
variables ?t is the vector of random
variables ct(.) is the cost function ft(.) is
the multivariate transition function ?t is the
set of constraints
7SDP Optimal Value Function
- For stages t through T
- Recursive formulation
- SDP solves backwards (t T,, 1) to determine
the optimal value functions Ft(xt) in all stages.
8SDP Statistical Modeling Process
9Wastewater Treatment System
- Research with Jining Chen (Tsinghua U, China),
Bruce Beck (U of Georgia), and Julia
Tsai (Purdue U) - Decision-Making Framework (DMF)
- Evaluate current and emerging technologies,
including usage at a unconventional level. - Provide most promising directions for future
research. - Liquid treatment line up to 11 levels.
- Solid treatment line up to 6 levels.
- Advantages over Typical Approach
- Comprehensive vs. isolated evaluation of
processes. - Exploration of many potential levels of
treatment.
10(No Transcript)
11(No Transcript)
12Wastewater Treatment DMF
- Objective Minimize capital and operating costs
- size (land area, volume) odor emissions.
- Maximize robustness global desirability.
- Complication Dependencies between some
technologies. - Stages Levels of the treatment system (T 17)
- State xt entering Level t Liquid and solid
attributes. - Decision ut in Level t Technology process
employed. - Constraints Attribute quantities exiting each
Level. - Random Variables Performance and cost of the
technologies. - Transition xt1 xt ? quantity removed
quantity added.
13Wastewater Treatment DMF
- Monitored attributes are the 20 SDP state
variables
14Dependencies Between the Technologies
Diagnostic state variables indicate what
technology has been utilized in an earlier level.
15DMF Solution (Economic Cost)
16DMF Solution (Economic Cost)
17DMF Solution (Economic Cost)
18DMF Solution (Economic Cost)
19DMF Solution (Economic Cost)
20Controlling Ozone Pollution
- Research with Michael Chang (Georgia Tech),
Melanie Sattler (UTA), and John Priest
(UTA) - Decision-Making Framework (DMF)
- Explore the necessary emissions reductions over
time and space to prevent an ozone exceedance. - Identify reductions in emissions that will lead
to new control strategies. - Provide information and guidance to government
decision makers for creating new control
strategies - Advantages over Typical Approach
- Comprehensive approach vs. trial and error.
- Dynamic and focused vs. static control strategies.
21Ozone Pollution DMF
- Objective Minimize cost to avoid ozone
exceedance days. - Complications
- Ozone is not directly emitted, but is formed by
the reaction of emitted pollutants, such as
nitrogen oxides (NOx). - State variables are temporally and spatially
related. - Stages Hours or groups of hours.
- State xt at the beginning of Stage t Ozone, NOx
at different locations and at hours prior to
stage t. - Decision ut in Stage t Reductions of NOx by
location. - Constraint Ozone exceedance limit ? adds penalty
cost - Random Variables Uncertainty in how Ozone, NOx
change over time and space.
22Ozone Pollution DMF
- Transition Atmospheric Chemistry Module
- E1 NOx emissions between 6am?9am
- E2 NOx emissions between 9am?12pm
23Atlanta Urban Airshed Model (UAM)
- U.S. EPA (1990)
- UAM An advanced three-dimensional, photochemical
air quality grid model that encompasses a
160?160 km square region containing the
metropolitan area. - Spatial modeling grid 40?40
- Point sources 102
- Temporal modeling 24 hours
- July 29 - August 1, 1987 episode
- Computational Issues one UAM simulation run gt
one hour - Impractical to implement within SDP optimization
- Instead develop a metamodel of the UAM.
24Atmospheric Chemistry Metamodel
25Atlanta Modeling Domain
1-15 were Georgia
Power
26Preliminary Metamodeling
- Transition Function Metamodel Setup
- All 5?5 grid regions and 15 significant point
sources - Reduce 24 hours ? five 3-hour time periods
- time period 0 from 4 a.m. to just before 7 a.m.
- time period 1 from 7 a.m. to just before 10
a.m. - time period 2 from 10 a.m. to just before 1
p.m. - time period 3 from 1 p.m. to just before 4 p.m.
- time period 4 from 4 p.m. to just before 7 p.m.
- Total NOx over each time period for each grid
region and point source ? 5(2515) 200
predictor variables - Maximum hourly-averaged ozone at the 4 stations
and in time periods 1, 2, 3, 4 ? 16 response
variables
27Preliminary Metamodeling
- Transition Function Metamodel Fit
- Experimental Design A 500-point Latin hypercube
design determined variations in NOx emissions
from zero up to the nominal value. - Statistical Model Multiple linear regression
models at each of the 4 Atlanta monitoring
stations for maximum ozone in time period p 1,
2, 3, 4, as a function of - maximum ozone at the 4 Atlanta monitoring
stations in time periods earlier than p - total NOx emissions over the 25 grid regions and
at the 15 point sources in time periods p and
earlier - Achieved regression R2 gt 0.90 for all except
sites Tucker and S. Dekalb in time period 2.
28State and Decision Variables
- Initial State Space for Stage t For all time
periods lt t, - maximum ozone at the 4 Atlanta monitoring sites
- total NOx emissions over the 25 grid regions and
at the 15 point sources - ? Stage 1 44 dimensions
- Stage 2 88 dimensions
- Stage 3 132 dimensions
- Stage 4 176 dimensions
- Initial Decision Space for Stage t For time
period t, - total NOx emissions over the 25 grid regions and
at the 15 point sources - ? 40 dimensions
29State and Decision Variables
- Key Only those variables included in the
metamodels need to be maintained in the state and
decision spaces. - Reduced State Space ? Stage 1 17 dimensions
- Stage 2 25 dimensions
- Stage 3 23 dimensions
- Stage 4 19 dimensions
- Reduced Decision Space ? Stage 1 17 dimensions
- Stage 2 9 dimensions
- Stage 3 9 dimensions
- Stage 4 3 dimensions
- SDP Dimensionality The dimension of the SDP
optimization problem is the largest state space
dimension ? 25
30Emissions Reductions Decisions for Atlanta
31Last Period SDP Optimal Value Function
32Metamodel Refinement
- Multivariate Adaptive Regression Splines (MARS)
- Employed new robust MARS algorithm that favors
lower-order terms. - Maximum ozone model for S. Dekalb in time period
2. - ? R2 98.6
- Maximum ozone model for Tucker in time period 2.
- ? R2 91.5
- Resulting metamodels were clearly curvilinear.
- Undesirable chemical reaction in which high NOx
leads to low ozone. - However, SDP transition functions must be
monotonic in order to guarantee convexity for the
optimization.
33Maximum Ozone at S. Dekalb in Period 2
34Restructured S. Dekalb Metamodel
35On-going and Future Work
- DMF and SDP
- Alternate experimental designs / statistical
models - Parallel computing
- Wastewater Treatment
- Verify technology selections via simulation
- Consider other objectives (odor emissions,
robustness) - Ozone Pollution
- Complete SDP solution for Atlanta
- Simulate solution via Atlanta Urban Airshed Model
- Refine transition function metamodels
- Transfer Atlanta DMF to Dallas/Fort-Worth
36References
- Chen, V. C. P. (1999). Application of
Orthogonal Arrays and MARS to Inventory
Forecasting Stochastic Dynamic Programs.
Computational Statistics and Data Analysis, 30,
317-341. - Chen, V. C. P., D. Ruppert, and C. A. Shoemaker
(1999). Applying Experimental Design and
Regression Splines to High-Dimensional
Continuous-State Stochastic Dynamic Programming.
Operations Research, 47, 38-53. - Chen, V. C. P., Chen, J., and Beck, M. B. (2000).
Statistical Learning within a Decision-Making
Framework for More Sustainable Urban
Environments. In Proceedings of the Joint
Research Conference on Statistics in Quality,
Industry, and Technology, Seattle, Washington. - Chen, V. C. P. (2001). Measuring the Goodness
of Orthogonal Array Discretizations for
Stochastic Programming and Stochastic Dynamic
Programming. SIAM Journal of Optimization, 12,
322-344.
37References
- Chen, V. C. P., J. C. C. Tsai, E. K. Lee, and E.
L. Johnson (2001). A Decision-Making Framework
for Evaluating Wastewater Treatment
Technologies. In Proceedings of the 5th Annual
Green Chemistry and Engineering Conference,
Washington, D. C. - Tsai, J. C. C. (2002). Statistical Modeling of
the Value Function in High-Dimensional,
Continuous-State Stochastic Dynamic Programming.
PhD dissertation, Georgia Institute of
Technology. - Chen, V. C. P., K.-L. Tsui, R. R. Barton, and J.
K. Allen (2003). A Review of Design and
Modeling in Computer Experiments. Handbook in
Statistics Statistics in Industry, 22, 231-261,
Khattree, R. and Rao C. R. (eds.), Elsevier
Science, Amsterdam. - Chen, V. C. P., Z. Yang, M. E. Chang, and T. E.
Murphy (2003). A Decision-Making Framework for
Studying Ozone Pollution in Urban Atlanta. In
Proceedings of the 7th Annual Green Chemistry and
Engineering Conference, Washington, D. C.
38References
- Tsai, J. C. C., V. C. P. Chen, B. S. Boardman
(2003). Flexible Implementation of Multivariate
Adaptive Regression Splines in Solving
Continuous-State Stochastic Dynamic Programming.
In Proceedings of the 2003 IE Research
Conference, Portland, OR. - Tsai, J. C. C., V. C. P. Chen, J. Chen, and M. B.
Beck (2004). Stochastic Dynamic Programming
Formulation for a Wastewater Treatment
Decision-Making Framework. Annals of Operations
Research, Special Issue on Applied Optimization
Under Uncertainty, to appear. - Cervellera, C., V. C. P. Chen, and A. Wen (2004).
Optimization of a Large-Scale Water Reservoir
Network. In Proceedings of the 2004 IE Research
Conference, Houston, TX. - Yang, Z., V. C. P. Chen, M. E. Chang, and T. E.
Murphy (2004). Formulation of a Decision-Making
Framework for Studying Ozone Pollution in Urban
Atlanta. In Proceedings of the 2004 IE Research
Conference, Houston, TX.
39References
- Chen, V. C. P., K.-L. Tsui, R. R. Barton, and M.
Meckesheimer (2002). Design, Modeling, and
Applications of Computer Experiments. IIE
Transactions on Quality and Reliability, in
review. - Tsai, J. C. C., E. K. Lee, V. C. P. Chen, and E.
L. Johnson (2003). Parallelization of the MARS
Value Function Approximation in a Decision-Making
Framework for Wastewater Treatment. Journal of
Statistical Computation and Simulation, in
review. - Cervellera, C. and V. C. P. Chen (2003). Neural
Network and Regression Spline Value Function
Approximations for Stochastic Dynamic
Programming.'' Computers and Operations
Research, in review. - Yang, Z., V. C. P. Chen, M. E. Chang, T. E.
Murphy, and J. C. C. Tsai (2004). Mining and
Modeling for a Metropolitan Atlanta Ozone
Pollution Decision-Making Framework. IIE
Transactions on Quality and Reliability, Special
Issue on Data Mining, in review.