A Hybrid Stochastic-Interval Analytic Hierarchy Process Approach for Prioritizing the Strategies of Reusing Treated Wastewater

This paper proposes a hybrid stochastic-interval analytic hierarchy process (SIAHP) approach to address uncertainty in group decision making by integrating interval judgment, probabilistic distribution, lexicographic goal programming, and Monte Carlo simulation. A case study related to wastewater treatment plant (WWTP) effluent reuse was conducted to demonstrate the feasibility of the proposed approach. Four candidate alternatives including city moat landscaping, municipal reuse, industrial reuse, and agricultural irrigation were evaluated by five experts according to technical, economic, and environmental criteria. The results suggest that industrial reuse (0.18–0.3) is more preferred over municipal reuse (0.16–0.25) or agricultural irrigation (0.17–0.26) in most replications. The final score of city moat landscaping ranges from 0.11 to 0.31 which indicates a great divergence of expert opinions. It can be concluded that choosing industrial reuse seems to give the best overall account of technical, economic, and environmental concerns.The proposed SIAHP approach can aid group decision making by accommodating linguistic information and dealing with insufficient information or biased opinions.


Introduction
As a major factor influencing social and economic development, the global water shortage has been exacerbated due to population growth, urbanization, climate change, and industrialization [1][2][3][4].The reuse of treated wastewater from wastewater treatment plants (WWTPs) has been gaining significant attention as a potential solution to cope up with the increasing water stress.Possible beneficial applications include, but are not limited to, toilet flushing, irrigation, groundwater recharge, fire protection, cooling water, landscaping, boiler feed water, and some working procedures of stocking, pulping, and paper making [5][6][7].It is of importance that conventional WWTPs are commonly designed to remove suspended solids and bulk organics rather than pathogenic microorganisms.Therefore, treated wastewater may still contain significant levels of fecal coliforms and many other pollutants (e.g., heavy metals) that need to be removed prior to reuse.Friedler et al. highlighted that improper planning of wastewater reuse may expose people to hazardous pathogenic microorganisms, posing a serious threat to public health [6].Kon and Watanabe argued that the promotion of reusing treated wastewater can be compromised by many factors such as the need for posttreatment, the concerns associated with water quality and functionality, and the economic efficiency [8].Former works have shown that the selection of reuse schemes depends not only on their technical, economic, and environmental feasibility but also mainly on public support, in other words, the decision makers who represent the interests of society.
Decision makers usually use multicriteria decision making (MCDM) methods to aid the decision making process [9][10][11][12][13][14][15].A widely used one is the analytic hierarchy process (AHP), a prominent tool for making decisions in situations involving multiple objectives [12].Many investigations have been undertaken in the past to apply the AHP to problems with high complexity and uncertainty, especially in the environmental sector [9,[16][17][18][19][20].Despite its numerous applications, there are critiques in the literature about using inconvenient crisp values and not including uncertainty analysis.Rosenbloom reported that alternatives are not statistically distinguishable if their scores are too close [21].That is to say, the AHP is not capable of addressing statistical interactions or feedback dependency between different hierarchies of a decision making problem [22,23].Another major drawback lies on its inability to translate the imprecision of the decision maker's perceptions [24].The exact nine-point AHP scale used in pairwise comparisons may fail to capture the imprecision or vagueness in the mind of respondents [25].
In response to these drawbacks, there have been many research attempts on incorporating interval uncertainty into the AHP.It has been well accepted that decision makers are more comfortable and confident to give interval judgments rather than to evaluate pairwise comparisons using single numeric values.Islam developed a lexicographic goal programming (LGP) approach in estimating weight vectors from interval pairwise comparison matrices (IPCMs) [5].Wang et al. stated that fuzzy inputs can also be addressed by IPCMs using -level sets and the extension principle [26].Chandran et al. presented a method to estimate the weights of IPCMs based on linear programming [27].Yu et al. advanced the LGP method by using multiplicative constraints to cope with the inherent deficiencies.On the other hand, to reflect the statistical characteristics of the traditional AHP, pairwise comparison elements were suggested to be viewed as random variables and computed via Monte Carlo simulation [9].The triangular distribution is the most commonly used distribution for modeling expert judgment in the AHP [28,29].It is advantageous over normal distribution and lognormal distribution due to its bounded nature.However, it may place too much emphasis on the most likely value [22].It is possible to overcome this disadvantage by using the uniform distribution, which is not much affected by the lack of information.To date, there has been no study specifically investigating how to handle both interval and probabilistic uncertainty in the AHP.Therefore, the objective of this study is to address such issue by integrating interval judgments, probabilistic distribution, the LGP method, and Monte Carlo simulation as a hybrid stochastic-interval analytic hierarchy process (SIAHP) approach.A case study of prioritizing the strategies of reusing treated wastewater from a WWTP in the city of Shuangcheng, China, was carried out to verify the feasibility and efficiency of the proposed approach.

Methodology
The proposed SIAHP approach uses expert knowledge as interval judgments and aggregates the lower and upper bounds into two independent arrays.For each nondiagonal element in the upper portion of each IPCM, the uniform distribution is assumed on both arrays such that random lower and upper bounds can be generated using Monte Carlo simulation.The lexicographic goal programming (LGP) method is then employed to estimate the priorities of each IPCM.The detailed steps are summarized as follows.
Step 1. Define the problem and structure the hierarchy from the top (i.e., goal) to the bottom (i.e., decision alternatives).
Step 2. Construct IPCMs for alternatives (or criteria) on each hierarchy level in terms of their performance against the ones on the level immediately above.Note that the IPCMs are reciprocal matrices which means judgments are required only for the upper portion: Step 3. Expert opinions are collected as intervals through questionnaires, interviews, surveys, and round-table discussion (2).The upper and lower bounds should not exceed the standard preference scale where 1 represents equal importance and 9 stands for absolute importance where  and  are the lower and upper bounds of interval judgments, respectively.
Step (e.g., normal, beta) that require more precise knowledge can be used.
Step 5.For each aforementioned upper portion element (e.g., ]), generate uniformly distributed pseudorandom numbers and map them to the ranges of the two arrays (i.e., the lower and upper bounds) predefined in Step 4. Note that the ranges may intersect and hence the random upper bound must be greater than or equal to the random lower bound.The lower portion of each IPCM can be determined by taking the reciprocal values of the corresponding upper portion.At the end of this step, a series of random IPCMs are obtained such that their non-diagonal elements are random intervals.
Step 6.The lexicographic goal programming (LGP) method is used to estimate the weight vector [23],  = ( 1 , . . .,   ), of each IPCM of rank .The interval judgments contained in the IPCMs can be defined by the following inequality: where   and   are the randomly generated lower and upper bounds of each element, respectively,   and   are real numbers representing the weights of the th and th entries, and deviation variables   and   are nonnegative real numbers.The weight vector can be derived by minimizing the inconsistency of the upper portion of the interval comparison matrix, in other words, a summation of all deviation variables   and   as stated in the following: Step 7. As with the traditional AHP, the inconsistency of judgements needs to be examined in order to ensure the validity of the random IPCMs.In this approach, the summations of all deviation variables  are expected to be less than or equal to 0.1; otherwise, the IPCM should be declined and regenerated.
Step 8.The final scores of each decision alternative can be calculated by aggregating the weights throughout the hierarchy as follows: where  (Score) is the final performance score for the th alternative,   is the merit of the th alternative with regard to the th criterion on the preceding level,  is the number of criteria, and   is the normalized weight of the th criterion against the goal.
Step 9. Repeat Steps 5 through 8 using Monte Carlo simulation for a number of replications.The overall scores of decision alternatives can be obtained as probability density functions.As one of the major cities in the province, the city of Shuangcheng has a population of 830 thousand with 640 thousand people dependent on agriculture.An advanced wastewater treatment plant associated with a drainage system was built in 2009 and has been operated since then to treat sewage and wastewater from individual residences, businesses, schools, and so forth.This wastewater treatment plant has a daily treatment capacity of 30,000 tonnes, which is sufficient for population growth and any future expansion demand.However, some important issues with regard to the operation of the treatment plant have been emerging, such as how to optimize the operating parameters to adapt to varying environmental conditions, how to effectively use the treated effluent, sludge, and other byproducts, and particularly how to efficiently operate this plant in a sustainable manner.In this case study, the local government invited experts to evaluate wastewater reuse options in order to explore the possibility of saving natural water resources in the future.

