69 5
69 5
69 5
Méthodologie de
la cartographie des risques
naturels liés
aux mouvements de terrain
A. FARES
M. ROllET 1 (Région de Taounate .. Rif, Nord du Maroc)
~~
E
Cette étude présente une méthodologie de la
~BROQUET cartographie des risques naturels liés aux mouvements
~ de terrain à l'échelle du 1/25000, qui utilise un mode de
Laboratoire de géologie ~~ cartographie numérique basé sur une simulation du
structurale et appliquée ~ risque en fonction de l'importance et de la typologie des
1~ place Leclerc instabilités de versant.
25000 Besançon Les facteurs constituant les éléments du risque pour cette
cartographie sont relatifs à la topographie, à la lithologie
et à la géomorphologie. Ce choix des facteurs nous
permet de délimiter le cadre propice aux déclenchements
des mouvements de terrain. L'évaluation du risque
résultant de l'action conjointe des facteurs a été effectuée
suivant deux méthodes, la méthode arithmétique et la
méthode probabiliste, aboutissant à deux cartes
synthétiques dont une seule, fondée sur la méthode
probabiliste servira comme carte des risques naturels;
facile à lire, elle est donc utilisable par les différents
décideurs de l'aménagement du territoire.
1
Mots clés: Cartographie, géologie, simulation,
mouvement de terrain, risque, méthodologie,
quantification.
i;I::
~~JW
Ct.')..- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
1
~ :~i
Methodical test of natural risks
relative to ground movements
cartography
(Taounate area - Rif, North Marroco)
63
REVUE FRANÇAISE DE GÉOTECHNIQUE
W69
4e tnmestre 1994
et une fréquence variables dans le temps et dans
l'espace.
Introduction Nous avons fait le choix de cartographier à l'échelle
du 1/25 000 les facteurs permanents afin de délimiter
Le Rif forme un domaine d'étude riche et varié des les secteurs qui présentent le cadre favorable aux
instabilités de versants. Ces mouvements de terrain, au déclenchements des mouvements de terrain.
niveau du site de Taounate (Rif central), gênent l'amé-
nagement de la ville du même nom et constituent un 1
risque, d'où l'intérêt d'une cartographie à caractère
prévisionnelle pour définir les sites les plus exposés aux Établissement de la carte
instabilités des versants.
Historiquement, la cartographie des risques natu-
des risques naturels
rels a commencé pour la première fois au Maroc dans
les années soixante, en même temps que les cartes géo- :11
techniques de Casablanca, de Fès, de Tanger... En 1968,
A. Millies-Lacroix établit pour l'ensemble du Rif une Méthodologie (Fig. 1)
carte prévisionnelle des mouvements de terrain à une Pour aboutir à une carte des secteurs susceptibles
échelle de 1/1 000 000. Les facteurs utilisés dans cette d'être le siège des mouvements de terrain, nous avons
cartographie sont la lithologie au sens large, les hau- retenu les facteurs liés à la topographie, à la nature et la
teurs moyennes de précipitation annuelle, la topogra- structure du matériau et à la géomorphologie. Ces fac-
phie, le couvert végétal et l'action anthropique. La teurs sont cartographiés séparément pour obtenir les
région de Taounate apparaît dans des zones comprises cartes de facteurs du risque. La superposition de ces
entre les degrés 6 et 9, correspondant aux susceptibili- cartes fait ressortir les relations entre la lithologie, la
tés maximales de la carte prévisionnelle. pente et les indices d'instabilité. Une étude basée sur
En France, les études de cartographie des risques les calculs de rapports de surfaces nous permet de fon-
naturels ont commencé à partir du début des années der une simulation du risque que nous avons graduée
soixante-dix, par le Bureau de Recherche Géologiques en cinq niveaux dans une échelle numérique (1, 2, 3, 4
et Minières (BRGM) et le laboratoire de géologie appli- et 5). Ceci nous conduit à transcrire les cartes de fac-
quée de l'Université de Grenoble, dans une cartogra- teurs en cartes numériques appelées cartes de simula-
phie commune intitulée Plan Zermos cc zones exposées tion du risque. La carte de l'interprétation du risque
aux risques des mouvements du sol et du sous sol)) délimite des secteurs où figurent le niveau du risque
(W. Chazan, 1973). Les objectifs de ce programme n'ont assigné à chaque facteur selon des triplets.
pas été totalement atteints (P. Antoine, 1977), puisque L'évaluation du risque résultant de l'action conjointe
la conception du risque reste une notion subjective des facteurs a été effectuée suivant deux méthodes: la
interprétée différemment selon les auteurs des cartes méthode arithmétique et la méthode probabiliste abou-
Zermos réalisées jusqu'à maintenant. tissant à deux cartes synthétiques dont une seule ser-
Nous avons repris les préoccupations de la carto- vira comme carte des risques naturels à la suite de la
graphie Zermos et défini une méthodologie en utilisant comparaison des deux démarches.
le mode de cartographie numérique basée sur la simu-
lation des facteurs à partir des observations et des don- fi
nées de surface. Notre méthode s'appuie sur des
moyens faciles à mettre en œuvre et peu coûteux, à Les cartes de facteurs
savoir: une carte topographique à l'échelle du 1/25 000,
une carte géologique levée à l'échelle du 1/25 000, des
__ photographies aériennes et quelques essais géotech-
niques. A partir de ces données, nous avons élaboré
une série de cartes qui s'enchaînent dans un ordre Carte des pentes
logique pour aboutir à la carte des risques naturels liés
-
aux mouvements de terrain de la région de Taounate. Nous avons établi une carte des pentes à partir de la
carte topographique à l'échelle du 1/25 OOO.Elle pré-
sente six classes de pente de [0-2,5 0/0], [2,5-5 0/0], [5-
10 0/0], [10-20 0/0], [20-40 0/0] et [40-100 0/0].
Les facteurs déterminant
la genèse des mouvements
de terrain Carte des faciès lithologiques
La genèse d'un mouvement de terrain résulte de Elle a été établie à partir de la carte géologique à
l'action simultanée de plusieurs facteurs. Certains sont l'échelle du 1/50000 (G. Suter, 1967) et de levés nou-
permanents et forment le cadre propice créant les veaux. Elle présente une notice explicative qui classe
conditions indispensables aux mouvements, à savoir: les faciès suivant les caractéristiques géologiques,
la topographie, la lithologie et la géomorphologie; les hydrogéologiques et géotechniques. Nous y avons
autres dynamiques ont un rôle déclenchant et sont annexé également la carte de distribution de la fractu-
représentés par le climat, l'hydrologie, l'évolution de la ration dans les massifs gréseux qui présente en plus
dynamique actuelle des mouvements de terrain (A. Mil- des rosaces de fracturation, le pendage et le degré
lies-Lacroix, 1981), l'action de l'homme, la végétation et d'hétérogénéité en fonction de la dimension des blocs
N° 69
4e trimestre 1994
CARTES DES FACTEURS
Carte des pentes Cart~ lit~ologique . Carte géomorphologique
~ superposItIon et sImulatIon ~
~
..
Indice du risque (Il ------.
Ip (1, 2~ 3~ 4) IL (1, 2~ 3~ 4~ 5) IG(1, 2~ 3~ 4~ 5)
1 Cartes de simuTation du risque 1
, liées à chaque facteur (CSRF) ,
CSRL (IL) CSRG (IG)
- . superposition
..
Carte de l'interprétation du risque
~-
,
(IpILI G)
..
Evaluation du risque global
~
Évaluation arithmétique Évaluation probabiliste
sommation conversion de l'indice du risque
(Ip + IL + IG) en Probabilité (P)
..
I~P
(3 à 14)
•
Carte de la sommation du risque
..
P=PpXPLXP G
Carte des probabilités d~ apparition spatiale du risque
..
(1 à 1000/0)
Carte des classes du risque Carte des degrés du risque
[ (3, 4~ 5); (6, 7~ 8); (9~ 10~ 11); (12~ 13~ 14)] [(0 à 100/0); (10 à 30); (30 à 50); (50 à 1000/0)]
-------------. ..
Comparaison et choix
Carte des degrés du risque
~
•
Carte des risques naturels liés aux mouvements de terrain
Degré d'instabilité des classes des pentes des marnes bleues de la région de Taounate.
Instability degree of blues marIs slopes classes of Taounate area.
:1'
Carte des indices d'instabilité Les cartes analytiques
Elle a été établie à partir de l'étude des photogra-
phies aériennes de plusieurs missions espacées dans le
temps (1951, 1960, 1979 et 1986) complétée par les
observations sur le terrain concernant les dommages
occasionnés. En dehors des ravinements et des bad- Quantification du risque
lands qui sont des phénomènes actuellement actifs, la
cartographie des mouvements de masse a englobé sys- Dans la nature, chaque faciès lithologique présente
tématiquement les formes actives où le risque est cer- une pente limite en dessous de laquelle le versant
tain et les formes anciennes stabilisées ou en période n'engendre aucune instabilité. Nous avons procédé à
de rémission provisoire. une étude des surfaces basée sur la superposition sec-
65
REVUE FRANÇAISE DE GÉOTECHNIQUE
W69
4e trimestre 1994
teur par secteur de la lithologie, des pentes et des
indices d'instabilité; nous avons ainsi défini pour
chaque faciès les degrés d'instabilité en fonction des Les cartes de simulation du risque
classes de pentes (rapport des surfaces instables sur la
surface totale de la classe considérée). Cette confronta- A chacun des facteurs du risque, nous avons attri-
tion des cartes de facteurs nous permet de situer bué un indice (If) dans une échelle numérique en rap-
approximativement la pente limite de chaque faciès: port avec son importance relative dans la genèse des
instabilités. Cette simulation nous permet de traduire
Les marnes bleues ne présentent pas de pente supé-
les cartes de facteurs en cartes de simulation du risque.
rieure à 40 % (Tab. 1). A partir des pourcentages des
surfaces instables, on peut caler approximativement la • Carte de simulation du risque lié à la pente
pente limite vers 10 % et plus précisément entre 5 et La détermination cartographique de la pente limite
10%. naturelle pour chaque matériau rend possible l'utilisa-
Le pourcentage des surfaces instables élevé de la tion directe de la carte des pentes dans une cartogra-
classe [10, 20] traduit l'héritage d'anciennes coulées et phie des risques naturels. Nous pouvons comparer les
correspond à des terrains relativement altérés (Tab. II). pentes des versants (Pv) aux pentes limites (Pl) des dif-
La pente limite des marnes schisteuses ne pourra avoir férents faciès et attribuer un indice de risque (I ),
une valeur supérieure à 20 % que si la roche est saine. comme suit: P
Contrairement aux marnes bleues et aux marnes - si Pv est inférieur à Pl, la pente du versant est stable
schisteuses, les grès n'ont pas de pente limite (Tab. III). et l prend la valeur minimum de : 1;
Les instabilités correspondant aux différentes classes - stpv est sensiblement égal à Pl, l est égal à : 2;
dépendent d'autres facteurs tels que la structure et - si Pv est supérieur à Pl, lest égaf à : 3;
l'altération. - si Pv est très supérieur à Pl, l prend la valeur maxi-
Il ressort de cette confrontation, que toutes les mum de: 4. P
pentes inférieures à 10 % sont relativement stables et Le tableau IV montre les différents indices attribués
qu'à partir de 10 % la stabilité des terrains dépend de la en fonction des degrés d'instabilité des classes de
nature lithologique, de la structure ou de l'altération. pente:
Degré d'instabilité des classes des pentes des marnes schisteuses du Sénonien de la région de Taounate.
Instability degree of Senonian schistous marIs slopes classes of Taounate area.
5à 100/0 0,433 0 0
10 à 200/0 3,733 t4124 37,83 Coulées anciennes
20 à 400/0 3,2715 1,3185 40,3
>à400/0 0,5308 0,2455 46,25 Niches de départ et ravinement
Degré d'instabilité des classes des pentes des grès du Tortonien supérieur de la région de Taounate.
Instability degree of Upper Tortonian sandstones slopes classes of Taounate area.
5 à 100/0 0,2314 0 0
10 à 200/0 2,5422 0,4453 17,5 Glissement de terrain
20 à 400/0 7,9492 0,4242 5,4
> à 400/0 1,6726 0,1403 8,38 Glissement et éboulement
20 à 400/0 (-) 4 3 3 3 1 3
40 à 1000/0 (-) (-) 4 4 4 2 (-)
66
REVUE FRANÇAISE DE GÉOTECHNIQUE
N° 69
4e tnmestre 1994
• Carte de simulation du risque lié à la lithologie -la pente est favorable à l'instabilité: 3;
Nous pouvons faire une première approximation en - le matériau est très sensible aux mouvements de ter-
comparant le degré d'instabilité des différents faciès rain: 4;
lithologiques (Tab. V). - les risques manifestés ou à craindre sont les glisse-
Compte tenu du degré et du type d'instabilité, de la ments de terrain et les éboulements: 5.
nature et de l'épaisseur de la couverture d'altération, et
de la structure; nous avons traduit la carte des faciès
lithologiques en une carte de simulation du risque lié à I:évaluation du risque globa~
la lithologie, basée sur les indices du risque attribués
aux faciès suivants: et les cartes synthétiques
- les molasses: 5; La transformation des cartes de facteurs en cartes
-les marnes bleues saines et altérées, le complexe bré- numériques nous permet d'évaluer le risque naturel
chique gypsifère, les marnes schisteuses altérées, les suivant deux méthodes: la méthode arithmétique et la
grès très fracturés à pendage conforme à la pente, les méthode probabiliste.
grès à couverture d'altération issue d'un mouvement
de masse: 4;
- les marnes schisteuses, les grès à pendage conforme
à la pente: 3; Évaluation arithmétique
- les grès subaffleurants, les schistes gréseux: 2;
- les terrasses anciennes: 1. Cette méthode a été utilisée par A. Millies-Lacroix
• Carte de simulation du risque lié à la géomorphologie (1968). L'action simultanée des trois facteurs de la carte
L'évaluation du risque peut être faite en fonction de l'interprétation du risque est représentée par la
des masses déplacées, nous avons donc délimité les somme des différents indices dans chaque secteur.
secteurs où domine chaque type d'instabilité en leur L'exemple du triplet l I l précédent c'est-à-dire (3, 4, 5)
fournit une somme d~ 12~
J
Degré d'instabilité des terrains de la région de Taounate en fonction de leur nature lithologique.
Instability degree of Taounate area formations aeeording to their lithologie nature.
67
REVUE FRANÇAISE DE GÉOTECHNIQUE
N°69
4e trimestre 1994
La carte de l'interprétation du risque.
The map of the risk interpretation.
68
REVUE FRANÇAISE DE GÉOTECHNIQUE
N° 69
4e trimestre 1994
o lKm
1
69
REVUE FRANÇAISE DE GÉOTECHN~?~~
4e trimestre 1994
La carte des degrés du risque au contraire voit cette
imprécision levée, car elle exprime le risque par une
Évaluation probabiliste probabilité d'apparition spatiale, et par conséquent la
lecture de la carte directe est facilitée, donc accessible à
Nous avons gradué le risque de 0 à 1 et converti un large public.
les indices du risque en terme de probabilité partielle
(Tab. VI). • Rapports de la carte des risques avec la carte géolo-
gique
En conséquence, l'évaluation du risque est estimée
égale au produit des probabilités des différents facteurs. La carte des classes du risque et la carte des degrés
du risque présentent les mêmes tendances dans les
PR (0/0) == Pp x Pl x Pg x 100 secteurs à haut risque et des recoupements quand les
PR est la probabilité d'apparition spatiale du risque. secteurs sont classés avec un risque douteux. Cepen-
Lorsque tous les facteurs sont à leur maximum dant les limites de la carte des classes du risque suivent
d'intensité (PR) est de 100 0/0, valeur pour laquelle la en quelque sorte les contours géologiques, ce qui limite
conception du risque est certaine sinon évidente: l'intérêt de la démarche arithmétique à l'échelle locale.
PR (0/0) == 1 x 1 x 1 x 100 == 100 0/0 Cette concordance fait défaut quand on superpose la
carte des degrés du risque avec la carte géologique. Il
Par contre, si les facteurs sont à leur minimum paraît possible alors d'adapter cette approche probabi-
d'intensité, la probabilité est réduite à 1 0/0. liste de la quantification du risque à d'autres régions
PR (0/0) == 0,25 x 0,2 x 0,2 x 100 == 1 0/0 présentant des problèmes similaires et ayant des carac-
• Carte des probabilités d'apparition spatiale du risque téristiques géologiques différentes, et c'est là qui réside
Les facteurs du risque considérés ne varient que tout l'intérêt de l'approche probabiliste par rapport à
dans l'espace et la probabilité d'apparition spatiale l'approche arithmétique.
écarte donc toute prévision dans le temps. • Méthode choisie pour l'établissement de la carte des
A partir de la carte de l'interprétation du risque, risques naturels
nous pouvons établir directement la carte de la proba- La comparaison des deux méthodes nous a conduit
bilité d'apparition spatiale du risque qui présente à abandonner la méthode arithmétique et à retenir la
23 secteurs avec des probabilités allant de 1 à 100 0/0. méthode probabiliste, car elle met en jeu l'interaction
• Carte des degrés du risque (Fig. 4) des facteurs et donne une idée relativement précise sur
Nous avons ramené les pourcentages de probabilité leur action conjointe par la probabilité d'apparition
de la carte précédente à quatre niveaux du risque crois- spatiale du risque. Nous utiliserons donc la carte des
sant: degré du risque faible ou négligeable (1 à 100/0), degrés du risque comme composante fondamentale
degré du risque moyen (10 à 30 0/0), degré du risque pour la carte des risques naturels.
élevé (30 à 50 %) et degré du risque très élevé (50 à
1000/0).
lp 1 2 3 4 Il 1 2 3 4 5 19 1 2 3 4 5
70
Pp 0,25 0,5 0,75 1 P1 0,2 0,4 0,6 0,8 1 Pg 0,2 0,4 0,6 0,8 1
~~- ./
~ t
~; \ ~~
~ "'" 'l
~ J..........'"
u.u..L
~lKm
1
1
1 1
71
REVUE FRANÇAISE DE GÉOTECHN~C?~~
Antoine P. (1977) - Réflexions sur la carto- aux mouvements de terrain; application visionnelle des mouvements de terrain
graphie Zermos et bilan des expériences à l'aménagement de la ville de Taounate dans le Rif. Mines et géologie, 27, pp. 45-
en cours. Bull. BRGM, III, 1-2, pp. 9-20. (Rif, Maroc). Thèse Univ. Besançon, 54.
Chazan W. (1973) - Le plan Zermos. Identi- 177p.
fication des zones exposées aux risques Millies-Lacroix A. (1981) - Classification
Humbert M , Jeannette A. (1962) - Les
liés aux mouvements du sol et du sous- des talus et versants instables. Bull. Lia-
cartes géotechniques: quelques prin-
sol. Symposium national Sol et Sous-sol son Labo. Ponts et Chaussées, Spécial X.
et sécurité des constructions, Cannes cipes de réalisation. R. Géogr. Maroc,
BRGM,1. n° 1-2, pp. 67-74.
Suter G. (1967) - Carte géologique de
Farès A. (1994) - Essai méthodologique de Millies-Lacroix A. (1968) - Les glissements Taounate-Aïn Aïcha à 1/50 000. Not. et
la cartographie des risques naturels liés de terrain, présentation d'une carte pré- Mém. du Serv. GéoI. Maroc, 186 bis.
72 * L'ensemble des références bibliographiques peut être trouvé dans la thèse de A. Fares, 1994.