CN107679317B - A kind of processing method of one-dimensional open channel hydrodynamic model pumping plant inner boundary - Google Patents
A kind of processing method of one-dimensional open channel hydrodynamic model pumping plant inner boundary Download PDFInfo
- Publication number
- CN107679317B CN107679317B CN201710898145.1A CN201710898145A CN107679317B CN 107679317 B CN107679317 B CN 107679317B CN 201710898145 A CN201710898145 A CN 201710898145A CN 107679317 B CN107679317 B CN 107679317B
- Authority
- CN
- China
- Prior art keywords
- pumping plant
- studied
- moment
- flow
- net lift
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Civil Engineering (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Architecture (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of processing methods of one-dimensional open channel hydrodynamic model pumping plant inner boundary, are related to.The method:Obtain the characteristic curve q of every water pump assembly flow and pumping plant net lift in pumping plant to be studiedi~Ht;Obtain the total flow and net lift functional relation Q of k moment pumping plant to be studiedk~Ht,k;Combination coefficient sequence of the pumping plant inner boundary in entire calculation interval;Obtain the compatible equations of one-dimensional open channel hydrodynamic model;Linearization process is carried out to pumping plant boundary compatible equations;The Hydrodynamic Process of open channel where obtaining pumping plant to be studied.The method of the invention improves versatility and the flexibility of one-dimensional open channel hydrodynamic model, promotes the precision of One-dimensional Numerical Simulation.
Description
Technical field
The present invention relates to open-channel flow numerical simulation technology field more particularly to a kind of one-dimensional open channel hydrodynamic model pumping plants
The processing method of inner boundary.
Background technology
Pumping plant is one of most important control building in river network in plain areas, in able to resisting flood and drought, agricultural irrigation, across basin
Water transfer, plumbing etc. play particularly important effect.The water level and flow of flow is adjusted in pumping plant, has broken river network system
The continuity of waterpower relationship.The numerical simulation of open-channel flow is planning, formulates the important evidence of pumping plant scheduling scheme.Therefore, it grinds
Studying carefully processing method of the pumping plant boundary in one-dimensional open channel hydrodynamic model has important value.
In one-dimensional open channel hydrodynamic model (such as single river channel, tree-type pipe network, looping river network), often by pumping plant node boundary
Be converted to flow boundary or water level boundary.Existing part research is using the condition of continuity and conservation of energy condition to pumping plant boundary
Generally changed, pumping plant lift need to be known quantity.Currently, there is no the pumping plant processing method directly against the operation of water pump assembly.Cause
This, the processing method on research pumping plant boundary is one dimensional numerical to improve versatility, the flexibility of one-dimensional open channel hydrodynamic model
The key points and difficulties problem of simulation.
Invention content
The purpose of the present invention is to provide a kind of processing methods of one-dimensional open channel hydrodynamic model pumping plant inner boundary, to solve
Certainly the versatility of one-dimensional river network model and flexibility are poor in the prior art, one obtained using one-dimensional open channel hydrodynamic model
The low problem of Dimensional Simulation of Air Motion result precision.
To achieve the goals above, the processing method of one-dimensional open channel hydrodynamic model pumping plant inner boundary of the present invention, institute
The method of stating includes:
S1 obtains the characteristic curve q of every water pump assembly flow and pumping plant net lift in pumping plant to be studiedi~Ht
According to pumping plant in the flow of each water pump assembly in pumping plant to be studied~lift performance curve and pumping plant to be studied
Necessary head curve of the installation system inquires into every water pump assembly flow qiWith pumping plant net lift HtCharacteristic curve qi~Ht, it is denoted as formula
(1);
I indicates that the number of the water pump assembly in work in pumping plant to be studied, n are the water in working in pumping plant to be studied
Pump assembly quantity, ai、bi、ciFor the characteristic curve q of water pump assembly flow and pumping plant net lift that number is ii~HtBe
Number;
S2 obtains the total flow and net lift functional relation Q of k moment pumping plant to be studiedk~Ht,k
It is bent according to the characteristic of every in running order water pump assembly in any one k moment pumping plant in calculation interval
Line qi~Ht, inquire into obtain the total flow Q in k moment pumping plant to be studiedkWith the net lift H of pumping plant to be studiedt,kFunctional relation Qk
~Ht,k, it is denoted as formula (2):
Ak、Bk、CkFor in the net lift functional relation Q of the total flow of k moment pumping plant to be studied and pumping plant to be studiedk~Ht,k
Coefficient;
S3, combination coefficient sequence of the pumping plant inner boundary in entire calculation interval
On the basis of S2, combination coefficient sequence of the pumping plant inner boundary in entire calculation interval of pumping plant to be studied is obtained,
It is denoted as (A (y), B (y), C (y));Y indicates any one moment in the calculation interval including moment k;
S4 obtains the compatible equations of one-dimensional open channel hydrodynamic model
On the basis of S2 and S3, by pumping plant to be studied any one moment of calculation interval y total flow~net lift
Relationship Qy~Ht,yBe converted to pumping plant to be studied the moment total flow~ga ge relation Qy~Zy, and formed and pumped with continuity equation
It stands the compatible equations on boundary;
S5 carries out linearization process to pumping plant boundary compatible equations;
S6, the Hydrodynamic Process of open channel where solving pumping plant to be studied
It is discrete to the progress of controlling channel equation, the pumping plant boundary compatible equations of linearisation and channel Saint-Venant equation are joined
It is vertical, channel Outer Boundary Conditions are added, complete system of linear equations is formed, solves the hydrodynamic force that this equation group can be obtained channel
Process.
Preferably, according to every water pump in running order in any one k moment pumping plant in calculation interval in S2
The characteristic curve q of uniti~Ht, inquire into the total flow Q in k moment pumping plant to be studiedkWith the net lift H of pumping plant to be studiedt,kLetter
Number relationship Qk~Ht,k, specifically according to following realizations:
S21, it is assumed that k moment pumping plant to be studied has the work of n platform water pump assemblies, in k moment each water pump assembly flow and water
The characteristic curve q of pump installation net lifti~HtRespectively:
S22 marks curve corresponding with each curve equation in formula (3) at one using flow as abscissa, with water pump
Device net lift is then to utilize willfully method in the rectangular coordinate system of ordinate, obtains the pump to be studied under the k moment net lift
The total flow stood;
S23 obtains the corresponding pumping plant total flow to be studied of various discrete net lift, obtains described in S21 and S22
To the total flow and net lift sequence in k moment pumping plant to be studied;
S24, the total flow and net lift sequence for being fitted k moment pumping plant to be studied obtain total stream at pumping plant k moment to be studied
Measure QkWith net lift Ht,kFunctional relation, i.e. formula (2).
Preferably, S3 is specially:The k moment in step S2 is changed to k ' moment of the calculation interval in addition to the k moment, then
Calculate the net lift functional relation Q in the total flow and pumping plant to be studied of k ' moment pumping plant to be studiedk~Ht,k′Curve coefficients;
Wherein, the k ' moment selected every time is all different, if the quantity at k ' moment is m, then
M=mc-mr+1 (4)
In formula, mcIndicate the adjustment number of the pump stations condition to be studied of calculation interval, mrIndicate pump to be studied in calculation interval
The number that reruns for operating mode of standing;
Obtaining each k ' moment is calculated curve coefficients, and the pumping plant inner boundary to be studied that obtains after combination is entirely calculating
The combination coefficient sequence (A (y), B (y), C (y)) of period.
Preferably, in S4, the compatible equations are:
ZoutFor pumping plant water outlet side water level to be studied, ZinFor pumping plant influent side water level to be studied;Q is pumping plant water delivery to be studied
Flow Q=Qin=Qout, QinFor pumping plant influent side flow to be studied, QoutFor pumping plant water outlet side flow to be studied;Pumping plant to be studied
Water delivery flow Q is the total flow Q at any one moment of calculation interval of pumping plant to be studiedy。
Preferably, in S5, linearization process is carried out to pumping plant boundary compatible equations using single order Taylor series expansions.
Preferably, in S6, controlling channel equation is carried out using Preissmann implied formats method discrete.
The method of the invention builds multicomputer using the pumping plant inner boundary in one-dimensional open channel hydrodynamic simulation as research object
The flow of pumping plant-net lift combined characteristic curve, and with continuity equation as inner boundary compatible equations, coupling channel waterpower control
Equation establishes the one-dimensional open channel hydrodynamic model for including pumping plant inner boundary, scientific basis is provided for engineering simulation and actual schedule
And technical support.
Compared with prior art, the present invention has the following advantages and beneficial effect:
(1) prior art introduces pumping plant and needs pump head in pumping plant compatible equations;And the pump that the present invention uses
Net lift of standing is directly related with water level, is easy to couple with channel equation, and facilitate acquisition and calibration.
(2) prior art is when coupling water pump assembly, an only fixed characteristic curve;And each water pump machine in the present invention
The characteristic curve of group can be adjusted switching according to actual operating mode.
(3) prior art is when handling pumping plant boundary, using the scheduling rule of pumping plant, but regulation goal and the energy that actually draws water
Power tends not to fit like a glove;And the present invention uses pumping plant combined characteristic curve processing pumping plant boundary, promotes One-dimensional Numerical Simulation
Precision.
To sum up, the method for the invention improves versatility and the flexibility of one-dimensional open channel hydrodynamic model, is promoted one-dimensional
The precision of numerical simulation.
Description of the drawings
Fig. 1 is the method for the present invention flow chart;
Fig. 2 is water pump assembly flow~net lift performance curve;
Fig. 3 is multicomputer pumping plant flow~net lift composite behaviour curve;
Fig. 4 is village rent land pumping plant water outlet side-preceding osiery pumping plant-low bank between fields head pumping plant water outlet side interval diagram;
Fig. 5 is preceding osiery pumping plant influent side water level process;
Fig. 6 is preceding osiery pumping plant water outlet side water level process;
Fig. 7 is preceding osiery pumping plant discharge process.
Specific implementation mode
In order to make the purpose , technical scheme and advantage of the present invention be clearer, below in conjunction with attached drawing, to the present invention into
Row is further described.It should be appreciated that the specific embodiments described herein are only used to explain the present invention, it is not used to
Limit the present invention.
Detailed description about herein described method:
S1 builds the flow q of any one water pump assembly i in one-dimensional open channel hydrodynamic model pumping plantiIt is net with pumping plant
Lift HtCharacteristic curve.
In the same pumping plant, the same time of running, the actual motion lift phase of the water pump assembly in pumping plant in each work
Together, it is equal to pumping plant lift.Under any rotation angle or rotating speed, the flow q of every water pump assemblyiWith actually raising for the water pump assembly
The characteristic curve q of journey Hi~H is expressed as:
In formula, αi、βi、γiFor the characteristic curve q of water pump assembly iiThree coefficients of~H;Three coefficient values with corner or
The difference of rotating speed and it is different, but three coefficients are when corner or rotating speed determine, are known quantity.
Pumping plant necessary head curve of the installation system HrIt is represented by:
In formula, HstFor pumping plant static lift;SiGo out to water pump assembly i and water pump assembly i to pumping plant for pumping plant intake pool
The drag parameter of the pipeline in pond is known quantity;HtFor pumping plant net lift;It is generally acknowledged that pumping plant water-in and water-out side atmospheric pressure
Power is identical, then there is Hst=Ht。
In the running operating point of water pump and unit, there are H=Hr, integrated (1) and formula (2) obtain formula (3):
Inquired into and can be obtained by formula (3):
Simplified style (4) can obtain separate unit water pump assembly flow qiWith pumping plant net lift HtCharacteristic curve:
In formula:ai=αi;bi=βi;ci=γi-Si.The characteristic curve q of water pump assembly ii~HtSuch as Fig. 2.Due to Ht=
Zout-Zin, ZoutFor pumping plant water outlet side water level to be studied, ZinFor pumping plant influent side water level to be studied, therefore, pumping plant net lift Ht
With pumping plant upstream and downstream water level Z to be studiedout、ZinIt is directly related, facilitate the structure and calibration of hydrodynamic model.
S2, structure pumping plant k moment pumping plant flows QkWith net lift Ht,kFunctional relation.
Assuming that k moment multicomputer pumping plants share work water pump assembly n platforms, each water pump assembly current time, that is, k moment runs
Operating mode flow~net lift characteristic curve is respectively:
It is same using flow as abscissa, with pumping plant net lift HtTo be identified in the rectangular coordinate system of ordinate
Each water pump assembly current time operating condition flow and the net characteristic curve of pumping plant net lift, as shown in figure 3, then utilizing
Willfully method, i.e., the flow on the characteristic curve of each water pump assembly under same net lift are laterally added, obtain pumping under the net lift
The total flow stood.According to the preset or known each discrete net lift needed, gradually calculate under each discrete net lift point
Flow~net lift sequence of pumping plant to be studied can be obtained in flow value, is fitted this sequence and obtains the flow Q at pumping plant k momentkWith it is net
Lift Ht,kFunctional relation is formula (7):
S3, combination coefficient sequence of the generation pumping plant inner boundary in entire calculation interval.
In S2, the total flow of k moment pumping plant to be studied and the net lift functional relation Q of pumping plant to be studied have been obtainedk~
Ht,k′Curve coefficients Ak、Bk、Ck;K ' at the time of the k moment at calculating moment in S2 is changed to calculation interval in addition to k, is then counted
Calculate the net lift functional relation Q in the total flow and pumping plant to be studied of k ' moment pumping plant to be studiedk~Ht,k′, obtaining the k ' moment waits for
Study the net lift functional relation Q of the total flow and pumping plant to be studied of pumping plantk~Ht,k′Curve coefficients Ak′、Bk′、Ck′;I.e.:By k
Be changed to successively in calculation interval after other moment, execute S2 respectively, obtain the total flow of the pumping plant to be studied of different moments with
The net lift functional relation Q of pumping plant to be studiedk~Ht,kCurve coefficients, by obtained all curve coefficients combine after obtain
Combination coefficient sequence (A (y), B (y), C (y)) of the pumping plant inner boundary to be studied in entire calculation interval.
It is illustrated with A (y):A (y) is a sequence, it has several numbers, and what y was marked is this number in sequence A (y)
In position.For example calculation interval is 0-Y, y can be any one moment between 0-Y, that is, any time.
S4, total flow~net lift relationship Q of the pumping plant to be studied at any one moment of calculation intervaly~Ht,yIt is converted to
Total flow~ga ge relation Q of the pumping plant to be studied at the momenty~Zy, and form with flow conservation equation the phase of pumping plant to be studied
Hold equation.QyFor pumping plant to be studied any one moment of calculation interval y total flow;Ht,yBe pumping plant to be studied calculate when
The net lift of section any one moment y, ZyFor pumping plant to be studied any one moment of calculation interval y water level;Y is indicated
Any one moment in calculation interval including moment k.
It can be obtained by S2 and S3, the total flow Q of any time y pumping plants to be studied in calculation intervalyIt is net with pumping plant to be studied
Lift Ht,yFunctional relation Qy~Ht,yIt is expressed as:
Ht,y=A (y)+B (y) Qy+C(y)Qy 2 (9)
Because pumping plant net lift is the water-head of pumping plant water-in and water-out to be studied side, therefore there are formula (10):
Ht,y=Zout,y-Zin,y (10)
Simultaneous (9), (10) two formulas can obtain pumping plant to be studied the moment total flow~ga ge relation Qy~Zy:
Zout,y-Zin,y=A (y)+B (y) Qy+C(y)Qy 2 (11)
In formula, Zout,yFor in moment y pumping plant water outlet side water levels to be studied, Zin,yFor in moment y pumping plant influent sides to be studied
Water level;
Pumping plant flow meets continuity equation, if ignoring the flow loss inside pumping plant to be studied, exists:
Qin=Qout (12)
In formula, QinFor pumping plant influent side flow to be studied, QoutFor water outlet side flow to be studied.
Formula (11) pumping plant flow~ga ge relation and formula (12) pumping plant water-in and water-out effluent magnitude relation are collectively as one as a result,
Tie up the consistency condition of pumping plant inner boundary in open channel numerical simulation.
S5 carries out linearization process using single order Taylor series expansions to pumping plant boundary compatible equations
In one-dimensional open channel numerical simulation, it is assumed that pumping plant influent side is s sections, and water outlet side is s+1 sections, then the pump in S4
Boundary compatible equations of standing are expressed as:
To each unknown quantity Q in formula (13)s, Qs+1, Zs, Zs+1, single order Taylor series expansions are used respectively:
It is known quantity with the variate-value that target on * is a upper cycle in formula (14);Band △ be current variable values relatively on
The increment of one loop variant values.
Formula (14), which is substituted into formula (13), to be obtained:
It transplants respectively to two item formulas of equation group (15), left side is the incremental form of each variable, and right side is known quantity shape
Formula, meanwhile, the quadratic term of increment is ignored in formula (15), then obtains pumping plant compatible equations (16):
S6, the Hydrodynamic Process of open channel where obtaining pumping plant to be studied
Using Preissmann formats to controlling channel equation carry out it is discrete, by the pumping plant node compatible equations of linearisation with
Channel Saint Venant equations simultaneousnesses, in addition channel entrance boundary condition and channel export boundary condition, form complete line
Property equation group, the Hydrodynamic Process that this equation group can be obtained open channel is solved with chasing method.
In the present embodiment, it is assumed that channel shares x section, and channel is divided x-1 canal section, includes 2x variable altogether.It adopts
With the discrete channel Saint Venant equations of Preissmann formats, and linearization process is carried out, channel jth section and jth+1
Controlling channel equation between section is represented by formula (17):
Each coefficient can be by the hydraulic elements table of known range survey data, known hydraulic parameters and last moment in formula (17)
Show.
Assuming that pumping plant influent side is set as s sections, water outlet side is set as s+1 sections, is pumped in pumping plant compatible equations (16)
It stands on the basis of inner boundary compatible equations, each coefficient can be expressed as in controlling channel equation (17):
Controlling channel equation (17), pumping plant inner boundary compatible equations (16), channel entrance boundary condition (23) and channel go out
Mouth boundary condition (24) amounts to 2x-2 equations, can form the closing Algebraic Equation set containing 2x equation.Because channel may
There are several small ditch sections.Each canal section has one group of controlling channel equation (17), in addition the pumping plant inner edge of any one small ditch section
Boundary's compatible equations (16), can form Algebraic Equation set.
If channel upstream entrance boundary uses flow boundary, such as formula (19), channel lower exit boundary uses water level boundary,
Such as formula (20).Then have:
Q1=Q (y) (19)
Z2x=Z (y) (20)
Q1In first section flow, Z2xIndicate the water level of the 2x section;2x indicates the 2x section of channel,
It is exactly the last one section, t indicates that water level is the sequence changed over time.
Single order Taylor series expansions are done respectively to formula (19) and formula (20), are obtained:
In formula,Headed by section (w+1) time Δt flow value;For the water of terminal surface (w+1) time Δt
Place value.To be unified with controlling channel equation, formula (19) and formula (20) are written as form:
c1ΔZ1+d1ΔQ1=e1 (23)
a2xΔZx+b2xΔQx=e2x (24)
Each coefficient is in formula:c1=0;d1=1;a2x=1;b2x=0;
Being write closed Algebraic Equation set as matrix form is:AX=D.
Using chasing method solution matrix, you can obtain the unsteady flow process containing pumping plant inner boundary open channel.Matrix Solving
Method belongs to prior art, and details are not described herein again.
Embodiment
The present embodiment is right according to description method of the present invention using osiery pumping plant before Miyun Reservoir storage project as inner boundary
It collects rent land pumping plant water outlet side-preceding osiery pumping plant-low bank between fields head pumping plant influent side section (see Fig. 4) and carries out unsteady flow simulation.
Miyun Reservoir storage project is the conveyance system of Middle Line Project For Transferring Water From South To North Beijing Section, 103 kilometers of overall length, with Huairou reservoir
For boundary, it is divided into forward and backward two sections.Wherein, Round Castle lake to Huairou section, altogether include 6 pumping plants, respectively collect rent land, preceding osiery, low bank between fields head,
Pumping plant on Xing Shou, Li Shishan, West Platform is connected with open channel between grade.Pumping plant unit is three standby with one, and preceding 5 grades of pumping plants use axis
Flow pump machine group, the 6th grade of pumping plant use mixed flow pump assembly.Pumping plant is the major regulatory building of system, and Numerical-Mode is carried out to system
It is quasi-, pumping plant can be instructed to dispatch, it is ensured that the safe operation of whole system.
Simulated time be 25 days 0 April in 2016 when 29 days 24 April in 2016, totally 5 days, 120 hours.Upstream boundary
To collect rent land pumping plant water outlet side discharge process, downstream is low bank between fields head pumping plant influent side water level process.Preceding osiery pumping plant separate unit water pump assembly
Flow qiWith pumping plant net lift HtCharacteristic curve expression formula be shown in Table 1.The build-up curve coefficient of osiery pumping plant before the simulation period
Sequence is shown in Table 2.Preceding osiery pumping plant hydrodynamic simulation result is shown in Fig. 5-7.
Each corner q of osiery water pump in pump station unit before table 1i~HtCurve representation formula
The build-up curve coefficient series of osiery pumping plant before table 2
By using above-mentioned technical proposal disclosed by the invention, following beneficial effect has been obtained:
(1) prior art introduces pumping plant and needs pump head in pumping plant compatible equations;And the pump that the present invention uses
Net lift of standing is directly related with water level, is easy to couple with channel equation, and facilitate acquisition and calibration.
(2) prior art is when coupling water pump assembly, an only fixed characteristic curve;And each water pump machine in the present invention
The characteristic curve of group can be adjusted switching according to actual operating mode.
(3) prior art is when handling pumping plant boundary, using the scheduling rule of pumping plant, but regulation goal and the energy that actually draws water
Power tends not to fit like a glove;And the present invention uses pumping plant combined characteristic curve processing pumping plant boundary, promotes One-dimensional Numerical Simulation
Precision.
To sum up, the method for the invention improves versatility and the flexibility of one-dimensional open channel hydrodynamic model, is promoted one-dimensional
The precision of numerical simulation.
The above is only a preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art
For member, various improvements and modifications may be made without departing from the principle of the present invention, these improvements and modifications are also answered
Depending on protection scope of the present invention.
Claims (6)
1. a kind of processing method of one-dimensional open channel hydrodynamic model pumping plant inner boundary, which is characterized in that the method includes:
S1 obtains the characteristic curve q of every water pump assembly flow and pumping plant net lift in pumping plant to be studiedi~Ht
According to pumping plant needs in the flow of each water pump assembly in pumping plant to be studied~lift performance curve and pumping plant to be studied
Head curve inquires into every water pump assembly flow qiWith pumping plant net lift HtCharacteristic curve qi~Ht, it is denoted as formula (1);
I indicates that the number of the water pump assembly in work in pumping plant to be studied, n are the water pump machine in working in pumping plant to be studied
Group quantity, ai、bi、ciFor the characteristic curve q of water pump assembly flow and pumping plant net lift that number is ii~HtCoefficient;
S2 obtains the total flow and net lift functional relation Q of k moment pumping plant to be studiedk~Ht,k
According to the characteristic curve q of every in running order water pump assembly in any one k moment pumping plant in calculation intervali
~Ht, inquire into obtain the total flow Q in k moment pumping plant to be studiedkWith the net lift H of pumping plant to be studiedt,kFunctional relation Qk~
Ht,k, it is denoted as formula (2):
Ak、Bk、CkFor in the net lift functional relation Q of the total flow of k moment pumping plant to be studied and pumping plant to be studiedk~Ht,kBe
Number;
S3, combination coefficient sequence of the pumping plant inner boundary in entire calculation interval
On the basis of S2, combination coefficient sequence of the pumping plant inner boundary in entire calculation interval of pumping plant to be studied is obtained, is denoted as
(A (y), B (y), C (y));Y indicates any one moment in the calculation interval including moment k;
S4 obtains the compatible equations of one-dimensional open channel hydrodynamic model
On the basis of S2 and S3, by pumping plant to be studied any one moment of calculation interval y total flow~net lift relationship Qy
~Ht,yBe converted to pumping plant to be studied the moment total flow~ga ge relation Qy~Zy, and form pumping plant boundary with continuity equation
Compatible equations;
S5 carries out linearization process to pumping plant boundary compatible equations;
S6, the Hydrodynamic Process of open channel where solving pumping plant to be studied
It is discrete to the progress of controlling channel equation, by the pumping plant boundary compatible equations of linearisation and channel Saint-Venant equation simultaneous, then
In addition channel Outer Boundary Conditions, form complete system of linear equations, the Hydrodynamic Process that this equation group can be obtained channel is solved.
2. method according to claim 1, which is characterized in that according to any one k moment pumping plant in calculation interval in S2
In every in running order water pump assembly characteristic curve qi~Ht, inquire into the total flow Q in k moment pumping plant to be studiedkWith
The net lift H of pumping plant to be studiedt,kFunctional relation Qk~Ht,k, specifically according to following realizations:
S21, it is assumed that k moment pumping plant to be studied has the work of n platform water pump assemblies, is filled in k moment each water pump assembly flow and water pump
Set the characteristic curve q of net lifti~HtRespectively:
S22 marks curve corresponding with each curve equation in formula (3) at one using flow as abscissa, with pumping plant
Net lift is then to utilize willfully method in the rectangular coordinate system of ordinate, obtains the pumping plant to be studied under the k moment net lift
Total flow;
S23 obtains the corresponding pumping plant total flow to be studied of various discrete net lift, obtains in k described in S21 and S22
The total flow of moment pumping plant to be studied and net lift sequence;
S24, the total flow and net lift sequence for being fitted k moment pumping plant to be studied obtain the total flow Q at pumping plant k moment to be studiedkWith
Net lift Ht,kFunctional relation, i.e. formula (2).
3. method according to claim 1, which is characterized in that S3 is specially:
The k moment in step S2 is changed to k ' moment of the calculation interval in addition to the k moment, is then calculated in k ' moment pump to be studied
The net lift functional relation Q of the total flow and pumping plant to be studied stoodk′~Ht,k′Curve coefficients;Wherein, when the k ' selected every time
Quarter is all different, if the quantity at k ' moment is m, then
M=mc-mr+1 (4)
In formula, mcIndicate the adjustment number of the pump stations condition to be studied of calculation interval, mrIndicate pumping plant work to be studied in calculation interval
The number that reruns of condition;
Obtaining each k ' moment is calculated curve coefficients, and pumping plant inner boundary to be studied is obtained in entire calculation interval after combination
Combination coefficient sequence (A (y), B (y), C (y)).
4. method according to claim 1, which is characterized in that in S4, the compatible equations are:
ZoutFor pumping plant water outlet side water level to be studied, ZinFor pumping plant influent side water level to be studied;Q is pumping plant water delivery flow Q to be studied
=Qin=Qout, QinFor pumping plant influent side flow to be studied, QoutFor pumping plant water outlet side flow to be studied;Pumping plant water delivery to be studied
Flow Q is the total flow Q at any one moment of calculation interval of pumping plant to be studiedy。
5. method according to claim 1, which is characterized in that in S5, using single order Taylor series expansions to pumping plant side
Boundary's compatible equations carry out linearization process.
6. method according to claim 1, which is characterized in that in S6, using Preissmann implied format methods to channel
Governing equation carries out discrete.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710898145.1A CN107679317B (en) | 2017-09-28 | 2017-09-28 | A kind of processing method of one-dimensional open channel hydrodynamic model pumping plant inner boundary |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710898145.1A CN107679317B (en) | 2017-09-28 | 2017-09-28 | A kind of processing method of one-dimensional open channel hydrodynamic model pumping plant inner boundary |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107679317A CN107679317A (en) | 2018-02-09 |
CN107679317B true CN107679317B (en) | 2018-08-21 |
Family
ID=61138388
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710898145.1A Active CN107679317B (en) | 2017-09-28 | 2017-09-28 | A kind of processing method of one-dimensional open channel hydrodynamic model pumping plant inner boundary |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107679317B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109992855B (en) * | 2019-03-19 | 2023-07-25 | 黑龙江省水利水电勘测设计研究院(水利部黑龙江水利水电勘测设计研究院、黑龙江省水利水电规划院) | Solving method of river network one-dimensional constant flow based on river network association matrix |
CN111783369B (en) * | 2020-07-22 | 2024-01-26 | 中国水利水电科学研究院 | Short-term multi-objective optimal scheduling method for multi-gate-group open channel water diversion project |
CN112301962B (en) * | 2020-10-21 | 2021-12-17 | 长江勘测规划设计研究有限责任公司 | Method for improving flow capacity of large diversion project without affecting normal water supply |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2624547C (en) * | 2005-10-03 | 2016-02-16 | Central Sprinkler Company | System and method for evaluation of fluid flow in a piping system |
CN102722611B (en) * | 2012-05-29 | 2014-09-10 | 清华大学 | Method for carrying out parallelization numerical simulation on hydrodynamic force conditions of river provided with cascade hydropower station |
CN105005707B (en) * | 2015-08-05 | 2017-08-08 | 武汉大学 | A kind of stage discharge relation of Heavenly Stems and Earthly Branches stream confluence determines method |
-
2017
- 2017-09-28 CN CN201710898145.1A patent/CN107679317B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN107679317A (en) | 2018-02-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107679317B (en) | A kind of processing method of one-dimensional open channel hydrodynamic model pumping plant inner boundary | |
Watts et al. | A soil‐water‐nitrogen model for irrigated corn on sandy soils | |
CN105678067B (en) | Urban waterway is controlled flood and drained flooded fields dynamic forecasting control method and system | |
CN107306765A (en) | A kind of water-fertilizer integral irrigation system and irrigation method | |
CN107633141A (en) | A kind of discrimination method of one-dimension mathematical model pumping station performance curve coefficients | |
Liu et al. | Integrated modeling of conjunctive water use in a canal-well irrigation district in the lower Yellow River basin, China | |
Barros et al. | Long-term water balances in La Violada Irrigation District (Spain): II. Analysis of irrigation performance | |
CN108575673A (en) | Drought-hit area crop irrigation fertilizing method based on same day Weather Forecast Information and system | |
CN106934534A (en) | Rattan knot melon Irrigation Project Design water effective utilization coefficients computational methods long | |
CN102359146A (en) | Selective inlet water temperature stabilizing device and water temperature value simulating and forecasting method | |
CN106498896A (en) | A kind of reservoir operation process analogy method being coupled with Sedimentation | |
CN110826805B (en) | Low-water-head cascade hydropower station medium-term optimization scheduling method considering water unevenness | |
CN109258044A (en) | A kind of field water fertilizer managing and control system and method based on biogas slurry | |
CN104916049B (en) | A kind of motor-pumped well water metering system | |
CN114819753B (en) | Cascade reservoir ecological scheduling method for promoting fish spawning | |
CN110232479A (en) | A kind of city flood control by reservoir regulation compensation optimizing dispatching method | |
Frederick | Paper 4. Climate change impacts on water resources and possible responses in the MINK region | |
CN108647907B (en) | Cascade pumping station open channel water diversion project water-filling method based on lock pump combined regulating | |
CN210597185U (en) | Multifunctional canal system test system for laboratory | |
Plusquellec | Improving the operation of canal irrigation systems | |
CN113435631B (en) | Flood forecasting method, flood forecasting system, readable storage medium and computing device | |
CN114757405A (en) | Irrigated area water resource balance optimization configuration method | |
CN112990757A (en) | Basin dry-running dry-water year ecological water transfer method | |
Shatanawi et al. | Irrigation system performance in Jordan | |
Allen et al. | Concepts for energy-efficient irrigation system design |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |