Journal of Applied Mathematics

Volume 2013 (2013), Article ID 731568, 14 pages

http://dx.doi.org/10.1155/2013/731568

Research Article

Mathematical Modeling for Water Quality Management under Interval and Fuzzy Uncertainties

J. Liu, Y. P. Li, and G. H. Huang

MOE Key Laboratory of Regional Energy Systems Optimization, Resources and Environmental Research Academy, North China Electric Power University, Beijing 102206, China

Received 4 October 2013; Accepted 29 November 2013

Academic Editor: X. S. Qin

Copyright © 2013 J. Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

In this study, an interval fuzzy credibility-constrained programming (IFCP) method is developed for river water quality management. IFCP is derived from incorporating techniques of fuzzy credibility-constrained programming (FCP) and interval-parameter programming (IPP) within a general optimization framework. IFCP is capable of tackling uncertainties presented as interval numbers and possibility distributions as well as analyzing the reliability of satisfying (or the risk of violating) system’s constraints. A real-world case for water quality management planning of the Xiangxi River in the Three Gorges Reservoir Region (which faces severe water quality problems due to pollution from point and nonpoint sources) is then conducted for demonstrating the applicability of the developed method. The results demonstrate that high biological oxygen demand (BOD) discharge is observed at the Baishahe chemical plant and Gufu wastewater treatment plant. For nonpoint sources, crop farming generates large amounts of total phosphorus (TP) and total nitrogen (TN). The results are helpful for managers in not only making decisions of effluent discharges from point and nonpoint sources but also gaining insight into the tradeoff between system benefit and environmental requirement.

1. Introduction

Water is one of the most essential constituents for the human life, which is crucial to various socioeconomic issues such as industrial production, agricultural activity, environmental protection, and regional sustainability. In recent years, especially in China, degradation of water quality due to point and nonpoint source pollutions has become one of the most pressing environmental concerns. According to the 2010 Report on Water Environmental Quality of China, approximately 33.8% of the monitored river water (204 rivers with 409 monitoring stations) is in the worst two categories of water quality classification system (i.e., no longer fishable and of questionable agricultural value); around 53.8% of the assessed lakes and reservoirs are subject to different degrees of eutrophication [1]. Under such a circumstance, water quality management is an essential task for preserving valuable water resources and facilitating sustainable socioeconomic development in watershed systems [2]. In fact, water quality planning efforts are complicated with a variety of uncertainties, which may be derived from the random characteristics of natural processes (i.e., precipitation and climate change) and stream conditions (i.e., stream flow, water supply, and point/nonpoint source pollution), the errors in estimated modeling parameters, and the vagueness of system objectives and constraints. In general, the system objectives are often associated with a number of socioeconomic and ecological factors such as economic return, environmental protection, and ecological sustainability, while the constraints are related to pollutant discharges, soil losses, resources availabilities, environmental requirements, and policy regulations. Moreover, these uncertainties may be further amplified by not only interactions among various uncertain and dynamic impact factors, but also their associations with economic implications of violated environmental requirements [3].

Fuzzy mathematic programming (FMP), based on fuzzy sets theory, can facilitate the analysis of system associated with uncertainties being derived from vagueness or fuzziness [4]. FMP method is suitable for situation when the uncertainties cannot be expressed as probability density functions (PDFs), such that adoption of fuzzy membership functions becomes an attractive alternative [5]. Previously, a number of FMP methods were developed for water quality management [5–15]. For example, Julien [6] investigated the application of a fuzzy possiblistic programming to address imprecise parameters which were represented by possibility distributions in water quality decision-making problems. Mujumdar and Sasikumar [8] formulated a fuzzy flexible optimization model for dealing with the system’s fuzzy goal and constraints in a water quality management problem. Nie et al. [12] proposed a fuzzy robust optimization model for water quality management of an agricultural system to deal with uncertainties expressed as fuzzy membership functions in both left- and right-hand-side coefficients (of the model’s constraints). Maeda et al. [13] employed a fuzzy flexible optimization model which involved fuzzy set theory to express vagueness in constraints and objectives in river water quality management problems. Liu et al. [15] developed a two-stage fuzzy robust programming model for water quality management to address fuzzy parameters which were represented by possibility distributions in the left- and right-hand sides of the constraints.

Generally, FMP methods can be classified into three categories in view of the forms of uncertainties: (i) fuzzy flexible programming, (ii) robust programming, and (iii) fuzzy possibilistic programming. In detail, fuzzy flexible programming can deal with decision problems under fuzzy goal and constraints; however, it has difficulties in tackling ambiguous coefficients of the objective function and constraints. Robust programming improves upon fuzzy flexible programming by allowing fuzzy parameters in the constraints to be represented by possibility distributions. However, the main limitations of this method remain in its difficulties in tackling uncertainties in a nonfuzzy decision space. In fuzzy possibilistic programming, fuzzy parameters that are regarded as possibility distributions are introduced into the modeling frameworks. It can handle ambiguous coefficients in the left- and right-hand sides of the constraints and in the objective function.

