Abridged English version
1 Introduction
Hard rock aquifers generally occupy the first 100 metres of crystalline rocks from the top [3,14] that have been subjected to weathering processes [18,19]. The classical weathering profile comprises the following layers, which have specific hydrodynamic properties, from top to bottom (Fig. 1):
- – a weathering cover, with a thickness of 0 to 5 m; this layer has high porosity and low permeability. When saturated, it constitutes the reserve of the aquifer;
- – fissured zone, constituted by fractured hard rock, with a depth-decreasing density of fractures; these fractures result from the weathering process itself [18,19]: the weathering of micaceous minerals induces expansion and cracks in the rock (in granite, for example, the fractures are parallel to the topography during the weathering phase); this layer provides the transmissive function in the aquifer;
- – a solid basement, permeable only locally, where deep tectonic fractures are present.
When the spatial distribution of these layers with their hydrodynamic properties is understood, methodologies can be developed both for the assessment of groundwater resources and modelling of groundwater flow at the catchment scale [15,16,19].
The objective of this paper is to determine the coherence between geological observations and hydrodynamic properties of the aquifer, especially the anisotropy of permeability in the fissured zone.
2 Weathering profile
The Maheshwaram watershed located 30 km south of Hyderabad (India) is the main study area of the Indo-French Centre for Groundwater Research (French Geological Survey/National Geophysical Research Institute). Granites constitute the basement of the whole basin and the weathering profiles are easily observable through many dug wells used by the farmers for irrigation (Fig. 2). Under a few decimetres of red soils, the weathered layers with a thickness of 0 to 5 m cover the fissured zone. A high density of horizontal fractures is observed in the fissured zone, as illustrated by the picture in Fig. 2. Vertical fractures with a tectonic origin are also present. Due to the strong exploitation of groundwater resources, water levels are far below ground level and the weathered layers are dry, while only the fissured zone is saturated. This allows us to specifically test this layer through pumping tests.
3 Vertical anisotropy of the fissured zone
Two pumping tests were carried out at constant discharge rates. Drawdowns in observation wells are interpreted by the Neuman [12] method, while drawdowns in the pumping wells are analysed using the theory of the horizontal fracture developed by Gringarten [5].
3.1 Neuman method at observation wells
The theory initially developed by Boulton [1] to interpret some unusual drawdown curves obtained in the observation wells (see Fig. 3) takes into account the notion of “delayed yield from storage in unconfined aquifers” [2]. It was improved by Neuman [11,12], who developed an analytical solution adapted to anisotropic unconfined aquifers where Kr is the radial permeability parallel to the aquifer extension and Kz is the vertical permeability. The Neuman solution, available as an abacus, gives reduced drawdowns in an observation well located at a radial distance r from the pumping well, , as a function of:
- – reduced time for type-A curves;
- – reduced time for type-B curves;
The application of this method (Table 1) to the observation wells IFP-1/1, IFP-1/2 and IFP-9/1 leads to the evaluation of transmissivities, storage coefficients (S) and specific yields (Sy). The results are consistent with other observations and methods of interpretation. The determination of KD requires the knowledge of the aquifer thickness b. Flowmeter measurements during injection tests in eight wells of the catchment have shown that the solid basement does not contain any conductive fracture (Fig. 4). Thus, the top of the solid basement was chosen as the bottom of the aquifer. With this assumption, the anisotropy ratio and the permeability can be computed and are given in Table 2. The results show a horizontal permeability 2 to 30 times higher than the vertical one. This result is consistent with the observation of horizontal fractures in dug wells.
Transmissivités et paramètres d'emmagasinement obtenus par ajustement des rabattements (TA : transmissivité obtenue par ajustement sur la courbe de type A, TB : transmissivité obtenue par ajustement sur la courbe de type B, TAB : moyenne géométrique de TA et TB)
Transmissivity and storage parameters obtained by fitting of drawdown (TA : transmissivity obtained by fitting on type A curve, TB : transmissivity obtained by fitting on type B curve, TAB : average of TA and TB)
Puits d'observation | Puits de pompage | r (m) | TA () | TB () | TA/TB (–) | TAB () | S (–) | Sy (–) |
IFP-1/1 | IFP-1 | 28 | 1,76×10−5 | 1,96×10−5 | 0,90 | 1,86×10−5 | 7,0×10−5 | 1,7 ×10−3 |
IFP-1/2 | IFP-1 | 27,5 | 1,71×10−5 | 1,76×10−5 | 0,97 | 1,74×10−5 | 3,7×10−5 | 1,5 ×10−3 |
IFP-9/1 | IFP-9 | 30,7 | 5,55×10−4 | 7,65×10−4 | 0,73 | 6,51×10−4 | 7,1×10−4 | 3,4 ×10−3 |
Perméabilités et rapport d'anisotropie déterminés aux puits d'observation par la méthode de Neuman
Permeability and anisotropy ratio at observation wells using Neuman's method
Puits d'observation | β (–) | r (m) | b (m) | Kr () | Kz () | KD (–) | 1/KD (–) |
IFP-1/1 | 1,00 | 28 | 21,8 | 8,5×10−7 | 5,2×10−7 | 0,606 | 1,7 |
IFP-1/2 | 0,20 | 27,5 | 21,8 | 8,0×10−7 | 1,0×10−7 | 0,126 | 8,0 |
IFP-9/1 | 0,60 | 30,7 | 7,3 | 9,0×10−5 | 3,0×10−6 | 0,034 | 29,5 |
3.2 Gringarten method at pumping wells
Vertical flowmeter profiles in wells IFP-1 and IFP-9 (Fig. 4) show that a few fractures are conductive, respectively three in F1/1, F1/2 and F1/3 and one in F9/1. Actually, at IFP-1, only the deepest one was saturated during the whole pumping test. Similarly, the only conductive fracture crossed by IFP-9 was also saturated during the pumping test. Moreover, the analysis by Neuman's method shows the existence of anisotropy due to horizontal fractures. Thus, the method developed by Gringarten [5,6] for a vertical well crossing a horizontal fracture in an anisotropic aquifer, applicable to the pumping well, is well adapted to the hydrogeological context of the IFP-1 and IFP-9 wells.
The interpretation is done through the fitting of observed drawdowns on theoretical curves of an abacus [5] giving the reduced drawdown in a pumping well for a fracture located at the centre of the aquifer, as a function of the reduced time. Knowing the geometry of the case (i.e. the thickness H of the aquifer), assuming the distance zf between the bottom of the aquifer and the fracture to be and using the value determined with Neuman's method for Ss, the hydrodynamic properties of the aquifer are evaluated (Table 3): Kr the horizontal permeability, Kz the vertical permeability and the radius rf of the horizontal fracture. With this method, an anisotropy of the aquifer is also shown in the pumping wells, with a horizontal/vertical ratio on the order of 5 to 20.
Perméabilités, rapport d'anisotropie et rayon de la fracture horizontale déterminés aux puits de pompage par la méthode de Gringarten ( ∗: Ss est la moyenne des coefficients d'emmagasinement spécifique déterminés aux puits d'observation par la méthode de Neuman)
Permeability, anisotropy ratio and radius of the horizontal fracture determined at the pumping wells using Gringarten's method ( ∗ : Ss is the average of specific storage coefficient determined for each site using Neuman's method at observation wells)
Paramètres connus | Paramètres déterminés par l'ajustement | |||||||||
Puits | H (m) | Fissure | zf (m) | zf/H (–) | (m−1) | HDG (–) | rf (m) | Kr () | Kz () | Kr/Kz (–) |
IFP-1 | 21,8 | F1/3 | 6,0 | 0,28 | 2,5×10−6 | 3 | 30,4 | 3,5×10−6 | 1,9×10−7 | 18,7 |
IFP-9 | 7,3 | F9/1 | 4,0 | 0,55 | 9,7×10−5 | 5 | 3,4 | 2,9×10−5 | 5,2×10−6 | 5,5 |
4 Conclusion
The geological observations on the granites of the Maheshwaram watershed (India) confirm the existence, as in many other areas of the World, of a high density of horizontal fractures in the fissured zone of the weathering profile. Measurements made with a flowmeter confirm that only a few of these fractures are conductive. The interpretation of pumping tests on several wells with observation wells regularly shows the existence of a vertical anisotropy of permeability: the horizontal permeability is 2 to 30 times higher than the vertical permeability. These results confirm the major influence of fissures with a weathering origin on the hydraulic parameters of hard-rock aquifers, whereas fissures with a tectonic origin have a less important role. The application of the Gringarten theory also makes it possible to determine the radius of horizontal fractures crossing the pumping wells.
1 Introduction
Les roches de socle (roches plutoniques et métamorphiques) ne sont aquifères qu'au sein de leur frange superficielle, sur les 100 premiers mètres environ sous la surface du sol [3,14], qui a été soumise à des processus d'altération supergène [18,19].
Le profil d'altération type comprend, de haut en bas, les horizons suivants, qui possèdent des propriétés hydrodynamiques spécifiques (Fig. 1).
- – Les altérites (ou saprolites) couvrent la roche « saine » sur une épaisseur variable (de 0, là où cet horizon est érodé, à plusieurs dizaines de mètres). Du fait de sa dominante argileuse (argilo-sableuse dans le cas de roches riches en quartz), ce compartiment se caractérise par une porosité d'interstices significative (au sein de granites porphyriques, la porosité efficace peut dépasser 5 %) et par une perméabilité relativement faible. Lorsqu'il est saturé en eau, ce compartiment assure une fonction capacitive de stockage des eaux souterraines au sein de l'aquifère.
- – L'horizon fissuré est constitué de roche dure, affectée, sur plusieurs dizaines de mètres d'épaisseur, par la présence de fissures, de fréquence décroissante avec la profondeur. Il est maintenant démontré que ces fissures, dont l'origine a longtemps été associée à la « décompression superficielle », résultent, elles aussi, du processus d'altération [18,19]. Leur genèse est due principalement à l'altération des minéraux phylliteux (biotite en particulier), dont le gonflement provoque cette fissuration. Au sein de roches à texture isotrope (granites, par exemple), la fissuration s'exprime préférentiellement le long de plans sub-horizontaux, perpendiculaires à la contrainte minimale. Par conséquent, ces plans sont généralement parallèles à la morphologie contemporaine de l'altération. Au sein de roches anisotropes non plissées (granites foliés, gneiss), la fissuration reste planaire, mais son orientation et son intensité sont déterminées par l'angle formé par la foliation et la surface topographique. Au sein de roches intensément plissées (schistes, micaschistes), la fissuration a une orientation généralement aléatoire, même à l'échelle de l'affleurement. Au sommet de l'horizon fissuré, la fréquence des microfissures devient voisine de la taille du grain de la roche (de l'ordre du millimètre), marquant la transition avec les altérites sus-jacentes. Une conséquence importante de ces modalités de genèse est que l'horizon fissuré du granite est parallèle aux paléosurfaces d'altération. Cet horizon assure la fonction transmissive de l'aquifère et est capté par la plupart des forages réalisés en région de socle.
Le substratum sain n'est perméable que très localement, là où il est affecté par des fractures d'origine tectonique, le plus souvent sub-verticales.
Ces (paléo)profils d'altération sont connus dans de très nombreuses régions du monde [13] : Amérique du Nord, Amérique du Sud, Afrique, Inde, Chine, Corée, Japon, Australie, Europe, y compris Europe du Nord [9,10,15,17], etc.
La compréhension, d'une part, des modalités de genèse de ces profils d'altération (conditions géomorphologiques, climatiques, géodynamiques, de durée) et, d'autre part, des conditions de leur éventuel démantèlement ultérieur par l'érosion, a permis de mettre au point une méthode de cartographie de l'épaisseur résiduelle de ces différents compartiments aquifères, qu'il est possible de mettre en œuvre à l'échelle de bassins versants de plusieurs dizaines à plusieurs milliers de kilomètres carrés de superficie [15,16,19]. Les cartes obtenues permettent de caractériser la variabilité spatiale des propriétés hydrogéologiques des roches de socle, à cette même échelle, et de réaliser, par exemple, des cartes prédictives des potentialités en eau souterraine d'une région donnée [8,19].
La disponibilité de ces éléments cartographiques ouvre également des perspectives en terme de mise en œuvre, au sein de ces aquifères spécifiques, des modèles de gestion de la ressource en eau souterraine classiquement utilisés au sein des milieux poreux (modèles « maillés » [7]). Les écoulements souterrains peuvent être modélisés au sein de chacune des couches dont la géométrie est ainsi définie, à l'échelle du bassin versant, seule pertinente pour la gestion des ressources en eau.
Dans ce cadre, les altérites peuvent aisément être assimilées à un milieu poreux. Des travaux sont, en revanche, nécessaires pour caractériser les propriétés hydrodynamiques de l'horizon fissuré, dans la perspective de prendre en compte l'hétérogénéité et l'anisotropie de ses propriétés, de manière pertinente, à l'échelle de la maille d'un modèle hydrogéologique. Le présent travail s'attache à la caractérisation de l'anisotropie de perméabilité au sein de la zone fissurée.
2 Caractérisation géologique du profil d'altération
Dans la perspective de la mise au point de la méthodologie de modélisation des aquifères de socle, présentée ci-dessus, des travaux de recherche sont mis en œuvre conjointement par le BRGM et le National Geophysical Research Institute, dans le cadre du Centre franco-indien de recherche sur les eaux souterraines (CEFIRES) d'Hyderabad (Inde).
Le bassin versant étudié autour du village de Maheshwaram, situé à 30 km au sud d'Hyderabad, d'une soixantaine de kilomètres carrés de superficie, est majoritairement composé de granites archéens à texture isotrope : granites porphyroı̈des à biotite et leucogranites.
Les profils d'altération y sont aisément observables grâce aux nombreux puits ouverts, creusés par les agriculteurs au début du siècle pour l'irrigation de leur terre (Fig. 2). Ces profils sont généralement tronqués par l'érosion, l'horizon fissuré apparaissant sous les altérites, dont l'épaisseur dépasse rarement 5 m. Une forte densité de fissures sub-horizontales est observée dans la zone fissurée, comme l'illustre la photographie de la Fig. 2. Des diaclases verticales, dont l'origine pourrait être tectonique, sont également présentes.
Du fait de la surexploitation des eaux souterraines, les niveaux piézométriques sont profonds sur le site de Maheshwaram. Par conséquent, sur la majorité du bassin versant, seul l'horizon fissuré de l'aquifère est saturé en eau, les altérites étant dénoyées. Cette particularité a été mise à profit pour réaliser des pompages d'essai destinés à caractériser les propriétés hydrodynamiques du seul horizon fissuré et, en particulier, pour déterminer si des analogies existent entre la structure géologique de cette couche (fissuration à dominante sub-horizontale) et ses propriétés hydrodynamiques.
3 Anisotropie verticale de la perméabilité au sein de la zone fissurée
Deux essais de pompage de moyenne durée (17 et 12 h), à débit constant, ont été réalisés, avec mesure des rabattements aux puits d'observation et de pompage, sur deux forages expérimentaux réalisés sur le bassin versant. Les rabattements aux puits d'observation ont été interprétés par la méthode de Neuman [12], tandis que les rabattements aux puits ont été analysés au moyen de la méthode de la fracture horizontale de Gringarten [5,6].
3.1 Méthode de Neuman aux puits d'observation
À titre d'exemple, sur un diagramme bi-logarithmique (Fig. 3), les courbes de rabattements observés aux puits d'observation IFP-1/1 et IFP-1/2, lors de l'essai de pompage au puits IFP-1 (voir Tableau 1 pour les détails de l'essai), montrent une forme complexe, difficilement interprétable par les méthodes classiques. Les courbes de rabattement sont composées de trois segments : le premier, pour les temps courts, avec des pentes relativement fortes, est suivi d'une période intermédiaire, durant laquelle une stabilisation des niveaux semble prendre place ; une troisième phase caractérise ensuite les temps plus longs et montre un nouvel accroissement des pentes.
La théorie initialement développée par Boulton [1] pour interpréter ce type de courbe de rabattement fait appel à la notion de « débit retardé issu de l'effet capacitif d'aquifères libres » – delayed yield from storage in unconfined aquifers [2]. Elle a été améliorée par Neuman [11,12] qui a développé une solution analytique de l'équation de la diffusivité adaptée aux aquifères libres anisotropes.
Neuman [12] considère que l'aquifère est libre et d'extension infinie. Il repose sur une couche horizontale imperméable. L'aquifère est homogène, mais anisotrope, et ses perméabilités principales Kr et Kz sont orientées parallèlement aux axes de coordonnées (r,z). Lorsqu'un débit constant Q est pompé dans un puits complet, l'eau provient, d'une part, de l'emmagasinement au sein de l'aquifère et, d'autre part, du drainage gravitaire à sa surface libre. La solution de Neuman, donnée sous forme d'abaque, exprime les rabattements réduits dans un piézomètre situé à une distance radiale r du puits de pompage, en fonction :
- – du temps réduit pour les courbes de type A ;
- – du temps réduit pour les courbes de type B ;
Les données expérimentales recueillies lors de l'essai de pompage à IFP-1 permettent des ajustements sur l'abaque de Neuman (Fig. 3) avec les mêmes valeurs du paramètre β, tant pour les courbes A que pour les courbes B, comme le requiert cette technique d'interprétation.
Les paramètres hydrauliques obtenus au moyen de ces ajustements sont donnés au sein du Tableau 1, pour les essais de pompage à IFP-1 et IFP-9.
Les valeurs très similaires obtenues dans chacun des puits d'observation, pour TA et TB, mettent en évidence la cohérence de l'interprétation de cet essai de pompage au moyen de la méthode de Neuman.
Les coefficients d'emmagasinement (S) présentent des valeurs réalistes pour ce type de milieu rocheux. Il en est de même de la porosité de drainage, dont les valeurs sont comparables à celles obtenues, tant de manière indirecte par mesures RMP [18,19] que dans un puits proche, par ajustement des fluctuations de niveaux piézométriques avec les précipitations au moyen d'un modèle global [4].
Le rapport d'anisotropie de l'aquifère (KD) dépend de b, r et β. Alors que les deux derniers paramètres sont bien connus, le premier est incertain dans ce type d'aquifère fracturé. Des essais au micro-moulinet, effectués dans huit puits du bassin versant, montrent que les fractures hydrauliquement actives sont situées dans l'horizon fissuré, c'est-à-dire dans l'horizon situé au-dessus du granite sain et compact. C'est notamment le cas dans les puits de pompage IFP-1 et IFP-9 (Fig. 4). Le sommet du granite sain a donc été retenu comme plancher de la nappe pour la détermination de b. Comme pour les méthodes classiques, l'incertitude sur la valeur de r au forage pompé rend néanmoins difficile toute interprétation quantitative des rabattements dans le puits de pompage.
Les résultats de l'interprétation, présentés dans le Tableau 2, indiquent une anisotropie du tenseur de perméabilité, en accord avec les observations géologiques : la perméabilité horizontale est systématiquement supérieure à la perméabilité verticale, d'un facteur allant de 2 à 30, ce qui peut être mis en regard de la structure de l'horizon fissuré, montrant majoritairement des fissures horizontales.
3.2 Méthode de Gringarten aux puits de pompage
Les profils au micro-moulinet (Fig. 4) suggèrent que peu de fractures sont actives dans les puits IFP-1 et IFP-9, respectivement trois (F1/1, F1/2 et F1/3) et une (F9/1). Lors de l'essai de pompage à IFP-1, étant donné le rabattement au puits et la faible valeur des pertes de charges quadratiques, seule la fracture profonde située à environ 31,5 m de profondeur (F1/3) n'a pas été dénoyée et est restée active durant tout l'essai. De façon similaire, la seule fracture active intersectée par IFP-9 (F9/1), localisée à 29 m de profondeur, est restée noyée durant l'essai de pompage. L'analyse par la méthode de Neuman effectuée ci-dessus suggère une anisotropie du milieu qui résulte de la présence de fractures horizontales. La méthode développée par Gringarten [6], pour un puits vertical traversant une fissure horizontale unique en milieu de perméabilité anisotrope, applicable au puits de pompage, est donc adaptée à la configuration hydrogéologique des puits IFP-1 et IFP-9.
La complexité de la solution analytique nécessite l'utilisation d'un abaque [5], dont les courbes expriment le rabattement réduit sDG dans le puits en fonction du temps réduit tDG pour différentes configurations géométriques représentées par le paramètre HDG, quand la fracture est au centre de l'aquifère (zf/H=0,5), avec
Les ajustements des rabattements observés sur l'abaque de Gringarten conduisent à des valeurs élevées de HDG, suggérant une forte anisotropie de perméabilité. Connaissant les paramètres géométriques H et zf, ainsi que le coefficient d'emmagasinement obtenu précédemment par la méthode de Neuman (moyenne de 2,5×10−6 à IFP-1/1 et IFP-1/2 et 9,7×10−5 à IFP-9), il est possible d'estimer le rayon de la fracture ainsi que les perméabilités horizontales et verticales de l'aquifère (Tableau 3).
Dans les deux cas, l'anisotropie verticale est évidente, avec des facteurs de l'ordre de 5 à 20, tandis que les rayons de fracture estimés (de 3 à 30 m) sont cohérents avec les observations de terrain.
4 Conclusion
Les observations géologiques effectuées sur les granites du bassin versant de Maheshwaram (Inde) confirment l'existence, comme dans de très nombreuses autres régions du monde, d'une forte densité de fractures horizontales au sein de l'horizon fissuré du profil d'altération. Elles montrent aussi la présence de quelques diaclases sub-verticales. Les mesures au micro-moulinet confirment qu'une faible fraction seulement de ces structures est hydrauliquement active.
L'interprétation de pompages d'essai sur plusieurs forages dotés de piézomètres montre systématiquement l'existence d'une anisotropie verticale de perméabilité : la perméabilité horizontale est de 2 à 30 fois supérieure à la perméabilité verticale. L'utilisation de la solution analytique de Gringarten permet en outre de déterminer le rayon des fractures horizontales intersectant les puits de pompage.
Ces résultats confirment le rôle majeur, à l'échelle du bassin versant, de la fissuration liée à l'altération, horizontale dans le cas de granites, dans l'acquisition de la perméabilité des roches de socle ; la fracturation de type tectonique joue un rôle subordonné à cette échelle.
Remerciements
Ce travail a été réalisé par les chercheurs du Centre franco-indien (BRGM–NGRI) de recherche sur les eaux souterraines. Les auteurs remercient le ministère français des Affaires étrangères et le Service de coopération et d'action culturelle de l'ambassade de France à New Delhi, qui a financé les voyages d'experts français en Inde. Le Centre franco-indien de recherche sur les eaux souterraines a par ailleurs bénéficié d'un financement du CNRS dans le cadre du Programme ACI « Eau et Environnement » et d'un financement du CEFIPRA (Centre franco-indien pour la promotion de la recherche avancée).