Goal, Decision Alternatives, and Evaluation Criteria.
The goal of this case study was to select the best alternative for WWTP effluent management and to provide the local government with a more systematic overview of sustainable development.Four reclaimed water reuse alternatives including city moat landscaping (CML), municipal reuse (MR), industrial reuse (IR), and agricultural irrigation (AI) were chosen based on the literature review and discussion with experts from local authorities and educational institutions (Figure 1).The city of Shuangcheng is surrounded by a city moat that uses a large amount of water every day to maintain water level and water quality.It is not just a landmark of the city's history but also a tourist attraction.Treated wastewater can be safely discharged into the moat to replenish its water supply and to improve water quality on a daily basis.and flushing toilets in public buildings and households.It should be noted that effluent from the WWTP usually needs to be further disinfected to such an extent that is safe from a microbiological point of view.Reclaimed water may also serve industrial processes (e.g., cooling, boiler feed, and process water) after proper treatment to prevent scaling, corrosion, bacterial fouling, and foaming issues.However, its feasibility and economic applicability vary from case to case.For example, the electroplating process requires water of high quality to wash circuit boards while tanneries can use low-quality water for washing hides.Lastly, agricultural irrigation using treated wastewater (reclaimed water) is becoming a common practice in arid and semiarid regions.Treated wastewater can be directly or indirectly (e.g., after disinfection) applied for irrigation depending on its quality and usage.As depicted in Figure 1, these four alternatives were evaluated based on technical, economic, and environmental criteria adopted from the literature and expert recommendations [7].

