"$%":! " # 1 2, e $. %. &!" ' ( (rkm ) (rkm ), )* (rkm 0+000) + (rkm ),, )* (rkm 0+000) " (rkm ), )* $ ' (,).

Size: px

Start display at page:

Download ""$%":! " # 1 2, e $. %. &!" ' ( (rkm ) (rkm ), )* (rkm 0+000) + (rkm ),, )* (rkm 0+000) " (rkm ), )* $ ' (,)."

Della Allison
5 years ago
Views:

1 - 1,! 1,2, "! 1, # 3 1 -,, 2,, 3 -,, "$%":! " # 1 2, e $ % &!" ' ( (rkm ) (rkm ), )* (rkm 0+000) + (rkm ),, )* (rkm 0+000) " (rkm ), )* $ ' (,)- &, &,,,", ) (/)!" # 1, # 2, ', ) 0!,!" # 1,!" # 1 # 2, *!" # 2 1,, ) & HEC-RAS :, 1 2, HEC-RAS MATHEMATICAL MODEL OF UNSTEDY FLOW IN IRON GATES RESERVOIRES - BASIS FOR ENERGY PRODUCTION ANALYSIS Z Kova&evi' 1, M Spasojevi' 1,2, A Ostoji' 1, S Mili' 3 1 Energoprojekt-Hidroinženjering, Belgrade, Serbia 2 Faculty of Civil Engineering Subotica, University of Novi Sad, Serbia 3 Energoprojekt-Entel, Belgrade, Serbia Abstract: This paper presents the one-dimensional model of unsteady flow in HPP Iron Gates 1 and 2 reservoirs, caused by the combination of both HPP regimes and natural waves entering reservoirs The model results are used as basis for energy production analysis The area modeled includes the Danube River between the gauging station Lom (rkm ) and the city of Novi Sad (rkm ), the Sava River between its mouth (rkm 0+000) and the city of Sabac (rkm ), the Tisa River between its mouth (rkm 0+000) and the city of Novi Becej (rkm ), as well as the mouth areas of smaller tributaries of the Danube (Tamis, Velika Morava, Mlava, Nera, Pek, Porecka, and Timok) and the Sava (Kolubara) The model includes the HPP Iron Gates 1 structure, as well as both HPP Iron Gates 2 structures, one on the Danube River, and the other on the Gogos side channel The model s domain is split into three subdomains, the first one covering area upstream from the HPP Iron Gates 1, the second one covering area between HPP Iron Gates 1 an HPP Iron Gates 2, and the third one covering area downstram from the HPP Iron Gates 2 The model s formulation, calibration, and verification are based on existing bathymetric, discharge, and free-surface elevation data The model is based on the HEC-RAS program Key words: 1-D unsteady-flow modeling, HPP Iron Gates 1 and 2 reservoirs, HEC-RAS program 1

2 1 # 1 $ (rkm ) ' (rkm ),, )* (rkm 0+000) " (rkm ), )* (rkm 0+000) + (rkm ) # 2 2 $, m ' (rkm ~ ), ) (~rkm ), $ #1 # 1 2!",,! " # 1 $, $ $! " # 2, #1 # 2, $! " # 2 $, $, $ " 0! ", " " # 1 2,!" $ 0 ) 30, # 2, * $- ", ) ) 2)!,! " # 1 2, " (" $ )! " )!" ), ) "! " (St Venant-ove ", ", ) $ ", ), ) * $ " $" " )$ $, " $ ( [1] [2]) 0 )$ ",! " [1] [2] $!, " " ' $ # 1 ( ',,), " " "$ ", "" " $$ )*, 3$!* )$ ",! " $ # 1 2 ( 1) [3] [4], ) "$ ()"$ )$ ", * $ ), * [1] [2] 2

3 #% 1 "( ) $* 1 2 +# " (*# 2 ) / $ $ $ # 1 # , " " 1998, $" " )*, # (), /-,, # ( )) )*, : )*, ,30/70,40 ; # 1 -, ,59 ; /-,, ,00 (42,50 ); # 2 -, (41,25/41,40 ); )*,, , +/- 30 ; % $$ $" $ $" ",!$ $ $,,!, * " : 1 - # 1; 2 - # 2 ; 3 - ' # 2 3,,) - 3 $",!,, " # 1 # 2 '!, % : (1) "$ $ 3

