Rapport de projet de fin d’étude

Planification longe terme de l’exploitation souterraine de la

mine de Tizert

Réalisé par : Tuteur industriel :

Samir EL OTEBY M. Youssef ADLAOUI

Tuteur académique :
Pr. Damien Goetz

Période du 12 avril 2021 au 12 septembre 2021

 Remerciement :

Avant de commencer le développement de ce travail de fin d’étude il me parait tout naturel de

remercier toutes les personnes qui ont contribués de près ou de loin au sucées de mon travail d’option
chez Managem.

Je tiens tout d’abord à remercier les plus beaux parents du monde, ceux qui m’ont soutenu toute au
long de mon parcours académique et bien plus avant durant les premières années de la grande
expérience, celle de la vie.

Je remercie énormément mon professeur encadrant Pr. Damien Goetz pour ses conseils qui m’ont été
enrichissants, pours sa serviabilité et pour sa disponibilité.

Je remercie M. Tarik SOUROR directeur du REMINEX Mine et Carrière pour m’avoir offert la chance
d’effectuer ce stage.

Je remercie mon encadrant industriel M. Youssef ADLAOUI pour ses conseils précieux et son
accompagnement durant toute la période du stage.

Je tiens à exprimer ma gratitude et ma reconnaissance à mon parrain de stage M. Mbacké DIOP Pour
le temps qu’il a consacré pour répondre à mes questions et pour m’aider à développer ce sujet. Je le
remercie pour sa bienveillance et sa gentillesse. Je tiens aussi à remercier l’équipe de l’ingénierie :
Mme Atika EL AOUAM, Mme Jamila OUABID , Mme Halima ESSAHRAOUI, M. Ismail ACHKARA,
M. Noureddine RAISS Pour leurs conseils et leurs aides.

Finalement, je remercie l’ensemble du cadre professoral de l’université Mohammed IV Polytechnique

en général et de l’école EMINES-School of industrial Management en particulier.

2
 Résumé

Le groupe Managem est une société marocaine qui opère dans l’activité minière et hydrométallurgie avec
une expertise de 90 ans dans l’extraction, la valorisation et la commercialisation des métaux industriels
et des métaux précieux.

Le groupe vise à exploiter un gisement cuprifère de la région de Tizert à travers sa filiale AKA GOLD MINING
(AGM), c’est un gisement situé sur le bord nord de la boutonnière précambrienne d’Ighrem avec des
ressources estimées à 62 millions de tonnes avec des teneurs de minerai qui atteignent 0,96% de cuivre
et 20g/t d’argent. Une telle ressource donnera lieu à une exploitation de 3.3 Mt de minerai par an.

Afin de pourvoir lancer le travail, la filiale de l’ingénierie et de la recherche REMINEX a fait une étude de
faisabilité en collaboration avec Met-Chem, un groupe canadien qui offre le service d’ingénierie de
l’activité minière, pour juger sur le potentiel économique du projet. Cependant, l’apparition d’une
extension dans la direction Nord-Ouest du gisement principal a fait relancer une étude d’orientation afin
de transformer les ressources en réserve ainsi que la détermination de l’approche d’exploitation
convenable aux paramètres morphologiques et géotechniques, une optimisation du séquençage de la
production est aussi nécessaire afin de récupérer le maximum de la richesse.

Ce travail est donc pour objectif de mener une étude d’orientation technique et économique de
l’extension, et ensuite de définir une planification long terme qui optimise le séquençage de la production.

3
 Abstract

Managem Group is a Moroccan company operating in the mining and hydrometallurgy business with a
90-years expertise in the extraction, valorization, and marketing of industrial and precious metals.

The group aims today to develop a copper deposit in the Tizert region (250 km south of Marrakech, 80 km
east of Agadir) through its subsidiary AKA GOLD MINING (AGM). The Tizert deposit is located on the
northern edge of the Precambrian buttonhole of Ighrem with resources estimated at 62 million tons of
ore with ore grades that reach on average 0.96% of copper and 20g/t of silver.

REMINEX, the engineering and research subsidiary of MANAGEM, carried out a feasibility study, in order
to evaluate the economic potential of the project, in collaboration with Met-Chem, a Canadian group that
offers engineering services for mining activities. However, the discovery of an extension in the
northwestern direction of the main deposit, with resources estimated at 21 million tons of ore at ore
grades that reach on average 0.99% of copper and 18g/t of silver, has prompted a new feasibility study to
transform these new resources into reserves, to determine the most suitable mining approach in regard
to the morphological and geotechnical parameters, and finally to optimize the production sequence over
the deposit.

This work is a contribution to the update of the feasibility study. It is concentrated on the mining aspects,
i.e. the selection of the mining method to be implemented in each area of the deposit, and the
optimization of the production sequence.

4
Table des matières

Résumé ......................................................................................................................................................... 3
Abstract ......................................................................................................................................................... 4
Liste des Figures: ........................................................................................................................................... 7
Liste des Tableaux : ....................................................................................................................................... 8
Chapitre I : Contexte Générale du projet...................................................................................................... 9
I. Présentation Générale de Managem : .............................................................................................. 9
1. Introduction : ................................................................................................................................ 9
2. Structure et Activité : .................................................................................................................. 10
II. Présentation de REMINEX : ............................................................................................................. 11
1. Introduction ................................................................................................................................ 11
2. Secteur d’activité ........................................................................................................................ 11
III. Présentation du gisement : ......................................................................................................... 12
1. Situation Géographique : ............................................................................................................ 12
2. Contexte Géologique : ................................................................................................................ 12
3. Synthèse lithologique :................................................................................................................ 13
4. Contexte hydrogéologique : ....................................................................................................... 13
Chapitre II : Choix de la méthode d’exploitation ........................................................................................ 14
I. Etude statistique de l’extension :.................................................................................................... 14
1. Etude statistique sur le pendage : .............................................................................................. 14
2. Etude statistique sur la puissance : ............................................................................................. 16
II. Etude géotechnique : ...................................................................................................................... 18
III. Sélection de la méthode d’exploitation : .................................................................................... 19
1. Introduction : .............................................................................................................................. 19
2. Remplissage des matrices : ......................................................................................................... 21
3. Calcul du score : .......................................................................................................................... 22
4. Résultat et synthèse : .................................................................................................................. 23
5. Description de la méthode d’exploitation retenue : .................................................................. 24
6. Application Shiny :....................................................................................................................... 24
Chapitre II : Optimisation des ressources exploitables ............................................................................... 27
I. Méthodologie :................................................................................................................................ 27

5
 II. Paramètres économiques: .............................................................................................................. 28
III. Paramètres techniques : ............................................................................................................. 31
1. La géométrie des tailles : ............................................................................................................ 31
2. Récupération minière : ............................................................................................................... 34
IV. Résultat d’optimisation : ............................................................................................................. 34
V. Infrastructure minière :................................................................................................................... 36
1. Préparation des réserves : .......................................................................................................... 36
2. Description de l’Infrastructure : .................................................................................................. 37
Chapitre III : Planification de la Production ................................................................................................ 39
I. Découpage en panneaux d’exploitation : ....................................................................................... 39
II. Optimisation du séquençage de la production:.............................................................................. 40
1. Introduction : .............................................................................................................................. 40
2. Formulation mathématique : ...................................................................................................... 41
3. Résultat : ..................................................................................................................................... 45
4. Mode d’exploitation retenue : .................................................................................................... 46
Conclusion : ................................................................................................................................................. 48
ANNEXES ..................................................................................................................................................... 49
ANNEXE CHAPITRE II : ................................................................................................................................. 50
Les matrices de comparaison par paire ...................................................................................................... 50
ANNEXE CHAPITRE II : ................................................................................................................................. 53
Le code de l’application Shiny ..................................................................................................................... 53
ANNEXE CHAPITRE III : ................................................................................................................................ 60
Les différentes modèle d’optimisation ....................................................................................................... 60

6
 Liste des Figures:
Figure 1:Lieux d’activité du groupe Managem. ............................................................................................ 9
Figure 2: Les mines de Managem au Maroc. .............................................................................................. 10
Figure 3: Situation géographique de la zone du projet Tizert. ................................................................... 12
Figure 4: Coupe géologique du gisement. .................................................................................................. 12
Figure 5: Cartographie du pendage de la zone 09. ..................................................................................... 14
Figure 6: Cartographie du pendage de la zone 35. ..................................................................................... 15
Figure 7: Boîte à moustaches du pendage. ................................................................................................. 15
Figure 8: Histogramme du pendage............................................................................................................ 15
Figure 9: Cartographie de la puissance de la zone 09................................................................................. 16
Figure 10: Cartographie de la puissance de la zone 35............................................................................... 16
Figure 11: Cartographie de la puissance de la zone 35............................................................................... 17
Figure 12: Boîte à moustaches de la puissance. ......................................................................................... 17
Figure 13: La répartition des méthodes d’exploitation de la zone 09. ....................................................... 23
Figure 14: La répartition des méthodes d’exploitation de la zone 09. ....................................................... 23
Figure 15: Schéma de la méthode Step Room and Pillar. ........................................................................... 24
Figure 16: Les zones exploitables avec la méthode chambre et pilier. ...................................................... 27
Figure 17: Répartition des ressources suivant la valeur du NSR de la zone 09. ......................................... 29
Figure 18: Répartition des ressources suivant la valeur du NSR de la zone 35. ......................................... 30
Figure 19: Dimensionnement des piliers en fonction de la profondeur. .................................................... 31
Figure 20: Chambre et pilier. ...................................................................................................................... 32
Figure 21: Différentes cas des ouvertures. ................................................................................................. 33
Figure 22: Résultat de l’optimisation MSO (Vue en plan)........................................................................... 35
Figure 23: Résultat de l’optimisation MSO (Vue en 3D). ............................................................................ 35
Figure 24: Résultat de la filtration des blocs............................................................................................... 36
Figure 25: Vue en 3D de la minéralisation et de l’infrastructure................................................................ 38
Figure 26: Vue en 3D de la minéralisation et de l’infrastructure................................................................ 38
Figure 27: Scénario 1 du découpage. Figure 28: Scénario 2 du découpage. ........... 39
Figure 29: Modélisation du scénario 1........................................................................................................ 42
Figure 30: Modélisation du scénario 2........................................................................................................ 42
Figure 31: Résultat d’optimisation pour le scénario de découpage 1
Figure 32: Résultat d’optimisation pour le scénario de découpage 2 ........................................................ 45
Figure 33: La répartition temporelle de l’exploitation des panneaux. ....................................................... 46
Figure 34: La simulation de la production sous MineSched ....................................................................... 47
Figure 35: Le développement de l'infrastructure de la simulation sous MineSched ................................. 47

