Scientific Bulletin of Politehnica University of Timişoara
Transactions on HYDROTECHNICS
Volume 60(74), Issue 1, 2015 Advance hydraulic modeling of irrigation channel CP10, in irrigation system Teba – Timisat, Romania, Timis county Robert F. Beilicci1 Erika Beilicci1
Abstract – Study case is situated in irrigation system perimeter system: Cruceni Foieni, Ionel, Otelec, Teba - Timisat. To solve theoretical problems of Ivanda, Sinmartinu Sarbesc, Sinmartinu Maghiar, movement of water in the channel CP10, it requires Sînmihaiu Roman partially Dinias. modeling of water flow in this case. Numerical modeling Irrigation system Teba - Timisat has a total area of was performed using the program MIKE11. Advanced 28,063 ha, is located in interfluves TIMIŞ - BEGA computational modules are included for description of flow over hydraulic structures, including possibilities to downstream of Timisoara, the basin of the rivers describe structure operation. The input data are: area Timis and Bega, is limited to the north of Bega plan with location of cross sections; cross sections navigable channel to the west of Serbia, south Timis topographical data and roughness of river bed; flood river and drainage systems Rudna - Giulvăz and discharge hydrograph. After simulation with MIKE 11 Caraci and east system Sag - Topolovăţ. result the water level in each cross sections.
CS60
CS58 Q=5,5mc/s CS60a
CS57 CS62 CS64 Sanmihaiu CS66
CI2/3 CS1 CS67 CP18
CT1 CI1 CP5/1 CI2/2
CS61 CS68
CS70 CP16 CI2/1
CS71a CS75
CS59 CS73 CP1/1a
CS55 Roman CI1 CP10
CP1/1b CS1 CT2 CI2/3a CT10 CT47b CS1c CS16 CS67d CP16/1 CS99 Key words: hydraulic modeling, irrigation channel CS63 CS65 CS74
CS17 CS63 CS69 CS71 CS72 CT7 CP1/4 CT5/1 CS67b
CS6d CS54
CP1 CT2c CS1b CT9
CP1/2 CT1a CI2/A
CT45 CI1/1 CT1b CT1c CS2 CT1d CS1f CS3 CS13 CP5/1 CS67a CS90b
CP1/5 CT1e CT9a CS81 CS90a
CT44 CP1/3 CT1f CT5a CT12 CS53
CT3 CT4 CT5 CP16/3
CS6e CS90
CT1g CS12
CT4a CS94 CP19/1 CP1/1 CS15 CT2a CT2b
CS13 CP5/1
CI1 CS81 CT8 CP17/1 CT12
CT11/2 CS80 CP14/4
BEGA CT11 CP14/2 CS79 CP14/1
CT12a/1 CS13d CS6a CS12
CS1a CS5 CP19 CS6c
CT9
CT8
CP14/3
CS5c CP13/1 Canal CS5b CP11 CT11
CT7 CP5/1
CS78 CT10 CT12 CP2/1 CS17a1 CS4
CI3 CS17a CT56
CP17 CS14b CS77 CP14/3 CS100 CS8 CT12a CT52
CI2a CT54
CT11/1 CS17a2 CS76
CS16 CS91 CS97a
CS13a CS14 CS101 CS6b CS14a CP6 CP14 CS5d
CI2 CS89a
CP2/5 CS19 CS79a CPE DINIAS
CP2/2 CP2/3 CS17 CS6 CT6 CT6/2
CS8 Sanmartinu-Maghiar CS13a CT6a CP7 CS96 CT13b CS5a CT63/1
CT13 CT6/1
CP2/1a CS78a CS10 CP10
CT6/0 CT55
CS9 CP7/2 CS77a CS88
CP2/4 CT12b CT53
CP2/5 CS18 CT59
CS8 CT15
CS9 CS7 CS15 CS12d CS7 CP2/1b CT14a/1 CP12 CS92
CI1a CS9 CP5/1 CS89
CS31a CT26 CS87a CP1/1
CP2/6 CT63 CS97b CT14a
Otelec CS12b
CT14/1
CS19 CS86a
CT58 CP7/1
CP1/1b
CP4/10 CS14
CT14 CP4/11
CP1/1a
CP11/1 CT60 CP4/9 CT27
CS38 CT57 CS92
CS51
CS22
CS86 CP10/1
CT23 CP2/6a CP4/8
CS38a CS85
CT11b CS48a
CS10b CP7/3 CT22
CS1 CS62 CT11b CT64 CP4/6 CS52
CP1/3 CT8/1 CT11b/1 CS11a
CT11c CP9 CP15 CP4/7 AMONTE CS6 CT26 CS26 CP5a CS47a
CS29 CS47 CP15 CS46
CS6 CS30 CS31
CS21 CS1c CS43h CS10c CT11d CS10 CS51 CS10a CP2/7 CP2/8 CP1/4 CS11 CS27 CS1d CT11d CT21 CP2/9
CP3 CL CP7-CP2 CS28
CT25a/3 CP4/5 CS48 CT11e CS43g CT65
CS11 CT11f CT13 CT7/1 CS20a
CT5 CT25b CS26 CP5 CS20 CP4/4a
CS1
CS40 CS1 Dinias
CS12 CS25 CS5 CT19 CS18
CS1/c CS49 CT8
CPA CT25a CS50
CT18
CT7
CS1e CS24
CS17a3 CS45 CT40
CS10 CT23g/1 CT11e CP2/11
CS25 CS43e CT25a/2 CT14 CP2/4 CS1c CP4/4
CT11g/1
CT13 CS38
CT17 CT11g CT2
CS29 CT1 CS1 CT3b
CT23g CS17a4
CS1b CS28 CS27a CT15 CT18a CS17b CT25a/1 CPA/2 CT29 CT13/1 CS38 CP2/12 CT3c/1
CP1 CS1f CT5a CS45
CT23f CT3c CPA/1 CS11a CT23e CT16 CT39/1
CT13 CS12b CS27
CP2/2 CS4 CS50 CT3d CPA/11 CT28 CT42
CS11a CP2/7
CT11h CT39
CP7 CS3 CS12 CS25 CS48a CP2/15a CS1b CS28 CT40a CP2/14a CP2/15 CT3f I. INTRODUCTION CT31 CT3e CS65 CS67 CT17 CS43c CS48 CT3a CP4/2 CS45 CS2 CP1/8 CP3 CP2/14 CT3 CS1 CT23d CS45 CS1a
CT25/2 CS9b
CT23b CT28 CP8/3 CPE Otelec CPE CS46 CT30 CS28 CT33
CP2/6 CT3g CS9a CT27 CT43
CS13 CT16/2 CS23 CT25/1 CS11
CS1a CT18
CS43b
CT12/2 CT12/6 CP1/9
CS22c CT41 CT12/4 CS9 CP2/4 CS36 CP8/2 CP4/1 CT33a/1 CT18/1 CS38
CT25 CS66
CS23 CP3/1 CT30/1
CS10
CT16b CP8/3 CT19
CS22b CP4 CT18/2
CP1/9a CT24 CS46 CT12 1081CT16 CS2a CT34
CP8/6 CS27 CT17a
CT18 CS1 CL CP3-CP4 CT39 CS12a CP4/14 CS30 CP1/10 CT28a CS68c CT33a CS22 CP8/6
CS22e CS6a CS27 CT20a CT12/1 CS45 CS19
CP2/5 CP8/1 CS69 CP4 CS28 CS2 CS20a CS22 CT16a
CS33 CP6 CT21 CT12/3
CS12 CT37
CT36 CS11
CP4/15 CP8 CT12/5 CT37a/1 CT20 CS17
CS64 CT22 CS4 CS43a
CP4/11 CT23a CS34b
CS8 CP4/9 CS20 CS5a
CS14 CS34b CS62b CS33 CS3 CT37a CS31 CS22d CT37a
CS17a CS57 CS7 CT21a
CT12/7 CP4/1 CP5 CS33b CPE OTELEC CS58b
CP7/1 CS58a CT4
CT1 CS7 CS62a
CT1b CS6 CPE OTELEC
CS15 CT24a CT37b CS59
CPA CT1c CS21 CP7/1b CS33a
CS17b CP7/3
CT37a CS34
CP2 CS32j CS6 CP4/14 CS17c CT15 CT23h CS32h CS2
CS8 CP5/2 CS58 CT38 CP3 CS5
CP4/11 CS32f
CS63 CT22/1 CT2 CP2/7 CT3 CP7/1a CS58c
CP4/3 CT39
CT37
CP4/5 CS62 CS61 CS6b
CS32k CP6/2 CP6/4 CS7 CPE OtelecCP4/8 CS16 CP5/3a Sanmartinu
CP6/1 CP4/18
CT9/1 CT4 CT5 CP4/18 CP4/6
CT9/3 CS16 CP5/4 CT12/7 CS15a CP4/10 V.Temesit CP4/13 Sarbesc
CS15 CP4/4 CS32c CP5/3
CT10/2 CT6 CP4/6 CP4/4 CP4/12 CS32b
CPA CT7
CS17
CT10/4 CS32a CS32d
CT8 CP13/1 CS7
CT10 CS2 CS18 CS14 CS31
CS8 CP15 CP5/5a CS32
CT7 CP4/2 CS8
CPA CT5 CS7 CS32e CT10/1 CP4/2 CS6 CP15a CS31a
CT10/5 CT10/3 CS12
CS5 CP13 CP15b CS12 CS6b
CT2/2
CS7 CS6c CS6a
Arrangements cuprinds in this area are: Teba - CP3/1 CS2 CL CP8-CP2 CT11b CS4 CT6 CT11/1 CT2 Ionel CS11 CS18 CT1
CT10 CT8
CT11 CS17 CP12/2 CT17 CT3
CT11/3 CS8 CT11/4 CT11 CS5 CS15 CS16 CS17a CS10
CT16 CS17
CT11c CT6
CS9 CS25 CT1/1 CS14
CT6 CT11a CT15
CS16 CT13
CS22 CS24 CS11 CS13
CT12 CT7 CT14 CS11
CT14 CS7 CT13
CS15 CS10 CS14 CT6 CS21 CT6/1 CS5 CT10 CT11
CP8 CS6
Timisat, Rudna - Giulvăz, Bociar and Caraci. They CS14 CT4 CP12 CP12/1
CS8 CT24b CS13 CS11 CS12 CS19 CT3 CS18 CT4 CP10 CS18 CS3 CS9 CT5 CS1 CT5a/1 CT2
CS4 CT5 CT24 CT24 CT5a CS9
CS2 CS18a
CT5b/1 CT2 CS1 CT3 V.Teba
CT5b CP10 CS16
CT5/1 CT1 CT25 CP11/3 CS18 SPR CS15
CS8 TEMESIT CS20 Ivanda CT2
CS19 CS7 CP12/1a CS13a
CS13 CT3 CS40
CT46/1 CT45 CS14 CP11 CS5 CT4 CP11/2a arranged a total surface of 43,335 ha. But because of CS2 CS10 CT11/2 CS4 V.Teba CS9a CT49 CL CP7-CP8 CP11/1 CS3
CS3 CP11 CS1 CT1a CP11/2 CS9b CS4 CP6 CS38 CS4a
CS1 CS5 CP9/1 CS1b CT 15 (CPE Cruceni) CS1c
CT13 CT11/1 CS35a
CT14a CT14 CS36 CT32 CT12 CL CP9-CP11 CS2 CS29 CS10 CS38 v.Timisat CS1a CT5/1
CS43 CT44/1 CS15a
CT11 CT44 CT1b
CT23 CT1 CT33 CS1 CT12 CT34 CS21 SPR CS6b CI5 CT15 CT20
CT47 CP6/1 OTELEC CS8 CT11 CT10 CT26 CT22 CT9 CS37 SUD CT10/1 CT5/2 CS23 CT20a CS6c CS14 CI4 CT10 CT8
CS12 CS15a
CT21 CT8 CT5/3 CS8 drainage arrangements can be made independently on CP11 CI3 CS21 CS16 CL CP9-CP11 CT48 CS13 CT16 CT9 CT5/4 CP6 CT7 CT7
CT4/1 CS28
CS19a CT5/5 CT13 CS14a CS14b CT12 CS15 CS17 CT20 CS16 CS17a CT7 CS13 CT8 CT12 CS13
CP9 CT15a CT15 CT5/6 CI2 CS35 CS7 CT19 CT9 CS25 CS10 CS34b CT16a CS12
CT6 CS14 CS9 CT11 CT4 CT5 CS15 CT16 CT10 CT37 CP2 CS11 Giulvaz CT11 CT17a CS22 CT6 CT13a CT18 CS32 CS34 CT40
CS20
CS33 CT14 CT17 CS6 CS5 CT13 CS19 CP5 CP4 CS18 CT16 -CP6 CS28 CT38 CS18 CS1 CP2 CT3
CS12 CS4 CT39 CS2 CS4 CT8 CS16 CL SPR CT2 CP2 CS34a each fitting generally most suggestive aplor planning CT29 SANT CS17 CS16 IN COT CS15
CT25 CT41 CT6/4 CT6
CT6a CT1 CS14
CS13 CS13a CT25a CT6/2 CT31 CS3
CS8 CS31a CS31b CS13 CP4 CP2/2 CS8a CS30 CT30 CS31 CS11a CT4 CP2 CS11 CP2/1 Foeni CS7 Cruceni) CS10 CS6 CS9 (CPE CT5 CS9a CI2 CS9 CP2 CS9 for evacuation in state border surrounding area is CP4 CS5 CS4 CT1/1 v.Timisat CT1
CT1 CS1
CS1 CP3-CP4 CS3
CS3 CT1 CL SPR CT1a BICA CP2 CT3 CS3
CP3 CS2 CP3 SPR MLACA
CT1 CT2 CS2
CS4
CT2/1 CT5/2 CS1 CT2 CT5 CT5
CT2 arranging Teba - Timisat. Cruceni CS4
CS1c CP1/7
CS1b CS9
CP1/6 CS10
CL CS7 CP1-CPE(R-G) CP14 CP1/5 CT4 CS2 CS6 CS8
CS11 CS12
CT4a CT8 CS3 CS5
CP1/3 CP1 CT5/1
The purpose of the performed works is remove legator CT5
CT1 CT5/2 Canal CT1a
CT6 CP1/2
CT2 CS1 CT1b
CT3 CS4 CT3 CS1a CS1
CS12 SP CRUCENI CT7/1
CP1/1
CT7 excess moisture from the surface, and soil in the r.TIMIS active profile on areas with specific requirements in order to adjust its moisture, creating appropriate Figure 1. Irrigation system Teba – Timisat conditions for making timely farming operations, with positive effects on production. The level of groundwater in the area is influenced by To achieve this goal were achieved partly works and several factors, namely: drainage ditches in the area of 28 063 ha in Teba- - Change the channel levels Bega and Timis river; Timisat drainage system, as follows: - Precipitation on surfaces located in the upstream - Main collector drain arrangement (regulation side of the perimeter of the system (planning); Timisat, Teba, Temesit) and connection of large - Rainfall in the period from October to May each channels or valleys extending the length of 65 km; year. - Completion of drainage works by thickening and Hydrogeological zoning made prior to execution, deepening the existing channel network in an area of based on observations made during the period 1973 - 22213 ha; 1984 drilling IMH network and IEELIF existence - New network of drainage channels on an area of field highlighted areas related to groundwater depths 5,850 ha; intervals as follows: - Closed horizontal drainage works cross type on an - 0.5-1.0 m: 974 ha area of 285 ha, located in a high groundwater level - 1.0-1.5 m 4643 ha trough north of town Ionel. - 1.5-2.0 m: 7406 ha From the administrative point of view, the drainage - 2.0-3.0 m: 9873 ha system Teba - Timisat is located in the western part - 3926 m over 3 ha of Timis County, including related land cadastral Following regionalizing technique, maximum levels territories Foieni, Giulvăz, Peciu Nou Uivar, recorded in the same period, it appeared that the Sanmihaiu Roman and Sag in the localities falling 1 Politehnica University Timisoara, Faculty of Civil Engineering, Department of Hydrotecnnical Engineering, Timisoara, George Enescu Str., no. 1/A, Timisoara, Romania, [email protected] 9 perimeter, 78% of the area in the spring has deep - Duration excess moisture in dry years 5 months; phreatic level at 0 - 1.00 m. - Duration excess moisture in wet years: 7-8 months; Groundwater flow direction is oriented NE - SW with General natural slope of the land is northeast to average gradient of 0.3%. southwest, being a 0.2-3%. Excess moisture in the arrangement Thebes - Timis Emissaries vacating the main river waters are appear differently on surfaces varying depending on navigable Timis and Bega channel. Timis River the orographic, soil and local hydrogeological comes gravitational allow downloads of inland waters in from rainfall and groundwater perimeter fueled partly about 70% of days of average, and for about 131 days navigable channel Bega. during the growing season. Climatic factors that influence the excess moisture Typical levels in Timis river discharge in the main situation were analyzed for Timisoara resort and section SP Cruceni are: partialy Dinias-Rauti points, resulting in: - 5% = 80.00 mdMB level (Q = 880 m / s) - Average yearly rainfall 639 mm and 540.4 mm - 1% = 80.75 mdMB level (Q = 1200 m / s) station to station Timisoara evil; - 0.3% = 81.30 mdMB level (Q = 1400 m / s) - Amplitude precipitation: 350-400 mm in dry years - Share dig = 81.43 mdMB canopy. or 500-850 mm in wet years;
Figure 2. Longitudinal profiles in channel CP10
In terms of chemical groundwater presents a 5 1 S 3 mineralization of 0.53 to 5.73 gr / l. Q = i (5) The concentration of salts is lower near permanent n 2 water courses in the area, Bega and Timis. In the P 3 Dinias - Sanmartinu Sarbesc, Ivanda, Foeni appear For the particular case of trapezoidal section is noted mineralization of 4.3 to 5.73 g. / L, which also led to that the area (S) and the wetted perimeter (P) can be secondary soil salinization. expressed as: S = bh+mh2 = β+mh2 (6)
II. MATERIAL AND METHODS 2 2 P = b+2 1+m h = β +2 1+m h (7) Movement of water in canals, rivers or partially filled Resulting in replacement: pipes are examples of free level movements. Uniform 1 β + m1,5+y motion is done free level canals, straight and Q = h2,5+y i (8) n 1,5+y prismatic (canals, galleries, gutters, ditches, etc.). β + m' Free level uniform motion is made in white artificial where the m' were noted: straight and prismatic (canals, galleries, gutters, ' 2 ditches, etc.). Analytical solutions are obtained by m = 2 1+m (9) hydraulic calculation formulas Chezy type by Numerical modelling was performed using the substitutions with Manning relationship with general program MIKE11. MIKE 11 is a professional equations. In general, technical applications currently engineering software package for the simulation of are using Chezy formula type: flows, water quality and sediment transport in V = C Ri (1) estuaries, rivers, irrigation systems, channels and other water bodies. MIKE 11 is a user-friendly, fully Associating continuity equation is obtained for flow: dynamic, one-dimensional modelling tool for the Q =VS = SC Ri (2) detailed analysis, design, management and operation Hydraulic radius is explicit and obtain of both simple and complex river and channel 3 systems. With its exceptional flexibility, speed and S S 2 user friendly environment, MIKE 11 provides a Q = SC i C i (3) complete and effective design environment for P P engineering, water resources, water quality Associating Chezy formula is obtained for flow: management and planning applications. The 1 S1,5+y Hydrodynamic (HD) module is the nucleus of the Q = i (4) n 0,5+y MIKE 11 modelling system and forms the basis for P most modules including Flood Forecasting, By replacing Manning relationship is obtained: Advection-Dispersion, Water Quality and Non-
10 cohesive sediment transport modules. The MIKE 11 Numerical modelling was performed with the HD module solves the vertically integrated equations program MIKE11. Initially modelled the existing for the conservation of continuity and momentum, situation with CP10 trapezoidal channel. Site plan in i.e. the Saint Venant equations. Applications related to the MIKE 11 HD module include: this situation is shown in Figure 5. Cross sections - Flood forecasting and reservoir operation through the channel as topographical surveys are - Simulation of flood control measures shown in Figure 6. - Operation of irrigation and surface drainage systems - Design of channel systems - Tidal and storm surge studies in rivers and estuaries
III. RESULTS AND DISCUSSIONS
Based on these theoretical foundations analytical calculations were performed for proposed pipeline and the existing channel. Topographic studies for existing channel result input the following data: b = 2 m, m = 2, i = 0,0006, n = 0,02 (concrete) The construction curve Q = Q (h) is presented in Figure 3:
Figure 5. Plan view with the network model for the existent trapezoidal channel CP10
Figure 3. The construction curve Q = Q (h)
The volume of water that can be stored in the channel is determined by the relationship: V (b mh)h Lcanal (10) Lm3007 the length of the channel is designed canal Respectively calculate the volume of water that can Figure 6. Cross section for the trapezoidal channel be stored in the drainage channel for different water CP10 depths curve was drawn next channel. The construction curve V = V (h) is presented in Figure 4: According to data entry or formulated boundary conditions, namely the upstream inflow at chainage 3007 are time variant inflow (Figure 7) and in the downstream at chainage 0 curve key for downstream section of the channel. After running the program MIKE11 was obtained through existing channel longitudinal profile, presenting water levels along the channel (Figure 8).
Figure 4. The construction curve Q = Q (h)
11
Figure 7. Time variant inflow.
Figure 8. Longitudinal profile in channel CP10. This study presents the application of a 1- hydroinformatic tools, AGREEMENT NO LLP- dimensional unsteady flow hydraulic model used for LdV-ToI-2011-RO-002/2011-1-RO1-LEO05-5329. the simulation of flow in rivers: the MIKE 11 model This project has been funded with support from the from the Danish Hydraulic Institute (DHI). European Commission. This publication MIKE 11 is the preferred choice of professional river [communication] reflects the views only of the engineers when reliability, versatility, productivity author, and the Commission cannot be held and quality are the keywords. responsible for any use which may be made of the information contained therein. ACKNOWLEDGMENT REFERENCES
[1] Henderson, F.M. (1966). Open Channel Flow. MacMillan Company, New York, USA. [2] David, I. Hydraulic I and II, Polytechnic Institute Traian Vuia This paper can be possible thanks to project: Timisoara, Romania, 1984 Development of knowledge centers for life-long [3] *** Archive ANIF, 2015 learning by involving of specialists and decision [4] *** Mike 11 User Guide, Danmark, pp 1-542, 2012 makers in flood risk management using advanced 12
1 Politehnica University Timisoara, Faculty of Civil Engineering, Department of Hydrotecnnical Engineering, Timisoara, George Enescu Str., no. 1/A, Timisoara, Romania, [email protected] 13