Fuzzy credibility-constrained programming (FCP) is a computationally efficient fuzzy possibilistic programming approach that relies on mathematical concepts (i.e., the expected value of a fuzzy number and the credibility measure) and can support different kinds of fuzzy members such as triangular and trapezoidal forms as well as enabling the decision maker to satisfy some constraints in at least some given confidence levels [16]. When the credibility value of a fuzzy event reaches 1, the fuzzy event will certainly occur; when the credibility value of a fuzzy event reaches 0, the fuzzy event will not occur. For example, with respect to water quality management, if the allowable total phosphorus (TP) discharges are 7.0, 7.6, and 8.2 kg/day, and the amount of actual discharge may be 7.8 kg/day, then the credibility degree of the event, wherein the total phosphorus discharge can satisfy the river’s self-purification capacity, would be 0.33; the credibility degree of its complement event (water pollution occur) would be 0.67. No feature of fuzzy sets would be missing by using credibility measure [16, 17]. However, the main limitation of FCP lies within its deterministic coefficients for the objective function, leading to potential losses of valuable uncertain information; besides, when many uncertain parameters are expressed as fuzzy sets, interactions among these uncertainties may lead to serious complexities, particularly for large-scale practical problems [18]. In fact, in water quality management problems, uncertainty is an inherent component of any economic analysis, particularly those (e.g., effluent trading programs) associated with environmental policy and project appraisal [19]. For example, one major characteristic of nonpoint source pollution that differs from point source pollution is imperfect knowledge about pollutant loadings; the crop productivity and cost benefit coefficients are easier to be presented as intervals than by membership functions. Interval-parameter programming (IPP) is an alternative for handling uncertainties in the model’s left- and/or right-hand sides as well as those that cannot be quantified as membership or distribution functions, since interval numbers are acceptable as its uncertain inputs [20].

Therefore, the objective of this study is to develop an inexact fuzzy credibility-constrained programming (IFCP) method for water quality management, through coupling fuzzy credibility-constrained programming (FCP) with interval-parameter programming (IPP). The main advantage of IFCP is that it can effectively handle uncertain parameters expressed as both fuzzy sets and interval values in the objective function and constraints. IFCP would not lead to serious complexities in its solution process, and it is applicable to large-scale practical problems. Then, the developed IFCP method is applied to a real-world case of water quality management of the Xiangxi River, which faces severe water quality problems due to point and nonpoint source pollution. The results obtained can help decision makers to generate alternatives for industrial production scale, water supply, cropped area, livestock husbandry size, and manure/fertilizer application rate, with consideration of river water quality management.

The paper will be organized as follows: Section 2 describes the development process of the IFCP; Section 3 provides a case study of river water quality management; Section 4 presents result analysis and discussion; Section 5 draws some conclusions and extensions.

2. Methodology

When coefficients in the constraints are ambiguous and can be expressed as possibility distributions, the problem can be formulated as a fuzzy credibility-constrained programming (FCP) model as follows:

subject to:

where is a vector of nonfuzzy decision variables, are benefit coefficients, are technical coefficients, and are right-hand-side coefficients. Some or all of these coefficients can be fuzzy numbers. Cr is the credibility measure which is firstly proposed and was widely used in many research areas [16]. Let be a fuzzy variable with membership function , and let be real numbers. The credibility measure can be defined as follows [17]:

Noteworthy, since and , the credibility measure can be defined as follows:

Similar to the probability measure,

Consider a triangular fuzzy variable since it is the most popular possibility distribution, the fuzzy variable fully determined by the triplet () of crisp numbers with whose membership function is given by

Based on this membership function, credibility of can be expressed by

The inverse function of the credibility measure is , when . Normally, it is assumed that a significant credibility level should be greater than 0.5. Therefore, (6) can be written as

Then, (7) can be transformed into a deterministic constraint as follows:

The fuzzy credibility-constrained programming (FCP) model can be formulated as follows:

subject to

Obviously, model (9a), (9b), and (9c) can effectively deal with uncertainties in the right-hand sides presented as fuzzy sets when coefficients in the left-hand sides and in the objective function are deterministic. However, in real-world optimization problems, uncertainties may exist in both left- and right-hand sides (of the constraints) as well as objective-function coefficients; moreover, the quality of information that can be obtained is mostly not satisfactory enough to be presented as fuzzy membership functions [2]. For example, economic return, pollutant discharge, and resources availability are easier to be expressed as intervals than membership functions [21]. Since interval-parameter programming (IPP) is useful for addressing uncertainties expressed as interval values in modeling parameters, it can be integrated into the FCP model to deal with uncertainties presented in fuzzy and interval formats. Then, an interval fuzzy credibility-constrained programming (IFCP) can be formulated as follows:

subject to

where the “−” and “+” superscripts represent the lower- and upper-bounds of interval parameters/variables, respectively. Then, a two-step solution method is proposed for facilitating computations of the IFCP model. The first submodel can be formulated as follows:

subject to

where  ; < 0; when ;= 1 when ; is the lower bound of the credibility level value. The optimal solutions of the first submodel would be   () and   . In the second step, the submodel corresponding to can be formulated:

subject to

The optimal solutions of model (11a), (11b), (11c), and (11d) would be and can be obtained. Through integrating the solutions of the two submodels, the solution for the objective-function value and decision variables can be obtained as follows:

3. Case Study

The Xiangxi River (which ranges in longitude from 110°25′ to 111°06′ E and in latitude from 30°57′ to 31°34′ N) is located at 40 km upstream of the Three Gorges Reservoir [22]. It is 94 km long with a catchment area of 3099 km², and its elevation generally ranges from 154 m to 3000 m. It is located in the subtropical continental monsoon climate zone, with an annual temperature of 15.6°C (from 1961 to 2004) and the long-term annual mean runoff depth of 688 mm [23]. Moreover, it is one of the rainiest centers in the west of Hubei province, with an average annual precipitation ranges from 900 mm to 1200 mm. The temporal precipitation of this basin is uneven, which varies largely among different seasons. For example, more than 41% precipitation occurs in June to August, and the rainfall in the spring, autumn, and winter seasons occupy 28%, 26%, and 5% of the total precipitation per year, respectively. There are plenty of mineral resources (i.e., phosphate ore, coal, pyrite, and granite), where the reserve of phosphorite is among the top three in China, which reaches 357 million (i.e., tonne) [1]. Relying on these advantages, the number of phosphorus mining companies and related chemical plants are increasing along the banks of the Xiangxi River. Besides, multiple crops such as rice, maize, wheat, citrus, tea, potato, and vegetable are cultivated in the catchment since the land-use patterns are diverse, and the tillable area land is approximately 294.5 km². In addition, pig, ox, sheep, and domestic fowl are the main live stocks in animal husbandry. The main pattern is scattered livestock breeding instead of large-scale standardized breeding.

