Scientific Programming

Volume 2017, Article ID 1680813, 14 pages

https://doi.org/10.1155/2017/1680813

## Inexact Multistage Stochastic Chance Constrained Programming Model for Water Resources Management under Uncertainties

^{1}College of Water Conservancy and Hydropower, School of Science, Hebei University of Engineering, Handan 056038, China^{2}College of Arts, Hebei University of Engineering, Handan 056038, China^{3}School of Economics and Management, Handan University, Handan 056038, China

Correspondence should be addressed to Minghu Ha; nc.ude.uebeh@uhgnimah

Received 25 December 2016; Accepted 23 May 2017; Published 20 June 2017

Academic Editor: Fabrizio Riguzzi

Copyright © 2017 Hong Zhang 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 order to formulate water allocation schemes under uncertainties in the water resources management systems, an inexact multistage stochastic chance constrained programming (IMSCCP) model is proposed. The model integrates stochastic chance constrained programming, multistage stochastic programming, and inexact stochastic programming within a general optimization framework to handle the uncertainties occurring in both constraints and objective. These uncertainties are expressed as probability distributions, interval with multiply distributed stochastic boundaries, dynamic features of the long-term water allocation plans, and so on. Compared with the existing inexact multistage stochastic programming, the IMSCCP can be used to assess more system risks and handle more complicated uncertainties in water resources management systems. The IMSCCP model is applied to a hypothetical case study of water resources management. In order to construct an approximate solution for the model, a hybrid algorithm, which incorporates stochastic simulation, back propagation neural network, and genetic algorithm, is proposed. The results show that the optimal value represents the maximal net system benefit achieved with a given confidence level under chance constraints, and the solutions provide optimal water allocation schemes to multiple users over a multiperiod planning horizon.

#### 1. Introduction

With the economic development and population growth, the demands for water resources increase every year. Due to unreasonable exploitation and utilization of water resources, water pollution, and extreme weather conditions, more and more countries and regions are faced with different degrees of water shortage. Making the best use of limited water resources premises the sustainable development of economy and society [1]. In the past few decades, many researchers have applied optimization techniques [2–7] to deal with uncertainties in a number of system components and their interrelationships within water resources systems [8]. Among them, inexact multistage stochastic programming (IMSP) is regarded as a significant method for water resources management. For example, Li et al. [9] proposed an interval-parameter multistage stochastic linear programming model which incorporates inexact optimization and multistage stochastic linear programming to manage water resources. Zhou et al. [10] developed a factorial multistage stochastic programming which is a hybrid methodology of factorial analysis and IMSP to analyze the potential interrelationships among a variety of uncertain parameters and their impacts on system performance for water resources management. Suo et al. [11] proposed an inventory-theory-based inexact multistage stochastic programming for water resources management through introducing an inventory theory into the framework of IMSP to provide reasonable transferring schemes associated with various flow scenarios. IMSP was also combined with other types of uncertainties, such as fuzziness [12–14], to solve the water resources management problems. However, the above models cannot assess the risks in water resources management systems, to which chance constrained programming (CCP) is an effective method to resort to.

CCP, pioneered by Charnes and Cooper [15], provides a method to handle uncertainties by specifying a confidence level at which it is desired that the constraint holds [16], which has been applied in several areas [17, 18]. In water resources management, using CCP is conducive to not only making decisions of water allocation but also gaining insights into the tradeoffs between economic objectives and some policy factors for a water resource manager [19]. For instance, Liu et al. [1] provided a factorial multistage stochastic programming with chance constraints to deal with the issues of constraint-violation risks and interactive uncertainties in water resources management systems. However, this kind of CCP models can only be solved in some special cases. Meanwhile, the models should firstly be converted into deterministic equivalents and then be solved by using solution methods of deterministic mathematical programming [16] in the solution process. Complex chance constrained programming models are hard to solve, such as nonlinear models which contain multiply distributed data [20]. In order to overcome the limitations, Liu [21] presented a framework of CCP with the assumption that not only uncertain constraints but also uncertain objective would hold at different confidence levels provided by the decision-maker as an appropriate safety margin and gave a hybrid algorithm to solve the CCP for more cases. It is noticeable that the choice of confidence levels makes the model flexible to more situations. The CCP has been applied to many aspects of system decision-making, such as project scheduling [22], capital budgeting [23], capacitated location-allocation [24], and redundancy-optimization [25]. Up to now, Liu [21]’s CCP has not been applied in the optimization of water resources management although this method is needed because of the complexity of the water resources management systems.

