# Multiobjective Optimization Problem of Multireservoir System in Semiarid Areas

Authors: Z. J. Chen; Z. J. Cheng; X. Q. Yan
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

With the increasing scarcity of water resources, the growing importance of the optimization operation of the multireservoir system in water resources development, utilization, and management is increasingly evident. Some of the existing optimization methods are inadequate in applicability and effectiveness. Therefore, we need further research in how to enhance the applicability and effectiveness of the algorithm. On the basis of the research of the multireservoir system’s operating parameters in the Urumqi River basin, we establish a multiobjective optimization problem (MOP) model of water resources development, which meets the requirements of water resources development. In the mathematical model, the domestic water consumption is the biggest, the production of industry and agricultural is the largest, the gross output value of industry and agricultural is the highest, and the investment of the water development is the minimum. We use the weighted variable-step shuffled frog leaping algorithm (SFLA) to resolve it, which satisfies the constraints. Through establishing the test function and performance metrics, we deduce the evolutionary algorithms, which suit for solving MOP of the scheduling, and realize the multiobjective optimization of the multireservoir system. After that, using the fuzzy theory, we convert the competitive multiobjective function into single objective problem of maximum satisfaction, which is the only solution. A feasible solution is provided to resolve the multiobjective scheduling optimization of multireservoir system in the Urumqi River basin. It is the significance of the layout of production, the regional protection of ecological environment, and the sufficient and rational use of natural resources, in Urumqi and the surrounding areas.


---

## Body

## 1. Overview

