Structuration profonde dans la région de Bir M’Chergua (Nord-Est de la Tunisie). Implications hydrogéologiques
Deep structure in the Bir M’Chergua region (North-East of Tunisia). Hydrogeological implications

A. Hamed Ferjani^1,2, R. Guellala^1,2, A. Amiri¹, A. Merzougui³, M.H. Inoubli¹

¹Unité de Recherche de Géophysique Appliquée aux Matériaux et Minerais, Département de Géologie, Faculté des Sciences de Tunis, Campus Universitaire, 1060 Tunis, Tunisie. Email: ferjaniasma2013@gmail.com, rihabguellala@gmail.com; ORCID ID: https://orcid.org/0000-0001-8307-9794, https://orcid.org/0000-0002-4680-291X, https://orcid.org/0000-0002-0972-8319, https://orcid.org/0000-0002-9329-2780

²Laboratoire de Géoressources, CERTE, Technopôle de Borj Cedria, Tunisie

³Commissariat Régional de Développement Agricole de Zaghouan, Tunisie. ORCID ID: https://orcid.org/0000-0002-1435-1261

IntroductionTOP

En Tunisie, l’importante pénurie en eau due à plusieurs années de sécheresse a touché énormément les activités économiques et domestiques; elle a été à l’origine des fréquentes coupures en eau enregistrées dans tout le pays. Ainsi, l’exploitation des eaux souterraines est devenue une nécessité pour combler le besoin en eau. Cependant, la réussite de cette stratégie demande une bonne compréhension de la structuration profonde et de son effet sur les caractéristiques des aquifères afin d’assurer le succès des forages à réaliser. Dans ce contexte, est réalisé le présent travail qui intéresse la région de Bir M’Chergua située au nord-est de la Tunisie (Fig. 1) à cinquante kilomètres de la capitale et qui constitue un important pôle économique comportant différentes activités industrielles, agricoles et touristiques. Par ailleurs, Bir M’Cherga est caractérisée par la station thermale du Jebel Oust connue mondialement pour ses soins thérapeutiques.

Ce travail mené en collaboration entre l’Unité de Recherches de Géophysique Appliquée aux Matériaux et Minerais, le Centre de Recherches et de Technologies des Eaux, le Commissariat Régional de Développement Agricole de Zaghouan et l’Office National des Mines, vise à établir une carte structurale dans la région de Bir M’Chergua qui servira comme un document de base qui permet d’orienter les programmes d’exploitation des aquifères dans cette région. Cependant, les traits tectoniques profonds responsables de la structuration d’un aquifère ne peuvent pas être déterminés seulement par les données de forages; c’est pour cette raison que nous avons eu recours à la méthode gravimétrique appliquée avec succès dans différentes études hydrogéologiques (Najine et al., 2006; El Goumi et al., 2010; Ouerghi et al., 2013).

Contexte géologiqueTOP

La région de Bir M’Chergua constitue un point de convergence des principales zones structurales de la Tunisie: la plateforme pélagienne (Sahel), la partie nord de l’Atlas (Atlas septentrional), l’Axe N-S et la faille de Zaghouan (Fig. 1a). Cette région est caractérisée par différents traits tectoniques (Fig. 1b); les plus importantes sont les anticlinaux de Jebel Oust, Jebel Aziz et Jebel Ez Zebbas, les synclinaux de Henchir el Morra-Jebel Es Sba, de Bled el Mengoub et de Khelidia –Ain el Djour et le faisceau de failles transverses de direction NW-SE qui tronque au Sud le dôme du Jebel Oust tout en mettant en contact les dépôts crétacés inférieurs et ceux d’âge Eocène (Fig. 2).

