Caractérisation géologique et

géomécanique des massifs rocheux

pour l’abattage à l’explosif

Jean-Alain Fleurisson (CESECO-Géosciences)

 Le massif rocheux est le seul paramètre qui ne
peut être changé, et toutes les opérations minières
doivent être adaptées à lui

Méthode d’exploitation et équipement

Forage

Fragmentation à l’explosif

Stabilité des pentes et des cavités souterraines

Concassage et broyage
 CARACTERISTIQUES DES MASSIFS ROCHEUX

•MILIEUX DISCONTINUS: à différentes échelles

Discontinuité : toute surface qui interrompt physiquement la continuité

de la matrice rocheuse et de ses propriétés mécaniques, hydrauliques et
thermiques

mm Matrice rocheuse Joints de grain

Microfissures

cm Échantillon Fissures, schistosité, foliation

dm Bloc Diaclases

m Massif rocheux Contacts lithologiques

Zone fracturée

dam Massif rocheux Joints de stratification

Failles
 CARACTERISTIQUES DES MASSIFS ROCHEUX

• MILIEUX HETEROGENES

– Hétérogénéité de formation
Massifs sédimentaires Alternance de bancs de natures différentes
Variations latérales de faciès et d’épaisseur

Massifs granitiques Plusieurs unités pétrographiques traversées de

filons
– Hétérogénéité d’altération

Massifs plutoniques et Profils d’altération météoriques très irréguliers en

métamorphiques fonction du degré de fracturation du massif et de
l’altérabilité de la matrice
 CARACTERISTIQUES DES MASSIFS ROCHEUX

• MILIEUX ANISOTROPES

– Anisotropie de formation

Massifs Alternance de roches aux lithofaciès très contrastés

sédimentaires Joints de stratifications aux faibles caractéristiques
mécaniques

– Anisotropie de fracturation

• MILIEUX A FAIBLE POROSITE


MATRICE ROCHEUSE

RESEAU DE DISCONTINUITES
 Principaux paramètres de la matrice

Paramètres pétrophysiques

Nature: Mineralogie et texture
 Choix de la méthode de forage

Propriétés sismiques & Micro fracturation
 Type d’explosif (impédance)
 Propagation des vibrations
 Efficacité du forage
 Propagation des micro-fissures
 Vitesse des ondes ultrasonores (Vp & Vs)

V 
IC (%) = 100 ∗  p
• Indice de continuité V* 
 p
Vp: vitesse mesurée sur échhntillon
Vp*: vitesse théorique calculée à partir de la composition minéralogique de
la matrice
CLASSE VALEURS DE DESCRIPTION
INDICE DE CONTINUITE IC DE LA ROCHE
IC1 IC > 90% Très grande continuité
IC2 75% < IC < 90% Grande continuité
IC3 50% < IC < 75% Continuité moyenne
IC4 25% < IC < 50% Continuité faible
IC5 IC < 25% Très faible continuité
Principaux paramètres de la matrice rocheuse

Paramètres pétrophysiques

Minéralogie et texture
 Choix de la méthode de forage

Micro fracturation (Indice de Continuité IC)
 Efficacité du forage
 Propagation des microfractures

Altération
 Usure des équipements

Masse volumique (1.5 to 3 g/cm3)
 Type d’explosifs (impedance)

Propriétés abrasives
 Efficacité du forage
 Usure des équipements
Test d’abrasivité du CERCHAR : test d’usure d’un foret
Norme française P94-430-1

CLASSE VALEURS DE DESCRIPTION

L’INDICE D’ABRASIVITE CERCHAR
CAI 1 > 4.0 Roche extrêmement abrasive
CAI 2 2.0 – 4.0 Roche très abrasive
CAI 3 1.0 – 2.0 Roche abrasive
 Test d’abrasivité LCPC:
test d’usure d’une palette en rotation
(Norme française P94-430-2)

CLASSE VALEUR D’ABRASIVITE ABR DESCRIPTION

ABR1 ABR > 2000 Roche très abrasive
ABR2 1500 < ABR < 2000 Roche abrasive
ABR3 1000 < ABR < 1500 Roche moyennement abrasive
ABR4 500 < ABR < 1000 Roche peu abrasive
ABR5 0 < ABR < 500 Roche très peu abrasive
 Test d’abrasivité Norvégien (DRI)
Trondheim university in Norway

S20 Test de Fragmentation SJ test d’usure d’un foret en rotation

