A Generalized Fuzzy Integer Programming Approach for Environmental Management under Uncertainty

1 Faculty of Engineering and Applied Science, University of Regina, Regina, SK, Canada S4S 0A2 2MOE Key Laboratory of Regional Energy Systems Optimization, Resources and Environmental Research Academy, North China Electric Power University, Beijing 102206, China 3 College of Environmental Science and Engineering, Xiamen University of Technology, Xiamen, Fujian 361024, China 4College of Urban Construction, Hebei University of Engineering, Handan, Hebei 056038, China


Introduction
Municipal solid waste (MSW) management is a priority for many developed and developing countries throughout the world.Effective planning of MSW is critical for supporting sustainable socioeconomic development in urban communities.However, extensive uncertainties may exist in many system components and impact factors.For example, waste generation rate within a city is related to many socioeconomic and environmental factors and exhibits various uncertain features.Such uncertainties and their interactions can lead to increased complexities in the related planning efforts and will affect consequent decision processes.Besides, these uncertainties may be further multiplied because many system components are of multiperiod, multilayer, and multiobjective features Li and Huang [1].Moreover, waste-disposal facilities in a MSW management system usually have overall cumulative or daily operating-capacity limits.Increasing waste generation rates, as a result of population explosion and economic development, lead to intensified conflicts with decreasing waste-treatment/disposal capacities.Therefore, it is desired that the above uncertain and dynamic complexities be reflected in efforts for identifying effective environmental management alternatives.
In the past decades, a number of inexact optimization techniques were developed to deal with uncertainties and dynamics in MSW management.They were mainly classified as fuzzy, stochastic, and interval mathematical programming (FMP, SMP, and IMP, resp.)[2][3][4][5][6].For example, Li and Huang [7] proposed an inexact two-stage mixed-integer linear programming (ITSMILP) method for solid waste management in the city of Regina, through incorporating interval linear programming (ILP), two-stage stochastic programming (TSP), and mixed-integer programming (MIP) within a general mathematical programming framework.However, IMP could merely deal with interval uncertainties without distributional information; SMP was inapplicable to large-scale problems due to its stringent requirements for information of probabilistic distributions.
Fuzzy mathematical programming (FMP), as a branch of fuzzy set theory, could generally deal with uncertainties expressed as fuzzy sets or fuzzy goals/constraints [8][9][10][11].Recently, various FMP methods were employed to deal with uncertainties in MSW systems [12][13][14][15].For example, Fan et al. [12] explored a fuzzy linear programming (FLP) method for dealing with uncertainties expressed as fuzzy sets that exist in the constraints' left-hand and right-hand sides and the objective function; however, this method was unable to reflect dynamic complexities related to capacity expansion schemes for waste-treatment facilities.Srivastava and Nema [13] proposed a fuzzy parametric programming model for identifying desired treatment/disposal facilities, planning waste management capacities, and allocating waste flows under uncertainty; however, the proposed method generated deterministic waste allocation schemes without provision of bases for supporting generation of multiple decision options corresponding to the uncertain system conditions.
Generalized fuzzy linear programming (GFLP) (or fully fuzzy linear programming (FFLP)) methods extended traditional FMP approaches through permitting uncertain information in the optimization process and resulting solutions.Recently, several GFLP (or FFLP) methods were proposed to deal with uncertain information in both parameters and decision variables [16][17][18][19][20][21].For example, Hosseinzadeh Lotfi et al. [18] developed a lexicography method to solve the FFLP problem and generate approximate solutions presented as fuzzy sets.However, previous studies on fuzzy variables in FMP problems mainly focused on some special types of fuzzy sets (such as symmetric, triangular, or trapezoidal fuzzy sets).Furthermore, some of them might lead to complicated intermediate models and thus were not applicable for large-scale problems.Fan et al. [19,20] proposed another kind of fuzzy programming (named generalized fuzzy linear programming (GFLP) method) to deal with fuzzy uncertainty in both parameters and variables, in which all fuzzy sets with known membership functions can be treated through defuzzification method (i.e., -cut method).However, the GFLP approach cannot reflect dynamic features in environmental management problems.Moreover, no previous study was reported on capacity expansion issues under fuzzy conditions through the generalized fuzzy optimization approach, where expansion schemes were desired under multiple scenarios and -cut levels.
As an extension of developed GFLP approach, a generalized fuzzy integer programming (GFIP) method would be proposed for MSW management under uncertainty.The proposed GFIP approach integrates the techniques of generalized fuzzy linear programming (GFLP) and mixed-integer programming (MIP) within an optimization framework.In detail, (i) the GFIP method can deal with uncertainties expressed as fuzzy sets with known membership functions, regardless of whether these functions are linear or nonlinear; (ii) the proposed GFIP method can allow uncertainties to be directly communicated into the optimization process and the resulting solutions; (iii) the GFIP method can reflect dynamics in terms of waste-flow allocation and facilitycapacity expansion; (iv) compared with other inexact mixedinteger programming approaches (e.g., ITSMILP by Li and Huang [7]), the GFIP can analyze the inherent interrelationship between the uncertainty of fuzzy parameters (i.e., -cut levels) and capacity expansion options of waste management facilities, and such analysis can help decision makers make tradeoffs between system reliability and system cost.Then, a case study will be provided to demonstrate applicability of the GFIP method to support dynamic analysis for MSW management under uncertainty.The results will be used for generating different decision alternatives under various system conditions and thus for helping identify desired waste management policies.