Xinjiang is the most north-west province in China. It is located in the inland semiarid area, which accounts for 3% of water resources in China. The average annual surface water runoff in Xinjiang is 88.2 billion cubic meters and 5146 cubic meters per capita, which is 2.25 times the average of China [1]. The average annual groundwater exploitation is about 25.1 billion cubic meters. In addition, the glaciers reserves in Xinjiang accounted for 50 percent of China’s total. Yet, Xinjiang is located in the hinterland of the Eurasian continent; the climate is dry; it is affected by seasonal factors easily. So, the spatial and temporal distribution of the water resource is extremely uneven [2].As the capital and the center of the cultural, economic of the Xinjiang Uygur autonomous region, the population of Urumqi city is nearly 5 million. The urban water resources are mainly provided by the Urumqi River basin, which belongs to the northern slope of Tianshan inland river system. It consists of the Urumqi River, Nanshan River, Dongshan River, Toutunhe River, and so forth. Urumqi River originates in the Glacier number 1 of Tianger peak, crosses the urban area of Urumqi city, and, finally, flows into Beishawo in the southern margin of Junggar basin. Its length is 214 km; drainage area is 5803 square kilometers; average annual runoff is 2.37 million cubic meters. Along the Urumqi River basin, there are more than 30 reservoirs as follows: Daxigou reservoir, Ulabo reservoir, Hongyanchi reservoir, Mengjin reservoir, and so on. It provides domestic water, farmland and green space irrigation water, industrial water, and aquaculture water for the residents in Urumqi city and the surrounding area. Meanwhile, it also plays an important role in flood control, drought resistance, and others [3].In inland arid and semiarid regions, such as Urumqi, the shortage of water resources has become a major restricting factor of the city’s development. With the rapid development of the city as well as the gradual deterioration of the ecological environment, severe water shortage has become the bottleneck of the constraints and limitations of social and economic development in Urumqi. In order to alleviate this problem, in the beginning of this century, “YinEJiU” project was implemented. Using hydraulic engineering, part of northern Xinjiang Irtysh water was injected into the Urumqi River. With the implementation of this project, a new water source is added into the Urumqi River basin. It alleviated the status of the serious water shortage in Urumqi to a certain extent. However, due to the limited increase of the water resources, it cannot fully meet the demand. Therefore, from the management point of view, in addition to the plan of reservoirs construction in the upstream, we can also apply the multiobjective optimization operation. In this way, we can alleviate the contradiction between the supply and the demand of the water resource at maximum and realize the optimal allocation of water resources in Urumqi river basin gradually.Since the early 1960s, the research on the optimization operation of the multireservoir system has been carried out in China. The research developed rapidly after the 1980s and got a large amount of achievements. Dong introduced the principle and contents of the dynamic programming optimization technique and its application in the hydro power station reservoir (group) scheduling, as well as the principle of the multiobjective decision. Feng established the theoretical foundation for the multiobjective optimum scheduling of the reservoirs. He also discussed the basic principle and method of the large-scale system decomposition and coordination systematically. Wang introduced the theory and application of the large-scale system optimization. Using the fuzzy identification mode, Li achieved the level characteristic​​of the sample and the contribution rate of the relevant factors in the evaluation process of the evaluation, which aims at the comprehensive evaluation of the reservoir water resources carrying capacity. At the same time, he achieved assessment of water resources carrying capacity in regional reservoir [4].In the recent decades, several researchers (Houck, Yakowitz, and Needham et al.) proposed methods applied in optimization operation of multireservoir system, such as linear programming, nonlinear programming, and dynamic planning. Although those optimization methods provide a way to resolve the problem of reservoir group scheduling, there are also some problems. First, some models need a lot to simplify; therefore, their simulation performance is poor. Second, some of the nonlinear methods are restricted in practical application because of the lack of the boundary conditions and more parameters. Moreover, those methods are hard to simulate the experience of scheduling knowledge. Consequently, the research on the optimization scheduling is more, and the practical application is less [5].In the process of the multireservoir system operation, the regulation and scheduling consist of the conflicting and the intereffecting multitargets. In a specific field, people often encounter the matching problem which requires multiple solutions to achieve optimization at the same time. It is the multiobjective optimization problem, MOP. The validation of MOP is whether the optimization target is more than one and requires the simultaneous processing. In engineering applications, the MOP is very significant. These engineering applications are very complicated and difficult. They are also one of the main research fields. Since the early nineteen sixties, the MOP has attracted more and more researchers' attention from different fields. For this reason, it has very important research value and practical significance to solve MOP [6].There are a large number of monographs focusing on the methods to resolve the MOP. In general, these books can be divided into two categories. One uses the classical optimization methods, and the other uses the methods based on evolutionary algorithm which simulates the biological evolution or activity behavior. Because evolutionary algorithm’s constraint conditions are low and solution quality is high, they have been widely applied in resolving the complex optimization problem. At present, the popular evolutionary algorithm to resolve the MOP is as follows: the multiobjective genetic algorithm [7], the artificial ant colony algorithm [8], the multiobjective particle swarm optimization algorithm [9], and the shuffled frog leaping algorithm [10].The shuffled frog leaping algorithm (SFLA) is a biological evolutionary algorithm based on swarm intelligence, which was proposed by Eusuff and Lansey in 2003. The SFLA has the characteristics of simple concept, fewer parameters, fast calculation speed, good global search ability, easy to implement, and so on. So at present, it is wildly applied to the fields of the network of water resources allocation issues, the function optimization, the multiproduct pipeline network optimization, the combinatorial optimization, the image processing, the multiuser detection, and so forth. Also, the method achieved good results. As a kind of optimization method, the SHLA has certain advantages in the optimization operation of the multireservoir system and has broad application prospects [11–13].The purpose of this paper is to study the operation parameters. Based on these parameters, we set up the operation mathematical model and use the SFLA to realize the optimization allocation of the multireservoir system in Urumqi River basin [14].


## 2. The Description of the MOP

Using words, the MOP can be described as an optimization problem, which consists ofD decision variable parameters, N objective functions, and m+n constraints. The decision variables, the objective function, and the constraint condition are function relation. In noninferior solution set, decision makers can only choose a satisfactory noninferior solution as the final solution according to the specific requirements of the problem. The mathematic form of MOP can be described as follows:
(1)min(y)=f(x)=[f1(x),f2(x),…,fn(x)],min(y)=f(x)=min(y)=n=1,2,…,Ns.t.gi(x)≤0,i=1,2,…,mhj(x)=0,j=1,2,…,kx=[x1,x2,…,xd,…,xD]xd-min≤xd≤xd-max,d=1,2,…,D,
where x is the decision-making relative to the amount of the D-dimensional, y is the target vector, g(x)≤0 is the ith inequality constraint, hj(x)=0 is the jth equality constraints, fn(x) is the nth objective function, X the decision space formed by the decision vector, gi(x)≤0 and hj(x)=0  determine the feasible solution domain, and xd-min and xd-max are upper and lower searching limits of each dimension vector [15].For optimal solution or noninferior optimal solution in the MOP, we can define the following.Definition 1.
For anyd∈[1,D] to meet xd*≤xd and exits d0∈[1,D], there is xd*≤xd0, then the vector x*=[x1*,x2*,…,xD*] dominates the vector x=[x1,x2,…,xD].
Iff(x*) dominates f(x), the following two conditions must be met:
(2)∀n,fn(x*)≤fn(x),n=1,2,…,N,∃n0,fn0(x*)<fn0(x),1≤n0≤N.
The dominating relation off(x) is consistent with the dominating relation of x.Definition 2.
The Pareto optimal solution cannot be dominated by any solution in the feasible solution set. Ifx* is a point in the searching space, x* is the noninferior optimal solution. If and only if exit x (in the feasible domain of the searching space), it makes the establishment of (x*)<f(x), n=1,2…,N.Definition 3.
Given a MOP off(x), if and only if for any x (in the searching space), there is fn(x)<fn(x*); f(x*) is the global optimal solution.Definition 4.
The solution set consisted of all the noninferior optimal solutions is called Pareto optimal set of the MOP, also called acceptable solution set or efficient solution set.


