CN104047257A - Environmental-protection-oriented ecological reservoir capacity determination method for single water supply reservoir - Google Patents

Abstract

The invention relates to an environmental-protection-oriented ecological reservoir capacity determination method for a single water supply reservoir, and belongs to the technical field of environmental protection and comprehensive utilization of resources. Firstly, according to river ecological flow requirements regulated by a river management department, optimization objectives of a reservoir ecological water utilization schedule are determined, optimization variables and constraint conditions are combined, a reservoir ecological water utilization schedule optimization model is established, a weight set with multiple objectives is constructed, and the multiple objectives are simplified to a single objective for solving; secondly, an optimal water supply scheme is determined according to analysis of the economic society and ecology water lacking rate; at last, the ecological reservoir capacity is calculated according to optimized economic society and ecology water supply output. The method overcomes the defect that ecological water supply is not taken into consideration in traditional water supply reservoir capacity setting and meets the demands of reservoirs for ecological water utilization scheduling in the future; the result is figured out according to the optimization model, more reasonable reservoir capacity can be obtained, and the method has the advantages of being clear in theoretical significance, high in maneuverability and easy to apply practically.

Description

A kind of ecological characteristic storages of reservoir of single service reservoir of Environment Oriented protection is determined method

Technical field

The ecological characteristic storages of reservoir that the present invention relates to a kind of single service reservoir of Environment Oriented protection is determined method, belongs to environmental protection and comprehensive utilization of resources technical field.

Background technology

Operation and the scheduling of tradition reservoir, be all control flood, supply water, the social function such as generating and Economic Function greatly turn to target.The scheduling of these reservoirs is used, and has changed natural process and the rule of River Hydrology, has influence on that fish etc. are hydrobiological perches existence, and river ecosystem is threatened.Along with becoming increasingly conspicuous of ecological environment problem, water demand of river ecosystem also just becomes task more and more important in reservoir operation, some scholars have proposed the concept (Dong Zheren of ecological dispatching one after another, 2007), and ecological flow is combined with reservoir operation, Reservoir Operation Scheme and rule (Yang Zhifeng, 2010 of considering Water Requirement have been developed; Dai Huichao, 2012; Liu Pan, 2013).

These methods are all from traditional storage capacity of reservoir and dispatching method angle, by letting out the ecological water yield under water demand for natural service requirement constraint in downstream river course, in concrete enforcement, easily produce the contradiction of water for economy and ecological water, are difficult to coordinate.In recent years there is scholar to propose reservoir ecological characteristic storages of reservoir concept (Liao Sihui, 2011), and it has been discussed to ensureing necessity (Yong Ting, 2013 of ecological flow; Lu Sun Yun, 2013), it determines that method is that the ecological limit supply line of setting is optimized, be still the constrained dispatch on traditional storage capacity and dispatching method, easily make the dry season only have ecological water supply, and human lives's production needs water breach larger, contradiction is still outstanding.In addition, the storage capacity setting in service of traditional Design of Reservoirs only has minimum capacity of a reservoir, utilizable capacity and storage capacity to determine method (Gu Shengping, 2009), according to real needs, be necessary ad hoc ecological characteristic storages of reservoir on the basis of traditional storage capacity, developing ecology storage capacity is determined method.At present, in service reservoir, also do not have definite computational methods of ecological characteristic storages of reservoir, reservoir is in service to be also difficult to distinct and to coordinate the water regulation of human social economy, flood control and Ecological Target.

Summary of the invention

The ecological characteristic storages of reservoir that the object of the invention is to propose a kind of single service reservoir of Environment Oriented protection is determined method; in conjunction with current ecological flow demand and Optimized model instrument; set up the ecological characteristic storages of reservoir of existing single service reservoir and determine method; make service reservoir performance river flow regulate and store and economic society water supply effect in; employ special storage capacity and carry out ecology water supply; ensure environment flow in river course, implementation environment protection and ecological friendly ecological dispatching.

The ecological characteristic storages of reservoir of the single service reservoir of the Environment Oriented protection that the present invention proposes is determined method, and its FB(flow block) as shown in Figure 1, comprises the following steps:

(1), according to the requirement of river management department, determine respectively the Water Requirement in each month of Ecological Control section of service reservoir downstream river wherein k is month, k=1, and 2 ..., 12, by the Water Requirement in each month as ecological flow target;

