GLQ3401/GLQ3651 : Première partie

Cours 5 : Estimation des ressources

par méthodes conventionnelles

1
 Objectifs
• Comprendre le fonctionnement des méthodes conventionnelles
d'estimation encore en usage en mines : polygone, triangle, inverse de
la distance;

• Être capable d’utiliser les méthodes conventionnelles d’estimation sur

des cas simples;

• Expliquer les forces et les faiblesses de chaque méthode d'estimation;

• Comprendre le concept de validation croisée et pouvoir l'appliquer;

• Calculer la densité théorique d'une roche à partir de ses principaux

minéraux et d'une analyse chimique.

GLQ3401: Géostatistique et géologie minières 2

GLQ3651: Géologie minière
 Plan du cours
1. Définitions des ressources et norme canadienne NI43-101
• Définition
• Rapport technique de ressources (NI43-101)
2. Estimation et influence de la densité
3. Notions de statistique
4. Méthodes conventionnelles d’estimation des ressources
• Polygone
• Triangle
• Inverse de la distance
5. Critique des méthodes conventionnelles
6. Validation croisée
7. Méthode des sections
8. Calcul de la densité théorique

GLQ3401: Géostatistique et géologie minières 3

GLQ3651: Géologie minière
 1. Définitions des ressources et norme canadienne NI43-101
Estimation des ressources

Définition : Détermination et définition, à partir d’un modèle de blocs, du

tonnage et de la teneur en minerais d’un dépôt géologique.
• Basé sur une interprétation géologique;
• Nécessite un traitement statistique des données;
• Écrêtage (capping)
• Normalisation (compositing)
• Interpolation des résultats;
• Géostatistique;
➢Méthodes conventionnelles (p. ex., polygone, triangle, inverse de la distance);
➢Krigeage et cokrigeage;
➢Krigeage d’indicatrices (Multiple indicator kriging– MIK)
➢Conditionnement uniforme (Uniform conditioning – UC)
➢Simulations géostatistiques

GLQ3401: Géostatistique et géologie minières 4

GLQ3651: Géologie minière
 1. Définitions des ressources et norme canadienne NI43-101
Exemple d’un modèle de blocs

GLQ3401: Géostatistique et géologie minières 5

GLQ3651: Géologie minière
 1. Définitions des ressources et norme canadienne NI43-101
Norme canadienne NI43-101

Définition : Le Règlement NI43-101 est une norme qui régit la manière dont les
entreprises divulguent au public au Canada des informations scientifiques et
techniques sur des projets miniers.

➢Norme canadienne pour les compagnies canadiennes

• Rapport technique sur les ressources importantes, préparé dans un format prescrit;
• La personne qualifiée doit être indépendante de l'entreprise et de la propriété (dans
certaines circonstances);
• Une entreprise est tenue d'utiliser la terminologie spécifiée définie par le CIM lors
de la divulgation des ressources minérales et des réserves minérales.

➢Respecter les conseils du CIM… (pas NI43-101… mais CIM…)

• CIM → Canadian Institute of Mining, Metallurgy and Petroleum
• Définitions des concepts :
https://mrmr.cim.org/media/1128/cim-definition-standards_2014.pdf

GLQ3401: Géostatistique et géologie minières 6

GLQ3651: Géologie minière
 1. Définitions des ressources et norme canadienne NI43-101
Personne qualifiée

Définition : Une personne physique qui remplit les trois conditions suivantes :

1. il s'agit d'un ingénieur ou d'un géoscientifique qui compte au

moins cinq ans d'expérience dans le domaine de l'exploration
minérale, de l'aménagement ou de l'exploitation de mines, ou de
l'évaluation de projets miniers, ou dans une combinaison de ces
domaines;

2. elle a une expérience pertinente à l'objet du projet minier et du

rapport technique;

3. elle est membre en règle d'une association professionnelle.

GLQ3401: Géostatistique et géologie minières 7

GLQ3651: Géologie minière
 1. Définitions des ressources et norme canadienne NI43-101
Ressources

Définition : Une ressource minérale est une concentration ou une occurrence de

matière solide d'intérêt économique dans ou sur la croûte terrestre
sous une forme, une teneur, et une quantité telles qu'il existe des
perspectives raisonnables d'extraction économique éventuelle.

L'emplacement, la quantité, la teneur, la continuité et les autres

caractéristiques géologiques d'une ressource minérale sont connus,
estimés ou interprétés à partir de preuves et de connaissances
géologiques spécifiques, y compris l'échantillonnage.

GLQ3401: Géostatistique et géologie minières 8

GLQ3651: Géologie minière
 1. Définitions des ressources et norme canadienne NI43-101
