Bulletin de.l’Institut scientifique, Rabat, section Sciences de la Terre, 2001, n°23, 1-6.

Carte du gradient géothermique au Maroc

Abdelkrim RIMI

Institut Scientifique, Département de Physique du Globe, Av. Ibn Batouta, B.P. 703 Agdal, Rabat, Maroc. e-mail : rimi@[Link]

Résumé. Une nouvelle compilation de données sur les gradients géothermiques au Maroc est présentée en utilisant 461 températures de
fond de puits pétroliers (BHT) corrigées, 22 profils thermiques mesurés dans des forages de profondeur entre 75 et 500 m et en intégrant
les valeurs des zones frontalières. Une carte de courbes d’égales valeurs de gradient géothermique montre les tendances macroscopiques
de la distribution. Pour la comparaison, les sources chaudes sont aussi ajoutées sans tenir compte de la profondeur d’origine de l’eau. Dans
le nord-est du Maroc, le gradient élevé de 35 à 50°C/km et jusqu'à plus de 65°C/km en Alboran et le sud-est de l’Espagne est attribué au
volcanisme récent et l’amincissement de la croûte terrestre alors que dans le Prérif, l’anomalie de 35-45°C/km reste locale et s’explique
par la nature marneuse des formations et par la fracturation intense favorisant une circulation de fluides. Sur la marge atlantique, la
tendance d’élévation du gradient en direction des îles Canaries est associée au phénomène de hot spot à l’origine de leur naissance.
Mots clés : Maroc, géothermie, carte, gradient géothermique, source thermale

Geothermal gradient map of Morocco.

Abstract. A new compilation of thermal gradient data in Morocco and bordering regions has been attempted by using 461 corrected
bottom hole temperatures (BHT), 22 thermal profiles in shallow boreholes 75-500 m depth. A contour map shows the main spatial
variations of the thermal gradient where thermal springs are superimposed for comparison. In northeastern Morocco, the 35-50°C/km
thermal gradient range which exceeds 65°C /km in the Alboran sea and southeastern Spain is linked to the recent volcanism and crustal
thinning, while the 35-45°C/km anomaly in the Prerif remains local and may be explained by the marly formations and also by dense
fractures which generate fluid circulation. On the Atlantic margin the trend of thermal gradient increase toward the Canary Islands is a
consequence of the hot spot phenomenon at the origin of their formation.
Key words: Morocco, geothermics, map, geothermal gradient, thermal spring

INTRODUCTION inverser le gradient de température. La cartographie du

gradient thermique à l'échelle régionale présente donc
Depuis 1983, plusieurs travaux ont essayé d’estimer les
l’inconvénient d'avoir une résolution très limitée ; aussi, le
valeurs du gradient géothermique au Maroc. Les méthodes
traçage de la carte des variations régionales du gradient
adoptées varient depuis les mesures directes dans les
thermique doit tenir en compte le fait que le gradient peut
forages (Rimi & Lucazeau 1987, Zaghloule 1999) à des
varier latéralement dans le cas des structures d’ef-
évaluations indirectes à partir d’échantillonnages de
fondrement ou d’intrusion de granite ou de sel. Pour éviter
géothermomètres chimiques, principalement SiO2, dans des
cet effet, nous avions estimé les gradients thermiques
sources ou des forages (Bahi et al. 1983, Ziyadi 1993,
moyens des grands domaines géologiques du Maroc en
Boukdir 1994, Benabidate 1994, Lahrach 1994, Zaghloule
reportant, sur un même graphe pour chaque domaine, toutes
1999, et Cidu & Bahaj 2000). La composition chimique à la
les températures mesurées et corrigées dans les forages
résurgence doit être considérée avec précaution ; en effet,
pétroliers en fonction de la profondeur (Rimi 1990).
lors de son ascension, l’eau peut soit provenir directement
et rapidement des niveaux profonds et chauds, dans quel
cas la composition chimique reflète celle des formations
METHODOLOGIE
réservoirs, soit subir des transformations et mélanges
successifs caractérisant les différentes formations aquifères Dans le présent travail, pour pouvoir cartographier la
parcourues. distribution régionale et diminuer l'effet des variations
spatiales du gradient géothermique, les critères de choix
Des essais d‘élaboration de cartes régionales ont été
varient selon le type de puits.
effectués par Rimi et al. (1998) pour le Rif, Lahrach (1994)
pour le Maroc oriental puis par Boukdir (1994) pour le Dans les forages pétroliers, les températures de fond du
bassin de Tadla. Plus récemment, Zaghloule (1999) a dressé puits (BHT pour bottom hole temperatures) corrigées ont
une première carte du gradient géothermique au Maroc, été regroupées par site ; dans le cas d’une densité élevée de
sous forme de zonage de gradient en se basant sur le calcul puits (Gharb et Prérif), un gradient thermique moyen a été
de gradient moyen de toutes les températures dans chaque obtenu à partir de la représentation graphique de toutes les
forage. Cependant, une difficulté dans cette procédure vient BHT pour chaque site. Dans les régions où les forages sont
du fait que le gradient de température est par définition lié à relativement dispersés, le gradient est déterminé pour
la profondeur et à la nature de la roche, et donc un gradient chaque puits par la différence entre la température mesurée
moyen risque de ne pas représenter l’état thermique des à la plus grande profondeur et à la surface du sol, ce qui
formations en raison des contrastes de conductivité évite ainsi les niveaux intermédiaires qui peuvent être le
thermique et des circulations de fluide qui pourraient même siège d’une convection liée à des circulations d’eau.
 Rimi – gradient géothermique du Maroc