Data Acquisition.
Both quantitative and qualitative data were collected through questionnaire survey and the literature review.Key numbers such as capital cost and operational cost were obtained from different sources as shown in Table 1.They were normalized using a linear value function as follows: where   is the normalized performance score of the th alternative in terms of a criterion,  is the number of alternatives, and   is the numerical value of the th alternative.It should be noted that quantitative data has either increasing (e.g., benefit) or decreasing (e.g., cost, risk) behavior.Therefore, positive scores (between 0 and 1) were given to those with increasing values, while negative scores (between −1 and 0) were assigned to the decreasing values.Normalized values for different alternatives with respect to reclaimed wastewater reuse are reported in Table 2.In the absence of available quantitative data, expert subjective rankings were used to evaluate the alternatives in terms of other criteria.Five experts from local authorities (Bureau of Environmental Protection) and educational institutions (professors and graduate students at Memorial University of Newfoundland) were asked to assess the hierarchy with interval judgments.Equation (7) and Table 3 show an example of expert assessments in regard to human risk.The number of Monte Carlo iterations used in

Results and Discussion
As shown in Figure 2, the alternative priorities with respect to some subcriteria are depicted in the form of probability density estimates using the kernel-smoothing method.It clearly shows that most density estimates have two major peaks, indicating that the optimization results obtained from the LGP method tend to be concentrated within a certain interval possibly due to the randomly generated interval judgments.A few density estimates appear to be smooth and continuous over a large span while many others are steep and bell shaped.This may be attributable to expert assessments which diverge on unfamiliar knowledge (e.g., the application of city moat landscaping) but converge on some commonly agreed upon principles.The overall alternative priorities with regard to the whole sets of technical, economic, and environmental criteria are synthesized in Figure 3. City moat landscaping and agricultural irrigation have the best technical performance in most replications mainly due to their ease of operation and maintenance.Their overlaps with other alternatives are not statistically significant (Figure 3(a)); however, their own overlap reveals itself with a strong correlation coefficient of (0.956).Municipal reuse is seemingly preferred over the other options in terms of economic concerns (Figure 3(b)).This is more or less in agreement with the reality that reusing reclaimed water for municipal purposes can save extra costs of posttreatment and therefore creates more economic benefits.The correlation coefficients between the alternatives are all above 0.950 which suggests that their performance is positively correlated such that the increase of one is associated with the increase of the others.This can be further interpreted that the alternatives are not well distinguishable in terms of economic considerations; therefore, the significance of economic criteria is attenuated as demonstrated in Figures 3(b) and 4(d).From the perspective of environmental conservation, industrial reuse dominates the others in more than 95% of the replications (Figure 3(c)).This seems to fit the way most people think because using treated wastewater in industrial sectors can hardly cause any environmental issues or raise human health concerns.
The probability density distributions of subcriterion weights are demonstrated in Figure 4 using the kernelsmoothing method.Applicability and water quality requirement are the most influential technical subcriteria, followed by reliability.Benefit absolutely dominates other two economic subcriteria with the confidence level of 100%.Environmental subcriteria are to some extent tangle with each other while water reservation manages to lead in 89.2% of the replications.These aforementioned subcriteria should be prioritized and further considered in the decision making process if the final alternative scores are too close.The final scores of each decision alternative are plotted in Figure 5.It appears that industrial reuse (0.18-0.3) is more preferred over municipal reuse (0.16-0.25) or agricultural irrigation (0.17-0.26) though the existence of overlaps may contradict this assertion.Further analysis indicates that industrial reuse takes the first or second places in over 88.6% of the replications while this number for the other two options is somewhere between 20% and 30% (Table 4).Another worthy mentioned finding is related to the final score of city moat landscaping which spans over a large range from 0.11 to 0.31 (Figure 5).This range is larger than those of other alternatives mainly because of the great divergence of expert opinions.In other words, some experts may favor city moat over the others (in 54.2% replications), whereas some experts may rank it as the least preferable solution (in 35.3% replications) (Table 4).Figure 6 confirms these observations by using box plot to graphically illustrate the minimum, lower quartiles, medians, upper quartiles, and maximum of alternative overall scores.It can be argued that choosing industrial reuse seems to give the best overall account of technical, economic, and environmental concerns.

