République Algérienne Démocratique et Populaire

Ministère de l’Enseignement Supérieur et de la Recherche Scientifique

Université BATNA

Faculté de technologie
Département d’Hydraulique

THÈSE
En vue de l’obtention du diplôme de

DOCTORAT EN SCIENCES

Option : Hydraulique

Thème

Caractérisation des paramètres hydrodynamiques de l’aquifère de

Tadjnant –Chelghoum Laid et impact de la pollution des eaux de
surface sur les eaux souterraines

Présenté par

KHEDIDJA ABDELHAMID
(Ingénieur d’État en Hydrogéologie. Magister en géologie appliquée)

Soutenue publiquement le

Devant le jury composé de :

KALLA Mehdi Président Professeur Université de Batna

BOUDOUKA Abderrahmane Rapporteur Professeur Université de Batna
FEHDI Chemseddine Examinateur Professeur Université de Tébessa
HOUICHI Larbi Examinateur Professeur Université de Batna
CHABOUR Nabil Examinateur Maitre de conférences Université de Constantine
 RESUME

Caractérisation des paramètres hydrodynamiques de l’aquifère de Tadjenant –

Chelghoum Laid et impact de la pollution des eaux de surface sur les eaux
souterraines

La région explorée fait partie de la zone occidentale du grand bassin versant de Kébir
Rhumel situé à l’Est Algérien. Du point de vue géologique la région étudiée est
caractérisée par une formation de type sédimentaire dont la prédominance des
affleurements correspond à des formations Mio-Plio-Quaternaires.
La présente étude a pour objectif essentiel d’étudier le comportement hydrodynamique de
la nappe alluviale et de déterminer l’impact de la pollution naturelle ou anthropique sur
la qualité physico chimique des eaux de la nappe superficielle de la zone Tadjenanet –
Chelghoum Laid.
Les résultats du modèle hydrodynamique en régime transitoire étalée sur une période
d’observation de 20 ans, montrent un bilan hydrogéologique de la nappe plus au moins
équilibré traduit par un faible rabattement atteignant les 6 m au maximum.
L’étude du chimisme de la nappe alluviale a permis de dire que le comportement
hydrochimique est caractérisé par une grande variabilité spatiale. La minéralisation
augmente progressivement de l'Ouest vers Est coïncidant avec le sens d’écoulement
principal des eaux souterraines. La concentration très élevée en nitrates est en liaison avec
l’activité agricole dans la région.
L’analyse des cartes de vulnérabilité a mis en évidence une zone de forte vulnérabilité
au centre la plaine, alors que les zones à moyenne vulnérabilité se localisent à l’Est sur la
périphérie de Oued Rhumel, le reste du terrain a une faible vulnérabilité.

Mots clés : Hydrodynamique, nappe, alluvion, impact, qualité, eau.

 ABSTRACT

Characterization of hydrodynamic parameters of the aquifer of Tadjenant-

Chelghoum Laid and impact of pollution of surface water on groundwater

The region explored is part of the western area of the great watershed Kebir Rhumel
located in eastern Algeria. From a geological point of view the study area is characterized
by formation of sedimentary outcrops which predominantly corresponds to Mio- Plio-
Quaternaryformations.
The main objective of this study is to study the hydrodynamic behavior of the alluvial
aquifer and determine the impact of natural and anthropogenic pollution on the physico-
chemical water quality of surface water in the area of Tadjenanet - Laid Chelghoum.
The results of the hydrodynamic model transient over a period of observation of 20 years
, show a water balance of the water more or less balanced results in low drawdown up to
6m. .
The study of the chemistry of the alluvial aquifer was concluded that the hydrochemical
behavior is characterized by a large spatial variability. Mineralization increases gradually
from west to east coinciding with the main direction of groundwater flow. Very high
nitrate concentration is connected with agricultural activity in the region.

The analysis of vulnerability maps showed an area of high vulnerability to center the
plain, while the average vulnerability areas are located to the east on the outskirts of Wadi
Rhumel, the rest of the field at a low vulnerability.

Keywords: Hydrodynamics, nappe, alluvium, impact, quality, water.

 ‫ملـــــخـــص‬

‫دراسة هيروديناميكية للمكمن المائى لمنطقة تاجنانت – شلغوم العيد و أثر التلوث‬
‫المياه السطحية على المياه الجوفية (شرق الجزائر)‬
‫كبير الرمال‬ ‫تقع المنطقة المدروسة في الجهة الغربية لحوض التصريف‬
‫‪ (Kebir (Rhumel‬و الذي يقع بالجهة الشرقية للجزا ئر‪.‬‬
‫من جهة نظر جيولوجية تتميز منطقة الدراسة بتكوينات روسوبية اغلبها عبارة‬
‫عن نتوءات التي ترجع حقبتها الى العصر الرباعي البليوساني والميوساني‪ ,‬أين تحيط‬
‫بالمنطقة تكوينات روسوبية تعود الى العصر االيوساني و الطباشيري‪.‬‬
‫الهدف الرئيسي من هذه الدراسة هو دراسة السلوك الهيدروديناميكي للمكمن‬
‫المائي الجوفي الغريني وتحديد أثر التلوث الطبيعي والبشري على نوعية المياه الفيزيائية‬
‫والكيميائية للمكمن المائي الجوفي السطحي للمنطقة الممتدة من تاجنانت الى شلغوم‬
‫العيد‪.‬‬
‫نتائج نموذج الهيدروديناميكي المطبق على مدى فترة المراقبة المقدرة ب ‪02‬‬
‫عاما‪ ،‬تبين وجود توازن الى حد ما في حصيلة الماء المغذي للمكمن و المستغل مع‬
‫تسجيل انخفاض في منسوب الماء الجوفي للمكمن يصل إلى ‪ 6‬أمتار كحد أقصى‪.‬‬
‫الدراسة الكيميائية لمياه المكمن الجوفي الغريني خلصت الى أن السلوك الكيميائي للمياه‬
‫المدروسة يتميز ب تغير كبير من منطقة الى اخرى و هذا بالتباين الكبير المسجل على‬
‫التركيز للعناصر الكيميائية‪ .‬الدراسة سمحت بتسجيل زيادة في الملوحة تدريجيا من‬
‫الغرب إلى الشرق وتزامن ذلك مع االتجاه الرئيسي لتدفق المياه الجوفية ‪ .‬أيضا نشير‬
‫الى تواجد نسب عالية جدا في عنصر النترات والذي يرتبط أساسا مع النشاط الزراعي‬
‫في المنطقة‪.‬‬
‫تفحص الخرائط المنجزة والخاصة بمدى قابلية المكمن المائي للتلوث تظهر‬
‫بوجود هناك حدة لتلوث هذا االخير في الجهة الوسطى للسهل مع قابلية متوسطة في‬
‫الجهة الشرقية القريبة من وادي الرمال الى أضعف حدة في باقي المنطقة ‪.‬‬

‫الكلمات المفتاحية ‪ :‬ديناميكية المياه‪ ،‬مكمن مائي‪ ،‬غرين‪ ،‬أثر‪ ،‬نوعية‪ ،‬الماء‪.‬‬
 Dédicace

Je dédie cette thèse à

A la mémoire de mon très cher père et ma très chère mère

A ma femme et mes enfants

A mes frères et sœurs
A toute ma famille
 REMERCIEMENTS

Je tiens dans un premier temps à remercier le Professeur BOUDOUKHA

Abderrahmane, de l’université de Batna d'avoir accepté de m'encadrer et suivre mon
travail de près avec sa rigueur scientifique et ses conseils qui m’ont bien orienté dans
l'élaboration de cette thèse. Je le remercie surtout pour sa disponibilité malgré ses
responsabilités. Qu'il trouve ici l'expression de ma profonde gratitude.
Mes remerciements vont également à KALLA Mehdi, Professeur à l'université de Batna
d’avoir accepter de présider mon jury de soutenance, aux membres examinateurs :
HOUICHI Larbi ; Professeur à l'université de Batna, FEHDI Chemseddine ; Professeur à
l’université de Tébessa et CHABOUR Nabil, Maître de conférence (A) à l’université de
Constantine d’avoir bien voulu assister et juger cette thèse.
Mes remerciements s’adressent aux personnels de laboratoire de la station d’épuration de
Batna et de laboratoire de chimie de l’université de Batna.
Je tiens à remercier le directeur du Laboratoire d’hydrologie et de géochimie de
Strasbourg Philippe Ackerer qui m’a bien accueilli dans son laboratoire et m’a permis de
contacter des chercheurs dans le domaine de géochimie et la modélisation des
écoulements souterrains et de surfaces, tels que Mr Gérard Schafer.
J’ai l’agréable tâche de témoigner ma grande reconnaissance aux différents
services techniques auprès desquels j’ai bénéficié d’un aimable appui, notamment ceux
qui m’ont facilité l’accès aux données de mesures et à la documentation, en particulier :
- A.N.R.H. de Constantine (service d’hydrologie et hydrogéologie) ;
- Direction de l’Hydraulique de la Wilaya de Mila ;
- Subdivisions d’hydraulique de Tadjenanet et de Chelghoum Laid.
Ma reconnaissance à tous les propriétaires des points d’eau qui m’ont autorisé à utiliser
leurs puits et forages pour effectuer les prélèvements d’eau et les mesures
hydrodynamiques.
Et enfin, je tiens à remercier tous les amis qui m’ont aidé et encouragé pour réaliser cette
thèse, je cite spécialement mes collègues du Département des sciences de la terre et de
l’univers de l'université de Batna Derias T, Belkheri L ,Briniss N et Djenba S pour les
bons moments passés ensemble et les échanges scientifiques ainsi que Mr Baghiani B, de
département de biologie pour son aide pendant la campagne d’analyses chimique des
eaux .
Mes remerciements les plus sincères vont également à ma famille pour son soutien moral
et sa patience pendant les moments difficiles dans la réalisation de ce travail. Grâce à eux,
les moments de découragement furent rares et brefs. Qu'ils trouvent ici le témoignage de
ma profonde gratitude.
 TABLE DES MATIERES

INTRODUCTION GENERALE

GENERALITES ET ETAT DES CONNAISSANCES N°

Page
Introduction………..……………………………………………………………...………..
Problématique………………………………………………………………………….........
Objectif et méthodologie………………………………………………………….………...
Caractères physico-géographiques……………………………………………….…………
Situation géographique………………………………………………………….….……..
Orographie………………………………………………………………………..….........
Réseau hydrographique …………………………………………………………..………
Géomorphologie…………………………………………………………………..………
Relief ……………………………………………………………………….….……….
Plaine ……………………………………………………………………….….….........
Végétation……………………………………………………………………….………...

CHAPITRE I : CADRE GEOLOGIQUE

Introduction ………………………………………………………………………..….........
Description générale des diverses unités géologiques ………………………………...…...
Tertiaire et Quaternaire………………………………………………………….…...........
Eocène méridional………………………………………………………………………..
Eocène centrale …………………………………………………………………………..
Eocène inférieur…………………………………………………………………..……….
Crétacé supérieur et ultra-supérieur…………………………………………….…............
Crétacé inférieur et moyen…………………………………………………….….............
Le Trias…………………………………………………………………………..………..
Lithostratigraphie de la région d’étude……………………… ….………………….............
Formations peu ou pas tectonisées ………………………………………........................
Quaternaire……………………………………………………………………..……….
a. alluvions récentes…………………………………………………………………..………..
b. Terrasses, formation de pente et Quaternaire indéterminé……………………… …............
c. Glacis polygéniques…………………………………………………………………………
d. Villafranchien : croûtes calcaires…………………………………………………………...
e. Sebkha……………………………………………………………………………………….
f. Sols salés……………………………………………………………………….…................
g. Formations dunaires…………………………………………………………….......... …….
Mio – Pliocène continental……………………………………………………... ...........
 Formations tectonisées………………………………………………………….................
Unité supérieure à matériel Eocène……………………………………………..............
Nappe du Djebel Sattor et unités Eocènes à Nummulites et à Huîtres ………................
a. Lutétien à Huîtres…………………………………………………………………...............
b. Yprésien et Lutétien à débris………………………………………………………………..
c. Yprésien et Lutétien à Nummulites…………………………………………………………
d. Sénonien supérieur Lumachellique……………………………………………………....…
Nappe néritique Constantinoise……………………………………………………....…
a. Aptien calcaire………………………………………………………………….…………...
b. Barrémien carbonaté………………………………………………………………..……….
c. Néocomien calcaire………………………………………………………………………….
d. Jurassique carbonaté………………………………………………………………..……….
Ensemble allochtone Sud-Sétifien …………………………………………….………..
a. Miocène marin transgressif………………………………………………………………….
b. Barrémien……………………………………………………………………….…..............
Trias exotique ou extrusif…………………………………………………….…............
Analyse structurale………………………………………………………………..…...........
Ensemble Allochtone Sud-Sétifien……………………………………………………….……….
Unités du trias diapyrique.. ………………………..……………………………….….............
Unités du Crétacé-Eocène………………………………………………………………………
Unités Mio-Pliocène transgressives……………………………………………….............
Description géologiques à partir des logs de forage et données géophysiques……………..
Conclusion…………………………………………………………………………………..

CHAPITRE II : CADRE HYDROCLIMATIQUE

Introduction……………………………………………………………….……….. ………
Caractéristiques climatiques de la région d’étude………………………………………..…
Etude des paramètres climatiques………………………………………… ……………….
Précipitations………………………………………...…………………………………....
Coefficient pluviométrique………………………………………………………………..
Les températures…………………………………………………………………………..
Diagramme pluviothermique……………………………………………………………..
Indice d'aridité…………………………………………………………… ……………....
Indice d'Euverte…………………………………………………………………………...
Etablissement du bilan hydrologique ……………………………………………………….
Analyse des différents éléments du bilan………………………………………………....
Précipitations .…………………………………………………………………….…….
Evapotranspiration……………………………………………………………………....
Evapotranspiration potentielle (ETP)………………………………………………...
Les réserves facilement utilisables (R.F.U)……………………………………..…….
Evapotranspiration réelle……………………………………………………………...
Méthode de Thornthwaite………………………………………………………..….
Méthode de L. Turc………………………………………………………….……...
Méthode de A. Coutagne………………………………………………………..…..
 Estimation du ruissellement et de l'infiltration……………………………………….....
Analyse morphométrique……………………………………………………….……….….
Introduction……………………………………………………………………………......
Variation de l’altitude………………………………………………………………..…....
Forme du bassin……………………………………………………………………..…….
Densité de drainage……………………………………………………………………......
Profil en long………………………………………………………………………..…….
Conclusion………………………………………………………………………………......

CHAPITRE III : GEOMETRIE DU COMPLEXE AQUIFERE

Introduction ……………………………………………….……………………………..….
Situation de l’étude………………………………….………………………….…….…..…
Répartition des sondages électriques…………………………………….……………..…...
Interprétation………………………………………………………..………………………….
Examen de quelques familles de S.E…………………………………...………….....…….
Sondage électrique à l’extrême-Ouest, coupe A…………………………………....….
Sondage électrique étalon B3…………………………………………………..…….…
Cartes de résistivité……………………………………………………….…………....…
Carte de résistivité AB = 200 m ………………………………………….............…....
Carte de résistivité AB = m………………………………………..………....…....
Carte de résistivité AB = 3 m…………………………………………….…................
Interprétation des Coupes géo électriques………………………………..........………...
Carte du substratum………………………………………………………………………….
Carte de résistance transversale …………………………………………………………....
Conséquences hydrogéologique et conclusion……………………………………….……..

CHAPITRE IV : HYDROGEOLOGIE

Introduction…………………………………………………………………….………...…
Aquifère du Mio-Plio- Quaternaire…………………………………………….………...…
Description ……………………………………………………………………...………...
Géométrie de l’aquifère Mio-Plio-Quaternaire ………………………………………..…
Piézométrie……………………………………………………………………………….…
Inventaire des points d’eau………………………………………………….…..…….…...
Conditions aux limites ………………………………………………………..………..…
Interprétation de la carte piézométrique (Mai 7)………………………………..………
Interprétation de la carte piézométrique (Sep 7)……………………………..……….....
Battement de la nappe……………………………………………………………….....……
Evolution de l’état piézométrique ……………………………………………………....…..
Interprétation de la carte piézométrique (Mai 2013) …………………………………..…
Evolution spatio-temporelle de l’état piézométrique de l’aquifère ……………....………
Carte d’égale épaisseur de la couche mouillée……………………………..……….…..…..
Les caractéristiques hydrodynamiques ……………………………………….………....….
Interprétation des données……………………………………………………………..….
 Courbe de remontée…………………………………………………………….………....
Répartition de la transmissivité et de la perméabilité……………………………………….
Carte de transmissivité ……………………………………………………………………
Carte de perméabilité…………………………………………………...……………...….
. Conclusion ……………………………………………………………………...……...….

MODELISATION
Introduction………………………………………………………………………...........…
Formulations mathématiques…………………………………………………………….….
Loi de Darcy ……………………………………………………………………………...
Equation de conservation de la masse………………………………………..…………...
Equations d’état…………………………………………………………….………….….
Equation de diffusivité………………………………………………………..….……….
Résolution par les différences finies…………………………………………….………...
Conditions initiales et conditions aux limites………………………………..…….……….
Conditions initiales……………………………………………………………….……….
Conditions aux limites………………………………………………………...….……….
Les étapes de l’élaboration du modèle………………………………………………………
. Introduction des données………………………………………………………..………...
Calage du modèle……………………………………………………………..…………..
Exploitation du modèle…………………………………………………….……………...
Application du modèle en 3d au système aquifère……………………………….…………
Introduction………………………………………………………………………………..
Objectifs…………………………………………………………………….………..……
Discrétisation du domaine…………………………………………………..…………….
Données initiales………………………………………………………………….………
Calage du modèle en régime permanent……………………………………….……….....
Optimisation du champ de la perméabilité……………………………………………...
. Bilan hydrogéologique issu du calage………………………………………………......
. Exploitation du modèle et simulation……………………………………………………….
Conclusion………………………………………………………………………………..…

CHAPITRE V : HYDROCHIMIE
Introduction………………………………………………………………......…….……….
Matériels et moyens d’analyse………………………………………………………...........
Caractéristiques physico-chimiques………………………………………….……………..
Caractérisation des eaux de l’aquifère……………………………………………............
Répartition spatiale des paramètres physico-chimiques………………………….............
Paramètres physiques……………………………………………………………...........
Potentiel d’hydrogène (pH)………………………………………………….………..
La conductivité électrique……………………………………………………............
Paramètres chimiques………………………………………………………….……….
Le calcium……………………………………………………………….….………..
Le magnésium……………………………………………………....……….………..
 Le sodium…………………………………………………………………..…………
Les bicarbonates………………………………………………………….….………..
Les chlorures…………………………………………………………………...........
Sulfates……………………………………………………………………....……….
Les nitrates…………………………………………………………….………………
Evolution spatiale des nitrates………………………………………….….………..
Influence de la profondeur de l’eau sur la teneur en nitrates…………….…………
Nitrites………………………………………………………………………………...
Ammonium………………………………………………………………….………...
Phosphates………………………………………………………………..…………...
Evolution spatiale des fores azotées…………………………………………..………...
Faciès chimique des eaux………………………………………………………….………..
Diagramme de SCHOELLER-BERKALLOF…………………………………………….
Diagramme de PIPER……………………………………………………….…………….
Répartition spatiale des faciès chimiques……………………………………….………...
Analyse statistique des données physico- chimiques……………………………..………...
Corrélation entre les éléments…………………………………………………..………...
Analyse en composantes principales………………………………………….….............
Choix de la méthode………………………………………………………...…………..
Résultats et interprétations……………………………………………...……………….
Equilibres chimiques et indices de saturation……………………………………………….
Degré de saturation d’une solution vis-à-vis d’un minéral……………………..…………
Etat des minéraux dissouts dans l'eau…………………………………………..…………
Aptitude des eaux à l’irrigation……………………………………………….......………...
Introduction………………………………………………………………….…................
Risque de sodicité et de salinité (Méthode de Richard)………………………..………...
Classification des eaux en fonction du pourcentage du sodium absorbé par l’eau………..
Etude de la potabilité des eaux……………………………………………….…..………...
Introduction……………………………………………………………………................
La salinité…………………………………………………………………………………
La dureté totale……………………………………………………………………………
Potabilité des eaux……………………………………………………………….............
Conclusion………………………………………………………………………….............

CHAPITRE VI : VULNERABILITE A LA POLLUTION DE LA NAPPE

ALLUVIALE

Introduction ……………………………………………………………………….………...
Les méthodes d’estimation de la vulnérabilité…………………………………….………..
Vulnérabilité de la nappe alluviale par la méthode DRASTIC…………………..…………
Présentation de la méthode DRASTIC………………………………………….………...
Description des paramètres et réalisation des cartes thématiques………………………...
Profondeur de la nappe (D)……………………………………………….…..………..
La recharge efficace (R)………………………………………………...…….………..
 Nature de la zone saturée (A)…………………………………………….................…
La nature du sol (S)………………………………… ………………………...……….
La topographie (T)……………………………………………………..………………
Nature de la zone non saturée (I)…………………………………………….………...
Conductivité hydraulique (C)………………………………………………..…………
Carte de synthèse (vulnérabilité à la pollution)…………………………………………...
. Vulnérabilité de la nappe alluviale par la méthode de l’indice de contamination….............
Présentation de la méthode…………………………………………………………….….
Carte de la vulnérabilité à la pollution………………………………………...………..…
Conclusion………………………………………………………………..……………..…..

CHAPITRE VII : IMPACT DES EAUX USEES SUR LA QUALITE

PHYSICOCHIMIQUEDES EAUX SOUTERRAINES ET EAUX DE SURFACE

Introduction…………………………………………………………………………….……
Pollution d’oued Rhumel………………………………………………………………...….
Différentes sources de pollution………………………………………………...……….….
Pollution agricole……………………………………………………………………....….
Pollution d’origine domestique ………………………………………………..............…
Pollution d’origine industrielle……………………………………………………...…….
Caractérisation de la pollution……………………………………………………………....
Sites d’échantillonnage……………………………………………………………...…….
Résultats des analyses physicochimiques…………………………………...…………....
Evolution de la pollution………………………………………………………………....…
Interaction eaux de surface- eaux souterraine…………………………………………….…
Conclusion………………………………………………………………………..............…
CONCLUSION GENERALE……………………………….……………………................…
RECOMMANDATION……………………………………………………..…....………….…
REFERENCES BIBLIOGRAPHIQUE…………………………….…………….................…
ANNEXE…………………………………………………………………..……...……………
 Liste des figures
N° de figure Titre Page
Figure.1 Localisation du sous bassin du Tadjenanet – Chelghoum Laid …………….… ………....

Figure. Localisation du barrage Hammam Grouz à l’aval de la plaine (image satellitaire)………….

Figure.3 Localisation de la zone d’étude sur le plan régional ………………………………………..

Figure.4 Réseau hydrographique de la région d’étude ………………………………………………

Figure 5. Schéma structural de la méditerranée occidentale. (M. Durand Delga, )…………….…

Figure.6 Carte géologique de la région d’étude- d’après [Link] ( 8 ) Extrait de la carte de Sétif
- …………………………………………………………………………………………...

Figure.7 Carte géologique du sous bassin de Tadjenanet - Chelghoum Laid Selon J .[Link] 1980.....

Figure.8 Coupes géologiques à travers la région d’étude ………………………………………….

Figure.9 Extrait de la carte tectonique de [Link] ………………………………………….

Figure.10Description lithologique à partir des logs des forages implantés dans le terrain d’étude……

Figure.11 Carte bioclimatique de l’Est Algérien (Cote M., 8) ……………………………………

Figure.12 Répartition saisonnière des précipitations. Station de Hammam Grouz………………….

Figure.13 Variation des précipitations moyennes mensuelles ……………………………………….

Figure.14 Variations des précipitations moyennes [Link] Hammam Grouz (1988/89 -

3 )……………………………………………………………………………………………..

Figure.15 Variation des températures moyennes mensuelles (1988 - )………………………..

Figure.16 Variations des températures moyennes interannuelles…………………………………….

Figure.17 Diagramme ombro-thermique ……………………………………………………………..

Figure.18 Abaque de l’Indice d’aridité annuel de DE Martone ………………………………………

Figure.19 Bilan hydrique simplifié de la station de Hammam Grouz (1988/1989 - 3 )……

Figure.20 Bloc diagramme du bassin versant de Tadjenanet – Chelghoum laid………………………

Figure.21 Courbe hypsométrique du sous bassin de Tadjenanet – Chelghoum Laid ………………..

Figure.22 Profil en long de l’amont de Oued Rhumel et ses affluents ……………………………..

Figure.23 Carte de position des sondages électriques dans la plaine de Tadjenanet-Chelghoum Laid..

Figure.24 Courbe du sondage électrique A5. ………………………………………………………...

Figure. Courbe du sondage électrique d'étalonnage B3……………………………………………

Figure.26 Carte de résistivité AB = m………………………………………………………….

Figure.27 Carte de résistivité AB = m …………………………………………………………..

Figure.28 Carte de résistivité AB = 3 m…………………………………………………………

Figure.29 Coupes géoélectriques A, C, D et E ………………………………………………………

Figure.30 Coupes géoelectriques B’ ,C’,D’ et E’ ……………………………………………………

Figure.31 Carte du substratum ……………………………………………………………………….

Figure.32 Carte de la résistance transversale …………………………………………………………

Figure.33 Coupe hydrogéologique …………………………………………………………………..

Figure.34 Carte d’inventaire des points d’eau ……………………………………………………...

Figure.35 Carte piézométrique. Mai 7 …………………………………………………………..

Figure.36 Carte piézométrique. Septembre 7 …………………………………………………...

Figure.37 Carte de battement de la nappe …………………………………………………………....

Figure 38. Carte piézométrique. Mai 3 ………………………………………………….......……

Figure 39. Evolution spatio- temporelle de la piézométrie selon le sens d’écoulement………………

Figure.40 Carte d’égale épaisseur de la couche mouillée ……………………………………………

Figure.41 Essai de pompage – courbes de remontée …………………………………………….......

Figure.42 Carte de transmissivité en 10- m s ……………………………………………………..

Figure.43 Carte de perméabilité en 10- m s ………………………………………………………..

Figure.44 Résolution par modèle des problèmes hydrogéologiques…………………………………..

Figure 45. Bilan des flux sur une maille (i)……………………………………………………………

Figure.46 Discrétisation de la zone d’étude ………………………………………………………...

Figure.47 Reconstitution de l’état piézométrique en régime permanent …………………………….

Figure.48 Carte de la perméabilité de la nappe issue du modèle …………………………………...

Figure.49 Bilan hydrique de la nappe alluviale de Tadjenanet – Chelghoum Laid issu du modèle en
régime permanant ……………………………………………………………………………………...

Figure 50. Reconstitution de l’état piézométrique en régime transitoire. Mai 3……………………………

Figure.51 Simulation d’une exploitation à long terme …………………………………………….

Figure.52 Evolution temporelle de la piézométrie (scénario ) …………………………………….

Figure.53 Bilan hydrique de la nappe alluviale de Tadjenanet – Chelghoum Laid issu du modèle en
régime transitoire …………………………………………………………………………………….

Figure.54 Carte piézométrique résultant du deuxième scénario……………………………………..

Figure 55. Carte piézométrique résultant du troisième scénario …………………………………..…

Figure.56 Carte d’inventaire des points d’eau ayant fait l’objet d’une analyse physico- chimique…...

Figure.57 Coefficient de variation des variables physico-chimiques mesurées des eaux souterraines..

Figure.58 Carte du pH ……………………………………………………………………………….

Figure Carte de la conductivité électrique ………………………………………………………..

Figure.60 Carte de calcium …………………………………………………………………………

Figure.61 Carte du magnésium ……………………………………………………………………..

Figure Carte de sodium …………………………………………………………………………

Figure.63 Carte des bicarbonates …………………………………………………………………

Figure Carte des chlorures ………………………………………………………………………

Figure Carte des Sulfates ………………………………………………………………………..

Figure.66 Cycle de l’azote dans le sol (Hill, ) ………………………………………………...

Figure.67 Carte des nitrates ………………………………………………………………………..

Figure.68 Evolution spatiale des nitrates selon l’écoulement …………………………………….

Figure.69 Relation nitrate – profondeur de l’eau …………………………………………………...

Figure.70 Mécanisme de la pollution de l’eau par les nitrates. Julien, 3 ……………………….

Figure.71 Evolution spatiale des formes azotées et les phosphates selon l’axe d’écoulement……..

Figure.72 Diagramme de Schoeller-Berkallof. Mai ………………………………………….

Figure.73 Diagramme de Piper. Mai 7…………………………………………………………

Figure.74 Répartition des faciès chimiques selon les rapports caractéristiques ………………….

Figure.75 Carte des faciès chimiques .Mai 2007 ………………………………………………...

Figure.76 Répartition spatiale des faciès chimique de la nappe alluviale .Mai 7……………..

Figure.77.A. Corrélation entre la conductivité électrique et quelques éléments chimiques majeurs...

Figure.77. B. Corrélation significative entre quelque élément chimique…………………………..

Figure.78 Cercle de corrélation et répartition des variables………………………………………..

Figure.79 Variation de l’indice de saturation des minéraux carbonatés et évaporitiques dans les eaux
souterraines de la nappe alluviale…………………………………………………………………….

Figure.80 Classification des eaux pour l'irrigation selon la méthode de SAR …………………….

Figure.81 Carte d’aptitude des eaux à l’irrigation selon la méthode de Richard .Mai 7………..

Figure.82 Classification des eaux pour l'irrigation selon la méthode de Wilcox …………………..

Figure.83 Carte d’aptitude des eaux à l’irrigation selon la méthode de Wilcox .Mai 7…………

Figure.84 Diagramme de potabilité des eaux de la nappe alluviale Tadjenanet- Chelghoum Laid –
Normes Algérienne selon Boudoukha et al. 1999 ……………………………………………………

Figure.85 Carte indicielle : profondeur de l’eau .Mai 7……………………………………….

Figure.86 Carte indicielle : Recharge efficace …………………………………………………...

Figure.87 Carte indicielle : Matériau aquifère …………………………………………………..

Figure.88 Carte indicielle : Type du sol …………………………………………………………..

Figure. 89 Carte indicielle : Pente du terrain ……………………………………………………..

Figure.90 Carte indicielle : Impact de la zone non saturée ……………………………………….

Figure.91 Carte indicielle : Conductivité hydraulique …………………………………………..

Figure.92 Carte de vulnérabilité à la pollution des eaux souterraines de la nappe alluviale de
Tadjenanet – Chelghoum Laid selon la méthode DRASTIC ………………………………………..
Figure.93 Carte de l’indice de contamination …………………………………………………….
Figure.94 Carte de vulnérabilité à la pollution selon la méthode de l’indice de contamination…...
Figure.95 Carte d’inventaire des points de prélèvement des eaux de surface. Mai2007 et les
différentes sources de pollution dans la région d’étude ……………………………………………..
Figure.96 Photos représentant les stations de prélèvements des eaux de surface …………………
Figure.97 Evolution spatiale des éléments de pollution des eaux de surface sur le long de Oued
Rhumel ………………………………………………………………………………………………
Figure.98 Interaction eau de surface –eau souterraine à l’aval d’Oued Rhumel ……..……………..
 Liste des tableaux
N° du tableau Titre Page
Tableau.1 Coordonnées de la station climatologique de Hammam Grouz…………………………

Tableau.2 Les précipitations moyennes mensuelles, saisonnières et annuelles. Station de Hammam

Grouz 1988/89 à 2003/04 ………………………………………………………………………...

Tableau.3 Les précipitations et températures moyennes annuelles à la station de Hammam Grouz

1988/89 à 2003/04 ……………………………………………………………………………………

Tableau.4 Coefficient pluviométrique………………………………………………………………

Tableau.5 Les températures moyennes mensuelles, saisonnières et annuelles à la station de

Hammam Grouz. 1988/89 à 2003/04 ………………………………………………………………...

Tableau.6 Indice d'Euverte à la station de Hammam Grouz. 1988/89 à 2003/04 ………………...

Tableau.7 Evapotranspiration potentielle mensuelle. Station de Hammam Grouz. 1988/89 à

3 ……………………………………………………………………………………………….

Tableau.8 Calcul de l'ETR par les différentes méthodes ………………………………………….

Tableau.9 Bilan hydrologique de la station de Hammam Grouz (1988-2004) selon la relation de

Thornthwaite ……………………………………………………………………………………….

Tableau.10 Répartition altimétrique du bassin versant de Tadjenanet – Chelghoum Laid…………

Tableau.11 Calage du SE B3 sur les données du forage de Chabet Rmali ………………………..

Tableau.12 Paramètres hydrodynamiques de l’aquifère superficiel …………………………….......

Tableau.13 Bilan hydrique calculé en régime permanant…………………………………………..

Tableau.14 Bilan hydrique calculé en régime transitoire ………………………………………….

Tableau.15 Caractéristiques statistiques des paramètres physico-chimiques des eaux la nappe

alluviale de Tadjnanet- Chelghoum Laid. Mai 7………………………………………………….

Tableau.16 Estimation de la qualité de l’eau selon le pH …………………………………...……..

Tableau.17 Répartition des teneurs en nitrates dans la région d’étude. Mai 7………………...

Tableau.18 Coefficient de corrélation linéaire selon Bravais- Pearson…………………………

Tableau.19 Matrice de corrélation des paramètres physico-chimiques des eaux souterraines de

l’aquifère Tadjenanet-Chelghoum Laid ……………………………………………………………

Tableau.20 Caractéristiques de l’ACP appliquée aux données physico-chimiques des eaux de la

nappe alluviale …………………………………………………………………………………….
Tableau.21 Corrélations variables – facteurs …………………………………………………...

Tableau.22 Indice de saturation (SI) des minéraux présents dans l'eau ………………………….

Tableau.23 Classification des eaux souterraines de la région de Tadjenanet-Chelghoum Laid selon

les valeurs de TDS (Davis et De Wiest, 7) ……………………………………………………..

Tableau.24 Classification de l’eau selon la dureté totale ………………………………………..

Tableau.25 Méthodes d’estimation de la vulnérabilité choisies ………………………………….

Tableau.26 Poids attribués aux paramètres DRASTIC (Aller et al, 1987) ……………………...

Tableau.27 Classes et notes des sept paramètres (Lallemand- Barèes .1994)................................

Tableau.28 Notation et pondération selon la profondeur ………………………………………..

Tableau.29 Notation et pondération du paramètre de la recharge ……………………………….

Tableau.30 Notation et pondération du matériau aquifère ……………………………………….

Tableau.31 Notation et pondération de la nature du sol …………………………………………

Tableau.32 Notation et pondération du facteur de la pente ……………………………………...

Tableau.33 Notation et pondération de la zone non saturée ……………………………………..

Tableau.34 Notation et pondération de la conductivité hydraulique …………………………….

Tableau.35 Intervalle des classes de l’indice de contamination ………………………………….

Tableau.36 Indice de contamination ……………………………………………………………..

Tableau.37 Résultats d’analyses chimiques des eaux de surface…………………………………

GENERALITES ET ETAT DES CONNAISSANCES
 GENERALITES

1. Introduction

Les besoins de l’homme en eau augmentent considérablement dans le temps et touchent

des domaines de plus en plus nombreux pour satisfaire ses exigences de confort et
d’agrément et pour répondre à une demande économique croissante.
Cependant, les ressources en eau sont inégalement réparties dans le monde. L’eau douce,
celle qui est essentielle à nos besoins, ne représente que 1 % du total des eaux présentes
sur Terre, le reste étant formé par les eaux salées des mers et des océans. L’eau douce est
donc un capital limité, renouvelable mais fragile, car menacé par une consommation
croissante et par de multiples pollutions.

L’eau utilisable par l’homme se présente sous 3 formes : l’eau de pluie, les eaux de
surface (les fleuves, les lacs, etc.) et les eaux souterraines. Dans les régions arides et semi
arides, il est nécessaire de chercher l’eau dans les profondeurs de la terre à l’aide de puits
ou de forages, parfois très profonds. Dans les pays à climat non désertique, où les
ressources en eau sont issues des eaux de pluie et de surface, elles sont captées à l’aide de
barrages et de retenues collinaires.
L’eau, cette ressource rare et vulnérable reste parmi les principales préoccupations de
notre pays, durant ces dernières décennies. La rareté de cette denrée nous a incité à
adopter les techniques d’investigation plus modernes pour sa protection, vu les
méconnaissances des potentialités aquifères, leurs exploitation et leur gestion
anarchiques.

2. Problématique
La région de ‘’Tadjenanet – Chelghoum Laid’’ est située en Algérie orientale dans les
hautes plaines sétifiennes à environ 60 km au Sud du chef-lieu de la wilaya Mila. Il
s’agit d’une plaine de superficie de 1000 km environ, de nature lithologique hétérogène,
caractérisée par un matériel détritique varié d’âge Mio-Plio- Quaternaire. Certaines
formations peuvent constituer un réservoir important d’eau souterraine (alluvions, sable,
calcaire lacustre, gravier).
Sa situation en zone semi-aride implique une évaporation prononcée à même de basculer
le bilan dans un déficit chronique, hormis les périodes caractérisées par des épisodes
pluvieux.

- -
La vocation agricole de la région depuis l’époque coloniale, lui donne une importance
considérable dans ce domaine. D’autre part, la dispersion de la majorité des populations
dans les zones rurales, impose une implantation tout aussi diffuse des puits d’exploitation
d’eau souterraine et leur éloignement des grandes villes, pousse la population à se
focaliser sur l’activité agricole.
La diversification des cultures et la croissance des vergers sur la plaine ainsi que
l’accroissement des populations ont eu pour conséquences une augmentation des besoins
en eau pour l’irrigation et pour les besoins domestiques. Par conséquent le nombre de
puits n’a cessé d’augmenter d’une année à l’autre avec une demande croissante à cause de
l’implantation de forages de profondeur moyenne ainsi que l’implantation anarchique des
puits à usage individuel. Ceci a conduit à une mauvaise gestion des potentialités
hydriques de la région. Cette situation exige une étude approfondie des possibilités
hydrauliques du système aquifère de la plaine de Tadjenanet – Chelghoum Laid, ainsi que
l’effet de l’agriculture, l’installation des différentes unités industrielles dans les deux
villes, Tadjenanet et Chelghoum Laid. Cette situation peut entrainer un risque de
pollution de nature diverse, agricole, industrielle et domestique des eaux souterraine ainsi
que des eaux de surface (barrage de Oued el Athmania situé à l’aval de la plaine). Cette
situation est aggravée par l’absence de stations d’épurations dans la région et que tous les
rejets des eaux usées sont déversés dans l’oued du Rhumel sont aucun traitement
préalable.

Objectif et méthodologie
Notre tâche s’est limitée à une exploration sur terrain du réservoir aquifère superficiel de
la région de Tadjenanet –Chelghoum Laid, afin d’acquérir les données nécessaires, à cet
effet. On essayera à partir d’un traitement statistique et cartographique, de donner une
vision meilleure sur l’évolution des caractéristiques hydrodynamiques et chimiques des
eaux souterraines de la région. On abordera également des techniques géophysiques,
hydrodynamiques et hydrochimiques, permettant l’appréciation de la vulnérabilité de ce
système aquifère.
L’agriculture dans notre région d’étude représente la première activité socio-économique
et la plus consommatrice des ressources hydriques. Il est donc nécessaire de la prendre en
compte pour assurer une gestion efficace de l’eau. En effet la combinaison de la
cartographie géologique, géophysique, hydrogéologique et hydrochimique des données
récoltées sur terrain et leurs interprétations, peut constituer un excellent outil pour

- -
appréhender les grandes potentialités aquifères, leur répartition spatiale ainsi que leur
degrés de vulnérabilité aux diffèrent types de pollutions naturelle ou anthropique.
L’effet des changements climatiques a un impact sur le renouvellement des réserves
souterraines d’où la nécessité d’élaborer un modèle mathématique pour une bonne gestion
des ressources dans la plaine en fonction des variations climatiques qui affectent le bilan
hydrologique et en fonction de l’augmentation des besoins avec l’accroissement de la
population et des activités économiques. Pour atteindre ces objectifs, nous avons adopté
la méthodologie suivante :
- Dans une première étape : nous avons procédé à un inventaire complet de tous les points
d’eau (forage et puits), captant l’aquifère alluvial du Mio-Plio-Quaternaire ;
- Dans une deuxième étape : afin de définir l’évolution spatio- temporelle de la
piézométrie de l’aquifère en fonction des conditions climatiques, nous avons réalisé deux
compagnes piézométriques, l’une à fin de la période des hautes eaux (Mai 7), et
l’autre à la fin de la période de basses eaux (Septembre 2007).
Par ailleurs, nous signalons qu’une série d’essais de pompage a été réalisée durant la
période de Mai 2007.
Dans une dernière étape qui correspond à la dernière phase de la récolte des données:
nous avons effectué des prélèvements d’échantillon d’eau pour effectuer des analyses
chimiques des eaux souterraines de cet aquifère.
La dernière étape de notre travail, correspond à l’interprétation des données.
Le présent manuscrit se compose de six chapitres :
Dans le premier chapitre, nous présenterons le cadre géologique dans lequel nous avons
abordé la lithologie et la structure du terrain ;
Dans le deuxième chapitre, nous étudierons les paramètres hydroclimatiques afin de
quantifier les différentes composantes du bilan ;
Dans le troisième chapitre, nous étudierons la géométrie de l’aquifère par l’intermédiaire
de l’analyse de la prospection électrique qui a été effectuée dans la région ;
L’étude hydrogéologique a été abordée dans le quatrième chapitre, afin d’évaluer les
potentialités de cet aquifère par une modélisation de son comportement
hydrodynamique;
Dans le cinquième chapitre, nous traiterons le chimisme des eaux souterraines de cet
aquifère et son évolution spatiale à l’aide de la cartographie hydrochimique;
Dans le sixième chapitre, nous présenterons le concept de la cartographie de la
vulnérabilité à la pollution.
- -
Le dernier chapitre sera consacré à l’étude de l’impact des eaux usée sur la qualité
physicochimique des eaux souterraines et eaux de surface.
Finalement, nous donnerons une conclusion générale de cette étude, dans laquelle, nous
mettrons en évidence les résultats auxquels nous avons abouti, d’autre part nous
proposerons sur la base de ces résultats, des solutions adéquates pour préserver cette
ressource, ainsi qu’un programme de travail complémentaire pour une éventuelle étude
supplémentaire.
Caractères physico-géographiques
4.1. Situation géographique
La région explorée fait partie de la zone occidentale du grand bassin versant de Kébir
Rhumel situé à l’Est Algérien (fig. , 2 et 3). Cette plaine est située à 60 km environ au
sud ouest du chef lieu de la wilaya de Mila. Au Nord, elle est limitée par la commune de
Benyahia Abderrahmane, au Sud par la commune de Oueled Khlouf , à l’Est par Oued
Athmania et à l’Ouest par Bir El Arch. Cette région, a une superficie d’environ km
et administrativement, elle est rattachée à la wilaya de Mila.
La région est traversée par un cours d’eau appelé Oued Rhumel qui traverse la plaine
d’Ouest en Est ; il prend sa source vers 1 160 m dans les marges méridionales du Tell, au
Nord-Ouest de Bellaa. Il draine, suivant une orientation Sud-Ouest–Nord-Est, la région
Tadjenanet-Chelghoum Laid, puis franchit le Rocher de Constantine en s’encaissant
profondément (200 m de dénivellation) dans des gorges calcaires de [Link]
région est caractérisée par un climat semi-aride, avec un hiver pluvieux et froid et un été
sec et chaud.

Figure . Localisation du sous bassin du Tadjenanet – Chelghoum Laid

- -
Figure . Localisation du barrage Hammam Grouz à l’aval de la plaine
(Image satellitaire)

Figure . Localisation de la zone d’étude sur le plan régional

- -
 Orographie
La plaine est caractérisée par la présence de quelques massifs carbonatés tel que Djebel
Tnoutit (1170 m), Djebel Chebka (1044 m) et Maazoula El-Guablia (984 m). On
remarque également la présence de quelques collines réparties à travers la plaine.
4.3. Réseau hydrographique
De point de vue hydrographique la zone étudiée correspond à la haute vallée de l’Oued
Rhumel, bassin presque circulaire, fermé en aval par l’extrémité orientale de [Link].
Elle est drainée par oued Rhumel qui ne coule qu’à la faveur des précipitations à l’amont
du barrage. Il a donc un écoulement temporaire de direction Ouest – Est, sa pente est
régulière avec 83 m d’altitude à Tadjenanet, 720 m avant Hammam Grouz, et 760 m à
mi-distance aux environs de Chelghoum Laid. (Fig. ). A l’amont du barrage, les
ressources souterraines du bassin assurent à l’oued Rhumel, à l’aval du barrage, un
soutien en période d’étiage de 90 à 100 l/s en moyenne. (ANRH de Constantine).

Chelghoum Laid

Tadjenanet

Figure . Réseau hydrographique de la région d’étude

4.4. Géomorphologie
Du point de vue géomorphologique, on distingue les trois formes.
Relief
La région étudiée présente des reliefs importants avec une altitude comprise entre 1044
m et 7 m, caractérisant la partie Sud et Sud Est du terrain. Il s’agit des massifs

- -
carbonatés de Djebel Tnoutit, Djebel ed Dess , Djebel Grouz ,Djebel Tikouya et Djebel
Chebka.
Plaine
Sur le plan géologique, la plaine est caractérisée par une large dominance des formations
quaternaires (alluvions et terrasses, glacis souvent encroûtés) et Mio-Pliocènes (argiles et
marnes, avec niveaux de calcaires lacustres, grès et conglomérats). La majorité du terrain
étudié est caractérisé par une surface relativement plane dont la pente n’excède pas %,
avec des altitudes qui varient entre 7 m à l’Est et 930 m à l’Ouest. Cette disposition
topographique favorise le ruissellement des eaux de surface vers les dépressions de
l’Ouest vers l’Est.
4.5. Végétation
L’aspect général de la région est typique des hautes plaines sétifiennes à paysage
faiblement vallonné, de grandes étendues sans reliefs et de teinte généralement rougeâtre.
Région à climat beaucoup plus continental et semi-aride, est en grande partie cultivée par
des céréales. Seuls quelques lambeaux de forêts de chêne vert, à l’état de broussailles,
sont localisés sur les éminences calcaires. L’affleurement des formations carbonatées du
Crétacé sont en majeure partie dénudées. La région au printemps n’est qu’un immense
tapis vert, à l’automne après les labours il n’y a plus qu’une seule couleur, le brun
rougeâtre de l’argile interrompue par des taches claires ou grises de la croûte calcaire
récente.

- -
CHAPITRE I

CADRE GEOLOGIQUE
 CADRE GEOLOGIQUE

1. Introduction
L’étude géologique a pour but de préciser succinctement les terrains rencontrés et leur
contexte géologique. La feuille au . de Sétif fournit à l’exception de la chaîne
calcaire Kabyle et du socle Kabyle, et à partir de l’étude originale menée par Vila ( 8 ),
une série très complète des unités de l’orogène alpin d’Algérie orientale.
Du Nord au Sud et de haut en bas, cet édifice a pu être distingué (fig. 5):
- Le domaine des flyschs ;
- Le domaine tellien ;
- Le domaine néritique de la plate forme Sétifienne ;
- Les domaines parautochtones et autochtone hodnéen.

Zone d’étude

Figure 5. Schéma structural de la méditerranée occidentale.

(M. Durand Delga, 1969)

Au domaine des flyschs doivent être rattachées : la nappe numidienne, la nappe du flysch
mauritanien et la nappe du flysch Massylien.
Au domaine tellien appartiennent plusieurs unités : la nappe ultratéllienne, la nappe de
type Barbacha, la nappe de type Draa El-Arba, la nappe de type Beni-Abdellah, l’unité
supérieure à matériel Eocène , la nappe Bibanique , la nappe de Djemila et l’unité
inférieur à matériel Eocène, la nappe du Djebel Sattor, les unités à matériel Eocène et à
Huître et la nappe péni-tellienne.
Dans la plate forme allochtone, la nappe néritique Constantinoise est relayée vers l’Ouest

- -
par l’ensemble allochtone Sud sétifien. La région de Tadjenanet – Chelghoum Laid est
située à cheval entre le domaine tellien au Nord et hodnéen au Sud.
Du point de vue stratigraphique, les formations rencontrées sont en majeur partie
d’origine continentale et lacustre.
2. Description générale des diverses unités géologiques
L’examen des travaux de Vila 8 et Durozoy , ainsi que les logs stratigraphiques
des forages implantés dans le terrain d’étude, ont permis de visualiser les différentes
séries stratigraphiques. Les principaux affleurements rencontrés dans la région d’étude
sont de haut en bas:
- Les terrains récents du Quaternaire.
- le Miocène et le Pliocène,
- l’Eocène inférieur,
- le Crétacé supérieur et moyen
- le Trias
2.1. Tertiaire et Quaternaire
Le Miocène inférieur est formé de minces lambeaux d’argiles gréseuses se trouve tout à
fait au Sud de la zone étudiée. Les terrains du Miocène supérieur (Pontien), Pliocène et
Quaternaire, sont tous d’origine continentale. On n’y retrouve pas la continuité latérale
classique des terrains sédimentaires marins. Il s’agit essentiellement d’argiles rouges dans
lesquelles s’étendent des bancs de calcaires lacustres, grés et conglomérat, toujours
restreints en épaisseur et en étendue, avec changements latéraux de faciès rapides. Les
argiles, en certains endroits (Nord-Ouest, vers le Dj. Ed Dess, au Sud-est du [Link] et
autour des affleurements éocènes, contiennent de nombreux blocs ou cailloutis ou débris
de silex. Vers le Sud-est de la zone étudiée, contre le relief qui s’étend du [Link] au
[Link], nous avons constaté un redressement, jusqu’à la verticale des couches
tertiaires. En surface, il existe assez souvent une couche dure, calcaire englobant des
cailloutis, de quelques mètres au plus, appelée croûte calcaire.
2. 2. Eocène méridional
A la base se trouvent des calcaires massifs cristallins à nombreuses Nummulites (150m),
au sommet se trouvent des marnes jaunes.
2.3. Eocène centrale
Entre Tadjenanet et Chelghoum Laid, à Draa Mestaoua et Karouach à la base de
Thanétien, Yprésien et Lutétien inférieur, se trouve une assise de calcaire gris et de
calcaire marneux à silex (280 m), se débitant à l’affleurement, en blocs décimétriques. Au

- -
sommet (Lutétien supérieur) se trouvent des marnes jaunes (300 m), gypseuses parfois.
Durozoy, 1960.
2.4. Eocène inférieur
La lithologie de cette série est différente suivant que nous la considérons au centre de la
zone étudiée ou au Sud.
2.5. Crétacé supérieur et ultra-supérieur
Cette série est subdivisée par les géologues en de nombreux étages (Cénomanien
supérieur, Emschérien-Santonien-Auturien, Maestrichtien, Montien. Ces deux derniers
regroupés parfois sous le terme de Dano-Montien. Elle est essentiellement marneuse et
argileuse, et probablement très épaisse (centaine de mètres). Au Sud du dj Chebka, elle
est gypseuse. Le Maestrichtien est calcaire dans l’extrême Ouest de la zone. Ces niveaux
sont notés par Durozoy (1960), dans sa carte géologique de Tadjenanet, intercalés entre
les écailles des calcaires éocènes au voisinage de la RN 5 entre Tadjenanet et Chelghoum
Laid.
2.6. Crétacé inférieur et moyen
Le Crétacé constitue une unité épaisse de 500 à mètres (ou plus, si l’on considère les
calcaires jurassiques sous-jacents) de nature calcaire pour l’essentiel. Ces calcaires sont
massifs, cristallins ou dolomitiques, ils se débitent, à l’affleurement, en blocs dépassant
parfois plusieurs mètres cubes. Ils forment Dj. Grouz, qui correspond à un vaste
bombement de l’unité crétacée et le massif de Dj. Tnoutit.
2.7. Le Trias
Ce sont des argiles rouges gypseuses et salifères. Elles apparaissent au Nord vers le Dj Ed
Dess et au Sud et Sud-Est, en une série, d’affleurements, alignés Sud-Ouest à Nord-Est,
correspondant sûrement à un ‘’accident ‘’ tectonique profond. Les argiles triasiques
traversent les terrains plus jeunes du Crétacé et de l’Eocène en les découpant à
l’emporte-pièce et en les rebroussant suivant le style diapirique.
3. Lithostratigraphie de la région d’étude
L’étude lithostratigraphique de la région permet de diviser les formations de la région en
deux séries : (fig. 6, 7 et 8)
- Formations peu ou pas tectonisées ;
- Formations tectonisées.
3.1. Formations peu ou pas tectonisées :
De haut en bas on distingue :

- -
3.1.1. Quaternaire
Les formations Quaternaires occupent presque la totalité du terrain.
a. Alluvions récentes
Elles apparaissent sur tout le profil de Oued El-Mahri, ces dépôts sont formés d’une
banquette de limon gris et de graviers.
b. Terrasses, formation de pente et Quaternaire indéterminé
Ce sont des formations de pente correspondant à des terrasses dilacérées ou à des
alluvions anciennes sur lesquelles se sont développés des sols bruns peu ou pas calcaires.
Au Quaternaire indéterminé sont rapporté des placages dispersés de nature variée qui
occupent un grand espace sur les niveaux marneux.
c. Glacis polygéniques
Ils affleurent sur toute la périphérie des principaux reliefs (Djebel Tnoutit, Rokbet El-
Djemel). Ils correspondent à des surfaces très faiblement pentées près des plaines dont la
pente augmente progressivement vers la montagne, ils sont recouverts d’un matériel
clastique faiblement calibré et très fragmenté. Vers le bas, ils deviennent plus anguleux et
plus grossier.
d. Villafranchien : croûtes calcaires
Ce sont les formations les plus caractéristiques des hautes plaines sétifiennes, elles sont
constituées de bancs calcaires, des travertins, des croûtes calcaires, des carapaces
calcaires, de calcaire lacustre blanc et de limon avec une épaisseur de 50 m environ.
e. Sebkha
Ce sont des vases salées, souvent inondées pendant l’hiver et recouvertes par des croûtes
de sel en été, ils sont localisés au Sud et au Sud Est du terrain.
f. Sols salés
Il s’agit de limon salé à végétation typique de sol salsolacées correspondant à des sols
hydromorphes incultes, très mal drainés.
g. Formations dunaires
Ce sont des amas éoliens de sables calcareux et de poussière limoneuse mêlée à des
débris coquilliers fins (Gastéropodes) ou à des végétaux.

3.1.2. Mio – Pliocène continental

Il affleure sur une grande partie du terrain, ce sont des calcaires lacustres, des marnes
rouges, des limons et de sables rouges. L’ensemble de ces formations présente des faciès
très variables. En bordure des reliefs actuels, s’enrichissent en conglomérats rougeâtres

- -
mal cimentés, l’ensemble a une épaisseur de à 5 m environ.
3-2- Formations tectonisées
3-2-1- Unité supérieure à matériel Eocène
Il s’agit de lambeaux de séries avec les trois termes classiques de l’Eocène tellien, dont la
série commence parfois à la fin du Sénonien.
a- Yprésien et Lutétien :
Ce sont des calcaires marneux schistosés à patine blanche, à cassure noire, à silex noires
et à Globigérines.
b- Marnes fini-Crétacées et Paléocènes à boules jaunes.
C’est une série de marnes noires d’au moins 8 m d’épaisseur. Sa base est tronquée par
un contact anormal.
3-2-2- Nappe du Djebel Sattor et unités Eocènes à Nummulites et à Huîtres
a- Lutétien à Huîtres
Il s’agit d’une centaine de mètres de marnes jaunâtres, très riches en intercalations
Lumachellique à Ostrea multicostata.
b- Yprésien et Lutétien à débris
Ils forment une barre calcaire débutant par des niveaux phosphatés glauconieux et
comportant des silex à son tiers inférieur. D’abondant débris de Gastéropodes, de
Lamellibranches et d’entroques rendent l’aspect du faciès très granuleux.
c- Yprésien et Lutétien à Nummulites
Des calcaires identiques contiennent des Nummulites dans la partie Est de la carte au Sud
Est de Tadjenanet.
d- Sénonien supérieur Lumachellique
Il est marqué par le développement de gros bancs Lumachelliques à Ostréidés associés à
des marnes jaunes et à des calcaires sparitiques, parfois glauconieux. Il affleure au Sud du
terrain.
3.2.3. Nappe néritique Constantinoise
a. Aptien calcaire
Il débute par une barre à Rudistes de 15 m et se poursuit par 150 m de calcaire bien lité à
Salpingoporella.
b. Barrémien carbonaté
Le Barrémien franc débute par des calcaires plus ou moins dolomitiques à algues :
Salpingoporella Muelhbergü.
Au sommet, le Barrémien supérieur a livré : Palorbitolina lenticulaire avec quelques très

- -
rares niveaux marneux.
c. Néocomien calcaire
Il est formé par des dolomies cristallines, souvent bréchoides de couleur claire et mal
litée.
d. Jurassique carbonaté
La série attribuée au Jurassique supérieur, comprend des dolomies des micrites en
alternance avec des dolomies et calcaires dolomitiques.
3.2.4. Ensemble allochtone Sud-Sétifien
a. Miocène marin transgressif
Le Miocène débute par des conglomérats rougeâtres qui ravinent largement son
substratum. Il affleure au Sud du Rokbet El-Djemel où les niveaux marneux sont surtout
visibles.
b. Barrémien
Le Barrémien est formé par une série complexe composée de calcaire, de dolomies et de
grés avec de minces intercalations argileuses ou marneuses à Djebel Tnoutit.
3.2.5. Trias exotique ou extrusif
Il est formé par des argiles varicolores, gypses, cargneules, des glaçons calcaro-
dolomitiques. Il affleure au Nord et au Sud-est du terrain d’étude.

3. Analyse structurale
3.1. Ensemble Allochtone Sud-Sétifien
A quelques exceptions, l’ensemble allochtone Sud-sétifien est constitué par toute une
série d’écailles en relais constituant le trait dominant de la géologie structurale de la
région. Cet ensemble possède une série surtout carbonatée du Lias au Cénomanien.
La masse chevauchante de Djebel Tnoutit apparaît comme un élément de l’ensemble
allochtone Sud – sétifien. Vers l’Est son contact de base va rejoindre le front des masses
allochtones de la nappe néritique Constantinoise .Fig. .Vers l’Ouest cette surface est
visible jusqu’à la terminaison occidentale du Djebel Zdimm laquelle s’enfonce en tunnel
sous la nappe tellienne de la feuille de Ain Taghrout.

- -
Figure Carte géologique de la région d’étude- d’après [Link] (1980)
Extrait de la carte de Sétif -

- -
Figure 7. Carte géologique du sous bassin de Tadjenanet - Chelghoum Laid
Selon J .[Link] 1980

Figure 8. Coupes géologiques à travers la région d’étude

- -
 Figure 9. Extrait de la carte tectonique de Keiken, 1962

Un aperçu de la région nous permet de distinguer deux types de tectoniques :

-Une tectonique anté-miocène de style cassant qui favorise la karstification ;
-Une tectonique post-miocène de style souple.
La première se manifeste dans la partie Est et la deuxième dans la partie Ouest en
accentuant la tectonique précédente.
Suivant les caractéristiques structurales et stratigraphiques des terrains, on peut les
subdiviser en trois unités principales :
. unités du Trias diapir
. unités du crétacé –éocène.
3. unités du Moi-pliocène transgressives
3-2-Unités du Trias diapirique
De faciès gypso-salin avec des lambeaux de calcaire, leur remontée a produit des
charriages, des déversements et des chevauchements
3-3-Unités du Crétacé-Eocène
Caractérisant les aires anticlinales de directions NE-SW avec un faciès néritique du
Crétacé moyen et inférieur. Ces massifs sont découpés par des failles longitudinales. Les
zones synclinales (zone d’effondrement) de facies marneux du Crétacé supérieur et

- -
Nummulitique, sont en plus déversées et faillés de direction NW- SE au voisinage des
massifs et EW aux axes des synclinaux.
Ces mouvements ne sont pas produits à la fin du Lutétien et sont accentués à la fin de
l’Oligocène
3.4. Unités Mio-Pliocène transgressives
Ces unités caractérisent la partie Ouest de la région représentée par les bombements
anticlinaux qui se sont produits à la fin du Miocène et accentués après le Villafranchien.
4. Description géologiques à partir des logs de forage et des données géophysiques
Plusieurs forages ont été réalisés dans la région et seuls quelques-uns ont fait l’objet d’un
suivi relativement correct. Nous allons essayer de faire leur description tout en faisant le
lien avec la géologie et la géophysique.
Un forage a été réalisé à 4 Km au Nord de Tadjenanet, de coordonnées Lambert
suivantes : x=795.70 km et y=322.80 km vers le Chabet Rmali (Fig. 10). Sa profondeur
est de 278 m. Durozoy (1960) donne les indications suivantes, mis à part argiles et
marnes qui constituent l’essentiel :
- de 26.50 à 32.50m, calcaire lacustres
- de 72 à 79m, calcaire lacustres
- de 216 à 241m, galets
- de à 78m, calcaire de l’Eocène inférieur.
La géophysique parait affirmer le non enracinement des calcaires éocènes centraux, tout
au moins vers le Sud. Une solution serait la suppression par érosion de toute l’épaisseur
des calcaires éocènes sur le pourtour des amas ‘’ centraux’’, puis dépôts des terrains
tertiaires et quaternaires (fig. 10). Une seconde solution pourrait être apportée par les
explications tectoniques récentes concernant le domaine tellien de l’Algérie orientale.
Depuis les découvertes de la S.N REPAL (1952), de vastes mouvements de
chevauchement du Crétacé, et de charriage (nappes et écailles) touchant l’Eocène, depuis
le Nord, ont été mis en évidence dans la région de Bordj Bou Arreridj et Sétif et, dans
celle de Constantine. Ce sont les forages profonds, l’étude des disharmonies de structure
et les groupements de microfaunes fossiles qui ont permis de dresser des cartes
tectoniques nouvelles (KIEKEN .1962). Les études géologiques indiquent :
- un Crétacé inférieur et moyen, calcaire, épais (500,1000 m ou plus), découpé
en compartiments polygonaux par de très nombreuses failles (fig. 9).
- Un crétacé supérieur marno-argileux, qui peut être très épais (plus de 1000m
pour le Sénonien vers Bordj Bou Arreridj).

- -
 - Un Eocène inférieur, moins épais, qui présente le faciès à Nummulites et un
faciès sans Nummulites, de tectonique extrêmement souple avec plis déversés
vers le Sud pour le deuxième.
- Le Mio-plio-quaternaire, constitue le toit de la formation éocène (fig. 10).
La tectonique du Crétacé est nettement différente de celle de l’Eocène ‘’central’’. Elle
peut indiquer des mouvements indépendants du second par rapport au premier, facilité par
la plasticité des argiles et des marnes du Crétacé supérieur. Un argument dans ce sens
nous est apporté par Durozoy ( ), qui constate l’identité de faciès entre la zone de
Sétif (Eocène sans Nummulites, nappé) et celle ‘’centrale’’ de notre zone. On peut
comprendre ainsi :
- les plissements en accordéon de l’Eocène central ;
- le voisinage, à moins de 10 Km, des deux faciès ;
- Le passage sans transition de l’Eocène ‘’central’’ au vaste domaine Sud argilo
marneux (fig. 10).
Nous pouvons aussi supposer en allant dans le même sens, de grand dérangement dans
l’ordre géologique des marnes et argiles sous-jacentes à l’Eocène central, ce qui peut
poser des questions sur la nature géologique véritable du substratum profond et sur
l’éventuel rôle des argiles triasiques.

Forage 1 Forage 2

- -
Forage 3 Forage 4

Forage 5 Forage 6

- -
 Forge7

Figure 10. Description lithologique à partir des logs des forages implantés dans le
terrain d’étude

5. Conclusion

La région de Tadjenanet – Chelghoum Laid est située à cheval entre le domaine tellien au
Nord et hodnéen au Sud. Du point de vue stratigraphique, les formations rencontrées sont
en majeur partie d’origine continentale et lacustre. Il s’agit essentiellement d’argiles
rouges dans lesquelles s’étendent des bancs de calcaires lacustres, grés et de
conglomérat, toujours restreints en épaisseur et en étendue, avec changements de faciès
rapides. Les argiles, en certains endroits (Nord-Ouest, vers le Dj. Ed Dess, au Sud-est du
[Link] et autour des affleurements éocènes), contiennent de nombreux blocs ou
cailloutis ou débris de silex.
Du point de vue tectonique la région d’étude est caractérisée par deux styles de
tectoniques :
-Une tectonique anté-Miocène de style cassant qui favorise la karstification ;
-Une tectonique post-Miocène de style souple.
La première se manifeste dans la partie Est et la deuxième dans la partie Ouest.
L’ensemble allochtone Sud-sétifien est constitué par toute une série d’écailles en relais
constituant le trait dominant de la géologie structurale de la région. Cet ensemble possède
une série surtout carbonatée du Lias au Cénomanien. La masse chevauchante de Djebel

- -
Tnoutit apparaît comme un élément de l’ensemble allochtone Sud-sétifien. Vers l’Est, son
contact de base va rejoindre le front des masses allochtones de la nappe néritique
Constantinoise, caractérise la masse chevauchante de Djebel Grouz.
Ces masses chevauchantes de nature carbonatées ont une grande capacité hydrique
contribuant à l’alimentation de la nappe alluviale. Les formations du Mio-Pliocène
couvrant d’importantes superficies et ils ont de bonnes potentialités de point de vue
hydrogéologique, dont leur l’épaisseur varie de à 5 m .Ils sont constitué par des
dépôts fluvio- lacustres, de limon, de sable, marne et argiles.

- -
CHAPITRE II

CADRE HYDROCLIMATIQUE
 CADRE HYDROCLIMATIQUE

1. Introduction

La connaissance des paramètres climatiques d'une région est indispensable pour

l'évaluation des différents termes du bilan hydrologique. L'alimentation de l'aquifère se
fait en grande partie par les précipitations météoriques. La région de Tadjenanet –
Chelghoum Laid appartient au domaine des hautes plaines sétifiennes qui se caractérisent
par des pluies moyennes annuelles de l'ordre de 400 mm. Du point de vue hydrologique la
région est située dans la haute vallée du bassin versant de Kébir Rhumel avec une
superficie de l'ordre de 1130 km . Dans cette approche hydroclimatique, nous avons
utilisé les données de la station climatique de Hammam Grouz .On note que cette station
se trouve juste à la limite Est de notre terrain, précisément à l'exutoire du bassin versant,
dont les coordonnées sont résumés dans le tableau 1.

Tableau 1. Coordonnées de la station climatologique de Hammam Grouz

Station Oued Coordonnées
Hammam Longitude Latitude Z (m)
Rhumel
Grouz 6° ’ E 3 ° ’N 73

Caractéristiques climatiques de la région d’étude

L’Algérie, du fait de sa situation géographique se trouve partagée en trois zones
climatiques distinctes (fig.11); la première à climat méditerranéen au Nord, la seconde à
climat semi-aride vers l’intérieur du pays (hauts plateaux) et enfin un climat aride qui
caractérise le grand Sahara.
La région d’étude fait partie de la seconde zone, marquée par les risques et les aléas. Dans
ce qui suit nous présentons un aperçu hydroclimatique, indispensable à la compréhension
des processus d’alimentation et de restriction des eaux souterraines, ce qui nous permettra
d’avancer dans la connaissance du comportement hydrologique des cours d’eau, et dans
la variation des réserves des eaux souterraines. Ainsi, la variation des paramètres
climatiques donne des indications sur l’évolution des concentrations des éléments
chimiques, en fonction des saisons, tant au niveau des eaux souterraines qu’au niveau des
eaux de surface.

- -
 Figure Carte bioclimatique de l’Est Algérien (Cote, 1998)

3. Etude des paramètres climatiques

Les différentes données climatiques sont celle recueillies au niveau de la station
pluviométrique de Hammam Grouz (Tableau 2). Une période de seize ans (1988/89 à
2003/04) permet de caractériser des conditions climatiques locales de la région d’étude
comme on le verra plus loin.

3.1. Précipitations
La pluie est un facteur climatique très important conditionnant l’écoulement saisonnier et
par conséquent le régime des cours d’eau ainsi que celui des nappes. La pluviosité
augmente avec l’altitude et diminue au fur et à mesure que l’on s’éloigne de la mer.

Tableau 2. Les précipitations moyennes mensuelles, saisonnières et annuelles

A la station de Hammam Grouz de 1988/89 à 2003/04

Mois Sep Oct Nov Dec Jan Fev Mar Avr Mai Juin Juil Août Année
P(mm) 5.203 702.. 552.1 02205 33233 7121. 57234 53263 04204 61255 52.3 67211 527273
Automne Hiver Printemps Eté
Saison 97.23 mm 129.70 mm 112.45 mm 32.88 mm

- -
L'étude des précipitations montre que les saisons les plus pluvieuses sont l'hiver avec
34.84% des pluies moyennes annuelles et le printemps avec 30.21%. L'automne
représente 26.12% et l'été avec seulement 8.83% des précipitations moyennes annuelles.
Fig.12.
La représentation moyenne mensuelle des précipitations sur une période de seize ans,
nous a permis de constater que le mois de Janvier est le plus arrosé avec une moyenne
mensuelle égale à 55.59 mm alors que le mois de juillet est le plus sec avec une moyenne
mensuelle égale à 3.89 mm (fig.13). La moyenne interannuelle des précipitations pour la
période de 1988 à 2004 est de 372.25 mm (Tab. 3).

40
35
Précipitations (%)

30
25
20
15
10
5
0
Automne Hiver Printemps Eté
Saisons

Figure 12. Répartition saisonnière des précipitations. Station de Hammam Grouz

Figure 13. Variation des précipitations moyennes mensuelles

- -
 Tableau 3. Les précipitations et températures moyennes annuelles
à la station de Hammam Grouz 1988/89 à 2003/04.

Année
P (mm) . 73.3 7.5 5 3.5 7. 5. 5.7 5
T (°C) 5.5 . 8 . 8 3.7 5.5 5.8 5. 8 5.
Année
P (mm) .8 7 . 5 5.5 .7 78. 8.8 5 .
T (°C) 5. 5. .8 . .3 . . 5.5

On constate que l'année la plus humide est 2002-2003, avec une moyenne annuelle égale
à 668.80 mm alors que l'année la plus sèche est 1992-1993 avec une moyenne annuelle
égale à 147.6 mm. Figure 14.

Figure 14. Variations des précipitations moyennes interannuelles.

Station Hammam Grouz (1988/89 -

3.2. Coefficient pluviométrique

Le coefficient pluviométrique est un paramètre climatique qui permet de déterminer les

années déficitaires et les années excédentaires. Il est définit comme étant le rapport de la
pluviométrie de l'année considérée à pluviométrie moyenne pour une station donnée.
Cp = P/Pa
Où
Cp : Coefficient pluviométrique ;
P : pluviométrie de l'année considérée (mm) ;
Pa : pluviométrie moyenne pour une période d'année (mm)
- -
Si
Cp > 1 : l'année correspondante est excédentaire.
Cp < 1 : l'année correspondante est déficitaire.
D'après les valeurs des coefficients pluviométriques représentées dans le tableau n °4, on
constate que la station de Hammam Grouz présente 08 années excédentaires et 08 années
déficitaires sur toute la période d’observation ( 88 8 - 3 ).

Tableau 4. Coefficient pluviométrique

Années
Cp 2.11 2.1. 7.70 7..1 2.12 2.66 7.71 7..1
Années
Cp 2.17 7.06 2.10 7.70 2.60 2.11 7.12 7..1

3.3. Les températures

La température est un élément très important du climat, elle joue un rôle déterminant dans
l'estimation du déficit d'écoulement qui intervient dans l'estimation du bilan
hydrologique. Elle est liée à la durée d’insolation, à l'altitude et aux conditions
climatiques locales de la région d'étude (Tab.5).
Tableau 5.
Les températures moyennes mensuelles, saisonnières et annuelles
A la station de Hammam Grouz. 1988/89 à 2003/04
Mois Sep Oct Nov De Ja Fe Ma Avr Mai Jui Juil Ao Ann
c n v r ût ée
T 02.6 71.1 70.. 1.1 6.1 1.6 72.7 70.7 71.1 00.1 01.. 06.7 71.1
(C°) . 2 0 0 1 1 0 0 6 2 1 7 .
Sais Automne Hiver Printemps Eté
on 76.17 1..1 7..00 01.00

L'étude de la variation moyenne mensuelle de la température montre que les mois les plus
froids sont ceux de Décembre, Janvier et Février avec des températures moyennes
mensuelles qui varient entre 6 et 8 °C. Par contre les mois les plus chauds sont ceux de
Juin, Juillet et Août avec des températures moyennes mensuelles qui varient entre 22 et
26 °C (fig.15).

- -
 Figure 15. Variation des températures moyennes mensuelles (1988 - )

La représentation inter–annuelle des températures sur une période de seize ans (1988/89 à
2003/04 ), montre que l'année la plus froide est celle de 1991/92 avec une valeur
moyenne annuelle de 13.76 °C, alors que l'année la plus chaude est celle de 1998/99 avec
une valeur moyenne annuelle de 16.84°C( fig.16).

Figure 16. Variations des températures moyennes interannuelles

3.4. Diagramme pluviothermique

Selon Bagnouls et Gaussen une période sèche est une période pendant laquelle les
précipitations totales du mois considéré sont inférieures ou égales au double de la
température. La confrontation des données pluviométriques et celle des températures sur

- -
le diagramme pluviothermique (fig.17) révèle que la période sèche s'étale sur quatre mois
environ, du mois de juin au mois de septembre alors que la période humide s’étale sur le
reste de l’année.

Figure 17. Diagramme ombro-thermique

3.5. Indice d'aridité

Cet indice permet d'identifier le type du climat caractérisant la région étudiée. La formule
de De Martonne permet d'apprécier ce paramètre.
I = p/( t +10)
Où

P : précipitation moyenne mensuelle (mm) ;

T : température moyenne mensuelle (°C).
Pour :
20 < I <30 : le climat est tempéré;
10 < I < 20 : le climat est semi-aride,
7.5 < I <10 : le climat est steppique;
5 < I < 7.5 : le climat est désertique ;
I<5 : le climat est hyper aride.

- -
Pour notre région, comme t = 15.43° C et p = 372.25 mm permettent de donner un
indice d'aridité de l'ordre de14.64, ce qui implique un climat semi-aride (fig.18).

Indice d'aridité
I = P/ T+10
Précipitation moyenne annuelle (mm)

Abaque de l'indice d'aridité

50
2600 Ecoulement abondant

2400

2200 40
Ecoulement permanant
2000
1800
30Zone tempérée
1600
Drainage exterieur
1400

1200 20Régime semi-aride

Ecoulement temporaire
1000

800
10Régime désertique
600 Ecoulement temporaire
Station de Hammam Grouz 5
400 Drainage interieur

200

0
10 20 30
Température moyenne annuelle (°C)

Figure Abaque de l’Indice d’aridité annuel de DE Martone

3.6. Indice d'Euverte

Le rapport des précipitations et les températures moyennes établis mensuellement selon
Euverte, permet de connaître l'évolution de l'humidité du sol (Tab.6). Il permet également
de définir plusieurs types de régimes:
- Un régime très sec où p/t <1 ;
- Un régime sec où 1<p/t<2;
- Un régime subhumide où 2<p/t<3;
Un régime humide où p/t >.
Où
P : précipitation moyenne mensuelle (mm);
T : température moyenne mensuelle (°C)

Tableau 6. Indice d'Euverte à la station de Hammam Grouz. 1988/89 à 2003/04

Mois Sep Oct Nov Dec Jan Fev Mar Avr Mai Jui Juil Août
P/T 62.1 6207 7225 1246 .233 5202 5275 5275 7254 4225 4261 420.

- -
L’application de cette formule permet de constater que la région d'étude est caractérisée
par un régime de type très sec pendant les mois d’été (Juin, Juillet et Août) et un régime
sec pendant les mois de Septembre et Octobre. Le régime Humide caractérise la saison
hivernale, le reste de l'année est de type sub-humide

4. Etablissement du bilan hydrologique

Le but principal du bilan hydrologique est de comptabiliser les apports et les pertes en eau
d'un bassin versant, autrement dit de quantifier les transferts d'eau issus des précipitations.
L'équation générale du bilan au niveau d'un bassin versant s'écrit de la manière
suivante :
P (mm) = Ex (mm) + ETR (mm) avec Ex = I + R ±V
Où
P : précipitation moyenne annuelle (mm);
Ex : excèdent moyen annuel (mm) ;
ETR : évapotranspiration réelle moyenne annuelle (mm) ;
I : infiltration (mm) ;
R : ruissellement (mm).
V : variation de la réserve, généralement très faible.

4.1. Analyse des différents éléments du bilan

4.1.1. Précipitations
Pour établir le bilan hydrologique, nous avons utilisés données des précipitations
mesurées dans la station de Hammam Grouz sur une période de seize ans (1988/89 à
3 ).
4.1.2. Evapotranspiration
Ce terme regroupe à la fois l'évaporation à partir des surfaces libres et la transpiration des
végétaux. C'est l'élément le plus important du bilan hydrologique. Ce processus
d'évapotranspiration à partir des sols dépend essentiellement de l'humidité, de la
température, de l'insolation, de la couverture végétale.
[Link]. Evapotranspiration potentielle (ETP)
Notre région appartient au domaine semi aride, on a jugé bon d'utiliser la méthode de
Thornthwaite, pour faire sortir le maximum de paramètres. Cette méthode exprime le
pouvoir évaporant au maximum de l'atmosphère au sol en fonction de la latitude, plus
exactement de la durée théorique de l'insolation dans la région. La formule de cette
méthode est la suivante
ETP (mm) = 16(10t / I)a
- -
 Avec:
 ETP : évapotranspiration potentielle mensuelle (mm)
 I : indice thermique annuel ;
.5
 I= i=1∑
i=
i avec i = (t / 5) avec i: indice thermique mensuel.
 t : la température moyenne mensuelle (°C)
 a = (1.6 I /100)+ 0.5
Ces valeurs de l'ETP doivent être corrigées à l'aide d'un coefficient k qui dépend de la
latitude de la station qui prend en compte la durée de l'insolation théorique tirée d'un
abaque.
ETP corrigée = ETP.K (Tab. 7).
Tableau 7. Evapotranspiration potentielle mensuelle à la station de Hammam
Grouz. 1988/89 à 2003/04

Mois Sep Oct Nov Déc Jan Fév Mar Avr Mai Jui Juil Aoû

t (°c)
i
ETP non corrigée(mm)

ETP corrigée (mm)

[Link]. Les réserves facilement utilisables (R.F.U)

Ce terme peut être utilisable à partir de la formule de Ricard.J 1979.
R.F.U =1/3 [Link].P avec, Da = 1.4, He = 30%, P = 40 cm (Djouahra 1983).
Où : Da = densité apparente du sol ;
He : humidité équivalente du sol (capacité de rétention);
P : profondeur de la terre parcourue par les racines (m).
D’où la R.F.U = 5 mm

- -
[Link]. Evapotranspiration réelle
L'évapotranspiration réelle est la quantité d'eau évaporée et transpirée réellement à partir
du sol, des végétaux et des plans d'eau dans les conditions naturelles du bassin versant.
Plusieurs formules et méthodes permettent de déterminer ce permettre.
[Link].1. Méthode de Thornthwaite
On peut évaluer l'évapotranspiration réelle par l'établissement du bilan hydrologique
selon deux cas.
a / Si P> ETP, alors ETR = ETP
Dans ce cas l'excèdent provenant des précipitations est emmagasiné dans le sol pour
alimenter la RFU jusqu'à sa saturation. Son surplus sera attribué à l'excèdent (water
surplus).
b / si P< ETP, alors, ETR = P + R.F.U (toute ou une partie de la R.F.U).
Si R.F.U disponible est insuffisant pour satisfaire l’ETR, alors l'ETR = P, et elle sera
inférieur à l’ETP, dans ce cas on aboutira un déficit agricole (Da).
[Link].2. Méthode de L. Turc

Turc a établi une formule qui permet de calculer l’ETR annuelle.

.5
ETR = P / (0.9+(P /L ))

Où :

ETR : évapotranspiration réelle moyenne annuelle (mm);

P : précipitation moyenne annuelle (mm);

L : paramètre dépendant de la température

L = 300+250t+0.05t3

t : température moyenne annuelle (°C)

Application :

P = 372.25, t =14.33°C , L= 869.43

ETR = 357.65mm, soit 96.07% des précipitations.

[Link].3. Méthode de Coutagne

Coutagne a défini une formule empirique permettant de calculer l'ETR qui est la suivante:
ETR= P-λP
Où :

- -
P: précipitation moyenne annuelle (mm);
λ: coefficient dépendant de la température λ = (0.8+0.14t)
Où
t: température moyenne annuelle (°C).
Cette formule n'est applicable qu'a la condition suivante :

8λ ≤ P ≤ λ d’où la validité de la formule de Coutagne .

-Application p=372.25mm, t= 15.43°C
λ= .3
.3 ≤ P≤ . 7
Donc : ETR= P-λP = .37- .3 ( .37) .
ETR = 323.45mm, soit 86.89% des précipitations.
Les valeurs de l'évapotranspiration réelle calculée par les différentes méthodes sont
récapitulées dans le tableau 8.

Tableau 8. Calcul de l'ETR par les différentes méthodes

Méthodes ETR (mm) % par rapport des

précipitations
Thornthwaite .10.10 00.27
L. Turc 357. 5 06.21
A .Coutagne .0..11 8 .8

4.1.3. Estimation du ruissellement et de l'infiltration

La lame d'eau ruisselée sera estimée à l’aide de la formule de Tixeront – Berkaloff qui
s'écrit de la manière suivante :
R = P3/3 (ETP)
avec:
R : ruissellement moyen annuel (mm);
P : précipitation moyenne annuelle (mm);
ETP : évapotranspiration potentielle moyenne annuelle (mm).
On a : P= 372.25mm, ETP = 843.85mm; ce qui donne un ruissellement de :
R= . 5mm, soit .5 % des précipitations d’où une infiltration :
I = 5.58mm, soit 1.50% des précipitations
Les différents paramètres du bilan sont récapitulés dans le tableau 9 et Fig. 19.

- -
 Tableau 9. Bilan hydrologique de la station de Hammam Grouz (1988-2004) selon la
relation de Thornthwaite.
Mois Sep Oct Nov Déc Jan Fév Mar Avr Mai Jui Juil Aoû TOTAL
Paramètre
t (°c)
ETP
corrigée
(mm)
P (mm)
VR - - - - - - - -

R.F.U (mm) - - - - -

ETR (mm)
EX (mm) - - - - - - - - -

DE (mm) - - - - - -

Figure 19. Bilan hydrique simplifié de la station de Hammam Grouz (1988/1989 -

- -
5. Analyse morphométrique

5.1. Introduction

Le sous bassin du Haut Rhumel est l’un des sept sous bassins du grand bassin versant
Kebir Rhumel, alimentant le barrage du Hammam Grouz à Oued El Athmania (fig.20).
Le sous bassin couvre une superficie de 1120 Km et il est limité :
- Au Nord les monts de Benyahia Abderrahmane et Bouhatem.
- Au Sud par les monts de Oueled khlouf et Ltaya.
- A l’Ouest par Dj Ouled Hammou et dj Tavoudach à coté de Bir ElArch .
- A l’Est par Dj Grouz et Dj Tikouya (Coté de Chelghoum Laid).
Ces principaux affluent sont O. Merdj el kebir, O. Boumrah, [Link],[Link]
,[Link], [Link] , [Link], constituant l’amont de oued Rhumel, le cours d’eau
principal.
Le sous bassin-versant peut-être décomposé en quatre zones hydrologiques distinctes :
-Le Rhumel-amont (de Bellaa à Chelghoum Laid) avec Oued El Mehri, comme affluent.
Il draine près de la moitié de l’aire du bassin au barrage ;
- L’Oued Dekri constitue, dans la partie Nord du bassin, le plus important affluent par ses
capacités d’écoulement car charriant les eaux de ruissellement, issues des collines
relativement arrosées du piémont Sud-tellien ;
-Le sous bassin, constitué par les affluents de la zone Sud-Est qui s’écoulent à partir
d’une ligne de djebels Chebka-Tikouya;
- La zone hydrographique « inactive », au Sud de Tadjenanet (piémont du djebel Tafrent).
En raison du double effet de l’évaporation et de l’infiltration, l’écoulement sur ces glacis
fortement encroûtés atteint difficilement le cours principal. Le bassin du haut Rhumel est
soumis en grande partie au climat semi-aride des Hautes Plaines.

- -
 Figure 20. Bloc diagramme du bassin versant de Tadjenanet – Chelghoum laid

Variation de l’altitude
L’étude des classes d’altitude par tranches de m (Tab. 10) a permis de définir une
superficie totale du bassin versant de l’amont de oued Rhumel d’ordre km , dont le
périmètre est d’ordre de km.
La répartition des altitudes du bassin versant étudié, nous a permis de réaliser la courbe
hypsométrique et le diagramme de fréquence altimétrique en fonction de la surface.
Fig.21.
A partir de la courbe, nous pouvons déterminer :
Altitude médiane : 880m
Altitude moyenne : 898.45m
Dénivelée simple (H5% - H95%) : 313 m

- -
Tableau 10. Répartition altimétrique du bassin versant de Tadjenanet – Chelghoum
Laid
Classe d’altitude Superficie Superficie Superficie Superficie
(m) (Km (%) cumulée (km ) cumulées (%)
<
-
-
-
-
-
-

Figure 21. Courbe hypsométrique du sous bassin de Tadjenanet – Chelghoum Laid

5.3. Forme du bassin

La forme du bassin se détermine par l’indice de compacité de Gravelius donnée par la

formule suivante :
Kc = 0.28.P/S
où P et S sont respectivement, le périmètre et la surface du bassin.
Pour notre bassin Kc = 1.17, donc Kc > 1.12 autrement dit le bassin est allongé.
5.4. Densité de drainage
La densité de drainage d’un bassin donné est définie comme la longueur totale des cours
d’eau de tout ordre sur la surface drainée.

- -
 Dd = L/ A
où L : longueur totale cumulées des drains permanents et temporaire (km) ;
A : surface du bassin (km ).
On distingue la densité de drainage totale, permanant et temporaire.
Pour notre bassin Dd = .7 km/km dont 0.61 km/km de densité de drainage
temporaire.
Ce qui explique que la majorité des cours d’eau sont généralement à sec à l’ exception
des périodes de crues à la faveur des pluies torrentielles.
5.5. Profil en long
D’après le profil tracé le long du cours d’eau principal, le Rhumel (fig. ), il prend
naissance à partir des monts de Bellaa sur une altitude de plus de 1100 m avec une
direction perpendiculaire à la surface topographique jusqu’à son confluence avec oued El
Mehri à l’aval de Tadjenanet, puis il change de direction au niveau de Draa Mastaoua au
centre de la plaine vers l’Est, où il rencontre oued Dekri à l’aval de Chelghoum Laid.
Depuis ce point, il maintient un parcours plus au moins rectiligne jusqu’à son
engagement dans le barrage de Hammam Grouz à une altitude de 740 m.

Figure Profil en long de l’amont de Oued Rhumel et ses affluents

6. Conclusion
L’étude climatologique de la région étudiée montre que le climat dominant est
typiquement semi-aride, de type continental caractérisé par une irrégularité annuelle et
mensuelle. L’analyse morphométrique de notre bassin versant montre que le bassin est
allongé avec la dominance d’une zones d’altitude comprise entre 8 et m, dont la

- -
densité de drainage totale est d’ordre de .7 km km . L’établissement du bilan hydrique
selon la formule de Thornthwaite sur une période de seize ans (1988/89 -2003/2004) à
la station de Hammam Grouz, montre que cette région est déficitaire, avec
respectivement, une moyen annuelle des précipitations et de température de l’ordre de
37 . 5 mm et 5. 3°C. L’évapotranspiration réelle représente % des précipitations,
alors que l’infiltration est très faible et ne représente que .5% des précipitations
(5.58 mm). Le bilan est caractérisé par un déficit d’écoulement qui s’étend du mois de
juillet au mois d’octobre, avec une moyenne annuelle de l’ordre 5 mm et par un
excédent égal à 29 mm, soit 8% des précipitations.
Ceci laisse supposer que l’alimentation de la nappe ne se fait pas uniquement par
infiltration efficace, mais qu’une bonne partie de l’eau d’alimentation provienne par les
formations carbonatées bordent la nappe.

- -
 CHAPITRE III

GEOMETRIE DU COMPLEXE AQUIFERE

 GEOMETRIE DU COMPLEXE AQUIFERE
1. Introduction
La prospection électrique permet de quantifier l'effet produit par un courant électrique
traversant le sous-sol. Il existe une grande variété de techniques utilisant les méthodes
électriques. Par ailleurs, dans un domaine envisagé, il existe de nombreuses techniques de
mesure. Les méthodes fondées sur la mesure du paramètre " résistivité ", (imaginées en
l9l2 par les frères Schlumberger), sont actuellement les plus répandues, plus développées
et les plus diversifiées.
La méthode utilisée dans le cas de Tadjenanet-Chelghoum Laid, est celle des sondages
électriques (S.E) de type Schlumberger. Une campagne de géophysique par SE a été
effectuée par la société Algérienne de géophysique (ALGEO) en 1977 dans la région de
Tadjenanet-Chelghoum Laid. Le but était l’exploration des ressources en eau pour les
besoins du complexe détergent (S.N.I.C.), localisé dans la zone industrielle de
Chelghoum Laid. Une étude géophysique antérieure a été réalisée par la compagne
générale de géophysique (CGG) en 1973 qui a concerné les hautes plaines sétifiennes et
qui a touché la région de Tadjenanet, afin d’étudier les possibilités aquifères de la région.
Pour caractériser le sous-sol par la mesure de la résistivité, ce paramètre peut varier de la
manière suivante : (Astier, 1971 ; Boudoukha, 2008)
- de 1 à quelques dizaines d'ohm.m pour les argiles et les marnes ;
- d'une dizaine à quelques centaines d'ohm.m pour les sables et les grès
marneux ;
- d'une centaine à plusieurs milliers d'ohm.m pour les calcaires et les roches
éruptives.

Situation de l’étude
La zone prospectée par ALGEO couvre 600 Km et s’étend de part et d’autre de
Chelghoum Laid (fig. 23).
- Vers l’Ouest, elle s’étend jusqu'à la limite administrative de la Wilaya de Mila
en passant par Tadjenanet ;
- Vers l’Est, elle s’étend jusqu'à Hammam Grouz.
La partie couverte par l’étude de la CGG en 73, elle s’étend de l’Ouldja à Bir el Arch à
la limite ouest du terrain étudié.

- -
3. Répartition des sondages électriques
Dans la zone étudiée, les sondages électriques ont été disposés suivant des profils
de direction Nord-Sud, espacés de km. Ces sondages étant distants d’environ km sur
ces profils. Les S.E., qui se trouvent dans la zone étudiée, forment les profils A, B, C, D,
E, F, G, H et I (ALGEO. 77) et A’, B’, C’ et D’ (CGG. 73). Fig. 3.

335 Dj. G
r ouz
ss
De
d H
.e
330 Dj J
F
f1
E
D
325
C
B c1 Chelghoum Laid
A
320

Tadjenanet
315

D' a
310 uy
iko
.T
B' Dj
C' A'
305
780 785 790 795 800 805 810 815 820 825 830

Légende Echèlle
A,B,C...Profils géoélectriques
. 0 2 4 km

. Sondagr électrique (ALGEO.1977)

Sondagr électrique (CGG.1973)
Oued
Ville
Limite du terrain d'étude

Figure 23. Carte de position des sondages électriques dans la plaine

de Tadjnanet-Chelghoum Laid.

4. Interprétation
4.1. Examen de quelques familles de S.E
Une courbe de sondage électrique représente la variation de la résistivité en fonction de la
demi-longueur des lignes AB. Elle traduit donc les variations de la résistivité apparente
selon la profondeur d’investigation. L’interprétation de cette courbe, permet de
déterminer les paramètres géoélectriques de la couche traversée par le courant. Les
difficultés de l’interprétation des données géophysiques sont liées à l’absence de forages
pouvant servir d’étalon, pour faire correspondre les résistivités aux terrains géologiques
présents.

- -
 Sondage électrique à l’extrême-Ouest, coupe A
Parmi les courbes de sondages électriques qui présentent une forme complexe, celle du
SE A5. Cette courbe montre la présence de cinq terrains où alternent des couches
résistantes et couches conductrices : 100, 10,50 et 40 ohm.m. Fig. 24. La couche
caractérisée par une résistivité de 50 ohm.m, correspondant à des calcaires lacustres ayant
une épaisseur de 23 m se trouvant à une profondeur de 40 m environ.

Figure 24. Courbe du sondage électrique A5.

4.1.2. Sondage électrique étalon B3

Un SE étalon en ligne AB= m, a été réalisé à côté d’un forage d’une profondeur
de 278 m dans la région de Chabet Rmali à l’Ouest de Tadjenanet aux coordonnées
suivantes : X = 800 km, Y = 317 km.
La courbe a une forme du fond du bateau correspond électriquement à un trois terrains.
L’analyse du log stratigraphique du forage, permet de voir qu’on a un empilement de six
terrains mais leurs natures ne leur permettent pas de les distinguer électriquement. Fig.
5. Le calage de ces données permet d’avoir les résistivités électriques des couches
correspondantes. Tableau 11.

- -
 Figure 25. Courbe du sondage électrique d'étalonnage B3.

Tableau 11. Calage du SE B3 sur les données du forage de Chabet Rmali

Résistivité ohm.m Lithologie Epaisseur (m) Profondeur (m)
7 alluvions
5 Marne et Calcaire lacustre 3
8 Argile sableuse 33
5 Calcaire lacustre 7
Marne
3 calcaire

4.2. Cartes de résistivité

Lors de cette compagne, les lignes d’injection AB étaient de m, m et 3

m ce qui correspond à des profondeurs d’investigation moyennes respectives égales à 35
m, 65 m et 500 m. Ces différentes lignes AB ont permis de tracer les cartes de résistivité
correspondantes.
4.2.1. Carte de résistivité AB = 200 m
Les résistivités égales ou supérieures à 150 ohm.m forment une tache au centre, au Sud et
au Sud-Ouest. Cette tranche de résistivité correspond aux calcaires marneux du Crétacé et

- -
de l’Eocène. Par ailleurs les valeurs plus faibles (voisines de 20 ohm.m), correspondent
aux formations du Mio-Plio-Quaternaire (argile, marne, sable et alluvions). Fig. 26.
335 Dj. G
rouz

Ohm.m
s
es
D
ed
330 j.
D 140

120
325

100

320
80

60
315

40
a
310 uy
i ko
.T 20
Dj
Dj. T
n outit
305 0
780 785 790 795 800 805 810 815 820 825 830

Echèlle

0 2 4 km

Figure 26. Carte de résistivité AB = 200 m

4.2.2. Carte de résistivité AB = 400 m

Cette carte permet de distinguer une répartition des plages de résistivité semblable à cette
de AB= m. Il s’agit toujours des calcaires marneux de Crétacé et de l’Eocène. Le
reste du terrain est occupé par une formation plus au moins conductrice caractérisant les
formations du Mio-Plio-Quaternaire, de résistivité inférieure à 20 ohm.m. Fig. 27.
335 Dj. G
r
Ohm.m

s ouz
es
D
ed
330 j. 180
D
160

325 140

120

320
100

80
315
60

ya 40
310 ou
Tik
. 20
Dj
Dj. T
n outit
305 0
780 785 790 795 800 805 810 815 820 825 830

Echèlle

0 2 4 km

Figure.27 Carte de résistivité AB = 400 m

- -
4.2.3. Carte de résistivité AB = 3000 m
L’examen de cette carte montre la variation de la résistivité à une profondeur moyenne
de 500 m. A ce niveau, la résistivité moyenne est de l’ordre 5 ohm.m le long de Oued
Rhumel. On retrouve également cette plage de résistivité au Nord-Est et à l’Ouest de Dj.
Tikouya. Au Sud-Ouest, les calcaires du Barrémien du massif de dj. Tnoutit plongent
jusqu’à 5 m de profondeur où ils sont caractérisés par une résistivité de 5 ohm.m.
L’enracinement des dj. Grouz et Ed Dess, n’apparait pas en profondeur puisque la
région est caractérisée par de faibles valeurs de résistivité (<30 ohm.m). Cette carte
permet donc de délimiter en profondeur (5 m), les zones où la recherche d’eau est à
éliminer. Fig. 28.
335 Dj. G
r

Ohm.m
s ouz
es
D
ed
330 j.
D 120

110

325 100

90

80
320
70

60
315
50

40
a
310 uy
i ko 30
.T
Dj 20
Dj. T
n out i t 10
305
780 785 790 795 800 805 810 815 820 825 830

Echèlle

0 2 4 km

Figure.28 Carte de résistivité AB = 3000 m

4.3. Interprétation des coupes géoélectriques

Les coupes sont réalisées pour suivre l’extension latérale des différents horizons
aquifères selon la résistivité. Plusieurs coupes géoélectriques ont été tracées selon les
différents profils. Fig. 29 et 30. L’examen de ces coupes permet de tirer plusieurs
remarques.
-Coupe A :
Le niveau, entre 20 et 60 m de profondeur et de 50 ohm.m est bien visible, correspond
aux formations d’argiles sableuses et calcaires lacustres du Pliocène. Vers le Nord il
affleure une formation très conductrice de trias.

- -
-Coupe C et D :
Sur la coupe C, l’Eocène s’étend entre C et C11. Les couches conductrices sont
abondantes. Sur D, le sous ensemble Nord est bien distinct du sous ensemble Sud. Un
niveau résistant (d’après CGG 73), il se développerait vers le Sud et se subdivise en deux
niveaux, l’un résistant qui pourrait appartenir au Crétacé supérieur (Maestrichtien) et
l’autre conducteur qui correspond une formation triasique
-Coupe E :
La coupe E montre bien la formation de l’Eocène, avec des épaisseurs de à
mètres, dont la résistivité moyenne est de l’ordre de ohm.m. Le substratum remonte
vers le sud, mais reste profond.
-Coupe G :
La coupe représente des niveaux conducteurs (marnes et argiles). Le substratum est
relativement peu profond contre le flanc Sud-Est de dj. Grouz (ALGEO 1977).
Nous signalons l’affleurement des masses triasiques localisées dans les monts de Dj. Ed
Dess au nord-ouest et Dj Tikouya au Sud-est.
L’examen des coupes D’ et E’ localisés dans la partie sud-ouest de Tadjenanet
précisément le nord-ouest de El Oueldja , nous a permis de constater :
- La présence d’une masse très conductrice d’une résistivité inférieure à .m. Ces
masses conductrices sont très épaisses et correspondent à des formations marneuses ou
argileuses liées au Trais.
- La présence des masses résistantes épaisses qui ont été attribuées en général au Crétacé
inférieur (calcaires et dolomies). Ces masses résistantes (100 à 300 .m) surmontent
parfois une formation conductrice.

- -
Figure 29. Coupes géoélectriques A, C, D et E

- -
 Figure Coupes géoélectriques B’ ,C’,D’ et E’

D’après la morphologie du substratum, on a des formes synclinales, avec un empilement

de plusieurs couches successivement conductrices et résistantes. On distingue également

- -
la présence de contacts anormaux dus probablement à la présence côte à côte de
formations conductrices et de formations résistantes. La présence d’une formation
résistante dont la résistivité moyenne varie entre 20 et 50 .m, avec une épaisseur de
l’ordre de à 60 m attribuée à la formation de calcaires lacustres du Mio-Pliocène. Cette
formation est tabulaire et elle a un substratum conducteur.
5. Carte du substratum
L’interprétation des coupes géoélectriques permet de renseigner sur la
morphologie du substratum de l’aquifère des calcaires lacustres ce qui nous a permis de
tracer une carte de substratum. Fig. 31. Le substratum étudié est formé par une couche
conductrice dont la résistivité varie entre 8 et 10 .m correspondant à des argiles,
identiques sur toute la plaine.
Ainsi la carte du substratum montre que ce dernier remonte au partie ouest du
terrain près de Djebel Tnoutit et Djebel Ed Dess . Le substratum devient plus profond en
allant vers l’Est (d’environ m de profondeur). Il présente une légère convexité au
centre du terrain. Donc nous sommes en présence d’un substratum ayant une morphologie
irrégulière qui épouse approximativement la topographie.
335
Dj. G
rouz (m)
ss
De
d 860
.e
330 Dj
820

325 800

780
320
760

315 740

720
a
310 uy
ko
. Ti 700
Dj
690
305
780 785 790 795 800 805 810 815 820 825 830

Echèlle
0 2 4 km
Figure 31. Carte du substratum
6. Carte de la résistance transversale

Pour étudier la potentialité de l’aquifère superficiel, nous avons établis une carte de la
résistance transversale. Fig. 32. Cette dernière montre la variabilité de la résistivité de la
couche en fonction de son épaisseur. A partir de cette carte, on remarque qu’au centre et à
l’Est, une grande résistance transversale (8 m.m) donc de bonnes caractéristiques
hydrogéologiques. Par contre au Nord et au Sud de Tadjenanet, ainsi que sur le flanc Sud

- -
de Dj Grouz, on a de faibles valeurs de la résistance transversale (<2000 m.m),

traduisant la présence d’une formation de mauvaises caractéristiques hydrogéologiques.

335 Dj. G

2
rou

Ohm.m
s z
es
D
ed
330 j.
D
8000

7000
325
6000

320 5000

4000

315
3000

2000
ya
310 ou
Tik
. 1000
Dj
Dj. T
nout
305 it 0
780 785 790 795 800 805 810 815 820 825 830

Echèlle

0 2 4 km
Figure.32 Carte de la résistance transversale

7. Conséquences hydrogéologiques et conclusion

L’étude géophysique de la région nous amène aux conclusions suivantes :

De grandes zones peuvent être éliminées, du fait de l’absence d’horizon intéressants
(calcaires ou autre) et du fait de la trop grande profondeur du substratum.
- zone Sud
- zone Nord, excepté le début du profil D
- zone comprise entre l’Eocène ‘’central’’ et Dj. Grouz (quelques S.E, toujours
isolés, indiquent des bancs ‘’résistants’’ de faible étendue).
- La zone extrême Sud est contaminée par la présence du Trias triasique.
La zone extrême Ouest vers le profil A, elle présente un banc à 50 ohm.m peu profond
dont l’intérêt, s’il est aquifère, ne peut être que local. Le flanc Sud-Est de Dj Grouz offre
une situation qui pourrait être favorable de même que pour les calcaires de l’éocène ‘’
central ‘’ qui ne semblent pas aquifères.
Les cartes établis selon les différentes profondeurs, nous ont permis de constater la
présence de deux grandes formations géologiques. L’une a une résistivité moyenne
inférieur à 50 Ohm.m correspondant aux formations du Mio-Plio-Quaternaire (argiles,
alluvions, croutes calcaire, calcaire lacustre, alluvions) se trouvant sur la totalité du
terrain (plaine et terrasse de oued Rhumel), l’autre plus au moins résistante (résistivité

- -
sup. à 50 Ohm.m) caractérisant les formations carbonatées appartenant aux formations de
l’éocène et du Crétacé. Ces deux formations peuvent avoir un intérêt hydrogéologique.
Le substratum de l’aquifère est constitué par une formation conductrice (argile et
marnes), ayant une morphologie irrégulière, ce qui influencera sans doute le
comportement hydrodynamique de l’aquifère.

- -
CHAPITRE VI

HYDROGEOLOGIE
 HYDROGEOLOGIE

1. Introduction

Les différentes études géologiques, géophysiques et hydrogéologiques menées dans la

région (KHEDIDJA, 2001 ; CGG, 1973 ; BENSOUILAH, 1995) ont montré l'existence
de trois formations géologiques perméables permettant la formation de trois aquifères.
- Un aquifère dont la perméabilité est liée à un réseau de fissure caractérisé par une
circulation karstique dans les formations calcaires néritiques du Crétacé et de
l’ensemble allochtone su d Sétifien. (KHEDIDJA. 2001 ; ISSAADI. 8 ).
- Un aquifère superficielle dans les alluvions du Quaternaire développé au niveau
d’Oued Rhumel qui est en phase de déstockage.
- Un aquifère continu dans les formations fluvio- lacustres d’âge Mio-Pliocène.
Le premier s’étend sur toute la partie Nord Est de la région dans le massif carbonaté du
Crétacé de Djebel Grouz et la partie Sud-Ouest, dans le massif carbonaté de Djebel
Tnoutit. On retrouve également ce type d’aquifère dans les calcaires de l’Eocène qui
affleurent dans la partie centrale de la plaine, au Nord et au Sud.
Le deuxième aquifère est peu épais et il est étroitement lié au cours d’eau. En période
d'étiage, il alimente le cours d’eau, par contre en période de hautes eaux c’est l’inverse
qui se produit. Le long de Oued Rhumel, cette nappe est en partie alimentée par l’aquifère
des calcaires néritiques du Crétacé. Toutefois, il est intéressant de noter qu’une partie de
cette nappe alluviale a une forte minéralisation suite à une intense évaporation. Son
niveau piézométrique est relativement proche de la surface du sol, de 0 à 3m (FARAH,
), il a donc peu d’intérêt du point de vue hydrogéologique.
Par contre le troisième aquifère se trouvant dans les formations du Mio-Pliocène couvrant
d’importantes superficies et il a de bonnes potentialités. L’épaisseur de l’aquifère varie de
100 à 150 m et il est constitué par des dépôts fluvio- lacustres offrant habituellement une
coloration rougeâtres assez prononcée. Ces formations affleurent beaucoup plus sur la
périphérie du bassin versant. Il s’agit des calcaires lacustres, de marne rouge, de limon et
de sable rouge.
Dans cette partie on ne s’intéressera qu’au deuxième et troisième aquifère.

- -
 Aquifère du Mio-Plio- Quaternaire
Description
Cette nappe de grande extension, occupe la quasi-totalité de la région d’étude, elle s’étend
de Bir El Arch à l’Ouest jusqu'à Oued Athmania à l’Est en passant par Tadjnanet et
Chelghoum Laid. Cet aquifère s’est formé dans les grés, graviers, sables, alluvions et
calcaires lacustres dans lesquels on y trouve des passages d’argile et de marnes. Le
substratum est formé par des argiles et marnes du miocène.
L’alimentation de cet aquifère se fait à partir des précipitations par infiltration surtout
durant les périodes de crue de oued Rhumel, et par une alimentation latérale à partir des
formations carbonatées se trouvant sur les limites au niveau des Dj Dess, Tnoutit et
Grouz. Nous signalons que cet aquifère est interrompu par des accidents tectoniques ou
complètement érodés laissant ainsi apparaitre les calcaires de l’Eocène au centre du
terrain, ce qui est confirmé par la prospection électrique. Ceci est très visible au bord de
la route nationale5 (RN 5) entre Tadjnanet et Chelghoum Laid. Cet aquifère joue un rôle
très important dans l’alimentation en eau potable des habitants de la région d’étude.
Géométrie de l’aquifère Mio-Plio-Quaternaire
Le dépouillement des log lithostratigraphiques de quelques forages mécaniques
implantés dans la plaine d’étude nous a permis d’établir une coupe hydrogéologique. Fig.
33.
L’examen de la coupe hydrogéologique de direction Sud-Ouest - Nord Est, établie à partir
des logs stratigraphiques des forages a mis en évidence la présence d’un remplissage
Mio-Plio-Quaternaire, constituée par des dépôts très variés correspondant aux alluvions,
grés ,sable et calcaire lacustres, d’épaisseur allant jusqu’à m, reposant sur un
substratum imperméable formé de marnes et argiles. L’aquifère superficielle devient mince
près de Oued Rhumel suite aux contacts anormaux, qui laissent apparaitre les formations de
l’Eocène au centre du terrain. À la partie sud-ouest, l’aquifère du Mio- Plio- Quaternaire est
en contact direct avec la formation carbonatée du Crétacé par la mise en place d’un accident
tectonique. Selon l’orientation de la coupe, l’écoulement des eaux souterraines convergent
vers Oued Rhumel.

- -
 Figure 33. Coupe hydrogéologique

3. Piézométrie
La piézométrie schématise le fonctionnement de l’aquifère et son comportement
hydraulique. Elle permet principalement :
*D’analyser la morphologie de la surface piézométrique;
*De définir la direction des axes principaux d’écoulement et des axes de
drainage ;
*De définir les conditions aux limites.
Inventaire des points d’eau
Durant le mois de Mai de l’année 7, nous avons effectué un inventaire de tous les
points d’eau (puits et forages), implantés dans la plaine (annexe. ). En parallèle, nous
avons mesuré les niveaux piézométriques dans les différents points d’eau. 7 points d’eau
ont été recensés. Fig. 3 . La forte densité des points d’eau, nous a permis d’adopter un
maillage de l’ordre du Km . Les puits réalisés sont pour la plupart des ouvrages récents,
forés à la main et soutenus par du parpaing, ayant un diamètre qui varie entre 1 m et 4 m.
Ces ouvrages ne sont pas couverts et équipés par des pompes électriques. Tous ces points
d’eau sont destinés à l’alimentation en eau potable et à l’irrigation des cultures
maraîchères. Nous signalons que la majorité des puits sont imparfaits.

- -
 335 Dj. G
rouz
ss
De
d
.e 31 29 27
330 Dj
32 30 28 26
Vers Constantine
2524 9
33 10
325 36 Chelghoum Laid8 76
44 35 34
37 5
42 40 RN.515 12 4
434139 38 23 16 14
22 11 3
Tadjenanet 46 13 2
45 21 1718 1
320 Vers Sétif 49 48
56 47 19
60 20
57 55 50
61 51
54
53 52
315 626559
63
64

67
a
310 uy
i ko
.T
Dj
305
780 785 790 795 800 805 810 815 820 825 830

Légende
. Points d'eau
0
Echèlle

2 4 km
----- Oued
Ville
Limite du terrain d'étude

Route
Figure Carte d’inventaire des points d’eau

Les données piézométriques mesurées au niveau des 33 points d’eau retenus, nous ont
permis de tracer :
- Deux cartes piézométriques la première au début de la période de basses eaux (Mai
2007) et la seconde à la fin de la période de basses eaux (Septembre2007) ;
- Une carte d’égale épaisseur de la couche mouillée ;
- Une carte de battement.
Conditions aux limites
L’analyse de la carte piézométrique (fig. 35) montre que la nappe a :
- Une limite à flux entrant à partir des massifs carbonatés qui encadrent la plaine au Nord,
au Sud-Ouest, et au Nord ;
- Une limite étanche au Sud-Est du terrain caractérisée par la présence des formations
imperméables évaporitiques et marneuses. Le substratum est formé par des argiles et des
marnes.
- Une dépression piézométrique confondue avec la cuvette du barrage
Interprétation de la carte piézométrique (Mai 2007)
L’étude de la carte piézométrique (fig.35), montre que la nappe a une surface
piézométrique irrégulière, sa profondeur par rapport au sol est variable. L’examen de
l’allure des courbes isopièzes permettent de caractériser une dépression piézométrique à
- -
la partie avale de Oued Rhumel suite à l’abaissement du substratum selon la direction
Ouest -Est.
L’aquifère étudié est caractérisée par un écoulement convergeant vers l’est
caractérisé par des courbes isopièzes fermées de direction Ouest – Est qui coïncide avec
la morphologie du substratum. Les courbes isopièzes sont très serrées à l’Ouest et à l’Est
traduisant un gradient hydraulique moyen de l’ordre de 3%, ce qui implique une faible
perméabilité. Au centre et à l’extrême Est du terrain, les courbes sont plus espacées
traduisant un gradient hydraulique plus faible de l’ordre de %. Ce faible gradient est dû à
la bonne caractéristique du terrain.
L’alimentation se fait en grande partie par les formations carbonatées du Crétacé
inférieur du Djebel Tnoutit et Djebel Grouz ainsi que les formations éocènes de Draa
Mestaoua et Djebel ed Dess.
335

(m)
Dj. G
rouz
ss 950
De 940
d
.e 930
330 Dj 920
910
900
325 890
880
870
860
850
320 840
830
820
810
315 800
790
780
770
310 ya
ou 760
Tik 750
.
Dj 740
730
305 720
780 785 790 795 800 805 810 815 820 825 830

Légende Echèlle
0 2 4 km

800 Courbe piézométrique (m)

Ligne de courant
Oued
Ville

Figure 35. Carte piézométrique. Mai 2007

4. Interprétation de la carte piézométrique (Sep 2007)

Afin de voir l’évolution de la piézométrie de la nappe dans le temps, on a établi
une deuxième carte piézométrique durant la fin de la période des basses eaux pendant
laquelle l’exploitation est intensive, sous des conditions climatiques très défavorables.
Fig. 36. La morphologie piézométrique de la nappe est identique à celle du mois de Mai

- -
2000, ceci permet de dire que le comportement hydrodynamique de la nappe est stable
dans le temps. La différence entre les deux périodes n’est visible qu’au niveau de la côte
piézométrique comme le montre la carte de battement.

335

(m)
Dj. G
r ouz
ss
De 930
d 920
.e
330 Dj 910
900
890
880
325
870
860
850
320 840
830
820
810
315 800
790
780
770
310 ya 760
ou
Tik 750
Dj. .
Tn Dj 740
ou 730
tit
305 720
780 785 790 795 800 805 810 815 820 825 830

Légende Echèlle
0 2 4 km
800 Courbe piézométrique (m)
Ligne de courant
Oued
Ville

Figure .36 Carte piézométrique. Septembre 2007

5. Battement de la nappe

La confrontation des deux cartes piézométriques (Mai et Septembre 2007), nous a permis
de tracer la carte de battement de l’aquifère. Fig. 37. Selon la carte établie on a enregistré
un abaissement de la nappe de 1 à 8 m. Cet abaissement est dû à un excès de pompage,
sous un régime non influencé. Les fortes valeurs de battement ont été enregistrées au
centre du terrain ( à 8 m). Ceci est dû à la présence d’un nombre élevé de points d’eau
implantés dans cette zone. Par contre le reste du terrain n’a enregistré qu’un faible
battement qui varie entre et m surtout à la périphérie de la nappe là où l’alimentation
latérale est bien marquée.

- -
 335

(m)
Dj . G
rouz 8
ss
De 7.5
d
.e
330 Dj 7
6.5
6
325
5.5
5
4.5
320
4
3.5
3
315
2.5
2
a 1.5
310 o uy
ik 1
.T
Dj 0.5

305 0
780 785 790 795 800 805 810 815 820 825 830

Légende Echèlle
0 2 4 km
2 Courbe d'égale battement de la nappe (m)
Ville

Figure 37. Carte de battement de la nappe

6. Evolution de l’état piézométrique

Afin de suivre l’évolution spatio-temporelle du niveau piézométrique de la nappe

alluviale dans région d’étude, une campagne piézométrique a été réalisée durant la
période de Mai 2013.
6.1. Interprétation de la carte piézométrique (Mai 2013)

L’examen de la carte piézométrique du mois de Mai 3 (fig.38) nous a permis de

constater que la nappe alluviale garde la même structure des courbes isopièzes avec une
dépression piézométrique moins importante par rapport à celle de la période de Mai 2007.
Ceci s’est traduit par de faibles rabattement dans sa partie Est. Le gradient hydraulique est
devenue très important au centre, de l’ordre de % du soit à un fort écoulement ou à une
surexploitation de la nappe à cet endroit.

- -
 335

(m)
Dj. G
rouz
ss 930
De 920
d
.e
330 Dj 910
900
890
880
325 870
860
850
840
320 830
820
810
800
315 790
780
770
760
310 ya 750
ou 740
Tik
. 730
Dj 720
710
305
780 785 790 795 800 805 810 815 820 825 830

Légende Echèlle
0 2 4 km

800 Courbe piézométrique (m) Profil d'évolution

Ligne de courant
Oued
Ville

Figure 38. Carte piézométrique. Mai 2013

6.2. Evolution spatio-temporelle de l’état piézométrique de l’aquifère

L’allure générale des isopièzes entre la période de Mai 7 et celle de Mai 3 reste
presque identique, l’écoulement se fait généralement de l’Ouest vers l’Est. L’axe
principal d’écoulement convergent vers l’aval et coïncide avec la cuvette du barrage de
Hammam Grouz. Fig. 3 . L’examen de ce diagramme a permis de noter que la nappe est
influencée par la surexploitation à l’aide de pompage et la diminution de la pluviométrie
durant ces dernières années.
La nappe alluviale présente des forts abaissements du niveau piézométriques de l’ordre de
10 m dans sa partie centrale tandis que des faibles abaissements sont enregistrés dans la
partie Est, estimée de 2m.

- -
Figure 39. Evolution spatio- temporelle de la piézométrie selon le sens d’écoulement

Carte d’égale épaisseur de la couche mouillée

La superposition de la carte du substratum et la carte piézométrique (Septembre 2007),

nous a permis de tracer une carte d’égale épaisseur de la couche mouillée de l’aquifère.
Fig. 40. Cette carte montre que la zone la plus favorable à l’implantation d’éventuels
futurs forages se trouve à l’Ouest et à l’Est du terrain, dont l’épaisseur est atteint de 8 m.
Dans les régions centre et sud, l’aquifère possède une épaisseur très faible (inférieur à
m), ce qui implique l’irrégularité de la morphologie du substratum.
335 Dj
.G (m)
ss ro
uz
De 90
330 . ed
Dj
80

325 70

60

320 50

40
315
30
ya
ik ou 20
.T
310 D
j.T Dj
no 10
ut
it
305 0
780 785 790 795 800 805 810 815 820 825 830

Légende Echèlle
0 2 4 km
20 Courbe d'égale epaisseur de la couhe mouillée (m)

Ville

Figure Carte d’égale épaisseur de la couche mouillée

- -
8. Les caractéristiques hydrodynamiques
Pour mieux comprendre le comportement hydrodynamique de l’aquifère, une série
d’essais de pompages a été réalisée sur six puits. Le choix des points d’eau s’est fait selon
le type d’équipement et sa position dans la plaine. La durée de pompage varie entre 1
heure 40 mn et 3 heures 45 mn pour la descente et de 2 heures 15 mn à 3 heures pour la
remontée.
Interprétation des données
Lors du pompage on s’est aperçu que le régime permanent n’a pas pu être atteint et que
les puits tarissent au cours du pompage car, le débit de la pompe est beaucoup plus
supérieur au débit d’alimentation, ce qui nous a donc obligé d’utiliser uniquement les
données relatives à la remontée selon le modèle de [Link].
Q
s= W (u) avec :
T W (u) : fonction exponentielle de u
u : temps réduit

X S
u=
4T t
Où :
s : rabattement (m)
Q : débit de pompage (m /s)
T : transmissivité (m s)
t : temps écoulé à un instant donné depuis le début de pompage (s)
S : coefficient d’emmagasinement
X : distance entre le puits et le piézomètre m).

.2. Courbe de remontée

La courbe de la remontée est obtenue en portant sur un diagramme bi-

logarithmique, la fonction s’ = f (tp t’)

Où s’ : rabattement résiduel (m) ;

tp : temps de pompage (s) ;
t’ : temps de remontée (s).

Les courbes de ces essais de pompage (fig.41) montrent que le rabattement résiduel
diminue en fonction du temps. Au début, la remontée est rapide, puis elle décroît
progressivement

- -
 Puits n° 09 Puits n° 11
3,0 0,9
2,6 0,7

Rabattement résiduelle (m)

2,2
0,5
Rabattement résiduelle (m)

1,8
1,4 0,3

1,0

0,1
0,6

0,2
1 4 7 10 40 70 100 400 1 3 5 7 9 20 40 60 80 100 300

Temps (s) Temps (s)

Puits n°23
Puits n ° 31
5,00

3,00
1,00
0,80
Rabattement résiduelle (m)

Rabattement résiduelle (m)

0,60
1,00
0,40 0,80
0,60

0,20 0,40

0,20
0,09
0,07
0,09
0,05 0,07
0,05
0,03
0,03
5 15 25 35 1 3 5 7 9 20 40 60 80 100 300

Temps (s ) Temps (s)

Puits n°50
puits n°43 5,0
4,5
2,00 4,0
3,5
3,0
2,5
Rabattement résiduelle (m)

0,90 2,0
Rabattement résiduelle (m)

0,70
1,5
0,50
1,0
0,30

0,5

0,10
0,08
0,06

0,04

3 5 7 9 20 40 60 80 100 1 3 5 7 9 20 40 60 80 100 300

Temps (s) Temps (s)

Puits n°57 Puits n°67

0,7

0,5
1,00
Rabattement résiduelle (m)

Rabattement résiduelle (m)

0,3 0,80
0,60

0,40

0,1 0,20

0,09
0,07

0,05
1 3 5 7 9 20 40 60 80 100 1 3 5 7 9 20 40 60 80 100

Temps (s) Temps (s)

Figure 41. Essai de pompage – courbes de remontée

- -
 L’application du modèle de Thiès a permis la détermination de la transmissivité pour
l’ensemble des essais effectues. Khedidja, 7. Les résultats de cette interprétation ont
été représentés dans le tableau 12.

Tableau 12. Paramètres hydrodynamiques de l’aquifère superficiel

Epaisseur
Transmissivité Perméabilité
N° du point X (km) Y (km) Z (m) - utile -
m s) m/s)
d’eau (m)

Répartition de la transmissivité et de la perméabilité

9.1. Carte de transmissivité
Le résultat des essais de pompage a permis de tracer une carte des transmissivités de la
région. Fig. . Cette carte montre que l’aquifère est caractérisé par une transmissivité
-
moyenne de l’ordre de .8. m /s dans la partie médiane et nord-est du terrain, ce qui
implique probablement la présence des passages argileux. Par contre à l’Ouest et au Sud-
-
Ouest, on a une transmissivité relativement plus élevée, de l’ordre de . m /s. Ceci est
du à la présence d’un matériel plus propre lié aux formations carbonatées de Djebel
Tnoutit et Djebel Ed Dess.,

- -
 335 Dj -4 2
.G En 10 m /s
ss ro
uz
De
2.6
d 2.5
330 .e
Dj 2.4
2.3
2.2
325 2.1
2
1.9
1.8
320
1.7
1.6
1.5
315 1.4
1.3
a
uy 1.2
Tiko 1.1
310 D .
j. Dj 1
Tn 0.9
ou
tit 0.8
305 0.7
780 785 790 795 800 805 810 815 820 825 830

Légende
Echèlle
-4 2 0 2 4 km
1 Courbe d'égale transmissivité ( en 10 m/ s )
Ville

Figure 42. Carte de transmissivité en 10- m /s

9.2. Carte de perméabilité

La carte de perméabilité (fig.43) montre une augmentation de ce paramètre du Nord vers
-
le Sud du terrain, avec une valeur maximale de l’ordre de . . m /s au Sud-ouest de
-
Tadjenanet, et une valeur minimale de l’ordre de . m /s au Nord du terrain au niveau
du puits n° 9 situé au nord-est de Chelghoum Laid caractérisant la présence d’un
matériau plus au moins fin.
335
Dj -4
En 10 m /s
.G
ss ro
uz
De 31 27 1.22
d
330 .e 28
Dj 32 1.2

9 1.1
33
325 36 76 1
44 3735 34
40
42 39
43 23 0.9
38 22 11
46 2
21 18 0.8
320 49 48
0.7
57 50
0.6
53
62
63
315 0.5

0.4
67 ya
ik ou 0.3
310 D . T
j. Dj 0.2
Tn
ou
tit 0.1

305 0.02
780 785 790 795 800 805 810 815 820 825 830

Echèlle
Légende 0 2 4 km
-4
2 Courbe d'égale perméabilité ( en 10 m / s )
Ville Oued

Figure 43. Carte de perméabilité en 10- m /s

- -
10. Conclusion
L’étude hydrogéologique de la région de Tadjenanet - Chelghoum Laid nous a permis de
dire que :
- Le remplissage du Mio-Plio-Quaternaire est la formation aquifère la plus exploitée
surtout par les agriculteurs (alluvions, sable et calcaire lacustre).
- L’aquifère est capté par un nombre élevé de puits et forages qui se répartissent sur
l’ensemble du système aquifère, ce qui a conduit à l’abaissement du niveau piézométrique
d’une manière alarmante atteint les 10 m.
- La morphologie piézométrique est plus au moins identique à celle de la surface
topographique et elle du substratum
- L’alimentation de la nappe se fait essentiellement par les limites par déversement des
eaux des formations Carbonatées limitrophes du Crétacé et de l’Eocène.
- L’aquifère présente des caractéristiques hydrodynamiques très variables, avec un
maximum sur la limite sud-ouest. Les zones ayant les meilleures caractéristiques
hydrodynamiques ne correspondent pas forcement aux zones ayant le maximum de
l’épaisseur utile.

- -
 MODELISATION

1. Introduction

La mise en œuvre d’un modèle nécessite la compréhension du fonctionnement

des aquifères et donc, requiert au préalable une description géométrique du
système hydrogéologique ; en particulier les unités hydrogéologiques dont il est
formé. Ceci nécessite obligatoirement une connaissance suffisante et
relativement complète du milieu physique et du système d’écoulement. La
réalisation d’un modèle numérique passe de ce fait par les étapes suivantes
(fig.44):
a. élaboration d’un modèle conceptuel,
b. construction de la géométrie du modèle,
c. modélisation des écoulements.
Le modèle conceptuel correspond à une hypothèse sur le mode de
fonctionnement d’un système ou d’un processus. Cette hypothèse peut être
exprimée quantitativement sous la forme d’un modèle mathématique. Les
modèles mathématiques sont des modèles abstraits qui assimilent les processus
à des équations, les propriétés physiques à des constantes ou des coefficients
intervenant dans ces équations, et les mesures de l’état du système à des
variables. La plupart des modèles hydrogéologiques utilisés aujourd’hui sont des
modèles mathématiques déterministes. Ils sont basés sur les principes de la
conservation de masse, de la force d’impulsion et de l’énergie et décrivent des
relations de causes à effets.

Pour toute modélisation de nappe de type déterministe, trois lois physiques

élémentaires sont strictement nécessaires : loi de conservation de masse, loi de
Darcy et les équations d’état.

Dans notre cas nous avons utilisé un modèle basé sur les différences finis, il
s’agit du Modflow.

- -
 Figure 44. Résolution par modèle des problèmes hydrogéologiques

2. Formulations mathématiques

2.1. Loi de Darcy

Pour les fluides incompressibles, l’équation de mouvement se simplifie sous la

forme tensorielle suivante :

U= -[Link] h

h : charge ou hauteur piézométriques,

K : tenseur de perméabilité,

U : Vitesse moyenne de DARCY.

Si l’on réduit le tenseur de perméabilité à ses composantes diagonales, en

choisissant les axes des coordonnés colinéaires aux directions principales du
tenseur K : Kxx, Kyy, Kzz.

- -
La relation se décompose de la manière suivante :
Ux = - kxx .  h /  x
Uy = - kyy .  h /  y
Uz = - kzz . h /  z ,

Avec :

Ux, Uy et Uz, les vitesses moyennes de filtration de Darcy.

2.2. Equation de conservation de la masse

L’équation de continuité exprime le principe général de la conservation de la

matière au sein d’un volume fermé, c’est à dire la variation des flux entrant et
sortant au sein d’un volume expérimental représentatif (VER) considéré. La
déformation du milieu étant négligée, l’équation s’écrit :

div (U)+  / t ()+q = 0

 : Masse volumique de l’eau,

U : Vitesse de filtration de DARCY,

t : temps

 : Porosité du milieu correspondant au rapport du volume des vides dans le

VER sur le volume total,

q : Débit volumique d’eau prélevé (ou apportée) par unité VER en chaque
point.

Equations d’état
Elles sont indispensables à l’application des deux lois précédentes car elles traduisent les
variations des paramètres mécaniques de l’eau en mouvement et de la matrice rocheuse
(masse spécifique, viscosité …) en fonction de la pression (Bonnet, 1982), elles s’écrivent
de la manière suivante :

d/ = .dp (eau)

dv / v = - .d = . Dp (matrice)

 : Masse volumique de l’eau,

- -
 : Coefficient de compressibilité de la matrice poreuse,

 : Coefficient de compressibilité de l’eau,

p: Pression interstitielle de l’eau,

 : Contrainte effective au sein du volume élémentaire représentatif (VER) de

volume v.

2.4. Equation de diffusivité

Les trois équations précédentes élémentaires de la circulation de l’eau en milieu

poreux (nappe libre) se combinent pour donner l’équation de diffusivité
suivante:

div (K grad h) = Ss. ( h /  t) + q

S: coefficient d’emmagasinement spécifique du milieu poreux sur tout le VER,

égale à .g (+w. )

En faisant une intégration et une linéarisation, nous obtenons l’équation de

diffusivité en nappe libre suivante :

2.5. Résolution par la méthode les différences finies

Cette méthode consiste à discrétiser le domaine à modéliser en mailles carrées

ou rectangulaires et à rechercher à définir sur chaque maille une fonction
constante (Hi) choisie par l’approximation de la charge hydraulique. Les calculs
se feront en admettant que la valeur de la charge hydraulique soit attribuée au
centre de la maille. On y procédera de la même façon pour les différents
paramètres de l’équation, en définissant sur chaque maille :

- une transmissivité Ti (ou la perméabilité Ki);

- un coefficient d’emmagasinement Si;

- un débit total algébrique Qi.

- -
Les fonctions d’approximation auront à satisfaire uniquement les données
locales (chaque maille), on obtient alors un système d’équations linéaires
définissant les valeurs de Hi au centre des mailles.

Sachant que les calculs à l’intérieur d’une maille se font en fonction des nœuds
environnants, on isole une maille du domaine (fig.45) avec ses quatre mailles
voisines, que l’on désigne par N (nord ), S (sud ), E ( est ) et W( ouest ) :

* Le principe de continuité qui implique la conservation du débit d’eau entrant

algébriquement par les limites de la maille se traduit par :
Qn + Qs + Qe + Qw = Qi + Qem

Avec Qemi : débit emmagasiné dans la maille i.

Figure 45. bilan des flux sur une maille (i).

* La loi de Darcy, permet de définir chaque composante du débit entrant en

fonction de la transmissivité et du gradient hydraulique. Elle s’écrit :
Qn = Tn . a .(Hn –Hi) / a = Tn. .(Hn –Hi)

Tn : Transmissivité de l’aquifère entre la maille I et sa voisine dans la direction

Nord. Hi et Hn : Approximations de la charge respectivement sur les mailles I et
N.

* L’équation d’état fournit l’expression du débit emmagasiné :

Qemi = a . Si .( dHi/dt) dt

Si : coefficient d’emmagasinement sur la maille i.

- -
Tous ces calculs permettent d’obtenir :

Tn ( Hn – Hi ) + Ts ( Hs – Hi ) + Te ( He – Hi ) + Tw ( Hw – Hi ) = Qi + a Si .dHi/dt

Si l’on applique la même équation à chaque maille du modèle, on obtiendra un

système différentiel linéaire du premier ordre à « n » équations, à « n » fonctions
inconnues du temps (Hi), que nous pouvons résoudre. Pour simplifier la
présentation, on adoptera la notation matricielle suivante :

[T]. H = Q + a . [S] .dh/dt

En définissant les vecteurs : H et Q, et les matrices T et S.

3. Conditions initiales et conditions aux limites

3.1. Conditions initiales
Elles sont définies par les valeurs de la charge « h » connue en tout point du
domaine modélisé à un temps de référence t = 0.
3.2. Conditions aux limites
Les conditions aux limites concernent les règles d’échange des flux d’eau ou de
matière entre le domaine modélisé et le milieu extérieur. La détermination de ces
conditions est nécessaire pour résoudre l’équation de la diffusivité. On connaît
quatre types de conditions aux limites.
a) Limite à potentiel imposé (condition de Dirichlet)
Tout au long de cette limite, la charge hydraulique est indépendante des
conditions de circulation de la nappe. C’est le cas du contact d’une nappe avec
un plan d’eau libre (mer, lac, oued)
b) Limites à flux impose (condition de Newman)
Le débit par unité de surface est constant ou nul le long d’une limite à flux
imposé. On distingue les limites à flux nuls (contact de l’aquifère avec une
formation imperméable) et les limites à flux imposé non nul qui sont
matérialisées par les limites du domaine ou par les ouvrages de captage.
c) Condition de Fourrier
Définit le cas d’un oued dont le fond est colmaté et qui est en communication
avec une nappe libre (drainage ou alimentation).

d) Conditions de surface libre