Currently, water quality problems due to point and nonpoint source pollution discharges become more and more challenging in this catchment. Main point sources include five chemical plants (i.e., GF, BSH, PYK, LCP, and XJLY), six phosphorus mining companies (i.e., XL, XH, XC, GP, JJW, and SJS), and four wastewater treatment plants (WTPs) (i.e., Gufu, Nanyang, Gaoyang, and Xiakou), while four agricultural zones (AZ1 to AZ4) are the main nonpoint sources due to the application of manure/fertilizer. These point and nonpoint sources scatter along a length of about 51 km river stretch which is segmented into five reaches, and the reaches are marked as I to V. The main water quality problems include (i) the immoderate discharge of high-concentration phosphorus-containing wastewater and industrial soil wastes (i.e., chemical wastes, slags, and tailings) far exceed what can be decomposed by self-purification (according to the field investigation, 23.89 of phosphorus enters downstream of the Xiangxi River); (ii) high potential for generating soil erosion and surface runoff due to the special geography and heavy rainfall (i.e., the average erosion modulus reaches 6,488 t/(km²·a) in Xiangshan County in this catchment); (iii) large amounts of nutrient pollutants (in terms of phosphorus and nitrogen) in livestock wastewater and wastes (from pig, ox, sheep, and domestic fowl breeding) are drained into the river by direct discharge or in rainfall. In such a circumstance, decision makers should seek to develop a sound pollution control plan to ameliorate the current situation of the water environment since it is infeasible and technical impossible to ensure zero emission of pollutants.

In this study, the planning horizon is one year. Moreover, since some crops should be grown in dry season, while some other crops should be cultivated in wet season, two periods are chosen to cover the planning horizon. The first period is from June to October (i.e., dry season), and the second period is from November to May of the next year (i.e., wet season). The objective is to maximize the net system benefit subject to the environmental requirements under uncertainty over the planning horizon. Policies in terms of the related human activities (i.e., industrial, municipal, and agricultural activities) and the pollutant discharges (from fifteen point and four nonpoint sources) are critical for ensuring a maximum system benefit and a safe water quality [24]. Based on field investigations and related literatures, biological oxygen demand (BOD), total nitrogen (TN), and total phosphorus (TP) are selected as water quality indicators [1, 25]. To develop the local economy in a sustainable manner, pollutant discharge should be controlled by setting the thresholds for TP, TN, and BOD discharge in each reach [26]. However, human-induced imprecision in acquiring these thresholds (i.e., lack of available data and biased judgment) make it more complicate.

On the other hand, uncertainties in the study system include the following: (a) cost of wastewater treatment, manure, and fertilizer purchase are associated with many uncertain factors, which are expressed as interval numbers (e.g., an interval of RMB¥/t is denoted as cost of manure purchase of AZ1 in dry season); (b) the BOD and TP treatment efficiencies of wastewater in WTPs and chemical plants are related to operating conditions of the treatment facilities, which cannot be obtained as deterministic numbers (e.g., an interval of is denoted as treatment efficiencies of wastewater in Xiakou WTP); (c) nonpoint source losses of nitrogen and phosphorus from agricultural zones fluctuate dynamically due to variability in soil erosion (corresponding to solid-phase nitrogen) and surface runoff (corresponding to dissolved nitrogen) (e.g., average soil loss from AZ1 planted with citrus in dry season would be  t/ha), and runoff from AZ1 planted with citrus in dry season would be  mm); (d) the amount of fertilizer and manure applications may vary with the soil fertility to meet the nutrient demands of each crop (i.e., nitrogen and phosphorus); (e) energy and digestible protein demands of human and animals are determined by crops’ yield (e.g., yield of citrus planted in AZ1 during dry season would be  t/ha). Therefore, based on the IFCP method developed in Section 2, the study problem can be formulated as:

 subject to(1)wastewater treatment capacity constraints:

(2)BOD discharge constraints:

(3)nitrogen discharge constraints:

(4)phosphorus discharge constraints:

(5)soil loss constraints:

(6)fertilizer and manure constraints:

(7)energy and digestible protein constraints:

(8)production scale constraints:

(9)total yield of crops:

(10)planning area constraints:

(11)nonnegative constraints:

