CN105089003A - Reservoir flood routing calculating method - Google Patents

Abstract

The invention relates to the field of hydraulic and hydro-power engineering and discloses a reservoir flood routing calculating method in order to solve the problems of low calculation efficiency, high time consumption and low calculation precision existing in the traditional technology. The reservoir flood routing calculating method comprises the steps of reading related information of a known reservoir, and calculating the average reservoir inflow within a certain time interval; finding the slope k=S'(W0) of a reservoir capacity and discharge capacity curve at the reservoir capacity W0 according to the initial W0 within the time interval; directly calculating the final reservoir capacity W1 within the time interval according to the formula shown in the specification; directly calculating the final discharge capacity O1 of the reservoir within the time interval according to the formula as shown in the specification; finding a reservoir level and capacity curve to obtain the final water level of the reservoir within the time interval; calculating the water level, reservoir capacity and discharge capacity within the next time interval by taking the final water level, reservoir capacity and discharge capacity calculated within the time interval as the initial water level, reservoir capacity and discharge capacity within the next time interval until finishing the calculation within all time intervals; and outputting the calculated water level series, reservoir capacity series and discharge capacity series. The reservoir flood routing calculating method is free of trial calculation, simple and convenient in calculation process, relatively high in precision, high in solving speed, easy to realize by programming and suitable for reservoir flood routing calculation in the hydraulic and hydro-power engineering.

Description

A kind of storage routing for reservoir method

Technical field

The present invention relates to Hydraulic and Hydro-Power Engineering field, be specifically related to a kind of storage routing for reservoir method.

Background technology

Storage routing for reservoir is very important in Hydraulic and Hydro-Power Engineering prospective design, operational management and the work on basis, the master tool determined Hydraulic and Hydro-Power Engineering scale He carry out reservoir scientific dispatch, therefore improve and improve adjust flood calculate precision and computational efficiency significant.

The normally known Flood process of reservoir of storage routing for reservoir, water level storage-capacity curve, discharge capacity curve, and the initial water level of reservoir, and then inquire into discharge flow process and reservoir level process, the fundamental equation inquired into is water balance equation, see formula (1), its implication is reservoir inflow and deducts the variable quantity that earial drainage (outbound) flow equals storage capacity.

$I\left(t\right)-O\left(t\right)=\frac{dW}{dt}---\left(1\right)$

Wherein: O=q (Z) (2)

W＝V(Z)(3)

O＝S(W)(4)

In formula:

I (t) ... reservoir inflow process;

O (t) ... earial drainage (outbound) discharge process;

W ... reservoir filling storage capacity, represents the water retention capacity in reservoir;

Z ... reservoir level;

Q (Z) ... discharge capacity and reservoir level relation are monotone increasing curve;

V (Z) ... reservoir capacity and reservoir level relation is monotone increasing curve;

S (W) ... discharge capacity and reservoir capacity relation, can be converted to by formula (2), (3), be all monotone increasing curve.

It is exactly solving above-mentioned equation in fact that the tune flood of reservoir calculates, and because above-mentioned formula (2), (3), (4) are non-explicit expression, strict economics analysis method therefore cannot be adopted to solve.At present for solving the problem of storage routing for reservoir, the method ills method, iterative method, numerical solution etc. of main application, wherein graphical method needs mapping, be unfavorable for that computer programming realizes, applied less at present, and iterative method and numerical solution are applied more at present, are now described below:

Iterative method:

First iterative method carries out difference simplify processes to equation (1), the step-length that design is calculated is Δ t, then formula (1) is variable is changed to differential type (5), bring formula (2) and formula (3) into formula (5), iterative method design formulas (6) can be obtained.Due to warehouse-in process I (t) and initial water level Z _t=0known, t is given for calculation interval Δ, then every single-step iteration tentative calculation calculation procedure is as follows: first hypothesis period Mo water level is then bring the right of formula (6) into, try to achieve last storage capacity according to formula (3), the last water level of calculating can be tried to achieve then judge whether be less than allowable error, if then illustrate that hypothesis period Mo water level is the last water level of this period calculating, otherwise, again suppose period Mo water level, repeat said process and calculate, till supposing that the difference of last water level with the last water level of calculating is less than allowable error.Flood is adjusted to calculate according to above-mentioned computational methods, from initial time, calculate until terminate, reservoir capacity process W (t) after the tune flood finally can asked, water level process Z (t), earial drainage (outbound) discharge process O (t) etc. by the period.