Du point de vue stratigraphique, les affleurements géologiques dans le secteur d’étude s’étendent du Trias au Quaternaire (Fig. 2a). Les dépôts triasiques caractérisant Jebel Aziz et Jebel Zebbas sont représentés par du gypse, des grès, des calcaires et de dolomies (Fig. 2b). La série jurassique identifiée dans les massifs d’Oust et Aziz est composée de calcaires avec des intercalations marneuses (Fig. 2b). Les dépôts liasiques constituent le cœur du dôme du Jebel Oust. Le Crétacé est marqué par des alternances marno-calcaires (Fig. 2b) qui couvrent une large partie de la feuille de Bir M’Chergua (Fig. 2a). L’Eocène est représenté par des calcaires intercalés d’argiles riches en gypse (Fig. 2b). La série oligocène épaisse de 275 m (Jauzein, 1957) est formée de grès, sables et argiles. Les sédiments mio-pliocènes sont composés de calcaires, argiles sableuses, grès et conglomérats. Le Quaternaire est caractérisé par une croûte calcaire, d’argiles à intercalations sableuses, des galets et d’alluvions éoliennes argilo-gypseuses (Fig. 2b).

Contexte hydrogéologiqueTOP

Les données de forages d’eau relatives à la région de Bir M’Chergua mises à notre disposition par le Commissariat Régional de Développement Agricole de Zaghouan montrent que cette région est caractérisée par différents aquifères dont l’âge s’étend du Jurassique au Quaternaire. Onze forages représentatifs sont sélectionnés et mentionnés sur la carte géologique (Fig. 2a).

Au niveau du Jebel Oust, les calcaires jurassiques à la base du forage 7 montrent des résistivités allant de 100 à 200 Ohm^.m (Fig. 3). Cependant, les eaux jaillissant de ce forage (Fig. 2a) sont caractérisées par une température de l’ordre de 50^° C et une salinité élevée de 16 g/l. Cette contradiction impose une question sur l’origine des eaux thermales du Jebel Oust.

Les dépôts crétacés qui affleurent largement dans la région de Bir M’Chergua renferment plusieurs aquifères composés de calcaires fissurés à différente perméabilité. Ceux d’âge barrémien et campanien –maestrichtien présentent de meilleures caractéristiques (Table 1) et sont ainsi les plus sollicités.

Les calcaires éocènes constituent un important aquifère particulièrement exploité dans la partie sud-est de la région de Bir M’Chergua. Etant donné la grande variabilité du taux de fracturation de ces calcaires, différents résultats sont obtenus au niveau des forages qui les traversent; le forage 9 (Fig. 1b) est négatif alors que le forage 8 situé à proximité (Fig. 1b) est productif avec un débit de 8,57 l/s et une salinité de l’eau de 1,66 g/l. Dans ce forage, les calcaires éocènes sont caractérisés par une transmissivité de 0,81 10⁻⁴ m²/s et des valeurs de résistivité déterminées à partir des enregistrements diagraphiques qui atteignent 200 Ohm^.m (Fig. 4).

Dans la zone nord-est de Bir M’Chergua, ce sont plutôt les aquifères oligocène et quaternaire qui sont captés. Au niveau du forage 1 (Fig. 1b), l’aquifère oligocène composé essentiellement de sables et de grès est caractérisé par une bonne qualité de l’eau (salinité de l’eau: 0,8 g/l) et un débit relativement fort de 18 l/s. Peu profond, l’aquifère alluvial quaternaire est exploité par des puits de surface. Il s’étend vers l’est à la plaine de Mornag.

Données et méthodologieTOP

650 mesures gravimétriques couvrent le secteur d’étude selon une maille presque régulière d’un point par km². Elles sont fournies par l’Office National des Mines et acquises en 2009 avec un gravimètre Scintrex CG3 capable d’offrir une précision de 10⁻³ mGals (Gabalda et al., 2003).

Le positionnement des stations gravimétriques est effectué par un système GPS bi-fréquentiels de type Leica. La projection adoptée est Lambert Nord en utilisant l’ellipsoïde de Clark 1880 et le datum Carthage.

L’anomalie de Bouguer a été calculée en tenant compte de toutes les corrections (correction à l’air libre, correction de Bouguer) et en adoptant une densité de réduction de 2,36 g/cm³ déterminée par de mesures directes sur des échantillons obtenus des affleurements des formations géologiques représentatives de la zone d’étude et de mesures indirectes par la méthode de Nettleton et celle de triplets (Schoeffler, 1975; Nely, 1986).

Afin d’exploiter au mieux les mesures gravimétriques et en extraire le maximum d’informations sur l’organisation du sous-sol de la région de Bir M’Chergua, différentes techniques qui ont généré différentes cartes ont été appliquées en utilisant le logiciel Geosoft Oasis Montaj 6.4. Geosoft est une société de développement des logiciels géophysiques et géologiques basée à Toronto, Canada et fondée en 1982 pa Ian MacLeod et Colin Reeves.