where : chemical plant; Gufu (GF); Baishahe (BSH); Pingyikou (PYK); Liucaopo (LCP); Xiangjinlianying (XJLY); : agricultural zone; ; : main crop; citrus; tea; wheat; potato; rapeseed; alpine rice; second rice; maize; vegetables; : phosphorus mining company; Xinglong (XL); Xinghe (XH); Xingchang (XC); Geping (GP); Jiangjiawan (JJW); Shenjiashan (SJS); : livestock; pig; ox; sheep; domestic fowls; : town; Gufu; Nanyang; Gaoyang; Xiakou; : time period; dry season; wet season; : length of period (day); : net benefit from chemical plant during period (RMB¥/t); : production level of chemical plant during period (t/day); : benefit from water supply to municipal uses (RMB¥/m³); : quantity of water supply to town in period (m³/day); : average benefit for per unit phosphate ore (RMB¥/t); : production level of phosphorus mining company during period (t/day); : average benefit from livestock (RMB¥/unit); : number of livestock in the study area (unit); : yield of crop planted in agricultural zone during period (t/ha); : average benefit of agricultural product (RMB¥/t); : planting area of crop  in agricultural zone during period (ha); : cost of manure collection/disposal in agricultural zone during period (RMB¥/t); : cost of purchasing fertilizer in agricultural zone during period (RMB¥/t); : amount of manure applied to agricultural zone with crop during period (t); : amount of fertilizer applied to agricultural zone with crop during period (t); : wastewater generation rate of chemical plant during period (m³/t); : wastewater treatment cost of chemical plant during period (RMB¥/m³); : wastewater discharge amount at town during period (m³/m³); : cost of municipal wastewater treatment (RMB¥/m³); : irrigation quota for crop in zone during period (m³/ha); : water supply price (RMB¥/m³); : water consumption of per unit production of chemical plant during period (m³/t); : BOD concentration of raw wastewater from chemical plant in period (kg/m³); : BOD treatment efficiency in chemical plant during period (%); : allowable BOD discharge for chemical plant in period (kg/day); : BOD concentration of municipal wastewater at town during period (kg/m³); : BOD treatment efficiency of WTPs at town during period (%); : allowable BOD discharge for WTPs at town during period (kg/day); : nitrogen content of soil in agricultural zone planted with crop (%); : average soil loss from agricultural zone planted with crop in period (t/ha); : runoff from agricultural zone with crop in period (mm); : dissolved nitrogen concentration in the runoff from agricultural zone planted with crop in period (mg/L); : maximum allowable nitrogen loss in agricultural zone   during period (t/ha); : tillable area of agricultural zone during period (ha); : amount of manure generated by livestock [t/(unit·day⁻¹)]; : amount of manure generated by humans [t/(unit·day⁻¹)]; : total rural population in the study area during period (unit); : manure loss rate in period (%); : wastewater generation of per capita water consumption during period [m³/(unit·day⁻¹)]; : dissolved nitrogen concentration of rural wastewater during period (t/m³); : maximum allowable nitrogen loss from rural life section in period (t); : phosphorus concentration of raw wastewater from chemical plant in period (kg/m³); : phosphorus treatment efficiency in chemical plant in period (%); : amount of slag discharged by chemical plant in period (kg/t); : slag loss rate due to rain wash in chemical plant during period (%); : phosphorus content in slag generated by chemical plant in period (%); : allowable phosphorus discharge for chemical plant in period (kg/day); : phosphorus content of soil in agricultural zone planted with crop (%); : average soil loss from agricultural zone planted with crop in period (t/ha); : runoff from agricultural zone with crop in period (mm); : dissolved phosphorus concentration in the runoff from agricultural zone planted with crop in period (mg/L); : maximum allowable phosphorus loss in agricultural zone during period (t/ha); : dissolved phosphorus concentration of rural wastewater during period (t/m³); : maximum allowable phosphorus loss from rural life during period (t); : phosphorus concentration of municipal wastewater at town in period (kg/m³); : phosphorus treatment efficiency of WTP at town in period (%); : allowable phosphorus discharge for WTP at town in period   (kg/day); : wastewater generation from phosphorus mining company in period (m³/t); : phosphorus concentration of wastewater from mining company in period (kg/m³); : phosphorus treatment efficiency in mining company (%); : amount of slag discharged by mining company during period (kg/t); : phosphorus content in generated slag (%); : slag loss rate due to rain wash (%); : allowable phosphorus discharge for mining company during period (kg/day); : maximum allowable soil loss agricultural zone in period (t/ha); : nitrogen volatilization/denitrification rate of fertilizer in period (%); : nitrogen volatilization/denitrification rate of manure in period (%); : nitrogen content of fertilizer (%); : phosphorus content of fertilizer (%); : nitrogen content of manure (%); : phosphorus content of manure (%); : nitrogen requirement of agricultural zone with crop during period (t/ha); : phosphorus requirement of crop in agricultural zone during period (t/ha); : net energy content of crop (Mcal/t); : net energy requirement of livestock (Mcal/unit); : net energy requirement of human beings (Mcal/unit); : digestible protein content of crop (%); : digestible protein requirement of livestock (t/unit); : digestible protein requirement of human beings (t/unit); : minimum crop production requirement for crop in period (t); : capacity of wastewater treatment capacity (WTPs) (m³/day); : capacity of wastewater treatment capacity (chemical plants) (m³/day).

Table 1 provides the net benefits of chemical plant productions, phosphorus mining company, and municipal water supply. According to the crops’ growth periods, dry season crops include wheat, potato, rapeseed and alpine rice; second rice, maize and vegetable are identified as main crops during the wet season; citrus and tea grow over the entire year. Table 2 shows the net benefits of each crop in every AZ and yields of each crop. To guarantee the stream water quality, wastewater treatment measures have to be adopted at each point source. Based on the local environmental regulations, a safe level of water quality must be guaranteed to protect aquatic life and maintain aerobic condition in the stream system [2]. Thus, the BOD and TP loading amount would be controlled strictly. However, the imprecision of the allowable BOD and TP discharge could introduce uncertainties in the water quality management. For modeling purpose, the vagueness of the allowable BOD and TP loading amount are encoded by triangular fuzzy membership functions. Figure 1 presents the fuzzy set with triangular membership function. The minimum, maximum, and most likely values (, , and ) that define these fuzzy sets are estimated according to previous research regarding water quality monitoring and environmental capacities, as tabulated in Table 3. Moreover, by setting acceptable interval credibility levels, the constraints can be at least basically satisfied and at best practically satisfied. Lower-bound of the interval numbers would be no less than 0.5, and upper- bound of the interval numbers would be no more than 1 [16]. Lower-bound of credibility level represents a situation when the decision makers are optimistic about this study area, which may imply a higher risk of violating the river’s self-purification capacity. Conversely, upper-bound of credibility level corresponds to a situation when the decision makers prefer a conservative policy that could guarantee that the river’s self-purification capacity be satisfied.

