International Journal of Geosciences
Vol.3 No.5(2012), Article ID:24988,13 pages DOI:10.4236/ijg.2012.35112

Relationships between Hydraulic Parameters of the Nubian Aquifer and Wells in El Shab Area, South Western Desert, Egypt (A Case Study)

Maged El Osta

Geology Department, Faculty of Science, Damanhour University, Damanhour, Egypt

Email: drmagedelosta.edu.alex@hotmail.com

Received August 15, 2012; revised September 14, 2012; accepted October 12, 2012

Keywords: Nubian Sandstone Aquifer; Transmissivity; Specific Capacity; Well Loss; Entrance Velocity; Well Efficiency

ABSTRACT

For proper water management in the new reclaimed areas, hydraulic parameters of both aquifer and wells related to transmissivity, specific capacity, well loss, formation loss, and water entrance velocity, as well as the relationship between these parameters are the main target after construction of production wells. In El Shab area, the Nubian Sandstone aquifer has a large range of transmissivity (from 483.12 to 1489.24 m2/day) and, also, specific capacity (from 203 to 486.32 m2/day). Relationship of specific capacity and transmissivity is constructed and the established empirical equations can be used to predict the transmissivity of the Nubian aquifer in all new proposed sites for well drilling at which the specific capacity measured without performing pumping tests. On the other hand, the drilled wells in El Shab area exhibit relatively high well losses (25%). The causes for high well losses (entrance velocity (Vn) through water well screen and the distance from the point of water entrance in the well to the point of intake in the pump) are discussed and the relationships are constructed, which seem a positive linear correlation. Relationships between well losses constant for 30 wells with transmissivity and specific capacity, are constructed. These relationships are useful for estimating hydraulic characteristics that are needed for the designs of wells and well fields and for preliminary water-resources management.

1. Introduction

El Shab area, lies in south of the Western Desert and is considered as natural extension of Darb El Arbaein project (Figure 1). In Darb El-Arbaein, the groundwater is the only water resource. The aquifer system starts from Paleozoic-Mesozoic to Upper Cretaceous sandstone rocks overlaying the basement rocks with confined condition. They irrigate 12,000 feddans (1 feddan ≈ 1.038 acres) and help 16 villages for farmers along about 400 km from Paris town towards the Egypt-Sudan border. The pumping wells are located in three parts: the northern, the central, and the southern part of Darb El-Arbaein. The northern part extends 90 km to the south from Paris town and has an area of 90 km2. The central part extends 80 km to the south of the northern part and has 120 km2 in area. The southern part (El Shab area) extends 200 km to the south of the central part and has 170 km2 in area. The Darb El-Arbaein is located in an arid region where the absolute maximum air temperature is 48.6˚C in May and minimum air temperature is 2˚C in February. The annual rainfall is <1.1 mm, and the total monthly rainfall is 0.3 mm. The annual evaporation is about 172 mm, and the maximum evaporation rate in June is about 21.32 mm. The evaluating of water resources and their use is main target for proper management in this new reclaimed area. The purpose of ground-water development in this arid region, is to bring about an additional supply of water, whether it be used for agricultural, domestic, industrial, or other purposes. An understanding of the various factors causing excessive well losses is important because savings in well design and construction and operation can be made by increasing the efficiency of a well and thus preventing much needless waste of materials pumping energy.

In this paper, an attempt is made to evaluate the factors causing excessive well losses in El Shab area and to relate the wells and aquifer hydraulic characteristics to find out the relationship between these hydraulic parameters. These relationships are useful for estimating hydraulic characteristics that are needed for the designs of wells and for preliminary water resources management.

Figure 1. Location map of El Shab area, South Western desert, Egypt.

2. Geological and Hydrogeological Conditions