Les cartes d’anomalie de Bouguer, d’anomalie résiduelle, de la dérivée verticale et du gradient horizontal sont calculées.

La carte d’anomalie de Bouguer reflète toutes les variations de la densité en profondeur (Schoeffler, 1975; Blakely, 1996). Vu que notre objectif est l’étude des aquifères, nous allons nous focaliser essentiellement sur les hétérogénéités de masse proches de la surface.

La technique de prolongement vers le haut (Jacobsen, 1987; Zeng et al., 2007; Jacoby & Smilde, 2009) permet de filtrer les structures de faibles longueur d’ondes proches de la surface laissant apparaître les structures de plus grandes longueurs d’ondes (c’est un filtre passe bas); l’anomalie régionale est engendrée par des sources profondes caractérisées par une large extension horizontale alors que l’anomalie résiduelle résulte de sources situées près de la surface.

La dérivée verticale (Blakely, 1996) permet d’accentuer les anomalies gravimétriques de courtes longueurs d’ondes et de renforcer l’amplitude des structures géologiques proches de la surface. Cette technique tend à réduire le phénomène de coalescence des anomalies (Dubois et al., 2011).

La technique du gradient horizontal est très utilisée pour localiser les discontinuités en subsurface à partir de données gravimétriques (Fedi & Florio, 2001). Les maxima du gradient horizontal coïncident aux positions de ces discontinuités (Dubois et al., 2011).

La technique de déconvolution d’Euler, appliquée aux cartes magnétiques (Thompson, 1982) et gravimétriques (Reid et al., 2003) permet d’estimer la localisation des sources d’anomalies ainsi que leur profondeur.

Les solutions de la déconvolution d’Euler dépendent des paramètres suivants: l’indice structural (IS), la fenêtre (F) et la tolérance (T).

Dans notre étude, nous avons retenu les valeurs correspondant à IS = 0, F= 8 km et T= 20 % qui ont donné le meilleur regroupement des solutions. Par ailleurs, IS = 0 est généralement utilisé pour identifier les failles (Reid et al., 2003)

Résultats et discussionsTOP

Carte de l’anomalie de BouguerTOP

La carte de l’anomalie de Bouguer (Fig. 5) ne montre que des valeurs positives allant de 0 à 25 mGal malgré l’abondance des terrains quaternaires à faible densité dans la région de Bir M’Chergua (Fig. 1).

L’anomalie la plus élevée caractérisée par une géométrie subcirculaire couvre le dôme du Jebel Oust qui correspond à une structure anticlinale essentiellement crétacée à cœur liasique (Jauzein, 1957). Vers l’ouest, au niveau de la plaine à remplissage quaternaire située au nord du Jebel Mouine est identifiée une anomalie à forte amplitude qui peut être expliquée par l’effet d’un terrain dense plus profond. En effet, in situ, le forage d’eau 2 (Fig. 1b) montre sous une couche d’argiles épaisse de 16 m, des alternances marno-calcaires du Barrémien. Elle peut être interprétée également comme une structure soulevée qui engendre généralement, à l’inverse d’une structure affaissée, une réponse positive (Prieto, 1996). Une anomalie d’ampleur comparable est localisée à proximité au niveau de l’affleurement crétacé du Jebel Oued El Hessi.

La mise en évidence des structures hautes qui comportent des dépôts calcaires perméables est un résultat très intéressant pour l’investigation des eaux souterraines dans la région de Bir M’Chergua.

Au nord –ouest du secteur d’étude, une anomalie à amplitude élevée caractérise Jebel Fedja qui correspond à un anticlinal constitué par les calcaires du Campanien, du Maestrichtien et de l’Eocène. A l’extrémité sud-est, Jebel Bou Khouf exprime aussi une forte réponse qui reflète sa structure anticlinale et sa composition en calcaires campaniens. Il s’agit d’un anticlinal secondaire du synclinal de Khélidia – Ain El Djour qui occupe toute la partie orientale de la feuille de Bir M’Chergua (Fig. 1b). Cet anticlinal s’amortit rapidement vers le nord (Jauzein, 1957).