Table 1: Net benefits from each production.

Period

Net benefits from chemical plant (RMB/t)

GF
718.6, 834.5
743.75, 876.2

BSH
1291.5, 1499.7
1336.7, 1574.8

PYK
743.0, 862.8
769.0, 905.9

LCP
1324.2, 1537.5
1370.6, 1614.4

XJLY
1524.4, 1770.0
1578.7, 1858.5

Net benefits from water supply (RMB/m³)

Gufu
39.1, 42.1
43.4, 47.7

Nanyang
29.0, 31.2
32.2, 35.4

Gaoyang
35.3, 38.0
39.2, 43.1

Xiakou
32.1, 34.6
35.7, 39.3

Net benefits from phosphorus mining company (RMB/t)

XL (Xinglong, )
150, 173
147, 180

XH (Xinghe, )
126, 145
130, 156

XC (Xingchang, )
135, 155
135, 162

GP (Geping, )
144, 166
150, 180

JJW (Jiangjiawan, )
137, 158
141, 169

SJS (Shenjiashan, )
140, 164
145, 175

Table 2: Crop yields and net benefits.

Agricultural zone

AZ1
AZ2
AZ3
AZ4

Yields (t/ha)

    Citrus
10.3, 12.6
9.8, 11.9
8.4, 10.2
9.2, 11.2

    Tea
0.1, 0.2
0.1, 0.2
0.3, 0.6
0.2, 0.4

    Wheat
1.4, 2.1
1.3, 1.7
1.2, 2.5
2.3, 3.1

    Potato
2.4, 3.2
2.0, 2.8
1.9, 2.8
2.8, 3.8

    Rapeseed
1.6, 1.7
1.3, 1.5
1.3, 1.7
1.7, 1.8

    Alpine rice
6.9, 9.6
6.7, 7.9
4.4, 6.6
5.7, 7.0

    Second rice
7.4, 10.2
7.3, 8.5
4.7, 7.0
6.2, 7.4

    Maize
2.9, 3.9
2.7, 2.9
2.5, 2.6
4.0, 4.7

    Vegetable
21.6, 24.0
10.4, 13.4
21.8, 24.2
26.2, 29.1

Net benefits (RMB/t)

    Citrus
1082, 1350
1010, 1260
1126, 1405
1203, 1502

    Tea
18000, 18900
19125, 20081
20170, 21178
22500, 23625

    Wheat
1031, 1242
972, 1169
945, 1134
1010, 1218

    Potato
915, 1001
827, 897
893, 969
846, 917

    Rapeseed
2435, 2934
2191, 2640
2305, 2778
2490, 3000

    Alpine rice
1105, 1214
1047, 1153
1084, 1178
1149, 1265

    Second rice
1159, 1274
1098, 1209
1137, 1236
1205, 1327

    Maize
1295, 1451
1267, 1422
1230, 1378
1308, 1467

    Vegetable
2360, 2691
1613, 1865
1927, 2197
1845, 1950

Table 3: Allowable BOD and TP discharges.

Allowable BOD loading discharge

from chemical plant (kg/day)

GF
0.76
0.78
0.82
0.85
0.88
0.92

BSH
437.2
442.97
444.265
448.605
451.33
454.24

PYK
7.25
7.83
7.53
7.89
7.81
7.95

LCP
95.43
97.23
97.82
102.83
100.21
108.43

XJLY
48.34
49.43
49.635
50.645
50.93
51.86

Allowable BOD loading discharge

from WTP (kg/day)

Gufu
145
143
150
148
155
153

Nanyang
7
6
10
9
13
12

Gaoyang
38
36
40
38
42
40

Xiakou
28
30
30
32
32
34

Allowable TP loading discharge

from chemical plant (kg/day)

GF
1.18
1.29
1.31
1.43
1.44
1.57

BSH
393.73
485.94
419.14
514.27
444.55
542.59

PYK
80.23
99.57
84.87
105.75
89.51
111.92

LCP
385.59
470.16
416.21
497.96
446.83
525.76

XJLY
259.81
330.52
276.46
355.02
293.1
379.52

Allowable TP loading discharge

from WTP (kg/day)

Gufu
13.00
13.50
15.22
15.68
17.43
17.85

Nanyang
0.36
0.42
0.55
0.6
0.74
0.78

Gaoyang
1.82
1.88
2.03
2.27
2.23
2.65

Xiakou
2.05
2.34
2.25
2.49
2.45
2.64

Allowable TP loading discharge

from phosphorus mining company (kg/day)

XL
157.70
143.31
175.60
192.45
193.49
241.59

XH
89.40
104.19
99.51
120.53
109.62
136.87

XC
108.20
110.74
120.48
138.25
132.75
165.76

GP
146.20
140.03
162.78
181.99
179.36
223.95

JJW
89.80
104.49
99.97
121.01
110.13
137.52

SJS
86.60
102.19
96.41
117.40
106.21
132.61

Allowable TP and TN loading discharge

from rural life (kg/day)

TP
21.5
20.5
24.45
24.4
27.4
28.3

TN
102.1
106.0
107.5
107.0
112.8
118.0

Figure 1: Fuzzy set and credibility.