Ressources minérales : les trois types
Définition :
Présumée: partie d'une ressource minérale dont la quantité et la teneur sont estimées sur la
base de preuves géologiques et d'échantillonnages limités. Les preuves géologiques
sont suffisantes pour impliquer, mais non pour vérifier la continuité géologique et la
teneur.
Indiquée : partie d'une ressource minérale pour laquelle la quantité, la teneur, les densités, la
forme et les caractéristiques physiques sont estimées avec suffisamment de
confiance pour permettre l'application de facteurs de modification de façon
suffisamment détaillée pour soutenir la planification de la mine et l'évaluation de la
viabilité économique du gisement.

Mesurée : partie d'une ressource minérale pour laquelle la quantité, la teneur, les densités, la
forme et les caractéristiques physiques sont estimées avec suffisamment de
confiance pour permettre l'application de facteurs modificateurs afin de soutenir la
planification détaillée de la mine et l'évaluation finale de la viabilité économique du
gisement.
GLQ3401: Géostatistique et géologie minières 9
GLQ3651: Géologie minière
 1. Définitions des ressources et norme canadienne NI43-101
Ressources minérales versus réserves

Définition : Une réserve minérale est la partie économiquement exploitable d'une ressource
minérale mesurée et/ou indiquée.

• Dois être démontrée par une étude de préfaisabilité ou une étude de

faisabilité;

• Les études démontrent que, au moment de la déclaration, l'extraction

peut être raisonnablement justifiée.

GLQ3401: Géostatistique et géologie minières 10

GLQ3651: Géologie minière
 1. Définitions des ressources et norme canadienne NI43-101
Ressources minérales versus réserves : la clarification

Résultats
d’exploration

Ressources Réserves
minérales minérales

Augmenter le niveau de Présumées

confiance et de
Indiquées Probables
connaissance géologique
Mesurées Prouvées

Dépends de certains facteurs liés à l'exploitation minière,

au traitement, à la métallurgie, aux infrastructures, à
l'économie, au marketing, à la législation, à
l'environnement, à la société et au gouvernement.

GLQ3401: Géostatistique et géologie minières 11

GLQ3651: Géologie minière
 1. Définitions des ressources et norme canadienne NI43-101
Ressources minérales versus réserves : la clarification

Résultats Technique de la
d’exploration passe d’estimation
Ressources Réserves
minérales minérales

Augmenter le niveau de Présumées

confiance et de
Indiquées Probables
connaissance géologique
Mesurées Prouvées

Dépends de certains facteurs liés à l'exploitation minière,

au traitement, à la métallurgie, aux infrastructures, à
l'économie, au marketing, à la législation, à
l'environnement, à la société et au gouvernement.

GLQ3401: Géostatistique et géologie minières 12

GLQ3651: Géologie minière
 1. Définitions des ressources et norme canadienne NI43-101
Contenue obligatoire d’un rapport technique NI43-101
1. Page titre 16. Vérification des données
2. Table des matières 17. Terrains adjacents
3. Résumé 18. Essais de traitement des minerais et essais
4. Introduction et mandat métallurgiques
5. Mise en garde 19. Estimation des ressources minérales et des
6. Description et emplacement du terrain réserves minérales
7. Accessibilité, climat, ressources locales, 20. Autres données et informations pertinentes
infrastructures et géographie physique 21. Interprétations et conclusions
8. Historique 22. Recommandations
9. Contexte géologique 23. Références
10. Types de gîtes minéraux 24. Data
11. Minéralisation 25. Règles supplémentaires pour les rapports
12. Travaux d’exploration techniques sur les terrains de l’aménagement et
13. Forage sur les terrains en production
14. Méthode d’échantillonnage et approche 26. Illustrations
15. Préparation, analyse et sécurité des
échantillons
GLQ3401: Géostatistique et géologie minières 13
GLQ3651: Géologie minière
 2. Estimation et influence de la densité
Estimation : Comment obtenir un modèle de blocs ?

GLQ3401: Géostatistique et géologie minières 14

GLQ3651: Géologie minière
 2. Estimation et influence de la densité
Estimation

𝑡4
tv 𝑡5
𝑛
∗
𝑡2
𝑡0 = ? 𝑡0 =෍ λ𝑖 𝑡𝑖
𝑖=1
𝑡1 𝑡3
Pour un bloc:
Intégrer les valeurs ponctuelles dans le bloc ou discrétiser le bloc

GLQ3401: Géostatistique et géologie minières 15

GLQ3651: Géologie minière
 2. Estimation et influence de la densité
Calcul de la densité théorique

𝑡1 , 𝑉1 , ρ1 𝑡2 , 𝑉2 , ρ2 𝑉1 ρ1 + 𝑉2 ρ2
ρ=
𝑉1 + 𝑉2

𝑉 = 𝑉1 + 𝑉2 𝑡1 𝑉1 ρ1 + 𝑡2 𝑉2 ρ2
𝑡=
ρ=? 𝑉1 ρ1 + 𝑉2 ρ2
𝑡 =?