Les anomalies à faible amplitude couvrent la zone sud-ouest du secteur d’étude: la plaine quaternaire de Bir M’Cherga (Bled Smindja et plaine de Depienne) et la cuvette endoréique de Bled Boucha. Elles s’étendent aussi sur une grande partie de la zone septentrionale; au niveau des synclinaux de Bled el Mengoub (nord –ouest) et de Khelidia (nord-est) et à Henchir el Mengoub (centre) qui fait partie de la plaine alluviale d’oued Chafrou.

La carte d’anomalie de Bouguer reflète l’expression de la majorité des structures en surface et celles qui sont enfouies. Toutefois elle révèle une réponse similaire pour Jebel Ez Zebbas et Henchir El Morra malgré que le premier soit un anticlinal et le second est un synclinal. Ils expriment une anomalie à moyenne amplitude (couleur verte). De surcroit, cette carte ne montre pas une anomalie spécifique à l’anticlinal jurassique du Jebel Aziz; sa réponse est influencée par celles des plaines quaternaires qui l’entourent.

Cartes prolongéesTOP

Pour une meilleure caractérisation des différentes structures géologiques proches de la surface, nous avons réalisé différents prolongements vers le haut (Fig. 6) afin de soustraire une composante régionale. Des modifications significatives sont enregistrées pour un prolongement de 3 km (Fig. 6a); les fortes anomalies au niveau du Jebel Fedja et la plaine quaternaire au nord du Jebel Mouine disparaissent indiquant le caractère peu profond des sources perturbatrices à l’origine de ces réponses gravimétriques. Aussi l’anomalie forte caractérisant le dôme du Jebel Oust commence à se fusionner avec celle de Jebel Bou Kouf. La carte prolongée à une distance de 5 km (Fig. 6b) montre l’installation d’un gradient gravimétrique croissant de l’ouest vers l’est. A un prolongement de 10 km (Fig. 6c), les courbes isogales sont caractérisées par un comportement lisse mais avec l’existence d’ondulations. A 20 km (Fig. 6d), ces ondulations sont complètement atténuées et la carte résultante représente la composante régionale du secteur d’étude. Elle révèle un gradient gravimétrique positif orienté SW-NE; un excès de masse au nord-est et un déficit de masse au sud-est. Cette tendance peut être expliquée par l’amincissement de la croûte en allant du sud vers le nord de la Tunisie, mis en évidence à partir de la sismique profonde dans le cadre du programme de Geotraverse (Buness et al., 1992) et confirmé par la gravimétrie (Jallouli & Mickus, 2000).

Carte de l’anomalie résiduelleTOP

La soustraction de la composante régionale de l’anomalie de Bouguer a donné la carte résiduelle de la région de Bir M’Chergua (Fig. 7) qui montre des valeurs allant de -11à 12 mgal.

Comparée à la carte d’anomalie de Bouguer (Fig. 5), cette carte montre l’atténuation de l’anomalie positive du dôme du Jebel Oust qui peut s’expliquer par l’effet des dépôts alluvionnaires d’Oued Meliane. En revanche, elle révèle l’accentuation des anomalies positives situées au nord du Jebel Mouine, à Jebel Oued el Hessi et à Jebel Fedja. En outre, elle décèle des anomalies négatives bien définies à l’emplacement des dépressions quaternaires de Khélidia, Henchir el Mengoub qui reflètent la faible densité de leur remplissage. L’anomalie négative de la plaine de Bir M’Chergua (Bled Smindja et plaine de Depienne) et celle de la cuvette endoréique de Bled Boucha sont séparées par une anomalie positive à l’aplomb des affleurements crétacés supérieurs et éocènes de K. El Feurcha el kartine qui peut correspondre à une zone soulevée dont la ligne de crête constituerait une ligne de partage des eaux souterraines. Cette zone pourrait aussi jouer le rôle d’une barrière entre les eaux saumâtres de la Sabkha et les aquifères de la plaine de Bir M’Chergua.