In our previous relevant research, [26] extended Liu [21]’s CCP from probability space to Sugeno measure space, in which Sugeno measure was one of representative nonadditive measures. Although the theoretical basis of the CCP was discussed in detail, we did not apply the model to water resources management. Paper [27], which managed water supply risk by constructing a model of water option pricing, did not refer to water resources allocation and assess the risks by using CCP. As an extension of the previous efforts, this study aims to develop an inexact multistage stochastic chance constrained programming (IMSCCP) model, which incorporates the CCP proposed by Liu [21] into the IMSP framework for water resources management. The model can tackle uncertainties in the objective and constraints which present themselves as probability distributions and interval with multiply distributed stochastic boundaries, analyze various policy scenarios when the promised policy targets are violated, and deal with the issue of risks presented as stochastic constraints and stochastic objective with predetermined confidence levels. A hypothetical case study of water resources management within three planning periods is given to demonstrate the applicability of the method. Moreover, a hybrid algorithm incorporating stochastic simulation, back propagation (BP) neural network, and genetic algorithm (GA) is proposed to solve the model. The developed method obtains the results of the model in which the optimal value represents the maximal net system benefit achieved with a given confidence level subject to chance constrains.

#### 2. Methodologies

##### 2.1. Multistage Inexact Stochastic Programming Model

It is important for a water resources manager to allocate water to multiple users such as municipal, industrial, and agricultural sectors from an unregulated reservoir over a multiperiod planning horizon in an optimized way, which could be formulated as maximizing the expected economic revenue based on the water allocation in the region over the planning horizon [28]. The manager promises a quantity of water to each user. If the promised water is delivered, net benefits will be brought to local economy; otherwise, the deficient water must be obtained from alternative and more expensive sources or the demand must be curtailed, resulting in economic penalties on local economy [29]. Considering the randomness of the water flow and the dynamic feature of the long-term water allocation plans, the problem of water allocation can be formulated as the following scenario-based multistage stochastic programming model for water resources management under uncertainties [9]: where is the net benefit of the water allocation system ($); is the net benefit when per unit () of water is allocated to user in period ($/); is the fixed amount of water allocation target promised to user in period (); is the water flow in period which is a random variable (); is the amount of water shortage to user when the seasonal flow is (); is the loss when per unit () of water is not allocated to user in period ($/); is the expected value of a random variable; is the surplus water inflow in period () according to ; is the amount of maximum allowable allocation for user (); is the number of water users; is the index of water user for , with for the municipality, for the industrial user, and for the agricultural sector. In the model, is the decision variable and is the objective function. It is observed that model (1) reflects nonanticipativity, since a decision must be made in each stage without the knowledge of realizations of random variables in future stages.

Let take values of with probability for scenarios in each planning period (), where is the probability of occurrence for scenario in period and . Then, we have ,, where is the sum of scenarios in period ; is the amount of water allocation shortage to user when the scenario occurs in period . Then model (1) can be reformulated as follows [9]:

Considering that the uncertainties exist in variables and coefficients, the fixed values of the parameters, such as , , , and , cannot be determined exactly. Thus, the inexact multistage stochastic programming (IMSP) model, which introduces interval parameters into model (2), is proposed as follows [9]:

##### 2.2. Inexact Multistage Stochastic Chance Constrained Programming Model