4 ; (2) ), "$ ; (3) ) 31 +#% %$" /$01"#1!$1 $ # 1-1 (',,,,- &,",,&,,,),/ ), # 2-2 ' # 2-3, ) $ $ )* " / $ )* 1:25000,, * 4 1 ) "$ $ 3 " 3 1 )* 315 "$ 2005/2006 & # "$ 2005/2006 & 3 - # 2 82 "$ / ) )* "$ 2005/06 # 2 ), -)* 3, -)* (-() ",! " )*,&-! # 2, 2 # -/, ) 30 0! ", - (, ) * $ 33 *2 1( $"2( * 4 $ # 1 2, %" $-) (%&3) 0 "$,! # 1 # 2 : $,,, 3 " $ , 1 2, $ $ "$ 4

5 +# 1 #3 *2 0 4$($5 6$11/ # 6$ # $% 0$ *"% " r%m *"% r%m m*# 0*1-3 1 ) Qz m *1-3 1 ) $* 2 Zh/Qh ) $* 2- ~ *1-2 1 ) $* 2 Zh/Qh ) $* 2- ~ *1-2 1 ) $* 1 Zh/Qh ) $* 1- ~ *1-1 1 ) $* 1 Zh/Qh ) $* 1- ~ *1-1 1 ) Qd * *1-1 ) Qd *1-1 ) Qd 1" *1-1 #+$ ) Qd -# -$* *1-1 m8 ) Qd 8m *1-1 $ ) Qd " *1-1 # ) Qd -$" *1-1 $ ) Qd *1-1 % ) Qd *1-1 $% ) Qd # # 2 #3 *2 0 %#+$2! $4%2! (*# + $% 0$ *"% " $%m *"% 1!+#3 $4# m*# 1 ) Zd -! ) Zd "# ) Zd ) Zd ) Zd -6$* ) Zd * ) Zd #1%m ) Zd m ) Zd ) $* 1 Zh ) Qd m*$ ) Zd m*$ ) $* 1 Zh $ (m) ) $* 1 Zh # ) $* 1 Zh +$ ) $* 1 Zh # ) $* 1 Zh/Qh $* 1- ~ ~ ) $* 1 Zh/Qh $* 1- ~ ~ ) $* 2 Zh #* ) $* 2 Zh -$0 #1% ) $* 2 Zh #" ) $* 2 Zh/Qh $* 2- ~ ) $* 2 Zh/Qh $* 2- ~ ) $* 2 Zh 8 m% ) Zd "" ) Zd #4" ) Zd $ ) $ /*$*1(%6 ("("%6 (*# 0 (#2! "5 "( 3 " # 1 2 HEC-RAS 41, (USA, COE - Hydrologic Engineering Center) [6, 7, 8 9] 5

6 42 0+$ %$%"$ "1/ $"2! ($*) 0 %#+$2! (*# 4 ) ) $ # 1 # 2, $,!$ " " )! " # 1 2,, "$, / $ # 1 2, ( ) $, # 2, %&3 43 *2 0 4$($5 6$11/ #, %#+$2! $4%2! (*# $ "$ $, ($% $("$ $("$ %#+$2! (*#! " # 1-1, $, " $, " 5-15, $ 4 / $ )* Manning- n ( Strickler- S tr ) $ 3 1 n = m -5 s (S tr = 23-50), " 3 "!, ) " ", "! ), ", )*,, )* ) * $, * $ $ $ 4 $ )* $ $, (, 2)! " # 2 $ # 1, $ # 2 ) 3 2 " n = 0025 m -5 s (S tr = 38) 3 3, # 2, n = m -5 s, n = 0026 m -5 s (S tr = 38) 6