Solution Method of GFIP Model through
Stepwise Interactive Algorithm.If the parameters and variables in model (1a)-(1e) are triangular fuzzy numbers, several methods can be applied to solve the model, such as the lexicography method proposed by Hosseinzadeh Lotfi et al. [18] and the methods of Fan et al. [12] and Kumar et al. [21].However, when the parameters of model (1a)-(1e) are expressed through other kinds of fuzzy numbers, the above methods are not applicable.Consequently, in this study, a new method named stepwise interactive algorithm (SIA) will be proposed to solve model (1a)-(1e).This algorithm is based on computational principles related to fuzzy intervals [25][26][27][28] (see Appendix section).The detailed proof of the solution algorithm can be found in Fan et al. [20].The inherent idea of the stepwise interactive algorithm is based on the design of experiment, in which the optimization model would be considered as an experiment with the -cut levels being the inputs and the optimal solutions being the outputs.The detailed procedures of the SIA method include (i) discretizing the membership function grade of fuzzy parameters into a set of -cut levels; (ii) converting the GFIP problem into an inexact mixedinteger linear programming (IMILP) problem under each cut level; (iii) solving the IMILP problem through an interactive algorithm; and (iv) approximating the membership function for decision variables through statistical regression methods.Compared with the previous methods, SIA can allow uncertainties to be directly communicated into the optimization process.Moreover, it will not lead to complex intermediate submodels and thus lead to a relatively low computational requirement.This is meaningful when the SIA method is applied to solve large-scale management models.Finally, the proposed SIA method can generate solutions expressed as fuzzy sets.
Since the parameters in model (1a)-(1e) are expressed as fuzzy sets, these parameters will be defuzzified before the model is solved.Various defuzzification methods have been proposed to convert fuzzy sets into crisp sets, including cut, max-membership principle, centroid, weighted average, mean-max membership, center of sums, center of largest, and first of maxima or last of maxima methods.In this study, the -cut would be applied to defuzzify the fuzzy parameters in model (1a)-(1e) due to its popularity and ease of implementation.The concept of -cut is important in reflecting the relationship between fuzzy sets and crisp sets.Each fuzzy set can be uniquely represented by all of its cuts.As stated by Kreinovich [29], fuzzy data processing is computable for -cuts but, in general, not computable for membership functions.Consequently, the fuzzy parameters and decision variables in model (1a)-(1e) are defuzzified through the -cut method instead of their membership functions.Through the -cut method, the fuzzy parameters and decision variables in model (1a)-(1e) will be converted into the related fuzzy intervals.The optimization model with interval parameters can then be transformed into deterministic submodels, which can be solved through ordinary solution methods (e.g., simplex method).Therefore, before solving model (1a)-(1e), a set of -cut levels (i.e.,  ] will be appointed firstly to cut model (1a)-(1e).Then an inexact mixed-integer linear programming (IMILP) model can be formulated as follows: ≥ 0  = 1, 2, . . ., = interval continuous variables, (  ) ∈ () ±  (1) , = 1, 2, . . .,  ( < ) , = interval integer variables, (  ) ∈ () ±  (1) , =  + 1,  + 1, . . ., , where , and (  ) ±  (1) are fuzzy intervals under  (1) .
]. Fuzzy intervals under other -cut levels also have similar expressions.Furthermore, an interval number ( ± ) can be defined as Model (2a)-(2e) shows the formulation of interval mixedinteger linear programming (IMILP) method with all parameters expressed as interval numbers.The IMILP model was developed through introducing the concept of interval analysis into a mixed-integer linear programming framework [3].It allowed uncertainties to be directly communicated into the optimization processes and resulting solutions and did not lead to complicated intermediate models [3].
Formulas ((3a)-(3e)) to ((5a)-(5b)) show the detailed solution process of an IMILP model through the interactive algorithm (also named two-step method).Based on the interactive algorithm, the original IMILP model is firstly reformulated into two submodels corresponding, respectively, to its upper and lower bounds of objective function; the two submodels are then solved separately one after another [30].The sequence to solve two submodels is subject to the nature of objective function (max or min).For a maximized problem [i.e., model (2a)-(2e)], the submodel corresponding to the upper bound of the objective function is solved first, followed by solving the submodel corresponding to the lower bound of the objective function; besides, the optimal solutions from the first submodel should be used as constraints for the second submodel [30].
Based on solutions of model (2a)-(2e), we will select  (2) to  () in sequence and then formulate corresponding IMILP models as follows: where are the optimal solutions obtained from the IMILP model under  (−1) .Formula (6e) is proposed to reflect the property of the fuzzy number that Based on the interactive algorithm, model (6a)-(6f) will be converted into two submodels as follows.
Based on formulas ((2a)-(2e))-((9a)-(9b)), we can obtain a series of fuzzy interval solutions for model (1a)-(1e) under different -cut levels.Then we can approximate the membership function for continuous decision variables by statistical regression methods.In this procedure, the GFIP model is supposed to be an experiment, with -cut levels being its inputs (i.e., independent variables) and the lower and upper bounds of decision variables being its outputs (i.e., dependent variables).Take ( opt ) −  as an example, we can obtain a regression function between ( opt ) −  and  based on the fuzzy interval solutions.Such a regression function will be considered as the inverse function of the left shape function for x , denoted as  −1 (); then we can acquire the left shape function for x , expressed as ().In the same way, we can obtain the right shape function for x , expressed as ().