Deux conditions définissent une surface libre :

- -
- Condition de potentiel : c’est la pression atmosphérique qui s’exerce sur tous
points M de la surface libre. Exprimée en charge, on écrit : h = z
- Condition de flux : Si la nappe n’est pas alimentée par sa surface, le flux est
nul et il s’écrit :
dh / dn = 0

Si un flux traverse la surface libre de la nappe (alimentation ou évaporation). On a:

dh / dn = a

n : normale orientée vers l’extérieur.

a : flux entrant ou sortant (positif ou négatif).

Les étapes de l’élaboration du modèle

4.1. Introduction des données
Après avoir défini les conditions initiales et les conditions aux limites et procédé
à la discrétisation du domaine à modéliser en un certain nombre de mailles, nous
introduisons pour chaque maille les données physiques du système à savoir:

- La côte piézométrique de la carte qui servira de référence,

- La géométrie de l’aquifère (la côte du toit du substratum),

- Les paramètres hydrodynamiques T, K, S.

-les paramètres climatiques (la recharge efficace)

4.2. Calage du modèle

Le but de cette phase quantitative est de reconstituer la variation des paramètres
de l’aquifère en utilisant le modèle pour les affiner c’est à dire en minimisant les
écarts entre les réponses du modèle (calculées) et les réponses du système
(observées). Ce calage est important car il permettra de réaliser des simulations
sur le domaine étudié.