$W_{t+\mathrm{\Δ}t}=(\frac{I_{t+\mathrm{\Δ}t}+I_t}{2}-\frac{O_{t+\mathrm{\Δ}t}+O_t}{2})\mathrm{\Δ}t+W_t---\left(5\right)$

$W_{t+\mathrm{\Δ}t}=(\frac{I_{t+\mathrm{\Δ}t}+I_t}{2}-\frac{q\left(Z_{t+\mathrm{\Δ}t}\right)+q\left(Z_t\right)}{2})\mathrm{\Δ}t+V\left(Z_t\right)---\left(6\right)$

In formula:

Δ t ... for calculation interval, i.e. material calculation;

Time in formula (1), (2), (3) adopts subscript to represent in formula (5), (6); All the other parameter meanings are the same.

Numerical solution:

Numerical solution brings formula (4) into formula (1), obtains the differential equation about reservoir filling storage capacity W, then according to conservation storage W at the beginning of the known period _t=0=W ₀, form Solution of Initial Value Problems of First-order Ordinary Differential, see formula (7), then solved by numerical solution.

$\left\{\begin{array}{c}{W}^{\′}=\frac{dW}{dt}=I\left(t\right)-S\left(W\right)\\ W(t=0)=W_0\end{array}\right.---\left(7\right)$

The numerical solution of formula (7) has a lot, and for the mainly Fourth order Runge-Kutta of storage routing for reservoir, its design formulas is as follows:

$\left\{\begin{array}{c}{W}_{t+\mathrm{\Δ}t}=W_0+\frac{\mathrm{\Δ}t}{6}(K_1+2K_2+2K_3+K_4)\\ K_1=I\left(t\right)-S\left(W\right(t\left)\right)\\ K_2=I(t+\frac{\mathrm{\Δ}t}{2})-S\left(W\right(t)+\frac{\mathrm{\Δ}t}{2}{K}_1)\\ K_3=I(t+\frac{\mathrm{\Δ}t}{2})-S\left(W\right(t)+\frac{\mathrm{\Δ}t}{2}{K}_2)\\ K_3=I(t+\mathrm{\Δ}t)-S\left(W\right(t)+{\mathrm{\ΔtK}}_3)\end{array}\right.---\left(8\right)$

Its calculation procedure is: first according to initial water reservoir level, calculates initial reservoir capacity and discharge flow, the setup algorithm period, initializes current time t=0, then solves K successively according to formula (8) ₁, K ₂, K ₃, K ₄, finally try to achieve calculation interval end storage capacity W _{t+ Δ t}, according to period Mo storage capacity, then look into water level storage-capacity curve and storage capacity discharge capacity curve, period Mo water level and period Mo vent flow, with the initial value that the result of calculation at this period end is subsequent period, circulation enters subsequent period and adjusts flood to calculate, until all calculation intervals calculate complete.

Due to the water level in storage routing for reservoir process, storage capacity, discharge capacity, relation between water process is very complicated, cannot be solved by analytic method, above-mentioned iteration trial and error procedure and numerical solution are two kinds of the most conventional at present calculation methods, but also there are some defects in above-mentioned two kinds of methods: each step of iterative method needs to carry out iteration tentative calculation, computational efficiency is low, calculate length consuming time, and computational accuracy is also limited by got allowable error, if allowable error obtains excessive, error of calculation then can be larger, if that gets is too small, tentative calculation number of times can be more, calculate meeting consuming time longer, in addition because computer exists rounding error, too small likely the causing that allowable error is got calculates situation about not restraining, fourth order Runge-Kutta numerical solution mathematically has Strict Proof, but as can be seen from formula (8), it often walks calculating needs inquiry four storage capacity earial drainage curves, carries out five step calculating, amount of calculation is comparatively large, and numerical solution has certain loss in computational accuracy.

Summary of the invention

Technical problem to be solved by this invention is: propose a kind of storage routing for reservoir method, solves storage routing for reservoir mode exists in conventional art low, the consuming time length of computational efficiency and the low problem of computational accuracy.

The present invention solves the problems of the technologies described above adopted scheme: a kind of storage routing for reservoir method, comprises the following steps:

A. read known Flood Routing through Reservoir relevant parameter, comprising: reservoir level storage-capacity curve, reservoir discharge capacity curve, reservoir reservoir inflow process, initial water reservoir level information;

B. according to reservoir level storage-capacity curve and reservoir discharge capacity curve, reservoir capacity discharge capacity relation curve is obtained;

C. setup algorithm period step delta t;

D. according to initial time reservoir level, look into reservoir level storage-capacity curve and reservoir discharge capacity curve, obtain reservoir initial storage and the initial discharge flow of reservoir, thus obtain water level Z at the beginning of the period of calculation interval ₀, storage capacity W ₀, discharge flow O ₀;

E. the average reservoir inflow in calculation interval step delta t;

F. according to storage capacity W at the beginning of the period ₀, checking storehouse holds discharge capacity curve at storage capacity W ₀slope k=S ' (the W at place ₀);

G. calculation interval end storage capacity W ₁;

H. calculation interval end reservoir discharge flow O ₁;

I. according to period Mo storage capacity W ₁, look into reservoir level storage-capacity curve, obtain period Mo reservoir level Z ₁=V ^-1(W ₁);

J. the period Mo water level Z this period calculated ₁, period Mo storage capacity W ₁with period Mo reservoir discharge flow O ₁as subsequent period period at the beginning of water level, storage capacity at the beginning of the period and discharge flow at the beginning of the period, then jump to step e, carry out the calculating of subsequent period, until all periods all calculate complete;

K. export the water level series of calculating, storage capacity is serial and discharge flow is serial.

Here " water level series " refer to water level at the beginning of the period of the day part calculated, period Mo water level; " storage capacity series " refer to storage capacity at the beginning of the period of the day part calculated, period Mo storage capacity; " discharge flow series " refer to reservoir discharge flow at the beginning of the period of the day part calculated, period Mo reservoir discharge flow.

Further, in steps A, read reservoir level storage-capacity curve W=V (Z), reservoir discharge capacity curve O=q (Z), reservoir reservoir inflow process I (t) and reservoir initial water level Z ₀.

Further, in step B, according to reservoir level storage-capacity curve and reservoir discharge capacity curve, the method obtaining reservoir capacity discharge capacity relation curve is:

O＝q(V ^-1(W))＝S(W)

Further, in step D, storage capacity W at the beginning of the period of calculation interval ₀=V (Z ₀); The initial earial drainage O of reservoir ₀=q (Z ₀).

Further, in step e, the average reservoir inflow according in following formula calculation interval step delta t:

$\stackrel{\‾}{I}=\frac{I(t+\mathrm{\Δ}t)+I\left(t\right)}{2}.$

Further, in step G, the method for calculation interval end storage capacity is:

$W_1=W_0+\frac{e^{k\mathrm{\Δ}t}-1}{{\mathrm{k\Δte}}^{k\mathrm{\Δ}t}}(\stackrel{\‾}{I}-O_0)\mathrm{\Δ}t$

Definition α=k Δ t is for adjusting big vast index, then above formula becomes:

$W_1=W_0+\frac{e^{\mathrm{\α}}-1}{{\mathrm{\αe}}^{\mathrm{\α}}}(\stackrel{\‾}{I}-O_0)\mathrm{\Δ}t$

Definition for adjusting big vast coefficient, above formula continues to be transformed to:

$W_1=W_0+\mathrm{\β}(\stackrel{\‾}{I}-O_0)\mathrm{\Δ}t.$

Further, in step H, calculation interval end reservoir discharge flow O ₁method be:

$O_1=\stackrel{\‾}{I}-\frac{1}{e^{k\mathrm{\Δ}t}}(\stackrel{\‾}{I}-O_0)=\mathrm{\α}\mathrm{\β}\stackrel{\‾}{I}+(1-\mathrm{\α}\mathrm{\β})O_0$

Further, in step I, calculation interval end reservoir level Z ₁method be:

Z ₁＝V ^-1(W ₁)。

Further, in step J, the calculating of subsequent period, initial water level, storage capacity and discharge flow adopt the period Mo result of calculation of this period, return step e cycle calculations, until all periods calculate complete.

The invention has the beneficial effects as follows: the present invention adopts integral analytic method directly to start with from the water balance differential equation, pass through theory deduction, directly Integration Solving is carried out to the water balance differential equation, obtain the universal equation of storage routing for reservoir, this method without the need to tentative calculation, calculate easy, solution efficiency is high, is easy to programming realization, and owing to being solve the direct integral of the water balance differential equation, computational accuracy is high.According to example test show, computational accuracy of the present invention and computational efficiency all than iterative method and runge kutta method high, compared with iterative method, computing time about 70.7% can be saved, compared with runge kutta method, can save computing time about 34.4%, computational efficiency promotes obviously.

Accompanying drawing explanation

Fig. 1 is storage routing for reservoir method flow diagram of the present invention.

Detailed description of the invention

In the present invention, storage routing for reservoir principle is as follows:

Known: reservoir reservoir inflow process I (t); Reservoir level storage-capacity curve W=V (Z); Reservoir discharge capacity curve O=q (Z); Reservoir initial water level Z ₀; Reservoir initial storage W ₀=V (Z ₀); The initial earial drainage O of reservoir ₀=q (Z ₀); Storage routing for reservoir period step-length is Δ t.

If adjusting flood to calculate initial time is t=0.

By the water balance differential equation:

$\frac{dW}{dt}=I\left(t\right)-O\left(t\right)---\left(9\right)$

According to reservoir level storage-capacity curve W=V (Z) and reservoir discharge capacity curve O=q (Z), because the two is all the monotonic increasing function about reservoir level Z, reservoir capacity discharge capacity relation therefore can be derived, that is:

O＝q(V ^-1(W))＝S(W)(10)

In calculation interval Δ t, because the general value of Δ t is less, therefore reservoir earial drainage storage capacity closes and ties up to calculation level W ₀near can be approximately straight line, if its slope is k=S ' (W ₀), then at W ₀near point, discharge flow O has:

O≈k(W-W ₀)+O ₀(11)

Formula (11) is substituted into formula (9):

$\frac{dW}{dt}+kW=I+{\mathrm{kW}}_0-O_0---\left(12\right)$

In calculation interval Δ t, reservoir inflow I gets period average, is calculated by following formula:

$\stackrel{\‾}{I}=\frac{I(t+\mathrm{\Δ}t)+I\left(t\right)}{2}---\left(13\right)$

Then formula (12) is written as:

$\frac{dW}{dt}+kW=\stackrel{\‾}{I}+{\mathrm{kW}}_0-O_0---\left(14\right)$

Because formula (14) the right is all known terms, according to the solution of First order linear non-homogeneous differential equation, integral and calculating is carried out to formula (14), can obtain:

$\begin{array}{c}W\left(t\right)=e^{-{\∫}_0^{t}kdt}\left({\∫}_0^t\right(\stackrel{\‾}{I}+{\mathrm{kW}}_0-O_0)e^{{\∫}_0^{t}kdt}dt+C)\\ =e^{-kt}\left({\∫}_0^t\right(\stackrel{\‾}{I}+{\mathrm{kW}}_0-O_0)e^{kt}dt+C)\\ =e^{-kt}(\stackrel{\‾}{I}+{\mathrm{kW}}_0-O_0)\frac{e^{kt}-1}{k}+e^{-kt}C\end{array}$

During by t=0, W=W ₀, then: constant term C=W ₀, can obtain thus:

$W\left(t\right)=W_0+\frac{e^{kt}-1}{{\mathrm{ke}}^{kt}}(\stackrel{\‾}{I}-O_0)---\left(15\right)$

For calculation interval end storage capacity:

$W_1=W_0+\frac{e^{k\mathrm{\Δ}t}-1}{{\mathrm{k\Δte}}^{k\mathrm{\Δ}t}}(\stackrel{\‾}{I}-O_0)\mathrm{\Δ}t---\left(16\right)$

Definition α=k Δ t is for adjusting big vast index, then formula (16) becomes:

$W_1=W_0+\frac{e^{\mathrm{\α}}-1}{{\mathrm{\αe}}^{\mathrm{\α}}}(\stackrel{\‾}{I}-O_0)\mathrm{\Δ}t---\left(17\right)$

Fixed for adjusting big vast coefficient, then formula (17) writing:

$W_1=W_0+\mathrm{\β}(\stackrel{\‾}{I}-O_0)\mathrm{\Δ}t---\left(18\right)$

Further, formula (15) is substituted into formula (11), and arrangement can obtain outbound discharge process:

$O\left(t\right)=\stackrel{\‾}{I}-\frac{1}{e^{kt}}(\stackrel{\‾}{I}-O_0)---\left(19\right)$

For calculation interval end storage outflow:

Further, asking the average storage outflow of reservoir in period Δ t, first to formula (16) integration in calculation interval Δ t, then divided by period Δ t, can obtain through deriving:

$\left\{\begin{array}{c}\stackrel{\‾}{O}=\frac{1}{\mathrm{\Δ}t}{\∫}_0^{\mathrm{\Δ}t}O\left(t\right)dt\\ =\frac{1}{\mathrm{\Δ}t}\[\stackrel{\‾}{I}\mathrm{\Δ}t+\frac{1}{k}(\stackrel{\‾}{I}-O_0)(\frac{1}{e^{k\mathrm{\Δ}t}}-1)\]\\ =\stackrel{\‾}{I}+\frac{1}{k\mathrm{\Δ}t}(\stackrel{\‾}{I}-O_0)(\frac{1}{e^{k\mathrm{\Δ}t}}-1)\\ =(1-\mathrm{\β})\stackrel{\‾}{I}+{\mathrm{\βO}}_0\end{array}\right.---\left(21\right)$

From above-mentioned derivation, in adjusting flood to calculate, average storage outflow is also not equal to and adopts in difference water balance formula (5), the whole story storage outflow arithmetic average it is more accurate that this also illustrates that algorithm proposed by the invention calculates.After calculating period Mo storage capacity according to formula (18), by reservoir level storage capacity relation W=V (Z), period Mo reservoir level can be tried to achieve.

By above-mentioned formula summarizing, the basic calculating formula of Flood Routing through Reservoir integral analytic method can be obtained:

$\left\{\begin{array}{c}{W}_1=W_0+\frac{e^{k\mathrm{\Δ}t}-1}{{\mathrm{k\Δte}}^{k\mathrm{\Δ}t}}({\stackrel{\‾}{I}}_1-O_0)\mathrm{\Δ}t\\ {\stackrel{\‾}{I}}_1=\frac{I_0+I_1}{2}\\ k=S^{\′}\left(W_0\right)\end{array}\right.---\left(22\right)$

The Computing Principle of a period step-length during above-mentioned algorithmic descriptions adjusts flood to calculate, the calculating of multiple step-length is exactly in fact the repeatedly use to above-mentioned formula.As shown in Figure 1, it comprises the following steps a complete storage routing for reservoir step:

1) reservoir level storage-capacity curve, reservoir discharge capacity curve, reservoir reservoir inflow process, initial water reservoir level information is read;

2) according to reservoir level storage-capacity curve and reservoir discharge capacity curve, reservoir capacity discharge capacity curve is generated;

3) setup algorithm period Δ t;

4) according to initial time reservoir level, look into water level storage-capacity curve and reservoir discharge capacity curve, obtain initial reservoir capacity and initial reservoir discharge flow, and in this, as calculation interval period at the beginning of water level Z ₀, storage capacity W ₀, discharge flow O ₀;