7
 Liste des Tableaux :

Tableau 1: RMR par domaine lithologiques ................................................................................................ 18

Tableau 2: Classification du système RMR ................................................................................................. 18
Tableau 3: L’importance des différents critères de choix........................................................................... 19
Tableau 4: Echelle de Saaty......................................................................................................................... 20
Tableau 5: Les valeurs du ratio randomisé selon la taille. .......................................................................... 20
Tableau 6: Paramètres économiques du NSR. ............................................................................................ 28
Tableau 7: Répartition des ressources suivant le NSR de la zone 09.......................................................... 29
Tableau 8: Répartition des ressources suivant le NSR de la zone 35.......................................................... 30
Tableau 9: Paramètres géométriques des tailles........................................................................................ 32
Tableau 10: Marge opérationnelle actuelle des blocs pour le scénario 1. ................................................. 40
Tableau 11: Marge opérationnelle actuelle des blocs pour le scénario 2. ................................................. 41
Tableau 12: Le séquençage de l’exploitation selon le mode retenue. ....................................................... 46

8
 Chapitre I : Contexte Générale du projet

I. Présentation Générale de Managem :

1. Introduction :
Le groupe Managem est un acteur marocain dans l’activité minière africaine, avec une histoire de 90 ans
dans l’extraction minière, la valorisation et la commercialisation des métaux de base et des métaux
précieux. Avec plusieurs compagnies minières, le groupe se considère comme le leader de l’exploitation
minière national en plus de sa présence dans 6 pays africains. Tout au long de son existence, la société ne
cesse de faire entrer la recherche dans son plan de développement.

Depuis 1928, Managem s’est implanté dans plusieurs sites en commençant par la découverte du gisement
de Cobalt dans la région de Bou-Azzar, cette découverte donnera naissance à la compagnie Tifnout
Tighanimine (CTT) spécialisée dans la production des concentrés du Cobalt. Tout au long de son histoire,
le groupe a lancé plusieurs projets dans le territoire nationale et internationale avec un arsenal de 20
sociétés d’exploitation minières qui opère actuellement dans 12 mines localisées en Afrique. Avec des
projets aurifère au Soudan et à la Guinée, le groupe vise à se positionner dans le marché de l’or et à
monter en puissance surtout avec l’aboutissement de son projet d’exploitation d’un gisement d’or au
Soudan et aussi la mise en évidence d’importante réserve en Guinée. Le groupe vise aussi de renforcer
sa postion dans le marché du Cobalt valorisé par la mise en developpement de la mine Pumpi de Cobalt
dans la République démocratique du Gongo.

Figure 1:Lieux d’activité du groupe Managem.

9
 2. Structure et Activité :
Managem regroupe plusieurs filiales d’exploitation et de service qui s’installe dans plusieurs pays
africains. Au Maroc, il dispose de 9 compagnies d’exploitations minière qui contribuent à la production de
4 métaux de base : Cuivre, Cobalt, Zinc, Plomb et un métal précieux : Argent. Ces mines sont distribuées
sur le territoire comme suit :

Figure 2: Les mines de Managem au Maroc.

Certains produits de ces mines sont traités par le service hydrométallurgie du groupe représenté par le
complexe hydrométallurgique de Guemassa et la compagnie Tifnout Tighanimine, Ces centres gèrent
plusieurs unités opérationnelles qui se divise selon deux filières : La première filière se concentre sur la
valorisation du Cobalt afin de produire le cathode du cobalt qui se considère comme un métal critique au
niveau internationale, la deuxième filière est destinée à la production des dérivées de Zinc notamment
l’oxyde de Zinc.

Afin d’accompagner son modèle de développement, le groupe dispose de 2 filiales de service : TECHSUB
et REMINEX. La première filiale se concentre sur les sondages et les travaux souterraines, REMINEX à son
tour divise son service selon 3 axes qui représente les phases de développements d’un projet minière en
commençant par l’exploration géologique jusqu’à l’accompagnement de la mise en œuvre de
l’exploitation avec le centre d’ingénierie des mines et carrières. Cette filiale dispose aussi d’un laboratoire
de recherche pour faire les essaies et les analyses physico-chimiques des différents produits du groupe.

10
 II. Présentation de REMINEX :
1. Introduction
La filiale de recherches et d’ingénieries se considère comme un acteur primordial dans le plan de
développement de Managem. Fondée en 1983, elle démarre son aventure par une première réalisation
industrielle : La mise en place d’un procédé de valorisation du cuivre oxydé de la mine de Bleida. En 1997,
le centre de recherche de REMINEX a posé le procédé industriel de la valorisation du cobalt en cathode
de cobalt de haute pureté, c’est l’invention qui poserait Managem sur un marché de haute criticité vu
l’utilisation de ce métal dans l’industrie hybride. En parallèle, plusieurs procédés industriels ont vu la
lumière notamment pour le Zinc, le Cuivre et le fer.

En étant le mélange de plusieurs profils, la filiale donne une forte importance au travail d’équipe pour
atteindre les objectifs fixés en mettant leurs compétences en commun, cette variété de profils donne
naissance à un collectifs de vocation professionnels qui peut intervenir selon le besoin sur plusieurs
stades. Le partage de connaissance fait partie du travail de REMINEX que ce soit avec les parties prenantes
du groupe notamment les collaborateurs des sociétés d’exploitation ou avec des séminaires organisés
avec les autres sociétés externes qu’exercent une activité minière ou en faveur des étudiants, ce partage
constitue un levé de développement pour le continent et un moyen pour REMINEX d’être à jour de toutes
nouveautés technologiques.

2. Secteur d’activité
L’activité de REMINEX s’étale sur tous les stades de développement d’un projet minier, en commençant
par l’exploration stratégiques dans les zones métallogéniques du Maroc et de l’Afriques pour trouver de
nouveaux gisements. En collaboration avec les sociétés de services miniers, REMINEX se préoccupe des
études de faisabilités des gisements potentiels en menant des études économiques er financières pour
juger sur la rentabilité économique ainsi que la faisabilité technique en prenant en considération des
facteurs de risques de marché. En outre, les ingénieurs miniers exploitent leurs expertises ainsi que leurs
machines pour optimiser la structure interne de la mine de telle façon à faire sortir le maximum de la
richesse possible. Que ce soit une mine à ciel ouvert ou une mine souterraine, la planification est donc
nécessaire pour bien définir les objectifs à court terme et à long terme. C’est le laboratoire qui prend le
relai afin de développer des procédés industriels pour valoriser le minerai extrait, et cela grâce à
l’expérience cumulé pendant des années et aussi à l’utilisation de récente technologies. Avant d’entamer
la production en pleine régime, REMINEX offre un essai de taille semi-industrielle afin de valider les
différents paramètres déterminés dans la phase d’étude du projet. En gros, la filiale intervient dans toutes
les phases d’études du projet jusqu’à la mise en service de l’unité avec des formations organisées en
faveurs des personnels.

11
 III. Présentation du gisement :
1. Situation Géographique :
Ce travail d’option prend comme sujet le gisement cuprifère de Tizert qui se situe sur la bordure nord de
la boutonnière précambrienne d’Ighrem qui appartient à la région de Taroudant qui se trouve à 80 Km à
l’Est d’Agadir, et environ à 250km au sud de Marrakech, il est bordé au sud par une succession des petites
dépressions er de vallées jusqu’à Tata. A l’ouest, il est bordé par l’Adrar de Lkest et au nord par les plateaux
qui donnent sur la plaine du Souss, il se continue vers l’est par les Idrarenes des Inda-ou-Nidif et des Inda-
ou-Zal. La figure suivante illustre la situation géographique de la région :

Figure 3: Situation géographique de la zone du projet Tizert.

2. Contexte Géologique :
L’Anti-Atlas constitue un segment de la chaine panafricaine où les terrains précambriens affleurent dans
les boutonnières disposées suivant une direction globale NE-SW par effet de l’érosion combinée à des
effets d’évènement tectoniques : compressifs, extensionnels et transcurrents. La boutonnière d’Ighrem
est constitué d’un socle Paléoprotérozoïque et d’une couverture Néoprotérozoïque. Selon le découpage
lithostratigraphique adopté dans l’Anti-Atlas en deux groupes : le supergroupe de l’Anti-Atlas et le
supergroupe de Ouarzazate. Le gisement présente une minéralisation stratiforme synsédimentaire avec
des niveaux minéralisés situés sur la série de base infracambrienne.

Figure 4: Coupe géologique du gisement.

12
 3. Synthèse lithologique :
Le gisement de Tizert comporte une minéralisation sédimentaire transgressive composé du haut vers le
bas :

➢ Calcaire rosâtre de l’adoudounien.

➢ Dolomite de Tamjout.

➢ Série de base composée de bancs gréseux de plusieurs dizaines de mètres et de

microconglomérats.

➢ « Petit calcaire » : bancs carbonatés (dolomie microcristalline) d’une dizaine de mètre de

puissance.

➢ Conglomérats verdâtres de base : en discordance avec le précambrien terminal, leur épaisseur

est inférieure à 10m.

➢ Conglomérats à éléments grossiers du précambrien terminal.