L’anomalie négative de Bled El Mengoub s’atténue alors qu’à proximité à Jebel Ez Zebbas s’individualise une anomalie positive malgré la présence du Trias gypseux à faible densité à l’affleurement. Des anomalies positives caractérisant des structures salifères sont observées dans différentes localités de la Tunisie septentrionale telles que la haute vallée de Mejerda, Jebel Cheid, Jebel Mourra et Lansarine (Benassi et al., 2006; Hamdi Nasr et al., 2009; Amiri et al., 2011; Sarsar Naouali et al., 2017). Ces réponses indiquent la présence des roches denses comme l’anhydrite et la dolomie en subsurface. Par ailleurs, des anomalies négatives sont les réponses d’autres corps salifères coexistants (Inoubli & Mansouri, 2006).

Carte de la dérivée verticaleTOP

La carte de dérivée verticale de la région de Bir M’Chergua (Fig. 8) révèle l’individualisation et l’accentuation de certaines anomalies par rapport aux réponses enregistrées au niveau de la carte d’anomalie de Bouguer et même au niveau de la résiduelle. L’anomalie positive du Jebel Fedja s’étend plus au Sud indiquant la continuité de cette structure sous la couverture quaternaire. Celle du Jebel Ez Zebbas est remarquablement accentuée avec une forme subcirulaire et des limites bien définies. L’anomalie négative de Henchir el Mengoub s’allonge vers l’Ouest à Henchir el Morra et vers l’Est à Henchir Sidi El Hattab el Abbassi. A sa limite Nord, à Sidi Zouaoui, s’individualise une anomalie positive dirigée WNW-ESE qui peut constituer une structure haute composée de calcaires éocènes affleurant à proximité.

L’anomalie positive du dôme du Jebel Oust est bien matérialisée suivant une direction NE-SW surtout après sa séparation de l’anomalie du Jebel Bou Kouf. Ses amplitudes maximales sont vers la partie méridionale; elles couvrent l’affleurement crétacé inférieur du Jebel Er Rihana et le massif jurassique du Jebel Oust. Elle est faible dans la partie centrale du dôme où les affleurements sont essentiellement composés de dépôts quaternaires.

L’extension horizontale et l’orientation E-W de l’anomalie positive du Jebel Mouine sont mieux caractérisées malgré que sa réponse globale reste coalescente avec celle de l’anomalie du Jebel Oued el Hessi.

Carte de la dérivée horizontaleTOP

L’interprétation des maxima de la dérivée horizontale a permis d’identifier les discontinuités structurales majeures à l’échelle du secteur d’étude (Fig. 9a). Les alignements dégagés s’arrangent essentiellement suivant les directions E-W, NE-SW et NW-SE (Fig. 9b).

Linéaments de direction NE-SWTOP

La direction NE-SW est représentée principalement par les linéaments F1, F2, F3, F4 et F5. F1 et F2 délimitent respectivement le dôme du Jebel Oust à l’est et à l’ouest. La technique de la déconvolution d’Euler (Fig. 10) confirment les tracés de ces linéaments et décèlent leur caractère profond avec des solutions qui dépassent dans la majorité des points 1000m pour F1 et 1500 m pour F2.

F3 qui tronque la partie septentrionale de ce dôme est aussi enraciné; sa profondeur estimée excède 1500 m. L’alignement F4 borde à l’est la structure du Jebel Fedja alors que F5 tronque le synclinal de Henchir el Morra-Jebel Es Sba. Jebel Fedja est traversé par les failles NE-SW de Jerifet et el Fedja qui constituent des segments de l’accident régional de Tunis-Ellès à jeu actuel inverse senestre (Ben Chelbi et al., 2008). F4 et F5 peuvent être également interprétées comme des tronçons de cet accident.

Dans la partie méridionale de la carte, un alignement NE-SW peu étendu délimite à l’ouest la zone affaissée de Bir M’Chergua (Bled Smindja et plaine de Depienne).

Les linéaments NE-SW mis en évidence sont parallèles à la partie septentrionale de la faille plurikilométrique de Zaghouan qui prolonge la dorsale tunisienne (Rolland, 1885; Solignac, 1927; Castany, 1951; Turki, 1980; Creusot et al., 1992; Morgan et al., 1998) sur 80 km depuis Jebel Boukornine de Hammam Lif jusqu’à Jebel Bargou au Sud-Ouest. Cette faille correspond à un alignement composé de trois accidents distincts (Turki, 1980) qui sont du nord -est au sud-ouest: le chevauchement de Zaghouan proprement dit (partie septentrionale), l’accident subvertical du Jebel El Menassir-El Kemkine (partie centrale) et la faille d’effondrement du Jebel Serdj (partie méridionale).