(2) set up a reservoir ecological water Scheduling Optimization Model, the object function of Optimized model is:

$\mathrm{MinZ}=\underset{i=1}{\overset{n}{\mathrm{\Σ}}}(w_1{W}_{\mathrm{Dlack}}^i+w_2{W}_{\mathrm{Ilack}}^i+w_3{W}_{\mathrm{Alack}}^i)+\underset{j=1}{\overset{m}{\mathrm{\Σ}}}{w}_4{W}_{\mathrm{elack}}^j$

In above formula: for the life water deficit of i contributing region of service reservoir, be industry and the urban water shortage amount of i contributing region, be the agriculture water deficit of i contributing region, for the ecological water deficit of j Ecological Control section of service reservoir downstream river, w ₁, w ₂, w ₃, w ₄for weight coefficient, by generating at random or manually, meet w ₁+ w ₂+ w ₃+ w ₄the contributing region sum of=1, n for being supplied water by service reservoir, m is the Ecological Control section sum of reservoir downstream river;

(3) set about the Ecological Control section of service reservoir, river, service reservoir place, contributing region, reservoir downstream river, the constraints of the Optimized model of conflux node and channel is as follows:

Service reservoir:

V(t+1)＝V(t)+W _{res_in}(t)-W _{res_out}(t)-W _{res_sup}(t)-L _res(t)

V _ds≤V(t)≤V _max

W _{res_out}(t)≤Q _{res_max}

W _{res_sup}(t)≤Q _dmax

River, service reservoir place:

W _{riv_in}(t)＝W _{riv_out}(t)+L _riv(t)

W _{riv_in}(t)≤Q _{riv_max}

W _{riv_out}(t)≤Q _{riv_max}

Contributing region:

W _sup(t)＝W _con(t)+W _ret(t)

W _sup(t)＝W _dem(t)-W _lack(t)

Ecological Control section:

W _eco(t)＝W _edem(t)-W _elack(t)+W _einc(t)

Node confluxes:

$W_{\mathrm{con}\_\mathrm{out}}\left(t\right)=\underset{l=1}{\overset{N}{\mathrm{\Σ}}}{W}_{\mathrm{con}\_\mathrm{in}}^l\left(t\right)$

Channel:

W _sup(t)＝(1-α)×W _{res_sup}(t)

Wherein, the period index that t is Optimized model, t=1,2 ..., T, T is the total time that model calculates, V (t), V (t+1) are respectively the water retention capacity of the service reservoir of t period and t+1 period, W _{res_in}(t), W _{res_out}(t), W _{res_sup}(t), L _res(t) be respectively the inbound traffics, letdown flow of the service reservoir of t period, to output and the loss amount of contributing region, V _dsfor the lower limit storage capacity of service reservoir, V _maxfor the upper limit storage capacity of service reservoir, Q _{res_max}for letting out ability, Q under service reservoir _dmaxfor the diversion ability of service reservoir, W _{riv_in}(t), W _{riv_out}(t), L _riv(t) be respectively the loss amount in river, upstream inbound traffics, downstream river outflow and river, place of the service reservoir of t period, Q _{riv_max}for the conveyance capacity of downstream river, W _sup(t), W _con(t), W _ret(t), W _dem(t), W _lack(t) be respectively the water yield, water requirement and the water deficit of output that t period contributing region obtains, water consumption, recurrence downstream river, W _eco(t), W _edem(t), W _elack(t), W _einc(t) be respectively t period service reservoir downstream river Ecological Control section ecological output, water requirement, water deficit and add a large amount of, W _{con_out}(t) be the outflow of node of confluxing the t period, represent conflux l branch's inbound traffics of node of t period, N is for flowing into the same node branch sum that confluxes, α is the loss coefficient of channel, span is the real number of 0-1, in above-mentioned constraints, the loss coefficient effluent flow management department of letting out the diversion ability of ability, service reservoir, the conveyance capacity of downstream river and channel under the lower limit storage capacity of service reservoir, the upper limit storage capacity of service reservoir, service reservoir provides;

(4) adopt linear programming method to solve, take 10 days as time step, according to river, the upstream inbound traffics W of the service reservoir in 30 years hydrology historical summaries _{riv_in}(t), solve the Optimized model that above-mentioned steps (2) and step (3) form, the service reservoir being optimized is to the output W of contributing region _{res_sup}and the ecological output W of service reservoir (t) _eco(t);