- -
4.3. Exploitation du modèle

C’est l’étape ultime où le modèle révèle son utilité, son intérêt principal est
d’établir des prévisions des réactions hydrodynamiques futures du système
aquifère à des scénarios d’aménagement éventuels (exploitation, sécheresse,
recharge …), ou faire des simulations de la propagation des polluants vers la
nappe.

5. Application du modèle en 3D au système aquifère

5.1. Introduction

Le modèle tridimensionnel MODFLOW (Mc Donald and Harbaugh, 1988) a été

utilisé pour simuler les écoulements souterrains et estimer la recharge de la
nappe. C’est un modèle à bases physiques, déterministe, capable de représenter
des écoulements laminaires monophasiques tridimensionnels dans des systèmes
multicouches. Il résout l’équation de diffusivité aux dérivées partielles de
l’écoulement des eaux souterraines en milieu poreux (combinaison de la loi de
Darcy et de l’équation de continuité) par la méthode des différences finies.

Le modèle permet d’estimer les ressources en eau à l’état actuel et de prévoir

leur changement dans la future par proposition des scénarios, qui peuvent être
choisis comme une stratégie de gestion de la ressource en eau sans trop de risque
de défaillance.

5.2. Objectifs
L’objectif est de bâtir un modèle numérique d’écoulement afin de tester le
modèle conceptuel. Le concept prévisionniste permettra une utilisation durable
de la ressource suite à une éventuelle surexploitation. L’objectif principal de
l’étude est la simulation des écoulements (simulation des directions des
écoulements des eaux souterraines et la distribution des charges).
Il s’agit de ce fait de :
- estimer la recharge (les apports à partir des flux latéraux),
- prévoir les rabattements futurs et les conséquences des surexploitations,
- implanter au mieux les ouvrages à venir,
- gérer rationnellement les débits afin d’optimiser l’exploitation de la nappe
(réaliser de nouveaux scénarios de prévision).