Geologically, the exposed rocks in Darb El-Arbaein range from Pre-Cambrian to Quaternary sediments (CONOCO [1]). The lithostratigraphic successions are divided into seven units, from base to top (CONOCO [1]; Fathy et al. [2]; Korany et al. [3]; Ghazal [4]; El Gammal [5]): 1) Pre-Cambrian basement; 2) Paleozoic-Mesozoic sandstone; 3) Lower Cretaceous; 4) Upper Cretaceous; 5) Paleocene; 6) Eocene; and 7) Quaternary. The Nubian sandstone is classified into Shab, Qusier, and Taref members (Issawi [6]). The Kiseiba Formation is equivalent to both Shab and Qusier members and Taref is equivalent to Taref member. The Dakhla shale represents the capping aquifer only in the central and northern Darb El-Arbaein. The sedimentary cover increases far in the northern part (more than100 m), while it decreases to the southern part of the northern Darb El-Arbaein. The lithostratigraphic sequence in the southern part of Darb El-Arbaein (El Shab area) is mainly composed of Kiseiba and Taref Formations with a little thickness variation (Figure 2).

Hydrogeologically, 30 production wells are drilled and constructed (year 2001) in the Nubian Sandstone aquifer by the Research and Ground Water Co. (REWA) in a project for developing the area of the New Valley Governorate (Figure 1). The regional hydrogeological section constructed along Darb E1 Arbaein area (Figure 3) shows that the Nubian aquifer in the concerned area is the unique source of fresh water and built of Kiseiba and Taref Formations of the Upper Cretaceous. The upper most layer of Kiseiba Formation acts as a confining layer. The aquifer thickness increases in the southern part, attributed to high concentration of faults to the north. The aquifer is composed of coarse-grained sandstone at the base that changes gradually to fine and medium upward. It directly overlies the basement rocks. The drilled depth of wells is running almost between 120 and 210 m (reached the basement complex). The details of drilling and well design are shown in Table 1, which includes the whole wells. After deciding the production well assembly, the pilot hole is reamed to proper diameter and the well assembly is set in position by placing the specified slotted section against the water zones. Gravel of proper size is shrouded between the well assembly and the walls of the borehole, and the well is washed thoroughly. The wells are developed by backwashing and airlifting. The lifting pump is sited for a depth between 25 and 42 m only, which may not be suitable for the expected amount of water withdrawn from the wells (200 m3/h/well and for 12 hours daily) as well as the expected head decline in wells. The estimated transmissivity of the Nubian aquifer from pumping tests analysis ranges between 483.12 m2/day

Figure 2. Lithostratigraphic sequence in El Shab area, South Western desert, Egypt.

Figure 3. A very generalized hydrogeological cross section along the three areas of Darb El Arbaein.

(well No. 62) and 1489.24 m2/day (well No. 89). This wide variation may be attributed to the rapid lateral facies changes as well as the variation in the thickness of the productive water bearing units.

3. Materials and Methods

The materials used in this paper were collected through carrying out two field trips in E1 Shab area in 2010 year, in which the following activities were performed: 1) Inventory of the existing water wells; 2) Carrying out, depth to water measurements from the ground surface in 30 water wells; 3) Registration of wells discharge. The other basic hydrologic materials were got from the Ministry of Water Resources and Irrigation, Groundwater sector, El Wadi El Gidid area during these field trips.

Table 1. Data for wells in El Shab area (Data source: ministry of water resources and irrigation, groundwater sector, El wadi El gidid area[7]).

These materials include collection of archival data such as well drilling reports, pump position, screen length, proposed operating systems for both groundwater supply and reclaimed area beside step drawdown test data of 30 drilled wells. In addition, the methods applied in this work consists of the following steps:

3.1. Calculating of Specific Capacity (Sc)

Specific capacity is a measure of the productivity of a well. It is determined from a test that involves pumping a well at a constant rate and measuring the resulting drawdown in water level. Specific capacity is the parameter most often provided by drillers to characterize the performance of a well, both because it is easy to be measured and because well owners easily understand how much water they can get from the well. In contrast, estimation of transmissivity (T) requires more rigorous, longer-term testing, which is expensive. This explains why (Sc) data are much more abundant than (T) data. As, transmissivity is the preferred parameters for many hydrogeologic purposes, it is common practice to calculate it from specific capacity data through the relationship established between these parameters. In mathematical terms, specific capacity (Sc) is defined as the pumping rate (Q) divided by the observed decline in hydraulic head in the well (∆hw) and has the same dimensions of transmissivity (m2/day) as follows:

(1)

For steady-state conditions, specific capacity of a well is a function of well radius, degree of aquifer penetration and transmissivity. It may also be increased by increasing the radius of the well and by increasing the percent penetration of the total saturated thickness of the aquifer. For transient conditions, it is also a function of time, well loss and storativity.

Well loss leads to miscalculation of specific capacity from what would be expected if specific capacity only reflected the productivity of the aquifer, as the total head decline in a well (∆hw) is the sum of the head loss in the aquifer (∆ha) at the well and the well loss (∆hL), which represents the part of the drawdown in a well that is due to turbulent flow near the well, through the well screen and in the well casing.

(2)

which is also:

(3)

where B is the well function and C is the well-loss constant (Jacob [8]; Eagon and Johe [9]). If there is no well loss (C equals 0) the specific capacity only reflected the productivity of the aquifer. If C is greater than 0, well loss (GhL) is nonzero, and the observed decline of hydraulic head in the well (Ghw) becomes larger leading to miscalculation of specific capacity of the well. Calculated specific capacity (Sc) for the 30 drilled well in El Shab area is shown in Table 1. The specific capacity, thus calculated, ranges between 203 m2/day in well No. 73 (low productivity) and 486.32 m2/day in well No. 89 (high productivity). This relatively wide variation may be due to the variation in the well loss of the productive water, which in turn is controlled by the placing of specified slotted section, total screen open area, gravel size between the well assembly and the walls of the borehole, and the position of lifting pump.

3.2. Analysis of Step-Test Data

The step drawdown pumping test is very important to determine the definite discharge and the efficiency of both aquifer and the wells. The drawdown in a well created due withdrawal of water, is usually made up of two types of head losses resulting from laminar flow in the formation (formation loss) and turbulent flow through the well screen and the well casing (well loss).

A step drawdown test data of 30 drilled wells in El Shab area were analyzed by the author using graphical method according to the following equation of Rorabaugh [10] (Kruseman and de Ridder [11]):

(4)

Where, S is the anticipated drawdown (m), Q is the discharge rate (m3/h), B is the formation loss coefficient (h/m2), C is the well loss coefficient (h2/m5).

By dividing the drawdown by the discharge rate of the well, the above equation changes as:

(5)

The ratio S/Q is the specific drawdown (reciprocal of specific capacity) of the well.

On the other hand, the well efficiency percentage () at any pumping rate can be also calculated by applying the following equation:

(6)

By plotting the specific drawdown of the well against discharge rate for several steps, the well loss and formation loss coefficients (C and B) were determined from the slope and intercept of the straight line plots respectively. Two of these graphical representations are shown in Figure 4. Based on these graphical representations the values of the well loss (CQ2), formation loss (BQ), well loss percentage (CQ2%) and the well efficiency percentage () are calculated for each step and listed in Table 2.

The obtained results indicate that, the well loss coefficient ranges between 2.9 × 10–5 (h2/m5) in well No. 63 and 1.57 × 10–4 (h2/m5) in well No. 69, while the formation loss coefficient varies from 0.0324 (h/m2) in well No. 75 to 0.0942 (h/m2) in well No. 72. Moreover, the estimated total drawdown (formation loss + well loss) by the graphical method is very close to the actual drawdown observed in the well. However, the formation loss percentage (well efficiency) varies from 52% (well No. 69) to 98% (well No. 70), with an average of about 75%. In contrast, the well loss percentage ranges between 2% (well No. 70) and 48% (well No. 69), with an average of about 25%. These results indicate that the drilled wells in El Shab area exhibit relatively high well losses. For the proper design of a well, attempts should be made to keep the well loss to a possible minimum. It is therefore, necessary to find out the factors that cause the excess well losses in the well drilled in El Shab area. These factors may be attributed to: 1) the partly clogged interstices of the gravel back around the screens of the wells; 2) the low area of screen openings; and 3) the high distance from the point of water entrance in the well to the point of intake in the pump (pump position).

3.3. Calculating of Entrance Velocity (Vn)

Among the many important parameters that take into account in the design and operation of water wells for

Table 2. Well and formation characteristics obtained from graphical method, in El Shab area.

     