Models (1), (2), and (3) do not readily assess the risks, and they only deal with uncertainties in the right hand side such as the water flow . It is difficult to handle uncertainties in both the left and right hand sides (i.e., , , and ) [1] which are presented as interval with stochastic normal distributed boundaries. In view of the above considerations, Liu et al. [1] combined chance constrained programming (CCP, initiated by Charnes and Cooper [15]) with IMSP to propose the following inexact multistage stochastic programming model with chance constraints to solve problems with the request that chance constraints should hold at least with prescribed levels of probability (i.e., confidence levels):

Model (4) can handle uncertainties presented as interval with normally distributed boundaries. However, uncertainties in the water resources management systems have more diverse forms of performance besides normal distribution, which suggests a need for models that can handle uncertainties presented as interval with multiple distributed boundaries. And the confidence levels only in constraints may be not enough to indicate the relationship between the economic objective and the system risk. Moreover, a water resources manager may want to obtain the maximum that the objective function achieves with a given confidence level subject to stochastic constraints with other confidence levels, which means that the confidence levels are not only in constraints but also in objective.

Liu [21] proposed a framework of nonlinear chance constrained programming with confidence levels occurring in constraints and objective and provided a stochastic simulation based genetic algorithm to solve the CCP. This CCP model can be formulated as follows:where is a decision vector; is a random vector; and are predetermined confidence levels for stochastic objective and stochastic constraint(s), respectively; denotes the probability of the event in . This programming aims to maximize that the objective function achieves with at least probability of ( is called* α*-optimistic value to ). is feasible if and only if the probability measure of the set is at least for .

By incorporating CCP [21] and multistage stochastic programming (model (2)), a multistage stochastic chance constrained programming model for water resources management could be formulated as follows:where and are predetermined confidence levels for stochastic objective and stochastic constraint(s), respectively. This programming aims to obtain the *α*-optimistic value to the objective function of . Obviously, model (6) gives predetermined confidence levels for constraints and objective, which indicates more comprehensive risk assessment in the water resources management systems.

Considering that the fixed values of the parameters in model (6) cannot be determined exactly, an inexact multistage stochastic chance constrained programming (IMSCCP) model is proposed as follows:

Let , where , , . Model (7) could be replaced by the following form:

Next, we can change model (8) into two submodels corresponding to the upper and lower bounds for the desired objective function value.

Firstly, we formulate submodel (9) corresponding to the upper bound of the objective function value; that is,where and are the decision variables and is the -optimistic value to the net system benefit . In submodel (9), the optimal solutions are denoted by and , and the optimal value is denoted by .

Secondly, submodel (10) corresponding to the lower bound of the objective function value can be formulated as follows: where is the decision variable and is the -optimistic value to the net system benefit . In submodel (10), the optimal solution is denoted by and the optimal value is denote by .

Thirdly, the real allocation of water to user when the scenario occurs in period can be calculated by .

Finally, the optimized interval solution for the decision variable and the real interval allocation of water can be obtained, while the optimized interval objective value can be generated through the variation of .

#### 3. Case Study

In this section, the IMSCCP model is applied to the water resources management systems. Just as the statement for the hypothetical problem described in Huang and Loucks [30], a water resources manager shoulders the responsibility of delivering water from an unregulated reservoir to three sectors, that is, municipality, industry, and agriculture, during three periods. All users want to know allocated water amount that they can expect over the three periods. If the allocated water satisfies the demand of the user, per unit water will gain net benefit for the local economy. Otherwise, the user will spend more to obtain water from other reservoirs or curtail their expansion plans, which means per unit water deficiency will gain penalty. Table 1 provides the water allocation targets in the three planning periods. Tables 2 and 3 present the distribution of the stream flows in the three periods. Table 4 shows the related economic data. Obviously, the boundaries of the data in Tables 3 and 4 are random variables satisfying different distribution forms. The objective is to obtain 90%-optimistic value to the net system benefit while the constraints of water availability hold with at least a probability of 95% over the planning horizon.