(5) random or artificial generation organized weight coefficient w more ₁, w ₂, w ₃, w ₄, meet w ₁+ w ₂+ w ₃+ w ₄=1, obtain weight set w ^Ω, repeating step (4), calculates the service reservoir of many group optimizations to the output W of contributing region _{res_sup}and the ecological output W of service reservoir (t) _eco(t), form one by the many groups of service reservoirs of optimizing the output W to contributing region _{res_sup}and the ecological output W of service reservoir (t) _eco(t) the Noninferior Solution Set W forming ^Ω(W _{res_sup}(t), W _eco(t));

(6) service reservoir of optimizing according to above-mentioned many groups is to the output W of contributing region _{res_sup}and the ecological output W of service reservoir (t) _eco(t), by following formula, calculate the Deficient Ratio R of human social economy _{s_lack}with ecological Deficient Ratio R _{eco_lack}:

$R_{s\_\mathrm{lack}}=(\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{W}_{\mathrm{dem}}\left(t\right)-\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{W}_{\mathrm{res}\_\sup }\left(t\right))/\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{W}_{\mathrm{dem}}\left(t\right)$

$R_{\mathrm{eco}\_\mathrm{lack}}=(\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{W}_{\mathrm{edem}}\left(t\right)-\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{W}_{\mathrm{eco}}\left(t\right))/\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{W}_{\mathrm{edem}}\left(t\right);$

As the Deficient Ratio R of human social economy _{s_lack}with ecological Deficient Ratio R _{eco_lack}value while equating (or some specific coordination values), from Noninferior Solution Set W ^Ω(W _{res_sup}(t), W _eco(t)), retrieve the Deficient Ratio R with human social economy _{s_lack}with ecological Deficient Ratio R _{eco_lack}(or some specific coordination values) corresponding W under the condition of equating _{res_sup}and W (t) _eco(t);

(7) according to above-mentioned W _{res_sup}and W _eco, utilize following formula to calculate ecological characteristic storages of reservoir V _eco:

$V_{\mathrm{eco}}=V_{\mathrm{active}}g(\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{W}_{\mathrm{eco}}\left(t\right)/(\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{W}_{\mathrm{res}\_\sup }\left(t\right)+\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{W}_{\mathrm{eco}}\left(t\right)))$

Wherein, V _activethe utilizable capacity of the current service reservoir that Wei Cong river management department obtains.

The ecological characteristic storages of reservoir computational methods of the single service reservoir based on environmental protection that the present invention proposes; first according to river ecological flow requirement and the social economy of river management department regulation, need water; comprehensive optimization aim and the constraints of determining the scheduling of reservoir ecological water; set up reservoir ecological water Scheduling Optimization Model; construct multiobject weight set; multiple target is changed into single goal to be solved; determine optimization water supply scheme that human social economy is equal with ecological Deficient Ratio or the optimization water supply scheme of a certain certain preference, carry out the calculating of service reservoir ecological characteristic storages of reservoir.The inventive method has made up the setting of traditional water supply reservoir capacity and has failed to consider the ecological drawback supplying water, on the basis of existing utilizable capacity, storage capacity, ecological characteristic storages of reservoir is set, take into full account the important function of ecological characteristic storages of reservoir in actual reservoir ecological dispatching, met the demand that following reservoir carries out ecological water scheduling; Next is according to reservoir ecological water Scheduling Optimization Model, by random weight, re-generate many group sights, carry out the analysis of human social economy and ecological Deficient Ratio, can access more rational ecological characteristic storages of reservoir, also can ensure better that service reservoir considers the enforcement of ecological scheduling rule, have advantages of that theory significance is clear and definite, workable, easily obtain practical application and popularization.

Accompanying drawing explanation

Fig. 1 is the FB(flow block) that the ecological characteristic storages of reservoir of the single service reservoir of the Environment Oriented protection that proposes of the present invention is determined method.

Fig. 2 is the schematic diagram of the service reservoir ecological water Scheduling Optimization Model set up of the present invention.