## 3. The Description of the Shuffled Frog Leaping Algorithm (SFLA)

The SFLA is a kind of cooperative evolutionary algorithm enlightened by natural biological mimics, which simulates the information exchange classified by subgroup in the process of the frog group looking for food. First, we generate a group of initial solution (population) from solution space randomly. Then, we divide the whole population into more than one subpopulation, in which the frogs execute subpopulation internal search according to a certain strategy. After the end of the number of searches in defined subpopulation, mix all the frogs, sort and divide the subpopulation again, and exchange information between each subpopulation globally. The sub-population internal search and the global information exchange continue alternating, until they satisfy the convergence condition or reach the maximum number of evolution. The following is mathematical model of the SHLA.


### 3.1. Dividing Subpopulation

Suppose that the number of frogs in the population (candidate solution) isN, the number of subpopulation is k, the number of candidate solutions is n. Xj=(xj1,xj2,…,xjD) represents the jth candidate solution, j(0≤j≤N), and D represents the dimensions of the candidate solution. For the initial population S=(X1,X2,…,XN) generated randomly, after sorted according to the fitness f(x) in descending, the first candidate solution is assigned to the first subpopulation, the second candidate solution is assigned to the second subpopulation,…, the kth candidate solution is assigned to the kth subpopulation, the (k+1)th candidate solution is assigned to the first subpopulation, and the (k+2)th candidate solution is assigned to the second subpopulation, repeat until N candidate solutions are distributed.


### 3.2. Internal Search in Subpopulation

Suppose thatXb is the candidate solution, in which fitness is best in the subpopulation, Xw is the candidate solution, in which fitness is worst in the subpopulation, and Xm is the candidate solution, in which fitness is best in the whole population. For every subpopulation, we do internal search. It is to update Xw of the subpopulation. The search strategy is as follows:
(3)X′=Xw+R×(Xb-Xw),
where X′ is the new solution produced by (3) and R is a random number between 0 and 1. If X′ is superior to Xw, then Xw=X′. Otherwise, Xb is replaced by Xm in (3) and repeats the search strategy. If X′ is not always superior to Xw, a random candidate solution is generated to replace Xw. We repeat the above steps until the number of searches is greater than the largest number of the largest subpopulation internal search.


### 3.3. Global Information Exchange

When all the internal updates of subpopulation are completed, we do the subpopulation dividing and internal search again, until it meet the termination conditions (converges to the optimal solution or reaches the maximum number of evolving generation).


## 3.1. Dividing Subpopulation

Suppose that the number of frogs in the population (candidate solution) isN, the number of subpopulation is k, the number of candidate solutions is n. Xj=(xj1,xj2,…,xjD) represents the jth candidate solution, j(0≤j≤N), and D represents the dimensions of the candidate solution. For the initial population S=(X1,X2,…,XN) generated randomly, after sorted according to the fitness f(x) in descending, the first candidate solution is assigned to the first subpopulation, the second candidate solution is assigned to the second subpopulation,…, the kth candidate solution is assigned to the kth subpopulation, the (k+1)th candidate solution is assigned to the first subpopulation, and the (k+2)th candidate solution is assigned to the second subpopulation, repeat until N candidate solutions are distributed.


## 3.2. Internal Search in Subpopulation

