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, , Timis county Robert F. Beilicci1 Erika Beilicci1

Abstract – Study case is situated in irrigation system perimeter system: Cruceni Foieni, Ionel, , 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Ş - 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 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, , 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 = β+mh2 (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 β + m1,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  mh)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