La quantité de métal dans un bloc est 𝑡𝑉ρ = 𝑡𝑀, 𝑀 la masse totale du bloc.
Pour une même teneur et un même volume, il y a plus de métal dans un bloc à forte densité

Souvent, les variations de densité peuvent être considérées comme

négligeables par rapport aux variations des teneurs ou des volumes
GLQ3401: Géostatistique et géologie minières 16
GLQ3651: Géologie minière
 2. Estimation et influence de la densité
Que faire lorsque l’épaisseur et/ou la densité varient :
Généralement, on estime la quantité de métal dans un volume donné et on divise
par le tonnage estimé pour ce même volume.
Soit : 𝑡𝑖 : Teneur au point i
λ∗𝑖 : Facteur de pondération dépendant de la méthode utilisée
𝑒𝑖 : Épaisseur au point i
𝜌𝑖 : Densité au point i
σ𝑖 λ𝑖 𝑡𝑖
À un point « 0 », on estime la teneur par : 𝑡0∗ =
σ𝑗 λ𝑗
Avec λ𝑖 , le facteur total de pondération :
𝝀𝒊 = 𝝀∗𝒊 quand « t » varie
λ𝑖 = λ∗𝑖 𝑒𝑖 quand « t » et « e » varient
λ𝑖 = λ∗𝑖 𝑒𝑖 𝜌𝑖 quand « t » , « e » et « 𝜌 » varient
GLQ3401: Géostatistique et géologie minières 17
GLQ3651: Géologie minière
 3. Notions de statistique
Moyenne arithmétique Écart-type des teneurs : Variance des teneurs
des teneurs :
𝑛 𝑛 𝑛 𝑛

𝑡ҧ = ෍ λ𝑖 𝑡𝑖 𝜎= ෍ λ2𝑖 𝑡𝑖 − 𝑡ҧ 2 = ෍ λ𝑖 𝑡𝑖 2 − 𝑡ҧ 2 𝜎 2 = ෍ λ2𝑖 𝑡𝑖 − 𝑡ҧ 2

𝑖=1 𝑖=1 𝑖=1 𝑖=1

Espérance : Variance centrée d’ordre 2 :

∞

𝔼[𝑋] = ෍ λ𝑖 𝑥𝑖 2
𝑉𝑎𝑟 𝑋 = 𝔼 𝑋−𝔼 𝑋
𝑖=1
Coefficient d’asymétrie: Coefficient Kurtosis :
3 4
𝑋−𝜇 𝑋−𝜇
𝛾1 = 𝔼 𝛽2 = 𝔼
𝜎 𝜎

GLQ3401: Géostatistique et géologie minières 18

GLQ3651: Géologie minière
 3. Notions de statistique

La distribution des teneurs se résume par certaines statistiques

La moyenne (m) : tendance centrale L’écart-type (σ) : la dispersion

autour de la moyenne

GLQ3401: Géostatistique et géologie minières 19

GLQ3651: Géologie minière
 3. Notions de statistique
La distribution des teneurs se résume par certaines statistiques

Coefficient
Asymétrie : moment d’ordre d’asymétrie
trois de la variable centrée
réduite =0.00

Positif : dérive vers la gauche

=0.44
Négatif : dérive vers la droite

=3.84

GLQ3401: Géostatistique et géologie minières 20

GLQ3651: Géologie minière
 3. Notions de statistique
Le type de distribution peut influencer sur la valeur du gisement (asymétrie)
Coefficient
d’asymétrie
Teneur de coupure : 5
=0.00

Diminution du profit et de la teneur

lorsque l’asymétrie augmente

=0.44

=3.84

GLQ3401: Géostatistique et géologie minières 21

GLQ3651: Géologie minière
 3. Notions de statistique
Stationnarité de la moyenne

Diviser les domaines

GLQ3401: Géostatistique et géologie minières 22

GLQ3651: Géologie minière
 3. Notions de statistique
Stationnarité de la moyenne (Stationnarité intrinsèque)

GLQ3401: Géostatistique et géologie minières 23

GLQ3651: Géologie minière
 3. Notions de statistique
Stationnarité (moyenne et variance)
Stationnarité du 2ième ordre Stationnarité intrinsèque

GLQ3401: Géostatistique et géologie minières 24

GLQ3651: Géologie minière
 4. Méthodes conventionnelles d’estimation des ressources
Méthode des polygones (plus proche voisin)

Principe : La teneur, 𝑡0∗ , estimée en un point, 𝑥0 , est égale à la teneur du

point connu le plus proche

→ On doit définir des polygones (polygones de Voronoï) à teneur constante

Le même principe est appliqué pour l’épaisseur