La minéralisation est composée de malachite et de chalcosite, ils sont imprégnés dans les fractures ou se
présentent sous forme de lentilles portés par les dolomites de Tamjout ou par les microconglomérats de
la série de base. Sous l’effet de la tectonique, la structure minéralisée adopte la forme d’un paléo relief
avec des épaisseurs variables et deux sens de pendage.

Malachite Cu(OH)₂ Chalcosite Cu₂S

4. Contexte hydrogéologique :
La région de Tizert est bien localisée dans le sud du Maroc, c’est une zone caractérisée par la rareté des
ressources en eaux souterraines et la faiblesse des précipitations. Elle est tectoniquement polyphasée et
à grande variété lithologique. En effet, le contexte hydrogéologique semble complexe à caractériser.

13
 Chapitre II : Choix de la méthode d’exploitation
Avant d’entamer toutes sorte d’étude concernant la méthode d’exploitation, il est indispensable de
démarrer par une étude statistique des paramètres morphologiques notamment la puissance et le
pendage afin de déterminer les paramètres caractéristiques pouvant conditionnés les limites techniques
de méthodes d’exploitation.

I. Etude statistique de l’extension :

1. Etude statistique sur le pendage :
Une étude statistique de pendage est menée sur l’extension NW en utilisant les données de la base de
données fournie sous forme de block model, la zone NW est constituée de deux zones superposées : La
zone 09 et la zone 35, les deux zones caractérisées par un pendage faible au milieu de la structure et qui
commence à s’accroitre en se rapprochant des bords. En générale, le pendage reste faible pour l’ensemble
du gisement et suivant la direction SW-NE avec quelques fluctuations de l’ordre de 60° à 75° vers les bords.
Les figures 5 et 6 représentent la cartographie du pendage pour les zones 35 et 09

Figure 5: Cartographie du pendage de la zone 09.

14
 Figure 6: Cartographie du pendage de la zone 35.

Figure 7: Boîte à moustaches du pendage.

La figure 7 illustrent la distribution des valeurs du pendage, Les boites illustrent la marge de distribution
des valeurs qui regroupent 50% des valeurs, donc 75% des valeurs sont inférieurs à 25°, on disait alors
que le gisement adopte une forme aplatie. La zone 09 présente un écart-type important en comparaison
avec celui de la zone 35. Les deux parties rencontrent des valeurs aberrantes qui sont généralement les
valeurs situées à la limite de la structure.

Figure 8: Histogramme du pendage.

On remarque la présence de 3 familles de valeurs dans la figure 8 avec une densité plus grande pour la
famille qui varie aux alentours de 20°, les deux autres familles adoptent une forme presque similaire.

15
 2. Etude statistique sur la puissance :
La puissance verticale constitue une vraie importance dans le processus de choix de la méthode
d’exploitation, une analyse statistique est donc nécessaire. Les deux figures suivantes 9 et 10 représentent
la cartographie de la puissance pour chaque zone, on remarque qu’un pourcentage important du
gisement adopte une puissance inférieure à 6m, au milieu la structure commence petit à petit d’avoir des
puissances entre 6 et 12 m avant d’atteindre la tranche maximale avec des puissance supérieur à 12m.
On remarque une brusque variation de puissance pour la zone 09 où on passe d’une puissance supérieure
à 12m à une puissance inférieure à 6m, ceci peut être expliqué par la présence d’un artefact due à une
erreur de modélisation. La zone 35 n’adopte pas la tranche des grandes puissances (Puissance supérieur
à 12m) avec la présence d’une petite tache de puissance entre 6m et 12m.

Figure 9: Cartographie de la puissance de la zone 09.

Figure 10: Cartographie de la puissance de la zone 35.

16
Les dernières cartographies illustrent en 2D la répartition spatiale de la puissance, mais cela ne suffit pas
pour avoir une idée de la répartition 3D de la puissance. Par conséquent, on fera l’histogramme de la
puissance par zone pour voir la densité de présence de chaque tranche.

Figure 11: Cartographie de la puissance de la zone 35.

Pour la zone 09, les formes des tranches des moyennes et grandes puissances se ressemblent vue que la
densité est maximale aux bords gauches des intervalles, alors que la densité de la tranche des petites
puissances est maximale juste au milieu de l’intervalle notamment aux alentours de 4m.

Figure 12: Boîte à moustaches de la puissance.

La figure 11 nous permet de dire que la zone 35 est caractérisé par une faible variabilité de la puissance
en comparaison par celle de la zone 09, avec la non-présence des valeurs aberrantes vu que la
distribution est plus ou moins continue.

17
Pour Synthétiser, l’extension NW présente une minéralisation quasi-plane avec une valeur moyenne du
pendage qui ne dépasse pas 25°, la cartographie du pendage illustre la présence des valeurs aberrante qui
sont peut-être dû à un artefact de modélisation. L’analyse statistique de la puissance a montré que le
gisement adopte une puissance verticale moyenne de l’ordre de 6m, sans oublier de citer qu’il y avait un
changement brusque de la puissance au centre de la minéralisation qui est dû peut-être à un accident
géologique.

II. Etude géotechnique :

Dans le moment actuel aucune étude géotechnique ni sondage n’est fait pour l’extension NW, Ces
données sont primordiales pour la détermination de la résistance des matériaux constituant le gisement.
Pour résoudre cela, on considère que la zone NW est une continuité du gisement principale, donc on
prend le même résultat des sondages déjà fait pour le gisement principal, ces résultats sont présentés par
la valeur du RMR, une classification des massifs rocheux pour estimer la résistance à la compression, qui
peut varier entre 0 et 100.

Tableau 1: RMR par domaine lithologiques

Domaine Min Mean Max Ecart-type
Dolomite 25 51 77 11
Série de base 26 62 83 8,6
Conglomérat 27 65 75 6,6

Le système RMR propose une classification pour évaluer le comportement du massif rocheux, cette
classification possède 5 intervalles :

Tableau 2: Classification du système RMR

Classification RMR
Très mauvaise qualité 0 – 20
Mauvaise qualité 20 – 40
Qualité moyenne 40 – 60
Bonne Qualité 60 – 80
Roche intacte 80 – 100

La minéralisation se trouve dans les domaines lithologiques Conglomérat et série de base, ils sont classés
alors comme des massifs rocheux de moyenne qualité et le conglomérat est classé comme une roche de
bonne qualité.

18
 III. Sélection de la méthode d’exploitation :
1. Introduction :
Le choix de la méthode d’exploitation constitue un problème compliqué au niveau de l’exploitation
minière, différents critères rendent difficile la tâche de choisir celle la plus adapté. Pour résoudre cette
problématique on fait appel à la méthode AHP (analyse multicritère hiérarchique), il s’agit d’une
technique de prise de décision multi-attributs, développée pour la première fois en 1980 par Thomas
Saaty. C’est un outil permettant de combiner des facteurs qualitatifs et quantitatifs dans la sélection d’un
processus.

Il est utilisé pour définir les priorités dans une situation problématique complexe et multicritère. Son
principal avantage est sa capacité à traiter des problèmes complexes et mal structurés, qui ne peuvent
généralement être manipulés par des modèles mathématiques rigoureux

Pour notre cas, il s’agit de comparer une paire d’alternative, en ce qui nous concerne les méthodes
d’exploitation, par rapport à chaque critère fixé (paramètres caractéristiques du gisement) afin de donner
un poids relatif qu’on va l’utiliser pour calculer un score en le multipliant par la valeur intrinsèque du
facteur morphologique qui représente le critère.

Les critères adoptés sont cités dans le tableau ci-après :

Tableau 3: L’importance des différents critères de choix.

Critère Importance
Fondamentale Haute Très haute
Pendage *
Puissance *

Pour arriver à déterminer le poids relatif de chaque critère par rapport à un autre, il est à déterminer la
matrice de comparaison par paire, c’est ce qui nous permettra de calculer la priorité de chaque méthode
pour chaque partie du minerai.

Les principales méthodes d’exploitation les plus adaptés au gisement, en vue des résultats de l’analyse
statistiques des paramètres caractéristiques de ce dernier, sont :

19
On compare les méthodes citées précédemment par rapport à chaque critère mentionné dans le tableau
2, et cela suivant l’échelle de Saaty qui varie entre 1 et 9 :

Tableau 4: Echelle de Saaty.

Poids Le jugement
1 Même préférence
3 Préférence modérée
5 Forte Préférence
7 Très forte Préférence
9 Extrême Préférence
2, 4, 6, 8 Jugements intermédiaires par rapport aux adjacents

La matrice de comparaison par paire et par rapport à un critère donné adopte la forme suivante :

R&P SR&P PR&P SNA C&F LT

R&P 1 . . . . .
SR&P . 1 . . . .
PR&P . . 1 . . .
SNA . . . 1 . .
C&F . . . . 1 .
LT ( . . . . . 1)

Chaque coefficient de la matrice indique le poids avec lequel on favorise la méthode d’exploitation
correspondante à celle de la ligne par rapport à la méthode correspondante à celle de la colonne.

Il reste qu’à vérifier la cohérence de ces coefficients en calculant le ratio de consistance, le ratio se
décompose d’un indice de consistance et un ratio randomisé qui dépend de la taille de la matrice de
comparaison :

Tableau 5: Les valeurs du ratio randomisé selon la taille.

Le ratio de consistance ne dépend que de la valeur propre maximale ainsi que le nombre de critères ;
𝜆 max − 𝑛
𝐶𝑅 =
𝑛−1
Si la valeur du ratio reste inférieure à 10% on dit que la matrice est cohérente.

20
 2. Remplissage des matrices :
Le remplissage des matrices est établi en exploitant le résultat d’un questionnaire distribué pour les
experts de REMINEX, on introduit des intervalles de définitions de pendage et de puissance pour chaque
matrice vue que le facteur de préférence d’une méthode par rapport à une autre dépend de la valeur de
ces deux paramètres.

