Titre de l'écran-titre (visionné le 26 février 2024) / Les activités minières étant appelées à se développer à de plus en plus grande profondeur, l'intensité des réponses sismiques aux sautages de segments de galerie de développement devrait augmenter et conséquemment présenter un aléa plus important. Cette sismicité induite représente divers risques sismiques qui pourraient affecter la sécurité du personnel minier, les opérations et la rentabilité de la mine. Les réponses sismiques sont d'intensité variable et leurs facteurs d'influence ne sont pas bien compris. Habituellement, la gestion de l'aléa sismique aux réponses sismiques aux sautages de développement est liée à des règles conservatrices de type règles générales. Ces règles sont parfois ajustées selon l'expérience du personnel de la mine, non quantifiable et difficile à transmettre, ou l'historique sismique. Ces règles n'intègrent pas ou peu les connaissances de terrain. Cette thèse présente une méthodologie pour intégrer quantitativement l'influence de propriétés géologiques, structurales et mécaniques du massif rocheux dans la compréhension et l'anticipation des risques sismiques associés aux sautages de développement de galerie. La méthodologie développée est adaptable à différents contextes et peut être mise à jour par la considération de variables supplémentaires, selon les données disponibles et les conditions du site minier. Cela est possible par la quantification systématique de caractéristiques structurales, géologiques et liées au massif rocheux à chaque segment de galerie de développement. Pour démontrer l'applicabilité de la méthodologie développée, des sautages de développement de la mine LaRonde ont été étudiés. Ces sautages couvrent trois niveaux à 3 km de profondeur et pour une période s'étendant sur 17 mois. La réponse sismique au sautage de 379 segments de galerie de développement a été délimitée et quantifiée. L'acceptabilité des données collectées géologiques, structurales et liées au massif rocheux a été établie pour interpréter 32 variables pour chaque segment de galerie. Puis, des analyses statistiques bivariées et multivariées ont permis de quantifier l'influence des variables géologiques et structurales définies sur les réponses sismiques. Une approche par analyses factorielles de données mixtes (AFDM) a permis une évaluation quantitative de l'impact des variables géologiques et structurales sur les réponses sismiques, ce qui était auparavant difficile à obtenir. Cela permet de mieux évaluer le risque et planifier en conséquence. Cela peut être inestimable pour les opérations minières lors du développement minier dans une zone sismiquement active, qui demande une planification importance à court et à long terme. Les résultats d'AFDM ont montré, pour le secteur à l'étude, que quatorze variables géologiques et structurales, interprétées à chaque segment de galerie influencent significativement l'intensité des réponses sismiques au sautage de développement. Des modèles prédictifs ont été développés afin d'améliorer la compréhension et l'anticipation de la variation de l'intensité des réponses sismiques au sautage de développement. Des modèles d'analyses statistiques multivariées prédictifs d'arbre de régression d'inférence conditionnelle et de forêts aléatoires ont été développés avec 75 % des 379 segments de galeries. Les 100 segments restants ont été utilisés pour valider la performance des modèles. La stratégie d'utilisation des modèles prédictifs de forêts aléatoires proposée fournit des données quantitatives sur la variabilité des aléas sismiques associés au sautage de développement sur lesquelles les gestionnaires peuvent se baser. Les modèles prédictifs développés peuvent être utilisés pour comprendre, gérer et communiquer l'aléa sismique lié au sautage des segments de galerie de développement. Ils fournissent une évaluation quantitative de l'aléa sismique permettant aux décideurs de sélectionner des seuils de critères de performance jugés acceptables, selon leur tolérance au risque. La combinaison de l'approche proposée avec les protocoles sismiques actuels utilisés sur différents sites miniers permet d'améliorer la gestion du risque sismique associé au sautage de développement. En effet, l'utilisation du modèle prédictif pour le secteur et la période étudiés a permis d'augmenter l'exactitude, la sensibilité et la précision pour anticiper une réponse sismique d'intensité élevée à un sautage de développement. / As mining activities are expected to develop at greater depth, the intensity of seismic responses to blasting of development drift should increase and consequently present a greater hazard. This induced seismicity represents various seismic hazards that could affect the safety of mining personnel, operations, and mines' profitability. Seismic responses intensity varies and the influencing factors causing that variability are not well understood. Usually, the management of seismic hazard to seismic responses to development blasting is linked to conservative blanket rules. These rules are sometimes adjusted according to the experience of mine personnel, which is unquantifiable and difficult to communicate, or the seismic history. These rules do not consider field knowledge or at a minimum. This thesis presents a methodology to quantitatively integrate the influence of the rock mass's geological, structural, and mechanical properties in understanding and anticipating seismic risks associated with development blasting. The developed methodology is adaptable to different contexts and can be updated to consider additional variables, depending on available data and mining site conditions. This is possible by the systematic quantification of structural, geological, and rock mass-related characteristics at each segment of development drift. To demonstrate the applicability of the developed methodology, development blasts from the LaRonde mine were studied. These blasts cover three levels at a depth of 3 km and for a 17 months-period. The seismic response to the blasting of 379 segments of development drift was delineated and quantified. The collected geological, structural, and rock mass data acceptability was established to interpret 32 variables for each drift segment. Then, bivariate, and multivariate statistical analyzes made it possible to quantify the influence of the defined geological and structural variables on the seismic responses. A factor analysis of mixed data (FAMD) approach allowed a quantitative assessment of the impact of geological and structural variables on seismic responses, which was previously difficult to obtain. This allows to better assess risk and plan accordingly. This is invaluable to mining operations when mining in a seismically active area, which requires significant short- and long-term planning. For the studied sector, the FAMD results showed that 14 geological and structural variables, interpreted at each drift segment, significantly influence the intensity of seismic responses to development blasting. Predictive models have been developed to improve the understanding and anticipation of variations in the intensity of seismic responses to development blasting. Predictive multivariate statistical analysis models of conditional inference regression trees and random forests were developed with 75% of the 379 drift segments. The remaining 100 segments were used to validate the performance of the models. The proposed strategy for using predictive random forest models provides quantitative data on the variability of seismic hazards associated with development blasting that managers can rely on. The developed predictive models can be used to understand, manage, and communicate the seismic hazards associated with development blasting. They provide a quantitative seismic hazard assessment, allowing decision-makers to select acceptable performance criteria thresholds, according to their risk tolerance. Combining the proposed approach with current seismic protocols used on different mining sites makes it possible to improve seismic risk management associated with development blasting. Indeed, the use of the predictive model for the sector and period studied made it possible to increase the accuracy, sensitivity, and precision to anticipate a high-intensity seismic response to a development blast.