CLASSE VALEUR D’ABRASIVITE ABR DESCRIPTION
ABR1 ABR > 2000 Roche très abrasive
ABR2 1500 < ABR < 2000 Roche abrasive
ABR3 1000 < ABR < 1500 Roche moyennement abrasive
ABR4 500 < ABR < 1000 Roche peu abrasive
ABR5 0 < ABR < 500 Roche très peu abrasive
Principaux paramètres de la matrice rocheuse

Paramètres Pétrophysique

Minéralogie et texture

Propriétés sismiques et microfracturation

Altération

Masse volumique

Propriétés abrasives

Paramètres mécaniques

Paramètres élastiques (module de Young & coefficient
de Poisson)

Résistance à la compression simple

Résistance à la traction brésilien
 Déformation et rupture de la roche
 Résistance à la compression simple
(Norme française P 94-420 and P 94-425)

Résistance à la compression simple

Vitesse de déformation 120 micrometer/minute

25

Contrainte axiale MPa

20

15 y = 1959,4x - 4,3221

10

0
4
Déformation radiale * 10

0 0,002 0,004 0,006 0,008 0,01 0,012 0,014 0,016 0,018

-5
y = -0.199x + 0.000003
-10
Module de Young = 1959,4 MPa
-15
Coefficient de Poisson = 0,2
Résistance à la compression simple = 21,75 MPa
-20
Déformation longitudinale

Contrainte - Déformation globale Contrainte-Déformation axiale locale Déformation radiale - Déformation axiale
Classes de Résistance à la Compression simple (RCS)
(AFTES, 2003)

CLASSE VALEURS DE RESISTANCE A DESCRIPTION DE LA

LA COMPRESSION SIMPLE ROCHE

RC 1 σc > 200 MPa Extrêment résistant

RC 2 100 MPa < σc < 200 MPa Très résistant

RC 3 50 MPa < σc < 100 MPa Résistant

RC 4 25 MPa < σc < 50 MPa Moyennement résistant
RC 5 5 MPa < σc < 25 MPa Tendre
RC 6 1 MPa < σc < 5 MPa Très tendre
RC 7 σc < 1 MPa Extrêment tendre
Effet de la taille et de l’anisotropie sur le RCS
 Essai de traction brésilien
(Norme française P94-422)

F 2F
F σt =
πDL
6F
σy =
σy πDL

τ Cercle de Mohr

σt σt

σy
σt 0 3σt σ
 Essai de compression entre pointe
(Essai Franklin)
(Norme française P94-429)

F
IS =
De2

F = Charge à la rupture
De = Distance entre pointes

De = D pour un essai selon le diamètre

= 4A/π pour un essai sur un bloc
20 Is50 < RCS < 27 Is50 A = surface minimale de la section plane
passant par les pointes
 Masse volumique, RCS, Résistance à la traction et
modules d’élasticité de quelques roches courantes

ROCHE Masse RCS Résistance Module de Coefficient

volumique (Mpa) à la traction Young de Poisson
(103 kg.m-3) (Mpa)
Granite 2.6 – 2.75 100 - 230 13 - 14 10 - 70 0.1 – 0.25
Basalte 2.6 - 2.8 80 - 120 9 30 - 80 0.15 – 0.25
Calcaire 2.2 – 2.7 40 - 120 5 - 10 10 - 80 0.25 – 0.30
Dolomie 2.4 – 2.8 100 10 10 - 60 -
Grès 2.1 – 2.8 30 - 90 1 - 15 10 - 60 0.20 – 0.30
Schistes 2.4 – 2.8 10 - 150 10 10 - 30 0.10 – 0.20

MATRICE ROCHEUSE

RESEAU DE DISCONTINUITES
 CARACTÉRISATION
DU RÉSEAU DE DISCONTUITÉS

Étude à l'échelle de l'ouvrage considéré

Méthodologie bien spécifique

Analyses statistiques en rapport avec

la variabilité naturelle des paramètres
géométriques et mécaniques
 TYPES et ORIGINE DES DISCONTINUITES
Discontinuité : toute surface qui interrompt physiquement la continuité de la matrice
rocheuse et de ses propriétés mécaniques, hydrauliques et thermiques

Types de Paramètres géologiques et géométriques

discontinuité
Joints de grains Faible cohésion
Microfissures Surfaces de discontinuités possibles
Diaclases Épontes planes, peu rugueuses - Extension dm à dam
Réseaux de 2 ou 3 familles
Couloirs de fracturation Assemblage complexe de discontinuités
Possibilité de structure planaire
Plans de schistosité Extension variable
Surfaces de très grande extension possibles par rupture localisée
de ponts rocheux
Contacts lithologiques Phénomènes d’altération différentielle

Plans de stratification Très grande extension - Ondulation faible

Présence de matériaux argileux
Failles Extension variable (m à km)
Matériaux de remplissage avec faibles caractéristiques mécaniques
Diaclase plane aux épontes légèrement altérées et enduits de
ferro-manganèse dans un granite
Diaclase avec matériau de remplissage silteux dans un granite
Ondulation d’une éponte de diaclase dans un granite
Discontinuité avec ponts rocheux et mécanisme de rupture plane
Plans de schistosité
Plans de schistosité
Joints de stratification et rupture en dièdre
 Joints de stratification dans un grès
(effet arrière sur un mur de prédécoupage)
 ANALYSE DU RÉSEAU DE DISCONTINUITÉS
1. MESURES DE DISCONTINUITÉS SUR LE TERRAIN

Définition de zones homogènes

Levés systématiques de discontinuités

Mesure des principaux paramètres (ISRM)

Orientation (direction ou azimut, pendage)

Paramètres Connectivité du réseau
Espacement
géométriques Anisotropie du massif
Extension ou longueur de trace visible

Ouverture des discontinuités

Nature du matériau de remplissage
Paramètres Altération des épontes Résistance au cisaillement
mécaniques Dilatance
Ondulation et rugosité
Présence d’eau
 METHODE DE LA LIGNE DE MESURES

• Orientation (direction ou azimuth, pendage)

• Espacement
• Extension ou la longueur de trace visible
• Ouverture des discontinuités, nature et épaisseur du matériau de
remplissage
• Ondulation et rugosité (JRC: Joint Roughness Coefficient)
• Altération des épontes (JCS: Joint Compressive Strength)
• Présence d’eau
 Mesure sur le terrain de
Profils de rugosité de discontinuité
Scléromètre (marteau de Schmidt) pour la mesure sur le
terrain de la résistance des épontes JCS
2. CLASSIFICATION EN FAMILLES DIRECTIONNELLES

Projection Stéréographique (Rosace des azimuts du vecteur
pendage)
2. CLASSIFICATION EN FAMILLES DIRECTIONNELLES

Projection Stéréographique (Isoligne de densité de pôles)
 2. CLASSIFICATION EN FAMILLES DIRECTIONNELLES

Indice de distance
r r
Dij = ni . nj

Regroupement des familles les plus proches k (Nk) et l (Nl) (valeur

maximale de Dkl) en une nouvelle famille m (Nm = Nk + Nl )
– Calcul du vecteur “normale indicative”
r r
r' N n + Nl nl
nm = ( nmx
' '
, nmy '
, nmz ) = k rk r
N k nk + N l nl
– Calcul du vecteur normal réel (somme des vecteurs normaux orientés ou
réorientés dans la direction de la normale indicative)
Nm
r r r
∑ ω . ni ni . nm'
r
nm = ( nmx , nmy , nmz ) = i =1 ω = r r ' = ±1
Nm
r ni . nm
∑ ω . ni
i =1

Calcul de l’azimut et du pendage du vecteur pendage de la famille m

2. CLASSIFICATION EN FAMILLES DIRECTIONNELLES

azimut pendage
Nombre d’éléments
35,0 60,0 35,0 60,0 35,0 60,0
1 1 1 33,5 61,5 33,5 61,5 33,5 61,5
32,0 63,0 32,0 63,0 32,0 63,0 2 2 2 35,7 62,6
1 1 1 3
40,0 65,0 40,0 65,0 40,0 65,0 40,0 65,0 40,0 65,0 40,0 65,0
1 1 1 1 1 1 63,4 84,4
272,0 73,0 272,0 73,0 6
1 1 271,0 71,5 271,0 71,5 271,0 71,5
270,0 70,0 270,0 70,0 2 2 2 269,0 72,6 269,0 72,6
1 1 3 3
265,0 75,0 265,0 75,0 265,0 75,0 265,0 75,0 265,0 75,0 81,8 59,1
1 1 1 1 1 9
155,0 32,0
1 157,4 31,0 157,4 31,0 157,4 31,0
160,0 30,0 2 2 2 159,1 32,3 159,1 32,3 159,1 32,3 159,1 32,3
1 3 3 3 3
162,0 35,0 162,0 35,0 162,0 35,0 162,0 35,0
1 1 1 1
2. CLASSIFICATION EN FAMILLES DIRECTIONNELLES
Famille Nombre Vecteur pendage moyen Azimut Pendage
éléments Azimut Pendage σθ µ σ µ σ
7 1 3 103.1 43.5 16.3 100.0 20.8 45.3 8.2
2 10 41.3 80.0 13.7 41.4 9.4 80.1 5.3
3 107 333.3 81.7 18.0 333.3 15.2 82.0 6.8
4 18 72.6 68.0 12.0 72.3 9.7 68.2 7.9
5 15 252.7 68.6 14.4 252.8 6.4 68.7 11.9
6 4 152.5 26.0 12.1 157.8 13.6 26.5 10.8
7 22 292.1 72.7 12.8 292.1 9.6 73.0 9.0
6 1 7 125.8 30.9 21.1 134.0 33.3 34.6 13.5
2 10 41.3 80.0 13.7 41.4 9.4 80.1 5.3
3 107 333.3 81.7 18.0 333.3 15.2 82.0 6.8
4 18 72.6 68.0 12.0 72.3 9.7 68.2 7.9
5 15 252.7 68.6 14.4 252.8 6.4 68.7 11.9
6 22 292.1 72.7 12.8 292.1 9.6 73.0 9.0
5 1 28 61.1 71.6 20.0 61.4 17.7 72.5 7.8
2 7 125.8 30.9 21.1 134.0 33.3 34.6 13.5
3 107 333.3 81.7 18.0 333.3 15.2 82.0 6.8
4 22 292.1 72.7 12.8 292.1 9.6 73.0 9.0
5 15 252.7 68.6 14.4 252.8 6.4 68.7 11.9
4 1 37 276.6 70.0 22.6 276.3 21.1 71.3 10.5
2 7 125.8 30.9 21.1 134.0 33.3 34.6 13.5
3 107 333.3 81.7 18.0 333.3 15.2 82.0 6.8
4 28 61.1 71.6 20.0 61.4 17.7 72.5 7.8
3 1 65 261.2 86.4 33.1 261.1 26.2 86.8 18.4
2 7 125.8 30.9 21.1 134.0 33.3 34.6 13.5
3 107 333.3 81.7 18.0 333.3 15.2 82.0 6.8
2 1 114 332.4 85.0 24.7 332.3 17.6 85.2 15.5
2 65 261.2 86.4 33.1 261.1 26.2 86.8 18.4
1 1 179 311.9 84.5 43.0 312.9 39.5 85.3 16.2
3. ANALYSE STATISTIQUE DES PARAMÉTRES GÉOMÉTRIQUES
 ANALYSE DE L’ORIENTATION
• DIRECTION ET PENDAGE

Variables non corrélées

Loi normale

• VECTEUR PENDAGE

Loi normale hémisphérique (loi de Fisher)
P(θ): probabilité que l’angle entre le vecteur pendage et le
vecteur pendage moyen de la famille est inférieur à θ

P (θ ) =
[1 − exp[ K (cosθ − 1)]]
[1 − exp(− K )]

K = Indice de concentration
 ANALYSE DE L’EXTENSION
• LONGUEURS DE TRACE A 2D

Loi exponentielle négative ou log-normale

Biais géométrique des mesures
Effet de censure
Effet de troncature

• PASSAGE DES LONGUEURS DE TRACE À 2D À

L’EXTENSION RÉELLE À 3D
Nécessite des hypothèses sur la forme des discontinuités ( par ex. disques)

Distribution des Quelconque Log normale Exponentielle

rayons de disques

E (l) 16 E(r 3 ) 16 5 16 λ
exp(µ + σ 2 )
3π E(r 2 ) 3π 2 π
E(l2) 3 E(r 4 )
3 exp( 2µ + 6σ 2 ) 36 λ2
E(r 2 )
 ANALYSE DE L’ESPACEMENT

• LOIS DE DISTRIBUTION

Loi normale ou log-normale

Cas des joints stratigraphiques des massifs sédimentaires, des
joints de refroidissement des colonnes basaltiques

Distribution exponentielle négative

Cas des massifs rocheux ayant subi plusieurs phases tectoniques
(espacement observé résulte de la superposition des distributions
successives avec réduction en taille des plus grands
espacements)
RELATION DENSITÉ - ESPACEMENT - EXTENSION
Densité volumique λv Nombre moyen de discontinuités par unité de
volume
Densité surfacique λs Nombre moyen de traces de discontinuité par unité
de surface
Densité linéaire λl Nombre moyen d’intersection entre discontinuité et
λl = 1/s ligne de mesure (s = espacement)

Hypothèses Position du centre des disques suit une loi de

Poisson de paramètre λv
( λ v V )n
P ( N = n) = [exp(− λ v V )]
Discontinuité = disque n!
Espacement suit une loi de distribution exponentielle
négative 1
P( S = s) = [exp(−s / λl )]
λl

Θ = angle entre la ligne de mesure et le vecteur

λ v = λ l /[π . cos θ . E ( r )]
2
normal moyen de la famille
α = angle entre la normale à la surface et le vecteur
λ v = λ a /[2 . sin α . E ( r )]
normal moyen de la famille
 MODÉLISATION GÉOMÉTRIQUE À 3D DES MASSIFS ROCHEUX
Modèle booléen de disques aléatoires - Logiciel SIMBLOC

SIMULATION DES DISCONTINUITÉS

Discontinuité = disque
Centre et rayon (extension)
Direction et pendage

Nombre de discontinuités de chaque famille

Processus de Poisson (densité volumique)

Génération des centres de disques

Loi uniforme
Volume de simulation + marge tampon

Tirage au sort des paramètres géométriques

Orientation et extension
 MODÉLISATION GÉOMÉTRIQUE À 3D DES MASSIFS ROCHEUX
Logiciel Simbloc
ÉTUDE DE LA CONNECTIVITÉ
Identification des intersections

Elimination des discontinuités dont le
degré de connectivité est inférieur à 3

1m

SIMULATED CONNECTED
DISCONTINUITIES DISCONTINUITIES
MODELISATION GEOMETRIQUE A 3D DES MASSIFS ROCHEUX
Modèle booléen de disques aléatoires - Logiciel SIMBLOC

IDENTIFICATION DES BLOCS

Sommets Intersection de 3 traces

Arêtes Segment entre 2 sommets

Facettes Surface fermée dont la

frontière est constituée
d’arêtes orientées

Blocs Critères topologiques

 BLOCOMETRIE IN SITU et
ANISOTROPIE DES BLOCS
F REQUENCY
3D BLOCK SIZE DISTRIBUTION

AV ERAGE SIZE 0.36 m

X = G ( 35% ) 1.14 m

X= G ( 50% ) 1.38 m

X= G ( 65% ) 1.59 m

BLOCK SIZE (m) STEREOPLOT OF M AXIM UM BLOCK AXIS

VOLUME / TOTAL VOLUME

BLOCK SIZE (m)

Logiciel SIMBLOC (CGI)

 INDICES GLOBAUX DE DENSITE DE
FRACTURATION

• ROCK QUALITY DESIGNATION (RQD)

RQD =
∑ (longueur de carottes ≥ 10 cm)
Longueur totale de la passe carottée

• INTERVALLE ENTRE DISCONTINUITÉ ID

Moyenne des intervalles découpés par

des discontinuités successives le long
d’une ligne de mesure
DISTRIBUTION DES LONGUEURS DE CAROTTES

Histogramme et courbe cumulée des longueurs de carottes

(d’après C. Louis, 1974)
 Influence de la direction du sondage
sur la caractérisation des discontinuités
CARACTÉRISATION GLOBALE DU MASSIF ROCHEUX

1. MÉTHODES GÉOPHYSIQUES DE SURFACE

Relation vitesse sismique - défonçabilité
CARACTÉRISATION GLOBALE DU MASSIF ROCHEUX

1. MÉTHODES GÉOPHYSIQUES DE SURFACE

CARACTÉRISATION GLOBALE DU MASSIF ROCHEUX

1. MÉTHODES GÉOPHYSIQUES DE SURFACE

CARACTÉRISATION GLOBALE DU MASSIF ROCHEUX

2. DIAGRAPHIES INSTANTANÉES
CARACTÉRISATION GLOBALE DU MASSIF ROCHEUX

2. DIAGRAPHIES DIFFÉRÉES

Méthodes Paramétres Unité Conditions

mesurée d’utilisation
Electriques Résistivité Ohm.m Trou non tubé rempli
d’eau ou de boue
Soniques ou Temps ou vitesse de ms ou m/s Trou non tubé vide ou
microsismiques propagation rempli d’eau
Radioactivité Nombre d’impulsion Trou vide ou non, non
naturelle pendant un temps tubé de préférence
donné
Neutron - neutron Teneur en eau Trou vide non tubé
volumique
Radar Temps ou vitesse de Nano sec Trou tubé plastique ou
propagation Km/s non, vide ou rempli
d’eau
Diagraphie microsismique
Diagraphie microsismique et de radioactivité naturelle
Tomographie sismique