On définit alors 3 intervalles de définitions par rapport au pendage et 2 intervalles de définition par
rapport à la puissance :

Chaque intervalle adopte une matrice de comparaison qu’on va exploiter pour déterminer les poids
relatifs par le calcul du vecteur propre correspondant à la valeur propre maximale.

Exemple :

Dans le cas où la puissance est inférieure à 20m, la matrice de comparaison prendra la forme suivante :

R&P SR&P PR&P SNA C&F LT

R&P 1 3 3 7 5 1
SR&P 1/3 1 1/2 7 5 2
PR&P 1/2 2 1 7 5 2
SNA 1/7 1/7 1/7 1 1/5 1/4
C&F 1/5 1/5 1/5 5 1 1/3
LT ( 1 1/2 1/2 4 3 1 )

21
Après calcul du vecteur propre de la valeur propre maximale en utilisant le langage de programmation R,
on trouve les poids relatifs des méthodes par rapport à une puissance inférieure à 20m :
R&P 25%
SR&P 19%
PR&P 31%
SNA 2%
C&F 6%
LT (14%)

Le ratio de consistance est bien inférieur à 10% donc la matrice est cohérente.

Le vecteur des poids relatif indique le degré de préférence de chaque méthode d’exploitation si la
puissance est inférieure à 20m.

De même on fera pour toutes les matrices de comparaison pour en tirer les poids relatifs (voir annexe)

3. Calcul du score :
Après avoir calculer les poids relatifs de chaque méthode pour les différents intervalles de définition, on
calcule le score final pour chaque block de taille (1.25x1.25xPuissance verticale) par la formule suivante :

𝑆𝑐𝑜𝑟𝑒(𝑀é𝑡ℎ𝑜𝑑𝑒) = 𝑃𝑜𝑖𝑑𝑖𝑛𝑡𝑒𝑟𝑣𝑎𝑙𝑙𝑒 𝑑𝑒 𝑙𝑎 𝑝𝑢𝑖𝑠𝑠𝑎𝑛𝑐𝑒 ∗ 𝑃𝑢𝑖𝑠𝑠𝑎𝑛𝑐𝑒𝑛𝑜𝑟𝑚𝑎𝑙𝑖𝑠é𝑒 +

𝑃𝑜𝑖𝑑𝑖𝑛𝑡𝑒𝑟𝑣𝑎𝑙𝑙𝑒 𝑑𝑢 𝑝𝑒𝑛𝑑𝑎𝑔𝑒 ∗ 𝑃𝑒𝑛𝑑𝑎𝑔𝑒𝑛𝑜𝑟𝑚𝑎𝑙𝑖𝑠é

On fait la normalisation des valeurs de puissance et de pendage pour garder le degré d’importance de
chaque facteur morphologique. Une fois le score est calculé, on affecte la méthode obtenant le score
maximal au bloc correspondant.

22
 4. Résultat et synthèse :

Figure 13: La répartition des méthodes d’exploitation de la zone 09.

Figure 14: La répartition des méthodes d’exploitation de la zone 09.

Le système présente deux cas de figures correspondant aux chaque zones de la structure, on trouve 4
méthodes d’exploitations pour la zone 09 et 3 méthodes d’exploitation pour la zone 35. La méthode des
chambres et piliers couvrent le plus grand pourcentage de la superficie vu la présence des faibles
pendages et la puissance reste insuffisante pour faire l’exploitation en sous-niveaux abattus. Pour les
zones aux bords, l’algorithme propose de continuer l’exploitation en utilisant le « Cut and Fill », cela est
justifié par l’augmentation du pendage. Le sous-niveaux abattu est utilisé principalement dans les zones
où la puissance dépasse 20m.

23
 5. Description de la méthode d’exploitation retenue :
La méthode des chambres et piliers (Room and Pillar) est celle retenue par l’étude de choix pour
l’exploitation en souterraine de l’extension NW. La méthode prend sa forme classique au milieu de la
minéralisation où le pendage reste faible pour permettre aux engins de rouler, on utilise une rampe dans
les zones caractérisées par un pendage élevé pour vaincre l’angle de roulage des engins fixée à 10 degrés
(voir figure 14). L’exploitation sera par étage avec remblayage dans les zones où la puissance dépasse la
hauteur maximale 3,5m.

Figure 15: Schéma de la méthode Step Room and Pillar.

6. Application Shiny :
Pour faciliter l’approche de la sélection de la méthode d’exploitation, une application Shiny à base du
langage de programmation R est établie, cela permettra aux experts de REMINEX de gagner en temps vue
qu’elle fournit une cartographie en 3D des méthodes d’exploitations convenables à la morphologie du
gisement. Cette application peut être hébergé sur WEB, pour profiter de la capacité de calcul des serveurs,
ou utilisé localement sur les machines de REMINEX.

24
 i. Interface de l’application :

L’interface de l’application possède différentes entrés, On trouve tout en haut le bouton pour faire entrer
la base de données qui est le bloc model sous format CSV, une fois on clique sur le bouton «charger les
données » les listes déroulantes s’actualise pour permettre à l’utilisateur de choisir ses coordonnées X,
Y,Z et la zone d’étude, ces données sont utilisées à la fois pour déterminer la méthode d’exploitation et
pour préparer le graphe en 3D qui montrera les résultats en cliquant sur le bouton « Lancer les calculs »

25
 ii. Résultat :
Le programme affiche la forme de la zone sélectionnée, Les couleurs du graphe représentent les méthodes retenues par l’algorithme AHP. Ce
résultat dépend principalement de la modélisation du gisement, pour cette raison on doit assurer que le bloc model ne contient aucun artefact
et que la modélisation représente correctement le contexte morphologique du gisement.

26
 Chapitre II : Optimisation des ressources exploitables
L’optimisation des ressources nécessitent la prise en considération de plusieurs paramètres techniques
et économiques pour but de déterminer les réserves disponibles.

I. Méthodologie :
Pour déterminer les parties du gisement qui constituera les réserves à partir des ressources potentielles
on utilise le module MSO (Mineable Shape Optimizer), L’utilisation d’un tel outil permet de concevoir la
forme et la position géométrique des ouvertures en prenant en considération la géométrie des tailles,
l’ouverture minimale minière, la densité de stérile, la valeur de coupure en NSR (Net Smelter Return), la
dilution de contact et le pendage.

Figure 16: Les zones exploitables avec la méthode chambre et pilier.

27
 II. Paramètres économiques:
Pour déterminer la rentabilité économique il faut calculer la valeur du NSR de chaque bloc minéralisé, ce
NSR représente le restant de la valeur économique après traitement d’une tonne du minerai sans tenir en
compte le coût minier, la charge de structure, etc. Ceci servira à déterminer un seuil à partir duquel
l’exploitation minière du bloc devient rentable.

D’après le chapitre précèdent, la méthode d’exploitation retenue est les chambres et piliers, le calcul de
la valeur de coupure est composé alors de 17,65 $/tonne pour le coût opérationnel, 10,85$/tonne pour
le traitement du minerai et 2,08$/tonne pour les charges de structures, ceci donne une valeur de coupure
de 30,11$/tonne. En fait toute exploitation d’une chambre d’exploitation avec une valeur NSR inférieure
à cette valeur de coupure est considérée non rentable et ne sera pas compté en réserve.

Il reste maintenant à calculer la valeur NSR d’un bloc minéralisé en utilisant les paramètres suivants :

Tableau 6: Paramètres économiques du NSR.

Données de calcul du NSR
Teneur équivalente Variable
Teneur payé 96% * Teneur équivalente
Cours Cuivre ($/TM) 6000
Prix Cuivre du concentré ($/TC) Cours Cuivre * Teneur Payée
Frais traitement du concentré 140
($/TC)
Frais de raffinage du concentré 69.4
($/TC)
Frais d’approche ($/TC) 36.22

TM (*) = Tonne de métal

TC (**)= Tonne de concentré

Pour calculer la valeur NSR de chaque bloc on utilise la formule suivante :

𝑁𝑆𝑅𝑏𝑙𝑜𝑐 = 𝑇𝑒𝑛𝑒𝑢𝑟é𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡𝑒 𝑏𝑙𝑜𝑐 ∗ 𝑃𝑜𝑢𝑟𝑒𝑐𝑒𝑛𝑡𝑎𝑔𝑒. 𝑃𝑎𝑦é ∗ 𝐶𝑜𝑢𝑟𝑠𝐶𝑢𝑖𝑣𝑟𝑒

− (𝐹𝑟𝑎𝑖𝑠𝑇𝑟𝑎𝑖𝑡𝑒𝑚𝑒𝑛𝑡 + 𝐹𝑟𝑎𝑖𝑠𝑅𝑎𝑓𝑓𝑖𝑛𝑎𝑔𝑒 + 𝐹𝑟𝑎𝑖𝑠𝐴𝑝𝑝𝑟𝑜𝑐ℎ𝑒 )

Après calcul, il est possible de cartographier le NSR sur l’ensemble des lentilles probable à exploiter avec
la méthode d’exploitation retenue :

28
Figure 17: Répartition des ressources suivant la valeur du NSR de la zone 09.

Tableau 7: Répartition des ressources suivant le NSR de la zone 09.

29
Figure 18: Répartition des ressources suivant la valeur du NSR de la zone 35.

Tableau 8: Répartition des ressources suivant le NSR de la zone 35.

30
III. Paramètres techniques :
1. La géométrie des tailles :
Afin d’optimiser les ressources pour tirer le tonnage économiquement rentable, on définit premièrement
le dimensionnement de la méthode d’exploitation retenue. L’extension NW se considère comme une
continuité de la partie NORD, alors on adopte la même étude conceptuelle géotechnique déjà fait pour la
partie NORD.