Espagne

e
tiqu Alb
ora
n
n d'
tla
r
Me
R if

nA
35

en
éa

oy
Oc
M eseta

M
A tla s
H au t

ie
A nti A tlas

r
gé
Al
30
B assin de T indouf

25

Mauritanie

-20 -15 -10 -5 0

Figure 1. Valeurs ponctuelles du gradient géothermique en ° C/km

Pour les sondages peu profonds, les critères de sélection appelée habituellement "analyse structurale" ou
sont la représentativité de l’emplacement du puits, "variogramme" (Fig. 3) ; elle décrit la moyenne des carrés
l’équilibre thermique et l’absence de mouvements des écarts entre des paires de données lorsque la distance
convectifs. Ainsi, 22 mesures de flux de chaleur (Rimi qui les sépare augmente. Enfin, pour éliminer les effets de
1999) dans des forages peu profonds ont été incorporées, bord de la carte qui réduisent énormément l’information sur
avec 75 mètres comme valeur minimale de la profondeur l’étendue des zones géothermiques, le calcul de la grille de
pour être incluses. Les autres critères sont la repré- traçage, intègre les données des zones frontalières : la mer
sentativité de l’emplacement du puits et l’absence de d’Alboran (Polyak et al. 1996), le sud espagnol (Albert
mouvements convectifs. Deux cartes sont présentées ; la Beltran 1979) et la marge atlantique (Hermann et al. 1977,
première (Fig. 1) donne les valeurs ponctuelles du gradient Villinger 1984). Les résidus entre les résultats de la grille
par forage ou par site ; la deuxième (Fig. 2) est une carte de calculée et les données originales sont représentés sur la
courbes d’isovaleurs de gradient géothermique qui montre figure 4.
les tendances macroscopiques de la distribution sur laquelle
sont superposées les sources chaudes. Pour établir la grille RESULTATS
de traçage de la carte nous avons utilisé le krigeage qui est
une méthode géostatistique d’interpolation bien adaptée aux Gradient géothermique
données irrégulièrement distribuées dans l’espace et qui La valeur du gradient thermique varie de 10°C/km à plus de
permet l’estimation des valeurs locales et de générer les 100°C/km (Figs. 1 et 2). La valeur moyenne, 35±11°C/km,
tendances générales en considérant l’organisation spatiale reste plus élevée que la moyenne globale qui est de l’ordre
des données. Les données spatiales sont caractérisées par le de 20°C/km à 25°C/km. Aussi, les histogrammes de la
comportement erratique typique des variables aléatoires et figure 5, montrent que l’intervalle 30-45°C/km est le plus
par la corrélation qui dépend de leurs coordonnées fréquent. Ceci pourrait être dû à ce que la plupart des
géographiques. L'étude de ces corrélations spatiales est déterminations ont été effectuées dans des zones anormale-