- -
L’hypothèse d’un régime permanent implique des entrées et des sorties égales et
des variations néglige de l’emmagasinement.
Les aquifères superficiels conditionnés par les infiltrations et l’évaporation sont,
suivant les saisons, en régime transitoire traduis par des fluctuations du niveau
piézométrique.
En revanche, pour les nappes profondes dont l’alimentation est assurée à partir
des entrées latérales, les conditions de régimes transitoires sont établies quand
les paramètres hydrauliques (perméabilité) sont assez faibles ou en situation
d’exploitation intense (rapport flux/stock).

5.3. Discrétisation du domaine

Le nombre et la taille des mailles dépendent de la précision attendue et de la

nature des données sources. L’aquifère a été discrétisé en 60 lignes et 100
colonnes. Les mailles sont carrées et ont pour dimension 500 mètres de côté.
Le modèle est ainsi constitué de 4000 mailles actives, totalisant une superficie
de 1000 km (fig.46).

5.4. Données initiales

Elles concernent :
L’Etat piézométriques de référence qui a permis de caler le modèle en régime
permanent. Il s’agit de la piézométrie du mois de Mai 7 ([Link].
IV). La carte en courbes en isopièzes obtenue, semble en effet la mieux
convenir à la reproduction d’un état stable.
Les perméabilités reparties dans la zone d’étude (fig. 3. Chapitre. IV),
La moyenne annuelle des précipitations enregistrée à la station de Hammam
Grouz, sur une période de seize ans (1988/89 à 2003/04) est de l’ordre de 37
mm,
La part de l’infiltration efficace, évaluée à partir de la méthode de
Thornthwaite est estimée de l’ordre de 7 mm an,
Les conditions aux limites données par la carte piézométriques, généralement
les écoulements des eaux souterraines d’un système aquifère sont régis par leurs
conditions aux limites.
Des flux entrant au niveau des mailles de la bordure Nord-Ouest , sud-ouest ,
nord-est et sud-est, traduisant l’apport latéral à partir de ces zones.

- -
 Figure Discrétisation de la zone d’étude

Les mailles auxquelles il n’a pas été affecté de valeurs de débit ou de charge sont
considérées comme étanches (partie nord et sud).
Dans le cas des fortes crues, la nappe alluviale possède un flux à potentiel
imposée en liaison avec Oued Rhumel dans sa partie avale. Dans la partie
centrale, l’affleurement des calcaires de l’Eocène plus au moins fissurée en
liaison avec la nappe superficielle favorise à l’infiltration directe des eaux de
l’oued.
5.5. Calage du modèle en régime permanent
Le calage du modèle doit permettre une synthèse de toutes les données du
domaine discrétisé et la détermination des données manquantes. La technique
consiste à ajuster l’état calculé à l’état réel en changeant certains paramètres
hydrodynamiques, tel que la perméabilité.
- Reconstitution de la piézométrie
L’examen de la figure (fig. 7) montre une bonne similitude entre la
piézométrie mesurée et celle calculée par le modèle avec des écarts inférieurs à
1 m. A cet effet le modèle reste valable et le décalage entre les niveaux
piézométriques calculés et mesurés dans certaines zones sont probablement dues
à l’imprécision des mesures, et aux conditions hydrogéologiques spécifiques.

- -
 Figure Reconstitution de l’état piézométrique en régime permanent

5.5.1. Optimisation du champ de la perméabilité

Les différents itérations du modèle ont permis d’ajuster la conductivité
hydraulique de l’aquifère dont les valeurs de la perméabilité s’échelonnent
entre 10- m/s pour les secteurs à très faible perméabilité et 10- m/s pour les
secteurs à forte perméabilité caractérisant la partie sud-ouest du terrain d’étude
(fig.48).

- -
 Figure 48. Carte de la perméabilité de la nappe issue du modèle

5.5.2. Bilan hydrogéologique issu du calage

La phase du calage en régime permanent a permis d’établir le bilan suivant
(tab.13).
- Apport par les calcaires du crétacé et de l’éocène: .3 m3/j
- Apport par les précipitations (la recharge) : 20 559 m3/j.
- Débit de pompage : 1250 m3.