In Fig. 2, the 1st, river, the 2nd, channel, the 3rd, the branch of confluxing, the 4th, reservoir, the 5th, the node that confluxes, the 6th, Ecological Control section, the 7th, contributing region agricultural water, the 8th, contributing region industry, city and living water use.

The specific embodiment

The ecological characteristic storages of reservoir of the single service reservoir of the Environment Oriented protection that the present invention proposes is determined method; its FB(flow block) as shown in Figure 1, the schematic diagram of service reservoir ecological water Scheduling Optimization Model as shown in Figure 2, in Fig. 2; the 1st, river; the 2nd, channel, the 3rd, the branch of confluxing, the 4th, reservoir; the 5th, node confluxes; the 6th, Ecological Control section, the 7th, contributing region agricultural water, the 8th, contributing region industry, city and living water use.The method comprises following steps:

(1) according to the requirement of river management department, (or adopt Hydrology mechanics method to calculate river minimum ecological flow, general desirable Average Annual Runoff 10%), determine respectively the Water Requirement in each month of Ecological Control section of service reservoir downstream river wherein k is month, k=1, and 2 ..., 12, by the Water Requirement in each month as ecological flow target;

(2) set up a reservoir ecological water Scheduling Optimization Model, the object function of Optimized model is:

$\mathrm{MinZ}=\underset{i=1}{\overset{n}{\mathrm{\Σ}}}(w_1{W}_{\mathrm{Dlack}}^i+w_2{W}_{\mathrm{Ilack}}^i+w_3{W}_{\mathrm{Alack}}^i)+\underset{j=1}{\overset{m}{\mathrm{\Σ}}}{w}_4{W}_{\mathrm{elack}}^j$

In above formula: for the life water deficit of i contributing region of service reservoir, be industry and the urban water shortage amount of i contributing region, be the agriculture water deficit of i contributing region, for the ecological water deficit of j Ecological Control section of service reservoir downstream river, w ₁, w ₂, w ₃, w ₄for weight coefficient, by generating at random or manually, meet w ₁+ w ₂+ w ₃+ w ₄the contributing region sum of=1, n for being supplied water by service reservoir, m is the Ecological Control section sum of reservoir downstream river;

(3) set about the Ecological Control section of service reservoir, river, service reservoir place, contributing region, reservoir downstream river, the constraints of the Optimized model of conflux node and channel is as follows:

Service reservoir:

V(t+1)＝V(t)+W _{res_in}(t)-W _{res_out}(t)-W _{res_sup}(t)-L _res(t)

V _ds≤V(t)≤V _max

W _{res_out}(t)≤Q _{res_max}

W _{res_sup}(t)≤Q _dmax

River, service reservoir place:

W _{riv_in}(t)＝W _{riv_out}(t)+L _riv(t)

W _{riv_in}(t)≤Q _{riv_max}

W _{riv_out}(t)≤Q _{riv_max}

Contributing region:

W _sup(t)＝W _con(t)+W _ret(t)

W _sup(t)＝W _dem(t)-W _lack(t)

Ecological Control section:

W _eco(t)＝W _edem(t)-W _elack(t)+W _einc(t)

Node confluxes:

$W_{\mathrm{con}\_\mathrm{out}}\left(t\right)=\underset{l=1}{\overset{N}{\mathrm{\Σ}}}{W}_{\mathrm{con}\_\mathrm{in}}^l\left(t\right)$

Channel:

W _sup(t)＝(1-α)×W _{res_sup}(t)

Wherein, the period index that t is Optimized model, t=1,2 ..., T, T is the total time that model calculates, V (t), V (t+1) are respectively the water retention capacity of the service reservoir of t period and t+1 period, W _{res_in}(t), W _{res_out}(t), W _{res_sup}(t), L _res(t) be respectively the inbound traffics, letdown flow of the service reservoir of t period, to output and the loss amount of contributing region, V _dsfor the lower limit storage capacity of service reservoir, V _maxfor the upper limit storage capacity of service reservoir, Q _{res_max}for letting out ability, Q under service reservoir _dmaxfor the diversion ability of service reservoir, W _{riv_in}(t), W _{riv_out}(t), L _riv(t) be respectively the loss amount in river, upstream inbound traffics, downstream river outflow and river, place of the service reservoir of t period, Q _{riv_max}for the conveyance capacity of downstream river, W _sup(t), W _con(t), W _ret(t), W _dem(t), W _lack(t) be respectively the water yield, water requirement and the water deficit of output that t period contributing region obtains, water consumption, recurrence downstream river, W _eco(t), W _edem(t), W _elack(t), W _einc(t) be respectively t period service reservoir downstream river Ecological Control section ecological output, water requirement, water deficit and add a large amount of, W _{con_out}(t) be the outflow of node of confluxing the t period, represent conflux l branch's inbound traffics of node of t period, N is for flowing into the same node branch sum that confluxes, α is the loss coefficient of channel, span is the real number of 0-1, in above-mentioned constraints, the loss coefficient effluent flow management department of letting out the diversion ability of ability, service reservoir, the conveyance capacity of downstream river and channel under the lower limit storage capacity of service reservoir, the upper limit storage capacity of service reservoir, service reservoir provides;