2
 Rimi – gradient géothermique du Maroc

30

20
géothermique et la valeur de la grille (°C/km)
Résidu entre la valeur observée du gradient

10

-10

Intervalle = 49.8, Moyenne = 0.57, erreur = 0.25

intervalle de confiance 95% = 0.48
-20 intervalle de confiance 99% = 0.64
Variance = 32.35, Ecart type = 5.7
Numbre de points utilises : 534

-30

0 100 200 300 400 500 600

numéro du point de mesure

Figure 4. Variation du résidu entre les valeurs déterminées du gradient géothermique et celles interpolées par le
krigeage pour établir les courbes d’égale valeur.

ment chaudes en liaison avec le volcanisme ou la tectonique sud-rifaines et le Prérif, montre également une structure
d’extension (Maroc oriental). Le long de la marge hydrogéologique complexe avec deux sens opposés de
atlantique, le gradient reporté d’après les données marines circulation d’eau (Ostapenko 1985) : les eaux d’infiltration
est élevé (40-60°C/km) (Hermann et al. 1977, Villinger à partir de la Meseta et des Rides sud rifaines en direction
1984) (Figs. 1 et 2) et pourrait être lié au rajeunissement de du nord-est et du nord ; les eaux d’expulsion coulant vers le
la croûte suite au volcanisme des Iles Canaries, dont les sud-ouest. Dans le Miocène du Gharb, le mouvement de
éruptions successives se sont effectuées en plusieurs cycles, l’eau est essentiellement conditionné par une migration
le plus ancien remontant à avant 35 Ma, alors que le plus verticale ascendante depuis les séries infra-nappes.
tardif date de 1 à 2 Ma (Staudigel & Schminke 1984). Les
générations des volcans des Canaries sont interprétées
comme étant la manifestation en surface de profonds
panaches mantelliques ascendants ou hot spots (cf. par ex.
Morgan 1971). Histogramme
On observe également une grande variation latérale au sein
120
d’une zone étroite telle que le Prérif (Fig. 2). Ceci est dû a
une combinaison de différents processus thermiques 100
complexes (Rimi et al. 1998). Les variations spatiales du
80
gradient géothermique ne coïncident pas parfaitement avec
Fréquence

celles du flux de chaleur ; en effet, si dans l’Oriental le 60

gradient et le flux géothermiques sont tous deux élevés, ce
n’est pas le cas dans le Prérif et la côte atlantique. La zone 40
qui s’étend le long du front de la Nappe prérifaine montre 20
une importante anomalie du gradient géothermique, de 40-
45°C/km, alors que le flux de chaleur demeure normal, 0
de 60 à 70 mW/m2 (Rimi 2001). Des variations rapides du
10
19
29
38
48
57
66
76
85
95
104
113

flux de chaleur sont observées aux limites séparant le bassin Gradient géothermique (°C/km)
du Gharb, la nappe du Prérif et les Rides sud-rifaines. De
telles perturbations peuvent être dues aux hétérogénéités
latérales de conductivité thermique (dômes de sel, Figure 5. Histogramme de la distribution des valeurs du gradient
proximité du socle granitique ou quartzitique et densité de géothermique au Marocl
fracturation). Cette région, qui couvre le Gharb, les Rides

3
 Rimi – gradient géothermique du Maroc

Tableau I. Distribution des propriétés des sources d'eau naturelles au Maroc

Tmin (°C) Tmax (°C) nombre pourcentage (%) ρmin (l/s) ρmax (l/s) nombre pourcentage (%)
< 2 7 17.9
20 30 21 37.5 2 10 11 28.2
30 40 21 37.5 10 50 14 35.9
40 50 11 19.6 50 250 5 12.8
> 50 3 5.4 > 250 2 5.1
RSmin RSmax nombre pourcentage (%) pHmin pHmax nombre pourcentage (%)
(g/l) (g/l)
< 0,6 11 23.4 < 6,5 8 19.5
0,6 1,2 12 25.5 6,5 7,5 27 65,9
1,2 2,5 8 17.0 > 7,5 6 14.6
2,5 5 10 21.3
> 5 6 12.8
Faciès chimique nombre pourcentage (%)
sulfuré 11 23.4
carbo-gazeux 12 25.5
chloruré 8 17.0
sulfaté 10 21.3
bicarbonaté 6 12.8