L’ensemble des mines – et particulièrement celles exploitant un gisement filonien ou de faible puissance – sont très sensibles à la problématique de la dilution minière (i.e ., bris hors-profil de l’éponte supérieure). En effet, puisque cette dilution est généralement de faible teneur, le matériel qui se retrouve dans les chantiers miniers peut ne pas permettre de compenser les coûts associés au minage et par conséquent, la rentabilité de l’opération peut en être affectée. Ceci est particulièrement vrai lorsque la teneur de la dilution se situe sous le seuil de rentabilité (i.e., teneur de coupure). Par ailleurs, lorsque des effondrements (dilution) surviennent pendant les activités de soutirage, ceci représente un risque important pour les équipements mobiles (e.g., chargeuse-navette). Nombreux sont les chercheurs qui ont travaillé à développer et à améliorer les abaques de stabilité et critères de conception de chantiers, mais leur efficacité à prédire le bris hors-profil est souvent faible. Néanmoins, ces outils demeurent très utiles pour les projets sur sites vierges (« greenfield projects »), mais dès lors que les conditions géomécaniques d’un site deviennent mieux comprises, il importe de corriger et d’optimiser ces critères de conception afin de mieux représenter la réalité de l’opération minière. Ce mémoire expose donc la méthodologie employée afin de quantifier le bris hors-profil provenant des différentes surfaces des chantiers (i.e., épontes, toit, murs) exploités par la méthode de chambre vide. Une comparaison entre les performances réelles et celles prédites par les modèles couramment utilisés dans l’industrie a été effectuée. À l’aide de méthodes statistiques univariées et bivariées, les paramètres affectant le bris hors-profil peuvent être déterminés objectivement. L’élaboration de nouveaux modèles statistiques permet également d’accroître la fiabilité des prédictions (i.e., le bris hors-profil) et fournit des outils permettant d’assister les ingénieurs dans la conception et le dimensionnement des chantiers miniers. La mine Éléonore de Goldcorp a servi de cas d’étude afin de démontrer la validité d’une telle démarche.