7 5 ) # " #1!%6 (*# 8(" % $%0 (*# 1 ($3 *"% (! " # 1-1 $ ' # 1, )* +, )* "!" - & (",,, ',', &,,) /) 1 ) "$ $ 52 $%0 %(1"$ $0#"" %#+$2! (*# / $ $, 2, ' & ( rkm ), "!$ # " ( 2)! " $, $, "!$ ( ), $ $ " KALIBRACIJA - VS DONJI MILANOVAC - Qsr = 5000 m 3 /s KALIBRACIJA - VS DONJI MILANOVAC - Qsr = m 3 /s KOTA NIVOA (mnm) MERENO MODEL KOTA NIVOA (mnm) MERENO MODEL VREME (h) VREME (h) #% 2 #+$2! (*# Q= m 3 /s: 1 * (($1/(*#- 1 $*1 "), %2! : 5 #12 53 $%0 %(1"$ $0#"" $4%2! (*# 3 # 1-1 ) " "! (, = 100 )), )! $) -$" 0,!* (-), ) * ", " ",!", $,! - " 2006 ), 7

8 ,, "$ ' &,, $ " # "! #% 3 $4%2! (*# : 1 * (($1/(*#- 1 $*1 "), %2! : 5 #12 3 " " $ 2006 ", " ", 4 ( $ $ ), 41,! $ $ #% 4 $4%2! (*# : 1 * ($0*# ($1-(*#; 1 $*1 "), "% 6*1 0 1 " * $"2! %2! : 5 #12 5 $ " $ (%=0,999), " " " y = x (- ) " 91 % 5 cm 8