5) according to the average reservoir inflow of formula (13) calculation interval ;

6) according to W at the beginning of the period ₀, checking storehouse holds discharge capacity curve at storage capacity W ₀slope k=S ' (the W at place ₀);

7) according to formula (16) calculation interval end reservoir capacity W ₁;

8) according to formula (20) calculation interval end reservoir discharge flow O ₁;

9) according to period Mo W ₁, look into reservoir level storage-capacity curve, obtain period Mo reservoir level Z ₁=V ^-1(W ₁);

10) the last water level, storage capacity and the discharge flow that this period are calculated as subsequent period period at the beginning of water level, storage capacity and discharge flow, then turn the 5th step, carry out the calculating of subsequent period, until all periods all calculate complete;

11) last, export the water level series of calculating, storage capacity is serial and discharge flow is serial.

Embodiment: according to document (Chen Shouyu. the numerical solution of storage routing for reservoir and program [J] thereof. Journal of Hydraulic Engineering, 1980, (2): 44-49.) example that provides verifies.

Certain water surface area of reservoir is constant F=40 × 10 ⁶m ², flood releasing structure is without lock sluice weir, width B=100m, elevation of weir crest C=0m, the discharge coefficient m=0.4 on weir, flood into reservoir flow Q=2000m ³/ s.Weir formula adopts according to above-mentioned example situation, adopt elevation spacing to be that 0.5m calculates reservoir level ~ storage capacity ~ discharge flow relation curve, see the following form 1.

