- Open Access
- Total Downloads : 159
- Authors : Sahar Munder Ressen, Prof. Dr. Saleh I. Khassaf
- Paper ID : IJERTV3IS080268
- Volume & Issue : Volume 03, Issue 08 (August 2014)
- Published (First Online): 13-08-2014
- ISSN (Online) : 2278-0181
- Publisher Name : IJERT
- License: This work is licensed under a Creative Commons Attribution 4.0 International License
Development of an Empirical Formula for Computing Sediment Loads in Upstream of Al-Hafar Regulator
Prof. Dr. Saleh I. Khassaf1 Sahar Munder Ressen 2
1CivilDepartment
Faculty of Engineering University of Basrah
2Civil Department Faculty of Engineering
University of Kufa
Abstract:- In this research, the sediments transport and estimating their amount have been studied in upstream of Al- Hafar regulator on Euphrates river which is located in the south of Iraq within Dhi-Qar governorate. Twenty one cross section were selected along the study reach (5 km) in order to study the hydraulic parameters and characteristics of sediments transported.
The study was divided into two parts :the practical part (field and laboratory works) and the statistical part. Samples of water sediment and bed material have been taken at each section using homemade Sampler, while the hydraulic parameters were measured using (Acoustic Doppler Current Profile) ADCP device .
The Empirical Formula for estimating amount of total sediment discharge in study area was developed using dimensional analysis technical and by help SPSS program , the determination coefficient of the Empirical Formula is (R2=0.9946) .
Finally, The average predicted annual total sediment discharge has been estimated through field measurements ,it were (120783) ton ,While The Empirical Formula gave the average annual sediment discharge was (126144) ton.
Keywords: Sediment Transport, Empirical Formula, Al- Hafar regulator, Euphrates river
1. INTRODUCTION
Sediments are small particles like sand, gravel, clay and silt. The water in a river has a natural capacity of transporting sediments. Man made structures in a river may change the sediment transport capacity over a longer part of the river, or locally .Sediment can be classified as deposited or suspended, deposited sediment is that found on the bed of the river or lake, suspended sediment is that found in the water column where it is being transported by water movements, so the material transported by water is in suspensions, rolling or sliding on the bed. The border line between bed load and suspension is certainly not well defined, because it is hard to imagine a particle rolling or sliding on the bed without at some time loosing contact with the bed, executing short jump[1].
WhereSediment transport is important in the fields of sedimentary geology, geomorphology, civil engineering and environmental engineering. Knowledge of sediment transport is most often used to determine whether erosion or deposition will occur, the magnitude of this erosion or deposition, and the time and distance over which it will occur. . In a channel the water flow erodes the available material in the banks and /or the stream bed until the flow is "loaded" with as much sediment particles as the energy of the stream will allow it to carry[2].A derived of formula under the condition of the study reach is very important to estimate the quantity of sediment transport ,so in this paper a development of anew formula for estimation the sediment transport of the study reach.
DESCRIPTION OF AL-HAFAR REGULATOR
Al-Hafar regulator channel is considered as one of the most important projects of irrigation in the region for the large areas that benefited from it and the technical applets and performance.
It is located in Suq Al- Shuyukh area(AL- Karma) (south-east) in Dhi-Qar governorate on Euphrates River. The regulator is constructed in1957and it is consists of seven steel sliding gates in addition to the navigation lock, each gate with a dimension (7×7)m. The designed flood discharge of the regulator is 500 m3/sec . The reach of study upstream of the regulator with length 5 Km .It is located between longitude E 46°34'19.77" to E 46°32'1.62" and latitude N 30°52'52.22" to N 30°52' 9.46
-
FIELD MEASUREMENT
Twenty one cross sections all located up stream of regulator were selected downward the water flow direction with different distance range(100-300)m between section and another as required by the nature of the study area and hydraulic conditions as shown in Fig. 1.The cross sections of river were observed by
taking reference point on the right side bank (with respect to the flow direction) from the reference point at each cross section transect sampling .Cross- section area , velocity ,width ,depth , and other hydraulic properties are determined by passing Acoustic Doppler Current profile ADCP( M9) device along the path of each section which was operated from an anchored boat see Fig. 2. Son Tek River Surveyor (SNR) program starts recording all information after interconnecting with ADCP device. A summary of data used in the study is shows in Table 1.Three samples of bed material were taken at ¼, ½ and
¾ of the width of stream at each cross section in order to conduct the size analysis distribution. The obtained samples were mixed together and part of the mixture was taken to the laboratory for analysis [3].The home- made sampler is called (Van- Veen grab sampler) as
shown in Fig. 3.
Al-Hafar Regulator
Fig. 1: The positions of cross-sections (C.S.) in the region of study, by Google Earth©
Fig. 2: ADCP( M9) adapter
-
Before sampling
-
After sampling
Fig. 3: Bed material sampler
Suspended sediment in river is sampled to determine the sediment concentration .Different sediment samplers and procedures are available for sampling in a stream. Because of the flow acceleration ,there is some effect on the accuracy of measurement, it is imperative that
completion of filtration, the filter paper was dried in oven and re-weighed .The difference between the two weights, divided by the volume of the sample, gives the concentration of the suspended sediment as follow[6]:
sampling process would disturb the flow as little as possible[1]. Depth integrating samplers continuously accumulate a representative sample from a stream vertical. The sampler is lowered slowly down to approximately 10 cm (3.9 in) from the top of the streambed. It is then raised to the water surface at the same speed [4] in the same location of bed material
. The sampler consists of a bottle with the capacity of one liter, an intake nozzle of (8 mm) in diameter. An air escape of (6mm) in diameter with long plastic tube and control valve to control the entering of water sediment mixture into the sampler [5], shown in Fig. 4.
Fig. 4: Suspended sediment homemade sampler
-
-
LABORATORY MEASUREMENTS
The laboratory work is needed to construct the grain size distribution curve of the bed material samples was drown for each sections as a results of sieve analysis and hydrometer test, and to obtain specific gravity value of the observed sediment samples, also to find sediment concentration.
4.1 Sediment Concentration Measurement
The concentration of suspended sediment was determined from the samples that were taken from specified locations of a study reach by using a filtration method Fig. 5.
A well- mixed sample is filtered using a suitable filter and the residue retained on the filter is dried in oven at a constant (103-105) co.Each filter paper was pre-dried for 60 minutes and weighted, and then it was clipped to the filter funnel and moistened with distilled water. A volume of 250 ml of the sample was measured into a graduated cylinder and poured through the filter, and all interior surface of the cylinder were washed out into the filter funnel with distilled water. After the
C= Concentrtion of suspended sediment in mg/l (P.P.m).
W1= Weight of dry filter paper in mg.
W2= Weight of dry filter paper + suspended sediment in mg.
V= Volume of sample (l).
Fig. 5: Filtration method
Sec. No
1
2
3
Qw(m3/sec)
34.62
34.23
39.79
V (m/sec)
0.316
0.158
0.32
GS
2.73
2.66
2.65
d50 (mm)
0.152
0.113
0.112
A (m2)
109.4
216.3
124.3
B (m)
55.01
107.44
70.83
Rh (m)
1.99
2
1.75
(m2/sec)
1.05×10-6
1.05×10-6
1.05×10-6
Ws (m/sec)
0.016668
0.009407
0.009199
U* (m/sec)
0.0312
0.0313
0.0293
Sec. No
4
5
6
Qw(m3/sec)
40.37
40.87
42.56
V (m/sec)
0.312
0.345
0.343
GS
2.69
2.73
2.7
Table 1: Primary data and parameter
.(1)
Sec. No
16
17
18
Qw
(m3/sec)
27.82
13.48
14.76
V (m/sec)
0.257
0.123
0.162
GS
2.69
2.67
2.68
d50 (mm)
0.113
0.127
0.147
A (m2)
108.2
109.3
91.1
B (m)
50.01
69.22
64.99
Rh (m)
2.16
1.58
1.40
(m2/sec)
1.06×10-6
1.20×10-6
1.21×10-6
Ws (m/sec)
0.009462
0.010445
0.013612
U* (m/sec)
0.0325
0.0278
0.0262
Sec. No
19
20
21
Qw
(m3/sec)
13.93
12.46
14.10
V (m/sec)
0.135
0.117
0.135
GS
2.7
2.69
2.72
d50 (mm)
0.123
0.11
0.1
A (m2)
103.6
106.2
104.2
B (m)
52.18
44.11
47.94
Rh (m)
1.98
2.41
2.17
(m2/sec)
1.21×10-6
1.21×10-6
1.22×10-6
Ws (m/sec)
0.009915
0.007970
0.006724
U* (m/sec)
0.0312
0.0344
0.0326
d50 (mm)
0.105
0.1
0.102
A (m2)
120.1
116.7
121.3
B (m)
47.77
51.05
60.02
Rh (m)
2.51
2.28
2.02
(m2/sec)
1.06×10-6
1.07×10-6
1.07×10-6
Ws (m/sec)
0.007543
0.008431
0.008958
U* (m/sec)
0.0351
0.0334
0.0315
Sec. No
10
11
12
Qw (m3/sec)
22.97
22.95
24.97
V (m/sec)
0.25
0.195
0.23
GS
2.65
2.66
2.7
d50 (mm)
0.1
0.098
0.1
A (m2)
92
117.5
108.8
B (m)
60.19
49.60
74.36
Rh (m)
1.53
2.37
1.46
(m2/sec)
1.07×10-6
1.07×10-6
1.07×10-6
Ws (m/sec)
0.007294
0.008096
0.007508
U* (m/sec)
0.0274
0.0341
0.0268
Sec. No
13
14
15
Qw (m3/sec)
26.36
27.40
27.42
V (m/sec)
0.256
0.243
0.234
GS
2.69
2.68
2.7
d50 (mm)
0.099
0.101
0.125
A (m2)
102.9
112.6
116.9
B (m)
56.72
51.39
46.65
Rh (m)
1.81
2.19
2.50
(m2/sec)
1.06×10-6
1.07×10-6
1.07×10-6
Ws (m/sec)
0.008454
0.007599
0.011436
U* (m/sec)
0.0298
0.0327
0.035
-
SEDIMENT DISCHARGE MEASUREMENT
The measurement of sediment discharge is essential to determine the quantity of sediment load to establish or check analytical or empirical sediment transport equations. Then sediment discharge can be calculated by multiplying the concentration with the flow
discharge [2]. As the following
Qs=C ×Q× 0.001 ..(2)
Table 2 illustrates the calculated sediment discharge resulting from Equation (2).
Table 2: Field sediment discharge for each section
The number of primary dimensions involved is 3, i.e., m = 3 (M, L, T). The total numbers of variables are
- p>Therefore, the number of -terms is 7
Thus, F {1, 2, 3, 4, 5, 6, 7} = constant .(5) The repeating variables are selected ( ,ws , Rh) i.e
,the first variable representing the fluid property ,the second representing the sediment characteristics and the third representing the hydraulic property. The results of the analysis are listed in Table 3.
Table 3 : parameters
1
2
3
4
5
6
7
parameter
Then, the equation can be expressed as the following:
No. Sec.
water discharge (m3/sec)
Average Concentration (g/m3)
Sediment discharge (Kg/sec)
1
34.62
165.45
5.728
2
34.23
155.5
5.323
3
39.79
167.825
6.678
4
40.37
164.36
6.635
5
40.87
171.52
7.010
6
42.56
168.33
7.164
7
22.77
125.29
2.853
8
23.58
132.5
3.124
9
23.82
130.25
3.105
10
22.97
150.33
3.453
11
22.95
135.64
3.113
12
24.97
140.53
3.509
13
26.36
147.25
3.881
14
27.40
144.44
3.957
15
27.42
141.83
3.889
16
27.82
145.45
4.046
17
13.48
101.35
1.366
18
14.76
118.25
1.745
19
13.93
112.55
1.568
20
12.46
95.78
1.193
21
14.10
114.37
1.613
= F ( ) ( )
..(6)
6. DEVELOPMENT OF EMPIRICAL FORMULA
The dimensional analysis is a good way in dealing with a
The following procedure was followed to reduce the number of -terms:
3 and 4 may be replaced by another 8
3 × 4 =8 =( ) × = .(7)
complex problem if it is Correctly applied. The principal use of dimensional analysis is to deduce from a study of the dimensions of the variables in any physical system certain
Using the same approach, 2 and 7 can be combined to evaluate 9
limitations on the form of any possible relationship between those variables. The method is of great generality and mathematical simplicity [7] The result of
= 9 =
=
.(8)
dimensional analysis depended on the most important variables are selected according to groups[8].
Thus, the functional relationship becomes:
= F
Group 1: Variables related to the characteristics of sediments
.(9)
(d50 ,ws ,s , QS)
Group 2: Variables related to geometric and hydraulic properties of stream channel .
(B, Rh ,V, ,,U*)
Using Buckingham's -theorem procedure as presented in [9], the variables used for the field and the laboratory work and their relationship are as follows:
QS = (d50, ws ,s ,B, Rh ,V, ,,U*) .(3) Or F(QS, d50, ws ,s ,B, Rh ,V, ,,U*) = constant
.(4)
The final form of the equation has to be determined by
the conducting of the regression analysis with help of SPSS program on the observed data. These data were divided into two groups: the first group including of 13 different sections was randomly selected for the derivation of the new equation while the second group (8 sections) was used for verification of Empirical Formula .The regression analysis was conducted and be found by using the following formula:
..(10)
The coefficient of determination of equation (10) was found to be equal (R2=0.9946). Fig. 6 shows a well accepted relationship between the predicted and the observed values of sediment discharge for 13 sections.
R² = 0.9946
-
COMPARISON USING STATISTICAL RELATIONS
Three methods are used in this research to evaluate the performance Empirical formula through comparing with measured values.
-
Mean Standard Error
A mean standard error was used in order to select the best formula since due to the high difference between predicted and measured sediment rates at various intervals [10].
8
7
6
5
4
3
2
1
0
Observed Values
| |
(11)
In which:
0 1 2 3 4 5 6 7
8
Predicted Values
Fig. 6: Observed and predicted values for (13) sections.
7. VERIFICATION OF THE PROPOSED FORMULA
To verify that the proposed formula has been used eight remaining sections of the sediment discharge measured and predicted a long the area of study reach. It demonstrates the coefficient of determination (R2 = 0.9948). This step provides an independent verification of the precision of the mentioned formula because non of the data are used to obtained this proposed formula.
is a Mean Standard Error; an observed sediment rate; is a predicted sediment load and is the number of the predicted values.In this method, a lower statistical criterion (close to zero) shows a higher accuracy in the model performance.
The Empirical Formula gave the Mean Standard Error ( ) equal to 5.35%. Thus it is produced quite good performance to estimate the amount of bed material load in the study area comparing with measured values.
-
Discrepancy Ratio
Discrepancy ratio[11] is defined as the ratio between computed and measured sediment loads. It was used as an error measure that is calculated as:
A well agreement between the observed and the predicted sediment discharge can be realized for eight
discrepancy Ratio (R) =
.(12)
sections in Fig. 7.
When the discrepancy ratio is equal to one (R=1)[10] for the value of the Empirical Formula that indicate the predicted value is identical to the measured value for reach of study .The discrepancy ratio is scheduled with the ranges (0.75-1.25) ,(0.5- 1.5) ,and (0.25-1.75). The results are 100% for all ranges.
-
Root Mean Squared Error
The root mean square error calculation is a well known and frequently used method of error analyss. It accurately depicts the magnitude of deviations of and estimated (measured or calculated) value from the actual value sought[12]. The RMSE has the same units as the measured and calculated data. Smaller values indicate better agreement between the measured and the calculated values[13][14].
8
7
R² = 0.9948
6
5
4
3
2
1
0
0 1 2 3 4 5 6 7 8
Predicted Values
Observed Values
Fig. 7: The Observed and predicted values for the eight sections
..(13)
In which: observed sediment rate, is predicted sediment load and is the number of predicted values.The Empirical Formula gave equal to 0.27% .
-
-
CONCLUSIONS
-
According to the results which are obtained by this study for twenty one cross sections on Euphrates river up-stream of Al- Hafar regulator ,in Al- Nasiriyah city the following points are concluded:
A good agreement was observed between the measured values and the computed values of the total sediment discharge for Empirical Formula which it has been developed in term of five dimensionless parameters:
-
Ongly.,E., , Sediment Measurements , chap. 13, In:"Water Quality Monitoring" , (Ed.), Design And Implementation Of Fresh Water Quality Studies And Monitoring Programmes, World Health Organization, Geneva,1996 { Ongly.,E., , Sediment Measurements , chap. 13, In:"Water Quality Monitoring" , (Ed.), Design And Implementation Of Fresh Water Quality Studies And Monitoring Programmes, World Health Organization, Geneva,1996}.
-
Jasem , H., Mohammed Estimation of Sediment Quantity up stream of Al-Abbasiya barrage in Euphrates River , MSc. Thesis, Department of Civil Engineering, University of Kufa, 2012.
-
Sonion ,A.A."The physical Basis of Dimensional Analysis " Department of Mechanical Engineering, MIT, Cambridge, MA 02139,2001.
( and
This formula is based on thirteen section of data and other eight sections are used to certify it. The Empirical Formula is applicable under the following conditions:
-
Flow velocity (0.117 0.345) m/sec.
-
Water discharge (12.46 42.56) m3/sec.
A verage concentration of suspended load (95.78 171.52) mg/l.
-
Median grain size (0.098 0.152) mm.
From through statistical methods used for comparison noted that the Empirical Formula proposed is good accurate in calculation the amount of sediment in the search area comparing with measured values. Where the average predicted annual total sediment discharge has been estimated through field measurements ,it were (120783) ton ,While the Empirical Formula gave the average annual sediment discharge was (126144) ton.
REFERENCES
-
Addab, H. F: "Estimation Of Sediment Quantity Of Al- Meshkab Regulator Channel", M.Sc. Thesis, Department of civil engineering. University of Kufa,, 2011.
-
Van Rijn, L. C.," Principles of Sediment Transport in Rivers, Estuaries and Coastal Seas", Aqua Publications, Amsterdam, the Netherlands, 1993.
-
Adegbola, A., A., and Olaniyan, O., S.," Estimation of Bed LoadTransport in River Omi, South Western Nigeria using Grain Size Distribution Data", Department of Civil Engineering, University of Technology, Ogbomosho, International Journal of Engineering andTechnology Volume 2 No. 9,ISSN 2049-3444, 2012.
-
Lisa Frilay."Methods of Measuring fluvial sediment" MSc Thesis,Center for Urban Environmental Research and Education, University of Maryland, Baltimore County, Technology Research Center 102, January 21, 2004.
-
-
-
Simons, D. B., and Senturk, F., Sediment Transport Technology, water Resources Publications, Fort Collins, Colorado,
USA,1977
-
Arora, K., R., "Fluid Mechanics, Hydraulics and Hydraulic Machines", 1705-B, NAI SAPAK, post Box No.1106, DELHI, 2007.
-
Hassanzadeh, H., Faiznia, S., Shafai, B., M. and Motamed, A., "Estimate of Sediment Transport Rate at Karkheh River in Iran Using Selected Transport Formulas", World Applied Sciences Journal 13 (2): 376-384, ISSN 1818-4952, © IDOSI Publications, 2011.
-
Yang, C. T., "Erosion and Sedimentation Manual", US Department of the Interior, 2006.
-
Johnson, L., L.," A Comparison of Methods for Estimating RMS Error: A "Brute Force" Approach Versus a Mathematically-Elegant Approach, as Applied to the Calculation of a Specific Retrieval Error for a Limb-Scanning Microwave Radiometer-Spectrometer" ,M. Sc. Thesis,the Faculty of the School of Engineering Air Education and Training Command AFIT/GAP/ENP/95D- 10, December ,1995.
-
Scheaffer, R., L., "Probability and Statistics for Engineer",brooks/Cole, USA, 2011.
-
Sadiq ,N., Mohammed" Evaluation of Sediment Transport up Stream of Al-Shamia Barrage"" MSc. Thesis, Department of Civil Engineering, University of Kufa, 2013.
-
Notations
Symbol |
Meaning |
Unit |
Total sediment load |
kg/sec |
|
Density of fluid |
kg/m3 |
|
Fall velocity of particle |
m/sec |
|
Hydraulic radius |
m |
|
Mean velocity |
m/sec |
|
Shear velocity |
m/sec |
|
Specific gravity |
— |
|
Width of river |
m |
|
Median grain size |
m |
|
Kinematics viscosity |
m2/sec |
|
Q |
water discharge |
m3/sec |
C |
The sediment concentration |
PPm |
S |
Density of sediment |
kg/m3 |