Il est clair que le stresse devient de plus en plus important quand on gagne de la profondeur, Une variation
de dimensionnement des piliers de l’exploitation est donc nécessaire pour rendre les ouvertures plus
stables. D’après l’étude géotechnique fait par REMINEX en collaboration avec Met-Chem, on retient les
valeurs de la largeur et de la longueur des piliers :

Figure 19: Dimensionnement des piliers en fonction de la profondeur.

Après analyse de la profondeur du gisement, il apparait que la profondeur de l’extension varie entre 220m
et 816 m Pour la zone 09 tandis que la profondeur de la zone 35 varie entre 280 m et 455 m. Alors les
piliers de la zone 09 prennent une longueur et largeur de 16m au tout début après plus on descend ils
commencent à avoir des tailles plus importantes de l’ordre de 31 m. Le dimensionnement des piliers de
la zone 35 prend des valeurs moyenne entre 16m et 22m pour la largeur et la longueur.
La hauteur des tailles reste à 3,5m quand la puissance est inférieure à cette valeur, quand la puissance
dépasse cette limite il faut remblayer les chambres pour les confiner afin d’augmenter la résistance.

31
Le tableau suivant résume la largeur et la longueur des tailles utilisées pour l’optimisation MSO :

Tableau 9: Paramètres géométriques des tailles.

Profondeur (m)
Méthode d’exploitation 200 - 300 300 – 400 400 - 500 500 - 600 600 - 800
Taille (m)
Pilier X = 16 X = 20 X = 22 X = 27 X = 31
Y = 16 Y = 20 Y = 22 Y = 27 Y = 31
Chambre X=7 X=7 X=7 X=7 X=7
Y=6 Y=6 Y=6 Y=6 Y=6

L'exploitation minière récupère rarement toutes les ressources présentes dans un gisement de minerai.
Une certaine quantité de stérile est généralement mélangée au minerai pendant l'extraction, ces stériles
mélangés au minerai sont appelés dilution et sont généralement exprimés en facteur de dilution (%). Elle
provient principalement des épontes de la minéralisation toit et mure, cela dit que le taux de dilution
dépend de la méthode d’exploitation utilisée. Pour la méthode des chambres et piliers, on trouve
généralement un taux de dilution de 5%.
Pour la présente méthode d’exploitation on distingue de type de dilution; dilution des chambres primaires
et des chambres secondaires :

Figure 20: Chambre et pilier.

32
En ce qui concernent l’ouverture des chambres primaires, on envisage 5 cas d’exploitation possible :

• Cas 1 : le traçage avec une récupération maximale de minerai

• Cas 2 : le traçage des chambres visant la prise de la moitié de la quantité de stérile du premier cas

• Cas 3 : traçage dans le minerai seul sans toucher le stérile (dilution = 0%)

• Cas 4 : traçage avec sélectivité au niveau du toit de la minéralisation

• Cas MSO : le cas utilisé dans l’optimisation des ressources sur MSO

Figure 21: Différentes cas des ouvertures.

La forme des chambres primaires génère une problématique entre le taux de dilution et le taux de
récupération du minerai, les Cas 1, Cas2 et Cas 4 proposent un bon taux de récupération du minerai mais
avec une dilution importante alors que le cas 3 est caractérisé par un taux de dilution nul avec une
récupération beaucoup moine bonne que les autres cas. En effet, le Cas MSO a été introduit pour trouver
un arrangement qui minimise le facteur de dilution avec une bonne récupération du minerai. Dans ce cas
on prend une dilution de 7% en moyenne pour la globalité du gisement.

33
Les chambres secondaires n’ont pas le même impact sur la dilution globale de l’exploitation vu qu’elle se
limite à des longueurs faibles (longueurs des piliers) ce qui rend possible le contrôle des limites de ces
chambres en suivant le pendage et la puissance local de la structure.
2. Récupération minière :
L’exploitation par chambre et pilier adopte en moyenne un taux de défruitement de 85% sans tenir en
compte la récupération plus tard des piliers. Pour notre cas, la géométrie des tailles dicte à ce que le taux
de défruitement commence de 83% et continue à diminuer avec la profondeur. La récupération des piliers
ne se tient pas en considération pour l’analyse actuelle mais cela n’exclue pas une récupération plus tard
des piliers dans la durée de vie de la mine

IV. Résultat d’optimisation :

En tenant en compte tous les paramètres cités précédemment et utilisés pour achever l’optimisation avec
le module MSO, on réalise l’optimisation de chaque zone de l’extension NW.

Le processus prend comme résultat les tailles qui regroupent un ensemble de blocs ayant une valeur NSR
supérieure à la teneur de coupure en NSR qui vaut 30$/T. Ces tailles seront examinées en prenant en
considération les paramètres géomécanique pour en tirer la sélection finale.

Les résultats de l’optimisation présentent une réserve de 12 639 871 Tonnes du minerai avec une teneur
moyenne de 0.92% de cuivre et 17,2 Agppm.

Tonnes(T) Densité NSR($/T) Cu OK(%) Agppm Cueq

Réserve 12 639 871 2,72 60,81 0,92 17,20 1,09

Les deux figures suivantes présentent la forme spatiale des tailles après optimisation :

34
Figure 22: Résultat de l’optimisation MSO (Vue en plan).

Figure 23: Résultat de l’optimisation MSO (Vue en 3D).

35
 V. Infrastructure minière :
1. Préparation des réserves :
Le résultat de l’optimisation présente des zones non exploitables à cause de leurs positionnements par
rapport à l’ensemble des tailles ou de ses géométries critiques. Cela est due de la traitante indépendante
des ouvertures par MSO sans avoir recours aux voisinages :

• Suppression des blocks qui ne justifie pas l’investissement de l’infrastructure

• Les chambres auxquelles est impossible d’accéder

La réserve finale présente la forme suivante :

Figure 24: Résultat de la filtration des blocs.

36
 2. Description de l’Infrastructure :
La conception de l’infrastructure minière est très importante dans le développement d’un projet minier,
elle est indispensable pour le calcul du CAPEX, le calcul de la distance de roulage et la planification. Le
sujet actuel se concentre sur la planification longue terme de la mine souterraine, l’infrastructure doit
donc en premier lieu assurer l’accès à la minéralisation ainsi que le déplacement des mineurs et des engins
tout en respectant des normes d’aérages, elle doit aussi contenir un système d’acheminement du minerai
du fond de la mine à la surface.

La structure du gisement favorise une voix de roulage d’extraction du minerai en dessous de la

minéralisation tandis que l’accès à la minéralisation est assuré par un puit qui sera positionner dans les
alentours du barycentre de la minéralisation pour minimiser les coûts du transport. En dessus de cette
voix de roulage de 50m on construira une galerie principale permettant le déplacement au sein de la
minéralisation avec des virages et pente douce pour respecter l’angle de roulage des engins fixée à 15%.
Le transport du minerai se fait par des cheminées liant la galerie principale à la galerie d’extraction par
des cheminées collectrices en exploitant l’effet de la gravité. Le système d’extraction s’occupera du
transport du minerai.

Des ouvrages pour réaliser des ateliers de maintenance et des bassins pour le système d’exhaure sont si
importants pour l’infrastructure. Ci-après une description des ouvrages principales :

Puit :
Un puit reliera la surface à la minéralisation pour permettre l’accès des mineurs et aussi le transport du
minerai au travers la salle de chargement.
Galerie d’exploitation :
Il s’agit d’une galerie de 2 386m avec un angle maximal de 15%, elle assure l’accès aux différentes zones
de minéralisation.
Galerie d’extraction :
La galerie est située 50m en dessous de la galerie d’exploitation en ayant des cheminées de récupération
du minerai tous les 150m.
Cheminées d’aérage :
Les cheminées d’aérages sont positionnées de telle sorte à assurer l’aire fraiche le long du traçage.

37
Figure 25: Vue en 3D de la minéralisation et de l’infrastructure.

Figure 26: Vue en 3D de la minéralisation et de l’infrastructure.

38
 Chapitre III : Planification de la Production
Après avoir déterminé les ressources optimisées en fixant les paramètres techniques et économiques
nécessaire, on arrive à l’étape de la planification où on fixe les objectifs à atteindre et les moyennes
nécessaires ainsi que le séquençage de la production. Ce séquençage joue un rôle important pour les
recettes long terme de la production, c’est pour cela qu’on utilisera une optimisation de la marge
opérationnelle par une variable décisionnelle désignant le séquençage de la production.

I. Découpage en panneaux d’exploitation :

L’extension NW du gisement de Tizert présente 11,12 Mt de réserve, la stratégie de l’entreprise vise à
exploiter ces réserves dans une durée de vie de 10 ans, ce qui donnera un rythme de production moyen
de 1,11 Mt/an. On découpera alors le gisement en panneaux d’exploitations de telle sorte à l’exploiter
pendant une année, ce qui donnera 10 Panneaux. Ces panneaux sont formés par les blocs qui ont des
propriétés morphologiques similaire en assurant à ce que le maximum nombre de panneaux sont
accessible par la galerie d’exploitation.

Le résultat du découpage donne 2 scénarios de découpage, on utilisera ces scénarios ci-après pour
déterminer la valeur actuelle nette de la marge opérationnelle pour les évaluer.

Figure 27: Scénario 1 du découpage. Figure 28: Scénario 2 du découpage.

39
 II. Optimisation du séquençage de la production:
1. Introduction :
La durée de vie du projet est supposée prendre une durée de 10 ans, le long de cette durée la valeur
actuelle de chaque panneau diminuera une année après l’autre par effet de l’actualisation, on définit cette
valeur actuelle nette par la formule suivante :
𝑉𝑎𝑙𝑒𝑢𝑟𝑎𝑛𝑛é𝑒 0
𝑉𝑎𝑙𝑒𝑢𝑟𝑎𝑛𝑛é𝑒 𝑖 =
(1 + 𝑎)𝑖
Où la constante a constitué le taux de l’actualisation, qu’on prendra égale à 10%