GLQ3401: Géostatistique et géologie minières 25

GLQ3651: Géologie minière
 4. Méthodes conventionnelles d’estimation des ressources
Méthode des polygones (plus proche voisin)
Méthode des polygones
100

90

80

70

60

50
y

40

30

20

10

0
0 20 40 60 80 100
x
GLQ3401: Géostatistique et géologie minières 26
GLQ3651: Géologie minière
 4. Méthodes conventionnelles d’estimation des ressources
Méthode des polygones (plus proche voisin)
Méthode des polygones
100

90

80

70

60

50
y

Le seul paramètre à
40 spécifier est la règle
30 de fermeture pour les
20 polygones externes
10
Souvent: distance max
0
0 20 40 60 80 100
x 27
GLQ3401: Géostatistique et géologie minières
GLQ3651: Géologie minière
 4. Méthodes conventionnelles d’estimation des ressources
Méthode des polygones (plus proche voisin)
Méthode des polygones Teneur
100
9
90
8
80
7
Design de la mine, 70

la teneur moyenne 60
6

est la moyenne 5
50
y

pondérée par les

« surfaces » 40 4

30 3

20 2

10 1

0 0
0 20 40 60 80 100

GLQ3401: Géostatistique et géologie minières x 28

GLQ3651: Géologie minière
 4. Méthodes conventionnelles d’estimation des ressources
Méthode des polygones (plus proche voisin)
Comment obtenir les polygones ?

1. Obtenir la triangulation de Delaunay;

Triangles les + équilatéraux possibles
2. Tracer les perpendiculaires au milieu des côtés des
triangles (médiatrices);
3. Identifier l’intersection des 3 médiatrices se rencontrant
au centre du cercle circonscrit au triangle;
4. Relier les points d’intersection de deux triangles
partageant un côté commun → Un côté du polygone.

GLQ3401: Géostatistique et géologie minières 29

GLQ3651: Géologie minière
 4. Méthodes conventionnelles d’estimation des ressources
Méthode des polygones (plus proche voisin)
Comment obtenir les polygones ?
Triangularisation de Delaunay

B B
A A

D D

C C

D est dans le cercle ABC. Le triangle D n’est pas dans le cercle ABC. Le triangle
ABC n’est pas Delaunay ABC est Delaunay
GLQ3401: Géostatistique et géologie minières 30
GLQ3651: Géologie minière
 4. Méthodes conventionnelles d’estimation des ressources
Méthode des polygones (plus proche voisin)
Comment obtenir les polygones ?
Tracer les médiatrices Méthode des polygones
100

90

80
Note: 70

Les côtés des polygones sont 60

les médiatrices des triangles 50

y
de Delaunay. 40

30

20

10

0
0 20 40 60 80 100
x
GLQ3401: Géostatistique et géologie minières 31
GLQ3651: Géologie minière
 4. Méthodes conventionnelles d’estimation des ressources
Méthode des polygones (plus proche voisin)
Comment obtenir les polygones ?
Identifier le point milieu Méthode des polygones
100

90

80

70

60

50

y
40

30

20

10

0
0 20 40 60 80 100
x
GLQ3401: Géostatistique et géologie minières 32
GLQ3651: Géologie minière
 4. Méthodes conventionnelles d’estimation des ressources
Méthode des polygones (plus proche voisin)
Comment obtenir les polygones ?
Relier les points d’intersection Méthode des polygones
100

90

80

70

60

50

y
40

30

20

10

0
0 20 40 60 80 100
x 33
GLQ3401: Géostatistique et géologie minières
GLQ3651: Géologie minière
 4. Méthodes conventionnelles d’estimation des ressources
Méthode des polygones (plus proche voisin)
Estimations ponctuelles et blocs
Principe : Tous les points dans un polygone reçoivent la teneur de la 15
histogramme des données

donnée associée au polygone

10

Méthode des polygones

Exemple : 100
5

Une grille régulière → polygones 80

0
60 1 1.5 2 2.5 3
sont des carrés

y
40 histogramme des blocs (estimés)
15

20
Estimés des blocs → même 10
0
distribution statistique que les 0 50
x
100

données 5

Est-ce réaliste d’un point de vue statistique ?

0
1 1.5 2 2.5 3

GLQ3401: Géostatistique et géologie minières 34

GLQ3651: Géologie minière
 4. Méthodes conventionnelles d’estimation des ressources
Méthode des polygones (plus proche voisin)
Estimations ponctuelles et blocs
• Peut bien estimer la proportion de points dépassant un seuil;
application en environnement
• Ne permet pas de bien déterminer quels points dépassent le seuil car méthode
peu précise;
• Ne permets pas de déterminer la proportion de blocs ni quels blocs dépassent
un seuil;
• Méthode présentant un fort biais conditionnel.

GLQ3401: Géostatistique et géologie minières 35