Figure 4. Graph showing the relation between the specific drawdown (S/Q) and discharge (Q).

proper well management, is entrance velocity (Vn) through water well screen. It is suggested that an entrance velocity of 0.03 m/second to 0.0762 m/second barely lifts sand 0.25 mm to 0.50 mm diameter, whilst a velocity of 0.018288 m/second barely lifts clay or silt (Smith [12]). It is readily acknowledged that entrance velocity is an integral component of well design, model studies clearly indicate that the range of a relatively acceptable entrance velocity is actually considerably 0.03 m/second (R.M.C [13]). The usual index, used to determine the tendency of flow to become laminar or turbulent, is the Reynolds number (NR) which is defined as follows:

(7)

where V is the velocity (m/min), d is the grain diameter, and υ is the viscosity of water.

Experimental work on uniform spherical sand grains influences that departure from laminar flow begins at value of NR between 1 and 10. Secondly, if the values of d and υ are nearly the same in a particular area, the increase in the value of V causes more turbulence resulting in higher well loss. The entrance velocity (Vn), is defined as the pumping rate of flow (Q) in m3/h divided by the working area of screen (Ar) in m2 as follows:

(8)

Using this equation, entrance velocities for the proposed discharging rate (200 m3/h/well) for all drilled wells except well N0. 61 (250 m3/h) in the concerned area were calculated and tabulated in Table 3. For well design in the study area, the working area of screen equals 0.05 m2 for each meter. The entrance velocities, thus estimated vary from 31.46 m/h (well No. 75) to 89.72 m/h (well No. 84) with an average of about 45.66 m/h, which is quite high.

4. Results and Discussion

As was stated before, the transmisivity (T) and specific capacity (Sc) depend on the well loss in a pumping well, which in turn resulting from: 1) the turbulent flow in the zone outside the well i.e. in the gravel pack, 2) the turbulent flow through the well screen (entrance velocity of water), and 3) the turbulent flow inside the well to the pump intake. So, in the present work, an attempt was made to find out the relationship between these hydraulic parameters for the purpose of well and aquifer proper management as follows.

4.1. Empirical Relationship between Transmissivity and Specific Capacity

Specific capacity is useful for estimating transmissivity in aquifers that have few good pump tests. One approach of estimating transmissivity from specific capacity is to use analytical equations such as those derived by Thomasson et al. [14], and Theis [15]. However, Razack and Huntley [16], and Huntley et al. [17] showed that the analytical solutions assume negligible well loss and a homogenous, isotropic, granular aquifer. Instead of using analytical equations, Razack and Huntley [16], and Huntley et al. [17] developed empirical relationships between transmissivity and specific capacity for each aquifer.

In the concerned aquifer, the relationships are usually derived from a first order linear regression using the values of transmissivity (from pumping) and measured specific capacity both in m2/day units listed in Table 1. The best-fitting line between these parameters indicates a positive relationship with moderate coefficient of determination, R2, equals 0.370929 as follows (Figure 5):

(9)

Table 3. Entrance velocity of water in wells, for the proposed discharging rate in El Shab area.

This established empirical equation can be used to predict the transmissivity of the Nubian aquifer in all new sites for well drilling at which the specific capacity measured without performing, longer-term pump testing, which is expensive, but with low potential error of prediction.

4.2. Empirical Relationship of Well Losses Constant with Transmissivity and Specific Capacity

The empirical relationships between well losses constant (C in d2/m5) for 30 drilled wells in El Shab area, with transmissivity (T) and specific capacity (Sc) both in m2/ day units with their best fit lines are shown in Figure 6.

The specific empirical equation for the transmissivity with very low coefficient of determination, R2, equals 0.000595893 is as follow:

(10)

The slope of best fit line for these two parameters (C and T) does not indicate a definite trend. However, a correction can be made on transmissivity values that can be estimated from the well loss constant obtained from graphical method of step-test data with the quite acceptable of the error percentages.

In contrast, the relationship between well losses constant and specific capacity shows a negative linear relation with low coefficient of determination, R2, equals 0.0657616 as follows:

Figure 5. Relationship between transmissivity and specific capacity, showing the best-fit line.