Table 1 example water level ~ storage capacity ~ discharge flow relation table

Adopt iterative method (allowable error gets 0.001m) respectively, method that runge kutta method and the present invention propose, adjust flood to calculate by programming to above-mentioned example, calculating achievement is in table 2.Compare for ease of each algorithm computational accuracy, adopt the material calculation consistent with document, take passages out the result of calculation of corresponding calculating time point, carry out precision comparison; Compare for ease of each algorithm computational efficiency, 10 times are computed repeatedly to often kind of algorithm, often all over above-mentioned example is computed repeatedly 10000 times, the computational efficiency of carrying out on average consuming time of then getting often time compares (note: this test is carried out on PC, and basic configuration is CPU:Intel Duo i53.30GHz; Internal memory: 8GBDDR31600MHz; Hard disk: 1TB7200 turns), the contrast situation of each algorithm computational accuracy and computational efficiency is in table 3.

Table 2 example calculating achievement table

The each algorithm computational accuracy of table 3 and computational efficiency comparison sheet

From above-mentioned example relatively, the method computational accuracy standard error that the present invention proposes is only 0.009m, higher than the 0.021m of runge kutta method and the 0.011m of iterative method, method exhibits excellent in computational accuracy that the present invention proposes is described.In computational efficiency, the method that the present invention proposes calculates consuming time than runge kutta method saving 34.4%, saves 70.7%, illustrate that the method that the present invention proposes has a distinct increment in computational efficiency than iterative method.