Linéaments de direction E-WTOP

La direction E-W est matérialisée essentiellement par quatre alignements: F6 », F7, F8 et F9. L’alignement majeur F6, d’importante extension horizontale, marque la limite entre les structures hautes de Jebels Oued El Hessi, Mouine, Er Rihane et les structures affaissées de Bir M’Chergua et Bled Boucha. Seul Castany (1951) a signalé le tracé de cet alignement. La déconvolution d’Euler (Fig. 10) révèle que la profondeur de F6 varie de 1000 à 1500 m.

L’alignement F7 correspond à la limite méridionale des structures affaissées de Bled El Mengoub, Henchir el Mourra et Henchir el Mengoub et correspond à des profondeurs comprises essentiellement entre 1000 et 1500 m.

Le troisième alignement F8 d’extension horizontale relativement faible affecte le massif jurassique du Jebel Oust. Ce massif est également orienté E-W malgré qu’il constitue le cœur du dôme du Jebel Oust ayant une direction générale NE-SW (Jauzein, 1957).

Le quatrième alignement F9 est mis en évidence au niveau de l’affleurement crétacé couvrant la partie sud-est du dôme du Jebel Oust (Fig. 9b). Il est parallèle à plusieurs accidents identifiés en surface.

Linéaments de direction NW-SETOP

La direction NW-SE caractérise des alignements relativement nombreux mais d’extension latérale plus faible; il s’agit de F10, F11, F12 et F13 F14, F15 et F16.

L’alignement F10 coïncide en surface avec le réseau de failles transverses tronquant le bord sud du dôme du Jebel Oust (Fig. 9b) et délimitant au nord le fossé de Bir M’Chergua (Bajanik et al., 1977; Kadri & Ben Haj Ali., 1999).

F11 ayant une géométrie en escalier se prolonge sur 5 km de Jebel Mekhma à Sidi Ahmed avec des profondeurs majoritairement supérieures à 1500 m (Fig. 10). Cet alignement marque la limite entre l’affleurement éocène du Jebel Bou Hadjaba et celui d’âge oligocène du Jebel Mekhma qui renferment des dépôts perméables composés respectivement de calcaires et de sables. Il peut constituer une faille butoir qui entrave l’écoulement entre les aquifères oligocène et éocène ou au contraire une zone de circulation préférentielle des eaux souterraines et favorise l’intercommunication de ces aquifères. F12 coïncide avec la limite méridionale de l’anomalie positive décelée à Sidi Zouaoui (Fig. 8). F13 coïncide avec l’accident déjà reconnu en surface qui délimite à l’ouest Jebel el Mengoub (Jauzein ,1957). «F14» borde au nord le synclinal de Bled el Mengoub –Henchir el Mila. F15 sépare les zones affaissées de Henchir El Mengoub et Henchir El Morra (Fig. 9b). La mise en évidence de cet alignement confirme en partie l’accident supposé sur la carte géologique qui s’étend du Jebel el Mengoub jusqu’à Henchir Bou Ferda (Fig. 9b). F16 marque la limite occidentale de l’anomalie négative couvrant la dépression de Bled Boucha (Fig. 8).

Dans la partie nord – ouest de la carte, sont mises en évidence des linéaments NW-SE d’extension relativement faible qui sont transversales à F4. Au sud-ouest, le linéament qui relie les segments de F6 limitant au sud Jebel Mouine et Jebel Oued Hessi explique la direction NW-SE de l’anomalie positive exprimée par la carte résiduelle (Fig. 7) et celle du gradient vertical (Fig. 8).

Synthèse et conclusionsTOP

L’interprétation qualitative et quantitative des anomalies gravimétriques a permis d’approcher au mieux l’organisation des structures profondes de la région de Bir M’Chergua qui sont généralement associées aux jeux de failles de différentes directions.