GLQ3651: Géologie minière
 4. Méthodes conventionnelles d’estimation des ressources
Méthode des triangles
Pour estimation ponctuelle et de blocs ∗
𝑎
𝑡13 = 𝑡1 + 𝑡 − 𝑡1
Il s’agit d’une interpolation linéaire 𝑎+𝑏 3

𝑡3 𝑐
𝑡0∗ = 𝑡13
∗
+ ∗
𝑡2 − 𝑡13
𝑐+𝑑
b
Pour trouver les valeurs de 𝑎, 𝑏, 𝑐 𝑒𝑡 𝑑, il faut
déterminer la droite (𝑡1 𝑡3 ) et la droite (𝑡2 𝑡0∗ ).
∗
𝑡13
c 𝑡0∗ L’intersection des deux droites fournit les
a + d coordonnées du point 𝑡13∗
.
𝑡1 𝑡2
Pythagore permet de déterminer toutes les
Blocs : moyenne des teneurs distances nécessaires.
ponctuelles estimées dans le bloc
GLQ3401: Géostatistique et géologie minières 36
GLQ3651: Géologie minière
 4. Méthodes conventionnelles d’estimation des ressources
Comparaison (polygones versus triangles) 100 100
11

0.8
0.8
90 90

80 80 0.6
0.6

70 70 0.4
0.4

1 60 0.2
0.2
100
100 60

0.8 50 50 00
90
90
Triangles 40 40 -0.2
-0.2
80
80 0.6
30 30 -0.4
-0.4
70
70 0.4 20 20 -0.6
-0.6

60
60 0.2 10 10 -0.8
-0.8
0 0
50
50 0 -1
-1
0 0 20 20 40 40 60 60 80 80 100
100

40
40 -0.2
1
100 100
30
30 -0.4
90 0.8
90
20
20 -0.6 80 0.6
80
10
10 -0.8 70 70 0.4

00 60 60 0.2
-1
00 2020 4040 6060 8080 100
100
Polygones 50 50 0

40 40
Réalité -0.2

30 30 -0.4

20 20
-0.6

10 10
-0.8
0 0
GLQ3401: Géostatistique et géologie minières 0 20 20 40 -1 37
0 40 60 60 80 80 100
100
GLQ3651: Géologie minière
 4. Méthodes conventionnelles d’estimation des ressources
Inverse de la distance
L’estimateur est de la forme :

𝑡𝑖
σ𝑖 𝑏 𝑛
∗ 𝑑𝑖 1/𝑑𝑖𝑏
𝑡0 = = ෍ λ𝑖 𝑡𝑖 𝑎𝑣𝑒𝑐, λ𝑖 =
1 1
σ𝑖 𝑏 𝑖=1 σ𝑗 𝑏
𝑑𝑖 𝑑𝑗
où 𝑑𝑖 est la distance entre le point à estimer et le 𝑖 ème point observé.

GLQ3401: Géostatistique et géologie minières 38

GLQ3651: Géologie minière
 4. Méthodes conventionnelles d’estimation des ressources
Inverse de la distance σ𝑖
𝑡𝑖
𝑏
Exemple 𝑑 𝑖
𝑡0∗ =
1
Estimation au point A avec b=2 σ𝑖 𝑏
𝑑𝑖

# distance teneur %
1 40 1
2 40 1
3 30 1.5
4 35 1.5
5 20 3

GLQ3401: Géostatistique et géologie minières 39

GLQ3651: Géologie minière
 4. Méthodes conventionnelles d’estimation des ressources
Inverse de la distance 𝑡𝑖
σ𝑖
Exemple ∗ 𝑑𝑖𝑏
𝑡0 =
1
Estimation au point A avec b=2 σ𝑖 𝑏
𝑑𝑖
# distance teneur %
1 40 1
2 40 1
3 30 1.5
4 35 1.5
5 20 3

(1/40 2 + 1/40 2 + 1.5/30 2 + 1.5/35 2 + 3/20 2 )

t= = 2.05%
(1/40 + 1/40 + 1/30 + 1/35 + 1/20 )
2 2 2 2 2

GLQ3401: Géostatistique et géologie minières 40

GLQ3651: Géologie minière
 4. Méthodes conventionnelles d’estimation des ressources
Inverse de la distance 𝑡𝑖
σ𝑖
Exemple ∗ 𝑑𝑖𝑏
𝑡0 =
1
Estimation au point A avec b=2 σ𝑖 𝑏
𝑑𝑖
# distance teneur %
1 40 1
2 40 1
3 30 1.5
4 35 1.5
5 20 3
Pour diminuer le nombre de données dans l’estimation, on a recours à un
voisinage.
P. ex., on conserve que les points dont la distance est inférieure ou égale à 35’
GLQ3401: Géostatistique et géologie minières 41
GLQ3651: Géologie minière
 4. Méthodes conventionnelles d’estimation des ressources

