MODELISATION DE L'AQUIFERE PLIO-QUATERNAIRE EN REGIME
PERMANENT ET BILAN
La modélisation hydrodynamique de l'aquifère
côtière du Cap-Bon est réalisée par l'application du
code MODFLOW . Le maillage adopté dans ce cas est un maillage
carré régulier de 250m de côté. La
discrétisation du domaine est composée de 104 colonnes et 208
lignes avec 6352 mailles de calcul, 1145 mailles à potentiel
imposé et 1299 mailles à débits imposés, soit au
total 8796 mailles actives (Fig.3). La dimension des mailles a
été choisie en tenant compte de la densité et de la
variabilité des données disponibles .
Figure 3: Représentation du maillage choisi
et des conditions aux limites considérées.
Les conditions aux limites adoptées dans la conception
du modèle sont : potentiel imposé constant égal à
zéro le long de la mer Méditerranée ; potentiel
imposé constant le long de la limite Nord constituée par l'oued
H'jarr à Kélibia (La charge hydraulique considérée
est égale à la différence entre le niveau topographique et
la hauteur du lit de l'oued. Dans le cas présent, la charge varie entre
2 et 3 mètres) et le long de la limite Sud de l'aquifère
constituée par l'oued El K'bhir à Béni Khiar (en adoptant
les même hypothèses que pour la limite Nord) ; débit
imposé constitué par les puits d'exploitation de la
région. La limite Ouest est une limite fictive dans laquelle on
considère l'écoulement latéral .
La distribution des perméabilités est
représenté dans la figure 4.
Figure 4 : Carte de la
conductivité hydraulique de l'aquifère
plio-quaternaire.
Pour le régime permanent, on va adopté un
coefficient d'infiltration estimé à 15% pour les affleurements
Plio-quaternaires dominants sur l'ensemble des mailles actives du domaine de
travail. La recharge actuelle de l'aquifère plio-quaternaire s'effectue
par l'infiltration :
- des ruissellements en amont essentiellement par les
puissantes séries gréseuses inférieures qui enserrent
l'Oligocène du jbel Abd Errahmen,
- des eaux de crues des oueds traversant la plaine, des pluies
excédentaires des années exceptionnelles.
L'évaporation directe est considérée non
nulle lorsque la profondeur de la nappe est inférieure à 10m. On
comptabilise un déficit pluviométrique d'environ 790mm/an.
Les prélèvements par pompages sont
supposés nuls en régime permanent (état
piézomètrique de1970).
Calage du modèle
L'état de référence retenu pour caler le
modèle en régime permanent est la piézomètrie de
1970 correspondant à l'état avant le début des pompages.
La carte (Fig.5) montre la distribution de la perméabilité
obtenue à l'issue du calage. On distingue principalement quatre zones
:
- au Sud d'oued H'jarr (Kélibia) où
l'aquifère circule dans les calcaires et les dolomies, la
perméabilité varie entre 10 et 12m/j,
- au voisinage de l'Oued Chiba et El K'hbir, la
perméabilité est faible, allant de 4 à 5m/j (les
dépôts d'alluvions limoneux peuvent expliquer ces faibles
valeurs),
- à proximité du barrage Chiba où la
l'aquifère circule dans les schistes, la perméabilité
varie entre 9 et 11m/j, avec 10m/j comme valeur prépondérante,
- le long de la bande côtière, elle est en
général de l'ordre de 8m/j, ce qui est légèrement
inférieur à l'ordre des perméabilités
mesurées.
Figure 5 : Carte de la perméabilité
après calage du modèle plio-quaternaire de l'aquifère
de la côte orientale.
Figure 6: Comparaison entre la carte
piézomètrique de 1970 observée et celle simulée par
la modèle en régime permanent.
Ce calage du modèle en régime permanent a permis
la reconstitution du bilan global de l'aquifère par évaluation de
ses différentes composantes (Fig.7).
ENTRÉES À LA NAPPE (+) :
Recharge par les oueds+barrage =
| 