On exploite les valeurs du NSR calculées précédemment pour chaque bloc pour calculer la marge
opérationnelle de l’exploitation de chaque panneau.

Tableau 10: Marge opérationnelle actuelle des blocs pour le scénario 1.

Marge Opérationnelle*
Block Tonnage (Mt)
(M$)

P3-1 1,5 50,8

P3-2 1,43 52,1

P4-11 1,32 38,8

P4-21 1,24 25,6

P4-12 1,46 39,2

P4-22 1,32 38

P1-1 1,3 64,3

P1-21 1,23 50,3

P1-22 1,26 31,3

P1-23 1,12 20,15

(*) = Marge opérationnelle de la production hors coût de transport

40
 Tableau 11: Marge opérationnelle actuelle des blocs pour le scénario 2.

Marge Opérationnelle*
Block Tonnage (Mt)
(M$)

P3-1 1,5 50,8

P3-2 1,43 52,1

P4-11 1,32 38,8

P4-21 1,24 25,6

P4-12 1,46 39,2

P4-22 1,32 38

P1-1 1,3 64,3

P1-21 1,23 35,3

P1-22 1,26 29,8

P1-23 1,12 36,8

(*) = Marge opérationnelle de la production hors coût de transport

2. Formulation mathématique :
Comme mentionnée précédemment, l’objectif de cette étude est de maximiser la marge opérationnelle
en suivant un séquençage optimal, cela se traduira mathématiquement en un problème d’optimisation.
Ce problème est constitué d’une fonction objectif et des contraintes à respecter le long du processus, ces
contraintes sont comme suit :

• Contrainte de continuité : Le séquençage de la production doit avoir une forme continue

• Contrainte de capacité et ressources : Un seul block ne doit pas être exploité plus qu’une fois, et
on ne peut pas exploiter plus qu’un block dans une période

• Contrainte de démarrage: Il faut garantir une production le plus tôt possible

• Contrainte d’accessibilité : Avant d’exploiter un block, il faut assurer l’accessibilité spatiale

41
Afin de traduire ces écriture en une forme mathématique, on modélise tout d’abord le gisement par une
forme matricielle qui présente des différentes liaisons d’accès aux panneaux :

Figure 29: Modélisation du scénario 1.

Figure 30: Modélisation du scénario 2.

La ligne verte, reliant les panneaux, représente la possibilité d’accès d’un panneaux par l’autre, ces
représentation seront utilisées par la suite pour déterminer la forme de la contrainte géométrique

Le séquençage des panneaux est défini par la variable binaire 𝑥𝑖𝑗𝑡 qui prendra la valeur vrai si le
panneau, situé sur la ligne 𝑖 et la colonne 𝑗, est exploité pendant la période 𝑡.

42
La formulation mathématique du problème d’optimisation prendra la forme suivante :

𝑉𝑖𝑗0
Max ∑𝑖 ∑𝑗 ∑𝑡 ∗ 𝑥𝑖𝑗𝑡
(1+𝑑)𝑡

Contrainte des ressources : ∑𝑡 𝑥𝑖𝑡 == 1 ∀ (𝑖, 𝑗) ∈ [1,10]𝑥[1,10]

Contrainte de la capacité ∑𝑖 ∑𝑗 𝑥𝑖𝑗𝑡 == 1 ∀ 𝑡 ∈ [1,10]

Contrainte d’accessibilité : 𝑥𝑖𝑗𝑡 − ∑𝑙,𝑚∈𝐵𝑖,𝑗 ∑𝑡−1

𝑘=1 𝑥𝑙𝑚𝑘 ≤ 0 ∀ 𝑖, 𝑗, 𝑡 ∈ [1,10]𝑥[1,10]𝑥[1,10]

Contrainte de continuité : 𝑥𝑖𝑗𝑡 − ∑𝑙,𝑚∈𝐵𝑖,𝑗 𝑥𝑙𝑚(𝑡−1) ≤ 0 ∀ 𝑖, 𝑗, 𝑡

Contrainte de démarrage : ∑𝑖,𝑗∈𝐵𝑃 𝑥𝑖𝑗0 == 1

𝐵𝑖 : 𝐿′ 𝑒𝑛𝑠𝑒𝑚𝑏𝑙𝑒 𝑑𝑒𝑠 𝑏𝑙𝑜𝑐𝑘𝑠 𝑒𝑛𝑡𝑜𝑢𝑟𝑎𝑛𝑡 𝑙𝑒 𝑏𝑙𝑜𝑐𝑘 𝑖

𝐵𝑃 : 𝐿′ 𝑒𝑛𝑠𝑒𝑚𝑏𝑙𝑒 𝑑𝑒𝑠 𝑏𝑙𝑜𝑐𝑘𝑠 𝑙𝑒𝑠 𝑝𝑙𝑢𝑠 𝑝𝑟𝑜𝑐ℎ𝑒 𝑎𝑢 𝑝𝑢𝑖𝑡
On remarque que le problème d’optimisation est non linéaire vu la présence de la composante de
l’actualisation ce qui rend la résolution plus compliquée. En fait, pour résoudre ce problème de non-
linéarité on fait stocker la valeur actualisée d’une année à l’autre pour chaque panneau dans une matrice
qu’on va faire appelle sans avoir besoin de faire le calcul au sein de l’équation elle-même.

On définit alors la matrice des valeurs actuelles de la marge opérationnelle de taille 10x10 :

Donc la valeur 𝑉𝑖𝑡 représente la valeur actuelle de la marge opérationnelle du panneau i pendant
l’année t.

43
La formulation mathématique devient linéaire et prendra la forme suivante :

Max ∑𝑖 ∑𝑡 𝑉𝑖𝑡 ∗ 𝑥𝑖𝑡

Contrainte des ressources : ∑𝑡 𝑥𝑖𝑡 == 1 ∀ 𝑖 ∈ [1,10]

Contrainte de la capacité ∑𝑖 𝑥𝑖𝑡 == 1 ∀ 𝑡 ∈ [1,10]

Contrainte d’accessibilité : 𝑥𝑖𝑡 − ∑𝑙∈𝐵𝑖 ∑𝑡−1

𝑘=1 𝑥𝑙𝑘 ≤ 0 ∀ (𝑖, 𝑡) ∈ [1,10]𝑥[1,10]

Contrainte de continuité : 𝑥𝑖𝑡 − ∑𝑙∈𝐵𝑖 𝑥𝑙(𝑡−1) ≤ 0 ∀ (𝑖, 𝑡) ∈ [1,10]𝑥[1,10]

Contrainte de démarrage: ∑𝑖,𝑗∈𝐵𝑃 𝑥𝑖𝑗0 == 1

𝐵𝑖 : 𝐿′ 𝑒𝑛𝑠𝑒𝑚𝑏𝑙𝑒 𝑑𝑒𝑠 𝑏𝑙𝑜𝑐𝑘𝑠 𝑒𝑛𝑡𝑜𝑢𝑟𝑎𝑛𝑡 𝑙𝑒 𝑏𝑙𝑜𝑐𝑘 𝑖

𝐵𝑃 : 𝐿′ 𝑒𝑛𝑠𝑒𝑚𝑏𝑙𝑒 𝑑𝑒𝑠 𝑏𝑙𝑜𝑐𝑘𝑠 𝑙𝑒𝑠 𝑝𝑙𝑢𝑠 𝑝𝑟𝑜𝑐ℎ𝑒 𝑎𝑢 𝑝𝑢𝑖𝑡
Dans notre système il y a que deux contraintes nécessaire à la production : Contrainte de capacité et
Contrainte des ressources ; tout autre contrainte permet à l’entreprise de pousser encore plus la
complexité du problème et définir le mode de fonctionnement dépendamment de sa stratégie. Pour cela,
on construira plusieurs scénarios d’exploitation pour chaque scénario de découpage, ces scénarios
constituent le mode d’exploitation retenue par l’entreprise notamment le démarrage de l’exploitation par
les panneaux les proches au puit ou d’imposer une forme d’exploitation continue.

Exemple :

Prenant l’exemple d’une exploitation qui ne tient en considération que les contraintes indispensable à
l’exploitation minière, ces contraintes sont résumées comme suit :

• Contrainte d’accessibilité

• Contrainte des ressources

• Contrainte de capacité

44
 3. Résultat :
Après avoir codifier les différentes mode d’exploitation pour les deux scénarios (Voir annexe), le logiciel d’optimisation CPLEX STUDIO donne les
résultats optimaux du séquençage :

305 300
300,6 297,08
299,5 295,9
300
295
294,76 290,9
295
Marge Ope (M$)

Marge Ope (M$)

290

290
286,5 284,03
285
285

280
280

275 275
A B C D A B C D
Mode d'exploitation(*) Mode d'exploitation(*)

Figure 31: Résultat d’optimisation pour le scénario de découpage 1 Figure 32: Résultat d’optimisation pour le scénario de découpage 2
Mode d’exploitation(*) :

• A : Une exploitation sans tenir compte de la contrainte de continuité ou de la contrainte de démarrage.

• B : Une exploitation qui démarre par l’exploitation des panneaux les plus proches au puit.

• C : Une exploitation continue dans l’espace d’une année à l’autre.

• D : Une exploitation qui démarre par les panneaux les plus proches au puit et qui adopte une forme continue.

45
 4. Mode d’exploitation retenue :
Durant une activité minière on s’attend à ce que l’avancement de l’extraction se fait d’une manière
continue d’un bloc à l’autre, cette attitude contribue à minimiser les déplacements des machines. En
outre, l’infrastructure est indispensable dans ce genre de travaux notamment les galeries d’extraction et
les cheminées d’aérage, alors il faut tout d’abord favoriser le démarrage de l’exploitation par les panneaux
les plus proches au puit. Pour respecter ces contraintes d’exploitation on retient le mode de
fonctionnement D du premier scénario de découpage. Ce séquençage donne une marge opérationnelle
totale actuelle de 286,5 M$ en démarrant l’exploitation par le panneau P3-1, le reste de l’évolution de
l’exploitation est résumé par le tableau suivant :