Exercices en équipe
1) Estimation par méthodes
conventionnelles

GLQ3401: Géostatistique et géologie minières 42

GLQ3651: Géologie minière
 4. Méthodes conventionnelles d’estimation des ressources

20
𝑍1 = 6𝑝𝑝𝑚

16

𝑍4 = 10𝑝𝑝𝑚
12 𝑍0 =?

8
𝑍3 = 7𝑝𝑝𝑚
𝑍2 = 2𝑝𝑝𝑚
4

2 4 6 8 10 12
GLQ3401: Géostatistique et géologie minières 43
GLQ3651: Géologie minière
 4. Méthodes conventionnelles d’estimation des ressources
Inverse de la distance
Distance anisotrope

𝑥 2 + 𝑎𝑦 2

si a > 1 , alors un écart unitaire selon y correspond à une distance plus

grande que le même écart en x

B A est plus près de B que de C,

pourtant il est logique de
Strates considérer C plus semblable à A
(que B) mettre a>>1 permet de
A C l’imposer.

GLQ3401: Géostatistique et géologie minières 44

GLQ3651: Géologie minière
 4. Méthodes conventionnelles d’estimation des ressources
Note:

• Les méthodes sont définies pour une estimation ponctuelle;

• Pour un bloc, faire la moyenne des estimations ponctuelles dans le bloc;

• Si l’épaisseur (et/ou densité) varie, habituellement il faut estimer

l’accumulation (𝑎∗ ) et l’épaisseur (𝑒 ∗ ) séparément et calculer la teneur
comme 𝑡 ∗ = 𝑎∗ /𝑒 ∗ .

GLQ3401: Géostatistique et géologie minières 45

GLQ3651: Géologie minière
 4. Méthodes conventionnelles d’estimation des ressources

Exercices en équipe
2) Estimation par la
technique d’accumulation

GLQ3401: Géostatistique et géologie minières 46

GLQ3651: Géologie minière
 5. Critique des méthodes conventionnelles
• Inverse de la distance → plus de flexibilité ;
• Triangles → bons résultats si les données sont abondantes et de qualité;
Par exemple, représenter une topographie.
Extrapolation hasardeuse.
• Polygonale → à éviter;
Imprécise et ne fournit pas de teneurs réalistes pour des blocs.
• Localement, polygonale (1 donnée), triangles (3 données) et inverse de la
distance (habituellement 5-50);
• Triangle, polygonale et inverse de la distance : interpolateur exact ;
Choix discutable lorsque les données sont entachées d’erreur
• Triangle et polygonale : surtout problèmes 2D (ex. veine; niveau d’une mine);
• Comparer les performances des méthodes par la technique de la validation
croisée.
GLQ3401: Géostatistique et géologie minières 47
GLQ3651: Géologie minière
 6. Validation croisée
Objectifs :
Permet de valider une méthode d’estimation
Déterminer les paramètres d’une méthode d’estimation (p. ex., exposant b);
Déterminer l’importance et l’orientation d’une éventuelle anisotropie;
Distance maximale utilisée pour sélectionner les données lors de l’estimation;
Déterminer le nombre de données à utiliser pour l’estimation;
Comparer des méthodes d’estimation.

Technique du ‘leave-one-out’ :
Consiste à retirer une observation pour ensuite l’estimer à
partir des autres données observées.
Cela est répété pour tous les points.
𝑒𝑖 = 𝑡 𝑥𝑖 − 𝑡 ∗ 𝑥𝑖

Statistique pour comparaison :

1 1
෍ 𝑒𝑖 ≈ 0 , min ෍ 𝑒𝑖 𝑜𝑢 min ෍ 𝑒𝑖2
𝑛 𝑛
𝑖 𝑖 𝑖
GLQ3401: Géostatistique et géologie minières 48
GLQ3651: Géologie minière
 6. Validation croisée
Comparaison visuelle :
• Histogramme des résidus et résidus normalisés
• Carte des résidus et résidus normalisés

Diviser en deux
zones distinctes ?

GLQ3401: Géostatistique et géologie minières 49

GLQ3651: Géologie minière
 6. Validation croisée
Comparaison visuelle :
Sol contaminé à l’arsenic (données de l’EPA)
Validation croisée
0.6

0.59

0.58
Moyenne |e|

0.57

0.56

0.55

0.54 b=2.5 procure en

moyenne l’erreur
0.53
absolue minimale
0.52
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
Exposant b

GLQ3401: Géostatistique et géologie minières 50

GLQ3651: Géologie minière
 6. Validation croisée
100

Comparaison visuelle : 90

80

70

60

50

40

30

20

10

0
0 20 40 60 80 100
1 1
100 100 100 100
0.890 0.8
90 90 90

