Stoch Environ Res Risk Assess (2014) 28:167–187 DOI 10.1007/s00477-013-0738-6
ORIGINAL PAPER
An interval-parameter mean-CVaR two-stage stochastic programming approach for waste management under uncertainty C. Dai • X. H. Cai • Y. P. Cai • Q. Huo Y. Lv • G. H. Huang
•
Published online: 6 May 2013 Springer-Verlag Berlin Heidelberg 2013
Abstract In this research, approaches of interval mathematical programming, two-stage stochastic programming and conditional value-at-risk (CVaR) are incorporated within a general modeling framework, leading to an intervalparameter mean-CVaR two-stage stochastic programming (IMTSP). The developed method has several advantages: (i) it can be used to deal with uncertainties presented as interval numbers and probability distributions, (ii) its objective function simultaneously takes expected cost and system risk into consideration, thus, it is useful for helping decision makers analyze the trade-offs between cost and risk, and (iii) it can be used for supporting quantitatively evaluating the right tail of distributions of waste generation rate, which can better quantify the system risk. The IMTSP model is applied to the long-term planning of municipal solid waste management system in the City of Regina, Canada. The results indicate that IMTSP performs better in its capability
C. Dai X. H. Cai (&) Q. Huo Department of Environmental Sciences, Peking University, Beijing 100871, China e-mail:
[email protected] C. Dai e-mail:
[email protected] Q. Huo e-mail:
[email protected] Y. P. Cai (&) Y. Lv State Key Laboratory of Water Environment Simulation, School of Environment, Beijing Normal University, Beijing 100875, China e-mail:
[email protected] Y. P. Cai G. H. Huang Institute for Energy, Environment and Sustainable Communities, University of Regina, 120, 2 Research Drive, Regina, SK S4S 7H9, Canada
of generating a series of waste management patterns under different risk-aversion levels, and also providing supports for decision makers in identifying desired waste flow strategies, considering balance between system economy and environmental quality. Keywords Conditional value-at-risk Mean-risk Interval parameter Uncertainty Solid waste management
1 Introduction Municipal solid waste (MSW) management issues have increasing bearing on aesthetics, socio-economy, human health, and amenity of many communities, states and nations around the world (Kollikkathara et al. 2009). Hence, an effective planning for MSW management is important for facilitating sustainable socioeconomic development. However, various complexities exist in such an MSW management system, including the collection techniques to be used, the levels of service to be offered, and the facilities to be adopted (Dai et al. 2011). Many system parameters such as waste generation rate, facility capacity, diversion goal, and waste treatment cost may appear uncertain and cannot be expressed as deterministic values. Moreover, such uncertainties may be further multiplied by the complex features of the system components, as well as by their associations with economic penalties if the promised targets are violated (Li et al. 2009). Furthermore, there is an inherent necessary to reflect and deal with system risks associated with such uncertainties (e.g., variations in waste generation rates). Therefore, incorporation of various uncertainties, economic penalties and system risks within a general mathematical programming framework is desired for supporting MSW management and planning under such complexities.
123
168
Recently, a number of inexact optimization techniques were developed for dealing with uncertainties in planning problems, including fuzzy, stochastic, and interval mathematical programming (Ahmed 2004; Kall and Mayer 2005; Schultz and Tiedemann 2006; Qin et al. 2007; Cheng et al. 2009; Xu et al. 2009; Tan et al. 2010a, b; Shao et al. 2011; Zhu and Huang 2011; Cai et al. 2011; Noyan 2012). Among them, two-stage stochastic programming is a typical stochastic mathematical programming method. It is effective for tackling optimization problems where an analysis of policy scenarios is desired and the model’s coefficients are random with known probability distributions (Chen et al. 2012). In two-stage stochastic programming, a decision is first undertaken before values of random variables are known and, then, after the random events have happened and their values are known, a second decision is made in order to minimize ‘‘penalties’’ that may appear due to any infeasibility. The fundamental idea behind two-stage stochastic programming is the concept of recourse, which is the ability to take corrective actions after a random event has taken place (Huang and Loucks 2000). Over the past decades, two-stage stochastic programming with the expected recourse function was widely explored (Huang and Loucks 2000; Seifi and Hipel 2001; Maqsood and Huang 2003; Ahmed et al. 2004; Guo and Huang 2009a). For example, Maqsood and Huang (2003) developed an inexact two-stage stochastic programming model for planning solid waste management systems, where this method could tackle uncertainties in waste-generation rates that can be presented as probability density functions. However, the traditional two-stage stochastic programming model is risk neutral in the sense that it considers optimization as an expected criterion (Ahmed 2004). Concretely, the risk-neutral two-stage stochastic programming model takes the expected system cost as the objective function without considering risk-averting issues which may lead to excessive waste when the waste generation rate is in an extra-high condition. Therefore, a risk-averting method quantifying the effects of the variability of random waste generation rate would generate more robust solutions compared to the existing risk-neutral approaches. Conditional value-at-risk (CVaR), having appealing features such as sub-additivity and convexity, can effectively quantify risks based on known probability distributions of random variables (Rockafellar and Uryasev 2000). This approach is enhanced from a value-at-risk model, which has become an essential tool for quantifying portfolio market risks (Hsu et al. 2011). The proposed CVaR model can not only handle the expected loss of an event under extreme conditions, but also calculate the associated risks through an optimization model. More recently, many studies incorporating CVaR into constraints of two-stage stochastic programming were
123
Stoch Environ Res Risk Assess (2014) 28:167–187
conducted. Most of them focused on the introduction of CVaR into objective functions of two-stage stochastic programming models [i.e., mean-CVaR two-stage stochastic programming (MTSP)] (Noyan 2012). For example, Ahmed (2006) explored the computational suitability of various mean-risk objective functions in addressing risk in stochastic programming models, where the absolute semi-deviation risk and the quantile deviation risk measures (i.e. CVaR) were considered. Schultz and Tiedemann (2006) introduced CVaR and integer programming into a general two-stage stochastic programming framework, and transformed it into an explicit mixed-integer linear programming model when the probability distribution is discrete and finite. Fabian (2008) proposed a decomposition framework for handling a two-stage stochastic model that contains CVaR objectives and constraints. Schultz (2011) formulated a mean-CVaR stochastic integer programming model, and explored the structural properties of these optimization problems and the corresponding solution algorithms. Noyan (2012) considered a MTSP model and employed it for supporting disaster management, where uncertain information in the MTSP model was addressed through a discrete scenario method. The proposed MTSP model can not only effectively analyze predefined policy scenarios and deal with uncertainties that can be expressed as probability density functions, but also reflect expected losses associated with many extreme events and analyze trade-offs between system costs and risks. However, no application of MTSP to waste management was reported. Besides, MTSP model required probabilistic distributions for uncertain parameters, whereas in many real-world cases the quality of available information was mostly not satisfactory enough to be presented as probabilities; when uncertainties can only be obtained as interval parameters rather than probabilistic distributions, the MTSP approach may become inapplicable. In comparison, the interval mathematical programming methods were effective for handling uncertain parameters that could not be quantified as distribution functions in the model’s left- and/or righthand sides as well as in the objective function (Tan et al. 2010a). In interval mathematical programming method, interval numbers are acceptable as its uncertain inputs. Although interval mathematical programming is effective in dealing with interval numbers in objective function coefficients and constraint parameters, it may become infeasible when the right-hand sides are highly uncertain. Moreover, interval mathematical programming cannot reflect the trade-off between system cost and risk generated from the serious waste generation condition, which is common in many real-world waste management problems. Therefore, one potential approach for better accounting for the uncertainties, economic penalties, and system risk at extreme probability levels is to incorporate the interval
Stoch Environ Res Risk Assess (2014) 28:167–187
169
mathematical programming and MTSP within a general optimization framework. This leads to an interval-parameter mean-CVaR two-stage stochastic programming (IMTSP) method. The objective of this study is to develop such an IMTSP method and apply it to MSW management. The IMTSP method will be able to handle uncertainties presented as intervals and probability distributions; moreover, it can also help control the system risk and avoid problems of excessive waste. The model will then be used for planning a real MSW management system in the city of Regina, Canada. The proposed method can generate optimal solution under variable risks, which can be used for facilitating the reflection of trade-offs between system costs and risks.
min E½f ðx; xÞ þ kCVaRb ½f ðx; xÞ : x 2 X ;
ð2Þ
where CVaRb denotes the conditional value-at-risk at confidence level b; k is a non-negative weight to trade-off expected cost with risk, which is specified by decision makers corresponding to their risk preferences. Birbil et al. (2008) indicated that CVaRb ðZ þ aÞ ¼ CVaRb ðZ Þ þ a for a 2 R and Z 2 Z; which Z is a linear space of F-measurable functions on a probability space ðX; F ; PÞ: Then, we have CVaRb ½f ðx; xÞ ¼ cT x þ CVaRb ½Qðx; xÞ; and E½f ðx; xÞ þ kCVaRb ½f ðx; xÞ ¼ cT x þ E½Qðx; xÞ þ k cT x þ CVaRb ½Qðx; xÞ
ð3Þ
¼ ð1 þ kÞc x þ E½Qðx; xÞ þ kCVaRb ½Qðx; xÞ: T
2 Methodology A general form of the two-stage stochastic programming model can be formulated as follows (Birge and Louveaux 1988): min E½f ðx; xÞ ¼ min cT x þ Ex2X ½Qðx; xÞ : x 2 X ; x
x
ð1aÞ where f ðx; xÞ ¼ cT x þ Qðx; xÞ is the total cost function of the first-stage problem [i.e. model (1a)], and Qðx; xÞ ¼ min qðxÞT yðxÞ : DðxÞyðxÞ y
þTðxÞx hðxÞ; y 2 Yg
ð1bÞ
is the second-stage problem. Here c 2 R ; X R ; Y Rn2 ; x is a random variable from probability space ðX; F ; PÞ with X ! Rk ; q : X ! Rn2 ; h : X ! Rm2 ; D : X ! Rm2 n2 and T : X ! Rm2 n1 ; E denotes the expectation operates; x and y are the vectors of first-stage and second-stage decision variables, respectively. Model (1) is a risk neutral two-stage stochastic programming model in the sense that it is concerned with the optimization of an expectation objective. A common approach to address risks is to consider a weighted meanrisk criterion, where a dispersion statistic is used as a proxy for risk (Birbil et al. 2008). Markowitz (1952) developed a classical mean–variance portfolio optimization model, where variance was used as the risk measure. The disadvantage of this mean–variance formulation is that it considers the under-and-over-performances equally. When a typical dispersion statistics such as variance is used as risk measures, the mean-risk approach may lead to inferior solutions (Noyan 2012). CVaR as a downside risk measure can be proposed to remedy this drawback (Ahmed 2006). Therefore, we have the following meanCVaR model: n1
n1
Obviously, the essential issue for obtaining the numerical solutions of model (2) is to calculate the minimization of CVaRb ½Qðx; xÞ firstly. Let Qðx; xÞ denote a loss function associated with decision vector x 2 X Rn and random vector x 2 Rm : For convenience, the underling probability of x will be assumed to have a density function pðÞ (Rockafellar and Uryasev 2002). Suppose that Qðx; xÞ is continuous in x and measurable in x; and that EfjQðx; xÞjg\1 for each x 2 X. For a confidence level b and a fixed decision variable x; the value-at-risk, denoted by VaRb ðxÞ, can be defined as8follows: 9 > > Z = < ð4Þ pðxÞdx b VaRb ðxÞ ¼ min a 2 R : > > ; : Qðx;xÞ a
where
R
Qðx;xÞ a
pðxÞdx denotes the probability of Qðx; xÞ
not exceeding a threshold a. CVaRb ðxÞ is defined as the expected value of the loss that exceeds VaRb ðxÞ; which can be presented as follows: Z 1 CVaRb ð xÞ ¼ Qðx; xÞpðxÞdx: ð5Þ 1 b Qðx;xÞ VaRb ðxÞ
The CVaR is a coherent risk measure (Ahmed 2006). The problem involved CVaRb ðxÞ is different to proceed because of its convoluted and implicit version. Rockafellar and Uryasev (2000) developed the following function: Z 1 Fb ðx; aÞ ¼ a þ fmaxf0; ½Qðx; xÞ agg 1 b x2Rm
pðxÞdx;
ð6Þ
where Fb ðx; aÞ is shown to be convex and continuously differentiable with respect to a; and minx CVaRb ðxÞ ¼ mina Fb ðx; aÞ . Through introducing an auxiliary variable z; the minimization of CVaRb is equivalent to the following model (Tong et al. 2010):
123
170
Stoch Environ Res Risk Assess (2014) 28:167–187
min a þ
1 Ex2Rm ½zðxÞ ; 1 b
x
ð7aÞ
subject to x 2 X;
ð7bÞ
zðxÞ Qðx; xÞ a;
ð7cÞ
zðxÞ 0:
ð7dÞ
By coupling model (1) and model (7) within the mean-CVaR model [i.e., model (2)], a mean-CVaR two-stage stochastic programming (MTSP) model can be formulated as follows: min f ¼ ð1 þkÞcT x þ E½Qðx; xÞ 1 Ex2Rm ½zðxÞ þk aþ 1 b
ð8aÞ
subject to zðxÞ Qðx; xÞ a
ð8bÞ
zðxÞ 0
ð8cÞ
a0
ð8dÞ
x2X
ð8eÞ
with T
Qðx; xÞ ¼ min qðxÞ yðxÞ
ð8fÞ
subject to DðxÞyðxÞ þ TðxÞx hðxÞ
ð8gÞ
y2Y
ð8hÞ
By letting random variables (i.e. x) take discrete values P ps ¼ 1Þ; xs at probability levels ps (s ¼ 1; 2; . . .v and the above MTSP model can be equivalently formulated as the following linear programming model (Noyan 2012): v X
ps ðqs ÞT ys s¼1 ! v 1 X þk aþ ps z s 1 b s¼1
min f ¼ ð1 þ kÞcT x þ
ð9aÞ
subject to Ts x þ Ds ys hs ;
s ¼ 1; . . .; v
x2X ys 0;
ð9bÞ ð9cÞ
s ¼ 1; . . .; v
zs ðqs ÞT ys a; a 2 R; zs 0;
123
s ¼ 1; . . .; v s ¼ 1; . . .; v
ð9dÞ ð9eÞ ð9fÞ
Obviously, model (9) can deal with uncertainties in the right-hand sides presented as probabilistic distributions when coefficients in the left-hand sides and the objective function are deterministic. However, in real-world optimization problems, the quality of information that can be obtained is mostly not satisfactory enough to be presented as probabilities (Li et al. 2009). Such complexities cannot be solved through model (9). Interval mathematical programming is effective in tackling uncertainties expressed as interval values with known lower and upper bounds but unknown distribution functions (Huang et al. 1992). Therefore, through incorporating interval mathematical programming and MTSP within a general optimization framework, an interval-parameter MTSP (or IMTSP) method can be formulated as follows: v X T T
min f ¼ ð1 þ kÞ c x þ ps q
ys s s¼1 ! v 1 X þ k a þ ps z s 1 b s¼1
ð10aÞ
subject to
T
s ¼ 1; . . .; v s x þ Ds ys hs ; T
zs qs ys a ; s ¼ 1; . . .; v a 2 R ; zs 0; s ¼ 1; . . .; v
ð10bÞ
x2X
ð10eÞ
ys 0;
s ¼ 1; . . .; v
ð10cÞ ð10dÞ
ð10fÞ
where superscript ‘‘±’’ means interval-valued feature; the ‘‘-’’ and ‘‘?’’ superscripts represent lower and upper bounds of an interval parameter/variable, respectively. The objective function value of model (10) includes the expected cost and risk cost of excess waste amount. Model (10) can not only effectively reflect random uncertainty but also guarantee solutions to be more stable and reliable. Model (10) can be transformed into two deterministic submodels that correspond to the lower and upper bounds of the desired objective function value. This transformation process is based on an interactive algorithm, which is different from the best/worst case analysis (Huang et al. 1992). Interval solutions associated with varying levels of constraint–violation risks can then be obtained by solving the two submodels sequentially. The submodel corresponding to the lower-bound objective function value ðf Þ can be firstly formulated as follows:
Stoch Environ Res Risk Assess (2014) 28:167–187
min f
k1 X
¼ ð 1 þ kÞ
c j xj
þ
þ
s¼1
!
n2 X
ps q js yjs þ
k1 n1
X X tjs Sign t xþ þ tjs þ Sign tþ x js j js j
þ c j xj
j¼k1 þ1
j¼1 k2 v X X
n1 X
171
k2
X djs Sign d yþ þ
þ ps q js yjs
js
j¼k2 þ1
j¼1
v 1 X þk g þ ps z s 1 a s¼1
js
j¼k2 þ1
k2 n2
X X qjs Sign q y þ qjs Sign q yþ zs js js js js j¼k2 þ1
8s 8s
g 0
ð11gÞ
0;
8s;
j ¼ 1; 2; . . .; k2
ð11hÞ
yþ js 0;
8s;
j ¼ k2 þ 1; k2 þ 2; . . .; n2
ð11iÞ
where x
j ; j ¼ 1; 2; . . .; k1 , are interval variables with positive coefficients in the objective function; x
j ; j ¼ k1 þ 1; k1 þ 2; . . .; n1 , are interval variables with negative coefficients; y
js ; j ¼ 1; 2; . . .; k2 and s ¼ 1; 2; . . .; v, are interval random variables with positive coefficients in the objective function; y
js 0; j ¼ k2 þ 1; k2 þ 2; . . .; n2 and s ¼ 1; 2; . . .; v, are interval random variables with neg ative coefficients. Solutions of g opt , xj;opt (j ¼ 1; 2; . . .; k1 ), xþ j;opt (j ¼ k1 þ 1; k1 þ 2; . . .; n1 ), yjs;opt (j ¼ 1; 2; . . .; k2 ), þ and yjs;opt (j ¼ k2 þ 1; k2 þ 2; . . .; n2 ) can be obtained through submodel (11). Based on the above solutions, the second submodel for f þ can be formulated ! as follows: k1 n1 X X þ cþ cþ min f þ ¼ð1 þ kÞ j xj þ j xj
j¼1
þ ps qþ js yjs þ
subject to
! ps qþ js yjs
j¼k2 þ1
v 1 X þk g þ ps z s 1 a s¼1 þ
js
8s
ð12cÞ
8s
ð12dÞ
gþ g opt xþ j xj;opt ;
ð12eÞ j ¼ 1; 2; . . .; k1
þ 0 x j xj;opt ;
j ¼ k1 þ 1; k1 þ 2; . . .; n1
ð12fÞ ð12gÞ
þ Solutions of gþ opt ; xj;opt (j ¼ 1; 2; . . .; k1 Þ; xj;opt (j ¼ k1 þ 1; þ k1 þ 2; . . .; n1 Þ; yjs;opt (j ¼ 1; 2; . . .; k2 Þ; and y js;opt (j ¼ k2 þ 1; k2 þ 2; . . .; n2 ) can be obtained through submodel (12). Through integrating solutions of submodels (11) and (12), interval solution for model (10) under a set of ps (s ¼ 1; 2; . . .; v) levels can be expressed as follows: h i h i þ þ x
8j; y
8j; s; j;opt ¼ xj;opt ; xj;opt ; js;opt ¼ yjs;opt ; yjs;opt ; h i
þ fopt ¼ fopt ; fopt . Figure 1 shows the flow chart for the for-
j ¼ k1 þ 1; k1 þ 2; . . .; n1
n2 X
þ gþ ; zs 0;
js
j¼k2 þ1
ð11eÞ
y js
j¼k1 þ1
js
þ 0 y js yjs;opt ;
xþ j 0;
s¼1
8s
ð11dÞ
ð11fÞ
þ
s
yþ js yjs;opt ;
j ¼ 1; 2; . . .; k1
j¼1
js
ð11cÞ
x j 0;
k2 v X X
js
j¼1
ð11bÞ
zs 0;
þ djs Sign dþ y hþ ;
k2 n2
X X qjs þ Sign qþ yþ þ qjs þ Sign qþ y zs
j¼k1 þ1
þ g ;
ð12bÞ
js
j¼k2 þ1
k2 n2
X X djs þ Sign d þ y þ djs Sign d yþ h 8s þ js js js js s
j¼1
þ
ð11aÞ
subject to k1 n1
X X tjs þ Sign tþ x þ tjs Sign t xþ js j js j
j¼1
j¼1 n2 X
!
j¼1
j¼k1 þ1
j¼1
!
!
ð12aÞ
8s;
j ¼ 1; 2; . . .; k2 8s;
j ¼ k2 þ 1; k2 þ 2; . . .; n2
ð12hÞ ð12iÞ
mulation of the IMTSP method and its method of solution.
3 Case study The City of Regina is located in southeast Saskatchewan and is the provincial capital with a population of approximately 179,000 within the city limits. The amount of waste generated by the residential sector is considerably smaller than the amount of waste generated by the industrial, commercial and institutional (IC&I) and construction and demolition (C&D) sectors, which contributes about 70 % of the waste compared with 30 % for the residential sector (Dai et al. 2012). In 2008, the residential sector generated a total of approximately 68,000 tonnes of solid waste and of this diverted about 11,000 tonnes (16.5 %) of material through recycling or returned for deposit beverage container program (City of Regina 2011). For the long-term planning exercise, a general waste generation rate of [1.00, 1.17] kg/capita/day is used for the residential sector of the City (Li and Huang 2009). Residential waste collection
123
172
Stoch Environ Res Risk Assess (2014) 28:167–187
Fig. 1 General framework of the IMTSP method
min
service is provided to 52,000 households (76 % of the total) in Regina, as part of the general city services that are funded by property taxes. Collection of solid waste from residences is the most costly and visible part in MSW management system. The average frequency of collection is around once per week. The average garbage collection/ transportation cost was [56, 71] $/tonne for manual collection or [29, 43] $/ton for the automatic one. Costs for waste collection and transportation are estimated based on the existing conditions in the collection areas; the average container size, collection frequency, collection mode (automatic and manual), and collection time (per load) were considered. In this study, a projected interval of [38, 47] $/ton is used for approximating the cost of waste collection and transportation (Guo and Huang 2009b). The main approaches for waste treatment include recycling, composting and landfill. Recycling is the high profile activity that the public most identifies in terms of sustainable development and waste minimization. Therefore, the city operates many recycling programs to encourage residents to reduce the amounts of waste that end up at the landfill (Li and Huang 2009). For example, a big blue bin project was started in 1990. It provides residents with an option to recycle paper and cardboard. The amount of paper and cardboard recycled was enough to make 88,524,588 cartons in 2008. Over the last three years, the revenue from the sale
123
{Ε ⎡⎣ f ( x, ω )⎤⎦ + λρ ⎡⎣ f (x,ω )⎤⎦ : x ∈ X }
of paper has declined from $40/tonne to $35/tonne in 2007 and 2008, while the operation cost of the program has increased (City of Regina 2011). A white metal goods recycling program was operated for collecting and selling the scrap metal. Between 2006 and 2008, this program has diverted approximately 2,700 tonnes of material and generated $246,000 in gross revenue. Potentially recyclable material constituted approximately 41.79 % of the waste stream, while the observed recycling rate of residential waste was only 8.72 % in Regina, indicating a high potential for further improvements (Dai et al. 2012). Composting is the conversion of organic waste into a soil product for use elsewhere (Li and Huang 2009). The city is operating a backyard composting program to naturally process grass clippings (by leaving them on the lawns). The waste stream generated in the city contained a high percentage of compostable materials, including 20,000 tonnes of yard wastes (25.64 %) and 20,900 tonnes of organic wastes (26.79 %); the potentially compostable material constituted approximately 52.43 % of the waste stream (City of Regina 2011). Composting on-site removes this valuable material from the waste stream, saves space in the landfill and reduces the use of fossil fuels for waste transport. Composting returns nutrients and moisture to the soil, improves soil consistency, and impedes weed growth, all without the use of herbicides, fertilizers, and excessive watering. Revenue can be obtained
Stoch Environ Res Risk Assess (2014) 28:167–187
173
Table 1 Waste generation rates under different probability levels Level of waste generation
Scenario
Period 1 (k = 1)
Period 2 (k = 2)
Period 3 (k = 3)
Probability
Probability
Probability
Waste generation rate (tonne/week)
Waste generation rate (tonne/week)
Waste generation rate (tonne/week)
Extra low
S1
0.012
[1038, 1167]
0.019
[1114, 1237]
0.005
[1171, 1288]
Very low
S2
0.05
[1167, 1295]
0.066
[1237, 1359]
0.049
[1288, 1405]
Low Low to medium
S3 S4
0.161 0.248
[1295, 1424] [1424, 1552]
0.192 0.23
[1359, 1481] [1481, 1604]
0.138 0.252
[1405, 1522] [1522, 1639]
Medium
S5
0.26
[1552, 1680]
0.223
[1604, 1726]
0.28
[1639, 1755]
Medium to high
S6
0.178
[1680, 1809]
0.157
[1726, 1848]
0.173
[1755, 1872]
High
S7
0.077
[1809, 1937]
0.081
[1848, 1971]
0.079
[1872, 1989]
Very high
S8
0.01
[1937, 2066]
0.029
[1971, 2093]
0.018
[1989, 2106]
Extra high
S9
0.004
[2066, 2194]
0.003
[2093, 2215]
0.006
[2106, 2223]
from compost produced by composting technology. The residue ([8, 10] % of original waste) from composting facilities will be finally disposed by landfills. The landfill is located in the northeast quadrant of the city at the corner of Fleet Street and McDonald Street. It is the only solid waste landfill in the study region and which occupies 97 ha, with an actual landfilling area of 60 ha. In 2008, Regina’s citizens generate about 1,000 kg per capita of MSW, more than 80 % of which is disposed in the landfill with less than 20 % diverted (City of Regina 2011). Landfill deliveries have increased from 2005 to 2008 with 478,474 tonnes in 2005 rising to 645,456 tonnes in 2008, where the IC&I sector and C&D sector produced 61 % of the total waste flow, and the residential sector generated approximately 39 % (Dai et al. 2012). The total number of vehicle transactions has increased over the last 4 years from 137,128 in 2006 to 178,097 vehicles in 2008, with an average load of [3.48, 3.62] tonnes, respectively. In Canada, establishment of waste diversion targets and relevant regulations were currently a growing trend. In March 2010, that city’s Council adopts, in principle, the Waste Plan Regina’s Residential Option, and set a waste diversion target of 34 % by 2020, and 48 % by 2025 (City of Regina 2011). The study time is 15 years (from 2012 to 2027), which is further divided into three 5-year periods. Table 1 presents the waste-generation rates and the associated probabilities of occurrence in the three planning periods, indicating that the waste-generation amounts are highly uncertain, presented by intervals with the associated probabilities. Tables 2 and 3 contain regular costs for allowable waste flows, operating costs for waste management facilities, and penalty costs for surplus waste flows, as well as revenues from waste management facilities over the three planning periods (Li and Huang 2009; City of Regina 2011; Dai et al. 2012). Table 4 shows the capacityexpansion options and the relevant capital costs for composting and recycling facilities (City of Regina 2011). It is
also indicated that the expansion costs decrease along with time. Since the planning problem under consideration is dynamic with multiple stages, discount factors are necessary for each period to obtain a total present value for the object function. In this project, all the cash flows are counted in year 2012 dollars. Because uncertainties and risks exist in a variety of system components and a link to the pre-regulated policy as formulated by local authorities is desired, the IMTSP method developed in Sect. 2 is considered to be suitable for tackling this type of management problem. The model includes continuous and binary decision variables. The binary variables represent the development or expansion options for waste-management-facilities in different peri
ods (i.e. BC
ck and BRrk ); their solutions can be used for answering the questions related to timing, sizing and siting Table 2 Costs and revenues for allowable waste flows Planning period k=1
k=2
k=3
Operating costs of waste-management-facilities, OP
ik ($/tonne) Landfill
[6.9, 13.02]
[5.83, 11.02]
[4.94, 9.32]
Composting
[16.09, 19.92]
[13.61, 16.85]
[11.52, 14.26]
Recycling
[46.74, 51.95]
[39.54, 43.95]
[33.45, 37.18]
Collection and transportation costs, TR
ik ($/tonne) To landfill
[24.52, 28.35]
[20.74, 23.98]
[17.55, 20.29]
To composting
[55.93, 55.93]
[47.32, 47.32]
[40.03, 40.03]
To recycling
[77.38, 77.38]
[65.47, 65.47]
[55.39, 55.39]
Revenues from waste-management-facilities, RE
ik ($/tonne) Composting
[3.83, 7.66]
[3.24, 6.48]
[2.74, 5.48]
Recycling
[34.48, 42.14]
[29.17, 35.65]
[24.68, 30.16]
Costs of residue disposal at the landfill,
FT
ik
($/tonne)
From composting
[1.29, 1.61]
[1.09, 1.36]
[0.92, 1.15]
From recycling
[1.13, 1.29]
[0.95, 1.09]
[0.81, 0.92]
123
174
Stoch Environ Res Risk Assess (2014) 28:167–187
Table 3 Costs and revenues for excess waste flows
Planning period k=1
k=2
k=3
Operating costs of waste-management-facilities, DP
ik ($/tonne) Landfill Composting
[13.79, 26.05] [26.05, 32.18]
[11.67, 22.04] [22.04, 27.23]
[9.87, 18.65] [18.65, 23.03]
Recycling
[79.69, 88.34]
[67.41, 74.74]
[57.04, 63.23]
Collection and transportation costs, DR
ik ($/tonne) To landfill
[36.77, 42.52]
[31.11, 35.98]
[26.32, 30.27]
To composting
[83.89, 83.89]
[70.98, 70.98]
[60.05, 60.05]
To recycling
[116.07, 116.07]
[98.2, 98.2]
[83.08, 83.08]
Revenues from waste-management-facilities, RM
ik ($/tonne) Composting
[3.83, 7.66]
[3.24, 6.48]
[2.74, 5.48]
Recycling
[34.48, 42.14]
[29.17, 35.65]
[24.68, 30.16]
Costs of residue disposal at the landfill, DT
ik ($/tonne) From composting
[1.93, 2.41]
[1.63, 2.04]
[1.38, 1.73]
From recycling
[1.69, 1.93]
[1.43, 1.63]
[1.21, 1.38]
for waste-management-facility development and/or expansion under uncertainty. The continuous variables represent the optimized waste flows from the city to the waste-management-facilities. Furthermore, the continuous variables include two subsets: those (the first-stage ones, X
ik ) that must be determined before the random variables (i.e. waste-generation rates) are disclosed, and those (the
) that will be determined after the second-stage ones, Yiks random variables are disclosed. The modeling formulation based on the IMTSP method for the city’s MSW management can be presented as follows: " 3 X 3 X
min f ¼ ð1 þ kÞ Lk X
ik TRik þ OPik i¼1 k¼1
þ
3 X 3 X
Lk X
ik FEi FTik REik
"
zs þ a
3 X 3 X
Lk Y
iks DRik þ DPik
i¼1 k¼1
þ
3 X 3 X
L k Y
iks
FE
i DTik
3 X 3 X 9 X
þ
þ ð 1 þ kÞ
ð13cÞ (Total landfill capacity constraints) þ
Y
2ks
CC þ
þ
r¼1 k¼1
CRE
rk BRrk
123
DCC
ck BCck ;
c¼1 k¼1
ð13dÞ
(Composting capacity constraints) 0
X
3k
þ
Y
3ks 0
RC þ
3 X k X
DRC
rk BRrk ;
r¼1 k¼1
ð13eÞ
8s; k; k ¼ 1; 2; 3
9 1 X ps z s þ k a þ 1 b s¼1
!
ð13aÞ subject to
3 X k X
8s; k; k ¼ 1; 2; 3
CCE
ck BCck
#
ð13bÞ
8s
0
p s Lk Y
iks FEi DTik RMik
c¼1 k¼1 3 X 3 X
8s
0
p s Lk Y
iks DRik þ DPik
3 X 3 X
;
i¼2
k¼1
X
2k
i¼2 k¼1 s¼1
"
#
(Auxiliary variable of CVaR constraints) " # 3 3 X
X
Lk X1k þ Y1ks þ FEi Xik þ Yiks LC ;
#
i¼1 k¼1 s¼1 3 X 3 X 9 X
i¼2 k¼1
i¼2 k¼1
þ
RM
ik
(Recycling capacity constraints)
X
2k þ Y2ks DCk Wks ;
8k; s
ð13fÞ
X
3k þ Y3ks DRk Wks ;
8k; s
ð13gÞ
(Constraints of waste flow to the composting and recycling)
Stoch Environ Res Risk Assess (2014) 28:167–187
175
Table 4 Capacity expansion options and the relevant capacity costs Expansion option
Total expanded capacity
Expanded capacity for residential
Expansion cost (106) k=1
k=2
k=3
Composting facility, DCC
ck (tonne/week) Option 1 (c = 1)
189
132
[1.92, 2.31]
[1.62, 1.95]
[1.37, 1.65]
Option 2 (c = 2)
385
270
[3.83, 4.61]
[3.24, 3.89]
[2.74, 3.29]
Option 3 (c = 3)
483
338
[4.79, 5.57]
[4.05, 4.86]
[3.43, 4.12]
Recycling facility, DRC
rk (tonne/week) Option 1 (r = 1) Option 2 (r = 2)
350 700
140 280
[4.08, 4.53] [8.17, 9.06]
[3.46, 3.83] [6.91, 7.67]
[2.92, 3.24] [5.85, 6.49]
Option 3 (r = 3)
875
350
[9.83, 11.33]
[8.32, 9.58]
[7.04, 8.11]
3 X
X
ik
þ
Y
iks
W
ks ;
8k; s
ð13hÞ
i¼1
(Waste disposal demand constraints) a0
ð13iÞ
zs 0;
8s
0 Y
iks
X
ik
ð13jÞ Xik;max ;
8i; k; s
ð13kÞ
(Non-negativity and technical constraints)
( BC
ck ( BR
rk
¼ 1;
if capacity expansion for composting facility is undertaken
¼ 0;
if otherwise
¼ 1; ¼ 0;
if capacity expansion for recycling facility is undertaken ; if otherwise
(Binary constraints) 3 X
BC
ck 1;
8k
ð13nÞ
BR
rk 1;
8k
ð13oÞ
c¼1 3 X
flow to facility i during period k ($/tonne), i = 1, 2, 3; DR
ik is cost of collecting and transporting excess waste flow from the city to facility i during period k ($/tonne) (the
second-stage cost parameter), where DR
ik TRik , i = 1, 2, 3;
OPik is regular operating cost of waste management facility i during period k ($/tonne), i = 1, 2, 3; DP
ik is operating cost of facility i for excess waste flow during period k ($/ tonne) (the second-stage cost parameter), where
DP
ik OPik , i = 1, 2, 3; FEi is residue flow from facility i to the landfill (% of incoming mass to facility i), i = 1, 2, 3;
r¼1
(Expansions for composting and recycling facilities may occur in any given time period) where f is mean-risk function value ($); k and b are risk parameters; Lk is length of time period k (week); k is time period, k = 1, 2, 3; c is the name of expansion option for composting facility, c = 1, 2, 3; r is the name of expansion option for recycling facility, r = 1, 2, 3; s is the level of waste generation, s = 1, 2,…, 9; i is type of waste management facility, i = 1, 2, 3, where i = 1 for the landfill, 2 for composting facility, and 3 for recycling facility; TR
ik is cost of collection and transportation for allowable-waste
;
8c; k
8r; k
ð13lÞ
ð13mÞ
FT
ik is disposal cost for allowable residues generated by facility i during period k ($/tonne), i = 2, 3; DT
ik is disposal cost for excess waste residues generated by facility i during period k ($/tonne) (the second-stage cost param
eter), where DT
ik FTik , i = 2, 3; REik is revenue generated by processing allowable waste flows in facility i during period k ($/tonne), i = 2, 3; RM
ik is revenue generated by processing excess waste flows in facility i during period k ($/tonne), (the second-stage revenue parameter) i = 2, 3; W
sk is residential waste-generation rate with probability ps in period k (tonne/week); ps is probability of waste generation rate (W
ks ) with level s (%);
Xik is allowable-waste flow to facility i during period k (tonne/week) (the first-stage decision variable), i = 1, 2, 3; Xik;max is maximum allowable-waste flow from the city to facility i during period k (tonne/week); Y
isk is amount by which the allowable-waste level is exceeded when the
123
176
waste-generation rate is W
sk with probability ps (tonne/ week) (the second-stage decision variable); LC is existing landfill capacity (tonne); CC is existing composting facility capacity (tonne/week); RC is existing recycling facility capacity (tonne/week); DC
k is waste flow transferred from composting facility in period k; DR
k is waste flow transferred from recycling facility in period k; DCC
ck is available expanded/developed capacity for composting facility with expansion option c in period k (tonne/week); DRC
rk is available expanded/developed capacity for recycling facility with expansion option r in period k (tonne/ week); BC
ck is binary decision variable for expanding composting facility with option c at the start of period k; BR
rk is binary decision variable for expanding recycling facility with option r at the start of period k; CCE
ck is capital cost of expanding composting facility by option c in period k ($/tonne); CRE
rk is capital cost of expanding recycling facility by option c in period k ($/tonne).
4 Result analysis In this study, there are two risk-related parameters: k and b, where three b levels (0.8, 0.9 and 0.99) were examined under each k level. Table 5 shows the solutions of wasteflow allocation over periods 1–3 under b = 0.99. An excessive waste flow would be generated if the allowablewaste-flow level is exceeded (i.e. excessive flow = generated flow - assigned quota). The waste-flow patterns (including allowable and excess flows) would vary dynamically due to temporal and spatial variations in wastegeneration/management conditions. When the level of waste-generation rate is medium to high in period 1, the optimized allowable flows to the landfill would be 1,188 tonnes/week under k = 0.5, 1397 tonnes/week under k = 1, and [1557, 1698] tonnes/week under k = 5. The optimal excessive flows would be [141, 283] tonnes/week under k = 0.5, [0, 142] tonnes/week under k = 1 and 0 tonnes/week under k = 5, respectively. Thus, the total waste flows to the landfill would be [1329, 1471], [1397, 1539], and [1557, 1698] tonnes/week under k = 0.5, 1 and 5, respectively. Figure 2 presents the optimal waste-allocation patterns from the city to the landfill, composting and recycling facilities under different risk parameters. The results indicate that a plan for waste-flow allocation would be related to both waste-generation rate and capacity-expansion scheme under varied risk parameters. Figure 3 shows the varying trend of the excessive waste flows allocated to the landfill, composting and recycling facilities under different k and b levels. In the case of excess waste, the allotment to the landfill should be assigned initially, and
123
Stoch Environ Res Risk Assess (2014) 28:167–187
then to the incinerating and composting facilities. Since the incinerating and composting facilities have higher regular and penalty costs for treating wastes, allotment of wastes to the landfill would be more economical. Moreover, the related outputs are affected significantly by the risk parameters. For example, the excessive waste flows to any waste management facilities would increase when k is fixed and b changes from 0.8 to 0.99 or b is fixed and k changes from 0.5 to 5. Therefore, decreasing the value of k and/or b means the more excess waste flow, which may lead to a higher level of system risk; conversely, increasing the value of k and/or b means the less excess waste flow, which may lead to a lower level of system risk. Table 6 presents the solutions of facility expansion schemes under different risk parameters. When b = 0.8 or 0.9, the expansion schemes for composting and recycling facilities would be constant even though the weight k vary between 0.5 and 5. In comparison, if b = 0.99, they would vary under different k levels. For example, when k = 0.5, a centralized composting facility would be expanded by an increment of 483 tonnes/week at the starts of periods 1 and 2 respectively, with 70 % of the capacity (i.e. 338 tonnes/ week) being dedicated to the residential waste since the IC&I (industrial, commercial and institutional) and rural sectors would not be included within the scope of this study. This facility should be developed with a capacity of 189 tonnes/week at the start of period 3. When k = 1 and 5, a centralized composting facility would be developed with a capacity of 483 tonnes/week at the starts of periods 1 and 2 respectively; then this facility should be expanded by an increment of 385 tonnes/week at the start of periods 3. Consequently, the total expanded capacities for the composting facility would be 1155, 1351 and 1351 tonnes/ week under k = 0.5, 1 and 5, respectively. Comparatively, the recycling facility would be expanded once with an increment of 350 tonnes/week at the start of period 1 when k = 0.5, and expanded with the same capacity at the start of period 2 when k = 1 and 5, respectively. Figure 4 shows the mean-risk function value of f for different risk parameters. The results indicate that the optimal mean-risk function would increase when the risk parameters k or b increase. For example, the mean-risk function value would increase from $[96, 122] 9 106 to $[747, 941] 9 106 when b is fixed at 0.8 and k increases from k = 0.5 to k = 10. Besides, the mean-risk function value would increase from $[96, 122] 9 106 to $[107, 133] 9 106 when k is fixed at 0.5 and b increases from b = 0.8 to b = 0.99, respectively. Figure 5 presents the value of CVaR for different k and b levels. Generally, CVaR would increase as b increases. When k is fixed at 5 and b changes from b = 0.8–0.99, the value of CVaR would increase from $[69.7, 88.4] 9 106 to $[90.8,
s=8
s=9
s=9
s=9
3
1
2
3
s=8
s=8
1
s=7
3
2
s=6
s=7 s=7
3
1 2
s=6
s=6
1
s=5
3
2
s=5
s=5
1
2
s=4
s=4
2
s=4
1
3
s=3
s=3
2
3
s=2
s=3
3
1
s=1
s=2 s=2
3
1 2
s=1
s=1
1
k=1
Level of waste generation
2
Facility (i)
Period
0.004
0.004
0.004
0.015
0.015
0.015
0.064
0.064 0.064
0.187
0.187
0.187
0.298
0.298
0.298
0.251
0.251
0.251
0.138
0.138
0.138
0.037
0.037 0.037
0.006
0.006
0.006
Probability
136
198
1,188
136
198
1,188
136
1,188 198
136
198
1,188
136
198
1,188
136
198
1,188
136
198
1,188
136
1,188 198
136
198
1,188
0
136
198
[232, 374]
0
179
[246, 387]
0
[283, 425] 0
0
0
[141, 283]
0
0
[0, 141]
0
0
0
0
0
0
0
0 0
0
0
1,188
272
396
[1420, 1562]
136
377
[1434, 1575]
136
[1471, 1613] 198
136
198
[1329, 1471]
136
198
[1188, 1329]
136
198
1,188
136
198
1,188
136
1,188 198
136
198
1,397
68
198
1397
68
198
1,397
68
1,397 198
68
198
1,397
68
198
1,397
68
198
1,397
68
198
1,397
68
1,397 198
68
198
68
198
[159, 301]
0
110
[174, 315]
0
[142, 284] 0
0
0
[0, 142]
0
0
0
0
0
0
0
0
0
0
0 0
0
0
0
Excess waste flow
Allowable waste flow
Optimized waste flow
Allowable waste flow
Excess waste flow
k = 1 (tonne/week)
k = 0.5 (tonne/week)
Table 5 Solution of the IMTSP method for continuous variables for b = 0.99
136
396
[1556, 1698]
68
308
[1571, 1712]
68
[1539, 1681] 198
68
198
[1397, 1539]
68
198
1,397
68
198
1,397
68
198
1,397
68
1,397 198
68
198
1,397
Optimized waste flow
68
289
[1557, 1698]
68
289
[1557, 1698]
68
[1557, 1698] 289
68
289
[1557, 1698]
68
289
[1557, 1698]
68
289
[1557, 1698]
68
289
[1557, 1698]
68
[1557, 1698] 289
68
289
[1557, 1698]
Allowable waste flow
68
107
[0, 1]
0
0
34
0
0 0
0
0
0
0
0
0
0
0
0
0
0
0
0
0 0
0
0
0
Excess waste flow
k = 5 (tonne/week)
[136, 135]
396
[1557, 1699]
68
289
[1591, 1732]
68
[1557, 1698] 289
68
289
[1557, 1698]
68
289
[1557, 1698]
68
289
[1557, 1698]
68
289
[1557, 1698]
68
[1557, 1698] 289
68
289
[1557, 1698]
Optimized waste flow
Stoch Environ Res Risk Assess (2014) 28:167–187 177
123
123
s=8
s=8
s=9
s=9
s=9
1
2
3
s=7 s=8
3 1
2
s=7
2
3
s=6
s=7
s=6
2
3
s=6
1
1
s=5
s=5
2
3
s=4
s=5
3
1
s=4
s=4
1
2
s=3
s=3
2
3
s=2
s=2 s=3
2
3 1
s=1
s=2
3
1
s=1
s=1
1
k=2
Level of waste generation
2
Facility (i)
Period
Table 5 continued
0.004
0.004
0.004
0.015
0.015
0.064 0.015
0.064
0.064
0.187
0.187
0.187
0.298
0.298
0.298
0.251
0.251
0.251
0.138
0.138
0.037 0.138
0.037
0.037
0.006
0.006
0.006
Probability
158
367
1,047
158
367
158 1,047
367
1,047
158
367
1,047
158
367
1,047
158
367
1,047
158
367
158 1,047
367
1,047
158
367
1,047
0
158
367
[60, 206]
0
367
0 [72, 218]
198
[94, 241]
0
0
[146, 292]
0
0
[0, 146]
0
0
0
0
0
0 0
0
0
0
0
1,047
316
734
[1107, 1253]
158
734
158 [1119, 1265]
565
[1141, 1288]
158
367
[1193, 1339]
158
367
[1047, 1193]
158
367
1,047
158
367
158 1,047
367
1,047
158
367
[1107, 1193]
158
367
[1107, 1193]
158
367
158 [1107, 1193]
367
[1107, 1193]
158
367
[1107, 1193]
158
367
[1107, 1193]
158
367
[1107, 1193]
158
367
158 [1107, 1193]
367
[1107, 1193]
158
367
158
367
[0, 60]
0
367
0 [12, 72]
142
[90, 151]
0
0
[86, 146]
0
0
0
0
0
0
0
0
0 0
0
0
0
0
0
Excess waste flow
Allowable waste flow
Optimized waste flow
Allowable waste flow
Excess waste flow
k = 1 (tonne/week)
k = 0.5 (tonne/week)
316
734
[1107, 1253]
158
734
158 [1119, 1265]
509
[1197, 1344]
158
367
[1193, 1339]
158
367
[1107, 1193]
158
367
[1107, 1193]
158
367
158 [1107, 1193]
367
[1107, 1193]
158
367
[1107, 1193]
Optimized waste flow
171
733
[1107, 1253]
171
733
171 [1107, 1253]
733
[1107, 1253]
171
733
[1107, 1253]
171
733
[1107, 1253]
171
733
[1107, 1253]
171
733
171 [1107, 1253]
733
[1107, 1253]
171
733
[1107, 1253]
Allowable waste flow
145
1
0
0
0
0 0
0
0
0
0
0
0
0
0
0
0
0
0
0
0 0
0
0
0
0
0
Excess waste flow
k = 5 (tonne/week)
316
734
[1107, 1253]
171
733
171 [1107, 1253]
733
[1107, 1253]
171
733
[1107, 1253]
171
733
[1107, 1253]
171
733
[1107, 1253]
171
733
171 [1107, 1253]
733
[1107, 1253]
171
733
[1107, 1253]
Optimized waste flow
178 Stoch Environ Res Risk Assess (2014) 28:167–187
s=8
s=8
s=9
s=9
s=9
1
2
3
s=7 s=8
3 1
2
s=7
2
3
s=6
s=7
s=6
2
3
s=6
1
1
s=5
s=5
2
3
s=4
s=5
3
1
s=4
s=4
1
2
s=3
s=3
2
3
s=2
s=2 s=3
2
3 1
s=1
s=2
3
1
s=1
s=1
1
k=3
Level of waste generation
2
Facility (i)
Period
Table 5 continued
0.004
0.004
0.004
0.015
0.015
0.064 0.015
0.064
0.064
0.187
0.187
0.187
0.298
0.298
0.298
0.251
0.251
0.251
0.138
0.138
0.037 0.138
0.037
0.037
0.006
0.006
0.006
Probability
158
433
1,033
158
433
158 1,033
433
1,033
158
433
1,033
158
433
1,033
158
433
1,033
158
433
158 1,033
433
1,033
158
433
1,033
0
158
433
[13, 164]
0
433
0 [20, 171]
302
[0, 151]
0
23
[128, 279]
0
0
[0, 151]
0
0
0
0
0
0 0
0
0
0
0
1,033
316
866
[1046, 1197]
158
866
158 [1053, 1204]
735
[1033, 1184]
158
456
[1161, 1312]
158
433
[1033, 1184]
158
433
1,033
158
433
158 1,033
433
1,033
158
433
[908, 1059]
158
558
[908, 1059]
158
558
158 [908, 1059]
558
[908, 1059]
158
558
[908, 1059]
158
558
[908, 1059]
158
558
[908, 1059]
158
558
158 [908, 1059]
558
[908, 1059]
158
558
158
446
0
0
446
0 7
302
0
0
0
151
0
0
0
0
0
0
0
0
0 0
0
0
0
0
0
Excess waste flow
Allowable waste flow
Optimized waste flow
Allowable waste flow
Excess waste flow
k = 1 (tonne/week)
k = 0.5 (tonne/week)
316
1,004
[908, 1059]
158
1,004
158 [915, 1066]
860
[908, 1059]
158
558
[1059, 1210]
158
558
[908, 1059]
158
558
[908, 1059]
158
558
158 [908, 1059]
558
[908, 1059]
158
558
[908, 1059]
Optimized waste flow
165
1,004
[908, 1059]
165
1,004
165 [908, 1059]
1,004
[908, 1059]
165
1,004
[908, 1059]
165
1,004
[908, 1059]
165
1,004
[908, 1059]
165
1,004
165 [908, 1059]
1,004
[908, 1059]
165
1,004
[908, 1059]
Allowable waste flow
151
0
0
0
0
0 0
0
0
0
0
0
0
0
0
0
0
0
0
0
0 0
0
0
0
0
0
Excess waste flow
k = 5 (tonne/week)
316
1,004
[908, 1059]
165
1,004
165 [908, 1059]
1,004
[908, 1059]
165
1,004
[908, 1059]
165
1,004
[908, 1059]
165
1,004
[908, 1059]
165
1,004
165 [908, 1059]
1,004
[908, 1059]
165
1,004
[908, 1059]
Optimized waste flow
Stoch Environ Res Risk Assess (2014) 28:167–187 179
123
Fig. 2 Optimal waste-flow allocation pattern for a landfill, b composting facility and c recycling facility
Stoch Environ Res Risk Assess (2014) 28:167–187
Optimized waste flow (tonne/week)
180 1800
β = 0.8 Lower
(a)
1600
β = 0.99 Lower
β = 0.8 Upper
β = 0.9 Upper
β = 0.99 Upper
1400 1200 1000 800
λ = 0.5 λ = 1
λ = 5 λ = 10 λ = 0.5 λ = 1
Period 1
Optimized waste flow (tonne/week)
β = 0.9 Lower
1200
(b)
β = 0.8 Lower
1000 800
λ = 5 λ = 10 λ = 0.5 λ = 1
Period 2
λ = 5 λ = 10
Period 3
β = 0.9 Lower
β = 0.99 Lower
β = 0.8 Upper
β = 0.9 Upper
β = 0.99 Upper
600 400 200 0
λ = 0.5 λ = 1
λ = 5 λ = 10 λ = 0.5 λ = 1
Period 1
λ = 5 λ = 10 λ = 0.5 λ = 1
Period 2
λ = 5 λ = 10
Period 3
Optimized waste flow (tonne/week)
380
(c)
340
β = 0.8 Lower
300
β = 0.99 Lower
β = 0.8 Upper
260
β = 0.9 Upper
β = 0.99 Upper
β = 0.9 Lower
220 180 140 100 60
λ = 0.5 λ = 1
λ = 5 λ = 10 λ = 0.5 λ = 1
Period 1
110.4] 9 106. In addition, CVaR decreases as k increases. For example, when b is fixed at 0.99 and k changes from 0.5 to 10, the value of CVaR would decrease from $[90.8, 110.4] 9 106 to $[80.9, 100.6] 9 106. This is because of the changing trade-off between the expectation total cost and the CVaR criterion. Figure 6 illustrates the total expected cost under different k and b levels. It is indicated that the total expected cost would increase in responding to the lager risk parameters k and/or b. This implies that a higher system cost would guarantee a lower system risk. Conversely, if
123
λ = 5 λ = 10 λ = 0.5 λ = 1
Period 2
λ = 5 λ = 10
Period 3
the decision maker aims towards an economic plan, a higher expected loss could be confronted. In the real-world application, decision makers may need to choose between a more risky solution with a lower system cost and a more conservative solution with a higher system cost. Figure 7 presents the recourse cost under different k and b levels. According to these results, increasing k leads to a more risk-averse policy with higher total expected costs and lower recourse costs in general. Therefore, a more conservative (i.e., averting risks) policy would lead to a lower excessive waste flow and a lower recourse cost.
Excess waste flow (tonne)
Stoch Environ Res Risk Assess (2014) 28:167–187
181
100000 Lower
80000
Upper
60000
(a) 40000 20000 0
β= 0.8
β= β= β= 0.9 0.99 0.8 λ = 0.5
β= β= β= 0.9 0.99 0.8 λ=1
β= β= β= 0.9 0.99 0.8 λ=5
β= β= 0.9 0.99 λ = 10
Risk parameter
Excess waste flow (tonne)
16000
(b)
Lower
12000
Upper
8000
4000
0 β= 0.8
β= β= β= 0.9 0.99 0.8 λ = 0.5
β= β= β= 0.9 0.99 0.8 λ=1
β= β= β= 0.9 0.99 0.8 λ=5
β= β= 0.9 0.99 λ = 10
Risk parameter Excess waste flow (tonne)
600 Lower
(c)
Upper
500 400 300 200 100 0 β= 0.8
β= β= 0.9 0.99 λ = 0.5
β= 0.8
β= β= 0.9 0.99 λ=1
β= 0.8
β= β= 0.9 0.99 λ=5
β= 0.8
β= β= 0.9 0.99 λ = 10
Risk parameter
Fig. 3 Excess waste flows to a landfill, b composting facility and c recycling facility for different risk parameters
5 Discussion In order to better reflect the effect of risk parameters k and b on the optimal results, sensitivity analysis is conducted in this research. Figure 8 shows the total expected cost under different k and b levels, where the upper and lower parts correspond to the upper and lower bounds of total expected cost, respectively. It appears that different combination of k and b levels would notably influence the value of total expected cost. For example, the total expected cost would be $[62.1, 77.5] 9 106 when b = 0.5 and k = 0.1, $[66.9, 82.8] 9 106 when b = 0.9 and k = 3, and $[77.1, 95.0] 9 106 when b = 0.99 and k = 10. A higher
confidence level (b) and a higher acceptable risk level (k) would give rise to a higher total expected cost. Conversely, a lower confidence level and a lower acceptable risk level would correspond to a lower total expected cost. Figure 9 presents the distribution of CVaR under different k and b levels. The specified b level of CVaRb quantifies the mean value of the worst (1 - b) % of the total costs. When b increases, CVaRb accounts for the risk of larger realizations. Thus, larger b values would lead to more conservative policies, which give higher weights to worse scenarios. The specified k level represents the risk preference of decision makers. When k = 0, decision makers would not consider the variability of the uncertain recourse costs based on a risk-neutral attitude. When k = 100, decision makers would seriously consider the variability based on a risk-aversive attitude. Moreover, increasing the value of k would increase the relative importance of the risk related terms and thus would lead to risk-averting policies. Sensitivity analysis results of Figs. 8 and 9 proclaim the trade-offs between costs and risk levels. A plan with low levels of the parameter k and/or the parameter b would correspond to lower system cost, implying that the decision makers have an optimistic attitude; however, it might be associated with a higher risk levels (i.e., excess of waste generation rate). Conversely, a plan with high levels of the parameter k and/or the parameter b would better resist from waste excess, and it would result in a lower risk of system failure. Therefore, the IMTSP encourages decision makers to assign the b and k to adjust risk control levels based on their preferences, and ensure management policies be made with reasonable consideration of both system stability and cost. Also, Figs. 8 and 9 indicate another fact that b and k have different sensitivity levels. When k is less than 0.5 and/or b is less than 0.7, the variation degree of total expected cost would increase gradually and the variation degree of CVaR decrease gradually, with an increase of k value. This implies that the risk constraint in effect is poor. When k changes from 0.5 to 3 and b is larger than 0.7, the variation degree of total expected cost and CVaR are both sharp. This implies that the risk constraint in effect is robust. When k is larger than 5, the values of total expected cost and CVaR under different b would be nearly constant, which indicates that the risk constraint is of significant effect. Thus, the k and b values in the CVaR constraint should be properly chosen to avoid the invalidation of risk control. The decision makers should incorporate their implicit knowledge (such as socio-economic conditions) and preferences about the risk and costs into the WSM management problems for generating more practical decisions. Uncertainty sources are important for MSW management systems. The random characteristics of various processes and conditions, the errors in acquiring the modeling
123
182
Stoch Environ Res Risk Assess (2014) 28:167–187
Table 6 Solutions of IMTSP method for binary variables under different risk parameters Facility
Total expanded capacity (tonne/week)
Composting
Period
b = 0.9
k = 0.5
k=1
k=5
k = 0.5
k=1
k=5
k = 0.5
k=1
k=5
132
k=1
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
270
k=1
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
483 189
338 132
k=1 k=2
[1, 1] [0, 0]
[1, 1] [0, 0]
[1, 1] [0, 0]
[1, 1] [0, 0]
[1, 1] [0, 0]
[1, 1] [0, 0]
[1, 1] [0, 0]
[1, 1] [0, 0]
[1, 1] [0, 0]
385
270
k=2
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
483
338
k=2
[1, 1]
[1, 1]
[1, 1]
[1, 1]
[1, 1]
[1, 1]
[1, 1]
[1, 1]
[1, 1]
189
132
k=3
[1, 1]
[1, 1]
[1, 1]
[1, 1]
[1, 1]
[1, 1]
[1, 1]
[0, 0]
[0, 0]
385
270
k=3
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[1, 1]
[1, 1]
483
338
k=3
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
350
140
k=1
[1, 1]
[1, 1]
[1, 1]
[1, 1]
[1, 1]
[1, 1]
[1, 1]
[0, 0]
[0, 0]
700
280
k=1
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
875
350
k=1
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
350
140
k=2
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[1, 1]
[1, 1]
700
280
k=2
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
875
350
k=2
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
350
140
k=3
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
700
280
k=3
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
875
350
k=3
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
[0, 0]
1102 782
747
597 481
425
187
171
133
124
165
405
400
534
510
600
CVaR ($ 10 6 )
90.8 886
Upper
800
100
983
941
Lower
96
99
107
131
137
152
β= 0.9
β= 0.99
β= 0.8
β= 0.9
β= 0.99
λ= 1
80
75.1
69.7
β= 0.8
λ = 0.5
90
89.3
(a) Lower bound 81.2
80.9
71.4
71.4
68.6
68.5
68.2
λ= 1
λ= 5
λ = 10
74.0
70
60
0 β= 0.8
β= 0.9 λ= 5
β= 0.99
β= 0.8
β= 0.9
β= 0.99
λ = 0.5
Risk parameter
λ = 10
Risk parameter
400 350
Fig. 4 Mean-risk function value of f for different risk parameters
(b) Upper bound
CVaR ($ 10 6 )
parameters, and the imprecision of the related system constraints are all possible sources of the uncertainties (Li et al. 2009). Waste generation rate arises as a direct consequence of human activities and consumption level of individual residents. However, the total population and consumption level may be of the stochastic characteristic in the planning horizon. Accordingly, waste generation rate varies from \1,000 tonne/week (i.e., the extra low scenario) to [2,000 tonne/week (i.e., the extra high scenario), which could be expressed as a random variable. Thus, waste generation rate could be regarded as the stochastic sources for two stage stochastic model. Moreover, the object of this research was to develop an inexact
123
b = 0.99
189
1000
200
b = 0.8
385
1200
122
Mean-risk function value ($ 10 6 )
Recycling
Expanded capacity for residential (tonne/week)
300 250 200
110.4
108.5
100.6
100.6
94.0
92.5
90.0
89.6
88.4
86.9
86.0
85.9
λ = 0.5
λ= 1
λ= 5
λ = 10
150 100 50
Risk parameter
Fig. 5 Value of CVaRb for different risk parameters
optimization model for allocating the waste flow generated by city residents. Thus, waste generation rate was one of the most important parameters to be considered for the
Stoch Environ Res Risk Assess (2014) 28:167–187
183
100
Total expected cost ($ 10 6 )
β = 0.8 β = 0.9
90
β = 0.99
80
70
60
λ = 0.5
λ=1
λ=5
λ = 10
λ = 0.5
λ=1
Lower
λ=5
λ = 10
Upper
Fig. 6 Total expected cost for different risk parameters
Recourse cost ($ 10 6 )
7 6
β = 0.8
5
β = 0.9 β = 0.99
4 3 2 1 0
λ = 0.5
λ=1
λ=5
λ = 10
λ = 0.5
λ=1
Lower
λ=5
λ = 10
Upper
Fig. 7 Recourse cost for different risk parameters
developed IMTSP model. On the other hand, the variability of random waste generation rate may result in the risk of economic loss from extreme events. If the extra high scenario of waste generation rate happens, waste disposal facilities might not have capacity to deal with the extra solid waste. The excess waste would become dangerous for city environment and human health. Thus, the decision makers have to take expensive measures, such as transporting it to neighboring city, to deal with the excess waste. CVaR as a kind of risk measure method is incorporated within the objective function of IMTSP method to describe the expected losses on extreme conditions. Therefore, the IMTSP can help decision makers analyze the trade-off between system cost and risk. To better reflect advantages of the proposed IMTSP method, ITSP is applied to the study case for comparison. The total expected cost of ITSP model would be $[61.34, 77.42] 9 106, which would be lower than the one of IMTSP. Solutions of binary variables are the same as for ITSP and IMTSP methods (k and b are fixed at 0.5 and 0.8, respectively). Comparatively, solutions of continuous variables are different for ITSP and IMTSP. Table 7 presents solutions of continuous variables for ITSP. Different from IMTSP, ITSP searches for the optimal solutions to
obtain the minimum cost without any consideration of risks. Consequently, it is incapable of analyzing trade-offs between system costs and risks (Maqsood et al. 2005). A lower level of allowable waste flow (i.e. X
ik;opt ) would be generated by ITSP to obtain minimum costs, leading to a higher level of excessive waste flow (i.e. Y
isk;opt ). Moreover, even though a lower cost would be obtained by ITSP compared with IMTSP, the obtained solutions would be limited in a number of aspects. First, for a basic municipal facility, such as a landfill, a fixed amount of waste flow normally should be guaranteed to ensure its good running; an overly high excess of waste flow, as suggested by the ITSP model, would cause serious consequences to such a facility. Second, the solutions of ITSP may lead to unbalanced allocation patterns among different MSW management facilities; this could provoke serious contradictions and harm implementation of waste management policies. Third, although the total expected cost is minimum, the large deficit amount at a high waste flow condition (i.e., Y
391;opt ) would occur at a probability level of ps = 0.004. If this really happens, a high second-stage penalty would have to be paid and high system infeasibility may need to be faced. On the other hand, solutions of ITSP seem to be exactly the same as those of IMTSP when the value of k is zero, which indicates that ITSP is a special or simplified case of IMTSP. Without consideration of interval mathematical programming, the MTSP method can be used to solve the study case. Its deterministic parameters can be derived by averaging the lower and upper bounds of intervals, such as the capacity of MSW management facilities. Accordingly, the solutions are fixed values rather than discrete intervals although their varying trends are similar to the solutions of IMTSP. Moreover, the solutions of MTSP are special cases in the solutions obtained from IMTSP. In such a case, the decision alternative would be restricted to a single solution. Meanwhile, the effectiveness and flexibility of the alternative would be reduced since impact of uncertainties is not considered. Robust risk analysis method (RRAM) was introduced by Chen et al. (2012), which incorporated interval mathematical programming, two-stage stochastic programming and variance within a general modeling framework. Thus, the main difference between IMTSP and RRAM is the quantitative risk measure for the secondstage random variable, where the former uses the CVaR and the latter uses the variance. Moreover, compared with IMTSP, the disadvantage of RRAM is that it considers the under-and-over-performances of the probability distribution equally. For example, the high waste generation rate may result in system risks, while the low waste generation rate is safe to the management system. Therefore, quantitatively evaluating the right tail of distributions of waste
123
184
(a) Lower bound
80
Total expected cost ($ 10 6 )
Fig. 8 Distribution of the lower and upper bounds of total expected cost under different k and b levels
Stoch Environ Res Risk Assess (2014) 28:167–187
75
70
65
60
λ = 0.01
λ = 0.1
β = 0.5
λ = 0.2
λ = 0.3
β = 0.6
λ = 0.5 β = 0.7
Total expected cost ($ 10 6 )
λ=2
β = 0.8
λ=3
λ=5
β = 0.9
λ = 10 β = 0.95
λ = 50
λ = 100
β = 0.99
(b) Upper bound
100
95
90
85
80
75
λ = 0.01
λ = 0.1
λ = 0.2
generation rate is necessary; however consideration the left tail of distributions must not only waste the expected cost but also be of no use for the system stability. One of the difficulties in planning a MSW management system is that those interval and stochastic uncertainties may be further multiplied by their associations with economic penalties if the promised targets are violated. Fortunately, two-stage stochastic programming and multi-stage stochastic programming can both effectively deal with such a problem. In this research, a time interval of 15 years is considered, which is divided into three 5-year periods. This implies that multistage stochastic programming is a suitable method, because it extends the two-stage stochastic programming by permitting modified decisions in each time stage based on the information of sequentially realized uncertain events (Li et al. 2009). Moreover, scenario-tree based multistage stochastic programming can not only reflect the correlations between the waste generation rates of the different time periods, but also model the decision process and define all possible scenarios. However, the main objective of this study is to control risks associated with the recourse actions and analyze trade-offs between expected costs and
123
λ=1
λ = 0.3
λ = 0.5
λ=1
λ=2
λ=3
λ=5
λ = 10
λ = 50
λ = 100
system risks. This objective can be achieved by the proposed two-stage stochastic programming associated with CVaR. Moreover, more scenarios are necessary for CVaR to quantify the risks from the extreme high waste generation rate. If a multi-stage stochastic programming is taken place the two-stage stochastic programming, the complicated scenario-tree would result in a concern upon the huge computational efforts. In other hand, to merge above advantage of multi-stage stochastic programming, future studies would be undertaken to investigate the performances of a hybrid model by integrating the multi-stage stochastic programming, CVaR and interval mathematical programming within a general framework.
6 Conclusions In this study, an interval-parameter mean-CVaR two-stage stochastic programming (IMTSP) method was developed for MSW management under uncertainty. This method incorporated interval mathematical programming, two-stage stochastic programming and CVaR within a general
Stoch Environ Res Risk Assess (2014) 28:167–187
185
(a) Lower bound
100
CVaR ($ 10 6 )
90
80
70
60 λ = 0.01
λ = 0.1
λ = 0.2
β = 0.5
λ = 0.3
β = 0.6
λ = 0.5 β = 0.7
λ=1
λ=2
β = 0.8
λ=3
λ=5
β = 0.9
λ = 10
λ = 50
β = 0.95
λ = 100
β = 0.99
(b) Upper bound
120
CVaR ($ 10 6 )
110
100
90
80
λ = 0.01
λ = 0.1
λ = 0.2
β = 0.5
λ = 0.3
β = 0.6
λ = 0.5 β = 0.7
λ=1
λ=2
β = 0.8
λ=3
λ=5
β = 0.9
λ = 10 β = 0.95
λ = 50
λ = 100 β = 0.99
Fig. 9 Distribution of the lower and upper bounds of the CVaR under different k and b levels Table 7 Solution of the ITSP model for continuous variables Period
To Landfill (i = 1) k=1
To composting (i = 2) k=2
To recycling (i = 3)
k=3
k=1
k=2
k=3
k=1
k=2
k=3
[881, 1032]
198
367
433
136
158
158
Allowable waste flow (X
ik ) (tonne/week) [1046, 1187]
[900, 1046]
Excess waste flow (Y
isk ) (tonne/week) s=1
0
0
0
0
0
0
0
0
0
s=2
0
0
0
0
0
0
0
0
0
s=3
0
0
0
0
0
0
0
0
0
s=4
[0, 1]
[0, 1]
[0, 1]
0
0
0
0
0
0
s=5
142
147
152
0
0
0
0
0
0
s=6
[283, 284]
293
280
0
0
23
0
0
0
s=7
[425, 426]
[372, 373]
21
0
67
433
0
0
0
s=8
388
219
172
179
367
433
0
0
0
s=9
[374, 375]
207
165
198
367
433
136
158
158
framework. In the developed IMTSP, uncertainties could be presented as probabilistic distribution functions and discrete intervals. Also, CVaR as a kind of risk measure method was
successfully incorporated within the objective function of IMTSP to describe the expected losses under extreme conditions. In general, the proposed IMTSP method has several
123
186
advantages: (i) its objective function simultaneously takes expected cost and system risk into consideration. Both discrete random variables and intervals can be used to reflect uncertainties. Therefore, the solution results would possess characteristics of CVaR, interval mathematical programming, and two-stage stochastic programming. The MTSP and ITSP model do not have such an advantage; (ii) it can be used for supporting quantitatively evaluating the right tail of distributions of waste generation rate, which can better quantify the system risk compared to RRAM method; (iii) it is useful for helping decision makers analyze the trade-offs between cost and risk and identify desired waste management strategies under complex uncertainties; and (iv) it is effective to penalize the second-stage costs that are above the expected values, as well as to capture the notion of risk in stochastic programming. The developed IMTSP was applied to the long-term planning of MSW management system in the City of Regina, Canada. The initial allocation target value could be used to reflect decision maker’s opinion on waste allocation, and the risk acceptance level (i.e. k) and the percentile level (i.e., b) could reflect decision maker’s preference upon system cost and risk aversion. Moreover, in the case of excess waste, the allotment to the landfill should be assigned initially, and then to the incinerating and composting facilities. The study results demonstrated that IMTSP method could generate a series of waste management patterns under different risk-aversion levels, gain indepth insights into the effects of uncertainties, and maintain a balance between economic cost and environmental protection. Thus, the IMTSP method was valuable for supporting: (i) the adjustment or justification of the existing waste flow allocation patterns, (ii) the long-term capacity planning for waste-management facilities, and (iii) the strong bases for selecting desired waste management policies. In choosing appropriate solutions for practical applications, the decision makers should incorporate their implicit knowledge (such as socio-economic conditions) and preferences about the risk and costs into the WSM management problems for generating more practical decisions. Although this study is the first attempt for the planning of waste management systems through the developed IMTSP method, the results suggested that it would be an effective tool for decision makers in conducting long-term MSW planning. However, the IMTSP also has much space for improvement. Multistage stochastic programming is a more suitable method for the multi-period planning, and it can reflect the correlations between the waste generation rates of the different time periods. Thus, it is an interesting work to integrate the CVaR, interval mathematical programming and multi-stage stochastic programming instead of two-stage stochastic programming. Moreover, discrete
123
Stoch Environ Res Risk Assess (2014) 28:167–187
distribution values of stochastic variables (i.e., waste generation rates) are different to estimate. The multi-uncertainty (i.e., probability intervals) can used to define the parameter of waste generation rates. Thus, many other uncertainty methods which can used to solve the multiuncertainty have potentials to be further integrated into an IMTSP framework for reflecting more complex conditions. Acknowledgments This research was supported by the National key basic research development planning project (No. 2010CB 428501) and by the National high technology research and development program (Nos. 2008AA06A415 and 2009AA06A41802). The authors are grateful to the editors and the anonymous reviewers for their insightful comments and suggestions.
References Ahmed S (2004) Mean-risk objectives in stochastic programming. http://www2.isye.gatech.edu/*sahmed/sprisk.pdf. Accessed 20 June 2012 Ahmed S (2006) Convexity and decomposition of mean-risk stochastic programs. Math Program 106(3):433–446 Ahmed S, Tawarmalani M, Sahinidis N (2004) A finite branch-andbound algorithm for two-stage stochastic integer programs. Math Program 100(2):355–377 Birbil SI, Frenk J, Kaynar B, Noyan N (2008) Risk measures and their applications in asset management. In: Gregoriou GN (ed) The VaR implementation handbook. The McGraw-Hill Companies, New York, pp 311–337 Birge JR, Louveaux FV (1988) A multicut algorithm for two-stage stochastic linear programs. Eur J Oper Res 34(3):384–392 Cai YP, Huang GH, Tan Q, Chen B (2011) Identification of optimal strategies for improving eco-resilience to floods in ecologically vulnerable regions of a wetland. Ecol Model 222(2):360–369 Chen C, Huang GH, Li YP, Zhou Y (2012) A robust risk analysis method for water resources allocation under uncertainty. Stoch Environ Res Risk. doi:10.1007/s00477-012-0634-5 Cheng GH, Huang GH, Li YP, Cao MF, Fan YR (2009) Planning of municipal solid waste management systems under dual uncertainties: a hybrid interval stochastic programming approach. Stoch Environ Res Risk Assess 23(6):707–720 City of Regina [Internet] (2011) Waste plan Regina report. http:// www.regina.ca/AssetFactory.aspx?did=3119. Accessed 3 Sept 2011 Dai C, Li Y, Huang GH (2011) A two-stage support-vector-regression optimization model for municipal solid waste management: a case study of Beijing, China. J Environ Manag 92(12):3023–3037 Dai C, Li Y, Huang GH (2012) An interval-parameter chanceconstrained dynamic programming approach for capacity planning under uncertainty. Resour Conserv Recycl 62:37–50 Fabian CI (2008) Handling CVaR objectives and constraints in twostage stochastic models. Eur J Oper Res 191(3):888–911 Guo P, Huang GH (2009a) Two-stage fuzzy chance-constrained programming: application to water resources management under dual uncertainties. Stoch Environ Res Risk Assess 23(3): 349–359 Guo P, Huang GH (2009b) Inexact fuzzy-stochastic mixed integer programming approach for long-term planning of waste management-part B: case study. J Environ Manag 91(2):441–460 Hsu C, Huang C, Paul Chiou W (2011) Effectiveness of copula-extreme value theory in estimating value-at-risk: empirical evidence from Asian emerging markets. Rev Quant Financ Acc. doi:10.1007/ s11156-011-0261-0
Stoch Environ Res Risk Assess (2014) 28:167–187 Huang GH, Loucks D (2000) An inexact two-stage stochastic programming model for water resources management under uncertainty. Civil Eng Environ Syst 17(2):95–118 Huang GH, Brian WB, Gilles GP (1992) A grey linear programming approach for municipal solid waste management planning under uncertainty. Civil Eng Syst 9(4):319–335 Kall P, Mayer J (2005) Stochastic linear programming: models, theory, and computation. In: International series in operations research and management science. Springer, New York Kollikkathara N, Feng H, Stern E (2009) A purview of waste management evolution: special emphasis on USA. Waste Manag (Oxf) 29(2):974–985 Li YP, Huang GH (2009) Inexact minimax regret integer programming for long-term planning of municipal solid waste management-part B: application. Environ Eng Sci 26(1):219–234 Li YP, Huang GH, Chen X (2009) Multistage scenario-based intervalstochastic programming for planning water resources allocation. Stoch Environ Res Risk 23(6):781–792 Maqsood I, Huang GH (2003) A two-stage interval-stochastic programming model for waste management under uncertainty. J Air Waste Manag Assoc 53(5):540–552 Maqsood I, Huang GH, Huang YF, Chen B (2005) ITOM: an intervalparameter two-stage optimization model for stochastic planning of water resources systems. Stoch Environ Res Risk 19:125–133 Markowitz H (1952) Portfolio selection. J Financ 7(1):77–91 Noyan N (2012) Risk-averse two-stage stochastic programming with an application to disaster management. Comput Oper Res 39(3): 541–559 Qin XS, Huang GH, Zeng G, Chakma A, Huang Y (2007) An interval-parameter fuzzy nonlinear optimization model for stream water quality management under uncertainty. Eur J Oper Res 180(3):1331–1357
187 Rockafellar RT, Uryasev S (2000) Optimization of conditional valueat-risk. J Risk 2:21–42 Rockafellar RT, Uryasev S (2002) Conditional value-at-risk for general loss distributions. J Bank Financ 26(7):1443–1471 Schultz R (2011) Risk aversion in two-stage stochastic integer programming. Stoch Program 150:165–187 Schultz R, Tiedemann S (2006) Conditional value-at-risk in stochastic programs with mixed-integer recourse. Math Program 105(2–3): 365–386 Seifi A, Hipel K (2001) Interior-point method for reservoir operation with stochastic inflows. J Water Res Plan Manag 127(1):48–57 Shao LG, Qin XS, Xu Y (2011) A conditional value-at-risk based inexact water allocation model. Water Resour Manag 25(9):2125–2145 Tan Q, Huang GH, Cai YP (2010a) Waste management with recourse: an inexact dynamic programming model containing fuzzy boundary intervals in objectives and constraints. J Environ Manag 91(9): 1898–1913 Tan Q, Huang GH, Cai YP (2010b) A superiority-inferiority-based inexact fuzzy stochastic programming approach for solid waste management under uncertainty. Environ Model Assess 15(5): 381–396 Tong XJ, Qi LQ, Wu F, Zhou H (2010) A smoothing method for solving portfolio optimization with CVaR and applications in allocation of generation asset. Appl Math Comput 216(6):1723–1740 Xu Y, Huang GH, Qin XS, Huang Y (2009) SRFILP: a stochastic robust fuzzy interval linear programming model for municipal solid waste management under uncertainty. J Environ Informatics 14(2):72–82 Zhu H, Huang GH (2011) SLFP: a stochastic linear fractional programming approach for sustainable waste management. Waste Manag (Oxf) 31(12):2612–2619
123