Tableau 12: Le séquençage de l’exploitation selon le mode retenue.

Année Année Année Année Année Année Année Année Année Année Année
1 2 3 4 5 6 7 8 9 10
Panneau P3-1 P3-2 P4-21 P4-11 P4-12 P1-1 P1-21 P4-22 P1-22 P1-23
exploité

Ces données serviront à construire un enchainement dans le temps de la production comme suit :

Figure 33: La répartition temporelle de l’exploitation des panneaux.

46
 Figure 34: La simulation de la production sous MineSched

Figure 35: Le développement de l'infrastructure de la simulation sous MineSched

47
 Conclusion :
Ce travail consistait à mener une étude de faisabilité technique et économique de l’extension, et ensuite
de définir une planification long terme qui optimise le séquençage de la production. Pour ce faire, une
étude statistique a été mise en place pour déterminer les caractéristiques morphologiques de l’extension
NW notamment le pendage et la puissance verticale. En utilisant le résultat de cette analyse et par
l’application de l’algorithme AHP, on a déterminé la méthode d’exploitation la plus adapté à l’extension,
cette technique est développée par le langage de programmation R. Les variables d’entré pour
l’algorithme sont les valeurs du pendage et de la puissance et ensuite les coordonnées de chaque bloc
pour construire la cartographie du gisement par méthodes d’exploitations, ces variables sont extraites
depuis le bloc model du gisement construit par les géologues. L’analyse a donné une cartographie où la
majorité des blocs adoptent la méthode des chambres et piliers. On trouve une présence du SNA et du
Cut and Fill vers les bords, cela est dû à un artefact dans le bloc model. Le résultat final opte pour
exploitation par les chambres et les piliers.

Pour optimiser en temps, une application Shiny est mise place pour REMINEX. Elle permet de faire tout le
calcul en interne et afficher le résultat sous forme de graphe coloré par la méthode d’exploitation
convenable à chaque partie du gisement.

L’optimisation des ressources a démarré en calculant la valeur du NSR pour chaque bloc afin de la
comparer avec la teneur de coupure en NSR de la méthode d’exploitation retenue. En se référant à la
zone la plus proche pour en tirer des synthèses géotechniques, car l’étude géotechnique n’est encore faite
pour l’extension NW, on tire les paramètres géométriques des tailles pour lancer le processus
d’optimisation à l’aide du module MSO disponible sur le logiciel Datamine. Le résultat donne des réserves
de 12 639 871 Tonnes du minerai avec une teneur moyenne de 0.92% de cuivre et 17,2 Agppm. Après
avoir supprimer les chambres inexploitable et qui ne justifient pas l’investissement, une infrastructure
était mise en place pour extraire la totalité des réserves, composée de 3 cheminées d’aérages, un puit
pour l’accès à la minéralisation, un système de transport composés par des camions et un skip au niveau
du puit.

Avant de démarrer la production il faut une planification de production qui assurera l’optimisation du
résultat financier durant toute la durée de vie de la mine fixée à 10 ans, cela est fait à l’aide d’une
optimisation linéaire pour le séquençage de l’extraction. Cette optimisation a donné un séquençage
favorisant dans le temps l’extraction du panneaux le plus riche tout en respectant les différentes
contraintes appliqués notamment les contrainte de la continuité, d’accessibilité et d’infrastructure. On
retient alors un scénario d’exploitation générant 286,5 M$ de marge opérationnelle avec un rythme
annuelle de production allant de 1,1 Mt jusqu’à 1,5 Mt.

Pour donner suite à ce travail, il est vivement recommandé d’intégrer un processus de découpage dans
l’algorithme d’optimisation pour étudier plusieurs cas de découpage afin d’en tirer le meilleur. Cela
permettra de trouver des scénarios optimaux avant même de fixer la forme de l’infrastructure de la mine.

48
ANNEXES

49
 ANNEXE CHAPITRE II :
Les matrices de comparaison par paire
Les matrices de comparaison par paire selon les intervalles du pendage :

• Si Le pendage dépasse les 60 degrés :

R&P SR&P PR&P SNA C&F LT

R&P 1 1/5 1/5 1/7 1/7 1/3

SR&P 5 1 1/2 1/5 1/7 3
PR&P 5 2 1 1/5 1/5 3
SNA 7 5 5 1 5 7
C&F 7 7 5 1/5 1 7
LT ( 3 1/3 1/3 1/7 1/7 1 )

Les poids relatifs par rapport à cette matrice sont :

R&P 4.5%
SR&P 10%
PR&P 17%
SNA 47%
C&F 15%
LT (5.5%)

• Si le pendage est supérieur à 20 degrés et inférieur à 60 degrés :

R&P SR&P PR&P SNA C&F LT

R&P 1 1/5 1/5 2 1/7 1/3

SR&P 5 1 1/2 5 1/7 3
PR&P 5 2 1 7 1/5 3
SNA 1/2 1/5 1/7 1 1/7 1/3
C&F 7 7 5 7 1 7
LT ( 3 1/3 1/3 3 1/7 1 )

Les poids relatifs sont :

R&P 4.2%
SR&P 13.7%
PR&P 18.6%
SNA 3.17%
C&F 52.3%
LT ( 7.3% )

50
 • Si le pendage est inférieur à 20 degrés :

R&P SR&P PR&P SNA C&F LT

R&P 1 9 9 9 9 1
SR&P 1/9 1 1 7 6 1/5
PR&P 1/9 1 1 7 6 1/5
SNA 1/9 1/7 1/7 1 1/2 1/7
C&F 1/9 1/6 1/6 2 1 1/6
LT ( 1 5 5 7 6 1 )

Les poids relatifs sont :

R&P 40%
SR&P 10%
PR&P 10%
SNA 2%
C&F 3%
LT (32%)

La matrice de comparaison par paire et par rapport à la puissance est :

• Si la puissance dépasse les 20 m :

R&P SR&P PR&P SNA C&F LT

R&P 1 3 1/2 1/7 5 1

SR&P 1/3 1 1/3 1/7 3 2
PR&P 2 3 1 1/5 5 3
SNA 7 7 5 1 7 5
C&F 1/5 1/3 1/5 1/7 1 1/3
LT ( 1 1/2 1/3 1/5 3 1 )

Les poids relatifs sont :

R&P 11%
SR&P 7%
PR&P 17%
SNA 52%
C&F 3%
LT ( 7% )

51
• Si la puissance est inférieure à 20 m :

R&P SR&P PR&P SNA C&F LT

R&P 1 3 3 7 5 1
SR&P 1/3 1 1/2 7 5 2
PR&P 1/2 2 1 7 5 2
SNA 1/7 1/7 1/7 1 1/5 1/4
C&F 1/5 1/5 1/5 5 1 1/3
LT ( 1 1/2 1/2 4 3 1 )

Les poids relatifs sont :

R&P 25%
SR&P 19%
PR&P 31%
SNA 2%
C&F 6%
LT (14%)

52
 ANNEXE CHAPITRE II :
Le code de l’application Shiny
## library
library(dplyr)

library(ggplot2)

library(expss)

library(shinycssloaders)

library(kableExtra)

library(rstudioapi)

library(shiny)
library(tidyr)

library(plotly)

library(intrval)

## Pendage supérieur à 60 degrés

Pendage_sup60 <- matrix([Link](c(1,1/3,1/3,1/7,1/3,1/2,3,1,1/2,1/5,1/3,3,
3,2,1,1/5,2,3,7,5,5,1,3,7,3,3,1/2,1/3,1,3,2,1/3,1/3,1/7,1/3,1)),6,6,byrow=TRU
E)
ev_Pendage <- eigen(Pendage_sup60)
values <- ev_Pendage$values
vectors <- ev_Pendage$vectors
a = [Link](vectors[,1])
b = [Link](values[1])
c = Pendage_sup60%*%a
c=c/b
c_pendage_sup60 = c/sum(c)
CR = (b - 6)/(5)

## Pendage inférieur à 60 degrés

Pendage_inf60 <- matrix([Link](c(1,1/5,1/5,2,1/7,1/3,5,1,1/2,5,1/7,3,5,2,
1,7,1/5,3,1/2,1/5,1/7,1,1/7,1/3,7,7,5,7,1,7,3,1/3,1/3,3,1/7,1)),6,6,byrow=TRU
E)
ev_Pendage <- eigen(Pendage_inf60)
values <- ev_Pendage$values
vectors <- ev_Pendage$vectors
a = [Link](vectors[,1])
b = [Link](values[1])
c = Pendage_inf60%*%a
c=c/b
c_pendage_inf60 = c/sum(c)
53
CR = (b - 6)/(5)

## Pendage inférieur à 20 degrés

Pendage_inf20 <- matrix([Link](c(1,7,7,9,9,1,1/7,1,1,7,6,1/5,1/7,1,1,7,6,
1/5,1/9,1/7,1/7,1,1/2,1/7,1/9,1/6,1/6,2,1,1/6,1,5,5,7,6,1)),6,6,byrow=TRUE)
ev_Pendage <- eigen(Pendage_inf20)
values <- ev_Pendage$values
vectors <- ev_Pendage$vectors
a = [Link](vectors[,1])
b = [Link](values[1])
c = Pendage_inf20%*%a
c=c/b
c_pendage_inf20 = c/sum(c)
CR = (b - 6)/(5)

## Puissance supérieur à 20m

Puissance_sup <- matrix([Link](c(1,3,1/2,1/7,5,1,1/3,1,1/3,1/7,3,2,2,3,1,
1/5,5,3,7,7,5,1,7,5,1/5,1/3,1/5,1/7,1,1/3,1,1/2,1/3,1/5,3,1)),6,6,byrow=TRUE)
ev_Puissance <- eigen(Puissance_sup)
values <- ev_Puissance$values
vectors <- ev_Puissance$vectors
a = [Link](vectors[,1])
b = [Link](values[1])
c = Puissance_sup%*%a
c=c/b
c_puissance_sup = c/sum(c)
CR = (b - 6)/(5)

