Fig. 3 Attenuation and lagging of the routed hydrograph. Inflow Q

Fig. 3 Attenuation and lagging of the routed hydrograph. Inflow Q
p
and outflow Q
O
hydrographs Q = f (time). The water entering the storage is hatched vertically; the
water leaving the storage is hatched horizontally
Горната фигура показва принципната схема на пълнене и празнене на язовира при навлизане на ВВ. При наличие на клапи, когато нивото стигне горния ръб на клапите те се отварят и изведнъж започва преливане с преливна височина равна на височината на клапите или на височината на частичното им сваляне. При сегментни затвори се предвижда само цялостното им повдигане. След това височината се

69 изменя в зависимост от количеството на притока. Не се предвижда повдигане на клапата по време на преливане, защото тази част на процеса, която е след достигането на преливния максимум, не е интересна за изчислението.
Този процес протича във времето, в което се изменят Q
приток
(t ), Q
прел
(t
) и
Q
осн.изп
(t). Q
приток
(t
) зависи от формата на ВВ. При триъгълна форма има два клона – на повдигане и спадане. При първия клон Q
приток
(t
)= К
1
..t + Q
осн
, а при спад Q
приток
(t )=
Q
мах-
-
К
2
.(t-t възх
) + Q
осн
Q
осн е нормалният приток в язовира преди ВВ. Коефициентите К
1
.(t
) и К
2
.(t
) са ъгловите коефициенти на възходящия и низходящия клон на триъгълника – К
1
.(t)=
Q
мах-
/ Т
възх.
и К
2
.(t)=Q
мах-
/ Т
низх
. Т
възх и Т
низх са времената на възход и спад на вълната. Те се задават като части от цялото времетраене на вълната.
Основното уравнение, което описва този процес е това за баланса между притока и разхода:
(Q
приток
(t ),

- Q
разход
(t ))dt = dV.
Това е обикновено диференциално уравнение с начални условия: t
о
=0, V=V
o
, Q
оприток
=0 или Q
приток
=
Q
осн
; Q
разход
(0 )= Q
о
, Q
прел
(0)=0 (6) h
о
=0 или h о
= h k
Z=Z
о
- нивото в езерото в началния момент.
Уравнението се решава в крайни разлики в явна схема, за да не са нужни итерации:
(Q
приток
(t n
)
– Q
разход
(t n
))* ∆t = ∆V
n
(7)
∆t=n часа, t n +1=
t n +
∆t ; ∆Z(t n+1
)= Z(V
n+1
)- Z(V
n
), V
n+1
= V
n
+∆V
n
Q
разход
(t n
) =
έmB
g
2
h(t n
)
3/2
+ Q
максекспл
, + Q
осн.изп
(t n
) h(t n
)= Z(V
n
)- Z
max
По описания изчислителен метод е съставен алгоритъм и са програмирани електронни таблици в Excell и са изчислени два примера. Те са описани в
Приложение 11. Таблиците могат да бъдат използвани за всякакви язовирни стени.