(4) adopt linear programming method to solve, take 10 days as time step, according to river, the upstream inbound traffics W of the service reservoir in 30 years hydrology historical summaries _{riv_in}(t), solve the Optimized model that above-mentioned steps (2) and step (3) form, the service reservoir being optimized is to the output W of contributing region _{res_sup}and the ecological output W of service reservoir (t) _eco(t);

(5) random or artificial generation organized weight coefficient w more ₁, w ₂, w ₃, w ₄, meet w ₁+ w ₂+ w ₃+ w ₄=1, obtain weight set w ^Ω, repeating step (4), calculates the service reservoir of many group optimizations to the output W of contributing region _{res_sup}and the ecological output W of service reservoir (t) _eco(t), form one by the many groups of service reservoirs of optimizing the output W to contributing region _{res_sup}and the ecological output W of service reservoir (t) _eco(t) the Noninferior Solution Set W forming ^Ω(W _{res_sup}(t), W _eco(t));

(6) service reservoir of optimizing according to above-mentioned many groups is to the output W of contributing region _{res_sup}and the ecological output W of service reservoir (t) _eco(t), by following formula, calculate the Deficient Ratio R of human social economy _{s_lack}with ecological Deficient Ratio R _{eco_lack}:

$R_{s\_\mathrm{lack}}=(\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{W}_{\mathrm{dem}}\left(t\right)-\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{W}_{\mathrm{res}\_\sup }\left(t\right))/\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{W}_{\mathrm{dem}}\left(t\right)$

$R_{\mathrm{eco}\_\mathrm{lack}}=(\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{W}_{\mathrm{edem}}\left(t\right)-\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{W}_{\mathrm{eco}}\left(t\right))/\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{W}_{\mathrm{edem}}\left(t\right);$

As the Deficient Ratio R of human social economy _{s_lack}with ecological Deficient Ratio R _{eco_lack}value while equating (or some specific coordination values), from Noninferior Solution Set W ^Ω(W _{res_sup}(t), W _eco(t)), retrieve the Deficient Ratio R with human social economy _{s_lack}with ecological Deficient Ratio R _{eco_lack}(or some specific coordination values) corresponding W under the condition of equating _{res_sup}and W (t) _eco(t);

(7) according to above-mentioned W _{res_sup}and W _eco, utilize following formula to calculate ecological characteristic storages of reservoir V _eco:

$V_{\mathrm{eco}}=V_{\mathrm{active}}g(\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{W}_{\mathrm{eco}}\left(t\right)/(\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{W}_{\mathrm{res}\_\sup }\left(t\right)+\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{W}_{\mathrm{eco}}\left(t\right)))$

Wherein, V _activethe utilizable capacity of the current service reservoir that Wei Cong river management department obtains.

Below introduce examples of implementation of the inventive method:

Certain service reservoir is positioned at upper river, and level of dead water is 2560m, and minimum capacity of a reservoir is 0.145 hundred million m ³, emerging sharp water level 2624m, corresponding adjusting storage capacity is 3.52 hundred million m ³.The definite water yield of Water Requirement Shi You watershed management department of this reservoir.For coordinating contributing region, middle reaches, river and downstream ecology water requirement, improve the ecological environment, adopt the inventive method to calculate ecological characteristic storages of reservoir, ecology supplies water downstream specially.

