CALCULATION OF LAND SUBSIDENCE CAUSED BY GROUNDWATER ABSTRACTION IN MỚ ĐÁ MINERAL WATER SOURCE, HOÀ BÌNH PROVINCE

NGUYỄN PHÚ DUYÊN

 Institute of Geosciences, VAST, 84 Chùa Láng Street, Hà Nội

Abstract: The calculation of land subsidence caused by groundwater exploitation in Mớ Đá mineral water source (Kim Bôi District, Hoà Bình Province) is carried out after Terzaghi’s theory. The value of land subsidence is relatively low and evaluated in range of from 0 to 23 cm. The value of land subsidence does not considerably affect the construction work in the study region.


I. INTRODUCTION

In recent years, it has been becoming clear that the extensive exploitation of groundwater resource has been bringing with it an undesired environmental side effect. At many localities in the world, groundwater pumpage from unconsolidated aquifer-aquitard system has been accompanied by significant land subsidence (Osaka: maximal land subsidence: 3 m, area of land subsidence: 190 km2, duration: 40 years; Tokyo: maximal land subsidence: 4 m, area of land subsidence: 190 km2, duration: 50 years; Mexico City: maximal land subsidence: 9 m, area of land subsidence: 130 km2, duration: 32 years; Santaclara (California): maximal land subsidence: 4 m, area of land subsidence: 650 km2, duration: 50 years; Las Vegas (Nevada): maximal land subsidence: 1 m, area of land subsidence: 500 km2, duration: 28 years; Water Supply factory (Pháp Vân, Hà Nội): maximal land subsidence 228.6 mm, from 1988 to 1994.

Pumping is from sand and gravel aquifers and the pumping always is accompanied by vertical leakage from the adjacent aquitards consisting of clay, which is composed of high-compressibility material. As a consequence, the pumping leads to compaction or consolidation. Consolidation is the gradual reduction in volume of a fully saturated soil, often of low permeability, due to drainage (for example: groundwater abstraction) of some of the pore water. The consolidation process continues until the original pore water pressure has completely disappeared, where it is replaced by an equivalent amount of intergranular stress in order to keep the total stress constant. The increase in intergranular stress causes the compaction.

In this paper, calculation of land subsidence in the study area is presented.

II. CALCULATION OF LAND SUBSIDENCE AFTER TERZAGHI’S THEORY

The simplest way of description is that of one-dimensional consolidation. Consolidation settlement is the vertical displacement of the surface corresponding to the volume change at any stage of the consolidation process. Consolidation settlement will result, for example, if a structure is built over a layer of saturated clay, or if the water table is lowered permanently in a stratum overlaying a clay layer.

The assumptions made in Terzaghi’s theory [5] of one-dimensional consolidation are:

- The soil is fully saturated;

- The soil particles are incompressible;

- Compression and flow are one-dimensional (vertical);

- Strains are small;

- The coefficient of permeability and coefficient of volume compressibility remain constant throughout the process;

- There is an unique relationship, independent in time, between void ratio and effective stress.

The main limitation of Terzaghi’s theory arises from the last assumption. Experimental results show that the relationship between void ratio and effective stress is not independent in time, i.e., there is hysteresis.

Suppose that the layer is subjected to a decrease in groundwater head h.

According to the Terzaghi’s theory, the vertical settlement L of the clay layer can be calculated by the formula:

 (1)

where: c - coefficient of consolidation.

 P - change in effective stress in the solid material (in kPa), but on other hand:

 P = *g*h

 - density of water (in kg/m3)

g = 10 m / s2 - gravity constant

h - drawdown (in metre)

L - compression of the considered layer (in metre)

L - thickness of the layer (in metre)

P - initial effective stress which is calculated by the following expression: weight of the overlaying formations + 1/2 own weight - water pressure

Drawdown (h) and coefficient of consolidation (c) are used according to Geological Party 47 [4].

III. PUMPING ACTIVITY IN THE STUDY AREA

Mớ Đá mineral water source is situated in Vĩnh Đồng and Hạ Bì communes, in the northwest of Kim Bôi townlet, stretching along NW-SE direction. Area of the source is about 9 km2. A net of pumping wells was put into operation, withdrawing mineral water for treatment, domestic purposes, etc.. Most of the wells are concentrated around the Kim Bôi Rest Home, consisting of: LK.14, LK.12, LK.13, LK.7, LK.11, LK.8 and LK.5B. Besides the central part, there are LK.3 and LK.6 in the NW and LK.4 in the north of the Mớ Đá mineral source. Most of the wells pump mineral groundwater from Cẩm Thuỷ (P3 ct) water-bearing formation, at a depth from 140 to 150 m. The water head in these wells rises 10 m above surface. At the outcrops ĐL.133 and ĐL.145, mineral groundwater of Cò Nòi (T1 cn) Formation was revealed. According to the Geological Party 47, just flow discharge of the ĐL.133 (4.86 l/s) can satisfy the total demand of mineral water for local use and bottling (about 110 m3/day).

Pumping activity of the well LK.7 gives best quality of the water, because there is a filter clay layer above, that can prevents pollution, but after long time of exploitation can causes land subsidence. It is clear that, the land subsidence can occur in Quaternary deposits in a stratigraphical column of the well. Thickness of the Quaternary deposits is 25 m [4], in which thickness of upper part is 7 m, mainly consisting of clay, sandy clay, clayish sand; the lower part is 18 m of thickness, consisting of spilite boulder, pebble, gravel, etc..

IV. CALCULATION OF LAND SUBSIDENCE AT SOME PUMPING WELLS

1. Calculation of land subsidence in the well LK.7

Parameters of the wells LK.7, LK.3 and LK.9 are illustrated below (Fig.1).

As mentioned above, Quaternary deposits of the well LK.7 is 25 m thick and consists of 7 m of clay layer in the upper part and 18 m of spilite boulder, pebble, gravel, etc. in the lower part.

- Land subsidence of 7 m of clay layer in the upper part:


Figure 1. Stratigraphical column of the wells LK.7, LK.3 and LK.9

Substituting L, c,  P, P into (1), one gets:

 (2)


Where: c- coefficient of consolidation of the clay in the study area, determined according to Geological Party 47. Value of the coefficient is in the range of from 0.017 to 0.08, depending on organic composition and “c” is taken as 0.035 in average.

Density of the clay is 1650 kg/m3.

Land subsidence of upper 7 m clay layer is approximately 21 cm.

- Land subsidence of 18 m of spilite, pebble, gravel, etc. in the lower part:

Coefficient of consolidation of the spilite, pebble, gravel, etc. is taken from literatures and is  0.006.


 

Substituting L, c,  P, P into (1), one gets:

 (3)


Land subsidence of lower 18 m spilite, pebble, gravel, etc. layer is approximately 1,6 cm.

Total land subsidence in LK.7 is about 23 cm.

2. Land subsidence in well LK.3

The well LK.3 is artesian, so that the water overflows into land surface. That is why there is no drawdown and no subsidence as well.

3. Land subsidence in well LK.9

According to Fig.1, Quaternary deposits of the well LK.9 is 20 m thick and consists of two layers. The upper one is 5 m thick and mainly consits of clay. The other one is 15 m thick and consists of spilite, pebble, gravel, etc..

Total land subsidence in the LK.9 is determined similarly as in LK.7 and equal 19 cm.

Distribution of land subsidence isolines of the study area is presented in Fig.2.


Figure 2. Distribution of land subsidence isolines
Mớ Đá mineral water source, Kim Bôi District, Hoà Bình Province


 V. CONCLUSION

Value of land subsidence caused by groundwater abstraction in the study area depends on geological conditions and drawdown. The calculated value varies in the range of from 0 to 23 cm. Consequence of the land subsidence does not considerably affect construction work in the area.

REFERENCES

1. De Smedt F., Nguyen Phu Duyen, Kapongo F., 1995. Land subsidence of the Jabotabeek Basin, Indonesia. Intergrated Project Design. Free University of Brussel, Belgium.

2. Lê Quý An, Nguyễn Công Mẫn, Nguyễn Văn Quý, 1997. Soil mechanics. High and Secondary School Publisher, Hà Nội.

3. Lomtadze V. D, 1978. Engineering Geology-Petrology. High and Secondary School Publisher, Hà Nội.

4. Phạm Văn Công (Editor), 1989. Preliminary geological survey report on Mớ Đá mineral water source, Hà Tây Province.

5. Terzaghi K. and R. B. Peck, 1967. Soil mechanics in an engineering practice. II edition, John Wiley & Sons, New York.