these - Université de Poitiers

2%. Des valeurs de transmissivité plus faibles à l’intérieur des terres par rapport à la zone
côtière justifient cette différence de gradient.
Dans ce travail une simulation du comportement de l’interface de la nappe avec l’eau de mer
a été abordée par la méthode d’éléments aux limites. Pour un scénario de recharge de 15 E6
m 3 /an et 90% de prélèvement soit 13.5 E6 m 3 /an et une profondeur de l’interface à 50 m à la
distance de 5 km de la côte, une remontée de l’interface de 1 m/an a été estimée. La vitesse de
remontée de l’interface diminue progressivement et s’arrête au bout de 40 ans, l’épaisseur de
la zone saturée se stabilisant à 25 m environ.
Dans le cadre de sa thèse sur les propriétés géométriques et hydrodynamiques des
aquifères en milieux volcaniques fissurés sous climat aride, Jalludin, (1993) a élaboré un
premier modèle en différences finies, avec un maillage de 500x500 m à l’aide du logiciel
DIFINE. Les conditions aux limites adoptées dans ce travail sont (figure 2-15) :
i) au Nord et à l’Est, les limites de l’aquifère correspondent au contact avec la mer. Un
potentiel imposé h=0 a été appliqué sur ces limites.
ii) au Nord-Ouest dans la région d’Arta et au Sud-est dans la région de Loyada, des lignes de
courant dessinées à partir de la piézométrie observée ont été prises comme limites à flux nul
du modèle.
iii) à l’Ouest, au contact avec les formations basaltiques de Dalha, une relation hydraulique
entre les deux aquifères est considérée. Des conditions aux limites à potentiels imposés ont été
fixées à partir des données piézométriques disponible sur l’aquifère basaltique de Dalha.
iv) au Sud une interface qui représentait le contact entre les basaltes du Golfe et les basaltes
Somali, d’après les cartes géologiques de l’époque, (les dernières études géologiques ont
montré que cette interface ne se trouve pas à cet endroit) a été considérée comme une limite à
potentiel imposé. Une continuité hydraulique est considérée entre ces unités basaltiques et la
limite a été construite sur la base de la piézométrie connue dans les basaltes Somali et les
basaltes du Golfe sur des points d’eau au voisinage de l’interface.
Le régime permanent de l’aquifère a été calé dans le but de vérifier la cohérence des données
existantes et des conditions aux limites appliquées. Ce travail a permis de dégager des zones
de transmissivité allant de 10 -1 à 10 -3 m²/s, de discuter les relations hydrauliques avec les
formations adjacentes. L’infiltration dans les lits des oueds a été prise en compte comme
recharge de l’aquifère.
Il est probable que les conditions aux limites à charges imposées considérées au contact avec
les basaltes Dalha dans cette modélisation en régime permanent restreignent la
représentativité de ce modèle.
48