- -
Tableau 13. Bilan hydrique calculé en régime permanant
Entrées Sorties Variations de réserve Stockage
)
( m /j ) ( m /j ) (m /j (m /j )
-

Les résultats du bilan hydrique obtenus par le modelé, ont permis de révéler que le
système d’écoulement de la nappe en régime permanant est dominé en entrée par la
recharge latéral , et en sortie par le pompage intensif, dont la variation de réserve est
estimée de 3.65 m3. Le bilan hydrique fait ressortir une recharge latérale totale de
5 7.55 m3, dont % représente l’apport par les précipitations (la recharge
efficace).Fig.49.

Figure 49. Bilan hydrique de la nappe alluviale de Tadjenanet – Chelghoum Laid

issu du modèle en régime permanant

6. Exploitation du modèle et simulation

Afin de caractériser l’état piézométrique en régime transitoire nous avons

procéder premièrement au calage du modèle. Il s’agit de la piézométrie du mois
de Mai 2013 (fig.38 .Chapitre. IV). La technique consiste à ajuster l’état calculé
à l’état réel.
L’examen de la figure (fig.5 ) montre une bonne similitude entre la piézométrie
mesurée et celle calculée par le modèle avec des écarts inférieurs à 2 m. A cet
effet le modèle reste valable et le décalage entre les niveaux piézométriques
calculés et mesurés dans certaines zones sont probablement dues à l’imprécision
des mesures, et aux conditions hydrogéologiques spécifiques.

- -
 Figure Reconstitution de l’état piézométrique en régime transitoire Mai

Apres la calibration du modèle, et à fin de l’exploiter en régime transitoire pour

définir l’état future de la nappe et ses variation hydrodynamiques dans le temps,
la simulation en ce régime transitoire a été effectuée sur une période de 20 ans.

Trois scénarios ont été proposés :

er
- scénario : ce scénario sert à tester les capacités d’emmagasinement sur une période
de 20 ans, qui correspondent à une exploitation à long terme. Les débits sortants pris en
compte sont ceux de l’année de référence et sont maintenus constants. En gardant le
nombre d’ouvrage de captage et leurs débits d’exploitation fixe. Dans ce scénario, deux
cas ont été pris en compte : le premier cas sans recharge et le deuxième avec recharge. Le
modèle ne semble pas être influencé par ce paramètre. Ce qui explique l’importance de
l’apport latérale.
L’examen de la carte piézométrique obtenue avec ce scénario (fig.51), montre
que la nappe garde sensiblement la même allure piézométrique et les mêmes
directions d’écoulements, avec des rabattements estimés de à m à l’ouest et
de 5 à 6 m aux parties centre et Est du terrain d’étude . (fig.5 ).

- -
 Figure Simulation d’une exploitation à long terme

Figure .52 Evolution temporelle de la piézométrie (scénario 1)

La reconstitution du bilan global de l'aquifère en régime transitoire par évaluation de ses

différentes composantes est représenté dans le tableau suivant (tab.14). Il dépend de la
recharge et du régime de l'exploitation, ces deux paramètres présentent une forte
irrégularité.

- -
 Tableau.14 Bilan hydrique calculé en régime transitoire
Entrées Sorties Variations de réserve Stockage
( m /j ) ( m /j ) (m /j ) (m /j )
-

L’examen de la figure (fig.53) nous a permis de constater que le bilan hydrique de la nappe
après une période d’exploitation de ans est influencé par la recharge latérale et du régime
de l'exploitation, ces deux paramètres présentent une forte irrégularité.

Figure.53 Bilan hydrique de la nappe alluviale de Tadjenanet – Chelghoum Laid

issu du modèle en régime transitoire
ème
- scénario : consiste à maintenir le même nombre de puits de pompage et
de multiplier par dix leurs débits de pompage.
La simulation étalée sur 20 ans, nous a permis de constater que cette situation
provoque des grandeurs de rabattement d’ordre de à 3m dans la partie ouest et
de 5 à 9m sur toute la partie centre et Est de la nappe. Fig. 54.

- -
 Figure.54 Carte piézométrique résultant du deuxième scénario.

ème
- scénario : implantation de nouveaux forages sur les limites
La simulation étalée sur 20 ans, nous a permis de constater que le modèle ne semble pas
être influencé par cette proposition et la nappe garde la même allure du premier scénario.
Fig .55.

Figure 55. Carte piézométrique résultant du troisième scénario

- -
7. Conclusion

La caractérisation hydrodynamique de la nappe alluviale de Tadjenanet -

Chelghoum Laid en régime permanent simulé par MODFLOW, a permis
d’identifier les principaux paramètres du bilan global de la nappe. Les
simulations en régime transitoire effectuées sur une période de 20 ans, ont
permis de constater que la nappe subirai un rabattements, atteignant les 6 m,
étalé sur la partie centre et Est de la nappe avec des faibles grandeurs d’ordre à
2 m enregistrés à la partie ouest . A cet effet, nous proposons, l’implantation
d’éventuel ouvrage de captage aux extrémités des massifs carbonatés près des
limites de la région d’étude à fin d’éviter une telle situation entrainant la
pollution de la nappe surtout dans sa partie centrale, sachant que cette partie est
connue par sa vocation agricole.

- -
CHAPITRE V

HYDROCHIMIE
 HYDROCHIMIE

1. Introduction

Tout travail hydrochimique consiste à recueillir des échantillons d’eau, tout en respectant
les techniques de prélèvement prescrites. De procéder à un certain nombre d’analyse pour
avoir une idée sur l’identité physico-chimique de l‘eau et son mode de répartition spatiale
dans le milieu environnant. La composition chimique des eaux souterraines dépend en
grande partie du type de roche encaissante, de l'évolution hydrodynamique des
écoulements de la nappe, de la profondeur de cette dernière et de la zone climatique. Elle
donne de précieux renseignements sur les conditions de gisement et le mouvement des
eaux souterraines.
Ce travail permet d’aborder :
- La nature et la qualité physico- chimique des eaux de la région
- De montrer les conséquences de la qualité des eaux souterraines sur le
risque de dégradation des terres agricoles.
- Les problèmes de potabilité de l’eau dans la région.
-L’impact des facteurs anthropiques sur la qualité physicochimique des eaux
de surface et eaux souterraine.

Matériels et moyens d’analyse

La campagne de prélèvements d’échantillons a été effectuée durant la fin de la période de

hautes eaux (Mai 2007). Ces prélèvements ont été effectué selon une répartition et une
fréquence qui dépend de la disponibilité des ouvrages, du type des formations
géologiques et de la profondeur des eaux. Trente et trois (33) échantillons ont été pris au
niveau des points d’eau (puits et forage). Fig. 5 .
Les analyses chimiques ont été effectuées au laboratoire de la station d’épuration de
Batna et dans un laboratoire privé d’analyse de sol et des eaux d’irrigation. Les éléments
chimiques qui ont été dosés sont : Ca , Mg ,Na+ ;K+,HCO3-,SO -,Cl-,NO3-,NO -

,NH +,PO 3-. Les méthodes de dosage utilisées sont : la colorimétrie, la volumétrie et la
spéctrophotométrie (RODIER, 1996). La température (T), la conductivité électrique (CE)
et le pH ont mesurés sur terrain à l'aide d'un appareil portatif multi paramètres.

- -
 335 Dj. G
rouz
ss
De 27
d 31
.e
330 Dj 32 28

9 Vers Constantine
33
36 76
325 44 35 34 Chelghoum Laid
37
4240
43 39 23 RN.5
38 22 11
46 2
Tadjenanet 21 18
49 48
320 Vers Sétif
57 50
53
62
63
315

67

310 ya
ou
Tik
.
Dj
305
780 785 790 795 800 805 810 815 820 825 830

Légende
. Points d'eau Profil d'évolution
0
Echèlle

2 4 km
----- Oued
Ville
Limite du terrain d'étude

Route

Figure Carte d’inventaire des points d’eau ayant fait l’objet d’une analyse physico-
chimique

Caractéristiques physico-chimiques

Caractérisation des eaux de l’aquifère

La composition chimique des eaux souterraines dépend en grande partie du type de roche
encaissante, du temps de séjour, de l’évolution hydrodynamique, de la profondeur du
niveau statique et du climat. La qualité de l’eau peut également être influencée par
l’activité de l’homme dans le domaine de l’urbanisation, de l’industrialisation et surtout
les différentes activités agricoles.
Un autre facteur qui a un effet considérable sur la qualité des eaux, surtout dans les
régions arides et semi-arides, qui est l’augmentation de la salinité des eaux souterraines
surtout pour les aquifères superficiels, résultant du phénomène de l’évaporation (Drever,
7). Les résultats des analyses sont regroupés dans le tableau 15. La dispersion autour
de la moyenne (m) de ces éléments s’est faite à l’aide de l’écart type (σ) et du coefficient
de variation (CV= σ/m), on estime que la série est homogène lorsque CV <50%.

- -
 Tableau 15. Caractéristiques statistiques des paramètres physico-chimiques des
eaux la nappe alluviale de Tadjnanet- Chelghoum Laid. Mai2007
Paramètres T pH CE Ca Mg Na+ K+ HCO -
SO -
Cl- NO -
NO -
NH +
PO -

Statistiques (°c) (µs/cm) mg/l mg/l mg/l mg/l mg/l mg/l mg/l mg/l mg/l mg/l mg/l

Max . 8. 3 3 3 . 7 . 3 .3 5 . 33 . 5 8 . . 5 .5 .
Min 5. 7. 7 5 5. 8 3 . .3 . 88.75 3 . . .
Moy 7.3 7. 5 8.3 7 . 7 7 . 5 .7 5 . 3 7 . 3 .58 . . . 5 . 5
Ecart type . .3 3 8. 3 .37 8. .5 7. 8. 7 . 8 . 5 . . 8 .
CV (%) . . 7.3 . 55.7 3 .3 37. 37.88 . 7 . 7 .87 5 7 3.
Norme .5-8.5 5 5 5 5 . .
OMS

Selon les résultats d’analyses la chimie des eaux, la nappe alluviale est caractérisée par
une grande variation des concentrations des éléments chimiques : Na+ et Cl- (26-339
3--
mg/l), SO -, Ca , Mg , et HCO3- (29-1100 mg/l), NO3- (13-486 mg/l), NH +, PO ,
NO2- (0.01-2.2 mg/l). On remarque que l'écart par rapport à la moyenne est considérable
pour certains éléments tels que le magnésium, les sulfates et les nutriments. Fig. 57.
La grande variabilité des composés azotés telle que les nitrates est due probablement à
l’utilisation des engrais chimiques et organiques dans l’agriculture et à la décomposition de la
matière organique (Kehew, 2000). (Groupe 1). En ce qui concerne les sulfates qui
proviennent de la dissolution des formations gypseuses, les concentrations élevées sont
limitées dans l'espace et ne concernent pas toute la plaine. (Groupe 2).
La déviation du reste des paramètres par rapport à la moyenne est de l'ordre de 40%. Ceci
peut être expliqué par une répartition assez régulière de la concentration de ces paramètres
dans l'espace, excepté pour le magnésium où les fortes concentrations sont liées aux sulfates.
(Groupe 3). Concernant les paramètres physiques, ils ont une répartition identique sur
l’ensemble du terrain. (Groupe ).

Figure 57. Le coefficient de variation des variables physico-chimiques mesurées des

eaux souterraines
- -
3.2. Répartition spatiale des paramètres physico-chimiques
La répartition spatiale des paramètres physico-chimiques, permet de mettre en évidence
l’évolution de chaque élément et d’estimer l’origine de minéralisation.
3.2.1. Paramètres physiques
Potentiel d’hydrogène (pH)

La détermination du pH constitue une mesure de la concentration des ions H+ dans l'eau.

Ce paramètre caractérise l’équilibre physico-chimique et dépend des facteurs multiples
dont l’origine de l’eau. Le pH détermine l’acidité, l’alcalinité et la neutralité des
solutions, il est lié à la nature du terrain. En général, la modification du pH dans l’eau est
liée à la perte du gaz carbonique. Il convient aussi de signaler que la valeur du pH est
fonction de la matière organique thermostable dissoute ou non dans l’eau. Généralement
les valeurs de pH, diminuent pendant les périodes pluvieuses, cette diminution peut être
attribuée au gaz carbonique qui charge les eaux d’infiltration à travers les différentes
couches du sol. Les minéraux rencontrés seront dissous sous l’influence du H CO3
(Drever, 1997).
A partir du pH, on peut classifier les eaux de la manière suivante (tab.16)

Tableau Estimation de la qualité de l’eau selon le pH

pH Qualité de l’eau

pH < 5 Acidité forte

pH = 7 Neutralité (eau pure)

7< pH < 8 Neutralité rapprochée (eau de surface)

5.5< pH >8 Majorité des eaux souterraines

pH > 8 Alcalinité, évaporation intense

D’après la carte de pH (fig.58) les eaux souterraines de la région, ce paramètre varie

entre 7.2 et 8.63. La totalité des points échantillonnés sont conformes aux normes à
l’exception du point d’eau n° situé à l’extrême nord où le pH est de l’ordre de 8. 3.
L’évolution spatiale du pH montre une légère variation de l’Ouest vers l’est.

- -
 335
Dj. G
rouz
ss
De
d
.e
330 Dj
8.5

325

320 8

315

7.5
a
310 uy
iko
. T
Dj
305 7.1
780 785 790 795 800 805 810 815 820 825 830

Echèlle
0 2 4 km

Figure 58. Carte du pH

La conductivité électrique

La conductivité électrique mesure la capacité de l’eau à conduire le courant électrique. Sa

variation permet de suivre l’évolution de la salinité globale de la chimie de l’eau. La
classification des eaux en fonction de la conductivité électrique se présente de la manière
suivante :

- Conductivité égale à 0.05 (μS/cm): eau déminéralisée

- Conductivité de 10 à 80 (μS/cm): eau de pluie
- Conductivité de 80 à 100 (μS /cm) : eau peu minéralisée
- Conductivité de 300 à 500 (μS /cm) : eau moyennement minéralisée
- Conductivité de 1000 à 3000 (μS /cm) : eau saline
- Conductivité supérieure à 3000 (μS /cm) : eau de mer.
La répartition de la conductivité électrique montre une augmentation dans le sens
d’écoulement de l’Ouest vers l’Est. Le maximum est observé au centre avec 3610 μS/cm.
Fig. 59.

- -
 ( s /cm)
335
Dj. G
rouz
ss
De 3600
d
.e
330 Dj
3400

325 3000

2600
320
2200

315 1800

1400
a
310 uy
ko
. Ti 1000
Dj
305 700
780 785 790 795 800 805 810 815 820 825 830

Echèlle
0 2 4 km

Figure 59. Carte de la conductivité électrique

3.2.2. Paramètres chimiques

Le calcium

Le calcium, élément principal de la dureté de l’eau, est un métal alcalino-terreux extrêmement

répandu dans la nature et en particulier dans les roches calcaires sous forme de carbonates.
Dans les roches sédimentaires, il se présente sous forme de carbonates (calcite, dolomite
et aragonite). (Derver 1997). Le calcium fait partie de quelques minéraux évaporitiques à
savoir le gypse et l’anhydrite et aussi dans le minéraux argileux sous forme de Ca-
montmorillonite (Hem, 1985).
L’examen de la carte de la répartition du calcium (fig. ), montre que la teneur la plus élevée
ont est enregistrée au centre de la plaine, tandis que les faibles teneurs sont observées au
sud-est de la zone d’étude, notamment prés de Djebel Tnoutit et au nord de Tadjnanet. Ceci
est en étroite liaison avec l’évolution de la conductivité électrique.

- -
 335
Dj. G
rouz
ss (mg/l)
De
d
.e 300
330 Dj

260
325

220

320
180

315
140

a
310 uy 100
iko
. T
Dj
60
305
780 785 790 795 800 805 810 815 820 825 830

Echèlle
0 2 4 km
Figure 60. Carte de calcium

Le magnésium

Le magnésium est un élément très répandu dans la nature ( . % de l’écorce terrestre),

élément aussi principal de la dureté de l’eau. Dans les roches sédimentaires ses origines
sont comparables à ceux du calcium, car il provient de la dissolution des formations
carbonatées à fortes teneurs en magnésium comme la magnésite et la dolomite (Drever
7).
L’étude de la carte de la répartition spatiale du magnésium (fig.61) montre une répartition
similaire à celle du calcium. Les concentrations varient entre 32.49 mg/l et 172.44 mg/l.
Les teneurs les plus élevées sont observées au niveau des points d’eau localisés le long de
Oued Rhumel, au centre et au sud de la région d’étude. Ces fortes concentrations peuvent
être attribuées au contact direct de ces eaux avec les argiles gypseuses et les formations
carbonatés et dolomitiques limitrophes.

- -
 335
Dj. G
rouz
ss (mg/l)
De
d
.e
330 Dj
160

325
120

320

80

315

a 40
310 uy
iko
. T
Dj
305 30
780 785 790 795 800 805 810 815 820 825 830

Echèlle
0 2 4 km

Figure 61. Carte du magnésium

Le sodium

Cet élément, très abondant (environ 26g/kg de la croûte terrestre soit environ 2.8%) ne se
rencontre pas naturellement à l’état natif mais toujours associé à d’autres éléments
chimiques (Chlorures, sulfates….). Le principal minéral contenant le sodium est le sel
gemme (chlorures de sodium). Des teneurs anormalement élevées peuvent provenir du
lessivage ou de la percolation de terrains salés ainsi que de l’effet de la salinité marine.
Au niveau de la région d’étude les teneurs de sodium varient entre un minimum de .35
mg/l dans la partie sud-ouest du terrain et un maximum de 143.8 mg/l dans la partie Est.
La carte de la variation spatiale de cet élément (fig.62), montre que les plus grandes
concentrations se localisent dans le secteur est vers le centre. Ces teneurs sont attribuées
probablement des mécanismes chimiques différents liés à l’échange de base qui est
influencés par le pompage et le sens d’écoulement. Les parties ouest et sud-ouest se
caractérisent par de faibles teneurs avec un minimum d’ordre de mg l au niveau du
point d’eau n° 3.

- -
 335
Dj. G
rouz
ss (mg/l)
De
d 140
.e
330 Dj

120
325

100

320
80

315
60

a
310 uy 40
iko
.T
Dj
305 20
780 785 790 795 800 805 810 815 820 825 830

Echèlle
0 2 4 km

Figure 62. Carte de sodium

Les bicarbonates

Ils proviennent de la dissolution des formations carbonatées (calcaire, dolomie) par des
eaux chargées en gaz carbonique (CO ) selon l'équation suivante:
CaCO + H O + CO → HCO - + Ca
Les concentrations en bicarbonates dans les différents échantillons analysés, varient entre
un minimum de 122 mg/l au niveau du puits n°63 et un maximum de 512.4 mg/l au
niveau du puits n° . L’évolution spatiale des bicarbonates (fig. 3), montre une
augmentation des concentrations de cet élément dans la partie est qui coïncide avec la
zone de convergence des eaux souterraines de l’aquifère. Ces fortes teneurs sont dues
probablement à la dissolution des calcaires lacustres de Moi-Pliocène au Nord et au Nord-
Ouest de Tadjnanet, les formations carbonatées limitrophes d’âge Eocène, constituant le
massif de Djebel Tikouya au sud-est et Draa Mestaoua ainsi que les formations d’âge
Crétacé qui caractérisent le massif de Djebel Grouz.

- -
 335
Dj. G (mg/l)
r ouz
ss 520
De
d
.e 480
330 Dj
440

325 400

360

320 320

280

315 240

200

310 ya 160
ou
ik
.T
Dj 120

305 118
780 785 790 795 800 805 810 815 820 825 830

Echèlle
0 2 4 km
Figure 63. Carte des bicarbonates

Les chlorures
Très répandus dans la nature, généralement sous forme de sels de sodium (NaCl), de
potassium (KCl) et de calcium (CaCl ), les ions chlorures (Cl- ), la grande partie se trouve
dans les océans, constituent environ . 5 % de la lithosphère. De façon générale, l’ion
chlorure est présent dans toutes les eaux à des concentrations différentes. Dans les eaux
souterraines, les teneurs en ions chlorures peuvent atteindre quelques grammes par litre
au contact de certaines formations géologiques.
Dans la région, la teneur varie entre un minimum de 88.75 mg/l enregistré au niveau du
puits n°3 et un maximum de 33 . 5 mg l enregistré au niveau du puits n° 3. L’examen
de la carte des chlorures (fig. 64), montre des concentrations élevées dans la partie
centrale et sud-est du terrain identique à elle de sodium. Plus au Nord et au Nord-Ouest
l’eau est moins chargée. La répartition des teneurs des chlorures se corrèle bien avec
celles de sodium, ce qui confirme l’origine salifère de ces deux éléments.

- -
 335
Dj. G
rouz
ss (mg/l)
De
d
.e 320
330 Dj

280
325

240

320
200

315
160

a
310 uy 120
iko
. T
Dj
80
305
780 785 790 795 800 805 810 815 820 825 830

Echèlle
0 2 4 km

Figure 64. Carte des chlorures

Sulfates

Composés naturels des eaux, ils sont liés aux cations majeurs : calcium, potassium et
sodium. Ils proviennent essentiellement de la dissolution du gypse (CaSO4, 2H2O). Une
autre source de sulfates peut être attribuée à l’infiltration de sulfure utilisé comme
fertilisant ou comme pesticide .le gypse se dissout très rapidement en contact avec l’eau,
la réaction de dissolution de gypse est représentée par l’équation suivante (Freeze et
Cheery, 1979):
CaSO , 2H O Ca+ + + SO -
+ 2H O
L’examen de la carte des sulfates (fig. 5), nous a permis de constater que les eaux sont
très riches en sulfates, les fortes concentrations se localisent à l’Est, et au centre avec des
concentrations qui dépassent la norme, (1100mg/l) enregistrée au niveau du puits n°9. La
plus faible teneur est d’ordre de mg l enregistrée dans le puits n° 3. Les fortes teneurs
sont liées à l’effet de la dissolution des formations évaporitiques.

- -
 335
Dj. G
rouz (mg/l)
ss
De 1100
d
.e
330 Dj 1000

900
325
800

700
320 600

500

315 400

300
a
310 uy 200
iko
. T
Dj 100

305 0
780 785 790 795 800 805 810 815 820 825 830

Echèlle
0 2 4 km

Figure 65. Carte des Sulfates

[Link]. Les nitrates

Les nitrates sont abondamment répandus dans le sol, dans la plupart des eaux et ils sont
nécessaires à la synthèse des végétaux. Soluble dans l’eau, ils se retrouvent naturellement en
faible concentration dans les eaux souterraines et les eaux de surface et sont accentuées par
l'activité humaine qui accélère le processus d'enrichissement en éléments nutritifs par l'apport
d'engrais sur les sols cultivés. L’azote dans les sols provient principalement de la fixation de
l’azote atmosphérique, du dépôt de l’azote contenu dans les précipitations, des amendements
organiques et des apports d’engrais chimique.
Les apports d’azote sous forme organique dans les sols se font par amendements anthropiques
(fumier, etc…), restitution au sol des résidus de culture et ou restitution des déjections
animales. La constitution du stock d’azote minéral dans le sol est issue des apports d’engrais
de synthèse, sous forme d’ammonium, de nitrates et d’urée.
Le cycle de l’azote dans le sol est très complexe, il est représenté sous forme simplifiée dans
la figure 66.

- -
 Figure 66. Cycle de l’azote dans le sol Hill,

L’examen de la carte (fig. 7) montre que les zones les plus exposées à la pollution par les
nitrates se situent dans la partie centrale de la plaine, dû principalement à la nature
alluvionnaire des terrains qui favorise l’immigration des éléments azotés vers la zone
saturée sous l’effet de l’intensification des activités agricoles et l’excès de pompage. La
forte teneur est enregistrée au niveau du puits n°23 avec une teneur de 486.41mg/l. Les
faibles teneures sont enregistrées dans la partie sud-ouest du terrain vu la profondeur
élevée de la nappe ainsi que la nature argileuse du toit de la nappe, ce qui protège la
nappe contre les infiltrations des ions nitrates. La faible teneur est enregistrée dans le
puits n° par une concentration d’ordre 3mg l.
335
Dj. G
rouz mg/l
ss
De 440
d
.e
330 Dj 400

360
325
320

280
320 240

200

315 160

120
a
310 uy 80
iko
.T
Dj 40

305 0
780 785 790 795 800 805 810 815 820 825 830

Echèlle
0 2 4 km

Figure 67. Carte des nitrates

- -
La variation des teneurs en nitrates dans la région d’étude se répartit en classes (tab.
7).

Tableau Répartition des teneurs en nitrates dans la région d’étude Mai

Classe Intervalle Nombre %

NO >200mg/l
100<NO <200mg/l
50<NO <100mg/l
NO <50mg/l

Ce tableau montre que 78.78% des points d’eau répartit sur la région possèdent des
teneurs supérieur aux normes selon l’OMS fixée à 5 mg l, par contre la reste des points
d’eau présente des teneurs inférieures aux normes.

Evolution spatiale des nitrates

L’évolution spatiale des nitrates selon le long d’un profil Ouest-Est montre que les teneurs
augmentent d’une manière alarmante. Le pic est atteint au puits avec 8 mg l et
mg l à l’extrême Est. Cette situation est en liaison avec la présence d’une immense
surface agricole où l’utilisation irrationnelle des engrais chimiques est excessive. Fig. 8.

Figure Evolution spatiale des nitrates selon l’écoulement

Influence de la profondeur de l’eau sur la teneur en nitrates

L’évolution de la concentration des nitrates en fonction de la profondeur de l’eau (fig.

) montre que les points d’eau dont le niveau piézométrique est plus proche à la surface
du sol, possèdent de fortes teneurs en nitrates. Ceci est dû à la facilité des ions nitrates
d’être transportés par l’eau à travers la zone non saturée. Ceci est surtout marqué pour les

- -
points d’eau se trouvant dans la partie Est et centrale. La figure suivante montre les
différents mécanismes de pollution de l’eau par les nitrates. Fig. 7 .

Figure 69. Relation nitrate – profondeur de l’eau

Figure Mécanisme de la pollution de l’eau par les nitrates Julien,

[Link]. Nitrites
Les nitrites représentent la forme la moins oxygénée et la moins stable de l’azote, elle
représente la forme du passage entre les nitrates et l’ammonium, c’est une forme toxique.
Elle est due à l’effet de l’oxydation de la forme ammonium. Les analyses montrent que
les eaux souterraines de la région renferment des teneurs très faibles (0 à 0.05mg/l).

- -
[Link]. Ammonium
L’ammonium est la forme d’azote la plus toxique. Sa présence dans l’eau est liée soit aux
rejets urbains soit aux rejets industriels. Les concentrations observées sont dans la région,
sont inférieures à la norme et varient entre 0 à 0.56mg/l.
Phosphates
Les phosphates peuvent se trouver sous différentes formes oxydées, méta HPO 3, pyro
H P O7 et ortho H3PO . En milieu aqueux, les formes méta et pyro tendant vers la forme
ortho pour des pH de 5 à 8. Les réactions chimiques entre les différentes formes de
phosphate sont données par :
3-
PO + 3H+  HPO --
+ 2H+  H PO --
+ H+  H3PO
Les concentrations observées au niveau de la région sont dans l’ensemble inférieures à la
norme et varient entre 0 à 2.2 mg/l.
3.2.3. Evolution spatiale des formes azotées
L’examen de la figure d’évolution des différentes forme azotées selon le profil ouest est
de l’écoulement, nous a permis de constater que l’apparition des nitrites et l’ammonium
dans l’eau débute à partir d’un seuil d’une teneur moyenne en nitrates supérieur à
mg l dans la partie centrale du terrain d’étude, les faibles concentrations dans l’eau de ces
élément est dû aux conditions de réduction moins développées de la forme nitraté. Les
phosphates deviennent plus au moins important dans la partie Est, lié à l’utilisation des
engrais chimique de type NPK et aux rejets domestique et industrielle à l’aval de la ville
Chelghoum Laid.(fig.71)

Figure 71. Evolution spatiale des formes azotées et les phosphates selon l’axe
d’écoulement

- -
4. Faciès chimique des eaux

A partir des calculs des quantités en réaction (conversion des concentrations en méq/l) et
la classification des ions (formule ionique) des données hydrochimiques de 33 points
d’eau prélevés dans la plaine de Tadjenanet Chelghoum Laid, nous a permis de
déterminer quatre faciès : classification de Stabler
• Le faciès sulfaté calcique
• Le faciès chloruré calcique.
• Le faciès bicarbonaté calcique.
• Le faciès chloruré magnésien.
• Le faciès sulfaté magnésien
Dont le faciès sulfaté calcique est le plus dominant avec 36 % des échantillons analysés,
30 % le faciès chloruré calcique, Le faciès bicarbonaté calcique, chloruré magnésien et
sulfaté magnésien représentent respectivement 24 % ,6 % et 3 % des échantillons
analysés.
D’après cette répartition préliminaire des faciès chimique des eaux il parait que
l'influence des formations évaporitiques et carbonatées est bien marquée par la dominance
des sulfates, chlorures et les bicarbonates selon le pôle cationique et par calcium et
magnésium selon le pôle anionique.
Pour apprécier le faciès chimique des eaux, plusieurs techniques de traitement basées sur
l’établissement des diagrammes de classification ont été utilisées tels que celui de
Schoeller-Berkallof et Piper, ainsi que des méthodes de traitement statistiques, afin de
mieux cerner les particularités des eaux étudiées.

4.1. Diagramme de SCHOELLER-BERKALLOF

Ce type de représentation, établie par Schoeler et repris par Berkallof, nous portons sur un
semi logarithmique les quantités en réactions des anions et des cations sur l’axe des
ordonnées et sur les abscisses les éléments chimiques à pas régulier, et ceci pour chaque
point d’eau de prélèvement.
Le diagramme tracé pour la période de Mai 2007 (fig.72), confirme principalement un
faciès sulfaté calcique, permettant de dire que la minéralisation des eaux est liée aux ions
-
sulfatés (SO ) et au calcium (Ca ). Secondairement on observe un chloruré calcique.
Généralement les graphiques qui sont marquées par des fortes concentrations et qui ont en
première position sur la figure sont liés à la présence des terrains salifères ou
évaporitiques caractérisent les puits se trouvant à l’est et au centre de la région d’étude.
- -
Les nitrates présentent aussi des teneurs qui dépassent les normes témoignant ainsi
l’apport anthropique vers la nappe.

Figure 72. Diagramme de Schoeller-Berkallof. Mai 2000

- -
4.2. Diagramme de PIPER

Ce diagramme est un autre mode de représentation graphique et de classification d’une

eau analysée. De forme triangulaire, il est subdivisé en trois zones de faciès (fig.73).
Les nuages de points concentrés dans une zone représentent pour les différents
échantillons la combinaison des éléments cationiques et anioniques.
Les résultats d'analyse portés sur le diagramme de Piper confirment l'origine commune de
ces eaux. Elles se regroupent en un seul pôle influencé par la dissolution des formations
évaporitiques. Cette dissolution a générée plusieurs éléments chimiques, notamment les
sulfates, les chlorures, le calcium, le magnésium. Qu’en allant de l’ouest vers l’est les
eaux ont une tendance vers le pole anionique sulfaté. L’ensemble des points d’eau
représentés sur le triangle relatifs aux cations se mettent au milieu (partie mixte) ont une
tendance vers le pole calcique, excepte les points d’eau n°53et 7 situés à l’ouest du
terrain d’étude.

Figure 73. Diagramme de Piper. Mai 2007

Un autre moyen pour mieux apprécier l'origine des faciès caractéristiques des eaux de la
région est l'étude des rapports caractéristiques: SO /HCO3, Ca/Mg et SO /Cl, Ca/Mg,
afin d’avoir une idée sur la répartition et l’évolution des faciès chimiques.

- -
L'examen des deux relations des rapports caractéristique nous a permis de révéler la
dominance des ions sulfatés, chlorures et bicarbonatés avec le pôle cationique calcique,
ce qui explique l’origine de ces eaux est de nature diverse (fig.7 ).

Figure 74. Répartition des faciès chimiques selon les rapports caractéristiques

4.3. Répartition spatiale des faciès chimiques

L'étude de la répartition spatiale des faciès n'est qu'en réalité que l'étude de la répartition
des paramètres chimiques caractérisant ces faciès sur toute la plaine. Les paramètres
dominants des eaux souterraines sont les sulfates, les chlorures et le calcium. D’après
cette carte, nous constatons que les eaux passent d’un faciès chloruré et bicarbonaté à un
faciès sulfaté. Fig. 75. La coexistence de ces faciès est en liaison avec la diversité

- -
géologique des terrains. La présence de formations carbonatées et de formations
salifères triasiques permet à l’eau d’avoir un faciès variable selon la nature du terrain.

335
Dj. G
rouz
ss Bicarbonaté calcique
De
d
.e
330 Dj

Sulfaté magnésien
325

320 Chloruré magnésien

315
Chloruré calcique

a
310 uy
iko
.T Sulfaté calcique
Dj
305
780 785 790 795 800 805 810 815 820 825 830
Axe d'evolution Echèlle
0 2 4 km

Figure 75. Carte des faciès chimiques .Mai 2007

L'étude de l'évolution spatiale des faciès chimiques des eaux de la nappe en liaison avec
la conductivité électrique (fig.7 ) nous a permis de constater que selon l’axe
d’écoulement, le faciès chloruré calcique et bicarbonaté calcique appariaient à l’ouest
sont devenu sulfaté calcique au niveau du centre du terrain d’étude à partir des grandeurs
de conductivité électrique supérieur à 1500 S/m, dont elle atteint des valeurs très
élevées d’ordre de 3 S/m, ce qui est expliqué par l’accélération du phénomène de
dissolution et l’effet d’échange de base entre les alcalins et alcalino- terreux.

- -
 Nord Ouest Sud Est

Figure 76. Répartition spatiale des faciès chimique de la nappe alluviale .Mai 2007

5. Analyse statistique des données physico- chimiques

L’analyse multidimensionnelle est une méthode statistique de description d’un tableau de

données. Elle synthétise et classe les données et révèle les facteurs qui sont à l’origine de
l’évolution des variables .Cette méthode constitue donc un outil de choix en hydrochimie,
elle permet de remplacer une variable par une combinaison linéaire de plusieurs variables
(Plagnes 2000).

- -
5.1. Corrélation entre les éléments
Pour une meilleure caractérisation des paramètres chimiques et leurs évolutions, on a
procédé à l’analyse statistique par corrélation des différents paramètres analysés durant
la période d’investigation (anions majeurs, cations majeurs, pH, conductivité et éléments
azotés).
La matrice de corrélation est une matrice carrée caractérisée par un coefficient de
corrélation qui peut être utilisé pour distinguer la relation entre deux paramètres
(Wuensch, ).
Les différentes relations entre les paramètres chimiques par l’estimation du coefficient de
corrélation peuvent être déterminées sur la base du coefficient critique théorique tiré de la
table de Bravais-Pearson en fonction du nombre de couples mis en jeux. (tab . 8) Sur
cette base le coefficient critique est égal à . .
Une matrice de corrélations a été réalisée pour l’ensemble des prélèvements d’eau
analysés à l’aide du logiciel STATISTICA Version 7. (tab . ).
Tableau. 18 Coefficient de corrélation linéaire selon Bravais- Pearson
v/e v/e
. 877 . . 5 .3887 . 555 .5 85
. . 5 . 8 .3783 . 38 .5 55
.8 5 .8783 . 3 3 .3 87 . 3 .5 3
.7 3 .8 .88 .35 8 . 7 .
. .75 5 .83 .3 33 .38 . 5
. 5 .7 7 .7887 . .3 . 3
.58 . .7 8 . 7 .3 .38
.5 . 3 .7 55 . 573 .3 .3578
.5 . . 85 . 8 . 875 .338
. 73 .575 . 58 . 3 . 73 .3 8
. 7 .55 . 33 . 8 . 5 . 8
. 575 .53 . . 5 . 3 . 737
. .5 3 .5 3 . 8 . 7 . 5 5
. 5 . 73 .57 . 7 . 5 .
. . 8 .5577 . 38 . . 3
. . 83 .5 5
v: le degré de liberté. v =n-p-1 avec n: le nombre d'observations,
p: le nombre de variables explicatives. e : le risque d'erreur.

- -
 Ceci nous permet de regrouper les corrélations de la manière suivante :
Groupe 1 : Les paramètres avec un coefficient de corrélation supérieur à 0.9, ont une très
grande relation linéaire.
Groupe 2 : Les paramètres avec un coefficient de corrélation compris entre (0.80-0.90),
ont une grande relation linéaire.
Groupe 3 : Les paramètres qui ont un coefficient de corrélation entre (0.70-0.80), ont une
relation linéaire signifiante.
Groupe 4 : Les paramètres avec un coefficient de corrélation entre (0.60-0.70), ont une
bonne relation linéaire.
Groupe 5 : Les paramètres qui ont un coefficient de corrélation entre (0.40-0.60), ont une
Relation moyenne à acceptable.
Groupe 6 : Les paramètres avec un coefficient de corrélation inférieur à 0.40, ont une
faible relation linéaire.

Tableau 19. Matrice de corrélation des paramètres physico-chimiques des eaux

souterraines de l’aquifère Tadjenanet-Chelghoum Laid

Variable pH Cond Ca Mg Na K HCO SO Cl NO NO NH PO

pH
Cond -

Ca -

Mg -

Na
K
HCO
SO -

Cl -

NO -

NO -

NH -

PO -

Le coefficient de corrélation est couramment utilisé pour estimer la relation entre deux
variables. C'est un outil statistique simple pour montrer le degré de liaison entre deux
variables .L’examen de la matrice de corrélation, nous a permis d’établir quelques
corrélations significatives entre les différents éléments (fig.77).

- -
- Cond → Ca , Mg , Na+, K+, HCO -, SO - , Cl-, NO -, NO -,NH +,PO ---.
- Ca → Cond, Mg , Na+, K+, HCO -, SO - ,Cl-, NO -,NH +.
- Mg → Cond , Ca ,Na+, K+ , HCO -¯, SO - ,Cl-, NO ¯,NH +.
- Na+ → Cond , Ca , Mg ,K+, HCO -, SO - ,Cl-, NO -,NH +.
- K+ → Cond , Ca , Mg , HCO -, SO - , Cl-, NO ¯,NH +.
- Cl- → Cond, Mg , Na+, K+ , HCO -¯, SO - , NO -,NH +.
-
- SO → Cond, Ca , Mg , Na+, K+, HCO -, Cl- ,NO -,NH +,PO ---..
- HCO → Cond, Ca , Mg , Na+, K+, SO - , Cl-.
¯

- NO ¯ → Cond, NO -, NH +, PO ---.
- NO - → Cond ,Ca , Mg +, Na+, K+, SO - , Cl-, NO ¯, NO ¯,NH +,PO ---.
- NH + → Cond ,Ca , Mg , Na+, K+, SO - , Cl-, NO ¯, NO ¯,NH +,PO ---.
- PO --- → Cond, SO , NO ¯, NO ¯, NH +.

Figure.77.A. Corrélation entre la conductivité électrique et quelques éléments

chimiques majeurs

- -
 Figure.77. B. Corrélation significative entre quelque élément chimique

La conductivité électrique des différents échantillons, montre une bonne corrélation avec
tous les éléments (R> 0.40). (fig. 77 .a). Ces fortes corrélations indiquent la contribution
directe de ces éléments dans la salinité des eaux souterraines.

- -
Les relations significatives entre Ca et HCO3, Ca et Mg, Ca et Na, Ca et SO4, Ca et Cl

(fig.77.b), témoignent de la diversité des faciès chimiques des eaux souterraines de la

région.

5.2. Analyse en composantes principales

Le but de l’analyse statistique (ACP) est de faire ressortir les différents points d’eau qui
possèdent des concentrations élevées en éléments chimiques par rapport à l’ensemble des
points d’observation. Il sera ainsi possible d’avoir une idée sur le comportement
hydrochimique et les sites de pollution.
5.2.1. Choix de la méthode

La méthode regroupe le traitement de description et de comparaison simultanée des

différentes variables et l’analyse des données multidimensionnelles (statistique
descriptive).
Ce que nous cherchons à mettre en évidence, ce sont les liens éventuels qui pourraient
regrouper deux ou plusieurs variables au cours de leur évolution.
5.2.2 Résultats et interprétations
Le critère de Kaiser (Kaiser, 1960) a été appliqué afin de déterminer le nombre total de
facteurs significatifs. Selon ce critère, seuls les facteurs ayant une valeur propre
supérieure ou égale à 1 seront acceptés comme sources possibles de variance dans les
données (tab.20). Trois composantes principales (CP) ont été retenues et les résultats
montrent que trois CP représentent 70 % de la variance totale, ce qui est relativement bon
et peut être utilisée pour identifier les principales variations dans l’hydrochimie. (tab.21)
A % d’erreur le coefficient de corrélation critique de Pearson est de . . Sur cette base,
CP1 qui représente 55 % de la variance, est déterminé positivement par (Cond, Ca ,
Mg , Na+, K+, HCO3 -, Cl-, SO -
, NH et NO ). C’est le facteur de la contamination
par les éléments évaporitiques et carbonatés. Le CP2, qui représente 14 % de la variance
totale est déterminé par une opposition entre HCO3- et les éléments NO3- et PO ---
. Ce
dernier facteur est donc celui de la pollution agricole. Tous les éléments sont corrélés
négativement ave le pH. . La CP3 explique 9 % de la variance totale et il est déterminé
par le reste des variables (pH). Donc CP1, CP2 et CP3 sont supposés être représentatifs
du processus d’acquisition du chimisme de l’eau lié à l’interaction eau-roche et à la
pollution par l’activité agricole. (fig.78 A).
La projection des individus nous a permis de révéler trois familles de qualité d’eau, la
première et la deuxième correspondent aux zones de forte minéralisation qui sont localisées
- -
aux parties entre et l’est de la région d’étude. La troisième famille caractérise les zones où
l’utilisation des produits de fertilisation des sols dans l’agriculture est intense (fig.78.B).

Tableau Caractéristiques de l’ACP appliquée aux données physico-chimiques

des eaux de la nappe alluviale
CP Valeurs % de la Variance Valeurs propres Variance
propres du CP cumulée Cumulée
7. 8 55.7 7. 8 55.7
. 8 3 .75 . 7 . 3
. 55 3 .35 .37573 7 .8 33

Tableau 21. Corrélations variables – facteurs

Variable Fact .1 Fact .2 Fact .3 Fact .4

pH . 75 . 883 .8 78 .3 85 3
Cond - .8 7 - . 5 58 . 5 8 . 5 8
Ca - .87 . 37 - . . 878
Mg - .7 7 .33 737 - .38 5 5 . 335
Na - .8878 . 5 . 5 . 75 5
K - .75 . 5 3 . 575 - . 7 53
HCO - .5 85 7 .5 . 885 - . 3 3
SO - .883 . 7 5 - . . 77 5
Cl - .857887 .33 88 - . 55 5 . 5 8
NO - .58 88 - . 38 . 7 - . 37
NO - .77873 - . 5 5 - . 35 . 7 5
NH - .7 - . 35 . 758 - . 3 58
PO - .57 5 - . 33 7 . 88 3 - . 7

- -
 Plan factoriel (1

Plan factoriel (1

A* Projection des variables sur le plan factoriel

- -
 B* Projection des individus sur le plan factoriel (1

Figure 78. Cercle de corrélation et répartition des variables

6. Equilibres chimiques et indices de saturation

Degré de saturation d’une solution vis-à-vis d’un minéral
L’indices de saturation exprime le degré d’équilibre chimique entre l’eau et le minéral
dans la matrice de l’aquifère et peut être considéré comme une mesure du processus de
dissolution et ou précipitation concernant l’interaction eau-roche (Drever, 1997).
Le degré de saturation peut être évalué à partir de l’équation suivante :

IS = log (KIAP/Ksp). Où :

KIAP = le produit d’activité ionique des ions.

Ksp = le produit de la solubilité du minéral.
IS = l’indice de saturation.
- Si le IS < . L’eau est sous-saturée en ce qui concerne un certain minéral, c'est-à-dire
l’eau est encore capable de dissoudre ce minéral spécifique.
- Si le IS . L’eau est sursaturée en ce qui concerne ce minéral. Et le minéral va
précipiter
à l’intérieur de l’aquifère.
- Si le IS = . L’eau est en état d’équilibre.

- -
L’objectif principal des indices de saturation est d’évaluer les activités ioniques sous les
conditions de température et des forces ioniques prévues.

6.2. Etat des minéraux dissouts dans l'eau

Dans le but d’expliquer les mécanismes d’acquisition de la minéralisation, on a utilisé le
programme PHREEQC-2 (version 2.10) pour déterminer les indices de saturation de
quelques minéraux spécifiques (la calcite, de l’aragonite, de la dolomite, du gypse, de
l’anhydrite et l’halite).Tab. .
-les minéraux carbonatés présentent des degrés de saturation différents. En effet, celui de
la calcite varie de - .5 à . 5, suivi par celui de la dolomite qui varie de - . 5 à et
enfin celui de l’aragonite qui varie de - .7 à .3 .
Si l’on admet que l’état d’équilibre se fait dans l’intervalle de - 0.5 à + 0.5, on peut dire
que les trois minéraux sont légèrement sous-saturés et proche à état d’équilibre, (fig.79.a)
dont les minéraux carbonatés tendent à précipiter surtout sous la forme de calcite et de la
dolomite
- Les minéraux évoporitiques montrent des degrés de saturation plus faibles que les
minéraux carbonatés pour l’ensemble des échantillons analysés (sous-saturés) (fig.79.b).
Le gypse enregistre des indices qui varient de - . à - . , suivi de l’anhydrite dont
l’indice varie de - . jusqu’à - . 3 (83% des puits sont sous-saturés) et enfin l’halite
dont les indices varient de - . 8 à - . (donc très sous saturés) , ce qui permet aux ces
éléments de se présenter dans l’eau à des concentrations élevées. Cette simulation
confirme l'influence des roches évaporitiques sur la salinité des eaux analysées et
l'importance de l’halite dans le chimisme.

- -
Tableau. 22 Indice de saturation (SI) des minéraux présents dans l'eau

Minéral Minéraux carbonatés Minéraux évaporitiques

Calcite Dolomite Aragonite Gypse Anhydrite Halite
Pts d’eau CaCO3 CaMg(CO3) CaCO3 CaSO 2H O CaSO NaCl
.3 .7 . - . - . 7 - . 8
. 3 .3 - . - . - . 3 - .
. .5 . - .7 - . 3 - . 5
. 3 . - . - . - . 3 - .
. 5 .8 .3 - . - . - .
.3 .7 . 7 - .5 - .7 - .
. 5 . .3 - .7 - . - . 5
. .3 - . 5 - .3 - . - .5
- . 7 - . - . - .5 - .7 - . 7
- . . - . 5 - . - .37 - .73
- . 8 - . 7 - .33 - . 3 - .35 - .57
. . 5 - . 3 - . 8 - . - .3
- . - . 7 - . 8 - .83 - . 5 - . 8
- . - . 5 - . 8 - . - .3 - .
- . - . 8 - . - . - . - .
- . 3 - . 8 - . 7 - . - . - .77
- . - . - . - .87 - . - .78
- . 3 - .35 - . 8 - . - . 8 - .
- . - . - . 8 - . 7 - . - . 3
- . 5 - . - . - . 7 - . - .73
- . - .3 - .3 - .8 - . 8 - .
- . - . - . - . 5 - . 7 - .7
- . - . 7 - . - . - . - .7
. 3 - . - . - .8 - . - . 3
- . - . - . - . - . - .5
. . . - . - .88 - . 7
- . - . - . - .7 - . 8 - .
- . - . 8 - . - . 7 - . - .5
- .5 - . - .7 - . - .8 - .
- . - . - .3 - . - . 8 - . 7
- . 7 - . 3 - . - . - . - . 5
- .35 - .7 - .5 - .3 - .58 - . 3
- . 3 - . 5 - .57 - .38 - . - .8
Moyenne - . 788 - . 5 3 - . 7 - . 3 - . 7 - .5 7
Min - .5 - .7 - .7 - . - . - . 8
Max . 5 . .3 - . - . 3 - .
Ecar type . 5 . 5 . 5 7 .55 8 .5 33 . 8

- -
 a. Minéraux carbonatés

b. Minéraux évaporitiques

Figure Variation de l’indice de saturation des minéraux carbonatés et

évaporitiques dans les eaux souterraines de la nappe alluviale

- -
 Aptitude des eaux à l’irrigation

7.1. Introduction

L'irrigation avec des eaux riches en sels peut entraîner la fixation de sodium par le
complexe adsorbant du sol. L'intensité du processus de salinisation dépend des
caractéristiques du sol, de la qualité des eaux utilisées, des conditions de leur emploi et en
particulier de l'efficacité du système de drainage. La sécheresse de la dernière décennie a
engendré un déficit agricole (démontré par les calculs du bilan hydrique) suite à une forte
évapotranspiration, influençant le taux de salinité des eaux. Afin d’étudier l’aptitude des
eaux souterraine de la nappe alluviale à l’irrigation, nous avons procédé à utiliser les
méthodes de classification de Richards et Wilcox, les plus fréquemment utilisées.

7.2. Risque de sodicité et de salinité (Méthode de Richard)

La salinité peut entraîner des effets nocifs conséquents en raison de la fixation du sodium
des sels de chlorures de sodium par les colloïdes du sol. Le sodium exerce alors une
action néfaste sur la végétation, de façon indirecte, en dégradant les propriétés physiques
du sol. De part cette action, les sols deviennent compacts et asphyxiants pour les plantes
(Todd, 8 ).
Lorsque la concentration des ions Na+ à l’état soluble dans le sol est importante, ces ions
remplacent fréquemment Ca++ dans le complexe absorbant. Une eau chargée en sels peut
provoquer cet effet. Le risque est déterminé à partir de la valeur du sodium absorbable, «
Sodium Absorption Ratio» (S.A.R). Pour une même conductivité, le risque est d’autant
plus grand que le coefficient est plus élevé.
Le SAR est défini par la relation suivante :

S.A.R = Na++ / (Ca++ + Mg++)

Où Na, Ca et Mg sont exprimés en méq/l.

Ces éléments chimiques sont utilisés en combinaison avec la conductivité électrique de
l’eau.
Le report des valeurs du SAR en fonction de la conductivité électrique sur le diagramme
de Richard (fig.80) nous a parmi de dire que les eaux souterraines de la région,
appartiennent à deux classes. Il s'agit de C3-S1 de qualité admissible et C4-S1 de qualité
médiocre.

- -
Donc on peut dire que les eaux souterraines de la région d’étude peuvent être utilisées
sans risque pour le sol avec un risque de salinisation élevé dans la partie Est, et un risque
de sodicité faible. L’eau convenant à l’irrigation des cultures tolérantes au sel, sur des sols
bien drainés .

Figure 80. Classification des eaux pour l'irrigation selon la méthode de SAR

Les eaux de qualité admissible occupent la presque totalité du terrain alors que les eaux
de qualité médiocre se localisent au centre du terrain et aux alentours de Chelghoum Laid
et au nord de Tadjenanet (Feidh Nafaa). Fig. 81.

- -
 Figure Carte d’aptitude des eaux à l’irrigation selon la méthode de Richard
Mai 2007

Classification des eaux en fonction du pourcentage du sodium absorbé par l’eau

(Méthode de Wilcox)

Généralement il est connu que le développement des plantes est très faible ou nul dans les
sols saturés en sodium. En effet, dans un sol salé imprégné de NaCl, la perméabilité reste
bonne tant que la solution saline percolante présente une charge élevée en électrolyte,
charge qui maintient les argiles floculées (Servant 1975).
Si l'eau d'irrigation est chargée en sodium, le nombre d'ions Na combiné avec le sol,
augmente, cependant une quantité équivalente en cations et essentiellement en Ca sera
déplacée. On assiste alors à une dégradation de la structure, avec un gonflement des
argiles ce qui entraine une imperméabilisation du milieu. Le processus de dessalement et
de drainage risque alors d'être complètement stoppé. Dans le cas où le calcium est le
cation dominant, l'échange se fait dans le sens inverse, créant ainsi un sol plus floculé et
plus perméable. L'apport de CaSO au sol peut par l'échange de base, améliorer la texture
du sol ainsi que son drainage.
La classification de Wilcox fondée sur la conductivité électrique et la teneur en sodium
dans l’eau, généralement connue en pourcentage du sodium, est définie par la relation :
Na%= Na/( Ca+Mg+Na+K) Où tous les éléments sont exprimés en méq/l.

- -
Cette classification définit cinq classes d’eau : Excellente, bonne, admissible, médiocre,
et mauvaise.
La représentation du pourcentage de sodium en fonction de la conductivité sur le
diagramme de Wilcox (fig.8 ), permet de constater que 8 .8 % des points d’eau ont une
qualité bonne à excellente pour l’irrigation, % ont une qualité médiocre et % ont une
qualité mauvaise. On remarque qu’on a une dégradation progressive de la qualité des
eaux en se déplaçant de l'Ouest vers l'Est. Fig. 83.

Figure 82. Classification des eaux pour l'irrigation selon la méthode de Wilcox

- -
 335
Dj. G
rouz
ss
De
d
.e
330 Dj

Chelghoum Laid
325

320 Tadjenanet

315

310 ya
ou
Tik
.
Dj
305
780 785 790 795 800 805 810 815 820 825 830
Légende
Eau de qualité excellente Echèlle
0 2 4 km
Eau de qualité bonne
Eau de qualité médiocre
Eau de qualité mauvaise

Figure Carte d’aptitude des eaux à l’irrigation selon la méthode de Wilcox

Mai 2007

8. Etude de la potabilité des eaux

8.1. Introduction

Le but de l’étude de la potabilité est de définir les critères auxquels doit répondre une eau
destinée à l’alimentation humaine. En ce qui concerne l’eau de boisson, cette aptitude est
déterminée à partir des concentrations des anions et des cations majeurs.
8.2. La salinité

Traditionnellement les eaux souterraines ont été souvent classées sur la base des valeurs
de TDS (Davis et De Wiest, 1966).
La classification la plus souvent appliquée pour les eaux souterraines non fraîches est
présentée dans le tableau 23.
Tableau 23. Classification des eaux souterraines de la région de Tadjenanet-
Chelghoum Laid selon les valeurs de TDS (Davis et De Wiest,

Qualité de l’eau TDS (mg/l) % Echantillons Nombre d’échantillons

Douce - 8 .8 8
Saumâtre - 5. 5 5
Salée -
Saumure

- -
D’après cette classification on remarque que les eaux souterraines de la région sont en
général des eaux douces ( 0-1000 mg/l) avec 84.84% des cas tandis que les eaux
saumâtres représentent 15.15%.
8.3. La dureté totale
La dureté totale est exprimée en degrés hydrotimétriques (°TH) ou degrés français (°F).
Un degré français est la dureté d’une solution contenant 10 mg/l de CaCO3. En toute
rigueur la dureté totale ou titre hydrotimétrique d'une eau correspond à la somme des
concentrations en cations métalliques à l'exception de ceux des métaux alcalins et de l'ion
hydrogène. L'usage a limité cette notion à la dureté due aux ions calcium (dureté
calcique) et aux ions magnésium (dureté magnésienne) et ne prend pas en compte les ions
fer, aluminium, manganèse et strontium dont certains sont plutôt considérés comme
interférents de l'analyse. La dureté totale ainsi définie peut être calculée par :
TH (°f) = Ca (mg/l) / 4,008 + Mg(mg/l) / 2,431 Où :

TH : Dureté totale en °F
Ca et Mg : concentrations de calcium et du magnésium en meq/L.
Les classes suivantes sont généralement adoptées pour caractériser la dureté de l’eau.
(tab.24).

Tableau 24. Classification de l’eau selon la dureté totale

Dureté totale (°F) Qualité de l’eau
TH< 8F° Eau très douce
8≤ TH< 5 Eau douce
5≤ TH< 3 Eau moyennement dure

TH≥ 3 Eau très dure

Selon Rodier (1996), les eaux de qualité acceptable peuvent avoir une dureté jusqu’à 5
°F. A partir des données des eaux analysées, nous constatons 5 .5 % des points d’eau
analysés représentent une qualité acceptable (dureté inferieur à 50f°), caractérisent la
partie centrale et ouest du terrain.
8.4. Potabilité des eaux
Le report des résultats d’analyses physio chimique sur le diagramme de potabilité proposé
par Boudoukha et al., (1999) pour les zones arides et semi arides (fig.84) nous a permis
de constater que la plupart des eaux issues de la nappe alluviale du Mio-Plio-Quaternaire
offrent une eau potable d’une qualité admissible dans la partie ouest du terrain, et
médiocre dans la partie centrale et Est

- -
 Ca++ Mg++ Na+ K+ Cl- SO4-- NO3-
-
Cond - 2
- -

NON POTABLE

NON POTABLE
(µs/cm)
6
- -
- -100 -
---
10000 7
- -
--
1000 9
- - - 1000
- - - -
11

P O T A B I L I T E
- - - - 1000
- - - - 1000 18
- - - - - - -
- - - - 21
- - - -- 100
- - - 22
- -
--
- 100
- - - 23
- - - -
--

MOYENNE
- 27
- - - -

MOYENNE
- - - 28
- 1000 - -
- - - 10
-
- -
31

- --
100 32
- - 100 --
- - - - - 33
-
- -
-
-
- - -
-
100
-- 100 - 34
- - - - -
- - - 35
- - -
- -
-
- - - - 10 36
- -

BONNE
BONNE
- - 10
- -
- -
- -
- -
- 100 -1 - -

ADULTE

ENFANT
Teneurs en mg/l In A. BOUDOUKHA et all . 1992

Ca++ Mg++ Na+ K+ Cl- SO4-- NO3-

-
Cond - 37
- -
NON POTABLE

NON POTABLE
(µs/cm)
38
- -
- -100 - 39
---
10000
- - 40
--
1000
- - - 1000
- - - -
42
- -
- -
-
-- 1000 P O T A B I L I T E 43
- --
1000
- - - - - 44
- - - -
- - - - -- 100 46
- - - - 48
-
- --
100 -
- - - - 49
-- - -
MOYENNE

- - - - 50
MOYENNE

- - - -
- 1000
- 53
-
- - - 10
-
- - 57
--
100
--
--
- 100
- - - - 62
-
- -
-
-
- - -
-
100
-- 100 - 63
- - - - - -
- - - 67
- -
- - -
- - - - 10
-
BONNE
BONNE

- - 10 -
- - - -
- -
- -
- 100 -1 - -
ENFANT
ADULTE

Teneurs en mg/l In A. BOUDOUKHA et all . 1992

Figure 84. Diagramme de potabilité des eaux de la nappe alluviale Tadjenanet-

Chelghoum Laid –Normes Algérienne selon Boudoukha et al. 1999

- -
9. CONCLUSION

L’étude du chimisme de la nappe alluviale de Tadjenanet -Chelghoum Laid nous a permis

de conclure que le comportement hydro chimique est caractérisé par une grande
variabilité dans l'espace. La minéralisation augmente progressivement de l'Ouest vers Est
coïncidant avec le sens d’écoulement principal des eaux souterraines. Les eaux sont
caractérisées par une salinité faible à l’ouest près des monts de Djebel ed Dess et Djebel
Tnoutitt, alors qu'elles présentent au centre et à l’Est, une minéralisation assez élevée près
de la ville de Chelghoum Laid sur le long du oued Rhumel.
La conductivité électrique présente des variations importantes d'une zone à l'autre et
oscille entre 745 μ.S/cm et 3600 μ.S/cm avec une moyenne de 1465 μ.S/cm.
La répartition spatiale des concentrations en éléments chimique a permis de sortir trois
principaux types de faciès chimique, le faciès sulfatée calcique, chloruré calcique et
bicarbonaté calcique avec un degré faible le faciès chloruré magnésien et sulfaté
magnésien, dont le faciès chloruré calcique et bicarbonaté calcique appariaient à l’ouest et
le faciès sulfaté calcique au niveau du centre du terrain d’étude où les eaux possèdent
des grandeurs très élevés de la conductivité électrique , ce qui est expliqué par
l’accélération du phénomène de dissolution des roches évaporitiques et l’effet d’échange
de base entre les alcalins et alcalino terreux. Leurs répartition de la nature lithologique et
les conditions d’alimentation de l’aquifère.
Analyse statistique des données physico- chimiques de l’ensemble des points d’eau a
permis de montrer l’existence deux origine de polluants chimiques de nature diverse au
niveau des eaux de la nappe alluviale avec des concentrations élevées d’éléments
chimiques , l’un est dû à l’effets d’une pollution naturelle par l’effet de la salinité d’où la
forte concentration des eaux en éléments chimiques de Ca++,Mg++ ,HCO3-,Cl-, SO4--.
Le deuxième type de polluants est d’origine anthropiques marquée par la forte
concentration en NO3- dû essentiellement à l’utilisation intensive des engrais chimique
dans l’agriculture. Le calcul de l’indice de saturation des minéraux carbonatés et
évaporitiques a montré que les premiers ont tendance précipiter avant les seconds, ceci a
permis aux éléments chimiques provenant de la dissolution des gypses d’acquérir des
teneurs importantes. Concernant l'aptitude de ces eaux à l'irrigation, nous constatons que
les eaux de la zone Ouest de la région d’étude peuvent être utilisées sans risque sur le sol,
avec un risque de salinisation aux parties centre et Est du terrain d’étude dont l’eau
convenant à l’irrigation des cultures tolérantes au sel, sur des sols bien drainés.

- -
 CHAPITRE VI

VULNERABILITE A LA POLLUTION DE LA
NAPPE ALLUVIALE
 VULNERABILITE A LA POLLUTION DE LA NAPPE ALLUVIALE

1. Introduction

La vulnérabilité d’une nappe à la pollution est une notion relative, non mesurable et sans
dimension. La précision de son évaluation dépend essentiellement de la nature, de la
qualité et de la fiabilité des données utilisées. Sa caractérisation est basée généralement
sur l’estimation d’un certain nombre de paramètres plus ou moins importants relatifs
notamment à la couverture de la nappe (sol et zone non saturée) et à la zone saturée
proprement dite.
Les cartes de vulnérabilité montrent quelles sont, suivant les terrains rencontrés en
surface et les réservoirs contenant la première nappe d’eau souterraine, les possibilités de
pénétration et de propagation des polluants dans ces terrains. La vulnérabilité à la
pollution doit être considérée sous trois aspects :
- L’introduction de polluants dans les nappes, c'est-à-dire le transit, selon un
trajet surtout vertical, de corps polluants, entraînés par les eaux d’infiltration ou
des fluides polluants depuis la surface du sol jusqu’à une nappe souterraine, à
travers le sol et des roches non saturées d’eau ;
- La propagation de la pollution : l’écoulement de l’eau ainsi polluée, entraînée
par le mouvement naturel de la nappe atteinte, plus ou moins rapidement, selon
les caractéristiques de l’aquifère et la pente ;
- La persistance plus ou moins prolongée : après la suppression du foyer de
pollution, la durée de séjour est fonction du renouvellement naturel de l’eau dans
l’aquifère et également de la nature du polluant.
L’élaboration de la carte de vulnérabilité est basée sur la superposition de plusieurs cartes.
La carte finale permet de visualiser une zonation selon le degré de vulnérabilité (fort,
moyen et faible).

Les méthodes d’estimation de la vulnérabilité

Les méthodes d’estimation de la vulnérabilité des eaux souterraines à la pollution sont

très nombreuses. La multiplicité des Méthodes provient aussi bien des objectifs
recherchés (la réglementation, les activités Humaines, les polluants), que du nombre
d’organismes qui développent leur propre méthode pour une utilisation spécifique à leur
pays. Il est difficile de sélectionner une méthode d’évaluation de la vulnérabilité. Diverses
Raisons peuvent motiver le choix d’une méthode par rapport à une autre. Dans la présente

- -
étude, deux méthodes ont été appliquées : la méthode DRASTIC et la méthode des
indices de contaminations (Kherici 1993). (Tab.25).

Tableau Méthodes d’estimation de la vulnérabilité choisies

Méthodes
Paramètres DRASTIC Indice de contamination

Topographie X
Epaisseur / Texture du sol X X
Caractéristiques de la zone non X
saturée
Profondeur à la nappe X X
Caractéristiques hydrogéologiques X X
Perméabilité X X
Concentration des éléments X
chimiques

3. Vulnérabilité de la nappe alluviale par la méthode DRASTIC

3.1. Présentation de la méthode DRASTIC

La méthode DRASTIC a été mise au point par l’Environmental Protection Agency (EPA)
aux Etats-Unis en 85 et Aller et al. en 87, afin d’estimer le potentiel de pollution et
d’évaluer la vulnérabilité verticale des eaux souterraines. Elle se base sur sept paramètres
qui sont :
[D] : Depth to grouwndwater (Distance à la nappe, épaisseur de la zone non saturée)
[R] : Recharge de l’aquifère
[A] : Aquifer media (nature de la zone saturée)
[S] : Soil media (nature du sol)
[T] : Topography (topographie, pente en %)
[I] : Impact of the vadose zone (nature de la zone non saturée)
[C] : Conductivity (perméabilité de l’aquifère).
Chaque paramètre est cartographié avec un index, appelé aussi note (n), allant
typiquement de 1 à [Link] facteur de pondération (p) est ensuite appliqué aux différents
critères afin de relativiser leur importance respective en termes de vulnérabilité.
Une fois les différentes classes définies et leurs notes attribuées, la méthode détermine
l’indice DRASTIC (Id) qui permet de caractériser le degré de vulnérabilité d’un secteur

- -
donné de la nappe. La vulnérabilité est d’autant plus importante que l’indice (Id) calculé est
élevé. Cet indice est défini de la manière suivante:
ID = [Link]+[Link]+[Link]+[Link]+[Link]+[Link]+[Link]
où
D, R, A, S, T, I, C sont les paramètres cités plus haut.
p : facteur de pondération accordé à chaque paramètre.
n : notation accordée à chaque paramètre
Le tableau ci-après contient les poids attribués à chacun des paramètres ci-dessus. Les poids
les plus élevés (5) sont attribués à la nature lithologique et à l’épaisseur de la zone non
saturée. Les poids les plus faibles (1 et 2) sont attribués respectivement à la pente et à la
nature lithologique du sol. (Tab.26).

Tableau 26. Poids attribués aux paramètres DRASTIC (Aller et al, 1987)

Symbole Paramètre Propriétés Poids

D Profondeur de la Plus cette profondeur est élevée, plus le 5

contaminant met beaucoup de temps pour
nappe atteindre la surface piézométrique.

R Recharge Véhicule principal pour le transport du

contaminant. Plus cette recharge est
grande, plus le risque de contamination est
élevé.
Caractérisée par la granulométrie des
terrains saturés. Elle intervient dans le
piégeage du polluant qui peut s’échapper
A Lithologie de au pouvoir d’absorption du sol. Plus la 3
l’aquifère granulométrie est fine, plus le piégeage du
polluant est grand
Plus le sol est riche en argile, plus
l’absorption des métaux lourds est
S type de Sol importante, et plus la protection des eaux
souterraines est grande.
Plus la pente des terrains est grande, plus
le ruissellement des eaux est important et
T Topographie par conséquent la contamination des eaux
souterraines est faible
Son impact est déterminé à partir de la
texture des terrains qui la constituent. La
I Zone non saturée percolation du polluant jusqu'à la surface
5
piézométrique est d’autant plus grande que
cette texture est favorable
Plus ce paramètre est grand, plus le
transfert du polluant est rapide.
C Perméabilité 3

- -
-Le système de notation appliqué :
Pour chacun des paramètres DRASTIC, une note est affectée sur une échelle entière allant de
1 (le moins important) à 10 (le plus important). Pour un paramètre donné, une note basse
implique que la nappe concernée n’est pas vulnérable à la pollution et vice-versa (tab. 27).

Tableau 27. Classes et notes des sept paramètres (Lallemand- Barèes 1994)

D : Profondeur de la nappe (m) R : Recharge (mm)

Valeurs (m) Notation Valeurs (mm) Notation
0 – .5 5.5
1.5 – .5 17.5 – 5.5 8
4.5 – 7 10 – 7.5
- 5 5 5- 3
5- 3 -5
.5- 3
3
A : Nature de la zone saturée S : Nature du sol
Description Notation Description Notation
Calcaire Mince ou
karstique absent
Sable et gravier 8 Sables
Grès massif Limons
sableux
Métamorphique Terre grasse 5
altéré
Métamorphique 3 Limons silteux 3
Massive shale Argiles
T : Topographie (pente) I : lithologie de la couche vadose
Valeurs (%) Notation Description Notation
0à2 Calcaire Karstique
2à6 Sable et gravier
6 à 12 5 Sable et gravier avec silt et argile 8
12 à 18 3 Grès
8 Calcaire
Silt et argile
C : Conductivité hydraulique
Valeurs (m/s) Notation
, . -
,7. -
- , . -
8
3 , . -5 - ,7. -5

,7. -5 - 3 , . -5

,7. -5 - ,7. -5

,7. -7
- ,7. -5

- -
3.2. Description des paramètres et réalisation des cartes thématiques

Pour chacun des sept paramètres utilisés par la méthode "DRASTIC", une carte thématique
est réalisée. Sur chacune de ces cartes, sont délimitées les zones caractérisées par un indice
de vulnérabilité partiel du paramètre correspondant.

3.2 .1. Profondeur de la nappe (D)

La profondeur de la nappe est un paramètre important, qui constitue un facteur Déterminant

de la vulnérabilité. Plus cette profondeur est élevée, plus le contaminant met du temps à
atteindre la surface piézométrique. Pour la nappe alluviale de Tadjenanet –Chelghoum Laid,
ce paramètre a été déterminé pendant les relevés Piézométriques du mois de Mai 2007. Le
calcul des indices pour ce paramètre a été fait selon le tableau 28.

Tableau 28. Notation et pondération selon la profondeur

Intervalle (m) Dn Dp [Link]

1.5 – .5
4.5 –
- 5
5- 3
3-3

Ces indices ont permis d’établir la carte thématique du paramètre profondeur (fig.85).
L’observation de la dite carte montre que les plus faibles indices sont répartis dans la partie
sud-ouest prés de Djebel Tnoutit. Leurs valeurs ne dépassent pas 15, alors que les plus
importants indices se situent dans la zones centrale et Est du terrain d’étude, le long du oued
Rhumel, avec des valeurs de 35 à 5.

- -
 335
Dj. G
r ouz
ss
De
d
.e
330 Dj

325

320

315

310 ya
ou
Tik
.
Dj
305
780 785 790 795 800 805 810 815 820 825 830

Légende Echèlle
km
3 - 4.5 I.D = 45 0 2 4
4.5 - 9 Ville
I.D = 35
9 - 15 I.D = 25 Oued
15 - 23 I.D =15
> 23 I.D = 10

Figure 85. Carte indicielle : profondeur de l’eau Mai

3.2 .2. La recharge efficace (R)

La recharge efficace joue un rôle déterminant dans le transfert de la lame d’eau depuis la
surface du sol jusqu’à l’aquifère sous-jacent. Elle concerne la quantité d’eau qui franchit la
surface du sol pour renouveler les réserves d’eau souterraine et entretenir le débit de
l’écoulement souterrain, après circulation dans l’aquifère. Les résultats du bilan hydrologique
montrent que la valeur de ce paramètre est estimée à 7 mm sur l’ensemble de la région
d’étude. (Tab. )

Tableau 29. Notation et pondération du paramètre de la recharge.

Intervalle (cm) Rn Rp [Link]

- 5

L’indice DRASTIC est de l’ordre de répartie sur la totalité de la plaine (fig.8 ).

- -
 335
Dj. G
rouz
ss
De
d
.e
330 Dj

325

320

315

310 ya
ou
Tik
.
Dj

305
780 785 790 795 800 805 810 815 820 825 830

Légende Echèlle
km
Reharge = 7 mm Indice R = 4 0 2 4
0

Ville Oued

Figure 86. Carte indicielle : Recharge efficace

Nature de la zone saturée (A)

Ce paramètre désigne la lithologie de l’aquifère. Son identification a été basée sur les logs
des forages et des puits creusés pour le captage de la nappe. Les rapports de prospections
géophysiques (CGG 1973, ALGEO 1978) et les coupes géo-électriques interprétatives ont été
également exploités Ces documents ont montré que la nature lithologique de notre aquifère
d’étude est constitué essentiellement des alluvions, graviers, sables, argiles sableuses, limons
et calcaire lacustre.
Les indices correspondant à ce paramètre ont été estimés selon la notation illustrée dans le
tableau N° 30.

Tableau 30. Notation et pondération du matériau aquifère

Nature du milieu aquifère An Ap [Link]

Alluvions et gravier 8 3
Sable argileux 3
Calcaire lacustre 3 8

- -
L’examen de la carte des matériaux de la zone saturée montre que les indices les plus
importants se situent à la périphérie de oued Rhumel et les plus faibles au centre et l’Est du
terrain d’étude (fig.87).
335
Dj. G
r ouz
ss
De
d
.e
330 Dj

325

320

315

a
310 uy
iko
.T
Dj

305
780 785 790 795 800 805 810 815 820 825 830
Légende Echèlle
km
Alluvions et graviers Indice A = 24 0 2 4
Ville

Calcaire lacustre Indice A = 18 Oued

Sable argileux Indice A = 12

Figure 87. Carte indicielle : Matériau aquifère

La nature du sol (S)

La nature de la surface du sol influe sur la pénétration des polluants vers l’aquifère. Plus le
sol est riche en argile, plus l’absorption des polluants est importante, et plus la protection des
eaux souterraines est grande. Ce paramètre joue un rôle important dans la détermination des
zones polluées.
Les informations recueillies à partir des rapports et des cartes pédologiques établis sur les
hautes plaines sétifiennes (A.N.R.H de Constantine ),ainsi que la consultation des études
agro-pédologiques ont permis de constater que la nature du sol de la région d’étude est de
type calcimagnésique, occupe la totalité de la plaine, avec texture fine et parfois grossière
(terre grasse) . À proximité des bordures, près des massifs carbonatés et au centre du terrain,

- -
le sol est de type peu évolué d’apport colluvial, caractérisé par une épaisseur très mince. Les
valeurs indiciaires ont été déterminées selon les différentes natures de sol (tab.31).

Tableau 31. Notation et pondération de la nature du sol

Nature du sol Sn Sp [Link]

Terre grasse sableuse
Terre grasse argileuse 3
Sol mince

L’examen de la carte des indices relatifs à la nature du sol (fig.88) montre que le sol est de
type calcimagnésique ave une texture grossière caractérisé par une dominance sableuse
repartie sur le long de oued Rhumel, représente un risque de vulnérabilité moyen .Le sol de
type peu évolué d’apport colluvial se répartit au abord des massif carbonatés entourés la
région d’étude représente un risque élevé de vulnérabilité. Le reste du terrain est caractérisé
par un sol plus au moins argileux de texture fine.

335
Dj. G
r ouz
ss
De
d
330 D j. e

325

320

315

a
310 uy
iko
.T
Dj
305
780 785 790 795 800 805 810 815 820 825 830
Légende Echèlle
km
Sol mince 0 2 4
Indice S = 20
Terre grasse sableuse Indice S = 12
Terre grasse argileuse Indice S = 6

Ville Oued

Figure 88. Carte indicielle : Type du sol

- -
 La topographie (T)
Plus la pente des terrains est grande, plus le ruissellement des eaux est important et par
conséquent la contamination des eaux souterraines est faible. Les valeurs de pentes ont été
estimées en se référant aux cartes topographiques au 1/50000 de Tadjenanet et Chelghoum
Laid, et elles sont présentés dans le tableau N ° 32.

Tableau 32. Notation et pondération du facteur de la pente.

Intervalle (%) Tn Tp [Link]

0– 3 3
2– 3 7

L’examen de la carte indicielle montre deux plages de pente. La première varie de à %

qui occupe la majorité de la plaine, tandis que la seconde varie de 2 à 6 %, caractérise les
zone de la limite topographique du terrain d’étude (fig.8 ).
335
Dj. G
r ouz
ss
De
d
.e
330 Dj

325

320

315

310 ya
ou
Tik
.
Dj
305
780 785 790 795 800 805 810 815 820 825 830
Légende Echèlle
km
0-2% Indice T = 30 0 2 4

2-6% Indice T = 27

Ville Oued

Figure 89. Carte indicielle : Pente du terrain

- -
 Nature de la zone non saturée (I)
L’impact de la zone non saturée est considéré comme étant un paramètre très important dans
l’application de la méthode DRASTIC avec un poids (Ip = 5). Ainsi cette partie du sol est le
siège de nombreux phénomènes tels que l’oxydation, l’infiltration, l’évaporation …etc. Son
fonctionnement, ressemble à celui du milieu aquifère, autrement dit que l’atténuation des
contaminants est liée principalement à la granulométrie des matériaux meubles ou au degré
de fracturation de la roche.
La consultation des différents documents géologiques et géophysique (cartes et coupes), nous
a permis de dire que le milieu non saturé est constitué par un matériau détritique de
granulométrie très hétérogène, repartie sur la totalité de la plaine d’étude (sable et gravier
avec limon et argile. Les indices correspondant à ce paramètre ont été estimés selon la
notation illustrée dans le tableau N° 33.
Tableau 33. Notation et pondération de la zone non saturée.

Intervalle (%) In Ip [Link]

Gravier et sable argileux 5 3
Argile sableux 5

Selon la carte des indices relative à l’impact de la zone non saturée (fig. ) , on constate que
la partie centrale et Est du terrain d’études , sur le long de oued Rhumel possède un indice
très élevé (Indice I = 30) .

- -
 335 Dj. G
rouz
ss
De
d
.e
330 Dj

325

320

315

310 ya
ou
T ik
.
Dj
305
780 785 790 795 800 805 810 815 820 825 830

Légende
Echèlle
Gravier et sable argileux Indice I = 30
0 2 4
Argile sableux Indice I= 20

Ville Oued

Figure 90. Carte indicielle : Impact de la zone non saturée

3.2 .7. Conductivité hydraulique (C)

La perméabilité de la couche aquifère nous renseigne sur la vitesse de propagation des
polluants dans la nappe. Plus ce paramètre est important, plus le transfert des polluants est
rapide. L’estimation des indices de ce paramètre a été faite selon les intervalles présentés
dans le tableau N° 34.
Tableau 34. Notation et pondération de la conductivité hydraulique

Intervalle (m/s) Cn Cp Cn. Cp

,7. -5
- ,7. -5
3
,7. -7
- ,7. -5
3 3

L’établissement de la carte de perméabilité de l’aquifère alluviale (fig. , Chapitre IV.)

A mis en évidence deux zones de conductivité hydraulique .La première zone est
caractérisée par une variation de la conductivité hydraulique entre varie entre 10- et 10-5 m/s
qui s’étale sur la partie sud-ouest du terrain d’étude ; Ce qui donne un indice DRASTIC 3. La

- -
deuxième zone se caractérise par des valeurs de conductivité hydraulique varient entre 10 -5 et
-
caractérise le reste du terrain (fig.91).
335
Dj
.G
ss ro
uz
De 31 27
330 . ed 28
Dj 32

9
33
325 36 76
44 35 34
37
4240
43 39 23
38 22 11
46 2
21 18
320 49 48

57 50
53
62
63
315

67 a
uy
iko
310 D j .T
j. D
Tn
ou
tit
305
780 785 790 795 800 805 810 815 820 825 830

Echèlle
Légende 0 2 4 km
-5 -4
4.7. 10 m/s < C < 12.2. 10 m/s Indice C = 6
-6 -5
2.10 m/s < C < 4.7.10 m/s Indice C = 3

Ville Oued

Figure 91. Carte indicielle : Conductivité hydraulique

3.3. Carte de synthèse (vulnérabilité à la pollution)

La superposition des sept cartes thématiques relatives aux paramètres DRASTIC, nous a
permis de calculer l’indice DRASTIC (Id) à partir de quelle nous pouvons d’obtenir une carte
de vulnérabilité de l’aquifère en question.
Selon la classification adoptée par le ministère de l’environnement du Québec trois classes
de degré de vulnérabilité différentes peuvent être distingués :
- 150 < I : vulnérabilité forte
- 100 < I < 150 : vulnérabilité moyenne ;
- I<100 : vulnérabilité faible.

- -
L’examen de la carte de vulnérabilité de l’aquifère étudié (fig. ), nous a permis d’identifier
un indice DRASTIC qui varie entre 55 et 124, e qui met en évidence une nappe à
vulnérabilité moyenne à faible à la pollution
-La zone de vulnérabilité moyenne, occupe la partie centrale et Est du terrain d’étude sur le
long de oued Rhumel dont l’indice de vulnérabilité varie de à 5 qui est due
essentiellement la faible profondeur de la surface piézométrique et à la nature de la couche
vadose. Elle est située dans les aires influencées par les activités agricoles et les rejets
urbains.
- la zone de faible vulnérabilité, se localise dans le reste du terrain d’étude. Son indice de
vulnérabilité varie de 55 à 100.

335
Dj. G
s rouz
s
De
d
.e
330 Dj

325 Chelghoum Laid

320 Tadjenanet humel

O ue d R

315

a
310 uy
ko
. Ti
Dj
305
780 785 790 795 800 805 810 815 820 825 830

Légende Echèlle
km
0 2 4
I.D : Indice DRASTIC

100 < I.D < 150 : Vulnérabilité moyenne Ville

I.D < 100 : Vulnérabilité faible Oued

Figure 92. Carte de vulnérabilité à la pollution des eaux souterraines de la nappe

alluviale de Tadjenanet – Chelghoum Laid selon la méthode DRASTIC

- -
 Vulnérabilité de la nappe alluviale par la méthode de l’indice de contamination

4.1. Présentation de la méthode

L’étude hydrochimique des eaux souterraine de l’aquifère alluviale, nous a conduit à dire que
les concentrations très élevés en éléments chimiques sont liés principalement au lessivage des
roches magasins et à l’activité humaine (épandage d’engrais, rejets urbains industrie ….etc.)
Afin d’étudier la contamination des eaux de la nappe, nous avons appliqué la méthode basée
sur le calcul d’un indice regroupant l’ensemble des éléments chimiques considérés comme
étant d’origine anthropique.
Cette méthode consiste à admettre des intervalles de classes pour chaque élément. La somme
des classes pour chaque élément donne un indice final, (Kherici. 3), qui est d’autant plus
élevé, que le point de prélèvement est contaminé et donc vulnérable à la pollution.(tab.35)
L’indice est défini comme suit :

IC= Σ([Na+]+ [SO -- ]+ [Cl- ]+ [NO -])

Tableau Intervalle des classes de l’indice de contamination

Classes
Na+ (mg/l) < 5 5 -5 5 -75 75 -

Cl- (mg) < 5 5 -5 5 -75 75 -

--
SO (mg) < 5 5 -5 5 -75 75 -

NO - (mg/l) < 5 5- 5 5-35 35- 5 5-55 55- 5 5-75 75-85 85- 5

D’après le tableau de calcul de cet indice, les forts indices caractérisent essentiellement la
partie centrale et la partie Est du terrain, liés essentiellement aux teneurs très élevées en
nitrates, quant aux faibles indices, ils sont localisés sur les bordures et la partie ouest du
terrain (tab.36).

- -
 Tableau.36 Indice de contamination

N° pts d’eau
IC
N° pts d’eau
IC

Figure Carte de l’indice de contamination

4.2. Carte de la vulnérabilité à la pollution

Le but de l’établissement de cette carte est de montrer la susceptibilité de la nappe aquifère

à la pollution, afin de mieux apprécier la fragilité des captages d’eau potable existants, de
mieux cerner les zones à risque lors du choix de nouveaux sites des captages, et de prévoir
les mesures nécessaire pour lutter contre les risques de pollution chronique ou accidentelle,
diffuse ou ponctuelle .
La combinaison des différentes cartes , de la perméabilité, de la profondeur du plan d’eau, de
la nature du sol, la pente du terrain et celle des indices de contamination (concentration des
éléments chimiques), nous a permis de dresser la carte de vulnérabilité de la nappe alluviale
sur la base du paramètre principale de l’indice de contamination (Kherici. 3).

- -
La carte finale issue de la méthode de l’indice de contaminations (fig. ), fait ressortir trois
secteurs de faible, de moyenne et de forte vulnérabilité.
D’après la carte, Les zones de faible occupe presque la surface totale de la plaine d’étude,
alors que les zones à forte vulnérabilité se localisent au centre du terrain et les zones à
moyenne vulnérabilité s’étalent du centre à l’Est du terrain d’étude sur le long de oued
Rhumel.
335
Dj. G
rouz
ss
De
d
.e
330 Dj

Chelghoum Laid
Chelghoum Laid
325

320 Tadjenanet

315

310 ya
ou
Tik
.
Dj

305
780 785 790 795 800 805 810 815 820 825 830
Légende
Echèlle
Faible vulnérabilité 0 2 4 km

Moyenne vulnérabilité

Forte vulnérabilité

Figure Carte de vulnérabilité à la pollution selon la méthode de l’indice de

contamination

5. Conclusion

L’étude de la vulnérabilité de la nappe alluviale basée sur l’établissement de la carte

d’estimation de la sensibilité de la nappe d’eau souterraine à la contamination de manière
générale selon la méthode d’estimation DRASTIC et de l’indice de contamination a mis en
évidence une zone de forte vulnérabilité au centre la plaine ; cela nous amène à conclure que
cette zone est déjà contaminée, alors que les zones à moyenne vulnérabilité se localisent au
centre, à l’Est sur la périphérie de oued Rhumel, le reste du terrain à une faible vulnérabilité.

- -
 CHAPITRE VII

IMPACT DES EAUX USEES SUR LA QUALITE

PHYSICOCHIMIQUE DES EAUX
SOUTERRAINES ET EAUX DE SURFACE
 IMPACT DES EAUX USEES SUR LA QUALITE PHYSICOCHIMIQUE
DES EAUX SOUTERRAINES ET EAUX DE SURFACE

1. Introduction

Le développement industriel et la concentration de la population, font que le rejeté des

différents polluants dans le réseau hydrographique, crée des situations critiques et néfaste sur
l’hygiène et la santé humaine. Le rejet des eaux usées domestiques, industrielles ou agricole,
ont évidemment des conséquences considérables sur l’environnement. Toute étude
concernant la qualité des eaux passe obligatoirement par la connaissance précise des
différentes sources de pollution potentielle. A cet effet nous avons établi un inventaire des
différentes sources susceptibles de contribuer à la pollution que ce soit des eaux de surfaces
ou souterraines.

Pollution d’oued Rhumel

Les facteurs de pollution des eaux de oued Rhumel sont multiples et liés essentiellement à la
vocation agricole de la région, la nature des différents rejets domestiques et industrielles.
Nous signalons que le déversement des eaux usées de la ville de Tadjenanet et Chelghoum
Laid se font au niveau d’oued Rhumel, ce dernier en période d’étiage se transforme en égout
à ciel ouvert. En période de déficit hydrique, le pompage intensif à la périphérie de oued,
engendrent localement des axes d’écoulement dans le sens de oued-nappe, transmettant
ainsi la pollution vers la nappe souterraine.
Nous signalons que les deux grandes villes dans notre région d’étude (Tadjenanet et
Chelghoum Laid) se situent le long d’oued Rhumel avec une population estimée plus de
. habitant. Les deux agglomérations possèdent des décharges non contrôlées où l’une
est située à l’amont de oued Rhumel (environ km au nord-ouest de Tadjenanet).
Devant cette situation environnementale alarmante et le manque des stations d’épuration, on
a essayé de contribuer à l’étude de l’impact des eaux usées sur la qualité des eaux de surface
et souterraines par le bais d’une analyse physicochimiques.

- -
3. Différentes sources de pollution

3.1. Pollution agricole

L’usage de plus en plus généralisé d’engrais pose actuellement un grave problème de

pollution agricole. Il est notamment induit par les activités agricoles par l’épandage des
fertilisants chimiques, en particulier les engrais azotés et les phosphates qui ont pour
conséquence un enrichissement des sols et des eaux en nitrates et en phosphates. Il faut
signaler également l’activité liée à l’élevage du bétail, les fosses à purins non étanches, le
fumier, les eaux de lavage d’étables …etc
La région de Tadjnanet-Chelghoum Laid est caractérisée par une agriculture de type
céréaliculture et maraîchères qui nécessite une irrigation permanente surtout pendant la
période de basses eaux. Cette dernière se fait par des forages et puits de faible profondeur, et
à partir de pompage des eaux usées non traitées, s’écoulant dans l’oued.
Pollution d’origine domestique
Il s’agit d’une pollution aussi diversifiée qu’importante. La vie domestique contribue, pour
une part significative, à la dégradation de la qualité des eaux de surfaces et par la suite à la
qualité des eaux souterraines, nous distinguons la pollution domestique et la pollution
biologique, par les différents rejets d’eau (eaux ménagères, fosses septiques, station
d’essence, abattoir, déchets hospitaliers…). Le déversement des rejets domestiques dans la
région, se fait directement dans l’oued Rhumel sans contrôle.
Pollution d’origine industrielle
L’industrie a également une part prépondérante dans la pollution d’oued Rhumel, en
particulier le complexe des détergents et produits d’entretien (ENAD) de Chelghoum Laid,
les unités industrielles installés le long de la route N°5 en amont de oued Rhumel en
particulier les unités de vente et fabrication de trait à soudé et limonaderies.
4. Caractérisation de la pollution
Sites d’échantillonnage
Afin d’étudier l’impact des eaux usées sur la qualité physico-chimique des eaux souterraines
et de surface (eau de barrage Hammam Grouz), nous avons procédé à des prélèvements le
long de oued Rhumel selon une direction d’écoulement des eaux de surface et souterraine,
Ouest Est. Le premier prélèvement a été effectué à l’aval de la ville de Tadjnanet à Karir, le
- -
deuxième a été effectué à l’aval de la ville de Chelghoum Laid à Mzara, et le dernier
prélèvement a été effectué au niveau de la cuvette du barrage Hammam Grouz. Fig. 95 et 96.

335 Dj. G
rouz
ss
De
d
.e
330 Dj
Vers Constantine Barrage

325 Chelghoum Laid S3

S1 RN.5 S2
Tadjenanet el
320 Rhum
Vers Sétif Oued

315

a
310
o uy
ik
.T
Dj
305
780 785 790 795 800 805 810 815 820 825 830

Légende
Echèlle
Eau usée
----- Oued
0 2 4 km
Route Eau de barrage Hammam Grouz
Ville
Zone industrielle
Décharge public

Figure Carte d’inventaire des points de prélèvement des eaux de surface. Mai2007
et les différentes sources de pollution dans la région d’étude

- -
Station 1. Karir (Amont de oued Rhumel ) Station 2. MZARA (Aval de oued Rhumel)

Station 3. Barrage de Hammam Grouz Complexe des détergents - Chelghoum

Figure 96. Photos représentant les stations de prélèvements des eaux de surface

4.2. Résultats des analyses physicochimiques

Les éléments chimiques qui ont été dosés sont : chlorures (Cl), sulfates (SO ), nitrates (NO3),
nitrites (NO ), phosphates (PO ), fer (Fe), cuivre (Cu) et demande chimique et biochimique
en oxygène (DCO et DBO5). Les paramètres physiques qui ont été mesurés sont : la
conductivité, pH et le degré hydrotimétrique (TH).
D’après les résultats d’analyse regroupés dans le tableau n° 37, on constate que les
concentrations dépassent la norme telle que les sulfates et les chlorures avec la présence des
éléments métallique toxique (Fe et Cu). Ces fortes teneurs caractérisent une pollution des
eaux, traduisant par de fortes valeurs de la DCO et DBO5

- -
 Tableau Résultats d’analyses chimiques des eaux de surface

Cond DCO /
Station pH TH(°f) Cl SO NO NO PO4 Fe Cu Cd DCO DBO
(µ[Link]) DBO
S1
7.5 53 . 3 . . 3 . . 53 . 3 5 55 .
(Karir )
S2
8 7. 8 53 .5 .7 . 5 . 3 . 8 . 5
(Mzara)
S3
77 8. 5 5 8 3 . 5 . . 7 . 3 8 .
(Barrage)
Norme .5 -
5 3 5 5 . 5 . . 5
(OMS) 8.5
Teneurs en mg/l

5. Evolution de la pollution

D’après l’évolution spatiale des éléments chimique, nous constatons une légère élévation de
la salinité des eaux de surface de l’amont vers l’aval, puis elle diminue au niveau du barrage
Hammam Grouz sous l’effet du phénomène de la dissolution. Fig. 7. Pour le reste des
éléments, on s’aperçoit qu’on a une augmentation de tous les éléments ce qui montre une
forte pollution en liaison avec la nature des effluents. Les fortes valeurs de DCO et DBO5,
montrent la présence d’une charge polluante importante mobilisée à partir des eaux usées
rejetées dans l’oued Rhumel vers les eaux de barrage de Hammam Grouz. Cette constatation
a été déjà faite dans les études antérieures (Boudoukha et Boulaarak, 2013).
Le rapport de la demande chimique en oxygène et la demande biochimique pendant 5 jours
nous a permis d’apprécier la dégradabilité de la matière organique (AFNOR .1971).
Il est proche de 1.5 pour les eaux vanne, de 2 pour les eaux usées urbaines, et entre 3 et 5
pour les effluents issus des stations d’épuration et des zones industrielles. Les analyses des
eaux de la région, montrent le rapport DCO/DBO5 est légèrement supérieur à 2, ce qui
explique l’absence de la matière organique non biodégradable et confirme l’arrivée d’une
pollution urbaine.

- -
 350
Ouest Est
300
Concentration (mg/l)

250 NO3
NO2
200
PO4
150
Fe
100 Cu
50 DCO

0 DBO5
S1 (Karir) S2 (Mzara) S3 (Barrage)
Stations de prélèvement

Figure.97 Evolution spatiale des éléments de pollution des eaux de surface sur le long
de Oued Rhumel

6. Interaction eaux de surface- eaux souterraine

L’analyse des résultats physicochimiques des eaux souterraines et des eaux de surface, nous a
conduit à constater que la salinité des eaux de surface est semblable pratiquement aux
valeurs moyenne de celle des eaux souterraines de la nappe alluviale.
En terme des élément de pollution, les fortes concentration des éléments azoté trouves dans
les eaux de la nappe ne sont expliqués, que par leurs origines anthropiques par l’effet directe
de l’infiltration des eaux de surface chargés en sels (rejets domestique et industriel) ou par le
lessivage des terrain dans la zone non saturée (engrais azotés), ce processus est

- -
accentusurtout dans la période d’étiage, caractérisant principalement la partie centrale et Est
sur le long de oued Rhumel , ou il y a un pompage intense de la nappe .(fig.98)

Figure.98 Interaction eau de surface –eau souterraine à l’aval d’Oued Rhumel

8. Conclusion

D’après les concentrations élevées des indicateurs de pollution, nous constatons que les eaux
de surface montrent un apport conséquent en matières polluantes d’origine naturelle par le
lessivage des terrains géologiques ,et d’origine anthropique par les différentes nature de
rejets déversée à l’amont de l’oued, surtout de nature urbaine.

- -
CONCLUSION GENERALE
 CONCLUSION GENERALE

Au terme de e travail, nous pensons avoir contribué à une assez bonne connaissance
des paramètres intervenant dans l’étude de la caractérisation hydrodynamique et
hydrochimique de la nappe alluviale de Tadjenanet –Chelghoum Laid, ainsi que leurs
impacts environnementaux.

Cette étude a permis de dévoiler et maîtriser les phénomènes qui se posent dans la
région d’étude, et d’assurer une meilleur exploitation et une bonne protection de la nappe
contre la pollution.

La région explorée fait partie de la zone occidentale du grand bassin versant de Kébir
Rhumel situé à l’Est Algérien.
Du point de vue géologique la région étudiée est caractérisée par une formation de
type sédimentaire dont la prédominance des affleurements correspond à des formations Mio-
Plio-quaternaires. La description lithologique a montré une hétérogénéité dans les formations
(alluvions, argiles, sables, grès, limons, et les évaporites riches) caractérisant l’aquifère de la
région de Tadjenanet-Chelghoum Laid, qui issue d’une eau fortement minéralisée. On note
également la présence des formations du crétacé et Eocène de nature carbonatés entourant la
région d’étude, avec quelque affleurement d’éocène au centre du terrain. Ces formations à un
intérêt hydrogéologique important, contribuant à l’alimentation de la nappe alluviale. Nous
signalons l’affleurement des formations triasiques au contact de l’Eocène, caractérisent les
monts de ed Dess et Tikouya au nord-ouest et sud-est consécutive.
L’analyse des données climatiques de la région, a mis en évidence un climat de type
semi-aride, avec des étés chauds et secs et des hivers froids et humides. Les précipitations et
les températures moyennes annuelles sont de l’ordre de .8 ml an et 5. 3°C
respectivement.
Le bilan hydrique est dans son ensemble est déficitaire, avec une moyen annuelle des
précipitations d’ordre de 37 . 5mm et température moyenne annuelle de 5. 3°C.
L’évapotranspiration réelle représente % des précipitations, alors que l’infiltration est
très faible et ne représente que 1.5% des précipitations.
- -
 L’étude géophysique exécuté dans la région d’étude a permis d’identifier la
géométrie et l’extension des différentes formations géologiques susceptibles d’être aquifères.
L’examen des coupes géoélectriques et des cartes d’iso résistivités ont mis en évidence
l’irrégularité morphologique du substratum, ainsi que l’identification de la géométrie des
trois aquifères, le premier dans les formations du Mio-plio-Quaternaire qui est l’objet de
notre thème d’étude, le deuxième dans les formations carbonatées d’Eocène et le troisième
caractérise les formations carbonatés du crétacé.

L’examen de la carte piézométrique, réalisée en Mai 7, indique un écoulement de

direction principal Ouest- Est coïncidant avec l’axe d’écoulement des eaux de surfaces de
oued Rhumel. L’alimentation de la nappe, est assurée essentiellement par les formations
carbonatées limitrophes du Djeble ed Dess,Tnoutit, Grouz et Tikouya. La variation
spatiotemporelle de la surface piézométrique est expliquée par le changement des paramètres
climatiques et aux effets dû à l’excès de pompage.
L’interprétation des données de pompage d’essai, a montré que les zones favorables à
l’implantation d’éventuel ouvrages d’eau, et s’avèrent intéressantes du point de vue
potentialité aquifère se localisent au sud-ouest, près des massif carbonatés de Djebel Tnoutit.
Pour une meilleur gestion et exploitation des ressources en eau de la nappe alluviale
de Tadjenanet –Chelghoum laid et de prévoir leur évolution future au cour et à long terme,
on a essayé de pratiquer une approche par modélisation numérique tridimensionnelle, en
régime transitoire fondé sur la méthode des différences finies. Le modèle ayant pour finalité
l’optimisation de l’exploitation des réserves, permettront de placer une limite d’exploitation
et de choisir les meilleurs sites potentiels d’exploitation.
Les résultats du modèle hydrodynamique en régime transitoire étalant sur une
période d’observation de ans, correspondent à la surexploitation de la nappe par
une multiplication des débits pompés, montre un bilan hydrique de la nappe plus au
moins équilibré traduit par un faible rabattement atteignant les 6 m, étalé sur la
partie centre et Est de la nappe avec des faibles grandeurs d’ordre à m
enregistrés à la partie ouest. A cet effet, nous proposons, l’implantation d’éventuel
ouvrage de captage aux extrémités des massifs carbonatés entourés la région

- -
d’étude à fin d’éviter une telle situation entrainant la pollution de la nappe surtout
dans sa partie centrale, sachant que cette partie est connue par sa vocation agricole.

L’étude du chimisme de la nappe alluviale de Tadjenanet -Chelghoum Laid nous a

permis de conclure que le comportement hydrochimique est caractérisé par une grande
variabilité dans l'espace. La minéralisation augmente progressivement de l'Ouest vers Est
coïncidant avec le sens d’écoulement principal des eaux souterraines. Les fortes
minéralisations observées sont attribuables à une mise en solution au niveau des sols ou dans
l’aquifère d’une part, et à l’activité anthropique d’autre part, liés principalement à
l’utilisation irrationnelle des engrais chimique dans l’agriculture. Cette situation a été
confirmée par les résultats de l’analyse en composantes principales.
La cartographie hydrochimique a montré une zonalité des concentrations d’éléments
dosés ; où on constate une élévation de la salinité d’amont en aval dans le sens de
l’écoulement.
Le calcul de l’indice de saturation des minéraux carbonatés et évaporitiques a montré
que les premiers ont tendance précipiter avant les seconds, ceci a permis aux éléments
chimiques provenant de la dissolution des gypses d’acquérir des teneurs importantes.
Concernant l'aptitude de ces eaux à l'irrigation, nous constatons que les eaux de la
zone Ouest de la région d’étude peuvent être utilisées sans risque sur le sol, avec un risque de
salinisation aux parties centre et Est du terrain d’étude.
L’étude de l’impact des eaux usées sur la qualités des eaux de surfaces du barrage
Hammam Grouz ,et d’après les concentrations élevées des indicateurs de pollution
organique , il semble que les eaux de surface du barrage Hammmam Grouz montrent un
apport conséquent en matières polluantes d’origine naturelle par lessivage des terrains
géologiques ,et d’origine anthropique par les différentes nature de rejets déversée à l’amont
de l’oued, ainsi que une pollution agricole dû au épandage en masse des engrais chimique
sur des surface étendue localisée en amont du barrage.
Dans le but d’évaluer la vulnérabilité de la nappe alluviale de la région de Tadjenanet-
Chelghoum Laid à la pollution, et pour préserver cette ressource, on a adapté à établir deux
cartes de vulnérabilité relative aux indices DRASTIC et de contaminations calculés. Ces

- -
dernières tiennent compte de plusieurs facteurs, lithologique, piézométrique et
hydrochimique.
L’analyse des deux cartes de vulnérabilité a mis en évidence une zone de forte
vulnérabilité au centre la plaine ; cela nous amène à conclure que cette zone est déjà
contaminée, alors que les zones à moyenne vulnérabilité se localisent du centre à l’Est sur la
périphérie de oued Rhumel, le reste du terrain à une faible vulnérabilité.

Finalement, il est évident que ces synthèses sont relatives à une période d’observation
très limitée, et pour une étude très approfondie, nous proposons de l’effectuée à travers un
suivi spatiotemporel, par la mise au point des modèles hydrodynamique et hydrochimique
étalant sur des périodes d’investigation et d’exploitation diverses, qui vont renseigner, non
seulement, sur des informations sur l’origine des eaux, mais aussi sur l’état du réservoir et la
qualité de ses eaux.

- -
 RECOMMANDATIONS

Selon les articles relatifs à la loi de 2003 suivi par la lettre (b), le code des eaux a pour but
d’assurer la protection des eaux contre la pollution, le gaspillage, la surexploitation et
prévenir les effets nuisibles (Art 1).A cet effet un certain nombre de recommandations peut
être proposés, afin de lutter contre la dégradation de la qualité des eaux souterraines et de
surfaces :
- Arrêter, immédiatement, les nouveaux creusements de puits dans les secteurs où la nappe
alluviale est vulnérable à la pollution, surtout au niveau de la partie centrale et Est, et
modérer les pompages dans les puits déjà existants ;
- Adopter de nouvelles techniques d’irrigation (le goutte à goutte par exemple) plus
économiques de l'eau (pas de pertes par évaporation, fuites,…) que le système gravitaire, le
plus répandu dans la région
- Raisonner les apports en engrais azotés pour éviter les sur-fertilisations qui sont très
répandues dans la région ;
- Maîtriser le stockage et l’utilisation des déjections animales ;
- Eviter l’implantation des ouvrages d’AEP le long des berges de oued Rhumel, afin d’éviter
le flux imposé oued -nappe, ce qui amènerait une modification de la qualité des eaux de la
nappe.
- Réaliser un traitement préalable des rejets avant qu'ils atteignent l’oued et par la suite le
barrage, par la mise en place des stations d’épuration à l’aval des villes de Tadjenanet et
Chelghoum Laid.

- -
REFERENCES BIBLIOGRAPHIQUES
 REFERENCES BIBLIOGRAPHIQUES

Anonyme :
A.N.R.H : Agence Nationale des Ressource Hydrique
Fiche technique des forages. [Link].
Carte topographique et géologique de Tadjenanet, n°[Link] 1/ 500000.
Carte topographique de Chelghoum Laid, n° 119. Echelle 1/ 50000.
Carte géologique de Sétif ,1/ 200000.
Carte et rapports pédologiques des hautes plaines Sétifiennes .A .N.R.H de Constantine.
AFNOR ,(1971). Détermination de la demande biochimique en oxygène (D.B.O). Norme T
90-103. A.F.N.O.R., Paris, 5 p.
ALLER .L et al , (1987). DRASTIC: A standardised system for evaluating groundwater
pollution potential using hydrologic settings. U. S. EPA, Robert S. Kerr Environmental
ResearchLaboratory, Ada, OK, EPA/600/2-85/0108, 163 pp.
A.L.G.E.O (Société Algérienne de Géophysique), (1977). Etude par prospection
géophysique dans la région de Chelghoum Laid.
Astier. J.L, (1971) .Géophysique appliquée à l’ hydrogéologie. Edition Masson&Cie. Paris.
Bakalowicz . M, (1979). Contribution de la géochimie des eaux à la connaissance de l’aquifère
karstique et de la karstification. Thèse d’état Université P. et M. Curie Paris VI.
Barroin .G, 1990. La pollution des eaux par les phosphates. La recherche, 221.p 620- 7.
Bensouilah .S, (1995). Contribution à l’étude hydrogéologique des hautes plaines sétifiennes
dans le cadre de la haute vallée de l’oued Rhumel en amont d’oued
Athmé[Link].184p.
Bonnet. T, (1982). Méthodologie des modèles de simulation en hydrogéologie. Thèse de
doctorat d‟état en science (géologie appliquée). INP lorraine – NANCY. document BRGM
N°34.
Boudoukha. A, Djaffer. R, Retima. A, Djabri. L et Souag. M, (1999). Proposition de
normes de potabilité des eaux domestiques pour l’Algérie du Nord. Eudil Info-EAU. Edit.
Engineering Environment Consult. n° 24. 4-7.

- -
Brinis. N, (2011). Caractérisation de la salinité d’un complexe aquifère en zone aride cas de
l’aquifère d’El-Outaya, region Nord-Ouest de Biskra .Algerie. Thèse de doctorat, Univ.
Mohamed Kheider – Biskra, 191p.
Castany. G, (1963). Traité pratique des eaux souterraines. Dunod Edit. Paris, 641p.
Castany. G, (1968). Prospection et exploitation des eaux souterraines. Ed. Dunod Paris
France 683p.
Castany. G, (1982). Principes et méthodes de l’Hydrogéologie. Dunod Edit. Paris, 33 p.
Chabour. N, (2008): Hydrogéologie des domaines de transition entre l'Atlas saharien et la
plateforme saharienne à l'Est de l'Algérie. Thèse de doctorat d’état .Université mentouri,
Constantine. 176p.
C.G.G (Campagne générale de Géophysique), (1973). Etude par prospection géophysique
des hautes plaines Sétifiènnes. DHW de Mila –Algérie. 53p.
Code de l’eau, Loi83.17 du 16 juillet [Link] officiel de la république
Algérienne.
Davis. S. N, R. J. M. De Wiest , Hydrogeology. New York, Wiley, 463 p.
Djabri. L, (1996). Mécanismes de la pollution et vulnérabilité des eaux de la seybouse.
Origines géologiques, industrielles, agricoles et urbaines. Thèse de Doctorat d’État. Université de
Annaba 280p.
Debieche .T. H, (2002) : Evolution de la qualité des eaux (salinité, azote et métaux lourds) sous
l’effet de la pollution saline, agricole et industrielle. Thèse de doctorat. Univ. De franche comté.
199p.
Djouhra.R, (1983). Contribution à l’étude hydrogéologique de la région de Ain- Azel, mémoire
d’ingéniorat .USTHB Bab Elzouar, p.
Drever .J .I, (1997). The geochemistry of natural waters. 3rd ed. New York – Prentice-Hall
INC. 1997; 436 pp.
Duchaufour. P,( 2001). Introduction à la science du sol. Sol, Végétation, Environnement.
Dunod.
Durand. D.M, (1969) .Mise au point sur la structure du Nord-Est de la
Berbérie .Serv.géol .Algérie,N.S.,Bull.n°39,p.89- 3 .

- -
Durozoy. G, Les ressources en eau des massifs calcaires crétacés de la région de
Constantine, étude d'hydrogéologie appliquée, direction de l'hydraulique et de l'équipement
rural, service des études scientifiques, Alger.
Duchaufour. P, (2001) : Introduction à la science du sol, végétation, environnement. 6ème
édition, Dunod, Paris, pp.286- .
Farah. A.S, (1991). Etude du comportement d'un oued en zone méditerranéenne semi- aride
d'Afrique du nord et ses causes naturelles et anthropiques, l'oued Rhumel, Algérie .Thèse de
doctorat,Univ. Orléans, 271p.
Fehdi. CH, (2006). Apport hydrogéochimique et isotopique dans la caractérisation des
mécanismes d‟acquisition de la salinité des eaux souterraines du complexe aquifère 3
Morsott- el Aouinet (N.E. Algérien).thèse de doctorat es science. Université Badji Mokhtar-
Annaba.188p.
Freeze. A.R, Cherry J.A, (1979). Groundwater. Prentice-Hall, Inc., Englewood Cliffs, NJ,
604 p.
HILL .M, (1991). Nitrates and nitrites in food and water. Ellis Horwood, 194 pp.; New
York.
Hem. J. D, (1985). Study and interpretation of the chemical characteristic of natural water.
US Geological Survey. Water-Supply. Paper 2254, 263 pp ; Weinheim.

ISSAADI. A, ( Etude hydrogéologique des massifs du Guerion et Fortars. Thèse,3ème

cycle, Univ. Haouari Boumedienne, Alger.
Julien. T, (2003). Action N moins de nitrate dans l’eau. Office fédérale de l’agriculture
(suisse) [En ligne] 2003 Disponible sur : [Link]
Kaiser. H. F, (1960). The application of electronic computers to factor analysis.
Educational and Psychological Measurement, 20, 141- 5 .
Kieken. M, (1962).Esquisse tectonique de l’Algérie du nord .Présentation de la carte
tectonique au 100.000ème .Pub . Carte .gèol. Algerie. N° série .Bul.n° 3 .
Kehew.A.E, (2000). Applied chemical hydrogeology, Prentice-Hall international (UK)
limited, London, 368p.
Khedidja. A, (2001). Etude de la vulnérabilité de l’aquifère des calcaires lacustres de la
région de Bir EL-Arch ([Link] Sétif) .thèse de magister. Université de Tébessa.87 p
- -
Khedidja.A, Boudoukha.A, (2013). Risk assessment of agricultural pollution on
groundwater quality in the high valley of Tadjenanet – Chelghoum Laid (Eastern Algeria)
Publication au Journal International. Science and engineering. Desalination and water
treatment /Vol.51, January 2013, pp.292- 3 .
Kherici .N, (1993) .Vulnérabilité à la pollution chimique des eaux souterraines d’un système
de nappes superposées en milieux industriel et agricole (Annaba – la Mafragh) NE de
l’Algérie. Thèse de Doctorat Es-Sciences Université de Annaba 170p.
Konikow.L.F , Grove. D.B, (1977). Derivation of Equations Describing Solute Transport in
Ground Water. U.S. Geol. Survey Water-Res. Inv. 77-19: 30 pp.
Konikow.L.F, Reilly.T.E, (1998). Groundwater Modeling In: The Handbook of
Groundwater Engineering [J.W. Delleur, ed.], CRC Press, Boca Raton 20:1- . .
L’Organisation Mondiale de la Santé OMS Les lignes directrices de l'OMS en
ce qui concerne la qualité de l'eau potable, mises à jour en 2006 sont la référence en ce qui
concerne la sécurité en matière d'eau potable.
Lacroix. M, (1991). Impact de l’irrigation sur un aquifère alluvial. Dynamique du système
Basse Durance : nitrates et isotopes (cartographie), piézométrie (modélisation
mathématique). Thèse Doct. Univ. 161 p.
Lallemand-Barrès. A, (1995). Méthodes de dépollution des eaux souterraines.
Manuels et Méthodes. Orléans. BRGM. pp 60102.
Macko. S.A, Ostrom. N.E, (1994). Pollution studies using stable isotopes. In :Lajtha K. And
Michener R., Stable isotopes in ecology. Blackwell,Oxford, pp.45- .
McDonald .J .M, Harbaugh .A.W, (1988).A modular three-dimensional finite-difference
groundwater flow model, techniques of water resources investigations of the U.S Geological
Survey Book,6, pp.586.
Mebarki. A, (1982). Le bassin du Kébir Rhumel .Ressources en eau et aménagement en Algérie.
Thèse 3 ème cycle, université de NancyII.303p.
Marsily De. G, (2004). Cours d’hydrogéologie. Université pierre & marie curie. 226p.
Piper. A. M, (1944). A graphic procedure in geochemical interpretation of water analysis.
Trans. Amer. Geophys. Union, Vol 25 (6): pp. 914-928; Richmond, Va.
Plagnes .V, (2000). Structure et fonctionnement des aquifères karstiques. Caractérisation par
la géochimie des eaux. Documents du BRGM 294.
- -
Richards. L .A, (1954). Diagnosis and improvement of saline and alkali soils, Agric.
Handbook 60; USDA, Washington D.C., 160p.
Ricard. J, (1979). Evaluation des ressources d'un aquifère d'après l'analyse du bilan en eau :
application au système aquifère de Tebessa-Morsott (Algérie) Université des sciences et
techniques du Languedoc de Montpellier 2. Thèse : Sciences de l'eau
Rodier. J, (1996). L’analyse de l’eau, eaux naturelles, eaux résiduaires, eau de mer. 8ème
édition, Dunod, Paris.
Rouabhia Aek, 2006 : vulnérabilité et risques de pollution des eaux souterraines de la nappe
des sables miocènes de la plaine d’El Ma el Abiod Nord-Est Algerien. Thèse de doctorat.
Université Badji Mokhtar- Annaba.197p.
Schoeller. H, (1962). Les eaux souterraines. Hydrologie dynamique et chimique.
Recherche, exploitation et évaluation des ressources Masson et Cie. Editions, 642 pages.
Schoeller. H, (1962). Les eaux souterraines: Géochimie, eaux du gypse et des terrains
salifères , Masson et Cie édition Paris IV, 362 pages.
Servant. J, (1975). Contribution à l’étude des sols halomorphes. L’exemple des sols sales du
Sud et du Sud-Ouest de la France. Thèse de doctorat d’Etat, Université de Montpellier.
S.N REPAL (Société Nationale de Recherche et d’exploitation des Pétroles en Algérie,
Région sud telliennes et atlas saharien, 19th Int. Géol.
Stumm. W, Morgan. J.J, (1996). Aquatic chemistry ,3rd [Link]- Interscience, New
York, 1022p.
TODD. K, (1980). Groundwater hydrology. Seconde Edition. New York, [Link] & Sons,
510 pages.
Travi. Y, Mudry. J, (1997). Méthode pour l’évaluation et la gestion du risque nitrate dans
les aquifères de socle de la zone sahélienne d’Afriques de l’[Link]éologie : 3-21p.
Thornthwaite .C. W, (1980). An approach toward a rational classification of climate.
Geogr. Review, 38 : 55- .
VILA. J.M, (1977). Carte géologique de l'Est d'Algérie, feuille n° N-0-3-4 de Sétif au 200
000ème, Service de la carte géologique d'Algérie.
VILA. J.M, (1980). La chaîne alpine de l’Algérie orientale et des confins Algéro-Tunisiens.
Thèse de Doctorat- es -sciences, Université Pierre et Marie curie, Paris VI.

- -
Wilcox. L.V, . The quality of water for agricultural use. US Dept Agriculture
TechBull 962, Washington DC.
Wilcox.L.V,(1 ).Classification and use for irrigation waters. USDA,
Circ.969,Washington, D.C., 119 p.
Wuensch, K. L. (2001). Principal component analysis. [document disponible en ligne].
[Link]

- -
ANNEXE
 Annexe : Inventaire des points d’eau (puits, forages) Mai 2007

N° Points Côte piézométrique

X(km) Y (km) Z (m) Prof (m) NS (m)
d’eau (m)
8 7. 5 3 .85 8 8 .8 8 3.
8 7 3 .5 8 3 5 835
3 8 3.85 3 .3 78 ,5 77
8 . 3 .7 75 5 7 5
5 8 8.3 3 3.3 7 3 7. 788.8
8 7. 5 3 3.8 773 3. 7 .8
7 8 .5 3 . 77 5 3. 7 8.
8 8 5.7 3 .7 738 .7 7 8.3
8 .3 3 5.5 7 5 5. 7 .8
8 3. 3 5. 7 5 .8 7 .
8 .7 3 .3 7 8 3. 7 .
8 .7 3 . 77 . 7 5.8
3 8 . 5 3 .3 7 8 5 . 758.
8 . 3 . 773 5 . 7 .
5 8 3 . 5 78 8 7.8 77 .
8 7 3 .8 8 7 . 7 .8
7 8 8. 3 .8 785 .5 7 3.5
8 8 . 3 . 78 3 . 77 .
8 3 8.3 87 ,5 .7 8 7.3
8 . 3 7.5 837 8.7 8 8.3
8 .7 3 . 8 7 8.8 7 7.
8 . 3 8 8 7.5 8 .5
3 8 . 5 3 . 8 . 8 .8
8 8 3 5.3 7 7 75 8.5 778.5
5 8 7. 3 5.8 8 5 55 7 5
8 . 3 8. 8 7 .3 83 .7
7 8 . 3 .8 8 . 85 .8
8 8 5. 3 8. 5 8 . 83 .8

- -
 8 3. 5 3 .5 87 85 .7 85 .3
3 8 .35 3 8.3 8 7 3 8
3 8 .35 3 .5 88 5 . 87 .
3 7 .55 3 8. 87 3 873
33 8 3 .8 85 . 8 .8
3 8 . 3 .7 8 ,5 5.3 8 5.7
35 7 8.75 3 3. 85 5. 83 .
3 7 8. 5 3 3.8 8 . 837.
37 7 8.3 3 . 85 7 3. 837.
38 7 7. 3 .3 85 5. 838.
3 7 5. 5 3 . 87 5 . 855.8
7 5. 5 3 .5 88 33 8.3 8 .7
7 . 5 3 .75 87 3 7.5 85 .5
7 .35 3 . 88 3 3. 8 .8
3 7 3.75 3 . 5 8 3 . 87 .
7 . 5 3 3. 38 7. 3 .
5 7 3.5 3 .7 8 3 8 3.7 85 .3
8 3 . 83 75 .3 8 .7
7 8 . 5 3 8.5 8 .7 8 .3
8 7 7. 5 3 .3 85 5 . 8 .
7 3. 5 3 8.8 8 8 5 7.3 8 .7
5 7 7. 5 3 7. 8 3 83
5 7 . 5 3 .5 83 5 .8 8 3.
5 8 .7 3 .8 8 5 . 8 8.8
53 7 5. 3 5. 8 5 5 85
5 7 . 3 5.8 8 8 7 8. 85 .
55 7 3.85 3 . 5 87 .8 8 5.
5 7 3. 3 8.55 875 8 8.7 8 .3
57 7 . 5 3 .7 87 7.5 8 .5
58 7 . 5 3 5,5 888 5,5 3. 87 .8
5 78 .5 3 . 5 8 5 7 8.5 87 .5
78 3 7.75 8 .7 88 .3

- -
 788. 3 . 5 88
787.5 3 .55 8 . 88 .
3 78 . 5 3 .5 3 3 . 87 .
787. 5 3 3.7 5 875
5 788. 5 3 3 . 88 .
7 . 3 ,75 7 5. 875.
7 7 .35 3 . 5 88 7 5 87

- -