80 0.680 0.6 80
80

70 70 0.470 0.4 70

60 60 0.260 0.2 60

50 50 0 50 0 50

40 40 40
-0.2 -0.2 40

30 30 30
-0.4 -0.4 30

20 20 20 -0.6 20
-0.6

10 10 10 -0.8 10
-0.8

0 0 0 0
-1 -1
0 20 40 60 80 100 0 20 40 60 80 100 0 20 40 60 80 100 0 20 40 60 80 100

GLQ3401: Géostatistique et géologie minières 51

GLQ3651: Géologie minière
 6. Validation croisée
100

Comparaison visuelle : 90

80

70

60

Statistique : ෍ 𝑒𝑖
50

40

30
𝑖 20

10

0
0 20 40 60 80 100
1 1
100 100 100 100
0.890 0.8
90 90 90

80 0.680 0.6 80
80

70 70 0.470 0.4 70

60 60 0.260 0.2 60

50 50 0 50 0 50

40 40 40
-0.2 -0.2 40

30 30 30
-0.4 -0.4 30

20 20 20 -0.6 20
-0.6

10 10 10 -0.8 10
-0.8

0 0 0 0
-1 -1
0 20 40 60 80 100 0 20 40 60 80 100 0 20 40 60 80 100 0 20 40 60 80 100

Polygone Triangles Inv dist, b=1 Inv dist b=2

17 13 22 16
GLQ3401: Géostatistique et géologie minières 52
GLQ3651: Géologie minière
 7. Méthode des sections
Principe : Consiste à tracer une suite de sections (habituellement parallèles) à travers
un gisement, à déterminer les zones de minerais dans chaque section et à
compiler le volume total en utilisant les distances entre chaque section.

• Conserver l’approche géométrique pour définir l’enveloppe du gisement

en 2D (i.e. forages → définir l’enveloppe minéralisée sur chaque section)
• Passage section → 3D par « modélisation solide »
• Estimer les teneurs séparément dans un modèle de blocs (p. ex., par
krigeage)
• Intersection modèle de blocs et géométrie du gisement → teneur du
solide
Cette méthode est particulièrement indiquée dans le cas de
gisements en forme de veines ou de lentilles. Les sections parallèles
sont habituellement orthogonales au corps minéralisé.
GLQ3401: Géostatistique et géologie minières 53
GLQ3651: Géologie minière
 7. Méthode des sections – Déterministique (explicite)
Section S1 a y=0 apres interpolation
1

Points servant à définir le polygone

Coord z

0.5
Le tracé est interpolé en un grand
nombre de points
0
0 0.5 1
Coord x
Section S2 a y=1 apres interpolation
1
0.8
0.7
Coord z

0.5 0.6

0.5
0.4
0 2
0 0.5 1
0.8
Coord x 1 0.6
0.4
Section S3 a y=2 apres interpolation 0.2
1 0
Coord z

0.5
Les points d’interpolation de
deux sections consécutives
0
0 0.5 1 sont joints par des triangles
Coord x

GLQ3401: Géostatistique et géologie minières 54

GLQ3651: Géologie minière
 7. Méthode des sections – Déterministique (explicite)
Sept figures de cas simple

GLQ3401: Géostatistique et géologie minières 55

GLQ3651: Géologie minière
 7. Méthode des sections – Déterministique (implicite)

GLQ3401: Géostatistique et géologie minières 56

GLQ3651: Géologie minière
 8. Calcul de la densité théorique
Certaines mines calculent la densité du minerai à partir des analyses chimiques obtenues
Deux approches :
• formule empirique obtenue par régression;
• calcul basé sur la minéralogie déduite de l’analyse chimique
Exemple : Gisement de Cu
• Le Cu est contenu dans la chalcopyrite (ρ=4.2, la teneur en Cu dans chalcopyrite est de 35%);
• La chalcopyrite est le seul sulfure présent;
• Les autres minéraux ont une densité voisine de 3;

Déterminer la densité d’un minerai ayant 1% Cu ? 5% Cu ?

GLQ3401: Géostatistique et géologie minières 57

GLQ3651: Géologie minière
 8. Calcul de la densité théorique
Pour 1% Cu
1 % Cu → 1/0.35=2.86% chalcopyrite
100 g roche → 2.86 g chalco → volume de chalco = 2.86 g / 4.2 g/ cm3 → 0.68 cm3
→ 97.14 g gangue → volume de gangue → 97.14 g / 3 g/ cm3 → 32.38 cm3
Volume total = 33.06 cm3
Masse volumique théorique = 100g / 33.06 cm3 = 3.02 g/ cm3 → 𝝆 = 𝟑. 𝟎𝟐

