CN109002946A - A kind of " station of two libraries-two " system water resources optimal operation method of river and lake moisturizing - Google Patents

A kind of " station of two libraries-two " system water resources optimal operation method of river and lake moisturizing Download PDF

Info

Publication number
CN109002946A
CN109002946A CN201811216569.6A CN201811216569A CN109002946A CN 109002946 A CN109002946 A CN 109002946A CN 201811216569 A CN201811216569 A CN 201811216569A CN 109002946 A CN109002946 A CN 109002946A
Authority
CN
China
Prior art keywords
water
pumping plant
moisturizing
period
library
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.)
Granted
Application number
CN201811216569.6A
Other languages
Chinese (zh)
Other versions
CN109002946B (en
Inventor
龚志浩
程吉林
龚懿
蒋晓红
陈兴
汪靓
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yangzhou University
Original Assignee
Yangzhou University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Yangzhou University filed Critical Yangzhou University
Priority to CN201811216569.6A priority Critical patent/CN109002946B/en
Publication of CN109002946A publication Critical patent/CN109002946A/en
Application granted granted Critical
Publication of CN109002946B publication Critical patent/CN109002946B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Human Resources & Organizations (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Marketing (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Water Supply & Treatment (AREA)
  • Development Economics (AREA)
  • Game Theory and Decision Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Feedback Control In General (AREA)

Abstract

The invention belongs to Water Resources Irrigations to distribute technical field rationally, more particularly to a kind of " station of two libraries-two " system water resources optimal operation method of river and lake moisturizing, specific step is as follows: (1) seeking advice from by actual measurement or to irrigated areas administration department, collect the related data needed for calculating；(2) " station of two libraries-two " system water resources optimal operation model decomposition is high basin and its moisturizing pumping station system and low library and its moisturizing pumping station system two subsystems water resources optimal operation model by system decomposition；(3) subsystem optimizes；(4) system coordination, the final Optimized Operation scheme for obtaining " two libraries, two station " system.The present invention establishes corresponding water resources optimal operation mathematical model, it is solved on the basis of composition decomposition method using Dynamic Programming, it can get minimum water deficit in intake area in certain delivery period, corresponding each reservoir optimal water supply, abandon process water, and each pumping plant mentions process water, " station of two libraries-two " system Seasonal Water-Deficit is solved, the utilization rate of water resource is improved.

Description

A kind of " station of two libraries-two " system water resources optimal operation method of river and lake moisturizing
Technical field
The present invention relates to " station of two libraries-two " system water resources optimal operation methods of river and lake moisturizing, belong to Water Resources Irrigation Distribute technical field rationally.
Background technique
Currently since water resource spatial and temporal distributions are uneven, the socio-economic development in many areas has been restricted.The southeastern coastal areas Domestic Hills are in subtropical zone or monsoon climate of medium latitudes area, and flood season runoff is concentrated, and often result in flood；Remaining when Phase arid will lead to agricultural water difficulty without rain.In order to cope with the problem of water resource is unevenly distributed in time, locality labour The people build the storing facilities of waters such as reservoir, small reservoir and harvest flood season rainwater, need water requirement to meet other periods.Due to local mistake Border water resource is abundant, and this kind of water storage project is configured with toward contact draws water lift engineering, can mend library by water lift when necessary.So And there are also the higher reservoir of some topographies or small reservoirs since catchment area is smaller and water resource supplement of not passing by, in schedule periods Occur just needing to mend library from other large and middle reservoirs diversions when water shortage.
" station of two libraries-two " water resource system of river and lake moisturizing is the typical irrigation system of southern hilly area, although water source Sufficient condition, but since system does not unify reliable combined operating mechanism, the water resource of internal system cannot be adjusted reasonably Match, the contradiction for causing generation abandoning water, water shortage, moisturizing to coexist.
Summary of the invention
Abandoning water, water shortage, moisturizing caused by " station of the two libraries-two " water resource system of the present invention for river and lake moisturizing coexist Contradiction, establish corresponding water resources optimal operation mathematical model, on the basis of composition decomposition method use Dynamic Programming It solves, can get minimum water deficit in intake area in certain delivery period, corresponding each reservoir optimal water supply, abandons process water, with And each pumping plant mentions process water, solves " station of two libraries-two " system Seasonal Water-Deficit, improves the utilization rate of water resource.
The present invention program is as follows: a kind of " station of two libraries-two " system water resources optimal operation method of river and lake moisturizing.System It being made of two reservoirs and two pumping plants, two reservoirs include low library and high basin, respectively supply water to independent water user, two pumps It stands including low library moisturizing pumping plant and high basin moisturizing pumping plant, low library moisturizing pumping plant water lift from external river supplements low library, and high basin is mended Pump works supplements high basin from low library water lift, and system schematic is shown in Fig. 1.Specific step is as follows for system optimization dispatching method:
(1) it is seeked advice from by actual measurement or to irrigated areas administration department, collects the related data needed for calculating, specifically include:
1. irrigated area reservoir operation usually divide when number of segment T and corresponding day part length Ni；
2. represent reservoir day part in year carrys out water LS1,i, LS2,i(i=1,2 ..., T)；
3. representing the water requirement YS of irrigated area day part in year1,i, YS2,i(i=1,2 ..., T)；
4. the characteristic data of two reservoirs: reservoir initial storage V1,0, V2,0, storage capacity boundWith
5. the characteristic data of two pumping plants: design water lift flow Q1, Q2With total water lift amount data of actual year；
6. the water rights data of local water administrative sector distribution, refer mainly to system allows water lift total from the year of external river diversion Measure BZ1。
(2) " station of two libraries-two " system water resources optimal operation model decomposition is high basin and its small pump by system decomposition System of standing and low library and its moisturizing pumping station system two subsystems water resources optimal operation model.
1. high basin and its moisturizing pumping plant
Objective function:
In formula: F2For high basin and its day part water deficit quadratic sum in moisturizing pumping plant subsystem year；X2,iFor high basin and its benefit The water supply of i-th period of pump works subsystem；YS2,iFor high basin and its water requirement of moisturizing pumping plant the i-th period of subsystem.
Constraint condition:
Y2,i≤Q2×Ni (3)
V2,i=V2,i-1+LS2,i-X2,i+Y2,i-PS2,i-EF2,i (4)
V2,i m≤V2,i≤V2,i M (5)
In formula: LS2,iFor high basin period i two Phase flow amount；BZ2To allow water lift amount in the year of high basin moisturizing pumping plant； Y2,iThe water lift amount for being high basin moisturizing pumping plant within the i period；Q2For the design discharge of high basin moisturizing pumping plant；NiFor the length of period i； V2,iAnd V2,i-1Respectively reservoir storage of the high basin in i and i-1 period Mo；PS2,iFor high basin the i period abandoning water；EF2,iFor height Water loss of the library in the i period；For high basin the i period reservoir storage lower limit；It is high basin in the reservoir storage of i period Limit.
2. low library and its moisturizing pumping plant
Objective function:
In formula: F1For low library and its day part water deficit quadratic sum in moisturizing pumping plant subsystem year；X1,iFor low library and its benefit The water supply of i-th period of pump works subsystem；YS1,iFor low library and its water requirement of moisturizing pumping plant the i-th period of subsystem；Y2,iFor The water that high basin moisturizing pumping plant is extracted from low library within the i period.
Constraint condition:
Y1,i≤Q1×Ni (8)
V1,i=V1,i-1+LS1,i-X1,i+Y1,i-PS1,i-EF1,i (9)
V1,i m≤V1,i≤V1,i M (10)
In formula: BZ1To allow water lift amount, i.e., the water rights of local water administrative sector distribution low library moisturizing pumping plant year；Y1,iIt is low Water lift amount of the moisturizing pumping plant in library within the i period；Q1For the design discharge of low library moisturizing pumping plant；V1,iAnd V1,i-1Respectively low library is in i With the reservoir storage of i-1 period Mo；PS1,iFor low library the i period abandoning water；EF1,iFor low library the i period water loss； For low library the i period reservoir storage lower limit；For low library the i period the reservoir storage upper limit.
(3) subsystem optimizes
With the permission water lift amount BZ of high basin moisturizing pumping plant2For coordination variable, its practical year water lift total amount before being not optimised is enabled It is discrete by the step-length d for meeting system call required precision between the practical year total water shortage of high basin, to high basin and its moisturizing pumping plant Subsystem optimizes, and obtains each stage water deficit quadratic sum F in high basin year2, water supply process X2,i, abandon process water PS2,i And its moisturizing pumping plant mentions process water Y2,i；By the water lift amount Y of high basin moisturizing pumping plant2,iAs one of water requirement substitute into low library and The Optimized model of its moisturizing pumping station system obtains each stage water deficit quadratic sum F in low library year1, water supply process X1,i, abandon water Amount process PS1,iAnd its moisturizing pumping plant mentions process water Y1,i。
(4) system coordination
With stage water deficit quadratic sum F (F=F each in system year1+F2) it is dependent variable, allow to mention with high basin moisturizing pumping plant Water BZ2Regression analysis is carried out for independent variable, determines optimal F and BZ2The scheduling of value and its corresponding two reservoirs and pumping plant Process, the final Optimized Operation scheme for obtaining " two libraries, two station " system.
Preferably, the step (3) specifically includes the following steps:
(1) stage i=1:
f1(λ1)=min (X1-YS1)2 (11)
In formula, λ1For state variable, indicate preceding 1 water supply period reservoir yield, in correspondence feasible zone step-length d from It dissipates: λ1=0, W1,W2,…,∑LSi+BZ；To each discrete λ1, decision variable X1[0, λ in corresponding feasible zone1] press step-length d It is discrete, it determines correspond to each λ respectively1Value, optimal X1And its corresponding period minimum water deficit quadratic sum f1(λ1)。
Then, the 1st stage end pondage V1=V0+LS1-EF1-X1, not yet consider that pumping plant moisturizing or reservoir are abandoned at this time Water should be tested and be corrected:
1. if V1<V1 m, then Y1=min (V1 m-V1, Q × N1), PSi=0；
2. if V1>V1 M, then Y1=0, PS1=V1-V1 M；
3. if V1 m≤V1≤V1 M, then Y1=PS1=0；
4. revised pondage is V1'=V1+Y1-PS1。
By step 1.~4., correct and simultaneously determine the 1st stage end pondage, while can get corresponding reservoir and abandoning water Measure PS1Or pumping plant rate of water make-up Y1。
(2) stage i=2,3 ... T-1:
fi(λi)=min [(Xi-YSi)2+fi-1(λi-1)] (12)
In formula, state variable λiFor the preceding i period reservoir for water inventory, equally it is carried out respectively by step-length d discrete: λi=0, W1,W2,…,∑LSi+BZ；To each discrete λi, decision variable Xi[0, λ in corresponding feasible zonei] press step-length d It is discrete.
State transition equation:
λi-1=λi-Xi (13)
In formula: i=2,3 ..., T-1
To each discrete λi, by each discrete XiValue substitutes into the (X in formula (12) respectivelyi-YSi)2, by state transfer side Formula (13), to each discrete Xi, i-1 stage minimum f before searchingi-1(λi-1) value, it thus can get (Xi-YSi)2+min fi-1(λi-1), complete all of above discrete XiAfter optimizing, final can get meets min [(Xi-YSi)2+fi-1(λi-1)] require Preceding i period system minimum water deficit quadratic sum fi(λi) value and its corresponding each stage reservoir optimal water supply Xi。
Equally, it also needs to determine the i-th stage end pondage Vi, Vi=Vi-1+LSi-EFi-Xi, then test and repair Just
1. if Vi<Vi m, then Yi=min (Vi m-Vi, Q × Ni), PSi=0；
2. if Vi>Vi M, then Yi=0, PSi=Vi-Vi M；
3. if Vi m≤Vi≤Vi M, then Yi=PSi=0；
4. revised pondage is Vi'=Vi+Yi-PSi。
By step 1.~4., correct and simultaneously determine the i-th period end pondage, while can get corresponding reservoir and abandoning water Amount process PSiAnd pumping plant rate of water make-up process Yi。
(3) stage T:
fT(λT)=min [(XT-YST)2+fT-1(λT-1)] (14)
State variable λTFor the preceding T period reservoir for water inventory, equally it is carried out respectively discrete: λT=0, W1, W2,…,∑LSi+BZ；To each discrete λT, decision variable XT[0, λ in corresponding feasible zoneT] discrete by step-length d.
State transition equation:
λT-1=λT-XT (15)
To each discrete λT, by each discrete XTValue substitutes into the (X in formula (14) respectivelyT-YST)2, by state transfer side Formula (15), to each discrete XT, T-1 stage minimum f before searchingT-1(λT-1) value, it thus can get (XT-YST)2+min fT-1(λT-1), complete all of above discrete XTAfter optimizing, final can get meets min [(XT-YST)2+fT-1(λT-1)] require Preceding T period system minimum water deficit quadratic sum fT(λT) value, it finally obtains and meets the λTIt is required that the optimal water supply process of reservoir Xi, i=1,2 ..., T, corresponding reservoir abandoning process water PSi, i=1,2 ..., T, pumping plant rate of water make-up process Yi, i=1,2 ..., T and master mould objective function optimal value F*=min fT(λT)。
Preferably, the step (4) specifically includes the following steps:
1. with total permission water lift amount BZ of high basin moisturizing pumping plant2For coordination variable, enable its practical water lift before being not optimised total It is discrete by step-length d between amount and the practical year total water shortage of high basin, and model is substituted into, formula (1)-(5), to high basin and its small pump Subsystem of standing optimizes, and obtains each stage water deficit quadratic sum F in high basin year2, water supply process X2,i, abandon process water PS2,iAnd its moisturizing pumping plant mentions process water Y2,i。
2. by the benefit library water Y of high basin2iThe water resource for substituting into low library and its moisturizing pumping station system as one of water requirement is excellent Change scheduling model, formula (6)-(10) obtain each stage water deficit quadratic sum F in low library year1, water supply process X1,i, abandon water mistake Journey PS1,iAnd its moisturizing pumping plant mentions process water Y1,i。
3. with total permission water lift amount BZ of high basin moisturizing pumping plant2For coordination variable (independent variable), with each stage in system year Water deficit quadratic sum F (F=F1+F2) it is that dependent variable carries out regression analysis drafting F~BZ2Relation curve determines optimal F*With it is right The BZ answered2 *, the scheduling process of corresponding low library-high basin and its moisturizing pumping plant is searched, it is final to obtain the excellent of " two libraries, two station " system Change scheduling scheme, i.e. X1,i *, X2,i *, PS1,i *, PS2,i *, Y1,i *And Y2,i *。
Compared with prior art, the invention has the following advantages:
First, the hydrographic data of Reservoir-irrigated District actual measurement and the characteristic of actual reservoir and pumping plant has been fully utilized in the present invention Numerical value carries out the water resources optimal operation of " station of two libraries-two " system, overcomes the subjectivity and human intervention of routine dispactching method Stronger disadvantage precisely can easily determine optimal system call scheme, solve Seasonal Water-Deficit, improve water resource and utilize Rate.
Second, this method, without any particular/special requirement, and has comprehensively considered the requirement of regional water rights to data used, is easy to Southern hilly area is promoted the use of in similar "-two station of two libraries " Reservoir-irrigated District.
Detailed description of the invention
Fig. 1 is " station of two libraries-two " water resource system schematic diagram of river and lake moisturizing.
Fig. 2 is " station of two libraries-two " water resource system F~BZ2Regression curve.
Specific embodiment
With reference to the accompanying drawing, to this hair by taking " station of two libraries-two " water resource system of Eastern China province Hills as an example It is bright to be further described.
(1) it is seeked advice from by actual measurement or to irrigated areas administration department, collects the related data needed for calculating, specifically include:
1. schedule periods (water conservancy year) are divided into 20 periods, wherein flood season (September in June -), remaining pressed month by ten days point；
2. collecting reservoir carrys out water data, calculate determining reservoir day part carrys out water LS1,i=[115,45,11,32,2, 28,32,50,157,106,182,126,403,252,191,175,252,199,103,18] ten thousand m3, LS2,i=[33,55,1, 36,10,9,37,56,6,184,28,27,51,55,37,72,48,10,1,9] ten thousand m3(i=1,2 ..., 20)；
3. collecting the water requirement data of irrigated area day part, YS1,i=[23,15,18,17,23,26,30,31,13,458, 52,107,108,108,51,54,54,11,11,11] ten thousand m3, YS2,i=[59,79,17,16,23,27,31,35,46,460, 44,183,176,199,105,78,162,37,13,27] ten thousand m3(i=1,2 ..., 20)；
4. collecting the characteristic data of two reservoirs: reservoir initial storage V1,0=1,047 ten thousand m3, V2,0=629 ten thousand m3, on storage capacity Lower limitWith
5. collecting the characteristic data of two pumping plants: design water lift flow Q1=9936m3/ h, Q2=5400m3/ h, low library are mended The practical year water lift total amount of pump works is 3,600,000 m3It is 11,000,000 m with the practical year water lift total amount of high basin moisturizing pumping plant3；
6. determining that system allows water lift from the year of external river diversion according to the water rights data that local water administrative sector distributes Total amount, i.e., low library moisturizing pumping plant year allow water lift total amount [BZ1The m of]=3,600,0003。
(2) " station of two libraries-two " Large system optimization model decomposition is high basin and its moisturizing pumping station system by large system decomposition With low library and its moisturizing pumping station system two subsystems Optimized model.
1. high basin and its moisturizing pumping station system
Objective function:
Constraint condition:
V2, i=V2,i-1+LS2,i-X2,i+Y2,i-PS2,i-EF2,i (4)
457≤V2,i≤1393 (5)
2. low library and its moisturizing pumping station system
Objective function:
Constraint condition:
V1,i=V1,i-1+LS1,i-X1,i+Y1,i-PS1,i-EF1,i (9)
600≤V1,i≤1757 (10)
(3) subsystem optimizes
With the permission water lift amount BZ of high basin moisturizing pumping plant2For coordination variable, enable it in ten thousand m of feasible zone [1100,1260]3It is interior, 100,000 m of scheduling accuracy being accustomed to local water reservoir management department3It is discrete for step-length, to high basin and its moisturizing pumping plant subsystem into Row optimization, obtains each stage water deficit quadratic sum F in high basin year2, water supply process X2,i, abandon process water PS2,iAnd its moisturizing Pumping plant mentions process water Y2,i；By the water lift amount Y of high basin moisturizing pumping plant2,iLow library and its small pump are substituted into as one of water requirement The Optimized model for system of standing obtains each stage water deficit quadratic sum F in low library year1, water supply process X1i, abandon process water PS1,i And its moisturizing pumping plant mentions process water Y1,i。
With BZ2=1,100 ten thousand m3When for, by BZ2High basin and its moisturizing pumping plant system water priority scheduling of resource model are substituted into, Formula (1)~(5) solve the scheduling scheme for obtaining high basin and its moisturizing pumping station system using Dynamic Programming:
X2,i=[47,67,5,4,10,14,18,22,33,436,44,168,176,196,105,78,161,3 6,12, 26] ten thousand m3
PS2,i=[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] ten thousand m3
Y2,iTen thousand m of=[0,12,16,0,0,0,0,0,0,150,24,150,134,150,78,17,122,33,64,1 50]3
Corresponding F at this time2=2252.
By Y2,iLow library and its moisturizing pumping plant system water priority scheduling of resource model are substituted into as one of water requirement, formula (6)~ (10), the scheduling scheme for obtaining low library and its moisturizing pumping station system is solved using Dynamic Programming:
X1,i=[23,27,34,17,23,26,30,31,13,608,76,257,242,258,129,71,176,44,75, 161] ten thousand m3
PS1,i=[0,0,0,0,0,0,0,0,143,0,0,0,0,0,0,0,0,80,18,0] ten thousand m3
Y1,i=[0,0,0,0,0,0,0,0,0,114,0,47,0,0,0,0,0,0,0,199] ten thousand m3
Corresponding F at this time1=0
Finally, working as BZ2=1,100 ten thousand m3When, F=F1+F2=2252.
BZ2Corresponding objective function F when other discrete values is taken to see Fig. 2.
(4) big system coordination
With stage water deficit quadratic sum F (F=F each in system year1+F2) it is dependent variable, allow to mention with high basin moisturizing pumping plant Water BZ2Regression analysis is carried out for independent variable, draws F~BZ2Relation curve, as shown in Figure 2.Determine optimal F*=250 Hes BZ2 *=1,230 ten thousand m3And its scheduling process of corresponding two reservoirs and pumping plant, the final optimization for obtaining " two libraries, two station " system Scheduling scheme:
X1,i *=[92,52,47,17,30,53,39,32,13,592,75,257,242,258,127,69,177,43, 75,161] ten thousand m3
X2,i *=[59,79,17,16,22,26,30,34,45,459,43,168,175,196,104,77,161,36, 12,26] ten thousand m3
PS1,i *=[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,83,19,0] ten thousand m3
PS2,i *=[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] ten thousand m3
Y1,i *=[0,0,0,0,0,0,0,0,0,101,0,46,0,0,0,0,0,0,0,200] ten thousand m3
Y2,i *=[69,37,29,0,7,27,9,1,0,134,23,150,134,150,76,15,123,32,64,150] Ten thousand m3
The hydrographic data of Reservoir-irrigated District actual measurement and the intrinsic numeric of actual reservoir and pumping plant has been fully utilized in the present invention The water resources optimal operation of " station of the two libraries-two " system of progress, subjectivity and the human intervention for overcoming routine dispactching method are stronger The shortcomings that, it precisely can easily determine optimal system call scheme, solve Seasonal Water-Deficit, improve water utilization rate. This method, without any particular/special requirement, and has comprehensively considered the requirement of regional water rights to data used, is easy in southern hilly area It is promoted the use of in similar "-two station of two libraries " Reservoir-irrigated District.
The above is only a preferred embodiment of the present invention, it should be pointed out that: for the ordinary skill people of the art For member, various improvements and modifications may be made without departing from the principle of the present invention, these improvements and modifications are also answered It is considered as protection scope of the present invention.

Claims (3)

1. a kind of " station of two libraries-two " system water resources optimal operation method of river and lake moisturizing, " station of two libraries-two " of river and lake moisturizing System is made of two reservoirs and two pumping plants, and two reservoirs include low library and high basin, is respectively supplied water to independent water user, and two Seat pumping plant includes low library moisturizing pumping plant and high basin moisturizing pumping plant, and low library moisturizing pumping plant water lift from external river supplements low library, high Moisturizing pumping plant in library supplements high basin from low library water lift；
System optimization dispatching method the following steps are included:
(1) following information is obtained, is specifically included:
1. irrigated area reservoir operation usually divide when number of segment T and corresponding day part length Ni；
2. represent reservoir day part in year carrys out water LS1,i, LS2,i(i=1,2 ..., T)；
3. representing the water requirement YS of irrigated area day part in year1,i, YS2,i(i=1,2 ..., T)；
4. the characteristic data in low library, high basin: reservoir initial storage V1,0, V2,0, storage capacity boundWith
5. the characteristic data of low library moisturizing pumping plant, high basin moisturizing pumping plant: design water lift flow Q1, Q2With total water lift amount of actual year Data；
6. the water rights data of local water administrative sector distribution, referring mainly to system from the year of external river diversion allows water lift total amount BZ1；
(2) " station of two libraries-two " system water resources optimal operation model decomposition is high basin and its moisturizing pumping plant system by system decomposition System and low library and its moisturizing pumping station system two subsystems water resources optimal operation model；
1. high basin and its moisturizing pumping plant
Objective function:
In formula: F2For high basin and its day part water deficit quadratic sum in moisturizing pumping plant subsystem year；X2,iFor high basin and its small pump It stands the water supply of the i-th period of subsystem；YS2,iFor high basin and its water requirement of moisturizing pumping plant the i-th period of subsystem；
Constraint condition:
Y2,i≤Q2×Ni (3)
V2,i=V2,i-1+LS2,i-X2,i+Y2,i-PS2,i-EF2,i (4)
V2,i m≤V2,i≤V2,i M (5)
In formula: LS2,iFor high basin period i two Phase flow amount；BZ2To allow water lift amount in the year of high basin moisturizing pumping plant；Y2,iFor height Water lift amount of the moisturizing pumping plant in library within the i period；Q2For the design discharge of high basin moisturizing pumping plant；NiFor the length of period i；V2,iWith V2,i-1Respectively reservoir storage of the high basin in i and i-1 period Mo；PS2, iFor high basin the i period abandoning water；EF2,iIt is high basin in i The water loss of period；For high basin the i period reservoir storage lower limit；For high basin the i period the reservoir storage upper limit；
2. low library and its moisturizing pumping plant
Objective function:
In formula: F1For low library and its day part water deficit quadratic sum in moisturizing pumping plant subsystem year；X1,iFor low library and its small pump It stands the water supply of the i-th period of subsystem；YS1,iFor low library and its water requirement of moisturizing pumping plant the i-th period of subsystem；Y2,iFor high basin The water that moisturizing pumping plant is extracted from low library within the i period；
Constraint condition:
Y1,i≤Q1×Ni (8)
V1,i=V1,i-1+LS1,i-X1,i+Y1,i-PS1,i-EF1,i (9)
V1,i m≤V1,i≤V1,i M (10)
In formula: BZ1To allow water lift amount, i.e., the water rights of local water administrative sector distribution low library moisturizing pumping plant year；Y1,iFor low library benefit Water lift amount of the pump works within the i period；Q1For the design discharge of low library moisturizing pumping plant；V1,iAnd V1,i-1Respectively low library is in i and i- The reservoir storage of 1 period Mo；PS1,iFor low library the i period abandoning water；EF1,iFor low library the i period water loss；It is low Reservoir storage lower limit of the library in the i period；For low library the i period the reservoir storage upper limit；
(3) subsystem optimizes
With the permission water lift amount BZ of high basin moisturizing pumping plant2For coordination variable, its practical year water lift total amount and high basin before being not optimised is enabled It is discrete by the step-length d for meeting system call required precision between practical year total water shortage, to high basin and its moisturizing pumping plant subsystem It optimizes, obtains each stage water deficit quadratic sum F in high basin year2, water supply process X2,i, abandon process water PS2,iAnd its it mends Pump works mentions process water Y2,i；By the water lift amount Y of high basin moisturizing pumping plant2,iLow library and its moisturizing are substituted into as one of water requirement The Optimized model of pumping station system obtains each stage water deficit quadratic sum F in low library year1, water supply process X1,i, abandon process water PS1,iAnd its moisturizing pumping plant mentions process water Y1,i；
(4) system coordination
Using stage water deficit quadratic sum F each in system year as dependent variable, F=F1+F2, with the permission water lift amount of high basin moisturizing pumping plant BZ2Regression analysis is carried out for independent variable, determines optimal F and BZ2The scheduling of value and its corresponding two reservoirs and pumping plant Journey, the final Optimized Operation scheme for obtaining "-two station of two libraries " system.
2. the method according to claim 1, wherein the step (3) is by " station of two libraries-two " large system decomposition One-dimensional dynamic programming method is respectively adopted and optimizes for " the mono- station in single library-" subsystem, the specific steps are as follows:
(1) stage i=1:
f1(λ1)=min (X1-YS1)2 (11)
In formula, λ1For state variable, preceding 1 water supply period reservoir yield is indicated, it is discrete by step-length d in corresponding feasible zone: λ1 =0, W1,W2,…,∑LSi+BZ；To each discrete λ1, decision variable X1[0, λ in corresponding feasible zone1] discrete by step-length d, It determines respectively and corresponds to each λ1Value, optimal X1And its corresponding period minimum water deficit quadratic sum f1(λ1)；
Then, the 1st stage end pondage V1=V0+LS1-EF1-X1, not yet consider that water is abandoned in pumping plant moisturizing or reservoir, is answered at this time It tests and corrects:
1. if V1<V1 m, then Y1=min (V1 m-V1, Q × N1), PS1=0；
2. if V1>V1 M, then Y1=0, PS1=V1-V1 M；
3. if V1 m≤V1≤V1 M, then Y1=PS1=0；
4. revised pondage is V1'=V1+Y1-PS1；
By step 1.~4., correct and simultaneously determine the 1st stage end pondage, while can get corresponding reservoir and abandoning water PS1Or pumping plant rate of water make-up Y1；
(2) stage i=2,3 ... T-1:
fi(λi)=min [(Xi-YSi)2+fi-1(λi-1)] (12)
In formula, state variable λiFor the preceding i period reservoir for water inventory, equally it is carried out respectively by step-length d discrete: λi= 0,W1,W2,…,∑LSi+BZ；To each discrete λi, decision variable Xi[0, λ in corresponding feasible zonei] discrete by step-length d；
State transition equation:
λi-1=λi-Xi (13)
In formula: i=2,3 ..., T-1
To each discrete λi, by each discrete XiValue substitutes into the (X in formula (12) respectivelyi-YSi)2, by state transition equation formula (13), to each discrete Xi, i-1 stage minimum f before searchingi-1(λi-1) value, it thus can get (Xi-YSi)2+minfi-1 (λi-1), complete all of above discrete XiAfter optimizing, final can get meets min [(Xi-YSi)2+fi-1(λi-1)] require before I period system minimum water deficit quadratic sum fi(λi) value and its corresponding each stage reservoir optimal water supply Xi；
Equally, it also needs to determine the i-th stage end pondage Vi, Vi=Vi-1+LSi-EFi-Xi, it then tests and corrects,
1. if Vi<Vi m, then Yi=min (Vi m-Vi, Q × Ni), PSi=0；
2. if Vi>Vi M, then Yi=0, PSi=Vi-Vi M；
3. if Vi m≤Vi≤Vi M, then Yi=PSi=0；
4. revised pondage is Vi'=Vi+Yi-PSi；
By step 1.~4., correct and simultaneously determine the i-th period end pondage, while can get corresponding reservoir and abandoning water mistake Journey PSiAnd pumping plant rate of water make-up process Yi；
(3) stage T:
fT(λT)=min [(XT-YST)2+fT-1(λT-1)] (14)
State variable λTFor the preceding T period reservoir for water inventory, equally it is carried out respectively discrete: λT=0, W1,W2,…,∑ LSi+BZ；To each discrete λT, decision variable XT[0, λ in corresponding feasible zoneT] discrete by step-length d；
State transition equation:
λT-1=λT-XT (15)
To each discrete λT, by each discrete XTValue substitutes into the (X in formula (14) respectivelyT-YST)2, by state transition equation formula (15), to each discrete XT, T-1 stage minimum f before searchingT-1(λT-1) value, it thus can get (XT-YST)2+minfT-1 (λT-1), complete all of above discrete XTAfter optimizing, final can get meets min [(XT-YST)2+fT-1(λT-1)] require before T period system minimum water deficit quadratic sum fT(λT) value, it finally obtains and meets the λTIt is required that the optimal water supply process X of reservoiri, i =1,2 ..., T, corresponding reservoir abandon process water PSi, i=1,2 ..., T, pumping plant rate of water make-up process Yi, i=1,2 ..., T, with And master mould objective function optimal value F*=minfT(λT)。
3. the method according to claim 1, wherein with total permission of high basin moisturizing pumping plant in the step (4) Water lift amount BZ2For coordination variable, i.e. independent variable, with stage water deficit quadratic sum F, F=F each in system year1+F2, be dependent variable into Row regression analysis, the specific steps are as follows:
1. with total permission water lift amount BZ of high basin moisturizing pumping plant2For coordination variable, its practical water lift total amount and height before being not optimised is enabled It is discrete by step-length d between the practical year total water shortage in library, and model is substituted into, formula (1)-(5), to high basin and its moisturizing pumping plant subsystem System optimizes, and obtains each stage water deficit quadratic sum F in high basin year2, water supply process X2,i, abandon process water PS2,iAnd its Moisturizing pumping plant mentions process water Y2,i；
2. by the benefit library water Y of high basin2,iThe water resource optimization tune of low library and its moisturizing pumping station system is substituted into as one of water requirement Model is spent, formula (6)-(10) obtain each stage water deficit quadratic sum F in low library year1, water supply process X1,i, abandon process water PS1,iAnd its moisturizing pumping plant mentions process water Y1,i；
3. with total permission water lift amount BZ of high basin moisturizing pumping plant2For independent variable, it is with stage water deficit quadratic sum F each in system year Dependent variable carries out regression analysis, draws F~BZ2Relation curve determines optimal F*With corresponding BZ2 *, search corresponding low The scheduling process of library-high basin and its moisturizing pumping plant, the final Optimized Operation scheme for obtaining "-two station of two libraries " system, i.e. X1,i *, X2,i *, PS1,i *, PS2,i *, Y1,i *And Y2,i *。
CN201811216569.6A 2018-10-18 2018-10-18 Water resource optimization scheduling method for 'two-reservoir-two-station' system for replenishing water in rivers and lakes Active CN109002946B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811216569.6A CN109002946B (en) 2018-10-18 2018-10-18 Water resource optimization scheduling method for 'two-reservoir-two-station' system for replenishing water in rivers and lakes

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811216569.6A CN109002946B (en) 2018-10-18 2018-10-18 Water resource optimization scheduling method for 'two-reservoir-two-station' system for replenishing water in rivers and lakes

Publications (2)

Publication Number Publication Date
CN109002946A true CN109002946A (en) 2018-12-14
CN109002946B CN109002946B (en) 2021-06-25

Family

ID=64590210

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811216569.6A Active CN109002946B (en) 2018-10-18 2018-10-18 Water resource optimization scheduling method for 'two-reservoir-two-station' system for replenishing water in rivers and lakes

Country Status (1)

Country Link
CN (1) CN109002946B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110633849A (en) * 2019-09-10 2019-12-31 水利部交通运输部国家能源局南京水利科学研究院 Water resource configuration method based on heuristic large system decomposition coordination idea
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
CN111401750A (en) * 2020-03-18 2020-07-10 四川省水利科学研究院 Method for optimizing and allocating field water resources in reservoir irrigation area

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103981832A (en) * 2014-05-27 2014-08-13 重庆大学 Method and system for urban small river ecological water compensation scheduling
CN107742166A (en) * 2017-10-19 2018-02-27 扬州大学 More storehouse multiple station systems water resource optimal allocation methods in storehouse are directly mended under a kind of fully irrigation conditions
CN107798471A (en) * 2017-10-19 2018-03-13 扬州大学 More storehouse multiple station systems water resource optimal allocation methods of canal are directly mended under a kind of fully irrigation conditions

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103981832A (en) * 2014-05-27 2014-08-13 重庆大学 Method and system for urban small river ecological water compensation scheduling
CN107742166A (en) * 2017-10-19 2018-02-27 扬州大学 More storehouse multiple station systems water resource optimal allocation methods in storehouse are directly mended under a kind of fully irrigation conditions
CN107798471A (en) * 2017-10-19 2018-03-13 扬州大学 More storehouse multiple station systems water resource optimal allocation methods of canal are directly mended under a kind of fully irrigation conditions

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
史振铜等: "南水北调东线江苏段库群水资源优化调度方法研究", 《灌溉排水学报》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110633849A (en) * 2019-09-10 2019-12-31 水利部交通运输部国家能源局南京水利科学研究院 Water resource configuration method based on heuristic large system decomposition coordination idea
CN110633849B (en) * 2019-09-10 2022-05-20 水利部交通运输部国家能源局南京水利科学研究院 Water resource configuration method based on heuristic large system decomposition coordination idea
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
CN111401750A (en) * 2020-03-18 2020-07-10 四川省水利科学研究院 Method for optimizing and allocating field water resources in reservoir irrigation area
CN111401750B (en) * 2020-03-18 2020-11-17 四川省水利科学研究院 Method for optimizing and allocating field water resources in reservoir irrigation area

Also Published As

Publication number Publication date
CN109002946B (en) 2021-06-25

Similar Documents

Publication Publication Date Title
Lecina et al. Irrigation modernization and water conservation in Spain: The case of Riegos del Alto Aragón
CN105096004B (en) A kind of multi-reservoir supplies water transfer system real-time scheduling method
Dudley A single decision‐maker approach to irrigation reservoir and farm management decision making
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
CN109002946A (en) A kind of &#34; station of two libraries-two &#34; system water resources optimal operation method of river and lake moisturizing
CN106295893A (en) Fully directly mend single pumping plant list water reservoir system water resource optimal allocation method in storehouse under irrigation conditions
CN106327065A (en) Water resource optimization configuration method for single pumping station - single reservoir system for direct canal supplement under full irrigation condition
CN113407897A (en) Design method of distributed water circulation model based on multi-source mutual-aid water supply mode
CN102254237A (en) Method for predicting and dynamically regulating regional water resource supply-demand situations
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
Mahan et al. Market mechanisms and the efficient allocation of surface water resources in southern Alberta
CN116562583A (en) Multidimensional water resource supply and demand prediction method and system
CN115796381A (en) Actual runoff forecasting method based on improved Xinanjiang model
Hong et al. Analysis of changes in water allocations and crop production in the Zhanghe Irrigation System and District, 1966–1998
CN116402410B (en) Distributed water quantity and water quality configuration method and system
CN115496622B (en) MIKEBASIN water resource configuration model file rapid generation method
CN109829580A (en) Abundant irrigation conditions lower storage reservoir-benefit library pumping plant-benefit canal pumping station system water resource optimal allocation method
Mojid Emerging challenges and way forward for irrigated agriculture of Bangladesh
Lamaddalena et al. Participatory water management in Italy: case study of the Consortium “Bonifica della Capitanata”
Ramesh et al. Zero-one programming model for daily operation scheduling of irrigation canal
Taha et al. Projection of irrigation demand in Jordan: introduction to GIS-based model for the projection and management of irrigation demand
CN116128303A (en) Network risk assessment method for watershed water and soil resource system
Li et al. Water-saving Irrigation for Rice: Proceedings of an International Workshop Held in Wuhan, China, 23-25 March 2001
Miao et al. Irrigation Water-Saving Technologies to Adapt to Global Changes in the Yellow River Basin, China: A Hetao Case Study
Atarian A Climate, Land, Energy, and Water Nexus Assessment of Ethiopia

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
GR01 Patent grant
GR01 Patent grant