Determine object function and constraints, set up the ecological water Scheduling Optimization Model that comprises this service reservoir, construct multiobject weight set, multiple target is changed into single goal and solve.According to the analysis of economic society and ecological Deficient Ratio, determine and optimize water supply scheme.Economic society and ecological output according to optimizing, carry out the calculating of ecological characteristic storages of reservoir in the ratio that ecological output is shared, as follows:

V _eco＝V _activeg(W _eco/(W _{res_sup}+W _eco))

V wherein _ecoecological characteristic storages of reservoir, V _activethe utilizable capacity of the current service reservoir that Shi Cong river management department obtains, W _ecothe ecological output of service reservoir, W _{res_sup}it is the economic society output of service reservoir.

The ecological characteristic storages of reservoir that calculates this service reservoir is 0.64 hundred million m ³, account for 18.3% of utilizable capacity.This ecological characteristic storages of reservoir is relatively set and without the ecological lack of water situation under ecological characteristic storages of reservoir, as shown in the table:

Ecological lack of water situation	Without ecological characteristic storages of reservoir	There is ecological characteristic storages of reservoir	Changes in contrast
				Ecological water deficit (hundred million m 3)	1.96	1.55	0.41
Ecological Deficient Ratio	15.8％	11.6％	4.2％

As seen from table, and compare without ecological characteristic storages of reservoir result, the ecological characteristic storages of reservoir that utilizes the inventive method to calculate, the water deficit of ecological environment has reduced by 0.41 hundred million m ³, Deficient Ratio has reduced 4.2%, to improving downstream ecology water requirement, has positive effect.

Claims (1)

1. the ecological characteristic storages of reservoir of the single service reservoir of Environment Oriented protection is determined a method, it is characterized in that these computational methods comprise following steps:

(1), according to the requirement of river management department, determine respectively the Water Requirement in each month of Ecological Control section of service reservoir downstream river wherein k is month, k=1, and 2 ..., 12, by the Water Requirement in each month as ecological flow target;

(2) set up a reservoir ecological water Scheduling Optimization Model, the object function of Optimized model is:

$\mathrm{MinZ}=\underset{i=1}{\overset{n}{\mathrm{\Σ}}}(w_1{W}_{\mathrm{Dlack}}^i+w_2{W}_{\mathrm{Ilack}}^i+w_3{W}_{\mathrm{Alack}}^i)+\underset{j=1}{\overset{m}{\mathrm{\Σ}}}{w}_4{W}_{\mathrm{elack}}^j$

In above formula: for the life water deficit of i contributing region of service reservoir, be industry and the urban water shortage amount of i contributing region, be the agriculture water deficit of i contributing region, for the ecological water deficit of j Ecological Control section of service reservoir downstream river, w ₁, w ₂, w ₃, w ₄for weight coefficient, by generating at random or manually, meet w ₁+ w ₂+ w ₃+ w ₄the contributing region sum of=1, n for being supplied water by service reservoir, m is the Ecological Control section sum of reservoir downstream river;

(3) set about the Ecological Control section of service reservoir, river, service reservoir place, contributing region, reservoir downstream river, the constraints of the Optimized model of conflux node and channel is as follows:

Service reservoir:

V(t+1)＝V(t)+W _{res_in}(t)-W _{res_out}(t)-W _{res_sup}(t)-L _res(t)

V _ds≤V(t)≤V _max

W _{res_out}(t)≤Q _{res_max}

W _{res_sup}(t)≤Q _dmax

River, service reservoir place:

W _{riv_in}(t)＝W _{riv_out}(t)+L _riv(t)

W _{riv_in}(t)≤Q _{riv_max}

W _{riv_out}(t)≤Q _{riv_max}

Contributing region:

W _sup(t)＝W _con(t)+W _ret(t)

W _sup(t)＝W _dem(t)-W _lack(t)

Ecological Control section:

W _eco(t)＝W _edem(t)-W _elack(t)+W _einc(t)

Node confluxes:

$W_{\mathrm{con}\_\mathrm{out}}\left(t\right)=\underset{l=1}{\overset{N}{\mathrm{\Σ}}}{W}_{\mathrm{con}\_\mathrm{in}}^l\left(t\right)$

Channel:

W _sup(t)＝(1-α)×W _{res_sup}(t)