Mathematical Problems in Engineering
It is worth discussing the correlations between the final scores as this can help decision makers further interpret the experts' preferences.City moat landscaping and industrial reuse are in conflict with each other as their correlation coefficient is −0.858, implying that experts tend not to choose both options as their favorite (Figure 7(a)).The correlation coefficients between industrial reuse and municipal reuse, industrial reuse and agricultural irrigation, and municipal and agricultural irrigation are 0.701, 0.508, and 0.835, respectively, which agrees with the positive correlations shown in Figures 7(b)-7(d).That is to say, when either one of these three is favored by an expert, it is highly probable that the other two are also more or less favored.

Conclusions
This paper presents a hybrid stochastic-interval analytic hierarchy process (SIAHP) approach to address uncertainty issues in the decision making procedure by integrating interval judgments, probabilistic distributions, the lexicographic goal programming method, and Monte Carlo simulation.A case study of selecting the best effluent reuse strategy for a WWTP in the city of Shuangcheng, China, was carried out to verify the feasibility and efficiency of the proposed approach.Four candidate alternatives including city moat landscaping, municipal reuse, industrial reuse, and agricultural irrigation were evaluated by five experts from local authorities and educational institutions.Technical, economic, and environmental criteria and their associated sub-criteria were considered to investigate the performance of the alternatives.The results show that a few density estimates of alternative performance with regard to several subcriteria appear to be smooth and continuous over a large span while many others are steep and bell shaped.This may be attributable to expert assessments that diverge on unfamiliar knowledge but converge on some commonly agreed upon principles.City moat landscaping and agricultural irrigation have the best technical performance in most replications mainly due to their ease of operation and maintenance.Their own overlap is evident with a strong correlation coefficient of 0.956.Municipal reuse is seemingly preferred over the other options in terms of economic concerns while the correlation coefficients between all alternatives are above 0.950.This can be further interpreted that the alternatives are not well distinguishable in terms of economic considerations; therefore, the significance of economic criteria is attenuated.Industrial reuse seems to be the best option in terms of environmental reservation.Based on the final overall scores, industrial reuse (0.18-0.3) is more preferred than municipal reuse (0.16-0.25) or agricultural irrigation (0.17-0.26) in most replications.The final score of city moat landscaping ranges from 0.11 to 0.31 which is larger than those of other alternatives mainly because of the great divergence of expert opinions.In other words, some experts may favor it over the others (in 54.2% replications), whereas some experts may rank it as the least preferable solution (in 35.3% replications).It can be concluded that choosing industrial reuse seems to give the best overall account of technical, economic, and environmental concerns.The proposed SIAHP approach can aid group decision making by accommodating linguistic information and dealing with insufficient information or biased opinions.

Figure 2 :
Figure 2: Probability density estimates of alternative scores with regard to the subcriterion of (a) applicability, (b) human risk, (c) operational cost, and (d) water quality requirement.

Figure 3 :
Figure 3: Probability density estimates of the overall alternative scores with regard to (a) technical criteria, (b) economic criteria, and (c) environmental criteria.

Figure 4 :
Figure 4: Probability density estimates of the relative weights between each (a) technical subcriterion, (b) economic sub-criterion, (c) environmental subcriterion, and (d) criteria group.

Figure 5 :Figure 6 :
Figure 5: Probability density estimates of the final scores of decision alternatives.

Figure 7 :
Figure 7: Scatter plots of the final scores of (a) city moat landscaping versus industrial reuse, (b) municipal reuse versus industrial reuse, (c) agricultural reuse versus industrial reuse, and (d) agricultural reuse versus municipal reuse.

Table 1 :
Quantitative database for reclaimed wastewater reuse.

Table 2 :
Alternatives' normalized scores for quantitative criteria.

Table 3 :
Uniform distributions (in MATLAB) of the PCM with regard to human risk.

Table 4 :
Summary of the evaluation results of the SIAHP approach.