9 %! "#$$$$ % REZIDUALI (MODEL-MERENO) (m) KRIVA TRAJANJA REZIDUALA (MODEL-MERENO) TRAJNOST (%) #% 5 $4%2! (*# : %$#2! ($1/ 1 1 *+!1/ (#2!( 1 (*# ( 1 $*1 "), %$ "$!5 $0*#, %2! : 5 #12 3 $ " $,, $ "$ +# 3 11 "" "% $("$ $*1 " 1 ((*# ($1) O Mo Me / RMS RMSE Nash-Sutcliffe NSC R2 = 1- RMSE/RMS 0!,! m R = ) # " #1!%6 (*# 8(" % $%0 (*# 1 ($3 *"% (! " # 2-2 $ ' # 1, ('-) # 1 # 2 ', ) 1 ) "$ $ 62 $%0 %(1"$ $0#"" %#+$2! (#2! / $ $, 2, / ( ), $,!$ # 2 " ( 6) 9

10 4200 KALIBRACIJA - VS KLADOVO - Qsr = 5000 m 3 /s KALIBRACIJA - VS KLADOVO - Qsr = m 3 /s KOTA NIVOA (mnm) MERENO 4020 MODEL VREME (h) KOTA NIVOA (mnm) MERENO 4180 MODEL VREME (h) #% 6 #+$2! (*# Q= m 3 /s: 1 * (($1/(*#- 1 $*1 "), %2! : #* 63 $%0 %(1"$ $0#"" $4%2! (*# 3 # 2-2 ) " )!"! - " 2007 ) # 2 "! 53, $! (rkm ), ( 7) $ &! #% 7 $4%2! (*# : 1 * (($1/(*#- 1 $*1 "), %2! : -$0 #1% 10

11 "#'( & *)#"" & "#'" "#"( "#"" "#"( "#'" % *'#(" *'#"" *"#(" *"#"" $#(" '#""")! "#$$ "#'( " )""" *""" +""",""" '"""" ')""" '*""" '+""" $#"" $#"" $#(" *"#"" *"#(" *'#"" *'#(" *)#"" % #% 8 $4%2! (*# : 1 * ($0*# ($1-(*#) 0 1 " * $"2! %$#2! ($1/ 1 1 *+!1/ (#2!(, %2! : -$0 #1% 8 " $ (%=0,993), " " " y = x (-)!" 90 % 5 cm 4 $ " $,, $ "$ +# 3 11 "" "% $("$ $*1 " 1 ((*# ($1) O Mo Me / RMS RMSE Nash-Sutcliffe NSC R2 = 1- RMSE/RMS 0!,! m R = ) # " #1!%6 (*# 8(" % $%0 (*# 1 ($3 *"% " ", )*, ', %, +/- 30 cm 0 " / %2 % $ $ $ # 1 # , " " # 1 FE Portile De Fier 1998 " # 2-3, *, $ $" $ " ),, " ' )*, ' 2) 3!" * ), -)*, ) " $ $ '#, O ) )*, ', ", " 11

12 # 2, " ) 3 $ * km )*, * " (! " ', # 2, $ ' ('-) # 2 ', ('-) ) ), ( 72 $%0 %(1"$ $0#"" %#+$2! (#2! 3! 3, # 2, $ $, 1 2, $! " )* $ / $ $!,, ( 9) KALIBRACIJA - VS TIMOK - Qsr = 3400 m 3 /s KALIBRACIJA - VS TIMOK - Qsr = 7500 m 3 /s KOTA NIVOA (mnm) MERENO MODEL KOTA NIVOA (mnm) MERENO MODEL VREME (h) VREME (h) #% 9 #+$2! (*# Q= m3/s: 1 * (($1/(*#- 1 $*1 "), %2! : (%-8 73 $%0 %(1"$ $0#"" $4%2! (*# 3 # 2 ( 3 ) ) " , $, " 0 ) # 1 # 2 (!"* )),, -)*, $ "$ 2006, 3 - # 2, ) ,!* (-), ) * ", " ",!", $,! " - " ) " # 2, 12

13 "! 53, $!, -)* (rkm ), (#"" %- *#"" #"" )#"" '#"" "#"" )$#"" ),#"" /0"+ /0"+ 1"2 13"2 43"2 1"2 5"2 56"2 57"2 %/4"2 0"2 89"2 /0"2 #% 10 $4%2! (*# * *2 2007: 1 * (($1/(*#- 1 $*1 "), %2! : (%-8 "#)" %- %- % "#'( "#'" "#"( "#"" "#"( "#'" "#'( % (#"" *#"" #"" )#"" '#"" "#"" )$#"" "#$$$$! "#$$,( "#)" " '""" )""" """ *""" (""" +""" 2""",""" $""" '"""" ),#"" ),#"" )$#"" "#"" '#"" )#"" #"" *#"" (#"" % #% 11 $4%2! (*# : 1 * ($0*# ($1-(*#) 0 1 " * $"2! %$#2! ($1/ 1 1 *+!1/ (#2!(, %2! : (%-8 11 " $ (%=0,998), " " " y = x (-)!" 95 % 10 cm 4 ) " $ "$ 8 4!" $" " ) $ $" " $, $ $ *$ 13

14 ,,! $ 0 : #, ) ", ", $, * " ", $,!", ( " ), " "$! " # 1 2, ) $ " 9 : [1] "* # 2", 4, III, 1: "" " " ", -,, 1982 [2] *, &, " " $ $ #", 3 X4- " " 6 (/ III: "/)* *$ )"),, 20-22, 1984 [3] # 1 ", -,, 2011 [4] " $- * # 2", -,, 2012 [5] Cunge, JA, Holly, FM, Jr And Verwey, A, Practical Aspexcts Of Computational River Hydeaulics, Pitman Publising Ltd, London, 1980 [6] HEC-RAS (River Analysis System), Version 41,Application Guide, USCOE, USA, 2010 [7] HEC-RAS (River Analysis System), Version 41,Hydraulic Reference Manual, USCOE, USA, 2010 [8] HEC-RAS (River Analysis System), Version 41,Users Manual, USCOE, USA, 2010 [9] HEC-GeoRAS (GIS Tools for HEC-RAS), Version 4, Users Manual, USCOE, USA,

Modeling of One-Dimensional Unsteady Open Channel Flows in Interaction with Reservoirs, Dams and Hydropower Plant Objects

Journal of the Serbian Society for Computational Mechanics / Vol. 3 / No. 1, 2009 / pp. 154-181 UDC: 626.861:621.311.21 Modeling of One-Dimensional Unsteady Open Channel Flows in Interaction with Reservoirs,