Suppose thatXb is the candidate solution, in which fitness is best in the subpopulation, Xw is the candidate solution, in which fitness is worst in the subpopulation, and Xm is the candidate solution, in which fitness is best in the whole population. For every subpopulation, we do internal search. It is to update Xw of the subpopulation. The search strategy is as follows:
(3)X′=Xw+R×(Xb-Xw),
where X′ is the new solution produced by (3) and R is a random number between 0 and 1. If X′ is superior to Xw, then Xw=X′. Otherwise, Xb is replaced by Xm in (3) and repeats the search strategy. If X′ is not always superior to Xw, a random candidate solution is generated to replace Xw. We repeat the above steps until the number of searches is greater than the largest number of the largest subpopulation internal search.


## 3.3. Global Information Exchange

When all the internal updates of subpopulation are completed, we do the subpopulation dividing and internal search again, until it meet the termination conditions (converges to the optimal solution or reaches the maximum number of evolving generation).


## 4. Reservoir Group Optimization Dispatching of Urumqi River Basin

In the process of multireservoir system optimization operation of Urumqi River basin, we should consider the following five aspects of the requirements: first, to satisfy the demand for residential water; second, to satisfy the water demand of the industrial and the agricultural production; third, to satisfy the water demand of the ecological environment; fourth, to complete the industrial and the agricultural planning output target; last, the investment of water resources development is minimum, under the premise to satisfy the demand of water. We consider these five aspects, which are mutual compensation, to achieve the most satisfactory results, under the condition that satisfies the system performance [16].


### 4.1. The Objective Functions of Multireservoir System Optimization Operation

With the above, we take the following as the goal function: people living water consumption is the largest, the industrial and the agricultural production water consumption is the largest, ecological environment water consumption is the largest, the gross output value of industry and agriculture is the highest, and the investment of water resources development is minimum [17].According to the MOP, the objective function, the decision variables, and the constraints of the reservoir group optimization dispatching can be described as follows:(4)min/maxfm(x1,x2,…,xn),m=1,2,…,M;s.t.(1)gk(x1,x2,…,xn)≥0,k=1,2,…,M;(2)hl(x1,x2,…,xn)=0,l=1,2,…,L;(3)xi(L)≤xi≤xi(U),i=1,2,…,n,
where X⊂Rn is the decision space, x=(x1,x2,…,xn)∈X, and xi(L) and xi(H) are upper and lower bounds of variables, respectively. Further, we can define collection Ω as the feasible domain of (4):
(5)Ω={xi(L)x∣gk(x1,…,xn)≥0,hk(x1,…,xn)=0,Ωxi(L)≤xi≤xi(U)},Ω⊂X.Among the objective functions of multireservoir system optimization operation, the investment of water resources development is the minimizing objective function, and the others are the maximum objective functions. In order to maintain unity, by taking the negative of the water resources investment, we convert all the problems into the maximum objective functions MOP. The objective functions can be described as follows:(6)maxf1=∑t=1T∑iIa1iPGi(t)2+β1iPGi(t)+γ1i,maxf2=∑t=1T∑jJa2jPSj(t)2+β2jPSj(t)+γ2j,maxf3=∑t=1T∑kKa3kPWk(t)2+β3kPWk(t)+γ3k,maxf4=∑t=1T∑lLa4lPHl(t)2+β4lPHl(t)+γ4l,maxf5=∑t=1T∑mMa5mPUm(t)2+β5mPUm(t)+γ5m.For (6), I, J, K, L, and M denote the effects of five optimal goal influence factors, respectively. PGi(t) is the water consumption of the ith living water user in the tth time interval; PSi(t) is the water consumption of the jth industry and agriculture water user in the tth time interval; PWk(t) is the water consumption of the kth ecological water user in the tth time interval; PHl(t) is the output value of the lth production unit in the tth time interval; PUm(t) is the investment amount of the mth water resources development unit in the tth time interval. an*,βn*,γn*(n=1,2,3,4,5;*=i,j,k,l,m) is the influence coefficient of each factor; T is the scheduling cycle.


### 4.2. The Constraint Condition of Multireservoir System Optimization Operation

