CN104268645A - Industry water supplying quantity determining method with water supplying priority level taken into consideration in water resource allocation - Google Patents
Industry water supplying quantity determining method with water supplying priority level taken into consideration in water resource allocation Download PDFInfo
- Publication number
- CN104268645A CN104268645A CN201410500530.2A CN201410500530A CN104268645A CN 104268645 A CN104268645 A CN 104268645A CN 201410500530 A CN201410500530 A CN 201410500530A CN 104268645 A CN104268645 A CN 104268645A
- Authority
- CN
- China
- Prior art keywords
- water
- industry
- level
- mouthful
- water supply
- 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
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 440
- 238000000034 method Methods 0.000 title claims abstract description 85
- 238000013468 resource allocation Methods 0.000 title abstract 3
- JEGUKCSWCFPDGT-UHFFFAOYSA-N h2o hydrate Chemical compound O.O JEGUKCSWCFPDGT-UHFFFAOYSA-N 0.000 claims description 20
- 238000013461 design Methods 0.000 claims description 12
- 239000008400 supply water Substances 0.000 claims description 12
- 239000011159 matrix material Substances 0.000 claims description 7
- 208000005156 Dehydration Diseases 0.000 claims description 6
- 239000002245 particle Substances 0.000 claims description 6
- 230000001052 transient effect Effects 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 5
- 230000001133 acceleration Effects 0.000 claims description 4
- 230000005484 gravity Effects 0.000 claims description 4
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 claims description 3
- 241001672694 Citrus reticulata Species 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 7
- 238000007726 management method Methods 0.000 abstract description 6
- 238000009826 distribution Methods 0.000 description 6
- 238000012886 linear function Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000007599 discharging Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 230000002123 temporal effect Effects 0.000 description 5
- 238000011161 development Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 241000251468 Actinopterygii Species 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 238000012946 outsourcing Methods 0.000 description 2
- 241000196324 Embryophyta Species 0.000 description 1
- 238000006424 Flood reaction Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000008020 evaporation Effects 0.000 description 1
- 238000001704 evaporation Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 230000001932 seasonal effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000005641 tunneling Effects 0.000 description 1
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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses an industry water supplying quantity determining method with the water supplying priority level taken into consideration in water resource allocation. The industry water supplying quantity is determined according to the industry water supplying priority level in the future water resource allocation of appointed areas, an industry water supplying priority level three-layer structure is established through a wetted perimeter method, a line outer-wrapping method, a water balance method, an SOR method and various survey data. The water supplying quantity calculation of various water supplying industries is scientifically and effectively achieved according to the different water level limitation of water users of the water supplying industries and the different water supplying priority levels. According to the industry water supplying quantity determining method, the water supplying priority level of ecological water level control is introduced, the different water level limitation and the different water supplying priority levels of the water users in life, industries, ship locks, ecology, agriculture and the like, the water supplying quantities are independently determined, the increase of the utilization rate of the water resources is promoted, the water ecology is protected, and water resource management and allocation are more scientific and effective.
Description
Technical field
The present invention relates to a kind of Water Resources Allocation method, particularly relate to a kind of output defining method for different industries, belong to and produce the technical field such as runoff concentration simulation and supply demand, Water Resources Allocation, Hydraulic Projects colony dispatching, water conservancy construction, water operation.
Background technology
In recent years, for the development in pluralism of adjusted drainage system, implement efficiently to manage, Water Resources Allocation to external expansion by water quantity model, has been carried out water quantity and quality and has been combined the researchs such as configuration and Water Resource Adjustment & Distribution.
Along with the intensification that socioeconomic development and people are familiar with water, the domestic theoretical research to Water Resources Allocation develops into water resources development improvement as globality target from the water operation of initial " supply-decided model " and the equilibrium of supply and demand, attempt the reasonable allocation of water resources [Wang Hao realizing basin and region with continuable theory, trip marches. reasonable allocation of water resources research progress and progress [J]. and Journal of Hydraulic Engineering, 2008,39 (10): 1168-1175.].Wang Hao etc. propose " three-time balances " collocation method [Wang Hao based on sustainable utilization of water resource, Qin great Yong, Wang Jianhua, Deng. Yellow-huai river basin reasonable allocation of water resources [M]. Beijing: Science Press, 2003,10:21-159.], in basin and planning for regional water resources, obtain general applying.
In most of Water Resources Allocation research in the past, main employing lowest water level is controlling water level, water user supplies water simultaneously, manually determine water supply precedence according to management practice after determining Water Resources Allocation scheme, cause the river water operation problem that water resource macroscopic allocation disconnects mutually with the actual microcosmic demand of water user thus.
Application number is the application for a patent for invention of 201210241484.X, disclose a kind of water diversion project water resource optimal allocation dynamic model constructing system, wherein: system generally changes module can carry out abstract and generalization formation data file to practical problems, and this data file will be called in channel flow amount computing module and water user's water requirement statistical module; Channel flow acquisition module obtains channel flow data file; Water user's water requirement statistical module counts water user water requirement; Model construction module utilizes optimum theory, the objective function of Confirming model and every constraint condition, founding mathematical models, and mathematical model writes model file after confirming; Model calculation CMOS macro cell model result file; Model result processing module carries out overview display and process to model result.Adopt this invention to generate dynamic model, according to actual rainfall amount situation, Optimized model is revised in real time, thus improve the precision of model prediction.
Domesticly publish document and have " Base of Liaohe Plain underground water eco-grounderwater level and water resource optimal allocation are studied ", propose local ground watering eco-grounderwater level theory and the water quantity regulation computing method of ecology-oriented, but in such scheme, do not relate in Water Resources Allocation the industry output defining method considering water supply precedence.
Summary of the invention
Technical matters to be solved by this invention implements the work of industry water operation, solve the river water operation problem that water resource macroscopic allocation disconnects mutually with the actual microcosmic demand of water user, provide the water users such as one way of life, industry, ship lock, ecology, agricultural in Water Resources Allocation, consider the industry output defining method of water supply precedence.
The present invention is for solving the problems of the technologies described above by the following technical solutions:
Considering an industry output defining method for water supply precedence in Water Resources Allocation, for considering in the following Water Resources Allocation of appointed area that industry water supply priority ranking carries out industry output and determines, comprising the steps:
Step 1, according to waterworks intake head and the minimum water intaking depth of water of industrial self-contained water source intake head design, determine water water level Z at the bottom of the river cross-section of survey region
min, wherein, Z
min=E+H, E are bed level of the river, and H is the minimum water intaking depth of water, and unit is m;
Step 2, fracture morphology according to river course, adopt the river cross-section eco-grounderwater level Z of wetted perimeter method based on river course form and river course lowest navigable stage determination survey region
e, unit is m;
The Basic equation group of the current particle transient motion in the flat trapezoidal river course that step 3, basis are variable with water level and flow sets up survey region network of waterways node water balance system of equations, then adopts the current real time water level Z of SOR solution by iterative method river cross-section;
Step 4, all water user water requirement D of acquisition survey region
kj, water intaking channel length L
i, intake door ability A
kjdata, calculates each the accumulative of water supply industry k on every bar river course i and needs water
accumulative mouth door ability
pass through formula
determine the mouth door output of each water supply industry on every bar river course i; Wherein, j is water user's number, j=1,2 ... N, i are river course number, i=1,2 ... M, N, M are natural number;
Step 5, according to river course bed level of the river E, end water water level Z
min, eco-grounderwater level Z
eset up water supply priority ranking three-decker, define end water water level Z respectively
minabove water supply industry and eco-grounderwater level Z
eabove water supply industry, according to water supply industry residing level, current real time water level Z and Z in three-decker
minand Z
erelation, and the mouth door output of each water supply industry and industry allow the relation of the maximum flow utilized, and determine the output of water supply industry water supply priority ranking and correspondence thereof.
Further, consider the industry output defining method of water supply precedence in a kind of Water Resources Allocation of the present invention, described step 2 comprises the steps:
Step 201, according to river course first and last actual measurement big cross section data, draw its fracture morphology and classify, if this river course is natural river course, then fracture morphology is divided into parabolic type or W type, if this river course is artificial canal, then section configuration is divided into quadrilateral or approximate quadrangle form;
Step 202, to inquire into from the physical dimension of multiple river cross-section and discharge relation measured data, or calculate the relation curve of wetted perimeter and flow from single river channel one group of physical dimension and data on flows;
The position of the catastrophe point in step 203, foundation wetted perimeter and discharge relation curve, determines river channel ecology water level in conjunction with river cross-section typoiogical classification and river course lowest navigable stage.
Further, considering the industry output defining method of water supply precedence in a kind of Water Resources Allocation of the present invention, in described step 203, for the position of the catastrophe point in wetted perimeter and discharge relation curve, is adopt Slope Method or method of maximum curvature to ask for.
Further, in a kind of Water Resources Allocation of the present invention, consider the industry output defining method of water supply precedence, in described step 203, be specially in conjunction with river cross-section typoiogical classification and river course lowest navigable stage determination river channel ecology water level:
(1) if river cross-section form is parabolic type or W type, obtain its catastrophe point by method of maximum curvature or Slope Method, in both the water level corresponding to catastrophe point and river course lowest navigable stage, get maximal value as river channel ecology water level;
(2) if river cross-section form is quadrilateral or approximate quadrangular section, according to the eco-grounderwater level that the hydrobiological minimum existence depth of water in river course or river course wetted perimeter rate require, river course lowest navigable stage determines river.
Further, consider the industry output defining method of water supply precedence in a kind of Water Resources Allocation of the present invention, described in step 3, the Basic equation group of the current particle transient motion in flat trapezoidal river course is:
Wherein, Q is flow, and unit is m
3/ s; X is distance, and unit is m; B is that the water surface is wide, and unit is m; Z is water level, and unit is m; T is the time, and unit is s; q
lfor side enters outflow, unit is m
2/ s, becomes a mandarin as just, and it is negative for going out stream; U is flow velocity, and unit is m/s; A is discharge area, and unit is m
2; R is hydraulic radius, and unit is m;
wherein, sz represents river cross-section wetted perimeter, and unit is m; G is acceleration of gravity; N is channel roughness;
s is both sides slope coefficient sums, Δ Z
*for a upper period water level increment, Δ t is step-length computing time.
Further, consider the industry output defining method of water supply precedence in a kind of Water Resources Allocation of the present invention, the process of iteration of SOR described in step 3 is:
Wherein, x is water level, and k is iterations, i=1,2 ..., n, j=1,2 ..., n, j ≠ i, n is river course, region node total number; ω > 0 is relaxation factor;
represent the water level of i-th node after kth time iteration, a
iifor the matrix of coefficients of network of waterways node water balance system of equations is from diagonal coefficient left to bottom right; b
ifor the constant vector of network of waterways node water balance system of equations; a
ijfor the matrix of coefficients of network of waterways node water balance system of equations is except a
iicoefficient in addition.
Further, consider the industry output defining method of water supply precedence, water supply industry k={L in step 4, I, S, E, A} in a kind of Water Resources Allocation of the present invention, wherein L representative life, I represents industry, and S represents ship lock, and E represents ecology, and A represents agricultural.
Further, the industry output defining method of water supply precedence is considered in a kind of Water Resources Allocation of the present invention, the output of water supply industry water supply priority ranking and correspondence thereof is determined described in step 5, limit under comprising exsiccosis or partly limit the industry supplied water, restriction or part limit output, and the preferential industry supplied water.
Further, in a kind of Water Resources Allocation of the present invention, consider the industry output defining method of water supply precedence, if the mouth door output of life, industry, ship lock, ecology, agricultural is respectively WL
mouthful, WI
mouthful, WS
mouthful, WE
mouthful, WA
mouthful, step 5 comprises the steps:
Step 501, definition life, industry are supplied water successively at water water level of the average end
above, ship lock, ecology, agricultural supply water successively at average eco-grounderwater level
above;
If the current real-time mean water of step 502
lower than
then all industries do not supply water; If current real-time mean water
lower than
then be defined in average eco-grounderwater level
above water supply industry does not supply water; Otherwise water supply industry determines water supply order and output thereof according to the level residing for life, industry, ship lock, ecology, agricultural are in three-decker;
If step 503 is defined in water water level of the average end
above water supply industry exceedes for water inventory the maximum flow Q that industry allows utilization
after, then average eco-grounderwater level is defined in
above water supply industry does not supply water.
Further, the industry output defining method of water supply precedence is considered in a kind of Water Resources Allocation of the present invention, water supply industry described in step 502 determines water supply order and output thereof according to the level residing for life, industry, ship lock, ecology, agricultural are in three-decker, specific as follows:
The junior industry of step a, definition water allows the maximum flow Q utilized
after, the maximum flow Q utilized is allowed by the industry that water higher grade
first,
Wherein,
for the average bottom width in river course,
for the average end, river course is high,
for river course mean water,
for the average eco-grounderwater level in river course, S are side slope, L is channel length, and H is the minimum water intaking depth of water, and unit is all m, and Δ t is time step, and unit is min;
Step b, determine living water amount WL: work as WL
mouthful>Q
firsttime, make WL=Q
first, Q
first=0; Work as WL
mouthful≤ Q
firsttime, make WL=WL
mouthful, Q
first=Q
first-WL
mouthful;
Step c, determine indusqtrial water supply amount WI: work as WI
mouthful>Q
firsttime, make WI=Q
first, Q
first=0; Work as WI
mouthful≤ Q
firsttime, make WI=WI
mouthful, Q
first=Q
first-WI
mouthful;
Steps d, determine ship lock output WS: work as WS
mouthful>Q
aftertime, make WS=Q
after, Q
after=0; Work as WS
mouthful≤ Q
aftertime, WS=WS
mouthful, Q
after=Q
after-WS
mouthful;
Step e, determine ecological output WE: work as WE
mouthful>Q
aftertime, make WE=Q
after, Q
after=0; Work as WE
mouthful≤ Q
aftertime, make WE=WE
mouthful, Q
after=Q
after-WE
mouthful;
Step f, determine agriculture output WA: work as WA
mouthful>Q
aftertime, make WA=Q
after, Q
after=0; Work as WA
mouthful≤ Q
aftertime, make WA=WA
mouthful, Q
after=Q
after-WA
mouthful.
The present invention introduces eco-grounderwater level control concept, and determine the different water level restriction of the water users such as life, industry, ship lock, ecology, agricultural, different water supply priority level, determines output respectively.Based on the water supply precedence introducing eco-grounderwater level control, not only promote the lifting of water resource utilization efficiency, and protect Ecology, make water resources management and configure scientific and effective more.
The present invention adopts above technical scheme compared with prior art, has following technique effect:
(1) consider in the Water Resources Allocation that the present invention designs that the industry output defining method of water supply precedence is by adopting wetted perimeter method, outsourcing collimation method, water balance method, SOR method, utilize multiple survey data, set up industry water supply priority ranking three-decker, according to the different water level restriction of each water supply industry water user, different water supply priority level, scientificlly and effectively achieves each water supply industry water supply calculation separately.
(2) consider in the Water Resources Allocation of the present invention's design that restriction or part restriction are supplied water under the exsiccosis that the industry output defining method of water supply precedence obtains industry, restriction or part limit output, the preferential industry that supplies water and the achievement such as output, output spatial and temporal distributions, for planning engineering layout, the scheduling of lock station, water-saving society establishment, industrial pattern and structural adjustment provide decision-making foundation.
(3) the water supply computation process of the water supply industry water user that the industry output defining method of water supply precedence obtains, restriction water supply industry order and restriction output judged result, branch trade water supply spatial and temporal distributions is considered in the Water Resources Allocation that the present invention designs, precisely can apply in region, river course microcosmic Water Resources Allocation, improve Water Resources Allocation, water operation, scheduling, the efficiency of management evaluation and precision.
Accompanying drawing explanation
Fig. 1 is the industry output defining method process flow diagram considering water supply precedence in the Water Resources Allocation that designs of the present invention.
Fig. 2 is according to river course bed level of the river E, end water water level Z
min, eco-grounderwater level Z
eset up water supply priority ranking three-layered node composition.
Fig. 3 is the node schematic diagram in node water balance equation.
Fig. 4 is the typical river networks node schematic diagram in node water balance system of equations.
Fig. 5 is typical river course schematic diagram.
Fig. 6 be fortune Xi Zha in south on R9 river course (on) the actual measurement big cross section schematic diagram of 1980 years.
Fig. 7 is the water level ~ wetted perimeter graph of relation of south fortune Xi Zha.
Embodiment
Below in conjunction with accompanying drawing, technical scheme of the present invention is described in further detail:
As shown in Figure 1, in a kind of Water Resources Allocation, considering the industry output defining method of water supply precedence, for considering in the following Water Resources Allocation of appointed area that industry water supply priority ranking carries out industry output and determines, comprising the steps:
Step one, according to waterworks intake head and the minimum water intaking depth of water H of industrial self-contained water source intake head design, determine water water level Z at the bottom of the river cross-section of survey region
min(unit: m), wherein, Z
min=E+H, E are bed level of the river, and H is the minimum depth of water.
Step 2, fracture morphology (para-curve, W type, quadrilateral or approximate quadrilateral etc.) according to river course, adopt the river cross-section eco-grounderwater level Z of wetted perimeter method based on river course form and river course lowest navigable stage determination survey region
e(unit: m); It specifically comprises as follows step by step:
Step 201, according to river course first and last actual measurement big cross section data, draw its fracture morphology, and classify, be divided into natural river course fracture morphology parabolic type or " W " type, and artificial canal's section configuration quadrilateral or approximate quadrangle form;
Step 202, to inquire into from the physical dimension of multiple river cross-section and discharge relation measured data, or calculate the relation curve of wetted perimeter and flow from single river channel one group of physical dimension and data on flows;
The position of " catastrophe point " in step 203, foundation wetted perimeter and discharge relation curve, determines river channel ecology water level in conjunction with river cross-section Shape Classification and river course lowest navigable stage.
In step 203, the position for " catastrophe point " in wetted perimeter and discharge relation curve can adopt following mathematical method to ask for:
(1) Slope Method: select the point corresponding to △ Pw/ △ Q=1 to be " catastrophe point ", wherein Pw is wetted perimeter, refers to the circumferential length that on flow section, fluid contacts with solid wall surface, unit m, △ Pw is the variable quantity of wetted perimeter, and △ Q is the fluctuations in discharge amount corresponding with △ Pw.
(2) method of maximum curvature: under many circumstances, water level and wetted perimeter relation present relation or the repeatedly linear relationship of power function, and curve significantly " catastrophe point ", can not use method of maximum curvature to determine " catastrophe point " in this case.Curvature is the function of tangent angle and arc length.Draw curvature curve figure and ask for maximal value.Point corresponding to maximal value is " catastrophe point ".
In step 203, adopt the wetted perimeter method based on river course form and river course lowest navigable stage determination river channel ecology water level, be specially:
(1) if section configuration is parabolic type or " W " type, then there is stable wetted perimeter discharge relation curve in the river course of this type, obtain its " catastrophe point " by method of maximum curvature or Slope Method, in both the water level corresponding to " catastrophe point " and river course lowest navigable stage, get maximal value as river channel ecology water level;
(2) if river cross-section shape is quadrilateral or approximate quadrangular section, the discharge of river increases from 0, width and wetted perimeter increase fast, catastrophe point is produced when flow by chance floods riverbed, now, the water level that catastrophe point is corresponding does not also meet the requirement of eco-grounderwater level, in such cases, need to require and river course lowest navigable stage determines the eco-grounderwater level in river according to the hydrobiological minimum existence depth of water in river course or river course wetted perimeter rate.
Wherein, for the river course that section configuration is quadrilateral or approximate quadrangle form, the determination of its minimum existence depth of water, need to consider whether the depth of water that wetted perimeter method is tried to achieve can maintain the existence of fish from biological angle, with reference to external relevant achievement in research, select the respective standard of dividing work out season according to the biology needs of fish and the seasonal variety in river, its minimum existence depth of water gets 0.3m (i.e. 1ft), and the corresponding hydraulic tunneling feature reaching the critical conditions that water demand for natural service requires is that wetted perimeter rate accounts for 60% of main stem wetted perimeter.
Adopt wetted perimeter method, obtain the water level corresponding to " catastrophe point " that curvature is maximum, and judge whether its water level meets standard or the 60% wetted perimeter rate standard-required of minimum existence depth of water 0.3m.Be expressed as follows:
Z
e=max{Z
1,Z
2,Z
3}
In formula, Z
eriver channel ecology water level, Z
1for the water level calculated according to method of maximum curvature or the Slope Method of wetted perimeter method, Z
2for the river water level according to minimum existence depth of water standard or 60% wetted perimeter rate criterion calculation, Z
3for river course lowest navigable stage.
The Basic equation group of the current particle transient motion in the flat trapezoidal river course that step 3, basis are variable with water level and flow sets up survey region network of waterways node water balance system of equations, the wherein Basic equation group of the current particle transient motion in flat trapezoidal river course:
Wherein, Q is flow, and unit is m
3/ s; X is distance, and unit is m; B is that the water surface is wide, and unit is m; Z is water level, and unit is m; T is the time, and unit is s; q
lfor side enters outflow, unit is m
2/ s, becomes a mandarin as just, and it is negative for going out stream; U is flow velocity, and unit is m/s; A is discharge area, and unit is m
2; R is hydraulic radius, and unit is m;
wherein, sz represents river cross-section wetted perimeter, and unit is m; G is acceleration of gravity; N is channel roughness;
s is both sides slope coefficient sums, Δ Z
*for a upper period water level increment, Δ t is step-length computing time;
Adopt the current real time water level Z of SOR solution by iterative method river cross-section:
Wherein, x is water level, and k is iterations, i=1,2 ..., n, j=1,2 ..., n, j ≠ i, n is river course, region node total number; ω > 0 is relaxation factor;
represent the water level of i-th node after kth time iteration, a
iifor the matrix of coefficients of network of waterways node water balance system of equations is from diagonal coefficient left to bottom right; b
ifor the constant vector of network of waterways node water balance system of equations; a
ijfor the matrix of coefficients of network of waterways node water balance system of equations is except a
iicoefficient in addition.
Wherein, set up survey region network of waterways node water balance system of equations to comprise the steps:
Step 301. sets up node water balance equation (for typical node N1 in Fig. 3):
In Fig. 3, direction shown in arrow is the flow direction of definition.Be positive sign from first node-flow to end-node, otherwise be then negative sign.Q
1, Q
2, Q
3for river course 1,2,3 flow, Q
gfor weir flow amount, Q
4for directly joining the flow of node N1, can be effective rainfall (evaporation of rainfall deduction), turning over the water yield with water or pumping plant.
The flow flowed into or flow out N1 node is expressed as the linear function of this river first and last node water level:
Z
n1, Z
n2, Z
n3, Z
n4, Z
n5be respectively the water level of node N1, N2, N3, N4, N5; a
ij, b
ijfor matrix of coefficients, tried to achieve by (2), (3) of step 301.
Wherein, under free discharge pattern:
Q
g=C
1(Z
N1-H
0),
Under submerge discharging flow pattern:
Q
g=C
2(Z
N1-Z
N5),
Wherein, H
ufor water level on lock at the beginning of the period; H
dfor water level under lock at the beginning of the period; H
0for pocket floor elevation; B is the total span width in lock hole; ε is side direction contraction coefficient (getting 0.9); μ is free discharge coefficient (getting 0.35); φ is submerge discharging flow coefficient (getting 0.90); K is for opening a sluice gate coefficient; G is acceleration of gravity.
Node N1 water balance equation is:
A is node area, and unit is m
2, Δ t is step-length computing time,
for water level at the beginning of the period, it is the given value before each iteration.
In free discharge situation, formula (33) is substituted into (34) to obtain:
K
N1.1Z
N1+K
N1.2Z
N2+K
N1.3Z
N3+K
N1.4Z
N4=K
N1.R (35)
In formula,
k
n1.2=b
13, K
n1.3=b
23, K
n1.4=-a
33,
In submerge discharging flow situation, formula (33) is substituted into (34) to obtain:
K
N1.1Z
N1+K
N1.2Z
N2+K
N1.3Z
N3+K
N1.4Z
N4+K
N1.5Z
N5=K
N1.R (36)
In formula,
k
n1.2=b
13, K
n1.3=b
23, K
n1.4=-a
33, K
n1.5=C
2
K in formula
n1,1, K
n1,2, K
n1,3, K
n1,4, K
n1,5, K
n1, Rfor the coefficient in node water balance equation, all by coefficient that step (1) below, (2), (3) are tried to achieve.
In step 301, every bar river course first and last node 6 coefficient b are obtained by forward and reverse recursion
j1, b
j2, b
j3, a
j1, a
j2, a
j3, j representative calculates river course; For river course typical in Fig. 5, method is as follows:
(1). according to given value at the beginning of the period and selected time step-length and distance step size computation obtain six coefficient values below each section in every bar river course.
In formula, i is river cross-section numbering; J is at the beginning of calculation interval, and all leftover bits and piecess are designated as
person all represents and gets the average of i and i+1 section part functional value.
(2). node N1 is end-node (i.e. flow direction N1)
From the first section, manage this section flow Q
iexpress cost section water level Z
iwith first section water level Z
1linear function: Q
i=θ
i+ η
iz
i+ ν
iz
1
Y
1=C
i-1-η
i-1
Y
2=E
i-1η
i-1-F
i-1
W=Y
2-Y
1G
i-1
Recursion initial value is calculated by following formula:
W=E
1+G
1
From the first section, recursion obtains θ to last section backward successively
n, η
n, ν
n, then by θ
n, η
n, ν
ncorresponding assignment is to b
j1, b
j2, b
j3.
(3). node headed by node N1 (i.e. current wander about as a refugee N1)
From section second from the bottom, manage this section flow Q
iexpress cost section water level Z
iwith last section water level Z
nlinear function: Q
i=α
i+ β
iz
i+ ζ
iz
n
Y
1=C
i+β
i+1
Y
2=G
iβ
i+1+F
i
W=Y
2+Y
1E
i
Recursion initial value is calculated by following formula:
W=G
n-1+E
n-1
From section second from the bottom, recursion obtains α to first section forward successively
1, β
1, ζ
1, then by α
1, β
1, ζ
1corresponding assignment is to a
j1, a
j2, a
j3.
Step 302. sets up network of waterways node water balance system of equations (for typical river networks node in Fig. 4):
Be a simple network of waterways in Fig. 4, have 9 rivers, 11 nodes, 2 lock compositions.
Node 4,5,6 is water level frontier point, and water level is given value, can list following equation to each node:
Z in formula
1, Z
2z
11be 11 node water levels, Q
g1, Q
g2be respectively the lock flow of lock 1 and lock 2, the K of subscripting is the K in formula (35), (36).
Step 4, all water user water requirement D of acquisition survey region
kj, water intaking channel length L
i, intake door ability A
kjdata, calculates each the accumulative of water supply industry k on every bar river course i and needs water
accumulative mouth door ability
by formula W k
mouthful=∑
imin (D
ki, A
ki), determine that the life on every bar river course i, industry, ship lock, ecology, agriculture mouth door output are respectively WL
mouthful, WI
mouthful, WS
mouthful, WE
mouthful, WA
mouthful; Wherein, k classifies according to the water supply industry classification of life, industry, ship lock, ecology, agricultural, k={L, I, S, E, A}, wherein L representative life, and I represents industry, and S represents ship lock, and E represents ecology, and A represents agricultural; J is water user's number, j=1,2 ... N, i are river course number, i=1,2 ... M.
Step 5, according to river course bed level of the river E, end water water level Z
min, eco-grounderwater level Z
eset up water supply priority ranking three-decker, see Fig. 2, define end water water level Z respectively
minabove water supply industry and eco-grounderwater level Z
eabove water supply industry, according to water supply industry residing level, current real time water level Z and Z in three-decker
minand Z
erelation, determine the output (comprise restriction or part restriction are supplied water under exsiccosis industry, restriction or part and limit output, the preferential industry etc. supplied water) of water supply industry water supply priority ranking and correspondence thereof.Specific as follows:
Step I, definition life, industry are supplied water successively at water water level of the average end
above, ship lock, ecology, agricultural supply water successively at average eco-grounderwater level
above;
If the current real-time mean water of Step II
lower than
then all industries do not supply water; If current real-time mean water
lower than
then be defined in average eco-grounderwater level
above water supply industry does not supply water; Otherwise water supply industry determines water supply order and output thereof according to the level residing for life, industry, ship lock, ecology, agricultural are in three-decker;
If Step II I is defined in water water level of the average end
above water supply industry exceedes maximum flow Q for water inventory
after, then average eco-grounderwater level is defined in
above water supply industry does not supply water.
Described water supply industry determines water supply order and output thereof according to the level residing for life, industry, ship lock, ecology, agricultural are in three-decker, specific as follows:
The junior industry of step a, definition water allows the maximum flow Q utilized
after, the maximum flow Q utilized is allowed by the industry that water higher grade
first,
Wherein,
for the average bottom width in river course,
for the average end, river course is high,
for river course mean water,
for the average eco-grounderwater level in river course, S are side slope, L is channel length, and H is the minimum water intaking depth of water, and unit is all m, and Δ t is time step, and unit is min;
Step b, determine living water amount WL: work as WL
mouthful>Q
firsttime, make WL=Q
first, Q
first=0; Work as WL
mouthful≤ Q
firsttime, make WL=WL
mouthful, Q
first=Q
first-WL
mouthful;
Step c, determine indusqtrial water supply amount WI: work as WI
mouthful>Q
firsttime, make WI=Q
first, Q
first=0; Work as WI
mouthful≤ Q
firsttime, make WI=WI
mouthful, Q
first=Q
first-WI
mouthful;
Steps d, determine ship lock output WS: work as WS
mouthful>Q
aftertime, make WS=Q
after, Q
after=0; Work as WS
mouthful≤ Q
aftertime, WS=WS
mouthful, Q
after=Q
after-WS
mouthful;
Step e, determine ecological output WE: work as WE
mouthful>Q
aftertime, make WE=Q
after, Q
after=0; Work as WE
mouthful≤ Q
aftertime, make WE=WE
mouthful, Q
after=Q
after-WE
mouthful;
Step f, determine agriculture output WA: work as WA
mouthful>Q
aftertime, make WA=Q
after, Q
after=0; Work as WA
mouthful≤ Q
aftertime, make WA=WA
mouthful, Q
after=Q
after-WA
mouthful.
Consider in the Water Resources Allocation of the present invention's design that restriction or part restriction are supplied water under the exsiccosis that the industry output defining method of water supply precedence obtains industry, restriction or part limit output, the preferential industry that supplies water and the achievement such as output, output spatial and temporal distributions, for planning engineering layout, the scheduling of lock station, water-saving society establishment, industrial pattern and structural adjustment provide decision-making foundation.
The water supply computation process of the water supply industry water user that the industry output defining method of water supply precedence obtains, restriction water supply industry order and restriction output judged result, branch trade water supply spatial and temporal distributions is considered in the Water Resources Allocation of the present invention's design, precisely can apply in region, river course microcosmic Water Resources Allocation, improve Water Resources Allocation, water operation, scheduling, the efficiency of management evaluation and precision.
Embodiment:
Consider that the industry output defining method of water supply precedence is in the application process of reality, comprises the steps: in the Water Resources Allocation of the present invention's design
(1) investigation obtains waterworks intake head and industrial self-contained water source intake head design minimum water intaking depth of water H ∈ [0.5,1] of survey region, determines water water level Z at the bottom of the first section in the R9 river course of survey region
min is first=E
first+ H=2.61m, water water level Z at the bottom of last section
min end=E
end+ H=3.23m, other 202 river courses in like manner.
(2) according to fortune Xi Zha in south on R9 river course (on) the actual measurement big cross section data of 1980 years, draw its fracture morphology, as shown in Figure 6.
(3) water level ~ wetted perimeter relation curve of south fortune Xi Zha is drawn, as shown in Figure 7.
(4) according to water level ~ wetted perimeter relation curve, adopt method of maximum curvature to obtain " catastrophe point ", obtaining corresponding water level is 4.17m, the river water level 3.27m that 60% wetted perimeter rate standard is corresponding, R9 river course being transported lowest navigable stage on Dong Chuanzhazha, under lock is 0.8 meter, according to formula Z
e=max{Z
1, Z
2, Z
3the first node eco-grounderwater level in the R9 river course that calculates is Z
e=4.17m, in like manner can obtain R9 river course end-node eco-grounderwater level and other 202 river course first and last node eco-grounderwater level.
(5) step-length 15min and distance step-length 3km when selecting, according to first, present situation nineteen eighty-three type given value initial water level, initial flow at the beginning of 15 minutes, calculate R9 river course (N1=" N213 " in Fig. 5, N2=" N303 ", represents first and last node respectively) 6 coefficient values (in Fig. 5, n=11 represents section number) below each section.
In formula: i is river cross-section numbering; J is for calculation interval is at the beginning of 15 minutes, and all leftover bits and piecess are designated as
person all represents and gets the average of i and i+1 section part functional value.
(6) node N1 is end-node
From the first section 1-2, this section flowmeter is reached the linear function of this section water level and first section water level: Q
i=θ
i+ η
iz
i+ ν
iz
1
Y
1=C
i-1-η
i-1
Y
2=E
i-1η
i-1-F
i-1
W=Y
2-Y
1G
i-1
Recursion initial value is calculated by following formula:
W=E
1+G
1
From 1-2 section, recursion obtains θ to last section (n-1)-n backward successively
n, η
n, ν
n, then by θ
n, η
n, ν
ncorresponding assignment is to b
j1, b
j2, b
j3.
(7) node headed by node N1
From section second from the bottom (n-2)-(n-1), this section flowmeter is reached the linear function of this section water level and last section water level: Q
i=α
i+ β
iz
i+ ζ
iz
n
Y
1=C
i+β
i+1
Y
2=G
iβ
i+1+F
i
W=Y
2+Y
1E
i
Recursion initial value is calculated by following formula:
W=G
n-1+E
n-1
From section second from the bottom, recursion obtains α to first section 1-2 forward successively
1, β
1, ζ
1, then by α
1, β
1, ζ
1corresponding assignment is to a
j1, a
j2, a
j3.
(8) node water balance equation (Fig. 3 interior joint N1=" N213 " is the first node of river course R9) is set up
In Fig. 3, direction shown in arrow is the flow direction of definition.Be positive sign from first node-flow to end-node, otherwise be then negative sign.Q
1, Q
2, Q
3for river course 1,2,3 flow, Q
gfor weir flow amount, Q
4for directly joining node N
1flow, be given value.
The flow flowed into or flow out N1 node is expressed as the linear function of this river first and last node water level.
Wherein, under free discharge pattern:
Q
g=C
1(Z
N1-H
0),
Under submerge discharging flow pattern:
Q
g=C
2(Z
N1-Z
NS),
Node N1 water balance equation is:
In free discharge situation, formula (33) is substituted into (34) to obtain
K
N1.1Z
N1+K
N1.2Z
N2+K
N1.3Z
N3+K
N1.4Z
N4=K
N1.R (35)
In formula
k
n1.2=b
13, K
n1.3=b
23, K
n1.4=-a
33,
In submerge discharging flow situation, formula (33) is substituted into (34) to obtain
K
N1.1Z
N1+K
N1.2Z
N2+K
N1.3Z
N3+K
N1.4Z
N4+K
N1.5Z
N5=K
N1.R (36)
In formula
k
n1.2=b
13, K
n1.3=b
23, K
n1.4=-a
33, K
n1.5=C
2
K in formula
n1,1, K
n1,2, K
n1,3, K
n1,4, K
n1,5, K
n1, Rall try to achieve by step (7).
(9) network of waterways node water balance system of equations (be the part typical river networks node of survey region in Fig. 4, have 9 rivers, 11 nodes, 2 lock compositions) is set up
Node 4,5,6 is water level frontier point, and water level is given value, can list following equation to each node:
Z in formula
1, Z
2z
11be 11 node water levels, Q
g1, Q
g2be respectively the lock flow of lock 1 and lock 2.
(10) parameter value in SOR iterative formula is set:
(11) water D is needed according to all water users of survey region
kj, water intaking channel length L
i, intake door ability A
kjdeng field data, calculate each the accumulative of water supply industry k on R9 river course and other 202 river courses of survey region and need water
accumulative mouth door ability
by formula W k
mouthful=∑
imin (D
ki, A
ki), determine that the life of first the 15 minutes period of present situation nineteen eighty-three type on R9 river course and other 202 river courses of survey region, industry, ship lock, ecology, agriculture mouth door output are respectively WL
mouthful=156m
3, WI
mouthful=0m
3, WS
mouthful=591m
3, WE
mouthful=457m
3, WA
mouthful=104m
3.
(12) according to R9 river course bed level of the river E=2.42m, end water water level Z
min=2.92m, eco-grounderwater level Z
e=4.17m sets up water supply priority ranking three-decker, sees Fig. 2, definition end water water level Z
minabove water supply industry is life and industry, definition eco-grounderwater level Z
eabove water supply industry is ship lock, ecology, agricultural, according to water supply industry residing level, present situation nineteen eighty-three type first 15 minutes last water level Z=4.625m and Z in three-decker
minand Z
erelation, determine water supply industry water supply priority ranking for life, industry, ship lock, ecology, agricultural, correspondence output be respectively 156m
3, 0m
3, 591m
3, 457m
3, 104m
3, first, R9 river information system nineteen eighty-three type does not limit water supply industry and restriction output for 15 minutes.In like manner can try to achieve water supply industry water supply priority ranking and output, restriction water supply industry and the restriction output in other 202 river courses of survey region.
(13) using the node water level in present situation nineteen eighty-three type first 15 minutes last survey region 203 river course, section flow as second initial value at the beginning of 15 minutes, resume at step (5) carries out second and calculates for 15 minutes, the like, can in the hope of the water supply industry water supply priority ranking of present situation nineteen eighty-three type survey region 203 river course whole year and output, restriction water supply industry and restriction output.The output of life as annual in R9 river course, industry, ship lock, ecology, agricultural correspondence is respectively 0.05 hundred million m
3, 000,000,000 m
3, 0.21 hundred million m
3, 0.16 hundred million m
3, 0.62 hundred million m
3, unrestricted water supply industry and restriction output.
The industry output defining method of water supply precedence is considered in the Water Resources Allocation of Patent design of the present invention, by adopting wetted perimeter method, outsourcing collimation method, water balance method, SOR method, utilize multiple survey data, set up industry water supply priority ranking three-decker, according to the different water level restriction of each water supply industry water user, different water supply priority level, scientificlly and effectively achieve each water supply industry water supply calculation separately, limit under comprising exsiccosis or partly limit the industry supplied water, restriction or part limit output, the industry of preferential water supply and output, output spatial and temporal distributions etc., in region, can precisely apply in river course microcosmic Water Resources Allocation, improve Water Resources Allocation, water operation, scheduling, the efficiency of management evaluation and precision.
By reference to the accompanying drawings embodiments of the present invention are explained in detail above, but the present invention is not limited to above-mentioned embodiment, in the ken that those of ordinary skill in the art possess, can also makes a variety of changes under the prerequisite not departing from present inventive concept.
Claims (10)
1. consider an industry output defining method for water supply precedence in Water Resources Allocation, for considering in the following Water Resources Allocation of appointed area that industry water supply priority ranking carries out industry output and determines, it is characterized in that, comprise the steps:
Step 1, according to waterworks intake head and the minimum water intaking depth of water of industrial self-contained water source intake head design, determine water water level Z at the bottom of the river cross-section of survey region
min, wherein, Z
min=E+H, E are bed level of the river, and H is the minimum water intaking depth of water, and unit is m;
Step 2, fracture morphology according to river course, adopt the river cross-section eco-grounderwater level Z of wetted perimeter method based on river course form and river course lowest navigable stage determination survey region
e, unit is m;
The Basic equation group of the current particle transient motion in the flat trapezoidal river course that step 3, basis are variable with water level and flow sets up survey region network of waterways node water balance system of equations, then adopts the current real time water level Z of SOR solution by iterative method river cross-section;
Step 4, all water user water requirement D of acquisition survey region
kj, water intaking channel length L
i, intake door ability A
kjdata, calculates each the accumulative of water supply industry k on every bar river course i and needs water
accumulative mouth door ability
by formula W k
mouthful=∑
imin (D
ki, A
ki), determine the mouth door output of each water supply industry on every bar river course i; Wherein, j is water user's number, j=1,2 ... N, i are river course number, i=1,2 ... M, N, M are natural number;
Step 5, according to river course bed level of the river E, end water water level Z
min, eco-grounderwater level Z
eset up water supply priority ranking three-decker, define end water water level Z respectively
minabove water supply industry and eco-grounderwater level Z
eabove water supply industry, according to water supply industry residing level, current real time water level Z and Z in three-decker
minand Z
erelation, and the mouth door output of each water supply industry and industry allow the relation of the maximum flow utilized, and determine the output of water supply industry water supply priority ranking and correspondence thereof.
2. consider the industry output defining method of water supply precedence in a kind of Water Resources Allocation according to claim 1, it is characterized in that, described step 2 comprises the steps:
Step 201, according to river course first and last actual measurement big cross section data, draw its fracture morphology and classify, if this river course is natural river course, then fracture morphology is divided into parabolic type or W type, if this river course is artificial canal, then section configuration is divided into quadrilateral or approximate quadrangle form;
Step 202, to inquire into from the physical dimension of multiple river cross-section and discharge relation measured data, or calculate the relation curve of wetted perimeter and flow from single river channel one group of physical dimension and data on flows;
The position of the catastrophe point in step 203, foundation wetted perimeter and discharge relation curve, determines river channel ecology water level in conjunction with river cross-section typoiogical classification and river course lowest navigable stage.
3. in a kind of Water Resources Allocation according to claim 2, consider the industry output defining method of water supply precedence, it is characterized in that, in described step 203, for the position of the catastrophe point in wetted perimeter and discharge relation curve, be adopt Slope Method or method of maximum curvature to ask for.
4. consider the industry output defining method of water supply precedence in a kind of Water Resources Allocation according to claim 2, it is characterized in that, in described step 203, be specially in conjunction with river cross-section typoiogical classification and river course lowest navigable stage determination river channel ecology water level:
(1) if river cross-section form is parabolic type or W type, obtain its catastrophe point by method of maximum curvature or Slope Method, in both the water level corresponding to catastrophe point and river course lowest navigable stage, get maximal value as river channel ecology water level;
(2) if river cross-section form is quadrilateral or approximate quadrangular section, according to the eco-grounderwater level that the hydrobiological minimum existence depth of water in river course or river course wetted perimeter rate require, river course lowest navigable stage determines river.
5. consider the industry output defining method of water supply precedence in a kind of Water Resources Allocation according to claim 1, it is characterized in that, described in step 3, the Basic equation group of the current particle transient motion in flat trapezoidal river course is:
Wherein, Q is flow, and unit is m
3/ s; X is distance, and unit is m; B is that the water surface is wide, and unit is m; Z is water level, and unit is m; T is the time, and unit is s; q
lfor side enters outflow, unit is m
2/ s, becomes a mandarin as just, and it is negative for going out stream; U is flow velocity, and unit is m/s; A is discharge area, and unit is m
2; R is hydraulic radius, and unit is m;
wherein, sz represents river cross-section wetted perimeter, and unit is m; G is acceleration of gravity; N is channel roughness;
s is both sides slope coefficient sums, Δ Z
*for a upper period water level increment, Δ t is step-length computing time.
6. consider the industry output defining method of water supply precedence in a kind of Water Resources Allocation according to claim 1, it is characterized in that, the process of iteration of SOR described in step 3 is:
Wherein, x is water level, and k is iterations, i=1,2 ..., n, j=1,2 ..., n, j ≠ i, n is river course, region node total number; ω > 0 is relaxation factor;
represent the water level of i-th node after kth time iteration, a
iifor the matrix of coefficients of network of waterways node water balance system of equations is from diagonal coefficient left to bottom right; b
ifor the constant vector of network of waterways node water balance system of equations; a
ijfor the matrix of coefficients of network of waterways node water balance system of equations is except a
iicoefficient in addition.
7. consider the industry output defining method of water supply precedence in a kind of Water Resources Allocation according to claim 1, it is characterized in that, water supply industry k={L, I in step 4, S, E, A}, wherein L representative life, I represents industry, and S represents ship lock, and E represents ecology, and A represents agricultural.
8. in a kind of Water Resources Allocation according to claim 1, consider the industry output defining method of water supply precedence, it is characterized in that, the output of water supply industry water supply priority ranking and correspondence thereof is determined described in step 5, limit under comprising exsiccosis or partly limit the industry supplied water, restriction or part limit output, and the preferential industry supplied water.
9. consider the industry output defining method of water supply precedence in a kind of Water Resources Allocation according to claim 7, it is characterized in that, if the mouth door output of life, industry, ship lock, ecology, agricultural is respectively WL
mouthful, WI
mouthful, WS
mouthful, WE
mouthful, WA
mouthful, step 5 comprises the steps:
Step 501, definition life, industry are supplied water successively at water water level of the average end
above, ship lock, ecology, agricultural supply water successively at average eco-grounderwater level
above;
If the current real-time mean water of step 502
lower than
then all industries do not supply water; If current real-time mean water
lower than
then be defined in average eco-grounderwater level
above water supply industry does not supply water; Otherwise water supply industry determines water supply order and output thereof according to the level residing for life, industry, ship lock, ecology, agricultural are in three-decker;
If step 503 is defined in water water level of the average end
above water supply industry exceedes for water inventory the maximum flow Q that industry allows utilization
after, then average eco-grounderwater level is defined in
above water supply industry does not supply water.
10. in a kind of Water Resources Allocation according to claim 9, consider the industry output defining method of water supply precedence, it is characterized in that, water supply industry described in step 502 determines water supply order and output thereof according to the level residing for life, industry, ship lock, ecology, agricultural are in three-decker, specific as follows:
The junior industry of step a, definition water allows the maximum flow Q utilized
after, the maximum flow Q utilized is allowed by the industry that water higher grade
first,
Wherein,
for the average bottom width in river course,
for the average end, river course is high,
for river course mean water,
for the average eco-grounderwater level in river course, S are side slope, L is channel length, and H is the minimum water intaking depth of water, and unit is all m, and Δ t is time step, and unit is min;
Step b, determine living water amount WL: work as WL
mouthful>Q
firsttime, make WL=Q
first, Q
first=0; Work as WL
mouthful≤ Q
firsttime, make WL=WL
mouthful, Q
first=Q
first-WL
mouthful;
Step c, determine indusqtrial water supply amount WI: work as WI
mouthful>Q
firsttime, make WI=Q
first, Q
first=0; Work as WI
mouthful≤ Q
firsttime, make WI=WI
mouthful, Q
first=Q
first-WI
mouthful;
Steps d, determine ship lock output WS: work as WS
mouthful>Q
aftertime, make WS=Q
after, Q
after=0; Work as WS
mouthful≤ Q
aftertime, WS=WS
mouthful, Q
after=Q
after-WS
mouthful;
Step e, determine ecological output WE: work as WE
mouthful>Q
aftertime, make WE=Q
after, Q
after=0; Work as WE
mouthful≤ Q
aftertime, make WE=WE
mouthful, Q
after=Q
after-WE
mouthful;
Step f, determine agriculture output WA: work as WA
mouthful>Q
aftertime, make WA=Q
after, Q
after=0; Work as WA
mouthful≤ Q
aftertime, make WA=WA
mouthful, Q
after=Q
after-WA
mouthful.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410500530.2A CN104268645B (en) | 2014-09-25 | 2014-09-25 | Industry water supplying quantity determining method with water supplying priority level taken into consideration in water resource allocation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410500530.2A CN104268645B (en) | 2014-09-25 | 2014-09-25 | Industry water supplying quantity determining method with water supplying priority level taken into consideration in water resource allocation |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104268645A true CN104268645A (en) | 2015-01-07 |
CN104268645B CN104268645B (en) | 2015-06-10 |
Family
ID=52160165
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410500530.2A Active CN104268645B (en) | 2014-09-25 | 2014-09-25 | Industry water supplying quantity determining method with water supplying priority level taken into consideration in water resource allocation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104268645B (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107944603A (en) * | 2017-11-09 | 2018-04-20 | 中山大学 | Water resource optimal allocation newsboy's method based on water total amount control |
CN107958095A (en) * | 2016-10-18 | 2018-04-24 | 济南大学 | The parabola shaped open channel of 5/2nds powers of one kind and its optimal hydraulic cross-section |
CN108305006A (en) * | 2018-02-02 | 2018-07-20 | 浙江省水利河口研究院 | Reservoir usufruct of water resource amount primary distribution method based on the constraint of water supply fraction |
CN109767143A (en) * | 2019-03-14 | 2019-05-17 | 南京大学 | Water source method for tracing under multi-source water supply conditions |
CN109948235A (en) * | 2019-03-14 | 2019-06-28 | 南京大学 | Water resource scheduling and precision configuration method |
CN112862225A (en) * | 2019-11-27 | 2021-05-28 | 重庆师范大学 | Lake storage regulation excavation-based water shortage area cross-space-time scheduling method and system |
WO2021129039A1 (en) * | 2019-12-25 | 2021-07-01 | 华中科技大学 | Method and system for calculating cross-reservoir drainage basin river runoff |
CN113326632A (en) * | 2021-06-19 | 2021-08-31 | 南昌工程学院 | Optimization correction method for backward-thrust reservoir warehousing flow process |
CN113449890A (en) * | 2021-07-16 | 2021-09-28 | 中国水利水电科学研究院 | Efficient utilization and optimized allocation method for coal mine water area |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110852518A (en) * | 2019-11-15 | 2020-02-28 | 中国水利水电科学研究院 | Ecological water supplementing method for river estuary area of sediment-laden river |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101645114A (en) * | 2009-09-03 | 2010-02-10 | 北京师范大学 | Technology for analyzing ecological water demand of rivers based on district contribution method |
CN103034772A (en) * | 2011-10-10 | 2013-04-10 | 北京师范大学 | Method for solving contradiction between ecological water and agricultural water |
CN103544533A (en) * | 2012-07-12 | 2014-01-29 | 上海宝信软件股份有限公司 | Diversion project water resource optimal allocation dynamic model establishing system |
-
2014
- 2014-09-25 CN CN201410500530.2A patent/CN104268645B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101645114A (en) * | 2009-09-03 | 2010-02-10 | 北京师范大学 | Technology for analyzing ecological water demand of rivers based on district contribution method |
CN103034772A (en) * | 2011-10-10 | 2013-04-10 | 北京师范大学 | Method for solving contradiction between ecological water and agricultural water |
CN103544533A (en) * | 2012-07-12 | 2014-01-29 | 上海宝信软件股份有限公司 | Diversion project water resource optimal allocation dynamic model establishing system |
Non-Patent Citations (1)
Title |
---|
赵华清等,杨树滩等: "江苏省南水北调受水区需水计算模型初探", 《徐州工程学院学报(自然科学版)》, vol. 27, no. 2, 28 February 2012 (2012-02-28) * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107958095A (en) * | 2016-10-18 | 2018-04-24 | 济南大学 | The parabola shaped open channel of 5/2nds powers of one kind and its optimal hydraulic cross-section |
CN107958095B (en) * | 2016-10-18 | 2020-10-16 | 济南大学 | Bipartite quintuple parabolic open channel and hydraulic optimal section thereof |
CN107944603A (en) * | 2017-11-09 | 2018-04-20 | 中山大学 | Water resource optimal allocation newsboy's method based on water total amount control |
CN108305006A (en) * | 2018-02-02 | 2018-07-20 | 浙江省水利河口研究院 | Reservoir usufruct of water resource amount primary distribution method based on the constraint of water supply fraction |
CN109767143A (en) * | 2019-03-14 | 2019-05-17 | 南京大学 | Water source method for tracing under multi-source water supply conditions |
CN109948235A (en) * | 2019-03-14 | 2019-06-28 | 南京大学 | Water resource scheduling and precision configuration method |
CN109948235B (en) * | 2019-03-14 | 2023-05-09 | 南京大学 | Water resource scheduling and accurate configuration method |
CN112862225A (en) * | 2019-11-27 | 2021-05-28 | 重庆师范大学 | Lake storage regulation excavation-based water shortage area cross-space-time scheduling method and system |
WO2021129039A1 (en) * | 2019-12-25 | 2021-07-01 | 华中科技大学 | Method and system for calculating cross-reservoir drainage basin river runoff |
CN113326632A (en) * | 2021-06-19 | 2021-08-31 | 南昌工程学院 | Optimization correction method for backward-thrust reservoir warehousing flow process |
CN113449890A (en) * | 2021-07-16 | 2021-09-28 | 中国水利水电科学研究院 | Efficient utilization and optimized allocation method for coal mine water area |
Also Published As
Publication number | Publication date |
---|---|
CN104268645B (en) | 2015-06-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104268645A (en) | Industry water supplying quantity determining method with water supplying priority level taken into consideration in water resource allocation | |
CN106168991B (en) | A kind of tidal river network tide prediction method based on hydrodynamic simulation | |
Nagesh Kumar et al. | Optimal reservoir operation for flood control using folded dynamic programming | |
CN108764573B (en) | Inland arid region-oriented multidimensional equilibrium configuration system for water resources | |
CN109948235B (en) | Water resource scheduling and accurate configuration method | |
CN111046574B (en) | Flood control and drainage calculation method for lake and river gate pump system in plain lake region | |
CN105913146B (en) | Southern wet area water resource optimal allocation system | |
CN103544533B (en) | A kind of water diversion project water resource optimal allocation dynamic model constructing system | |
CN103106625A (en) | Reservoir, sluice and pump cluster combined saltwater tide control and scheduling method | |
CN108797511A (en) | A method of needing the Water Resources Allocation of water towards river channel ecology | |
CN113723024A (en) | Method for simulating 'stream' -river channel '-river mouth' distributed flood process suitable for coastal areas | |
Li et al. | Water resources optimal allocation based on large-scale reservoirs in the upper reaches of Yangtze River | |
CN108269199B (en) | Object-oriented small reservoir group time-space distribution type outflow calculation method | |
Mulat et al. | Evaluation of multi-storage hydropower development in the upper Blue Nile River (Ethiopia): regional perspective | |
CN110533326A (en) | A kind of water-related Eco-environmental Carrying Capacity estimates regulator control system | |
CN115470965B (en) | Method for rapidly determining and predicting tide branch estuary branch channel falling tide split ratio based on radial tide countermeasure mode | |
Tang et al. | A numerical study of hydrodynamic processes and flood mitigation in a large river-lake system | |
CN115860362A (en) | Multidimensional collaborative scheduling method for cascade reservoir group | |
Khidirov et al. | Management exploitation condition of Amu-Bukhara machine channel | |
CN106320256A (en) | Method for determining flux of sediment passing through turbine of pumped storage power station for heavily silt-carrying rivers | |
CN117674293B (en) | Long-term power generation optimal scheduling method and device for cascade hydropower station | |
Peng et al. | Positive effect of a canal system and reservoir group on the spatial-temporal redistribution of water resources in a pinnate drainage pattern | |
Ferdowsi et al. | A survey of PSO contributions to water and environmental sciences | |
CN116090744A (en) | Irrigation water allocation method, computer device and storage medium for small irrigation areas in hilly and hilly areas | |
Ikramova et al. | The reservoirs capacity assessment: the Tuyamuyun hydro complex in Khorezm region of Uzbekistan |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |