CN104268645A - Industry water supplying quantity determining method with water supplying priority level taken into consideration in water resource allocation - Google Patents

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

The industry output defining method of water supply precedence is considered in Water Resources Allocation

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:

$\frac{\∂Q}{\∂x}+B\frac{\∂Z}{\∂t}=q_L-{\mathrm{\δ}}_{*},$

$\frac{\∂Q}{\∂t}+2u\frac{\∂Q}{\∂x}+(\mathrm{gA}-{\mathrm{Bu}}^2)\frac{\∂Z}{\∂x}+\frac{{\mathrm{gn}}^2\left|u\right|Q}{R^{4/3}}=0,$

Wherein, Q is flow, and unit is m ³/ 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 ²/ 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 ²; 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:

$x_i^{(k+1)}=(1-\mathrm{\ω})x_i^{\left(k\right)}+\frac{\mathrm{\ω}}{a_{\mathrm{ii}}}(b_i-\underset{j=1}{\overset{i-1}{\mathrm{\Σ}}}{a}_{\mathrm{ij}}{x}_j^{(k+1)}-\underset{j=i+1}{\overset{n}{\mathrm{\Σ}}}{a}_{\mathrm{ij}}{x}_j^{\left(k\right)}),$

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 ₁,Z ₂,Z ₃}

In formula, Z _eriver channel ecology water level, Z ₁for the water level calculated according to method of maximum curvature or the Slope Method of wetted perimeter method, Z ₂for the river water level according to minimum existence depth of water standard or 60% wetted perimeter rate criterion calculation, Z ₃for 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:

$\frac{\∂Q}{\∂x}+B\frac{\∂Z}{\∂t}=q_L-{\mathrm{\δ}}_{*},$

$\frac{\∂Q}{\∂t}+2u\frac{\∂Q}{\∂x}+(\mathrm{gA}-{\mathrm{Bu}}^2)\frac{\∂Z}{\∂x}+\frac{{\mathrm{gn}}^2\left|u\right|Q}{R^{4/3}}=0,$

Wherein, Q is flow, and unit is m ³/ 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 ²/ 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 ²; 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:

$x_i^{(k+1)}=(1-\mathrm{\ω})x_i^{\left(k\right)}+\frac{\mathrm{\ω}}{a_{\mathrm{ii}}}(b_i-\underset{j=1}{\overset{i-1}{\mathrm{\Σ}}}{a}_{\mathrm{ij}}{x}_j^{(k+1)}-\underset{j=i+1}{\overset{n}{\mathrm{\Σ}}}{a}_{\mathrm{ij}}{x}_j^{\left(k\right)}),$

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 ₁, Q ₂, Q ₃for river course 1,2,3 flow, Q _gfor weir flow amount, Q ₄for 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:

$\left.\begin{array}{c}{Q}_1=b_{11}+b_{12}{Z}_{N1}+b_{13}+Z_{N2}\\ Q_2=b_{21}+b_{22}{Z}_{N1}+b_{23}{Z}_{N3}\\ Q_3=a_{31}+a_{32}{Z}_{N1}+a_{33}{Z}_{N4}\\ Q_g=f(Z_{N1},Z_{N5})\end{array}\right\}---\left(33\right)$

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 ₁(Z _N1-H ₀)， $C_1=\mathrm{\μ\ϵBk}\sqrt{2g}\sqrt{H_u-H_0}$

Under submerge discharging flow pattern:

Q _g＝C ₂(Z _N1-Z _N5)， $C_2=\mathrm{\φ\ϵBk}\sqrt{2g}\frac{H_d-H_0}{\sqrt{H_u-H_d}}$

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 ₀for 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:

$Q_1+Q_2-Q_3-Q_g+Q_4=A\frac{Z_{N1}-Z_{N1}^j}{\mathrm{\Δt}}---\left(34\right)$

A is node area, and unit is m ², Δ 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 ₁₃, K _n1.3=b ₂₃, K _n1.4=-a ₃₃, $K_{N1.R}=a_{31}-b_{11}-b_{21}-C_1{H}_0-Q_4-A\frac{Z_{N1}^j}{\mathrm{\Δt}};$

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 ₁₃, K _n1.3=b ₂₃, K _n1.4=-a ₃₃, K _n1.5=C ₂ $K_{N1.R}=a_{31}-b_{11}-b_{21}-Q_4-A\frac{Z_{N1}^j}{\mathrm{\Δt}},$

K in formula _n1,1, K _n1,2, K _n1,3, K _n1,4, K _n1,5, K _{n1, R}for 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.

$C_i=\frac{{\mathrm{\Δx}}_i}{2\mathrm{\Δt}}{B}_{i+1/2}^j$

$D_i=q_L\mathrm{\Δx}+C_i(Z_i^j+Z_{i+1}^j)-\frac{S}{2}\frac{{\left(\mathrm{\Δ}{Z}_{*(i+1/2)}\right)}^2}{\mathrm{\Δt}}\mathrm{\Δ}{x}_i$

$E_i=\frac{\mathrm{\Δ}{x}_i}{2\mathrm{\Δt}}-2u_{i+1/2}^j+\frac{g}{2}{\left(\frac{n^2\left|u\right|}{R^{4/3}}\right)}_i^j\mathrm{\Δ}{x}_i$

$G_i=\frac{\mathrm{\Δ}{x}_i}{2\mathrm{\Δt}}-2u_{i+1/2}^j+\frac{g}{2}{\left(\frac{n^2\left|u\right|}{R^{4/3}}\right)}_{i+1}^j\mathrm{\Δ}{x}_i$

$F_i={(\mathrm{gA}-{\mathrm{Bu}}^2)}_{i+1/2}^j$

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 ₁linear function: Q _i=θ _i+ η _iz _i+ ν _iz ₁

${\mathrm{\η}}_i=\frac{F_{i-1}{Y}_1-C_{i-1}{Y}_2}{W}$

${\mathrm{\ν}}_i=\frac{{\mathrm{\ν}}_{i-1}(Y_2+E_{i-1}{Y}_1)}{W}$

Y ₁＝C _i-1-η _i-1

Y ₂＝E _i-1η _i-1-F _i-1

W＝Y ₂-Y ₁G _i-1

Recursion initial value is calculated by following formula:

${\mathrm{\η}}_2=-\frac{C_1{E}_1+F_1}{W}$

${\mathrm{\ν}}_2=-\frac{F_1-C_1{E}_1}{W}$

W＝E ₁+G ₁

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

${\mathrm{\β}}_i=\frac{Y_2{C}_i+Y_1{F}_i}{W}$

${\mathrm{\ζ}}_i=\frac{{\mathrm{\ζ}}_{i+1}(Y_2-Y_1{G}_i)}{W}$

Y ₁＝C _i+β _i+1

Y ₂＝G _iβ _i+1+F _i

W＝Y ₂+Y ₁E _i

Recursion initial value is calculated by following formula:

${\mathrm{\β}}_{n-1}=\frac{1}{W}(C_{n-1}{G}_{n-1}+F_{n-1})$

${\mathrm{\ζ}}_{n-1}=\frac{1}{W}(C_{n-1}{G}_{n-1}-F_{n-1})$

W＝G _n-1+E _n-1

From section second from the bottom, recursion obtains α to first section forward successively ₁, β ₁, ζ ₁, then by α ₁, β ₁, ζ ₁corresponding 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 ₁, Z ₂z ₁₁be 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 ₁, Z ₂, Z ₃the 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.

$C_i=\frac{{\mathrm{\Δx}}_i}{2\mathrm{\Δt}}{B}_{i+1/2}^j$

$D_i=q_L\mathrm{\Δx}+C_i(Z_i^j+Z_{i+1}^j)-\frac{S}{2}\frac{{\left(\mathrm{\Δ}{Z}_{*(i+1/2)}\right)}^2}{\mathrm{\Δt}}\mathrm{\Δ}{x}_i$

$E_i=\frac{\mathrm{\Δ}{x}_i}{2\mathrm{\Δt}}-2u_{i+1/2}^j+\frac{g}{2}{\left(\frac{n^2\left|u\right|}{R^{4/3}}\right)}_i^j\mathrm{\Δ}{x}_i$

$G_i=\frac{\mathrm{\Δ}{x}_i}{2\mathrm{\Δt}}-2u_{i+1/2}^j+\frac{g}{2}{\left(\frac{n^2\left|u\right|}{R^{4/3}}\right)}_{i+1}^j\mathrm{\Δ}{x}_i$

$F_i={(\mathrm{gA}-{\mathrm{Bu}}^2)}_{i+1/2}^j$

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 ₁

${\mathrm{\η}}_i=\frac{F_{i-1}{Y}_1-C_{i-1}{Y}_2}{W}$

${\mathrm{\ν}}_i=\frac{{\mathrm{\ν}}_{i-1}(Y_2+E_{i-1}{Y}_1)}{W}$

Y ₁＝C _i-1-η _i-1

Y ₂＝E _i-1η _i-1-F _i-1

W＝Y ₂-Y ₁G _i-1

Recursion initial value is calculated by following formula:

${\mathrm{\η}}_2=-\frac{C_1{E}_1+F_1}{W}$

${\mathrm{\ν}}_2=-\frac{F_1-C_1{E}_1}{W}$

W＝E ₁+G ₁

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

${\mathrm{\β}}_i=\frac{Y_2{C}_i+Y_1{F}_i}{W}$

${\mathrm{\ζ}}_i=\frac{{\mathrm{\ζ}}_{i+1}(Y_2-Y_1{G}_i)}{W}$

Y ₁＝C _i+β _i+1

Y ₂＝G _iβ _i+1+F _i

W＝Y ₂+Y ₁E _i

Recursion initial value is calculated by following formula:

${\mathrm{\β}}_{n-1}=\frac{1}{W}(C_{n-1}{G}_{n-1}+F_{n-1})$

${\mathrm{\ζ}}_{n-1}=\frac{1}{W}(C_{n-1}{G}_{n-1}-F_{n-1})$