In dealing with the constraint condition, according to the objective function, we make multireservoir system operation satisfaction as constraint condition. In view of reservoir group optimal model, the constraint conditions are as follows.The water consumption demand of people living is mainly affected by the climate change and daily change. So, the living water consumption change and its rate of change are always in a certain range. As one of the important indexes which influence the reservoir group operation, under the premise of the whole objective optimization and reducing the consumption of water, as far as possible, the constraint conditions of the living water consumption can be described as follows:(7)PGmin<∑PG<PGmax,Δ∑PG<ΔPGmax,
where ∑PG is the living water consumption of the reservoir group system in unit interval (usually a day). Δ∑PG  is the variation of the living water consumption of the reservoir group system in unit interval. PGmin and PGmax are the minimum and the maximum living water consumption, respectively. ΔPGmax  is the largest living water quantity variation in unit interval [18].The industrial and agricultural water consumption is directly related to the production situation. In addition to this, the agricultural water consumption is also affected by the weather, the rainfall, the plant growth period, and other factors. The industrial water consumption is also affected by the environment temperature. As another important indicator of the reservoir group operation satisfaction, the constraint conditions of the industrial and agricultural water consumption can be described as follows:(8)PSmin<∑PS<PSmax,Δ∑PS<ΔPSmax,
where ∑PS is the industrial and the agricultural water consumption of the reservoir group system in unit interval (usually a day). Δ∑PS  is the variation of the industrial and the agricultural water consumption of the reservoir group system in unit interval. PSmin and PSmax are the minimum and the maximum industrial and agricultural water consumption, respectively. ΔPSmax is the largest industrial and agricultural water quantity variation in unit interval [19].The ecological environmental water consumption is the other important satisfaction index of the reservoir group operation. The constraint conditions of the ecological environmental water consumption can be described as follows:(9)PWmin<∑PW<PWmax,Δ∑PW<ΔPWmax,
where ∑PW is the ecological environmental water consumption of the reservoir group system within the unit interval (usually a day). Δ∑PW is the variation of the ecological environmental water consumption of the reservoir group system within the unit interval. PWmin and PWmax are the minimum and the maximum industrial and agricultural water consumption, respectively. ΔPWmax is the largest ecological environmental water quantity variation within the unit interval. The ecological environmental water consumption is influenced by the climate, environment, and external factors, and it relates to the living water consumption, the industrial, and agricultural water consumption.The gross output value of industry and agriculture is one of the most important satisfaction indexes of the reservoir group operation. We can describe it with the social gross output. So its constraint conditions can be described as follows:(10)∑PH≥PHmin,
where ∑PH is the gross output value of industry and agriculture in a unit of time (usually a month), and ΔPHmin is the minimum gross output value in unit interval.It is the necessary conditions for the investment of the water resources development to ensure the normal operation of the reservoir group. Under the general condition, the minimum investment maximizes the benefits. So, the constraint conditions of the investment of the water resources development can be described as follows:(11)∑PU≤PUmax,∑PUCount(PU)≥P¯Umin,
where ∑PU is the investment of the water resources development, PUmax is the maximum investment of the water resources development, Count(PU) is the water resources investment project number, and P¯Umin is the minimum social product of the individual water resources of the investment project.


### 4.3. Solving the MOP of the Reservoir Group Optimization Dispatching Using SFLA