Figure 6. Relationship between well losses constant with transmissivity and specific capacity, showing the best-fit lines.

(11)

This equation shows that, the productivity of a well increases as the well loss constant decreases. If (C) equals zero the specific capacity only reflected the productivity of the aquifer with no miscalculation.

4.3. Relationship of Well Loss Percentage for a Proposed Discharge Rate to Both Entrance Velocity and Distance from the Entrance Point to the Intake Point

In the concerned area, the factors which cause the main parts of well losses (%) are the entrance velocity of waters in wells (Vn in m/h) through the working area of screen and the distance from the entrance point to the intake point (D in m). The percentages of well losses estimated for a proposed discharge rate from the wells listed in Table 4 are plotted against the calculated entrance velocity for the same discharging rate as well as the distance from the pump to the first slotted pipe in each well (Table 3).

The scatter of points for nearly 85% of the wells on the graph for the relation between well losses (%) and entrance velocities (Figure 7), indicates a relatively positive and fitting relation between these factors with low coefficient of determination, R2, equals 0.0107456 as follows:

(12)

This equation shows that the well loss increases with higher values of entrance velocities and accordingly with the lower values of screen working area.

Moreover, the plot between Well losses (%) and the distance from the entrance point to the intake point (Figure 7), confirms the dependence of Well losses (%) with (D), with low coefficient of determination, R2, equals 0.0788977 as follows:

(13)

Equation (13) shows that the well loss factor increases as the distance from the point of entrance of water into the well to the pump intake increases.

Accordingly, it can be inferred that these might be the main causes for high well losses in the drilled wells constructed in El Shab area. Another factor contributing to the well loss is the turbulent due to the clogging of gravel pack itself by the formation materials, the pack becomes ineffective in its purpose and causes higher well losses.

Based on the above discussion and my experience in other areas, it is felt that this might form a considerable

Table 4. Well losses estimated for a proposed discharge rate from the wells, in El Shab area.

    

Figure 7. Relationship between well loss (%) with entrance velocity and distance from the entrance point to the intake point, showing the best-fit line.

part of the total well loss. For the well designer, this has important relevance when selecting the type and amount of well screen for a new well as well as the pump position.

5. Conclusions and Recommendations

In El Shab area, groundwater of the Nubian Sandstone aquifer is the unique source of fresh water. This aquifer has a large range of transmissivity (from 483.12 to 1489.24 m2/day) and, also, specific capacity (from 203 to 486.32 m2/day). It has been shown that, the percentage of well loss in the drilling wells varies from 2% to 48% with an average of about 25%. In contrast, the well efficiency varies from 52% to 98% (well No. 70), with an average of about 75%. For the present discharge rate, 10 wells having the well loss less than 20% (well efficiency <90%), and the other wells (20 wells) in the range of 21% - 48% (well efficiency between 80% and 50%). This leads to the understanding that in nearly 70% of the drilled wells in the area, excessive well losses are observed. An attempt was made to evaluate the factors causing excessive well losses so that remedial measures can be adopted in the design aspects of the drilled wells. The calculations of entrance velocities of water into the screen have shown that all values are failing in the range 31.46 to 89.72 m/h with an average of about 45.66 m/h, which is quite high. As such, the full open area of the slots, assumed for the calculations for entrance velocity, may not be available for the water to flow freely, which in turn causes excess well losses. Another factor contributing to the well loss is turbulent flow inside the well due to the high distance from the screen to the pump intake in the drilling wells (from 25 to 42 m). The established relationships between well losses constant with transmissivity and specific capacity indicate, indefinite clear trend relation for transmissivity with very low coefficient of determination, R2 (0.000595893), and a negative linear relation for specific capacity with a low coefficient of determination (0.0657616) respectively.

Based on the above results, the following recommenddations should be taken into consideration in the design of water wells for proper management:

1) Well should have better graded gravel packs. The angular pack materials should be avoided and consequently clogging of formation material in the same.

2) Maximizing the open area of screen than has adopted at present to minimize the entrance velocity of water into the well as well as the well loss percentage. Increasing the present screen working area by 20% (0.07 m2 for each screen meter) will decrease the entrance velocity by about 25% and hence the well loss percent by about 6% (Equation (12)).