## Puissance inférieur à 20m

Puissance_inf <- matrix([Link](c(1,3,1/2,7,5,1,1/3,1,1/2,7,5,2,2,2,1,7,5,
2,1/7,1/7,1/7,1,1/5,1/4,1/5,1/5,1/5,5,1,1/3,1,1/2,1/2,4,3,1)),6,6,byrow=TRUE)
ev_Puissance <- eigen(Puissance_inf)
values <- ev_Puissance$values
vectors <- ev_Puissance$vectors
a = [Link](vectors[,1])
b = [Link](values[1])
c = Puissance_inf%*%a
c=c/b
c_puissance_inf = c/sum(c)
CR = (b - 6)/(5)

## UI function
ui <- fluidPage(
sidebarLayout(
54
 sidebarPanel(
fileInput("df", "Merci d'entrer la base de données",
multiple = FALSE,
accept = c("text/csv")),
actionButton("set","charger les données", icon = icon("Data")),
div(style = "height:20px"),
selectInput("Zone","ZONE", choices = c("")),
selectInput("x","X", choices = c(" ")),
selectInput("y","Y", choices = c(" ")),
selectInput("z","Z", choices = c(" ")),

actionButton("go","Lancer les calculs", icon = icon("calendar"))

),
mainPanel(
div(plotlyOutput("plot", width ="100%", height = "100%"), align =
"center")
)
)
)

# Server logic

server <- function(input, output, session) {

options([Link]=200*1024^2)
v = reactiveValues(plot = NULL)
rv <- reactiveVal()
data <- reactive({
inFile <- input$df
if(){
data_puissance = [Link](inFile$datapath, stringsAsFactors = FAL
SE)
}
data_puissance
data_pendage<<-data_puissance
rv(data_pendage)
data_pendage$ZINC = [Link](data_pendage$ZINC)
data_pendage$TRDIP = [Link](data_pendage$TRDIP)
data_pendage = data_pendage[,]
data_pendage = data_pendage[,]

data_pendage$score_RP = 0
data_pendage$score_SRP = 0
data_pendage$score_PRP = 0
data_pendage$score_SNA = 0
data_pendage$score_CF = 0
data_pendage$score_LT = 0

55
 inf20degré_inf20m <<- subset(data_pendage, ZINC < 20 & TRDIP < 20)

inf20degré_sup20m <<- subset(data_pendage, ZINC >= 20 & TRDIP < 20)

inf60degré_inf20m <<- subset(data_pendage, ZINC < 20 & TRDIP < 60 & T

RDIP >= 20 )

inf60degré_sup20m <<- subset(data_pendage, ZINC >= 20 & TRDIP < 60 &

TRDIP >= 20 )

sup60degré_sup20m <<- subset(data_pendage, ZINC >= 20 & TRDIP >= 60 )

sup60degré_inf20m <<- subset(data_pendage, ZINC < 20 & TRDIP >= 60 )

rv(inf20degré_inf20m)
rv(inf20degré_sup20m)
rv(inf60degré_inf20m)
rv(inf60degré_sup20m)
rv(sup60degré_inf20m)
rv(sup60degré_inf20m)

inf20degré_inf20m$ZINC = (inf20degré_inf20m$ZINC - min(data_pendage$Z

INC))/(max(data_pendage$ZINC) - min(data_pendage$ZINC))
inf20degré_inf20m$TRDIP = (inf20degré_inf20m$TRDIP - min(data_pendage
$TRDIP))/(max(data_pendage$TRDIP) - min(data_pendage$TRDIP))

inf20degré_sup20m$ZINC = (inf20degré_sup20m$ZINC - min(data_pendage$Z

INC))/(max(data_pendage$ZINC) - min(data_pendage$ZINC))
inf20degré_sup20m$TRDIP = (inf20degré_sup20m$TRDIP - min(data_pendage
$TRDIP))/(max(data_pendage$TRDIP) - min(data_pendage$TRDIP))

inf60degré_inf20m$ZINC = (inf60degré_inf20m$ZINC - min(data_pendage$Z

INC))/(max(data_pendage$ZINC) - min(data_pendage$ZINC))
inf60degré_inf20m$TRDIP = (inf60degré_inf20m$TRDIP - min(data_pendage
$TRDIP))/(max(data_pendage$TRDIP) - min(data_pendage$TRDIP))

inf60degré_sup20m$ZINC = (inf60degré_sup20m$ZINC - min(data_pendage$Z

INC))/(max(data_pendage$ZINC) - min(data_pendage$ZINC))
inf60degré_sup20m$TRDIP = (inf60degré_sup20m$TRDIP - min(data_pendage
$TRDIP))/(max(data_pendage$TRDIP) - min(data_pendage$TRDIP))

sup60degré_sup20m$ZINC = (sup60degré_sup20m$ZINC - min(data_pendage$Z

INC))/(max(data_pendage$ZINC) - min(data_pendage$ZINC))
sup60degré_sup20m$TRDIP = (sup60degré_sup20m$TRDIP - min(data_pendage
$TRDIP))/(max(data_pendage$TRDIP) - min(data_pendage$TRDIP))

sup60degré_inf20m$ZINC = (sup60degré_inf20m$ZINC - min(data_pendage$Z

56
INC))/(max(data_pendage$ZINC) - min(data_pendage$ZINC))
sup60degré_inf20m$TRDIP = (sup60degré_inf20m$TRDIP - min(data_pendage
$TRDIP))/(max(data_pendage$TRDIP) - min(data_pendage$TRDIP))

for(i in 1:6){
inf20degré_inf20m[,21+i] = inf20degré_inf20m$ZINC*c_puissance_inf
[i] + inf20degré_inf20m$TRDIP*c_pendage_inf20[i]

inf20degré_sup20m[,21+i] = inf20degré_sup20m$ZINC*c_puissance_sup
[i] + inf20degré_sup20m$TRDIP*c_pendage_inf20[i]

inf60degré_inf20m[,21+i] = inf60degré_inf20m$ZINC*c_puissance_inf
[i] + inf60degré_inf20m$TRDIP*c_pendage_inf60[i]

inf60degré_sup20m[,21+i] = inf60degré_sup20m$ZINC*c_puissance_sup
[i] + inf60degré_sup20m$TRDIP*c_pendage_inf60[i]

sup60degré_inf20m[,21+i] = sup60degré_inf20m$ZINC*c_puissance_inf
[i] + sup60degré_inf20m$TRDIP*c_pendage_sup60[i]

sup60degré_sup20m[,21+i] = sup60degré_sup20m$ZINC*c_puissance_sup
[i] + sup60degré_sup20m$TRDIP*c_pendage_sup60[i]

Result <<- rbind(inf20degré_inf20m,inf20degré_sup20m,inf60degré_inf20

m,inf60degré_sup20m,sup60degré_inf20m,sup60degré_sup20m)
Result <<- rename(Result, c("Room & Pillar"= "score_RP","Step Room &
Pillar" = "score_SRP", "Post Room & Pillar" = "score_PRP","Sous-niveaux abat
tus" = "score_SNA","Cut and Fill" = "score_CF","Longe Taille" = "score_LT"))
Result$result <<- colnames(Result[,c(22,23,24,25,26,27)])[[Link](Res
ult[,c(22,23,24,25,26,27)], [Link] = "first")]
rv(Result)
})

observeEvent(input$set, {
rm(list = ls())
if([Link](input$df)){
showNotification("Veuillez entrer la BDD")
}
else{
showNotification("chargement des données")
if(exists('data_pendage')){
updateSelectizeInput(session, "Zone", choices = unique(data_p
endage$ZONE))
updateSelectizeInput(session, "x", choices = colnames(data_pe
ndage))
57
 updateSelectizeInput(session, "y", choices = colnames(data_pe
ndage))
updateSelectizeInput(session, "z", choices = colnames(data_pe
ndage))

}
else{
data()
updateSelectizeInput(session, "x", choices = colnames(data_pe
ndage))
updateSelectizeInput(session, "y", choices = colnames(data_pe
ndage))
updateSelectizeInput(session, "z", choices = colnames(data_pe
ndage))
updateSelectizeInput(session, "Zone", choices = unique(data_p
endage$ZONE))
}
}

})
observeEvent(input$go, {
if(exists('Result') & !(input$x == " ") &  & !(inp
ut$y == " ") & !(input$z == " ") & ){
showNotification("Lancement des calculs")
v$plot <- plot_ly(Result[Result$ZONE == input$Zone,], x = ~get(input$
x), y = ~get(input$y), z = ~get(input$z), color = ~result,type = 'scatter3d')

v$plot <- v$plot %>% layout(scene = list(aspectmode='manual',

aspectratio = list(x=2, y=
2, z=0.25)),
xaxis = list(title = 'X',
zerolinecolor = '#ffff',
zerolinewidth = 2,
gridcolor = 'ffff'),
yaxis = list(title = 'Y',
zerolinecolor = '#ffff',
zerolinewidth = 2,
gridcolor = 'ffff'),
zaxis = list(title = 'Y',
zerolinecolor = '#ffff',
zerolinewidth = 2,
gridcolor = 'ffff'))

}
if([Link](input$df)){
showNotification("Veuillez entrer la BDD")
}
if([Link](input$df) & ){
58
 showNotification("Veuillez entrer la BDD")
}

if((input$x == " ") |(input$y == " ") | (input$z == " ") ){

showNotification("Veuillez selectionner les variables")
}

})

output$plot <- renderPlotly({

if ([Link](v$plot)) return()
v$plot
})

## Run shiny app

shinyApp(ui = ui,server = server)

59
 ANNEXE CHAPITRE III :
Les différentes modèle d’optimisation
Optimisation en tenant en compte différentes scénarios d’exploitation :
[Link]

60