CN109829580A - Abundant irrigation conditions lower storage reservoir-benefit library pumping plant-benefit canal pumping station system water resource optimal allocation method - Google Patents
Abundant irrigation conditions lower storage reservoir-benefit library pumping plant-benefit canal pumping station system water resource optimal allocation method Download PDFInfo
- Publication number
- CN109829580A CN109829580A CN201910061144.0A CN201910061144A CN109829580A CN 109829580 A CN109829580 A CN 109829580A CN 201910061144 A CN201910061144 A CN 201910061144A CN 109829580 A CN109829580 A CN 109829580A
- Authority
- CN
- China
- Prior art keywords
- water
- reservoir
- pumping plant
- library
- benefit
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 319
- 238000005086 pumping Methods 0.000 title claims abstract description 217
- 238000000034 method Methods 0.000 title claims abstract description 107
- 238000003860 storage Methods 0.000 title claims abstract description 39
- 238000003973 irrigation Methods 0.000 title claims abstract description 27
- 230000002262 irrigation Effects 0.000 title claims abstract description 27
- 230000008569 process Effects 0.000 claims abstract description 61
- 230000003020 moisturizing effect Effects 0.000 claims abstract description 26
- 230000006735 deficit Effects 0.000 claims abstract description 8
- 238000013178 mathematical model Methods 0.000 claims abstract description 8
- 241000196324 Embryophyta Species 0.000 claims description 158
- 230000008901 benefit Effects 0.000 claims description 38
- 230000007704 transition Effects 0.000 claims description 18
- 230000008020 evaporation Effects 0.000 claims description 7
- 238000001704 evaporation Methods 0.000 claims description 7
- 238000005457 optimization Methods 0.000 claims description 7
- 238000011156 evaluation Methods 0.000 claims description 6
- 239000000203 mixture Substances 0.000 claims description 6
- 238000011160 research Methods 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 238000002360 preparation method Methods 0.000 claims description 3
- 239000008400 supply water Substances 0.000 claims description 2
- JEGUKCSWCFPDGT-UHFFFAOYSA-N h2o hydrate Chemical compound O.O JEGUKCSWCFPDGT-UHFFFAOYSA-N 0.000 claims 1
- 238000000354 decomposition reaction Methods 0.000 abstract description 4
- 238000013459 approach Methods 0.000 abstract description 2
- 230000008439 repair process Effects 0.000 description 4
- 241000372132 Hydrometridae Species 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000007726 management method Methods 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
Landscapes
- Feedback Control In General (AREA)
Abstract
Abundant irrigation conditions lower storage reservoir-benefit library pumping plant-benefit canal pumping station system water resource optimal allocation method, belong to Water Resources Irrigation and distributes technical field rationally, pass through the intake area water shortage situation under different water frequencies, library pumping plant is mended to reservoir moisturizing from single seat, there is single seat to mend canal pumping plant to feeder channel water delivery again simultaneously, establish one library of abundant irrigation conditions lower storage reservoir-benefit library pumping plant-benefit canal pumping plant, two station system water resource joint Optimized Operation mathematical model, it is solved using Approach by inchmeal, one library, two station system decomposition is obtained into the mono- single pumping plant subsystem Water Resources Allocation Model for mending library station subsystem Water Resources Allocation Model and directly mend canal of abundant irrigation conditions lower storage reservoir-, obtain minimum water deficit in intake area in certain delivery period, day part optimal water supply in corresponding reservoir delivery period, abandon process water, and mend library pumping plant, mend canal pump The optimal rate of water make-up process of day part of standing has most important theories meaning and practical application value.
Description
Technical field
The invention belongs to Water Resources Irrigations to distribute technical field rationally, be related to a kind of water resource optimal allocation method, specifically
Say be related to it is a kind of from 1 mend the moisturizing of library pumping plant single reservoir and 1 benefit canal pumping plant combine to intake area supply water a library
The two system combined traffic control methods in station.
Background technique
Current many areas restrict economic society and other field development since the spatial and temporal distributions of water resource are uneven.For
For the irrigated area of abundant irrigation conditions, under limited water resources total amount and water source project scale, the optimization of reinforced region water resource
Scheduling and management use the hydraulic engineering of irrigation system as unified entirety with turning to target with water comprehensive benefit maximum
And adjusting, so that it is played bigger effect with built engineering (such as reservoir and pumping plant combined dispatching are run), is to solve irrigated area to lack
The main path of water problems.One library, the two station system traffic control of " reservoir-benefit library pumping plant-benefit canal pumping plant " is water resources management one
Item content, it is hydraulic engineering technology that how reasonably to reach some target with system water scheduling of resource, which keeps system benefit best,
Middle urgent problem.
Summary of the invention
The purpose of the present invention is, since the spatial and temporal distributions of water resource are uneven, restricted for current many area economic society and
Other field development, and for the irrigated area of abundant irrigation conditions, under limited water resources total amount and water source project scale, area
It the Optimized Operation of domain water resource and manages not reasonable, proposes " reservoir-benefit library pumping plant-benefit canal pump under a kind of abundant irrigation conditions
Stand " two station system water resource optimization configuration method of a library reservoir can be given full play to by the method for the water resource optimal allocation
Pondage capacity, reduce the abandoning water of reservoir, the effective use of water resource can be better achieved, to reduce the water shortage of intake area
Amount.
Technical solution of the present invention: the two station system of a library of " reservoir-benefit library pumping plant-benefit canal pumping plant " under abundant irrigation conditions
Water resource optimal allocation method, which is characterized in that from 1 mend the moisturizing of library pumping plant single reservoir and 1 benefit canal pumping plant combine to
The co-supplying that intake area is supplied water, comprising the following steps:
(1) model construction
(1-1) with day part intake area water requirement in the 2 station system year an of library of " reservoir-benefit library pumping plant-benefit canal pumping plant " with
The least square of system total supply difference and be target, establish following objective function:
In formula: G be research object year in day part the difference for water requirement least square and；S is in research object year
The difference for water requirement of day part, ten thousand m3；N is the when number of segment divided in year；Segment number when i is (i=1,2 ..., N)；XiFor water
The water supply of i-th period of library, ten thousand m3；YiFor the rate of water make-up for mending the i-th period of canal pumping plant, ten thousand m3；DiFor the need of the i-th period of intake area
Water, ten thousand m3；Objective function is to accelerate to reduce between the system water supply amount and intake area water requirement using quadratic sum expression
Deviation；
Constraint condition is arranged in (1-2), for water inventory constraint condition, mends canal pumping plant in " reservoir-benefit library pumping plant " subsystem year
Year allows water lift total amount constraint condition, reservoir operation criterion constraint condition and reservoir capacity constraint condition；
(2) model solution
(2-1) data preparation was divided into N number of period for 1 year, and determines day part length；It is bent according to water level-storage capacity relationship
Line determines reservoir initial storage V0；Determine reservoir year for water inventory K, minimum capacity of a reservoir Vmin, the corresponding storage capacity V of flood controlP、
And utilizable capacity Vmin+Δ1；Measurement reservoir day part carrys out water LSi, evaporation with leakage EFi；Determine that mending library pumping plant year allows
Water lift total amount BK；Determine that mending canal pumping plant year allows water lift total amount BQ；Determine intake area day part under certain water frequency needs water
Measure Di(i=1,2 ..., N)；
(2-2) solves system model using Dynamic Programming successive approximation method, obtains day part reservoir optimal water supply mistake
Journey Xi, optimal abandoning process water PSi, mend the optimal rate of water make-up process Z of library Group of Pumping Stationi, and mend the optimal rate of water make-up process of canal pumping plant
Yi。
Constraint condition described in step (1-2) are as follows:
(1) " reservoir-Dan Buku pumping plant " system year is for water inventory constraint condition: in different level year, different fractions
In the case of, consider that intake area needs water requirement, the water that reservoir-benefit library pumping plant water supply project may provide；
In formula: K is the year of reservoir for water inventory, ten thousand m3；BK is to mend to allow water lift total amount, ten thousand m in library pumping plant year3；
(2) canal pumping plant year water lift total amount constraint condition is mended:
In formula: BQ is to mend to allow water lift total amount, ten thousand m in canal pumping plant year3；
(3) reservoir operation criterion constraint condition: according to reservoir water equilibrium equation,
Vi=Vi-1+LSi+Zi-PSi-EFi-Xi, (i=1,2 ..., N) (4)
(3-1) is lower than reservoir minimum capacity of a reservoir V when the i-th period end pondageminWhen, then the i-th period should be given by mending library pumping plant
Reservoir moisturizing, moisturizing to utilizable capacity (Δ more than reservoir minimum capacity of a reservoir1), it may be assumed that
Vi<VminWhen: Zi=Vmin-Vi+Δ1； (5)
This period reservoir abandons water PSi=0；
(3-2) when meeting with flood, the i-th period end pondage is greater than pondage corresponding to flood control
VPWhen, then the i-th period should carry out abandoning water to reservoir, abandon water to flood control by reservoir regulation limiting water level, it may be assumed that
Vi>VPWhen: PSi=Vi-VP； (6)
This period mends library pumping plant rate of water make-up Zi=0；
(3-3) is when the i-th period end pondage is between minimum capacity of a reservoir VminWith reservoir filling corresponding to flood control
Measure VPBetween, then the i-th period reservoir does not need to abandon water, mends library pumping plant and does not need moisturizing yet, it may be assumed that
Vmin≤Vi≤VPWhen: Zi=PSi=0； (7)
In formula: Vi、Vi-1Respectively reservoir i-th, the reservoir storage of i-1 period Mo, ten thousand m3；ZiFor the benefit library pumping plant of the i-th period
Rate of water make-up, ten thousand m3；LSi、PSi、EFiRespectively the i-th period of reservoir comes water, ten thousand m3, abandon water, ten thousand m3；Evaporation and leakage,
Ten thousand m3；Vmin、VPRespectively storage capacity corresponding to the minimum capacity of a reservoir of reservoir and flood control, ten thousand m3；
(4) reservoir capacity constraint condition: the pondage of day part should be between reservoir minimum capacity of a reservoir and flood control
Between corresponding storage capacity, it may be assumed that
Vmin≤Vi≤VP, (i=1,2 ..., N) (8)
The method that system model is solved using Dynamic Programming successive approximation method described in step (2-2) are as follows:
(1) it drafts and mends the initial rate of water make-up Y of canal pumping plant1, i(i=1,2 ..., N), it is desirable that meet and mend canal pumping plant year water lift total amount
Former intake area is regarded as one because mending at this point, deleting the benefit canal pumping plant subsystem from original system by constraint condition, i.e. formula (3)
A new intake area after the diminution of canal pumping plant moisturizing occurring area, the other compositions of system remain unchanged, i.e., sufficiently irrigate item for former
The conversion of the two station system Water Resources Allocation Model of a library of " reservoir-benefit library pumping plant-benefit canal pumping plant " under part, i.e. formula (1)~(8)
It is W for each stage water requirementi" reservoir-Dan Buku pumping plant " of the corresponding of no help canal pumping plant in new intake area of (i=1,2 ..., N)
System Water Resources Allocation Model:
Objective function:
Wherein: Wi=Di-Y1,i, due to day part Di、Y1,iIt is known that then WiIt is i.e. known；
Constraint condition:
(1-1) " reservoir-Dan Buku pumping plant " system year is for water inventory constraint condition: i.e. formula (2)；
(1-2) reservoir operation criterion constraint condition: i.e. formula (4)~(7)；
(1-3) reservoir capacity constraint condition: i.e. formula (8)；
(2) " reservoir-Dan Buku pumping plant " system water resource optimal allocation mathematical model, i.e. formula (2), (4)~(9) be with by
The minimum target of quadratic sum of the difference of pool day part water deficit, can divide by the stage of decision variable of day part reservoir yield
One-dimensional dynamic programming model obtain corresponding recurrence equation referring to one-dimensional dynamic programming evaluation principle are as follows:
(2-1) stage i=1:
g1(λ1)=min (X1,1-W1)2 (10)
In formula, state variable λ1Indicate preceding 1 water supply period reservoir yield, it can be discrete in corresponding feasible zone: λ1=
0,E1,E2,...,K+BK；To each discrete λ1, decision variable X1,1(reservoir yield) can be discrete in corresponding feasible zone, such as
00000 m3, 200,000 m3, 400,000 m3, 600,000 m3、…X1,1,maxDeng (X1,1,maxFor the 1st stage reservoir maximum water supply capacity), it should meet:
X1,1≥λ1；By X1,1Formula (10) are substituted into respectively, respectively obtain each discrete λ1When value, optimal X1,1And its corresponding g1(λ1)；
Then, according to formula (4), the 1st stage end reservoir capacity V1=V0+LS1-EF1-X1,1, not yet consider to mend library pump at this time
Stand moisturizing or reservoir abandons water, and Ying Caiyong formula (5)~(7) are tested and corrected:
(2-1-1) works as V1<Vmin, then consider to carry out moisturizing, Z to reservoir by mending library pumping plant1,1=Vmin-V1+Δ1, repair at this time
Positive storage capacity V1 *=Vmin+Δ1；
(2-1-2) works as V1> VP, then need to abandon dispatching requirement of the water to guarantee reservoir capacity, PS1,1=V1-VP, repair at this time
Positive storage capacity V1 *=VP；
(2-1-3) works as Vmin≤V1≤VP, then Z1,1=PS1,1=0, V at this time1 *=V1；
By step (2-1-1)~(2-1-3), corrects and determine the 1st stage end reservoir capacity V1 *, while can get and correspond to
Reservoir abandon water PS1,1Or mend library pumping plant rate of water make-up Z1,1；
(2-2) stage i=2,3 ..., N-1:
gi(λi)=min [(X1,i-Wi)2+gi-1(λi-1)] (11)
In formula, state variable λiFor the preceding i period reservoir for water inventory, equally carry out respectively discrete: λi=0, E1,
E2,...,K+BK；To each discrete λi, decision variable (reservoir yield X1,i) it is discrete ibid, and should meet:
State transition equation: λi-1=λi-X1,i (12)
In formula: i=2,3 ..., N-1
By each discrete X1,iValue substitutes into the (X in formula (11) respectively1,i-Wi)2, by state transition equation formula (12), to each
A discrete X1,i, search i-1 stage mingi-1(λi-1) value, thus to obtain (X1,i-Wi)2+min gi-1(λi-1), complete the above institute
There is discrete X1,iAfter optimizing, final can get meets min [(X1,i-Wi)2+gi-1(λi-1)] require optimal X1,iProcess and its right
The g answeredi(λi)；
Equally, it also needs to determine the i-th stage end storage capacity V using formula (4)i, then tested, corrected using formula (5)~(7)
Determine the i-th period end storage capacity Vi *, while can get corresponding reservoir and abandoning process water PS1,iWith benefit library pumping plant rate of water make-up process Z1,i
(i=2,3 ..., N-1), the same step of process (2-1-1)~(2-1-3)；
(2-3) stage N:
gN(λN)=min [(X1,N-WN)2+gN-1(λN-1)] (13)
In formula, state variable λNFor the top n period reservoir for water inventory: λN=K+BK；Decision variable (reservoir yield
X1,N) equally discrete in corresponding feasible zone；
State transition equation: λN-1=λN-X1,N (14)
It is final to obtain objective function optimal value G using step (2-2) the method1=gN(λN) and corresponding reservoir
Optimal water supply process X1,i, abandon process water PS1,i, mend library pumping plant rate of water make-up process Z1,i(i=1,2 ..., N)；
(3) X for being calculated step (1)1,iProcess substitutes into master mould, i.e. formula (1)~(8) as input condition, at this point, by
In reservoir yield X1,iWith benefit library pumping plant rate of water make-up Z1,iIt is known that can should " reservoir-Dan Buku pumping plant " subsystem from former system
It is deleted in system, former intake area is regarded as one because one after the diminution of " reservoir-Dan Buku pumping plant " co-supplying occurring area is new
Intake area, the other compositions of system remain unchanged；It can be by " reservoir-benefit library pumping plant-benefit canal pumping plant " under former sufficiently irrigation conditions
It is U that one library, two station system Water Resources Allocation Model, i.e. formula (1)~(8), which are converted into each stage water requirement,i(i=1,2 ...,
N the corresponding single benefit canal pumping plant Water Resources Allocation Model in new intake area):
Objective function:
Wherein: Ui=Di-X1,i, due to day part Di、X1,iIt is known that then UiIt is i.e. known；
Constraint condition:
Mend canal pumping plant year water lift total amount constraint condition: i.e. formula (3)；
This singly mends canal pumping plant Water Resources Allocation Model, i.e., formula (3), (15) are with the difference of intake area day part water deficit
The minimum target of quadratic sum, the one-dimensional Dynamic Programming mould that can divide as the stage of decision variable of canal pumping plant rate of water make-up is mended using day part
Type obtains corresponding recurrence equation still referring to one-dimensional dynamic programming evaluation principle are as follows:
(3-1) stage i=1:
g1(α1)=min (Y2,1-U1)2 (16)
In formula, state variable α1Canal pumping plant water lift amount is mended for preceding 1 period, a fixed step size can be pressed in corresponding feasible zone
It is discrete: α1=0, F1,F2,...,BQ；To each discrete α1, decision variable Y2,1(mending canal pumping plant rate of water make-up) can be feasible in correspondence
It is discrete in domain, such as 00,000 m3, 200,000 m3, 400,000 m3, 600,000 m3、…Y2,1,maxDeng (Y2,1,maxCanal pumping plant maximum is mended for the 1st stage to mention
Outlet capacity), it should meet: Y2,1≥α1；By Y2,1Formula (16) are substituted into respectively, respectively obtain each discrete α1When value, optimal Y2,1And its
Corresponding g1(α1)；
(3-2) stage i=2,3 ... N-1:
gi(αi)=min [(Y2,i-Ui)2+gi-1(αi-1)] (17)
In formula: state variable αiCanal pumping plant water lift total amount is mended for the preceding i period, is equally carried out respectively discrete: αi=0, F1,
F2,...,BQ；To each discrete αi, decision variable Y2,i(reservoir yield) is discrete ibid, and should meet:
State transition equation: αi-1=αi-Y2,i (18)
By each discrete Y2,iValue substitutes into the (Y in formula (17) respectively2,i-Ui)2, by state transition equation formula (18), search i-
1 stage min gi-1(αi-1) value, thus to obtain (Y2,i-Ui)2+mingi-1(αi-1), complete all of above discrete Y2,iAfter optimizing,
Final can get meets min [(Y2,i-Ui)2+gi-1(αi-1)] require optimal Y2,iProcess and its corresponding gi(αi)；
(3-3) stage N:
gN(αN)=min [(Y2,N-UN)2+gN-1(αN-1)] (19)
In formula, state variable αNFor the benefit canal pumping plant water lift total amount of top n period, αN=BQ；Decision variable Y2,N(mend canal pump
Stand rate of water make-up) it is equally discrete in corresponding feasible zone；
State transition equation: αN-1=αN-Y2,N (20)
Using step (3-2) the method, final obtain meets the αNIt is required that the optimal rate of water make-up process Y of benefit canal pumping plant2,i
(i=1,2 ..., N), and objective function optimal value G2=gN(αN)；
(4) by the Y of acquisition2,iProcess substitutes into master mould and to be again converted into " reservoir-Dan Buku pumping plant " system water resource excellent
Change configuration mathematical model, i.e. formula (2), (4)~(9), repeats step (1) and solve, obtain new reservoir yield process X3,i, abandon
Process water PS3,i, mend library pumping plant rate of water make-up process Z3,iAnd objective function optimal value G3；
(5) by the X of acquisition3,iProcess substitutes into master mould again and is converted into single benefit canal pumping plant Water Resources Allocation Model,
That is formula (3), (15) repeat step (2) and solve, obtain new benefit canal pumping plant rate of water make-up process Y4,iAnd objective function optimal value
G4；
(6) above step is repeated, the objective function optimal value obtained twice up to date meets | (Gn-Gn-1)/Gn|≤ε(ε
For iteration control precision), then iteration terminates；The G that will be obtained for the last time at this timenIt is corresponding each as objective function optimal value
Period reservoir optimal water supply process Xn-1,i, optimal abandoning process water PSn-1,i, mend the optimal rate of water make-up process Z of library pumping plantn-1,i, with
And mend the optimal rate of water make-up process Y of canal pumping plantn,iStrategy be original system water resource combined dispatching optimal policy.
The invention has the benefit that " reservoir-benefit library pumping plant-benefit canal under a kind of abundant irrigation conditions proposed by the present invention
The two station system water resource optimization configuration method of a library of pumping plant ", the water system are mainly made of 1 reservoir and 2 pumping plants,
Middle reservoir supplies water from 1 pumping plant to its moisturizing (i.e. benefit library pumping plant), then to intake area；Separately there is 1 pumping plant then to pass through feeder channel
It directly supplies water to intake area, thus constitutes the two station system of a library of " reservoir-benefit library pumping plant-benefit canal pumping plant ", combine to intake area
It supplies water.The present invention is directed to the intake area water shortage situation under different water frequencies, mends library pumping plant to reservoir moisturizing from single seat, while again
There is single seat to mend canal pumping plant to feeder channel water delivery, establishes under abundant irrigation conditions a library of " reservoir-benefit library pumping plant-benefit canal pumping plant "
Two station system water resource joint Optimized Operation mathematical models.Scheduling system is adjusted for specific year, supplies water and draws in known reservoir
Point when number of segment, initial storage, minimum capacity of a reservoir, utilizable capacity, the corresponding storage capacity of flood control, year for water inventory, it is each when
Section comes process water, evaporation and leakage process, and mending allows water lift total amount in library pumping plant year, and mending allows water lift total amount in canal pumping plant year,
And under the water requirement process condition of day part intake area, using Approach by inchmeal method for solving, two station system decomposition of a library is obtained
Single pumping plant subsystem of " reservoir-is mono- to mend library station " subsystem Water Resources Allocation Model and directly benefit canal under abundant irrigation conditions
Water Resources Allocation Model is respectively adopted one-dimensional dynamic programming and solves, and passes through gradually forcing for 2 subsystem model target values
Closely, minimum water deficit in intake area in certain delivery period, day part optimal water supply, abandoning water in corresponding reservoir delivery period are obtained
Process, and mend library pumping plant, mend the optimal rate of water make-up process of canal pumping plant day part.By to the system water resource optimal allocation, energy
The pondage capacity for giving full play to reservoir reduces the abandoning water of reservoir, and the effective use of water resource is better achieved, and realizes and sufficiently fills
The water resources optimal operation of the two station system of a library of " reservoir-benefit library pumping plant-benefit canal pumping plant " under the conditions of irrigating, method precision is reliable,
Meanwhile reducing and mending library pumping station operation energy consumption, achieve the purpose that Water Resources Irrigation is distributed rationally, improves irrigated area, society
And ecological benefits.
Detailed description of the invention
Fig. 1 is two station system of the library signal of " reservoir-benefit library pumping plant-benefit canal pumping plant " under sufficiently irrigation conditions of the invention
Figure.
Fig. 2 is " the mono- benefit library station of reservoir-" subsystem under the abundant irrigation conditions that two station system decomposition of a library of the invention obtains
Schematic diagram.
Fig. 3 is single pumping plant subsystem that canal is directly mended under the abundant irrigation conditions that two station system decomposition of a library of the invention obtains
Schematic diagram.
Specific embodiment
The present invention will be further explained below with reference to the attached drawings:
As shown in Figure 1-3, under a kind of abundant irrigation conditions " reservoir-benefit library pumping plant-benefit canal pumping plant " two station system of a library
Water resource optimal allocation method, which is characterized in that from 1 mend the moisturizing of library pumping plant single reservoir and 1 benefit canal pumping plant combine to
The co-supplying that intake area is supplied water, comprising the following steps:
(1) model construction
(1-1) with day part intake area water requirement in the 2 station system year an of library of " reservoir-benefit library pumping plant-benefit canal pumping plant " with
The least square of system total supply difference and be target, establish following objective function:
In formula: G be research object year in day part the difference for water requirement least square and；S is in research object year
The difference for water requirement of day part, ten thousand m3；N is the when number of segment divided in year；Segment number when i is (i=1,2 ..., N)；XiFor water
The water supply of i-th period of library, ten thousand m3；YiFor the rate of water make-up for mending the i-th period of canal pumping plant, ten thousand m3；DiFor the need of the i-th period of intake area
Water, ten thousand m3；Objective function is to accelerate to reduce between the system water supply amount and intake area water requirement using quadratic sum expression
Deviation；
Constraint condition is arranged in (1-2), for water inventory constraint condition, mends canal pumping plant in " reservoir-benefit library pumping plant " subsystem year
Year allows water lift total amount constraint condition, reservoir operation criterion constraint condition and reservoir capacity constraint condition；
(2) model solution
(2-1) data preparation was divided into N number of period for 1 year, and determines day part length；It is bent according to water level-storage capacity relationship
Line determines reservoir initial storage V0；Determine reservoir year for water inventory K, minimum capacity of a reservoir Vmin, the corresponding storage capacity V of flood controlP、
And utilizable capacity Vmin+Δ1；Measurement reservoir day part carrys out water LSi, evaporation with leakage EFi；Determine that mending library pumping plant year allows
Water lift total amount BK；Determine that mending canal pumping plant year allows water lift total amount BQ；Determine intake area day part under certain water frequency needs water
Measure Di(i=1,2 ..., N)；
(2-2) solves system model using Dynamic Programming successive approximation method, obtains day part reservoir optimal water supply mistake
Journey Xi, optimal abandoning process water PSi, mend the optimal rate of water make-up process Z of library Group of Pumping Stationi, and mend the optimal rate of water make-up process of canal pumping plant
Yi。
As shown in Figure 1-3, constraint condition in step (1-2) are as follows:
(1) " reservoir-Dan Buku pumping plant " system year is for water inventory constraint condition: in different level year, different fractions
In the case of, consider that intake area needs water requirement, the water that " reservoir-benefit library pumping plant " water supply project may provide；
In formula: K is the year of reservoir for water inventory, ten thousand m3；BK is to mend to allow water lift total amount, ten thousand m in library pumping plant year3；
(2) canal pumping plant year water lift total amount constraint condition is mended:
In formula: BQ is to mend to allow water lift total amount, ten thousand m in canal pumping plant year3；
(3) reservoir operation criterion constraint condition: according to reservoir water equilibrium equation,
Vi=Vi-1+LSi+Zi-PSi-EFi-Xi, (i=1,2 ..., N) (4)
(3-1) is lower than reservoir minimum capacity of a reservoir V when the i-th period end pondageminWhen, then the i-th period should be given by mending library pumping plant
Reservoir moisturizing, moisturizing to utilizable capacity (Δ more than reservoir minimum capacity of a reservoir1), it may be assumed that
Vi<VminWhen: Zi=Vmin-Vi+Δ1； (5)
This period reservoir abandons water PSi=0；
(3-2) when meeting with flood, the i-th period end pondage is greater than pondage corresponding to flood control
VPWhen, then the i-th period should carry out abandoning water to reservoir, abandon water to flood control by reservoir regulation limiting water level, it may be assumed that
Vi>VPWhen: PSi=Vi-VP； (6)
This period mends library pumping plant rate of water make-up Zi=0；
(3-3) is when the i-th period end pondage is between minimum capacity of a reservoir VminWith reservoir filling corresponding to flood control
Measure VPBetween, then the i-th period reservoir does not need to abandon water, mends library pumping plant and does not need moisturizing yet, it may be assumed that
Vmin≤Vi≤VPWhen: Zi=PSi=0； (7)
In formula: Vi、Vi-1Respectively reservoir i-th, the reservoir storage of i-1 period Mo, ten thousand m3；ZiFor the benefit library pumping plant of the i-th period
Rate of water make-up, ten thousand m3；LSi、PSi、EFiRespectively the i-th period of reservoir comes water, ten thousand m3, abandon water, ten thousand m3；Evaporation and leakage,
Ten thousand m3；Vmin、VPRespectively storage capacity corresponding to the minimum capacity of a reservoir of reservoir and flood control, ten thousand m3；
(4) reservoir capacity constraint condition: the pondage of day part should be between reservoir minimum capacity of a reservoir and flood control
Between corresponding storage capacity, it may be assumed that
Vmin≤Vi≤VP, (i=1,2 ..., N) (8)
As shown in Figure 1-3, the method solved using Dynamic Programming successive approximation method to system model in step (2-2) are as follows:
(1) it drafts and mends the initial rate of water make-up Y of canal pumping plant1, i(i=1,2 ..., N), it is desirable that meet and mend canal pumping plant year water lift total amount
Former intake area is regarded as one because mending at this point, deleting the benefit canal pumping plant subsystem from original system by constraint condition, i.e. formula (3)
A new intake area after the diminution of canal pumping plant moisturizing occurring area, the other compositions of system remain unchanged, i.e., sufficiently irrigate item for former
The conversion of the two station system Water Resources Allocation Model of a library of " reservoir-benefit library pumping plant-benefit canal pumping plant " under part, i.e. formula (1)~(8)
It is W for each stage water requirementi" reservoir-Dan Buku pumping plant " of the corresponding of no help canal pumping plant in new intake area of (i=1,2 ..., N)
System Water Resources Allocation Model:
Objective function:
Wherein: Wi=Di-Y1,i, due to day part Di、Y1,iIt is known that then WiIt is i.e. known；
Constraint condition:
(1-1) " reservoir-Dan Buku pumping plant " system year is for water inventory constraint condition: i.e. formula (2)；
(1-2) reservoir operation criterion constraint condition: i.e. formula (4)~(7)；
(1-3) reservoir capacity constraint condition: i.e. formula (8)；
(2) " reservoir-Dan Buku pumping plant " system water resource optimal allocation mathematical model, i.e. formula (2), (4)~(9) be with by
The minimum target of quadratic sum of the difference of pool day part water deficit, can divide by the stage of decision variable of day part reservoir yield
One-dimensional dynamic programming model obtain corresponding recurrence equation referring to one-dimensional dynamic programming evaluation principle are as follows:
(2-1) stage i=1:
g1(λ1)=min (X1,1-W1)2 (10)
In formula, state variable λ1Indicate preceding 1 water supply period reservoir yield, it can be discrete in corresponding feasible zone: λ1=
0,E1,E2,...,K+BK；To each discrete λ1, decision variable X1,1(reservoir yield) can be discrete in corresponding feasible zone, such as
00000 m3, 200,000 m3, 400,000 m3, 600,000 m3、…X1,1,maxDeng (X1,1,maxFor the 1st stage reservoir maximum water supply capacity), it should meet:
X1,1≥λ1；By X1,1Formula (10) are substituted into respectively, respectively obtain each discrete λ1When value, optimal X1,1And its corresponding g1(λ1)；
Then, according to formula (4), the 1st stage end reservoir capacity V1=V0+LS1-EF1-X1,1, not yet consider to mend library pump at this time
Stand moisturizing or reservoir abandons water, and Ying Caiyong formula (5)~(7) are tested and corrected:
(2-1-1) works as V1<Vmin, then consider to carry out moisturizing, Z to reservoir by mending library pumping plant1,1=Vmin-V1+Δ1, repair at this time
Positive storage capacity V1 *=Vmin+Δ1；
(2-1-2) works as V1> VP, then need to abandon dispatching requirement of the water to guarantee reservoir capacity, PS1,1=V1-VP, repair at this time
Positive storage capacity V1 *=VP；
(2-1-3) works as Vmin≤V1≤VP, then Z1,1=PS1,1=0, V at this time1 *=V1；
By step (2-1-1)~(2-1-3), corrects and determine the 1st stage end reservoir capacity V1 *, while can get and correspond to
Reservoir abandon water PS1,1Or mend library pumping plant rate of water make-up Z1,1；
(2-2) stage i=2,3 ..., N-1:
gi(λi)=min [(X1,i-Wi)2+gi-1(λi-1)] (11)
In formula, state variable λiFor the preceding i period reservoir for water inventory, equally carry out respectively discrete: λi=0, E1,
E2,...,K+BK；To each discrete λi, decision variable (reservoir yield X1,i) it is discrete ibid, and should meet:
State transition equation: λi-1=λi-X1,i (12)
In formula: i=2,3 ..., N-1
By each discrete X1,iValue substitutes into the (X in formula (11) respectively1,i-Wi)2, by state transition equation formula (12), to each
A discrete X1,i, search i-1 stage min gi-1(λi-1) value, thus to obtain (X1,i-Wi)2+min gi-1(λi-1), it is more than completion
All discrete X1,iAfter optimizing, final can get meets min [(X1,i-Wi)2+gi-1(λi-1)] require optimal X1,iProcess and its
Corresponding gi(λi)；
Equally, it also needs to determine the i-th stage end storage capacity V using formula (4)i, then tested, corrected using formula (5)~(7)
Determine the i-th period end storage capacity Vi *, while can get corresponding reservoir and abandoning process water PS1,iWith benefit library pumping plant rate of water make-up process Z1,i
(i=2,3 ..., N-1), the same step of process (2-1-1)~(2-1-3)；
(2-3) stage N:
gN(λN)=min [(X1,N-WN)2+gN-1(λN-1)] (13)
In formula, state variable λNFor the top n period reservoir for water inventory: λN=K+BK；Decision variable (reservoir yield
X1,N) equally discrete in corresponding feasible zone；
State transition equation: λN-1=λN-X1,N (14)
It is final to obtain objective function optimal value G using step (2-2) the method1=gN(λN) and corresponding reservoir
Optimal water supply process X1,i, abandon process water PS1,i, mend library pumping plant rate of water make-up process Z1,i(i=1,2 ..., N)；
(3) X for being calculated step (1)1,iProcess substitutes into master mould, i.e. formula (1)~(8) as input condition, at this point, by
In reservoir yield X1,iWith benefit library pumping plant rate of water make-up Z1,iIt is known that can should " reservoir-Dan Buku pumping plant " subsystem from former system
It is deleted in system, former intake area is regarded as one because one after the diminution of " reservoir-Dan Buku pumping plant " co-supplying occurring area is new
Intake area, the other compositions of system remain unchanged；It can be by " reservoir-benefit library pumping plant-benefit canal pumping plant " under former sufficiently irrigation conditions
It is U that one library, two station system Water Resources Allocation Model, i.e. formula (1)~(8), which are converted into each stage water requirement,i(i=1,2 ...,
N the corresponding single benefit canal pumping plant Water Resources Allocation Model in new intake area):
Objective function:
Wherein: Ui=Di-X1,i, due to day part Di、X1,iIt is known that then UiIt is i.e. known；
Constraint condition:
Mend canal pumping plant year water lift total amount constraint condition: i.e. formula (3)；
This singly mends canal pumping plant Water Resources Allocation Model, i.e., formula (3), (15) are with the difference of intake area day part water deficit
The minimum target of quadratic sum, the one-dimensional Dynamic Programming mould that can divide as the stage of decision variable of canal pumping plant rate of water make-up is mended using day part
Type obtains corresponding recurrence equation still referring to one-dimensional dynamic programming evaluation principle are as follows:
(3-1) stage i=1:
g1(α1)=min (Y2,1-U1)2 (16)
In formula, state variable α1Canal pumping plant water lift amount is mended for preceding 1 period, a fixed step size can be pressed in corresponding feasible zone
It is discrete: α1=0, F1,F2,...,BQ；To each discrete α1, decision variable Y2,1(mending canal pumping plant rate of water make-up) can be feasible in correspondence
It is discrete in domain, such as 00,000 m3, 200,000 m3, 400,000 m3, 600,000 m3、…Y2,1,maxDeng (Y2,1,maxCanal pumping plant maximum is mended for the 1st stage to mention
Outlet capacity), it should meet: Y2,1≥α1；By Y2,1Formula (16) are substituted into respectively, respectively obtain each discrete α1When value, optimal Y2,1And its
Corresponding g1(α1)；
(3-2) stage i=2,3 ... N-1:
gi(αi)=min [(Y2,i-Ui)2+gi-1(αi-1)] (17)
In formula: state variable αiCanal pumping plant water lift total amount is mended for the preceding i period, is equally carried out respectively discrete: αi=0, F1,
F2,...,BQ；To each discrete αi, decision variable Y2,i(reservoir yield) is discrete ibid, and should meet:
State transition equation: αi-1=αi-Y2,i (18)
By each discrete Y2,iValue substitutes into the (Y in formula (17) respectively2,i-Ui)2, by state transition equation formula (18), search i-
1 stage min gi-1(αi-1) value, thus to obtain (Y2,i-Ui)2+mingi-1(αi-1), complete all of above discrete Y2,iAfter optimizing,
Final can get meets min [(Y2,i-Ui)2+gi-1(αi-1)] require optimal Y2,iProcess and its corresponding gi(αi)；
(3-3) stage N:
gN(αN)=min [(Y2,N-UN)2+gN-1(αN-1)] (19)
In formula, state variable αNFor the benefit canal pumping plant water lift total amount of top n period, αN=BQ；Decision variable Y2,N(mend canal pump
Stand rate of water make-up) it is equally discrete in corresponding feasible zone；
State transition equation: αN-1=αN-Y2,N (20)
Using step (3-2) the method, final obtain meets the αNIt is required that the optimal rate of water make-up process Y of benefit canal pumping plant2,i
(i=1,2 ..., N), and objective function optimal value G2=gN(αN)；
(4) by the Y of acquisition2,iProcess substitutes into master mould and to be again converted into " reservoir-Dan Buku pumping plant " system water resource excellent
Change configuration mathematical model, i.e. formula (2), (4)~(9), repeats step (1) and solve, obtain new reservoir yield process X3,i, abandon
Process water PS3,i, mend library pumping plant rate of water make-up process Z3,iAnd objective function optimal value G3；
(5) by the X of acquisition3,iProcess substitutes into master mould again and is converted into single benefit canal pumping plant Water Resources Allocation Model,
That is formula (3), (15) repeat step (2) and solve, obtain new benefit canal pumping plant rate of water make-up process Y4,iAnd objective function optimal value
G4；
(6) above step is repeated, the objective function optimal value obtained twice up to date meets | (Gn-Gn-1)/Gn|≤ε(ε
For iteration control precision), then iteration terminates；The G that will be obtained for the last time at this timenIt is corresponding each as objective function optimal value
Period reservoir optimal water supply process Xn-1,i, optimal abandoning process water PSn-1,i, mend the optimal rate of water make-up process Z of library pumping plantn-1,i, with
And mend the optimal rate of water make-up process Y of canal pumping plantn,iStrategy be original system water resource combined dispatching optimal policy.
Claims (3)
1. abundant irrigation conditions lower storage reservoir-benefit library pumping plant-benefit canal pumping station system water resource optimal allocation method, which is characterized in that
The single reservoir and 1 benefit canal pumping plant for mending the moisturizing of library pumping plant from 1 combine the co-supplying to supply water to intake area, including following step
It is rapid:
(1) model construction
(1-1) is total with day part intake area water requirement in reservoir-benefit library pumping plant-benefit canal pumping plant 2 station system year of one library and system
The least square of water supply difference and be target, establish following objective function:
In formula: G be research object year in day part the difference for water requirement least square and；When S is each in research object year
The difference for water requirement of section, ten thousand m3；N is the when number of segment divided in year；Segment number when i is (i=1,2 ..., N)；XiFor reservoir i-th
The water supply of period, ten thousand m3；YiFor the rate of water make-up for mending the i-th period of canal pumping plant, ten thousand m3；DiFor the water requirement of the i-th period of intake area,
Ten thousand m3；Objective function is inclined between the system water supply amount and intake area water requirement in order to accelerate to reduce using quadratic sum expression
Difference；
Constraint condition is arranged in (1-2), and reservoir-benefit library pumping plant subsystem year, mending allowed in canal pumping plant year for water inventory constraint condition
Water lift total amount constraint condition, reservoir operation criterion constraint condition and reservoir capacity constraint condition；
(2) model solution
(2-1) data preparation was divided into N number of period for 1 year, and determines day part length；It is true according to water level-capacity curve
Determine reservoir initial storage V0；Determine reservoir year for water inventory K, minimum capacity of a reservoir Vmin, the corresponding storage capacity V of flood controlPAnd
Utilizable capacity Vmin+Δ1；Measurement reservoir day part carrys out water LSi, evaporation with leakage EFi；Determine that mending library pumping plant year allows water lift
Total amount BK；Determine that mending canal pumping plant year allows water lift total amount BQ；Determine the water requirement D of intake area day part under certain water frequencyi
(i=1,2 ..., N)；
(2-2) solves system model using Dynamic Programming successive approximation method, obtains day part reservoir optimal water supply process Xi,
Optimal abandoning process water PSi, mend the optimal rate of water make-up process Z of library Group of Pumping Stationi, and mend the optimal rate of water make-up process Y of canal pumping planti。
2. abundant irrigation conditions lower storage reservoir-benefit library pumping plant according to claim 1-benefit canal pumping station system water resource optimization is matched
Set method, which is characterized in that constraint condition described in step (1-2) are as follows:
(1) reservoir-Dan Buku pumping station system year is for water inventory constraint condition: at different level year, different fractions,
Consider that intake area needs water requirement, the water that reservoir-benefit library pumping plant water supply project may provide；
In formula: K is the year of reservoir for water inventory, ten thousand m3；BK is to mend to allow water lift total amount, ten thousand m in library pumping plant year3；
(2) canal pumping plant year water lift total amount constraint condition is mended:
In formula: BQ is to mend to allow water lift total amount, ten thousand m in canal pumping plant year3；
(3) reservoir operation criterion constraint condition: according to reservoir water equilibrium equation,
Vi=Vi-1+LSi+Zi-PSi-EFi-Xi, (i=1,2 ..., N) (4)
(3-1) is lower than reservoir minimum capacity of a reservoir V when the i-th period end pondageminWhen, then the i-th period should be by mending library pumping plant to reservoir
Moisturizing, moisturizing to utilizable capacity (Δ more than reservoir minimum capacity of a reservoir1), it may be assumed that
Vi<VminWhen: Zi=Vmin-Vi+Δ1； (5)
This period reservoir abandons water PSi=0；
(3-2) when meeting with flood, the i-th period end pondage is greater than pondage V corresponding to flood controlPWhen,
Then the i-th period should carry out abandoning water to reservoir, abandon water to flood control by reservoir regulation limiting water level, it may be assumed that
Vi>VPWhen: PSi=Vi-VP； (6)
This period mends library pumping plant rate of water make-up Zi=0；
(3-3) is when the i-th period end pondage is between minimum capacity of a reservoir VminWith pondage V corresponding to flood controlP
Between, then the i-th period reservoir does not need to abandon water, mends library pumping plant and does not need moisturizing yet, it may be assumed that
Vmin≤Vi≤VPWhen: Zi=PSi=0； (7)
In formula: Vi、Vi-1Respectively reservoir i-th, the reservoir storage of i-1 period Mo, ten thousand m3；ZiFor the benefit library pumping plant moisturizing of the i-th period
Amount, ten thousand m3；LSi、PSi、EFiRespectively the i-th period of reservoir comes water, ten thousand m3, abandon water, ten thousand m3；Evaporation and leakage, ten thousand
m3；Vmin、VPRespectively storage capacity corresponding to the minimum capacity of a reservoir of reservoir and flood control, ten thousand m3；
(4) reservoir capacity constraint condition: the pondage of day part should be corresponding between reservoir minimum capacity of a reservoir and flood control
Storage capacity between, it may be assumed that
Vmin≤Vi≤VP, (i=1,2 ..., N) (8).
3. abundant irrigation conditions lower storage reservoir-benefit library pumping plant according to claim 1-benefit canal pumping station system water resource optimization is matched
Set method, which is characterized in that the method solved using Dynamic Programming successive approximation method to system model described in step (2-2)
Are as follows:
(1) it drafts and mends the initial rate of water make-up Y of canal pumping plant1, i(i=1,2 ..., N), it is desirable that meet and mend the constraint of canal pumping plant year water lift total amount
Former intake area is regarded as one because mending canal pump at this point, deleting the benefit canal pumping plant subsystem from original system by condition, i.e. formula (3)
Moisturizing occurring area of standing reduce after a new intake area, the other compositions of system remain unchanged, i.e., will be under former sufficiently irrigation conditions
One library of reservoir-benefit library pumping plant-benefit canal pumping plant, two station system Water Resources Allocation Model, i.e. formula (1)~(8) are converted into each rank
Section water requirement is WiReservoir-Dan Buku pumping station system water the money of the corresponding of no help canal pumping plant in new intake area of (i=1,2 ..., N)
Source optimization allocation models:
Objective function:
Wherein: Wi=Di-Y1,i, due to day part Di、Y1,iIt is known that then WiIt is i.e. known；
Constraint condition:
(1-1) reservoir-Dan Buku pumping station system year is for water inventory constraint condition: i.e. formula (2)；
(1-2) reservoir operation criterion constraint condition: i.e. formula (4)~(7)；
(1-3) reservoir capacity constraint condition: i.e. formula (8)；
(2) reservoir-Dan Buku pumping station system water resource optimal allocation mathematical model, i.e. formula (2), (4)~(9) are each with intake area
The minimum target of the quadratic sum of the difference of period water deficit can divide one-dimensional using day part reservoir yield as the stage of decision variable
Dynamic programming model obtains corresponding recurrence equation referring to one-dimensional dynamic programming evaluation principle are as follows:
(2-1) stage i=1:
g1(λ1)=min (X1,1-W1)2 (10)
In formula, state variable λ1Indicate preceding 1 water supply period reservoir yield, it can be discrete in corresponding feasible zone: λ1=0,
E1,E2,...,K+BK；To each discrete λ1, decision variable X1,1(reservoir yield) can be discrete in corresponding feasible zone, and such as 0
Ten thousand m3, 200,000 m3, 400,000 m3, 600,000 m3、…X1,1,maxDeng (X1,1,maxFor the 1st stage reservoir maximum water supply capacity), it should meet:
X1,1≥λ1；By X1,1Formula (10) are substituted into respectively, respectively obtain each discrete λ1When value, optimal X1,1And its corresponding g1(λ1)；
Then, according to formula (4), the 1st stage end reservoir capacity V1=V0+LS1-EF1-X1,1, not yet consider to mend the moisturizing of library pumping plant at this time
Or reservoir abandons water, Ying Caiyong formula (5)~(7) are tested and are corrected:
(2-1-1) works as V1<Vmin, then consider to carry out moisturizing, Z to reservoir by mending library pumping plant1,1=Vmin-V1+Δ1, library is corrected at this time
Hold V1 *=Vmin+Δ1；
(2-1-2) works as V1> VP, then need to abandon dispatching requirement of the water to guarantee reservoir capacity, PS1,1=V1-VP, storage capacity is corrected at this time
V1 *=VP；
(2-1-3) works as Vmin≤V1≤VP, then Z1,1=PS1,1=0, V at this time1 *=V1；
By step (2-1-1)~(2-1-3), corrects and determine the 1st stage end reservoir capacity V1 *, while can get corresponding water
Water PS is abandoned in library1,1Or mend library pumping plant rate of water make-up Z1,1；
(2-2) stage i=2,3 ..., N-1:
gi(λi)=min [(X1,i-Wi)2+gi-1(λi-1)] (11)
In formula, state variable λiFor the preceding i period reservoir for water inventory, equally carry out respectively discrete: λi=0, E1,E2,...,K
+BK；To each discrete λi, decision variable (reservoir yield X1,i) it is discrete ibid, and should meet:
State transition equation: λi-1=λi-X1,i (12)
In formula: i=2,3 ..., N-1
By each discrete X1,iValue substitutes into the (X in formula (11) respectively1,i-Wi)2, by state transition equation formula (12), to each from
Scattered X1,i, search i-1 stage mingi-1(λi-1) value, thus to obtain (X1,i-Wi)2+min gi-1(λi-1), complete it is all of above from
Scattered X1,iAfter optimizing, final can get meets min [(X1,i-Wi)2+gi-1(λi-1)] require optimal X1,iProcess and its corresponding
gi(λi)；
Equally, it also needs to determine the i-th stage end storage capacity V using formula (4)i, then tested using formula (5)~(7), amendment determines
I-th period end storage capacity Vi *, while can get corresponding reservoir and abandoning process water PS1,iWith benefit library pumping plant rate of water make-up process Z1,i(i=
2,3 ..., N-1), the same step of process (2-1-1)~(2-1-3)；
(2-3) stage N:
gN(λN)=min [(X1,N-WN)2+gN-1(λN-1)] (13)
In formula, state variable λNFor the top n period reservoir for water inventory: λN=K+BK；Decision variable (reservoir yield X1,N)
It is equally discrete in corresponding feasible zone；
State transition equation: λN-1=λN-X1,N (14)
It is final to obtain objective function optimal value G using step (2-2) the method1=gN(λN) and corresponding reservoir it is optimal
Water supply process X1,i, abandon process water PS1,i, mend library pumping plant rate of water make-up process Z1,i(i=1,2 ..., N)；
(3) X for being calculated step (1)1,iProcess substitutes into master mould, i.e. formula (1)~(8) as input condition, at this point, due to water
Reservoir yield X1,iWith benefit library pumping plant rate of water make-up Z1,iIt is known that can delete the reservoir-Dan Buku pumping plant subsystem from original system
It removes, former intake area is regarded as one because of a new intake area after the diminution of reservoir-Dan Buku pumping plant co-supplying occurring area, is
Other compositions of uniting remain unchanged；It can be by the two station system of a library of former sufficiently irrigation conditions lower storage reservoir-benefit library pumping plant-benefit canal pumping plant
It is U that Water Resources Allocation Model, i.e. formula (1)~(8), which are converted into each stage water requirement,iThe new intake area of (i=1,2 ..., N)
Corresponding one single benefit canal pumping plant Water Resources Allocation Model:
Objective function:
Wherein: Ui=Di-X1,i, due to day part Di、X1,iIt is known that then UiIt is i.e. known；
Constraint condition:
Mend canal pumping plant year water lift total amount constraint condition: i.e. formula (3)；
This singly mends canal pumping plant Water Resources Allocation Model, i.e., formula (3), (15) are with the flat of the difference of intake area day part water deficit
Square and minimum target mends the one-dimensional dynamic programming model that canal pumping plant rate of water make-up can divide as the stage of decision variable using day part,
Still referring to one-dimensional dynamic programming evaluation principle, corresponding recurrence equation is obtained are as follows:
(3-1) stage i=1:
g1(α1)=min (Y2,1-U1)2 (16)
In formula, state variable α1Canal pumping plant water lift amount is mended for preceding 1 period, it can be discrete by a fixed step size in corresponding feasible zone:
α1=0, F1,F2,...,BQ；To each discrete α1, decision variable Y2,1(mending canal pumping plant rate of water make-up) can be in corresponding feasible zone
It is discrete, such as 00,000 m3, 200,000 m3, 400,000 m3, 600,000 m3、…Y2,1,maxDeng (Y2,1,maxCanal pumping plant maximum water lift energy is mended for the 1st stage
Power), it should meet: Y2,1≥α1；By Y2,1Formula (16) are substituted into respectively, respectively obtain each discrete α1When value, optimal Y2,1And its it is corresponding
G1(α1)；
(3-2) stage i=2,3 ... N-1:
gi(αi)=min [(Y2,i-Ui)2+gi-1(αi-1)] (17)
In formula: state variable αiCanal pumping plant water lift total amount is mended for the preceding i period, is equally carried out respectively discrete: αi=0, F1,
F2,...,BQ；To each discrete αi, decision variable Y2,i(reservoir yield) is discrete ibid, and should meet:
State transition equation: αi-1=αi-Y2,i (18)
By each discrete Y2,iValue substitutes into the (Y in formula (17) respectively2,i-Ui)2, by state transition equation formula (18), search i-1 rank
Section mingi-1(αi-1) value, thus to obtain (Y2,i-Ui)2+mingi-1(αi-1), complete all of above discrete Y2,iAfter optimizing, finally
It can get and meet min [(Y2,i-Ui)2+gi-1(αi-1)] require optimal Y2,iProcess and its corresponding gi(αi)；
(3-3) stage N:
gN(αN)=min [(Y2,N-UN)2+gN-1(αN-1)] (19)
In formula, state variable αNFor the benefit canal pumping plant water lift total amount of top n period, αN=BQ；Decision variable Y2,NIt (mends canal pumping plant to mend
Water) it is equally discrete in corresponding feasible zone；
State transition equation: αN-1=αN-Y2,N (20)
Using step (3-2) the method, final obtain meets the αNIt is required that the optimal rate of water make-up process Y of benefit canal pumping plant2,i(i=
1,2 ..., N), and objective function optimal value G2=gN(αN)；
(4) by the Y of acquisition2,iProcess substitutes into master mould and is again converted into reservoir-Dan Buku pumping station system water resource optimal allocation
Mathematical model, i.e. formula (2), (4)~(9) repeat step (1) and solve, obtain new reservoir yield process X3,i, abandon water mistake
Journey PS3,i, mend library pumping plant rate of water make-up process Z3,iAnd objective function optimal value G3；
(5) by the X of acquisition3,iProcess substitutes into master mould again and is converted into single benefit canal pumping plant Water Resources Allocation Model, i.e. formula
(3), (15) repeat step (2) and solve, obtain new benefit canal pumping plant rate of water make-up process Y4,iAnd objective function optimal value G4；
(6) above step is repeated, the objective function optimal value obtained twice up to date meets | (Gn-Gn-1)/Gn(ε is repeatedly to |≤ε
Generation control precision), then iteration terminates；The G that will be obtained for the last time at this timenAs objective function optimal value, corresponding day part
Reservoir optimal water supply process Xn-1,i, optimal abandoning process water PSn-1,i, mend the optimal rate of water make-up process Z of library pumping plantn-1,i, and mend
The optimal rate of water make-up process Y of canal pumping plantn,iStrategy be original system water resource combined dispatching optimal policy.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910061144.0A CN109829580A (en) | 2019-01-23 | 2019-01-23 | Abundant irrigation conditions lower storage reservoir-benefit library pumping plant-benefit canal pumping station system water resource optimal allocation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910061144.0A CN109829580A (en) | 2019-01-23 | 2019-01-23 | Abundant irrigation conditions lower storage reservoir-benefit library pumping plant-benefit canal pumping station system water resource optimal allocation method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109829580A true CN109829580A (en) | 2019-05-31 |
Family
ID=66861852
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910061144.0A Pending CN109829580A (en) | 2019-01-23 | 2019-01-23 | Abundant irrigation conditions lower storage reservoir-benefit library pumping plant-benefit canal pumping station system water resource optimal allocation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109829580A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110675281A (en) * | 2019-09-26 | 2020-01-10 | 扬州大学 | Water resource optimization scheduling method for 'long vine and melon-bearing' type multi-reservoir system in hilly area |
| CN112214885A (en) * | 2020-09-28 | 2021-01-12 | 西安理工大学 | Irrigation area crop irrigation water quantity optimal distribution method under insufficient irrigation condition |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106295893A (en) * | 2016-08-12 | 2017-01-04 | 扬州大学 | Fully directly mend single pumping plant list water reservoir system water resource optimal allocation method in storehouse under irrigation conditions |
| CN106327065A (en) * | 2016-08-12 | 2017-01-11 | 扬州大学 | Water resource optimization configuration method for single pumping station - single reservoir system for direct canal supplement under full irrigation condition |
| CN107491635A (en) * | 2017-07-26 | 2017-12-19 | 华北电力大学 | A kind of step reservoir water sand joint optimal operation method based on nested dimension-reduction algorithm |
| CN108197769A (en) * | 2017-10-19 | 2018-06-22 | 扬州大学 | Single library-multiple station systems water resource optimal allocation the method in library is directly mended under a kind of abundant irrigation conditions |
-
2019
- 2019-01-23 CN CN201910061144.0A patent/CN109829580A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106295893A (en) * | 2016-08-12 | 2017-01-04 | 扬州大学 | Fully directly mend single pumping plant list water reservoir system water resource optimal allocation method in storehouse under irrigation conditions |
| CN106327065A (en) * | 2016-08-12 | 2017-01-11 | 扬州大学 | Water resource optimization configuration method for single pumping station - single reservoir system for direct canal supplement under full irrigation condition |
| CN107491635A (en) * | 2017-07-26 | 2017-12-19 | 华北电力大学 | A kind of step reservoir water sand joint optimal operation method based on nested dimension-reduction algorithm |
| CN108197769A (en) * | 2017-10-19 | 2018-06-22 | 扬州大学 | Single library-multiple station systems water resource optimal allocation the method in library is directly mended under a kind of abundant irrigation conditions |
Non-Patent Citations (1)
| Title |
|---|
| 龚志浩 等: "山湖水库"一库两站"联合运行优化调度方法", 《排灌机械工程学报》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110675281A (en) * | 2019-09-26 | 2020-01-10 | 扬州大学 | Water resource optimization scheduling method for 'long vine and melon-bearing' type multi-reservoir system in hilly area |
| CN110675281B (en) * | 2019-09-26 | 2023-06-27 | 扬州大学 | Water resource optimization scheduling method for 'long vine melon' type multi-reservoir system in hilly area |
| CN112214885A (en) * | 2020-09-28 | 2021-01-12 | 西安理工大学 | Irrigation area crop irrigation water quantity optimal distribution method under insufficient irrigation condition |
| CN112214885B (en) * | 2020-09-28 | 2024-04-26 | 西安理工大学 | Irrigation area crop irrigation water quantity optimal distribution method under insufficient irrigation condition |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106295893A (en) | Fully directly mend single pumping plant list water reservoir system water resource optimal allocation method in storehouse under irrigation conditions | |
| CN105096004B (en) | A kind of multi-reservoir supplies water transfer system real-time scheduling method | |
| WO2019006733A1 (en) | Long-term joint peak regulation dispatching method for trans-provincial interconnected hydropower station cluster | |
| CN106327065A (en) | Water resource optimization configuration method for single pumping station - single reservoir system for direct canal supplement under full irrigation condition | |
| CN104063808B (en) | Trans-provincial power transmission cascade hydropower station group peak-shaving dispatching two-phase search method | |
| CN103745023A (en) | Coupling modeling method for hydropower station power generated output scheme making and optimal load distribution | |
| CN104268653A (en) | Cascade reservoir optimal scheduling method based on ESP | |
| CN106228276A (en) | Single pumping plant list water reservoir system water resource optimal allocation method of canal is directly mended under the conditions of insufficient irrigation | |
| CN104504455B (en) | A kind of lower GROUP OF HYDROPOWER STATIONS Long-term Optimal Dispatch method of step accumulation of energy control | |
| CN104047258A (en) | Environmental-protection-oriented ecological reservoir capacity determination method for combined scheduling reservoir group | |
| Li et al. | Water resources optimal allocation based on large-scale reservoirs in the upper reaches of Yangtze River | |
| CN106223394B (en) | The mono- water reservoir system water resource optimal allocation method of single pumping plant-in library is directly mended under the conditions of insufficient irrigation | |
| CN109829580A (en) | Abundant irrigation conditions lower storage reservoir-benefit library pumping plant-benefit canal pumping station system water resource optimal allocation method | |
| CN107798471B (en) | More libraries-multiple station systems water resource optimal allocation method of canal is directly mended under a kind of abundant irrigation conditions | |
| CN112581310A (en) | Cascade hydropower station group power generation index distribution method | |
| CN103093282A (en) | Hydropower station group terminal energy storage maximum short-term optimized dispatching method | |
| CN111864811B (en) | Photovoltaic capacity configuration method of water-light complementary system considering power grid output demand | |
| CN111915065A (en) | River dry season multi-target dynamic water resource optimal configuration system and method | |
| CN108197769A (en) | Single library-multiple station systems water resource optimal allocation the method in library is directly mended under a kind of abundant irrigation conditions | |
| CN109002946A (en) | A kind of " station of two libraries-two " system water resources optimal operation method of river and lake moisturizing | |
| CN114239992A (en) | Method for formulating dynamic flood limit water level of reservoir | |
| CN106524277A (en) | Multi-energy region energy supply system for winter heating | |
| CN116402410B (en) | Distributed water quantity and water quality configuration method and system | |
| CN105160443A (en) | Optimal complex reservoir group dispatching method based on extended linear quadratic Gaussian method | |
| CN107748930A (en) | Single storehouse multiple station systems water resource optimal allocation method of canal is directly mended under a kind of fully irrigation conditions |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190531 |
|
| RJ01 | Rejection of invention patent application after publication |