4. Result Analysis

Table 4 shows the solutions for industrial production and water supply during the two periods. The results show that significant variations in industrial production and water supply exist among different chemical plants, phosphorus mining companies, and towns. For chemical plants, the production scale of BSH (i.e.,  t/d (tonne/day) during period 1 and  t/d during period 2) would be larger than the other chemical plants (especially GF) because of its higher allowable BOD and TP discharges and net benefits. For water supply, more water would be delivered to Gufu (i.e.,  m³/day in period 1 and  m³/day in period 2) than those to the other towns due to its greater water demand, higher allowable BOD discharge, and higher economic return. For phosphorus mining company, the production scale of XL (i.e.,  t/d in period 1 and  t/d in period 2) would be larger than the other chemical plants because of its higher allowable TP discharge and net benefit.

Table 4: Solutions for industrial productions and water supplies.

Main point sources
Period

Chemical plant (t/day)

GF
25.64, 27.41
21.13, 27.02

BSH
536.23, 652.37
530.78, 629.19

PYK
81.71, 102.09
75.98, 80.92

LCP
379.24, 403.72
357.33, 453.48

XJLY
280.87, 313.47
320.42, 337.81

Phosphorus mining company (t/day)

XL
824.78, 1261.00
701.05, 1044.89

XH
326.00, 420.00
350.48, 400.00

XC
286.87, 400.00
260.38, 393.68

GP
534.99, 778.00
659.80, 778.00

JJW
425.99, 547.00
412.57, 547.00

SJS
510.53, 685.00
602.43, 685.00

Water supply (m³/day)

Gufu
11723.79, 14859.14
12264.89, 21739.14

Nanyang
762.29, 987.64
859.73, 1109.84

Gaoyang
2783.46, 3225.87
2962.32, 3226.25

Xiakou
1656.71, 2083.18
1592.84, 2156.09

The results for total crop areas and manure/fertilizer applications are listed in Table 5. The areas of citrus and tea would maintain the low levels over the planning horizon. This may be attributed to their pollutant losses as well as their low net energy and digestible protein contents (supplied for livestock). In period 1, wheat, potato, rapeseed, and alpine rice should be cultivated. The potato area (i.e.,  ha) would account for the largest one of the entire croplands. In period 2, wheat, potato, rapeseed, and alpine rice would be harvested and second rice, maize, and vegetable would be sown. The vegetable area (i.e.,  ha) would be the largest one among all croplands. The results indicate that the high levels of area planted with potato and vegetable are associated with their high crop yields, good market price, and low pollutant losses. In terms of manure and fertilizer application, their quantities would vary with crop areas. The results indicate that manure would be the main nitrogen and phosphorus source which satisfy the requirements of most crops due to its availability and low price to collect. Solutions demonstrate that domestic fowl husbandry would reach the highest level (i.e., unit) among all live stocks because it possesses more advantageous conditions than the other live stocks (i.e., higher allowable discharge, higher revenue parameter, and lower manure generation rate).

Table 5: Solutions for agricultural production.

Crops and live stocks
Growth period

Crop area (ha)

    Citrus
Whole year
35.7, 72.8

    Tea
Whole year
40.8, 53.4

    Wheat
Dry season
1112.6, 1894.2

    Potato
Dry season
1459.8, 1620.3

    Rapeseed
Dry season
296.5, 369.4

    Alpine rice
Dry season
206.5, 393.2

    Second rice
Wet season
351.2, 387.4

    Maize
Wet season
83.4, 97.8

    Vegetable
Wet season
1745.2, 1915.6

Manure application (10³ t)

    Citrus
Whole year
10.4, 19.3

    Tea
Whole year
8.2, 11.7

    Wheat
Dry season
2.9, 7.8

    Potato
Dry season
0.7, 2.8

    Rapeseed
Dry season
4.1, 4.7

    Alpine rice
Dry season
5.8, 6.3

    Second rice
Wet season
5.2, 5.7

    Maize
Wet season
0.8, 1.1

    Vegetable
Wet season
1.2, 1.4

Fertilizer application (t)

    Citrus
Whole year
0.7, 1.0

    Tea
Whole year
804.3, 893.6

    Wheat
Dry season
824.4, 917.3

    Potato
Dry season
1023.4, 1454.2

    Rapeseed
Dry season
2.4, 3.6

    Alpine rice
Dry season
39.6, 60.5

    Second rice
Wet season
292.7, 316.3

    Maize
Wet season
0.2, 0.3

    Vegetable
Wet season
1520.4, 1632.1

Size of livestock husbandry (10³ unit)

    Pig
Whole year
0.121, 35.462

    Ox
Whole year
0.688, 0.813

    Sheep
Whole year
21.430, 34.478

    Domestic fowl
Whole year
34.800, 132.561

Figures 2 and 3 present the amounts of BOD discharges from chemical plants and WTPs, respectively. The amount of BOD discharge is associated with a number of factors (e.g., production scale, wastewater generation rate, and wastewater treatment facility). The BOD discharge from BSH would be more than those from the other chemical plants, which would be  t (i.e., tonne) in period 1 and  t in period 2. The BOD discharge from GF would be the lowest among all chemical plants, with  t in period 1 and  t in period 2. Among all WTPs, Gufu wastewater treatment plant would discharge the highest BOD level, with  t in period 1 and  t in period 2; the BOD discharged from Nanyang WTP would be the lowest (i.e.  t in period 1 and  t in period 2).

Figure 2: BOD discharges from chemical plants. Note: the five chemical plants are abbreviated as their initials in this figure. GF: Gufu, ; BSH: Baishahe, ; PYK: Pingyikou, ; LCP: Liucaopo, ; XJLY: Xiangjinlianying, .