Pour 5% Cu
5 % Cu → 5/0.35=14.29% chalcopyrite
100 g roche → 14.29 g chalco → volume de chalco = 14.29 g / 4.2 g/ cm3 → 3.40 cm3
→ 85.71 g gangue → volume de gangue → 85.71 g / 3 g/ cm3 → 28.57 cm3
Volume total = 31.97 cm3
Masse volumique théorique = 100g / 31.97 cm3 = 3.13 g/ cm3 → 𝝆 = 𝟑. 𝟏𝟑

GLQ3401: Géostatistique et géologie minières 58

GLQ3651: Géologie minière
 8. Calcul de la densité théorique

Densité vs teneur en Cu
4.5
Note: 3.13  0.857*3+0.143*4.2=3.17

4
Densite

La variation de la densité n’est pas linéaire

en fonction de la teneur en Cu
3.5

3
0 10 20 30 40
Teneur en Cu, (%)

GLQ3401: Géostatistique et géologie minières 59

GLQ3651: Géologie minière
 8. Calcul de la densité théorique
Que faire si on a plus qu’un sulfure (p. ex., chalcopyrite et pyrite) ?
% Cu → % chalcopyrite
% S – (%S dans chalcopyrite) → % pyrite
Seulement si on suppose que S provient uniquement de la pyrite et de la chalcopyrite

Cas général → système d’équations linéaires 𝐴𝑥 = 𝑏

𝐴(𝑖, 𝑗) : teneur de l’élément « i » dans le minéral « j » connu par la formule chimique du minéral
𝑏(𝑖) : teneur de l’élément « i » dans la roche connu par l’analyse chimique
𝑥(𝑗) : teneur du minéral « j » dans la roche, x(n) est la à déterminer
gangue.
On a aussi la contrainte que σ 𝑥(𝑗) = 1

Pour réduire le nombre d’inconnues, on n’isole que les minéraux ayant

une densité nettement différente de la « gangue ».
GLQ3401: Géostatistique et géologie minières 60
GLQ3651: Géologie minière
 8. Calcul de la densité théorique
Exemple : Gisement de Cu et de Zn. Poids atomique et formule
• Le Cu dans la chalcopyrite (CuFeS2) et la chalcocite stœchiométrique :
(Cu2S);
35% Cu dans chalcopyrite
• Le Zn dans la sphalérite (ZnS);
• La roche contient aussi de la pyrite (FeS2);
80% Cu dans chalcocite
• Il y a en moyenne 2% de Fe dans la gangue, mais 67% Zn dans la sphalérite
pas de soufre. 35% S dans la chalcopyrite
20% S dans la chalcocite
Analyse → 6% Cu, 9% Zn, 10% S, 5% Fe 33% S dans la sphalérite
53% S dans la pyrite
30% Fe dans la chalcopyrite
47% Fe dans la pyrite

Cu-  0.35 00.80

.50 0 0 0  Chalcopyri te 0.06
Zn-   0 0 0.67 0 0   Chalcocite  0.09
   
S-  0.35 0.20 0.33 0.53 0   Sphalérite  = 0.10
     
Fe-  0.30 0 0 0.47 0.02  Pyrite   0.05
 1 1 1 1 1   Gangue   1 
GLQ3401: Géostatistique et géologie minières 61
GLQ3651: Géologie minière
 8. Calcul de la densité théorique
Solution : Calcul de la densité théorique :
Chalcopyrite   0 
 Chalcocite   7.7  Volume pour 100 g de roche
   
 Sphalérite  = 13.43 % chalcopyrite : 0 g/ 4.1g/cm3= 0 cm3
   
 Pyrite   7.9  chalcocite : 7.7/5.6 = 1.38 cm3
 Gangue   71.4 
    sphalérite : 13.43/4.1= 3.28 cm3
pyrite : 7.9/5.0= 1.58 cm3
Densité
Chalcopyrite : 4.1 gangue : 71.4/2.9= 24.62 cm3
Chalcocite : 5.6
Sphalérite : 4.1
Pyrite : 5.0 Volume total: 30.86 cm3
Gangue : 2.9
Masse volumique : 100 g/30.86 cm3 = 3.24 g/ cm3
Densité théorique: 3.24
GLQ3401: Géostatistique et géologie minières 62
GLQ3651: Géologie minière
 8. Calcul de la densité théorique
Effet de la porosité sur la densité :

Vvide
porosité  n =
Vroche + Vvide

M roche Vroche  théorique

 réelle = = = (1 - n) théorique
Vroche + Vvide Vroche + Vvide

GLQ3401: Géostatistique et géologie minières 63

GLQ3651: Géologie minière
 8. Calcul de la densité théorique

Exercices en équipe
4) Calcul de la densité
théorique et des teneurs
par une analyse chimique

GLQ3401: Géostatistique et géologie minières 64

GLQ3651: Géologie minière