Les anomalies positives du dôme du Jebel Oust, Jebel Fedja et Jebel Mouine correspondent à des anticlinaux résultant de la compression NW-SE du Miocène supérieur qui correspond à la phase paroxysmale de l’édification et l’installation de l’Atlas tunisien (phase atlasique) (Ben Ayed & Viguier, 1981; Bouaziz et al., 2002). Elles sont bordées par des accidents de direction NE-SW et E-W (Fig. 11). L’application de la compression N-S du Quaternaire ancien qui a donné la géométrie actuelle des structures atlasiques (Chihi, 1988; Mzali & Zouari, 2006) dans notre zone d’étude accentue les plis et réactive les accidents NE-SW et E-W en failles inverses avec une composante senestre pour les NE-SW.

Les anomalies négatives de Bir M’Chergua, Bled Boucha, Khelidia, Bled el Mengoub, Henchir el Mengoub et Henchir el Morra couvrent des dépressions d’accumulation de dépôts néogènes et quaternaires dont la plupart correspondent à des structures synclinales générées par la phase atlasique. Elles sont limitées par des linéaments orientés NE-SW, E-W et NW-SE (Fig. 11). Ces derniers répondent en failles normales à la compression N-S; les accidents tectoniques décrits en surface par Bajanik et al., 1977 et Kadri & Ben Haj Ali, 1999 comme étant les limites du fossé de Bir M’Chergua -Mograne sont les meilleurs exemples. Seul le réseau des failles traçant la bordure nord- orientale de ce fossé est déterminé dans la feuille de Bir M’Chergua; il correspond à F10 mis en évidence à partir de l’interprétation gravimétrique.

Le fossé de Bir M’Chergua est cité parmi les grabens transversaux au tracé de l’accident de Zaghouan (Turki, 1988).

A l’extrémité Sud –Ouest de la feuille de Bir M’Chergua, la faille inverse F6 et la faille normale F16 ont contribué conjointement à l’affaissement de la cuvette synclinale de Bled Boucha. La zone haute entre cette cuvette et la plaine de Bir M’Chergua correspond à la fois au compartiment soulevé de F16 à l’ouest et d’une faille inverse à l’est.

L’intersection de deux accidents profonds F6 et F16 identifiés au niveau du Jebel Aziz facilite l’ascension du matériel triasique.

Dans la partie orientale, le synclinal de Khelidia – Ain Djour est affaissé d’avantage sous l’effet du jeu inverse de F1 qui le sépare du dôme du Jebel Oust. Par ailleurs le jeu normal de F11 a effondré la partie nord-est de ce synclinal, caractérisée par des amplitudes plus faibles que la partie sud-ouest.

Dans la partie septentrionale, la dépression de Henchir El Halouan- Henchir El Mengoub est un compartiment affaissé à la faveur de trois accidents tectoniques: les failles normales F12 et F15 respectivement au Nord et à l’Ouest et la faille inverse F7 au Sud. L’effet conjugué de ces déformations a engendré l’approfondissement de cette zone et a favorisé son remplissage par des dépôts quaternaires à faible densité et peut ainsi expliquer l’importante anomalie négative qui la couvre. Le compartiment soulevé de la faille F12 coïncide avec la zone haute mise en évidence au niveau de Sidi Zouaoui. Le synclinal de Henchir El Morra correspond au compartiment soulevé de la faille F15. Les accidents E–W à NE–SW à jeu inverse F5 et F7 séparent ce compartiment et les structures hautes au Nord –Ouest, au Nord et au Sud.

L’anomalie négative de Henchir El Halouan- Henchir El Mengoub est très accentuée par rapport à celle de Henchir El Morra. Ceci peut être expliqué par une profondeur plus importante et par l’existence des sédiments quaternaires.

A l’extrémité Nord-Ouest de la carte, le jeu normal des failles F13 et F14 a générée une succession d’effondrements, traduite par une diminution progressive des réponses gravimétriques, depuis l’anticlinal du Jebel Fedja jusqu’au synclinal du Bled El Mengoub dont la limite Sud est tracée par la terminaison occidentale de l’accident majeur F7.

Le schéma structural de Bir M’Chergua ainsi établi (Fig. 11) montre des traits tectoniques extensifs et d’autres compressifs non reconnus auparavant. La mise en évidence de ces traits permet une meilleure connaissance de la structuration de l’Atlas septentrional en Tunisie.