La déformation due au fluage d'un massif rocheux autour d'un tunnel a été rencontrée fréquemment. Ce phénomène est évident où il y a des tunnels creusés dans la roche tendre, des masses rocheuses faible et fortement cisaillées, ou des massifs rocheux soumis à des contraintes in-situ élevées. La déformation due au fluage se produit fréquemment au moment d’excavation des tunnels longs où il y a des failles et des zones fracturées et cisaillées. Ce phénomène peut causer différents dommages sur des systèmes de soutènement en raison de la déformation excessive et des effondrements. La déformation excessive impose une ré-excavation de la section du tunnel, qui monte le coût supplémentaire, la durée de la réalisation du projet et le risque de la sécurité sur le projet. En plus, comme la stabilité de terrain est dans un état critique durant la ré-excavation, une petite négligence peut conduire à une grande caverne. Bien que la déformation de fluage est commune dans un massif rocheux à une faible résistance dans un tunnel très profond, mais ce phénomène a été observé dans des tunnels peu profonds. Une bonne compréhension des déformations causées par une excavation souterraine requiert la connaissance de l'interaction roche-support et l'interprétation des données de terrain. Auparavant, l’objet principal de la surveillance effectuée durant la construction du tunnel était des mesures de la pression au terrain imposé sur le revêtement du tunnel. Mais aujourd’hui, les méthodes modernes de construction de tunnel se concentrent sur la surveillance des déplacements pendant et après la construction. Afin de déterminer des déformations dans les tunnels, Panet et Sulem ont supposé que "Le tunnel a une section transversale circulaire et le milieu est homogène et isotrope, aussi le tunnel est suffisamment profond pour considérer que la distribution des contraintes est homogène". Mais dans le cas quasi réel, la distribution de la contrainte autour du tunnel est hétérogène et anisotrope. Dans cette étude, pour la modification des équations Panet et Sulem, certaines équations sont proposées en cas de matériau hétérogène et anisotrope pour généraliser le problème. La galerie de force motrice Seymareh a été considérée comme l’étude de cas. Celle-ci est une partie du conduit d’eau dans le projet de centrale électrique du barrage Seymareh. Ce projet est situé à l'ouest de l'Iran. Les données de surveillance de la galerie de force motrice sont collectées au moment de l’excavation du tunnel, et sont comparées avec les résultats de la modélisation numérique et de la solution analytique. Cette comparaison montre que les résultats des données expérimentales obtenues par la surveillance sont très proches des résultats de la solution analytique, mais il y a une différence entre les deux et la modélisation numérique. Il était prévisible, car l’effet d’autres activités comme l’excavation des tunnels verticaux n’est pas prise en compte dans l’analyse numérique et aussi dans la solution analytique. Il est évident que les autres activités comme l’excavation des tunnels verticaux et l'excavation du tunnel principal vers deux directions opposées, peuvent affecter sur les résultats de la surveillance. D'autre part, les données initiales utilisées dans l'analyse numérique et la solution analytique ne sont pas tout à fait exactes, car elles sont obtenues en tant que représentatives du massif rocheux de la région, mais pas pour une section particulière. Toutefois, le but de cette étude est le développement d'une solution analytique de la déformation dans les tunnels sur les conditions générales et la poursuite de cette étude pourra être plus développée. / The creep deformation of a rock mass around a tunnel has been encountered frequently. It is particularly common in tunnels excavated in soft rock, heavily sheared weak rock masses or rock masses subjected to high in-situ stresses. Creep deformation in fault and shear fractured zones are one of the frequently encountered difficulties in long tunnel construction, which tend to cause failure of supporting systems due to excessive deformation and cavern. Excessive deformation would necessitate re-mining of the tunnel cross section, thus imposing impacts such as extra cost, extended time schedule and safety risk on the project. Furthermore, as the ground stability is in critical condition during re-mining, the slightest negligence would lead to major cavern. Although creep deformation is common to extremely poor rock mass under high overburden in a tunnel alignment, but however this phenomenon is not limited to tunnels with high overburden. A good understanding of the deformations caused by an underground excavation requires simultaneously knowledge of the rock-support interaction and interpretation of field data. Formerly, the main purpose of the monitoring carried out during tunnel construction was to measure the ground pressures acting on the tunnel lining. Modern tunneling practice emphasizes the monitoring of the displacements occurring during and after the construction. Panet and Sulem for determining of deformations in tunnels have assumed that "The tunnel has a circular cross section and around the tunnel, the rock is homogeneous and isotropic and also the tunnel is deep enough to consider that the stress distribution is homogenous". But in almost real cases, the stresses distribution around the tunnel is not homogeneous and isotropic. In this study, for modification of the Panet and Sulem equations, some equations are proposed in case of nonhomogeneous and anisotropic for generalizing of the problem. Seymareh power tunnel which is considered as a case study is a part of the powerhouse waterways system of the Seymareh dam and hydroelectric power plant project. The project is located in west of Iran. The monitoring data of power tunnel which are collected during excavation of tunnel is compared with the results of numerical modelling and analytical solution results as well as. The results obtained from comparison show although the field data, which are collected through the monitoring, are very close to the analytical solution results (approximately), but there is a significant difference between both of them and numerical modelling results. It was predictable; because the influence of the other activities such as excavation of shaft and surge tank in the numerical analysis and also analytical solution are not considered. It is obvious that other activities such as excavation of shaft and surge tank and also excavation of mean tunnel from other direction which were under operation at the same time can effect on the results of monitoring. On the other hand, the initial data which are used in numerical analysis and analytical solution are not quite accurate; because they are extracted as a representative of the rock mass of region, not for a particular section. However the goal of this study is development of analytical solution of deformation in tunnels on general conditions and pursuit of the study could be leaded to more development in this field.