W＝G _n-1+E _n-1

From section second from the bottom, recursion obtains α to first section 1-2 forward successively ₁, β ₁, ζ ₁, then by α ₁, β ₁, ζ ₁corresponding 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 ₁, Q ₂, Q ₃for river course 1,2,3 flow, Q _gfor weir flow amount, Q ₄for directly joining node N ₁flow, 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.

$\left.\begin{array}{c}{Q}_1=b_{11}+b_{12}{Z}_{N1}+b_{13}+Z_{N2}\\ Q_2=b_{21}+b_{22}{Z}_{N1}+b_{23}{Z}_{N3}\\ Q_3=a_{31}+a_{32}{Z}_{N1}+a_{33}{Z}_{N4}\\ Q_g=f(Z_{N1},Z_{N5})\end{array}\right\}---\left(33\right)$

Wherein, under free discharge pattern:

Q _g＝C ₁(Z _N1-H ₀)， $C_1=\mathrm{\μ\ϵBk}\sqrt{2g}\sqrt{H_u^j-H_0}$

Under submerge discharging flow pattern:

Q _g＝C ₂(Z _N1-Z _NS)， $C_2=\mathrm{\φ\ϵBk}\sqrt{2g}\frac{H_d^j-H_0^j}{\sqrt{H_u^j-H_d^j}}$

Node N1 water balance equation is:

$Q_1+Q_2-Q_3-Q_g+Q_4=A\frac{Z_{N1}-Z_{N1}^j}{\mathrm{\Δt}}---\left(34\right)$

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 ₁₃, K _n1.3=b ₂₃, K _n1.4=-a ₃₃, $K_{N1.R}=a_{31}-b_{11}-b_{21}-C_1{H}_0-Q_4-A\frac{Z_{N1}^j}{\mathrm{\Δt}};$

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 ₁₃, K _n1.3=b ₂₃, K _n1.4=-a ₃₃, K _n1.5=C ₂ $K_{N1.R}=a_{31}-b_{11}-b_{21}-Q_4-A\frac{Z_{N1}^j}{\mathrm{\Δt}},$

K in formula _n1,1, K _n1,2, K _n1,3, K _n1,4, K _n1,5, K _{n1, R}all 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 ₁, Z ₂z ₁₁be 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:

$x_i^{(k+1)}=(1-\mathrm{\ω})x_i^{\left(k\right)}+\frac{\mathrm{\ω}}{a_{\mathrm{ii}}}(b_i-\underset{j=1}{\overset{i-1}{\mathrm{\Σ}}}{a}_{\mathrm{ij}}{x}_j^{(k+1)}-\underset{j=i+1}{\overset{n}{\mathrm{\Σ}}}{a}_{\mathrm{ij}}{x}_j^{\left(k\right)}),$ I=1,2 ..., n, n are the network of waterways, region node serial number; Relaxation factor ω=1.25, meet iterations K > 200 or adjacent 2 iteration water level value x _iabsolute error is less than 0.01, solves R9 river course and other 202 river course each node present situation nineteen eighty-three types of survey region first 15 minutes last water level Z=x (unit: m), then employing formula Q by step iteration during 15min _i=θ _i+ η _iz _i+ ν _iz ₁, Q _i=α _i+ β _iz _i+ ζ _iz _n, i is every bar river cross-section numbering, same river course same section simultaneous solution is obtained to present situation nineteen eighty-three type first 15 minutes last water level Z of each section in every bar river course _i(unit: m), flow Q _i(unit: m ³/ S).R9 river course first and last node present situation nineteen eighty-three type first 15 minutes last water level as tried to achieve is respectively 4.03m, 5.22m, mean water 4.625m.

(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 ³, WI _mouthful=0m ³, WS _mouthful=591m ³, WE _mouthful=457m ³, WA _mouthful=104m ³.

(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 ³, 0m ³, 591m ³, 457m ³, 104m ³, 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 ³, 000,000,000 m ³, 0.21 hundred million m ³, 0.16 hundred million m ³, 0.62 hundred million m ³, 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:

$\frac{\∂Q}{\∂x}+B\frac{\∂Z}{\∂t}=q_L-{\mathrm{\δ}}_{*},$

$\frac{\∂Q}{\∂t}+2u\frac{\∂Q}{\∂x}+(\mathrm{gA}-{\mathrm{Bu}}^2)\frac{\∂Z}{\∂x}+\frac{{\mathrm{gn}}^2\left|u\right|Q}{R^{4/3}}=0,$

Wherein, Q is flow, and unit is m ³/ 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 ²/ 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 ²; 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:

$x_i^{(k+1)}=(1-\mathrm{\ω})x_i^{\left(k\right)}+\frac{\mathrm{\ω}}{a_{\mathrm{ii}}}(b_i-\underset{j=1}{\overset{i-1}{\mathrm{\Σ}}}{a}_{\mathrm{ij}}{x}_j^{(k+1)}-\underset{j=i+1}{\overset{n}{\mathrm{\Σ}}}{a}_{\mathrm{ij}}{x}_j^{\left(k\right)}),$

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.