Figure 3: BOD discharges from WTPs.

Figure 4 shows TP discharges from point sources (i.e., chemical plants, WTPs, and phosphorus mining companies) and nonpoint sources (i.e., crop farming and agricultural life). In Figure 4, symbol “CP” denotes chemical plant; symbol “PMC” denotes phosphorus mining company; symbol “CF” denotes crop farming; symbol “AL” means agricultural life. For point sources, the chemical plants would be the major contributor to water pollution. The phosphorus pollutants can be discharged from wastewater and solid wastes (i.e., chemical wastes, slags, and tailings). The amount of TP discharge from chemical plant would be  t in period 1 and  t in period 2. Since the wastewater should be sluiced strictly according to the integrated discharge standards, the amount of TP discharge from WTPs would stay at a low level. Most of the TP would be from phosphorus-containing wastes (i.e., discharged directly and washed by rainfall). For nonpoint sources, the phosphorus pollutants from crop farming which can be generated though runoff and soil erosion (the latter would be a larger proportion) would be more than that from agricultural life. The amount of TP discharge from crop farming would be  t in period 1 and  t in period 2. This is associated with its high soil loss rate, low runoff, and low phosphorus concentration in the study area. Generally, TP discharge derives mainly from point sources, particularly from chemical plants. The results also indicate that the nitrogen pollutants would be generated by nonpoint sources (i.e., mainly from crop farming). TN discharge from citrus and tea can be neglected (i.e., nearly equal to aero). Figure 5 shows TN discharges from the other cropping areas. TN discharges from wheat and potato would be higher than those from rapeseed and alpine rice in period 1; TN discharges from vegetable would be higher than those from maize and second rice in period 2. This difference may be attributed to their planting areas, soil losses, runoff, and nitrogen concentration.

Figure 4: Solutions of TP discharge. Note: CP: chemical plant; WTP: wastewater treatment plant; PMC: phosphorus mining company; CF: crop farming; AL: agricultural life.

Figure 5: TN discharges from agriculture activities.

Net system benefit can be obtained from different industrial and agricultural activities, as shown in Figure 6. Chemical plants would be the major economic incoming source in the study area and could generate the highest revenue (RMB¥ ). Municipal water supply and phosphorus mining company would also make certain contribution to the economic development; their net benefits would be RMB¥ and RMB¥ , respectively. Livestock husbandry would bring the lowest benefit. Such an industry-oriented pattern may be related to the abundant mineral resources (particularly phosphorus ore) which can generate high economic return. Agricultural activities would make less contribution to the local economic development due to its topography which is not suitable for cultivation in large parts.

Figure 6: Net benefits from industrial and agricultural activities. Note: CP: chemical plant; MWS: municipal water supply; PMC: phosphorus mining company; CF: crop farming; LH: livestock husbandry.

5. Conclusions

In this study, an interval fuzzy credibility-constrained programming (IFCP) has been advanced for water quality management under uncertainty. This method integrates interval-parameter programming (IPP) and fuzzy credibility-constrained programming (FCP) techniques within a general optimization framework. Generally, the IFCP model has advantages in (1) handling uncertainties presented in terms of interval values and possibility distributions in the model, and (2) providing bases for determining optimal water quality management plans with desired compromises between economic benefits and environmental capacity-violation risks.

The developed model has been applied to a real-world case of planning water quality management in the Xiangxi River of the Three Gorges Reservoir Region. The objective is to maximize the net system benefit subject to the environmental requirements under uncertainty over the planning horizon. Pollutant discharges generated by various point and nonpoint sources were considered simultaneously. Interval solutions for production activities (i.e., industrial, municipal, and agricultural) and pollutant discharges (i.e., BOD, TP, and TN) under interval credibility levels have been generated by solving two deterministic submodels. The detailed results of related production scales and pollutant discharges can help identify desired water quality management schemes for developing the local economy in a sustainable manner. Some useful suggestions for the local economy development in a sustainable manner could be summarized: (i) advancing wastewater treatment technologies (e.g., tertiary treatment and depth processing technologies) to further improve pollutant removal efficiency; (ii) controlling the generation of phosphorus-containing wastes (from chemical plants and phosphorus mining companies) strictly in the production process and taking effective treatments and disposal measures to reach the goal of achieving TP abatement; (iii) taking control practices on soil erosion for reducing the transport of nitrogen and phosphorus pollutants to the river.

Although reasonable solutions and desired management policies have been obtained through the IFCP management model, there are still some extensive research works to be done. For example, the proposed IFCP method can deal with uncertainties expressed as fuzzy sets and interval numbers; however, the main limitations of the IFCP method remain in its difficulties in tackling uncertainties expressed as probabilistic distributions (stochastic uncertainties). Under such a circumstance, stochastic mathematical programming method is a suitable option to be introduced into the proposed IFCP method. Moreover, decision support regarding pollution management could be further provided by incorporating certain water quality simulation models into IFCP framework, which can effectively reflect dynamic interactions between pollutant loading and water quality.

Acknowledgments

This research was supported by the Natural Sciences Foundation of China (Grant nos. 51379075, 51225904, and 51109077), the 111 Project (B14008), and the Program for Innovative Research Team in University (IRT1127). The authors are grateful to the editors and the anonymous reviewers for their insightful comments and suggestions.

References

P. Du, Y. P. Li, and G. H. Huang, “Inexact chance-constrained waste-load allocation model for water quality management of Xiangxihe River,” Journal of Environmental Engineering, vol. 139, pp. 1178–1197, 2013.

Y. P. Li and G. H. Huang, “Two-stage planning for sustainable water-quality management under uncertainty,” Journal of Environmental Management, vol. 90, no. 8, pp. 2402–2413, 2009.

Y. P. Li, W. Li, and G. H. Huang, “Two-stage inexact-probabilistic programming model for water quality management,” Environmental Engineering Science, vol. 29, no. 7, pp. 713–725, 2012.

Y. P. Li and G. H. Huang, “Fuzzy-stochastic-based violation analysis method for planning water resources management systems with uncertain information,” Information Sciences, vol. 179, no. 24, pp. 4261–4276, 2009.

X. S. Qin, G. H. Huang, G. M. Zeng, A. Chakma, and Y. F. Huang, “An interval-parameter fuzzy nonlinear optimization model for stream water quality management under uncertainty,” European Journal of Operational Research, vol. 180, no. 3, pp. 1331–1357, 2007.

B. Julien, “Water quality management with imprecise information,” European Journal of Operational Research, vol. 76, no. 1, pp. 15–27, 1994.

C.-S. Lee and C.-G. Wen, “Fuzzy goal programming approach for water quality management in a river basin,” Fuzzy Sets and Systems, vol. 89, no. 2, pp. 181–192, 1997.

P. P. Mujumdar and K. Sasikumar, “A fuzzy risk approach for seasonal water quality management of a river system,” Water Resources Research, vol. 38, no. 1, pp. 51–59, 2002.

H.-W. Chen and N.-B. Chang, “Decision support for allocation of watershed pollution load using grey fuzzy multiobjective programming,” Journal of the American Water Resources Association, vol. 42, no. 3, pp. 725–745, 2006.

S. Karmakar and P. P. Mujumdar, “Grey fuzzy optimization model for water quality management of a river system,” Advances in Water Resources, vol. 29, no. 7, pp. 1088–1105, 2006.

Y. P. Li, G. H. Huang, S. L. Nie, and D. W. Mo, “Interval-parameter robust quadratic programming for water quality management under uncertainty,” Engineering Optimization, vol. 40, no. 7, pp. 613–635, 2008.

X. H. Nie, G. H. Huang, D. Wang, and H. L. Li, “Robust optimisation for inexact water quality management under uncertainty,” Civil Engineering and Environmental Systems, vol. 25, no. 2, pp. 167–184, 2008.

S. Maeda, T. Kawachi, K. Unami, J. Takeuchi, T. Izumi, and S. Chono, “Fuzzy optimization model for integrated management of total nitrogen loads from distributed point and nonpoint sources in watershed,” Paddy and Water Environment, vol. 7, no. 3, pp. 163–175, 2009.

T. Y. Xu and X. S. Qin, “Solving water quality management problem through combined genetic algorithm and fuzzy simulation,” Journal of Environmental Informatics, vol. 22, no. 1, pp. 39–48, 2013.

M. Liu, G. X. Nie, M. Hu, R. F. Liao, and Y. S. Shen, “An interval-parameter fuzzy robust nonlinear programming model for water quality management,” Journal of Water Resource and Protection, vol. 5, pp. 12–16, 2013.

Y. M. Zhang and G. H. Huang, “Optimal water resource planning under fixed budget by interval-parameter credibility constrained programming,” Engineering Optimization, vol. 43, no. 8, pp. 879–889, 2011.

M. S. Pishvaee, S. A. Torabi, and J. Razmi, “Credibility-based fuzzy mathematical programming model for green logistics design under uncertainty,” Computers and Industrial Engineering, vol. 62, no. 2, pp. 624–632, 2012.

Y. P. Li, G. H. Huang, and S. L. Nie, “Optimization of regional economic and environmental systems under fuzzy and random uncertainties,” Journal of Environmental Management, vol. 92, no. 8, pp. 2010–2020, 2011.

R. Brouwer and C. De Blois, “Integrated modelling of risk and uncertainty underlying the cost and effectiveness of water quality measures,” Environmental Modelling and Software, vol. 23, no. 7, pp. 922–937, 2008.

Y. R. Fan and G. H. Huang, “A robust two-step method for solving interval linear programming problems within an environmental management context,” Journal of Environmental Informatics, vol. 19, no. 1, pp. 1–9, 2012.

Y. Xu, G. H. Huang, and T. Y. Xu, “Inexact management modeling for urban water supply systems,” Journal of Environmental Informatics, vol. 20, no. 1, pp. 34–43, 2012.

Y. Chen, S. Y. Cheng, L. Liu et al., “Assessing the effects of land use changes on non-point source pollution reduction for the Three Gorges Watershed using the SWAT model,” Journal of Environmental Informatics, vol. 22, no. 1, pp. 13–26, 2013.

L. He, G. Q. Wang, and C. Zhang, “Application of loosely coupled watershed model and channel model in Yellow River, China,” Journal of Environmental Informatics, vol. 19, no. 1, pp. 30–37, 2012.

Y. P. Li and G. H. Huang, “A recourse-based nonlinear programming model for stream water quality management,” Stochastic Environmental Research and Risk Assessment, vol. 26, no. 2, pp. 207–223, 2012.

L. Ye, Q.-H. Cai, R.-Q. Liu, and M. Cao, “The influence of topography and land use on water quality of Xiangxi River in three Gorges Reservoir region,” Environmental Geology, vol. 58, no. 5, pp. 937–942, 2009.

F. A. Deviney Jr., D. E. Brown, and K. C. Rice, “Evaluation of bayesian estimation of a hidden continuous-time markov chain model with application to threshold violation in water-quality indicators,” Journal of Environmental Informatics, vol. 19, no. 2, pp. 70–78, 2012.