Case Study
A hypothetical municipal solid waste (MSW) management problem is used to illustrate the applicability of GFIP approach.The studied system includes three municipal cities.A planning horizon of 15 years is divided into three periods, with each one having a time interval of 5 years.Two types of facilities can be available for waste treatment/disposal.A landfill is considered in the proposed case due to its crucial role for MSW disposal in both developed and developing countries.For example, more than 54 percent of MSW was landfilled in the United States during 2009 [31], while 89.3 percent of the generated MSW (74.04 million tonnes) was landfilled in China in 2002 [32].The landfill is typically characterized as an overall capacity limit.Also, a waste-toenergy (WTE) facility, which can effectively minimize land depletion caused by landfilling, is employed to serve wastedisposal needs.It is characterized as a daily capacity limit.
In fact, a MSW management system involves several processes with socioeconomic and environmental implications, such as waste generation, transportation treatment, and disposal [33].Extensive uncertainties usually exist in these processes due to impacts of the economic development, population growth, and human activities.Moreover, probabilistic methods are not applicable to quantify these uncertainties when data are insufficient.Consequently, adoption of fuzzy set theory would be a potential alternative, especially when uncertainties can be consciously assumed by decision makers or experts.Furthermore, uncertain inputs in the MSW management system would lead to variations in the resulting solutions.Therefore, the GFIP method will be desired to reflect uncertain and dynamic complexities in the MSW management system and generate solutions expressed as fuzzy sets.
Table 1 shows related waste generation levels and cost coefficients, including waste generation rates in three cities, operation costs of two facilities, and transportation costs for shipping waste flows.These parameters are estimated as triangular fuzzy numbers with known most possible values, as well as left and right spreads.Table 2 presents capacity expansion options and related costs for waste disposal facilities.The total capacity of landfill is (1.8 × 10 6 , 0.2 × 10 5 , 0.3 × 10 5 ) tonne, which means the most possible capacity of landfill is 1.8×10 6 tonne, and the lower and upper bound is 1.78 × 10 6 and 1.83 × 10 6 tonne, respectively.The daily capacity of WTE facility is (390, 20, and 20) tonne/day, which means the most possible capacity and lower and upper bound is 390, 370, and 410 tonne/day, respectively.The WTE facility will generate residues of about 30% (10% as its left and right spread) of the incoming waste stream.The revenue from the WTE facility is approximately $20/tonne, with its left and right spreads being $2/tonne.In this study, all parameters are assumed to be triangular fuzzy numbers.The triangular fuzzy numbers are considered in this study because (i) the triangular form is the simplest type of fuzzy numbers, (ii) many other types of fuzzy numbers can be estimated through the triangular fuzzy numbers, and (iii) triangular membership function can provide the most important information for a fuzzy set: lower-bound value, upper-bound value, and the most possible value [34].Also, other kinds of fuzzy numbers can be treated through the proposed GFIP approach if their membership functions are known.
The problem under consideration is how to effectively allocate waste flows and choose appropriate capacity expansion options of waste-disposal facilities under a number of environmental, economic, and treatment/disposal constraints in order to minimize the overall system cost.A GFIP model can thus be formulated to solve this problem.
In this study, decision variable x represents the amount of waste flow from city  to waste-treatment facility  in period .The objective is to minimize the system cost through effectively allocating waste flows from three cities to two disposal facilities and choosing appropriate waste-disposalfacility options for excessive waste-disposal requirements.The constraints involve relationships between decision variables and waste generation/management conditions.Thus, a GFIP model can be formulated as follows: (1) Landfill capacity constraint: (2) WTE facility-capacity constraints: (3) Waste disposal demand constraints: (4) Nonnegativity constraints: (5) Nonnegativity and binary constraints: (6) Landfill expansion constraint: (7) WTE facility expansion constraints: where FE is the residue flow from WTE to landfill (% of incoming mass to WTE facility); FLC  is the capital cost of landfill expansion in period  ($); FT  is the transportation cost of waste flow from WTE to landfill in period  ($/tonne); FTC  is the capital cost of expanding WTE by option  in period  ($); ÕP  is the operating cost of facility  in period  ($/tonne); RE  is the revenue from WTE in period  ($/tonne); TE is the maximum capacity of WTE facility (tonne/day); TL is the capacity of landfill (tonne); TR  is the transportation cost for waste flow from city  to facility  during period  ($/tonne); WG  is the waste generation rate in city  during period  (tonne/day); ΔTC is the total amount of expansion capacity for landfill (tonne); ΔTC  is the amount of the th type of expansion capacity for WTE (tone/day); x is the waste-flow rate from city  to facility  in period  (tonne/day),  = 1, 2; j = 1, 2, 3;  = 1, 2, 3; Ỹ is the binary decision variable for landfill expansion at the start of period ; Z is the binary decision variable for WTE facility with expansion option  at the start period of ;  is the index for facility ( = 1 for landfill and  = 2 for WTE facility);  is the index for three cities ( = 1, 2, 3);  is the index for time periods ( = 1, 2, 3).
In model (10a)-(10i), the objective (i.e., formula (10a)) is to minimize the total cost of facility expansion and wasteflow disposal over the entire planning horizon, which will cover expenses of handling waste flows, charges of expanding facilities, and revenues from the WTE facility.Constraint (10b) specifies that the total amount of waste allocated to the landfill must not exceed its existing and expanded capacities.In this study, one year is assumed to have 365 days, and there are 5 years in each period.Consequently, the coefficient in constraint (10b) would be 1825 (i.e., 365 × 5).Constraint (10c) means that the actual daily waste flows shipped to the WTE facility should not exceed its existing and expanded capacities.Constraint (10d) indicates that, for each city in each period, the waste flows transported to the landfill and WTE must be not less than its waste-disposal demand in this period.This constraint also assumes that all solid wastes have to be shipped to a disposal site within a certain period after its generation, and no mass loss is incurred in the transportation process.Constraints (10f)∼(10i) define the binary variables related to capacity expansion decisions; constraint (10h) denotes that the landfill can only be expanded once within the entire planning horizon, and constraint (10i) means that the WTE can be expanded once in each period.
Figure 1 shows the schematic of the GFIP model for MSW management.Obviously, the GFIP model is an integration of generalized fuzzy linear programming (GFLP) and mixedinteger linear programming (MILP) methods.Each method has a unique contribution in enhancing the capability of GFIP in dealing with uncertainties and dynamic features in solid waste management.For example, fuzzy parameters can be handled by GFLP, and waste management facility expansion can be addressed by MILP.Besides, a stepwise interactive algorithm (SIA) is proposed for solving the proposed GFIP model, which can permit uncertainty to be directly communicated into the optimization process and resulting solutions.Through SIA, the developed GFIP model will firstly be converted into several IMILP submodels and then be further transformed into mixed-integer linear programming (MILP) submodels.Consequently, the computational complexity of the GFIP would be reasonable.For example, if  -cut levels are identified in solving the GFIP model,  IMILP submodels will be firstly generated.According to interactive algorithm, each IMILP submodel can be further converted into two MILP submodels; thus, the GFIP model will finally result in 2 MILP submodels with deterministic parameters.