Claims (9)

1. a storage routing for reservoir method, is characterized in that, comprises the following steps:

A. read known Flood Routing through Reservoir relevant parameter, comprising: reservoir level storage-capacity curve, reservoir discharge capacity curve, reservoir reservoir inflow process, initial water reservoir level information;

B. according to reservoir level storage-capacity curve and reservoir discharge capacity curve, reservoir capacity discharge capacity relation curve is obtained;

C. setup algorithm period step delta t;

D. according to initial time reservoir level, look into reservoir level storage-capacity curve and reservoir discharge capacity curve, obtain reservoir initial storage and the initial discharge flow of reservoir, thus obtain water level Z at the beginning of the period of calculation interval ₀, storage capacity W ₀, discharge flow O ₀;

E. the average reservoir inflow in calculation interval step delta t;

F. according to storage capacity W at the beginning of the period ₀, checking storehouse holds discharge capacity curve at storage capacity W ₀slope k=S ' (the W at place ₀);

G. calculation interval end storage capacity W ₁;

H. calculation interval end reservoir discharge flow O ₁;

I. according to period Mo storage capacity W ₁, look into reservoir level storage-capacity curve, obtain period Mo reservoir level Z ₁=V ^-1(W ₁);

J. the period Mo water level Z this period calculated ₁, period Mo storage capacity W ₁with period Mo reservoir discharge flow O ₁as subsequent period period at the beginning of water level, storage capacity at the beginning of the period and discharge flow at the beginning of the period, then jump to step e, carry out the calculating of subsequent period, until all periods all calculate complete;