Du point de vue énergétique, la profondeur pour atteindre la région. Les plus chaudes que 30°C jalonnent le front de
une température de 60°C est de 1000 m dans l’Oriental et le la Nappe prérifaine et les frontières entre les unités
Prérif, alors qu’il faudrait descendre à 2000 m pour arriver structurales majeures et résurgissent préférentiellement
à 100°C. entre les courbes isogradients 30 et 40°[Link]-1 (Fig. 2). En
revanche, les sources chaudes dans l'est du Rif et tout le
Maroc oriental sont associées à l'élévation du flux
Sources thermales
mantellique. Dans le Maroc central, la source Lalla Haya (T
Les propriétés physico-chimiques de la source sont un = 42°C), même si du point de vue structural, elle est liée au
indicateur qui peut être utilisé, après comparaison avec des massif granitique important de la région, le granite
données géologiques, dans l'exploration hydro- d’Oulmes, son origine demeure toujours météorique. Le
géothermique. Les températures observées donnent les Haut Atlas oriental recèle des sources thermo-minérales
premières évaluations de la profondeur de l'origine de l'eau s'alimentant dans le Lias calcaire. La plus chaude (T>52°C
thermale. La majeure partie de ces sources chaudes est avec un débit de 5 l/s) est située à Foum Zaâbel au nord
concentrée dans le nord du pays : Maroc central, Prérif, d'Errachidia. C'est une eau très minéralisée chlorurée
Sillon sud-rifain et Maroc nord-oriental (Fig. 2). sodique, sulfatée calcique et magnésienne, bicarbonatée et à
forte teneur en fluor (Saqualli 1970). Enfin, les deux
Les informations recueillies dans les rapports publiés sur sources les plus chaudes dans les provinces du Sud, Abeino
les propriétés physico-chimiques des sources thermo-
(T=42°C) et Timoulay (T=40°C) sont situées dans la zone
minérales du Maroc (température, débit, résidu sec, pH et d’anomalie élevée de flux géothermique entre les Iles
composition chimique, …), ont été vérifiées ; ainsi, 36 Canaries et le bassin de Tindouf (Rimi 1999).
points ont été éliminés sur 131, soit presque 30%, pour des
raisons inhérentes à une qualité insignifiante ou mauvaise
de l’eau, à un débit faible, à une température faible, et à des
CONCLUSIONS
données doutuses ou incomplètes.
Bien que la méthode utilisée dans cette étude pour estimer
En moyenne (Tabl. I), les sources les plus chaudes sont
les gradients géothermiques comporte des hypothèses
celles du Maroc nord-oriental dont le contexte est
simplifiées et qu’il y a encore plusieurs zones non couvertes
clairement lié au volcanisme et à la tectonique récents.
par les mesures pour définir la structure des gradients
Dans le Rif occidental, les sources d'eau sont en général
géothermiques, il n’en demeure pas moins que la
connectées à des systèmes de circulation profonde, et leur
distribution des gradients est compatible avec le cadre
montée artésienne est favorisée par les fractures actives de

4
 Rimi – gradient géothermique du Maroc

géologique. Dans le nord-est du Maroc, le gradient élevé est occidentale et l’Anti Atlas. Le sous-sol marocain dispose
associé au volcanisme, alors que dans le Prérif, il ne peut réellement de potentialités géothermiques, et les courbes
être expliqué que par la tectonique cassante et le réseau température / profondeur, sont un diagnostic des processus
dense de fractures favorisant une circulation de fluides qui de transfert de chaleur en sub-surface. La carte générale du
engendre un processus de convection. La tendance gradient géothermique met en évidence les zones où doit
d’élévation du gradient sur la marge atlantique en direction avoir une couverture détaillée, en utilisant des sondages peu
des îles Canaries est associée a la présence dans la région profonds avec un espacement de 2-3 km, nécessaire pour
d’un hot spot. Par ailleurs, de faibles valeurs sont observées l’investigation des traits du champ thermique local.
dans les zones tectoniquement stables telles que la Meseta