Result Analysis
In this study, a GFIP model is developed for supporting decision making in MSW management.A stepwise interactive algorithm (SIA) is proposed to solve the GFIP model.Based on SIA, six -cut levels (i.e., 0, 0.3, 0.5, 0.7, 0.85, and 1) would be considered.Under each -cut level, the fuzzy parameters presented in Tables 1 and 2 would be converted into corresponding fuzzy intervals, and model (10a)-(10i) would also be transformed into an inexact mixedinteger linear programming (IMILP) model.Tables 3 to 5 present waste-flow diversion schemes, capacity expansion options, and related system costs obtained through IMILP models under the selected -cut levels.The results indicate that waste-flow patterns and capacity expansion options would vary due to temporal and spatial variations in waste generation rates and waste management conditions.
In period 1, waste from city 1 would be initially shipped to WTE, with a flow amount of 195 tonne/day.For city 1, the WTE facility would be the first choice to serve its waste disposal requirement.Then the remaining waste would be  (q) , where allocated to the landfill, with a waste flow of (30, 25, and 25) tonne/day.The (30, 25, and 25) indicates a triangular fuzzy number, with 30, 5, and 55 as its most possible value and lower and upper bound, respectively.The fuzzy characteristic of the waste flow to the landfill indicates that the variation in the waste generation rate of city 1 would be handled through landfilling.The waste disposal scheme for city 2 is much different from that of city 1.The waste-treatment demand of city 2 can be satisfied through landfilling.Consequently, no waste would flow to WTE in this period.Conversely, all of the generated waste in city 3 would be shipped to WTE in spite of its variation in waste generation rate.The waste allocation schemes for three cities in period 2 would be similar to those in period 1.The majority of waste from city 1 would be allocated to WTE, with the residues being shipped to landfill.However, compared with the waste flows in period 1, more waste would be transported to WTE in period 2 as a result of temporally increasing waste generation rate.Moreover, the amount of waste allocated to WTE from city 1 would fluctuate within small intervals under low plausibilities ((225, 236.6) and (225, 229.9) tonne/day under  = 0 and 0.3, resp.).This is because the waste generation rate would vary within significant ranges under low plausibilities.Meanwhile, landfill would be the only choice to satisfy the waste-treatment demand of city 2, while all waste from city 3 would be shipped to WTE.
The waste-flow patterns would be changed significantly in period 3.All waste from the three cities would be delivered to the landfill due to its lower operation cost.In detail, waste flows shipped to the landfill from cities 1, 2, and 3 would be capacity of landfill achieves its lower bound and the waste generation rates reach their upper bounds).This means that the existing capacity of landfill would be sufficient to dispose of the solid waste under advantageous conditions (i.e., the capacity of landfill achieves its upper bound and the waste generation rates reach their lower bounds).However, when -cut level increases to 0.5, the landfill would require an expansion in period 2 under advantageous conditions due to the increase in the lower bounds of waste generation rates and the decrease in the upper bound of landfill capacity.Under  = 0.7, expansion option of the landfill would be similar to that under  = 0.3, except the option being applied in period 2. As shown in Table 4, when  ≥ 0.85, the landfill would be expanded in period 2 under both demanding and advantageous conditions.In brief, the results in Table 4 suggest that (i) expansion of the landfill in period 1 leads to low risk but high system cost and (ii) the landfill at least is expanded in period 2.
As shown in Table 5, the WTE facility would generally be expanded in periods 1 and 2. In period 1, option 1 (i.e., 150 tonne/day) would be applied for WTE, regardless of the impacts of uncertain inputs.However, in period 2, the capacity expanding options would be influenced significantly by system uncertainty.The WTE facility would be expanded with different options under different plausibilities in this period.When  = 0, option 1 is considered to be suitable for WTE under advantageous conditions, while option 3 is chosen under demanding conditions.As the -cut level varies between 0.3 and 0.7, options 1 and 2 are applicable for WTE under advantageous and demanding conditions.Furthermore, as shown in Table 5, the WTE facility would at least be expanded with option 2 in period 2.
Based on the waste-flow diversion schemes under selected -cut levels (as shown in Table 3), we can approximate their membership functions through statistical regression.Figure 4 shows the membership functions for the wasteflow schemes.It indicates that these membership functions can be well fitted based on the results in Table 3.However, for the variables indicating facility expansion options, it can hardly obtain their membership function since they are binary variables.But decision makers can still generate appropriate waste-flow diversion schemes and facility expansion options based on Figure 4 and Tables 4 and 5.For example, if a decision maker wants to identify the waste management policy under an -cut level of 0.6, the waste-flow patterns can be obtained through the membership functions in Figure 4 (e.g., ( 111 ) − 0.6 = (0.6 + 0.2)/0.04= 20, ( 111 ) + 0.6 = (2.2 − 0.6)/0.4= 40); the facility expansion options can be adopted based on results under  = 0.5 and 0.7 in Tables 4 and 5.
Table 2 also provides the total system costs (denoted as the objective function) obtained from the GFIP model under six -cut levels.The results suggest that different plausibilities of uncertain inputs lead to varied system costs.The lower bounds of the objective function correspond to advantageous conditions, while the upper bounds are associated with demanding conditions.When  = 0, the system cost would be $[6.09,8.02] × 10 9 .Under  = 0.3, the lower bound of system cost would be $6.18 × 10 9 and the upper bound would be $7.27 × 10 9 , leading to a fluctuating interval of $[6.18, 7.27] × 10 9 for system cost.However, the value of the objective function does not necessarily hold such a characteristic that as the -cut level increases, the lower bound increases and the upper bound decreases (as shown in Figure 5).For example, the lower bound of the objective function under  = 0.5 and 0.7 is $6.24 × 10 9 and $3.8 × 10 9 , respectively.This is due to the variation in capacityexpanding options under these two -cut levels.When  = 0.5, capacity-expanding for landfill is considered under both demanding and advantageous conditions; in contrast, the landfill would not be expanded under advantageous conditions when  = 0.7, which leads to a significant decrease of the objective function value.
Generally speaking, the conventional inexact mixedinteger linear programming (IMILP) method would consider uncertain parameters with known lower and upper bounds.For the IMILP method, it cannot consider any distributional information between the lower and upper bounds.In this study, the solutions of GFIP model under  = 0 are identical to the solutions obtained from the IMILP model.As shown in Tables 3 to 5, the solutions of IMILP model can only provide interval values for waste allocation schemes and neglect distributional information within these intervals.Particularly, when distributional information is available for uncertain inputs, the IMILP method can hardly reflect the relationships between uncertainty of the inputs and the resulting solutions.Conversely, the GFIP method can not only provide fluctuating intervals for waste-flow allocation schemes (i.e., lower and upper bounds) as well as corresponding capacity expanding options for waste-treatment facilities but also afford plausibilities for such waste management alternatives.