La coexistence de différents traits tectoniques dans la région étudiée offre une complexité à la géométrie du son système aquifère. Le plissement des séries sédimentaires, la réactivation des accidents NE-SW et E-W qui délimitent les bombements anticlinaux de cuvettes synclinales en failles inverses et le jeu normal des accidents NW–SE sont à l’origine de la compartimentation de ce système en structures soulevées et affaissées. Les structures soulevées peuvent constituer des zones de partage des eaux ainsi que des zones de recharge préférentielles des dépôts perméables particulièrement au niveau des affleurements de calcaires crétacés qui constituent les principaux réservoirs dans ces zones. Les zones affaissées renferment généralement un aquifère multicouche; les couches crétacées atteintes à une profondeur relativement importante sont surmontées par d’autres tertiaires et/ou quaternaires.

Les structures reconstituées peuvent ainsi expliquer le changement de la profondeur d’une formation perméable entre deux points d’eau voisins tel est le cas des forages 4 et 5 (Fig. 11). En effet, l’aquifère barrémien est recoupé uniquement à 11m de profondeur par le forage 5 situé au niveau du compartiment soulevé de la faille inverse F6 alors qu’il est atteint à 40 m au niveau du forage 5 implanté dans le compartiment affaissé de cette faille (Fig. 12). Cette faille est enracinée et peut aussi générée des variations de profondeur ainsi que d’épaisseur dans les séries sédimentaires plus profondes, non atteintes par les forages d’eau 4 et 5.

Ces structures peuvent aussi jouer un rôle dans les relations entre les différents aquifères caractérisant la région et influencer la circulation et les caractéristiques des eaux souterraines. Les eaux thermales du Jebel Oust constituent un très bon exemple.

En effet, un modèle gravimétrique 2.5D de direction NNW-SSE qui recoupe le dôme du Jebel Oust (Fig. 13) a permis de déceler d’avantage la géométrie de cette structure et son effet sur l’écoulement souterrain.

La modélisation a été réalisée en utilisant le logiciel GM-SYS (Geosoft), les contrastes de densité ont été calculés par rapport à la densité de réduction des mesures gravimétriques. Les géométries des couches ont été construites en intégrant les différentes coupes géologiques construites par Solignac, 1927 et Castany, 1951 et en adoptant des densités moyennes représentatives des formations géologiques de la zone d’étude (Table 2).

Le modèle gravimétrique a été ajusté en comparant la réponse gravimétrique calculée à celle extraite de la carte résiduelle calculée d’un prolongement de 5 km (Fig. 13a).

Le modèle obtenu (Fig. 13b) révèle une concordance avec les résultats de l’interprétation qualitative; la forte anomalie coïncide avec la présence de calcaires jurassiques et crétacés. La faible anomalie correspond à la réponse des dépôts quaternaires de faible densité.

Les séries jurassique et crétacée du dôme du Jebel Oust sont affectées par des accidents tectoniques d’importante profondeur qui dépasse par endroits 1500 m. cette structuration explique l’origine des eaux thermales à forte salinité captées par le forage 7, implanté au niveau du Jebel Oust (massif jurassique); les accidents mis en évidence favorisent une pénétration profonde des eaux de recharge ainsi que leur contact avec les terrains triasiques riches en sels.

La présente étude a abouti à des résultats qui peuvent contribuer énormément à une exploitation rationnelle des ressources en eau souterraine de la région de Bir M’Chergua. L’intégration d’autres méthodes géophysiques telle que la sismique réflexion qui est actuellement inaccessible dans cette région ne peut qu’ améliorer ces résultats et apporter plus des précisions au schéma structural établi. Par ailleurs cette étude peut servir d’exemple pour d’autres projets ayant pour objectif une meilleure gestion des eaux dans les régions arides dans le monde.

REMERCIEMENTSTOP

Les auteurs remercient Mr. Farhat Ben M’Barek, Directeur du Service Géophysique à l’Office National des Mines, de leur avoir fourni les mesures gravimétriques. Ils remercient également les deux reviewers anonymes pour leurs remarques constructives.

RéférencesTOP