Références

Albert Beltran J.F. 1979. El mapa español de flujas calorificas, Rimi A. 1990. Geothermal gradient and heat flow trends in
intento de correlacion entre anomalias geotermicas y Morocco. Geothermics, 19, 5, 443-454.
estructurales. Bol. Geol. Min., Madrid, XC-I, 36-48. Rimi A., Chalouan A. & Bahi L. 1998. Heat flow in the
Bahi L., El Yamine N. & Risler J.J. 1983. Linéaments westernmost part of the Alpine Mediterranean system (the
géothermiques au Maroc. C.R. Acad. Sci. Paris, 296, sér. II, Rif, Morocco). Tectonophysics, 285, 135-146
1087-1092. Rimi A. 1999. Mantle heat flow and geotherms for the main
Benabidate L. 1994. Contribution à l’étude hydrogéothermique du geologic domains in Morocco. Int. J. Earth Sci. (Geologische
Maroc nord-occidental (Rharb, Rides et Saïss). Thèse de Rundschau), 88, 458-466.
spécialité, Université de Sfax, Tunisie. Rimi A. 2001. Carte de densité de flux de chaleur et de gradient
Boukdir A. 1994. Contribution à l’étude géothermique du bassin géothermique dans le nord du Maroc. In : Symposium : La
du Tadla, Plateau des phosphates et Tassaout aval. chaîne rifaine dans son cadre méditerranéen occidental,
Application au réservoir calcaire du Turonien (Crétacé). Rabat 29-30 octobre, 2001.
D.E.S. Université Cadi Ayad, Fac. Sci. Marrakech, Maroc. Saqualli T. 1970. Les sources thermominérales du Prérif. Rapport
Cidu, R. & Bahaj S. 2000. Geochemistry of thermal waters from inédit n° 26, Ministère de l’Industrie et des Mines, SEGM,
Morocco. Geothermics, 29, 407-430. Meknès, Maroc. 101 p.
Hermann B.M., Langseth M.G. & Hobart M.A. 1977. Heat flow in Staudigel H. & Schminke H.U. 1984. The Pliocene seamount
the oceanic crust bounding West Africa. Tectonophysics, 41, series of La Palma/Canary Islands. J. Geophys., 96, 16083-
1-3, 61-78. 16100.
Lahrach A. 1994. Potentialités hydrogéothermiques du Maroc Villinger H. 1984. New heat flow values of the west of Morocco.
oriental. Doctorat de spécialité, Université de Sfax, Tunisie. Initial Reports DSDP, 79, 377-381.
Morgan W.J. 1971. Convection plumes in the lower mantle. Zarhloule, Y. 1999. Potentialités géothermiques du Maroc :
Nature, 230, 42-43. traitement intégré par les températures profondes et indices
Ostapenko S.V. 1985. Conditions hydrogéologiques et de surface. Doctorat d'Etat, Fac. Sci. Oujda, Maroc, 154 p.
perspectives pétrolières du bassin du Gharb-Prérif. Office Ziyadi R. 1993. Géologie appliquée à l’étude de l’environnement
National de la Recherche et d’Exploitation Pétrolière, Rabat, géothermique de la région de Nador (Rif nord- oriental,
Maroc, Rapport inédit, no 31389. Maroc). Doctorat de Spécialité, Université de Pau et des Pays
Polyak B.G., Fernandez M., Khutoskoy M.D., Soto J.I., Basov de l’Adour, France.
I.A., Comas M.C., Khain V. Ye, Alonso B., Agapova G.V.,
Mazurova I.S., Negredo A., Tochitsky V.O. , de la Lind J.,
Bogdanov N. A. & Banda E. 1996. Heat flow in the Alboran
Sea, western Mediterranean. Tectonophysics, 263, 191-218.
Rimi A. & Lucazeau F. 1987. Heat flow density measurements in Manuscrit reçu le 15 décembre 2001
northern Morocco. J. Afr. Earth Sci., 6, 835-843.