Conclusions
In this study, a generalized fuzzy integer programming (GFIP) method was developed for solid waste management under uncertainty.The developed GFIP could deal with uncertainties expressed as fuzzy sets that exist in the constraints' left-and right-sides and the objective function.A stepwise interactive algorithm (SIA) was proposed to solve the GFIP model and generate solutions expressed as fuzzy sets.The SIA firstly discretized the membership function grade into a set of -cut levels.Then the GFIP model was converted into a series of IMILP submodels.The interactive algorithm proposed by Huang et al. [3] was applied to solve the IMILP submodels and generate interval solutions under each -cut level.The membership functions for fuzzy continuous variables were finally obtained through statistical regression method based on those interval solutions.
The developed method was applied to a case of municipal solid waste management to illustrate the applicability of GFIP method.The solutions for binary variables, which were obtained as discrete fuzzy sets, provided different capacity expansion alternatives for waste-treatment facilities under different plausibilities of uncertain inputs.The solutions for continuous variables, which were expressed as fuzzy sets with known membership functions, provided optimal schemes for waste-flow allocations.These solutions were used for generating decision alternatives and thus helping decision makers to identify desired policies for MSW management under uncertainty.Although the case study in this study is just a hypothetical one, this case involved the main factors (e.g., transportation, operation, and capacity expansion costs) in MSW management.Consequently, the developed GFIP method can also be applied to real-world MSW management problems.Through the hypothetical (MSW) case study, the following advantages of the GFIP approach were presented: (i) GFIP could deal with fuzzy parameters with any kind of membership function shape; (ii) the solution process of GFIP would not generate complicated intermediate submodels; (iii) the computational requirement was reasonable; and (iv) the solutions of GFIP exhibited some distributional information, which was helpful for decision making.The above advantages indicated the usefulness and applicability of the developed GFIP approach in real-world MSW management problem.
Compared with the inexact mixed-integer linear programming (IMILP) method, uncertainties presented as fuzzy sets were incorporated within the GFIP's optimization processes.Solutions obtained from GFIP method specified some distributional information, which contained not only the lower and upper bounds with associated plausibility degrees but also the most possible values.
The developed GFIP could deal with various fuzzy sets with known membership functions, regardless of the shapes of these functions.However, it had difficulties in dealing with other kinds of uncertainties expressed as probabilistic distributions.Therefore, improvement for GFIP was further desired to enhance its capability of dealing with multiple uncertainties through incorporating interval programming or stochastic programming into the GFIP framework.Besides, the developed GFIP method cannot be used for evaluating the detailed effects of interactions among multiple uncertain inputs.Consequently, the factorial analysis method could then be integrated into the GFIP framework to address the above issue.Remark A.6.For any ã ∈  0 (), suppose two -cut levels  1 ,  2 are selected to cut ã; then these two -cuts can be formulated as: 1c) x = fuzzy continuous variables, x ∈ X ,  = 1, 2, . . .,  ( < ) (1d) x = fuzzy integer variables, x ∈ X,  =  + 1,  + 2, . . ., , MSW management under uncertaintyFuzzy parameters (e.g., waste generation, facility capacity, etc.)Permitting uncertain information in solutionsGeneralized fuzzy linear programming (GFLP) Generalized fuzzy integer programming (GFIP) model for MSW management Discretize membership grade into ( 1 ,  2 , . . ., , . . .,  q ) Reorder the -cut levels into Use  (i) to cut fuzzy parameters in GFIP model Convert the GFIP model into an inexact mixed integer linear programming (IMILP) model Interactive algorithm Lower-bound submodel (f − ) under  (i) -cut level Upper-bound submodel (f + ) under  (i) -cut level Obtain interval solutions for objective and decision variables under  (i) -cut level Generate membership function for continuous decision variables Generate optimal MSW management policies Stepwise interactive algorithm MSW disposal facility expansion Mixed integer linear programming (MILP) a sequence:  (1) ,  (2)

Figure 1 :
Figure 1: The schematic of the GFIP model for MSW management.

Figure 3 :
Figure 3: The lower bounds of waste flows under different -cut levels.

)Figure 5 :
Figure 5: The objective function values under different -cut levels.

Table 1 :
Waste generation levels and cost coefficients.

Table 2 :
Capacity expansion options and costs for landfill and WTE facilities.

Table 4 :
Capacity expanding options for the landfill under each -cut level.

Table 5 :
Capacity expanding options for the WTE facility under each -cut level.