Wherein, the period index that t is Optimized model, t=1,2 ..., T, T is the total time that model calculates, V (t), V (t+1) are respectively the water retention capacity of the service reservoir of t period and t+1 period, W _{res_in}(t), W _{res_out}(t), W _{res_sup}(t), L _res(t) be respectively the inbound traffics, letdown flow of the service reservoir of t period, to output and the loss amount of contributing region, V _dsfor the lower limit storage capacity of service reservoir, V _maxfor the upper limit storage capacity of service reservoir, Q _{res_max}for letting out ability, Q under service reservoir _dmaxfor the diversion ability of service reservoir, W _{riv_in}(t), W _{riv_out}(t), L _riv(t) be respectively the loss amount in river, upstream inbound traffics, downstream river outflow and river, place of the service reservoir of t period, Q _{riv_max}for the conveyance capacity of downstream river, W _sup(t), W _con(t), W _ret(t), W _dem(t), W _lack(t) be respectively the water yield, water requirement and the water deficit of output that t period contributing region obtains, water consumption, recurrence downstream river, W _eco(t), W _edem(t), W _elack(t), W _einc(t) be respectively t period service reservoir downstream river Ecological Control section ecological output, water requirement, water deficit and add a large amount of, W _{con_out}(t) be the outflow of node of confluxing the t period, represent conflux l branch's inbound traffics of node of t period, N is for flowing into the same node branch sum that confluxes, α is the loss coefficient of channel, span is the real number of 0-1, in above-mentioned constraints, the loss coefficient effluent flow management department of letting out the diversion ability of ability, service reservoir, the conveyance capacity of downstream river and channel under the lower limit storage capacity of service reservoir, the upper limit storage capacity of service reservoir, service reservoir provides;

(4) adopt linear programming method to solve, take 10 days as time step, according to river, the upstream inbound traffics W of the service reservoir in 30 years hydrology historical summaries _{riv_in}(t), solve the Optimized model that above-mentioned steps (2) and step (3) form, the service reservoir being optimized is to the output W of contributing region _{res_sup}and the ecological output W of service reservoir (t) _eco(t);

(5) random or artificial generation organized weight coefficient w more ₁, w ₂, w ₃, w ₄, meet w ₁+ w ₂+ w ₃+ w ₄=1, obtain weight set w ^Ω, repeating step (4), calculates the service reservoir of many group optimizations to the output W of contributing region _{res_sup}and the ecological output W of service reservoir (t) _eco(t), form one by the many groups of service reservoirs of optimizing the output W to contributing region _{res_sup}and the ecological output W of service reservoir (t) _eco(t) the Noninferior Solution Set W forming ^Ω(W _{res_sup}(t), W _eco(t));

(6) service reservoir of optimizing according to above-mentioned many groups is to the output W of contributing region _{res_sup}and the ecological output W of service reservoir (t) _eco(t), by following formula, calculate the Deficient Ratio R of human social economy _{s_lack}with ecological Deficient Ratio R _{eco_lack}:

$R_{s\_\mathrm{lack}}=(\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{W}_{\mathrm{dem}}\left(t\right)-\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{W}_{\mathrm{res}\_\sup }\left(t\right))/\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{W}_{\mathrm{dem}}\left(t\right)$

$R_{\mathrm{eco}\_\mathrm{lack}}=(\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{W}_{\mathrm{edem}}\left(t\right)-\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{W}_{\mathrm{eco}}\left(t\right))/\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{W}_{\mathrm{edem}}\left(t\right);$

As the Deficient Ratio R of human social economy _{s_lack}with ecological Deficient Ratio R _{eco_lack}value while equating (or some specific coordination values), from Noninferior Solution Set W ^Ω(W _{res_sup}(t), W _eco(t)), retrieve the Deficient Ratio R with human social economy _{s_lack}with ecological Deficient Ratio R _{eco_lack}(or some specific coordination values) corresponding W under the condition of equating _{res_sup}and W (t) _eco(t);

(7) according to above-mentioned W _{res_sup}and W _eco, utilize following formula to calculate ecological characteristic storages of reservoir V _eco:

$V_{\mathrm{eco}}=V_{\mathrm{active}}g(\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{W}_{\mathrm{eco}}\left(t\right)/(\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{W}_{\mathrm{res}\_\sup }\left(t\right)+\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{W}_{\mathrm{eco}}\left(t\right)))$

Wherein, V _activethe utilizable capacity of the current service reservoir that Wei Cong river management department obtains.