3) Screen materials should be corrosion resistant.

4) The distance from the pump intake to the point of water entrance into the well (screen position) has a bearing on the well losses. In order to minimize the well losses, the pump may have to be lowered to depth closer to the screen, taking into consideration the economics involved in the lift vis a vis a maximum limit of permissible well loss. Increasing the current distance between the pump and the screen by 20% will increase the well loss percent by 7% (Equation (13)).

REFERENCES

  1. CONOCO, “Geologic Map of Egypt, Scale 1:500,000,” Continental Oil Company, 1989.
  2. R. Fathy, M. Hefnawy and A. Abdel Hamid, “Contributions to the Hydrogeology Aspects of the Groundwater Aquifer in Southern Part of Darb El Arbaein Area, Southern Portion of Western Desert, Egypt,” Al-Azhar Bulletin of Science, Vol. 12, No. 2, 2001, pp. 175-194.
  3. E. Korany, G. Fathy and M. El Nagaty, “Contributions to the Hydrogeology of Nubian Sandstone Aquifer in the Middle Part of Darb El Arbaein, South Western Desert, Egypt,” Sedimentology of Egypt, Vol. 10, 2002, pp. 119- 143.
  4. A. Ghazal, “Hydrogeological Studies in Darb El Arbaein area, Southwestern Desert, Egypt,” Master Science, Geology Department, Cairo University, Cairo, 2002.
  5. N. El Gammal, “Hydrogeological Studies in Darb El Arbaein Area, South Egypt,” Master Science, Geology Department, Cairo University, Cairo, 2004.
  6. B. Issawi, “Geology of Southeastern Corner of Western Desert,” Ann GSE Cairo, Vol. 8, No. 3, 1973, pp. 25-30.
  7. Ministry of Water Resources and Irrigation (MWRI), “Hydrogeological Data of Water Wells in Darb El Arbaein Area, Western Desert, Egypt,” Internal Unpublished Report, 2010.
  8. C. E. Jacob, “Flow of Groundwater. Engineering Hydraulics,” John Wiley, New York, 1950.
  9. H. B. Eagon Jr. and D. E. Johe, “Practical Solutions for Pumping Tests in Carbonate—Rock Aquifers,” Ground Water, Vol. 10, No. 4, 1972, pp. 6-13. doi:10.1111/j.1745-6584.1972.tb02929.x
  10. M. J. Rorabaugh, “Graphical and Theoretical Analysis of Step Drawdown Test of Artesian Well,” Proceedings of the American Society of Civil Engineers, Vol. 79, No. 362, 1953, pp. 1-23.
  11. G. P. Kruseman and N. A. de Ridder, “Analysis and Evaluation of Pumping Test Data,” 2nd Edition, Wageningen Academic Publishers, Wageningen, 1990.
  12. E. R. Smith, “Screen Bore Construction,” Sydney, Water Research Foundation of Australia, Report No. 24, 1966.
  13. Roscoe Moss Company, “Handbook of Ground Water Development,” John Wiley and Sons, New York, 1990.
  14. H. J. Thomasson, F. H. Olmsted and E. R. LeRoux, “Geology, Water Resources, and Usable Groundwater Storage Capacity of Part of Solano County. CA. US,” Geological Survey Water Supply Paper 1464, 1960.
  15. C. V. Theis, “Estimating Transmissivity of Water—Table Aquifer from Specific Capacity of a Well,” Geological Survey Water Supply Paper 1536-I, 1963, pp. 332-336.
  16. M. Razack and D. Hantley, “Assessing Transmissivity from Specific Capacity in a Large and Heterogeneous Alluvial Aquifer,” Ground Water, Vol. 29, No. 6, 1991, pp. 856-861. doi:10.1111/j.1745-6584.1991.tb00572.x
  17. D. Hantley, R. Nommensen and D. Steffey, “The Use of Specific Capacity to Assess Transmissivity in FracturedRock Aquifers,” Ground Water, Vol. 30, No. 3, 1992, pp. 396-402. doi:10.1111/j.1745-6584.1992.tb02008.x