For evaluation of each frog (the solution), the following fitness function can be used:(12)f(x)=1J+α,
where J is the maximum value of the solutions, which is maximum satisfaction; α(0<α<0.1) is a constant. Obviously, the smaller the value of the objective function, the greater the fitness of frog.We generateL frogs (the solution of the problem) as the initial population; the update strategies of each population local depth search are as follows:
(13)U(q)=PW+S,S={min{int[rand(PB-PW)],Smax}Negativedition,Positivedirectionmax{int[rand(PB-PW)],-Smax}Negativection,Negativedirection,
where U(q) is the present moment solution, PW is the last moment solution, S is the moving distance on the components, rand is a random number between 0 and 1, and PB and PW are the optimal and worst solution of the fitness. The steps of the SFLA are as follows.Step 1.
Set parameters and initial frog population.Step 2.
Calculate the fitness of each frog in the population.Step 3.
Sort and divide the population according to fitness in descending.Step 4.
Update each population in part.Step 5.
Mix the updated population and replace the original group with new.Step 6.
If it meets the termination condition, output the optimal solution and finish operation; if not, go Step2.The program flow chart is shown in Figure1.Figure 1
The flow chart of SFLA.


### 4.4. The Improvement of the SFLA: The Weighted Variable Step SFLA

Through studying the basic principle of the SFLA and in order to solve the engineering problems using the SFLA, it is necessary to improve algorithm efficiency and search breadth. The method using the weighted variable step-size search is adopted. We call this method the weighted variable step SFLA.In the process of the basic SFLA, for each decision variable is the same value in rand. So it will limit the random of the decision variables to the optimal direction. This problem will be avoided, if each random decision variable generates different rand values. The step formula can be improved as follows:(14)S(i)={min{int[rand(i)(PB(i)-PW(i))],Smax}PositivedirectionionPositivedirectionmax{int[rand(i)(PB(i)-PW(i))],-Smax}PositivedirectionnNegativedirection.At the same time, in order to keep solution with better fitness and eliminate the solution with bad fitness, to each decision variable, rand is multiplied by a variable weight coefficient which relates to the fitness. In order to reflect the global search, steps do not only relate to present value, but also relate to before value. The step formula can be improved as follows:(15)S(i)={min{int[αf(x)rand(i)(PB(i)-PW(i))+βf(x)S(i-1)],Smax}PositivedirectionoPositivedirectionPositivedirectionmax{int[αf(x)rand(i)(PB(i)-PW(i))+βf(x)S(i-1)],-Smax}NegativedirectiNegativedirection.Negativedirection.αf(x)+βf(x)=1(0<αf(x),βf(x)<1).Figure2 is a decision effect using the improved SFLA, compared with the basic SFLA [20].Figure 2
The effect of the improved SFLA.In the process of local search, we can improve the search method; make improvements in the posterior half of the frog in the population, which is called 1/2 search method. Using this kind of method, the update is more efficient; the transformation can be ensured towards the high fitness direction; the search range can be widened.The optimization process flow of improved SFLA is shown in Figure3.Figure 3
The flow chart of improved SFLA.


### 4.5. The Greatest Satisfaction Transformation of the Reservoir Group MOP

There are five competing goals of the function in the reservoir group optimization dispatching MOP. We can use the fuzzy set theory to translate them into single-objective problem, using the descending half-line as their membership function [21, 22], shown in Figure 4.Figure 4
The objective membership function.The objective membership function can be described as follows:(16)μ(f(x))={1,f(x)<a,b-f(x)b-a,a≤f(x)≤b,0,f(x)>a,
where “a” is the optimization target of the single-objective function, “a+b” is the maximum acceptable values, and λ is the minimum membership variables in all the membership function, called the degree of satisfaction;
(17)λ=min{μ(f1(x),f2(x),f3(x))}.According to the principle of maximum degree of membership, we can translate the MOP into single-objective problem with satisfaction maximization [23].


## 4.1. The Objective Functions of Multireservoir System Optimization Operation

With the above, we take the following as the goal function: people living water consumption is the largest, the industrial and the agricultural production water consumption is the largest, ecological environment water consumption is the largest, the gross output value of industry and agriculture is the highest, and the investment of water resources development is minimum [17].According to the MOP, the objective function, the decision variables, and the constraints of the reservoir group optimization dispatching can be described as follows:(4)min/maxfm(x1,x2,…,xn),m=1,2,…,M;s.t.(1)gk(x1,x2,…,xn)≥0,k=1,2,…,M;(2)hl(x1,x2,…,xn)=0,l=1,2,…,L;(3)xi(L)≤xi≤xi(U),i=1,2,…,n,
where X⊂Rn is the decision space, x=(x1,x2,…,xn)∈X, and xi(L) and xi(H) are upper and lower bounds of variables, respectively. Further, we can define collection Ω as the feasible domain of (4):
(5)Ω={xi(L)x∣gk(x1,…,xn)≥0,hk(x1,…,xn)=0,Ωxi(L)≤xi≤xi(U)},Ω⊂X.Among the objective functions of multireservoir system optimization operation, the investment of water resources development is the minimizing objective function, and the others are the maximum objective functions. In order to maintain unity, by taking the negative of the water resources investment, we convert all the problems into the maximum objective functions MOP. The objective functions can be described as follows:(6)maxf1=∑t=1T∑iIa1iPGi(t)2+β1iPGi(t)+γ1i,maxf2=∑t=1T∑jJa2jPSj(t)2+β2jPSj(t)+γ2j,maxf3=∑t=1T∑kKa3kPWk(t)2+β3kPWk(t)+γ3k,maxf4=∑t=1T∑lLa4lPHl(t)2+β4lPHl(t)+γ4l,maxf5=∑t=1T∑mMa5mPUm(t)2+β5mPUm(t)+γ5m.For (6), I, J, K, L, and M denote the effects of five optimal goal influence factors, respectively. PGi(t) is the water consumption of the ith living water user in the tth time interval; PSi(t) is the water consumption of the jth industry and agriculture water user in the tth time interval; PWk(t) is the water consumption of the kth ecological water user in the tth time interval; PHl(t) is the output value of the lth production unit in the tth time interval; PUm(t) is the investment amount of the mth water resources development unit in the tth time interval. an*,βn*,γn*(n=1,2,3,4,5;*=i,j,k,l,m) is the influence coefficient of each factor; T is the scheduling cycle.


## 4.2. The Constraint Condition of Multireservoir System Optimization Operation

In dealing with the constraint condition, according to the objective function, we make multireservoir system operation satisfaction as constraint condition. In view of reservoir group optimal model, the constraint conditions are as follows.The water consumption demand of people living is mainly affected by the climate change and daily change. So, the living water consumption change and its rate of change are always in a certain range. As one of the important indexes which influence the reservoir group operation, under the premise of the whole objective optimization and reducing the consumption of water, as far as possible, the constraint conditions of the living water consumption can be described as follows:(7)PGmin<∑PG<PGmax,Δ∑PG<ΔPGmax,
where ∑PG is the living water consumption of the reservoir group system in unit interval (usually a day). Δ∑PG  is the variation of the living water consumption of the reservoir group system in unit interval. PGmin and PGmax are the minimum and the maximum living water consumption, respectively. ΔPGmax  is the largest living water quantity variation in unit interval [18].The industrial and agricultural water consumption is directly related to the production situation. In addition to this, the agricultural water consumption is also affected by the weather, the rainfall, the plant growth period, and other factors. The industrial water consumption is also affected by the environment temperature. As another important indicator of the reservoir group operation satisfaction, the constraint conditions of the industrial and agricultural water consumption can be described as follows:(8)PSmin<∑PS<PSmax,Δ∑PS<ΔPSmax,
where ∑PS is the industrial and the agricultural water consumption of the reservoir group system in unit interval (usually a day). Δ∑PS  is the variation of the industrial and the agricultural water consumption of the reservoir group system in unit interval. PSmin and PSmax are the minimum and the maximum industrial and agricultural water consumption, respectively. ΔPSmax is the largest industrial and agricultural water quantity variation in unit interval [19].The ecological environmental water consumption is the other important satisfaction index of the reservoir group operation. The constraint conditions of the ecological environmental water consumption can be described as follows:(9)PWmin<∑PW<PWmax,Δ∑PW<ΔPWmax,
where ∑PW is the ecological environmental water consumption of the reservoir group system within the unit interval (usually a day). Δ∑PW is the variation of the ecological environmental water consumption of the reservoir group system within the unit interval. PWmin and PWmax are the minimum and the maximum industrial and agricultural water consumption, respectively. ΔPWmax is the largest ecological environmental water quantity variation within the unit interval. The ecological environmental water consumption is influenced by the climate, environment, and external factors, and it relates to the living water consumption, the industrial, and agricultural water consumption.The gross output value of industry and agriculture is one of the most important satisfaction indexes of the reservoir group operation. We can describe it with the social gross output. So its constraint conditions can be described as follows:(10)∑PH≥PHmin,
where ∑PH is the gross output value of industry and agriculture in a unit of time (usually a month), and ΔPHmin is the minimum gross output value in unit interval.It is the necessary conditions for the investment of the water resources development to ensure the normal operation of the reservoir group. Under the general condition, the minimum investment maximizes the benefits. So, the constraint conditions of the investment of the water resources development can be described as follows:(11)∑PU≤PUmax,∑PUCount(PU)≥P¯Umin,
where ∑PU is the investment of the water resources development, PUmax is the maximum investment of the water resources development, Count(PU) is the water resources investment project number, and P¯Umin is the minimum social product of the individual water resources of the investment project.


## 4.3. Solving the MOP of the Reservoir Group Optimization Dispatching Using SFLA

For evaluation of each frog (the solution), the following fitness function can be used:(12)f(x)=1J+α,
where J is the maximum value of the solutions, which is maximum satisfaction; α(0<α<0.1) is a constant. Obviously, the smaller the value of the objective function, the greater the fitness of frog.We generateL frogs (the solution of the problem) as the initial population; the update strategies of each population local depth search are as follows:
(13)U(q)=PW+S,S={min{int[rand(PB-PW)],Smax}Negativedition,Positivedirectionmax{int[rand(PB-PW)],-Smax}Negativection,Negativedirection,
where U(q) is the present moment solution, PW is the last moment solution, S is the moving distance on the components, rand is a random number between 0 and 1, and PB and PW are the optimal and worst solution of the fitness. The steps of the SFLA are as follows.Step 1.
Set parameters and initial frog population.Step 2.
Calculate the fitness of each frog in the population.Step 3.
Sort and divide the population according to fitness in descending.Step 4.
Update each population in part.Step 5.
Mix the updated population and replace the original group with new.Step 6.
If it meets the termination condition, output the optimal solution and finish operation; if not, go Step2.The program flow chart is shown in Figure1.Figure 1
The flow chart of SFLA.


## 4.4. The Improvement of the SFLA: The Weighted Variable Step SFLA

Through studying the basic principle of the SFLA and in order to solve the engineering problems using the SFLA, it is necessary to improve algorithm efficiency and search breadth. The method using the weighted variable step-size search is adopted. We call this method the weighted variable step SFLA.In the process of the basic SFLA, for each decision variable is the same value in rand. So it will limit the random of the decision variables to the optimal direction. This problem will be avoided, if each random decision variable generates different rand values. The step formula can be improved as follows:(14)S(i)={min{int[rand(i)(PB(i)-PW(i))],Smax}PositivedirectionionPositivedirectionmax{int[rand(i)(PB(i)-PW(i))],-Smax}PositivedirectionnNegativedirection.At the same time, in order to keep solution with better fitness and eliminate the solution with bad fitness, to each decision variable, rand is multiplied by a variable weight coefficient which relates to the fitness. In order to reflect the global search, steps do not only relate to present value, but also relate to before value. The step formula can be improved as follows:(15)S(i)={min{int[αf(x)rand(i)(PB(i)-PW(i))+βf(x)S(i-1)],Smax}PositivedirectionoPositivedirectionPositivedirectionmax{int[αf(x)rand(i)(PB(i)-PW(i))+βf(x)S(i-1)],-Smax}NegativedirectiNegativedirection.Negativedirection.αf(x)+βf(x)=1(0<αf(x),βf(x)<1).Figure2 is a decision effect using the improved SFLA, compared with the basic SFLA [20].Figure 2
The effect of the improved SFLA.In the process of local search, we can improve the search method; make improvements in the posterior half of the frog in the population, which is called 1/2 search method. Using this kind of method, the update is more efficient; the transformation can be ensured towards the high fitness direction; the search range can be widened.The optimization process flow of improved SFLA is shown in Figure3.Figure 3
The flow chart of improved SFLA.


## 4.5. The Greatest Satisfaction Transformation of the Reservoir Group MOP

There are five competing goals of the function in the reservoir group optimization dispatching MOP. We can use the fuzzy set theory to translate them into single-objective problem, using the descending half-line as their membership function [21, 22], shown in Figure 4.Figure 4
The objective membership function.The objective membership function can be described as follows:(16)μ(f(x))={1,f(x)<a,b-f(x)b-a,a≤f(x)≤b,0,f(x)>a,
where “a” is the optimization target of the single-objective function, “a+b” is the maximum acceptable values, and λ is the minimum membership variables in all the membership function, called the degree of satisfaction;
(17)λ=min{μ(f1(x),f2(x),f3(x))}.According to the principle of maximum degree of membership, we can translate the MOP into single-objective problem with satisfaction maximization [23].


## 5. The Conclusion

For competitive multiobjective function of multireservoir system optimization operation in the inland river basin of semiarid areas, using fuzzy set theory, we convert the MOP of the multireservoir system into single-objective problem of maximum satisfaction. The only solution of fuzzy reasoning makes each objective function achieve optimization [24]. Using this method, we simplify the water resources optimization scheduling problem. It is advantageous for the MOP of the inland river basin of arid areas to accelerate the engineering process [25].For the requirements of the hydrology resource and the reservoir scheduling optimization analysis, there are more and more new theories and new technologies applied to the field of water resources. In the example of Urumqi River basin, we resolved the MOP of multireservoir system optimization operation in the semiarid areas by the SFLA. The research can promote the optimal scheduling and scientific management of water resources in the semiarid inland areas [26]. With the development of the city, the problem of rational management is increasingly important. It is of great significance for the development of the cities to resolve the MOP of the water resource optimal scheduling.In the process of research on the optimization operation of the multireservoir system, besides the appropriate optimization algorithm, the processing of the input information is essential. The input information of the system is random in essence. So the system can be considered as a nonlinear stochastic system. At present, the latest research of these issues focuses on theH∞ filtering [27, 28] and the distributed filtering [29]. The issue will be discussed in the late papers.



---

*Source: 354206-2013-05-07.xml*