K. export the water level series of calculating, storage capacity is serial and discharge flow is serial.

2. a kind of storage routing for reservoir method as claimed in claim 1, it is characterized in that, in steps A, read reservoir level storage-capacity curve W=V (Z), reservoir discharge capacity curve O=q (Z), reservoir reservoir inflow process I (t) and reservoir initial water level Z ₀.

3. a kind of storage routing for reservoir method as claimed in claim 1, is characterized in that, in step B, according to reservoir level storage-capacity curve and reservoir discharge capacity curve, the method obtaining reservoir capacity discharge capacity relation curve is:

O＝q(V ^-1(W))＝S(W)。

4. a kind of storage routing for reservoir method as claimed in claim 1, is characterized in that, in step D, and storage capacity W at the beginning of the period of calculation interval ₀=V (Z ₀); The initial earial drainage O of reservoir ₀=q (Z ₀).

5. a kind of storage routing for reservoir method as claimed in claim 1, is characterized in that, in step e, and the average reservoir inflow according in following formula calculation interval step delta t:

$\stackrel{\‾}{I}=\frac{I(t+\mathrm{\Δ}t)+I\left(t\right)}{2}.$

6. a kind of storage routing for reservoir method as claimed in claim 1, is characterized in that, in step G, the method for calculation interval end storage capacity is:

$W_1=W_0+\frac{e^{k\mathrm{\Δ}t}-1}{{\mathrm{k\Δte}}^{k\mathrm{\Δ}t}}(\stackrel{\‾}{I}-O_0)\mathrm{\Δ}t$

Definition α=k Δ t is for adjusting big vast index, then above formula becomes:

$W_1=W_0+\frac{e^{\mathrm{\α}}-h}{{\mathrm{\αe}}^{\mathrm{\α}}}(\stackrel{\‾}{I}-O_0)\mathrm{\Δ}t$

Definition for adjusting big vast coefficient, above formula continues to be transformed to:

$W_1=W_0+\mathrm{\β}(\stackrel{\‾}{I}-O_0)\mathrm{\Δ}t.$

7. a kind of storage routing for reservoir method as claimed in claim 1, is characterized in that, in step H, and calculation interval end reservoir discharge flow O ₁method be:

$O_1=\stackrel{\‾}{I}-\frac{1}{e^{k\mathrm{\Δ}t}}(\stackrel{\‾}{I}-O_0)=\mathrm{\α}\mathrm{\β}\stackrel{\‾}{I}+(1-\mathrm{\α}\mathrm{\β})O_0.$

8. a kind of storage routing for reservoir method as claimed in claim 1, is characterized in that, in step I, and calculation interval end reservoir level Z ₁method be:

Z ₁＝V ^-1(W ₁)。

9. a kind of storage routing for reservoir method as claimed in claim 1, is characterized in that, in step J, the calculating of subsequent period, initial water level, storage capacity and discharge flow adopt the period Mo result of calculation of this period, return step e cycle calculations, until all periods calculate complete.