Ce travail de thèse a consisté à établir une base de données d'éboulements rocheux la plus exhaustive possible (pour des volumes supérieurs à 0,1 m3) pour une paroi rocheuse active dominant l'agglomération grenobloise, en utilisant des données de scanner laser et de photographie. Le site d'étude est une falaise calcaire de plusieurs kilomètres en bordure du massif de la Chartreuse, constituée de deux barres de morphologies différentes, l'une en calcaire lité, l'autre en calcaire massif.Les nuages de points issus du scanner laser permettent de reconstituer la falaise et les compartiments éboulés en 3D. Les données de surface des falaises, ainsi que des informations sur la localisation, les dimensions, le mécanisme de rupture propre à chaque compartiment ont été analysées pour caractériser l'évolution morphologique des deux falaises. Il apparait que la falaise inférieure, dont la morphologie dépend fortement de la fracturation et de l'érosion torrentielle des marnes sous-jacentes, présente une fréquence d'éboulement 22 fois plus importante que la falaise supérieure, de morphologie et pente régulière. De plus, dans la falaise inférieure, le taux d'érosion est 4 fois plus élevé entre 900 et 1000 m d'altitude, qu'entre 1000 et 1100 m. Cela montre que le régime d'érosion de cette falaise est transitoire, alors que celui de la falaise supérieure pourrait être permanent.Les éboulements rocheux détectés ont également été datés par un suivi photographique pendant 2,5 ans. Un suivi quasi-continu (1 photo toutes les 10 min), avec un objectif grand angle a permis de dater 214 éboulements de plus de 0,1 m3. Un suivi mensuel, avec un téléobjectif, a permis de dater 854 éboulements de plus de 0,01 m3. L'analyse de ces deux bases de données montre que la fréquence d'éboulements rocheux peut être 7 fois plus grande lors d'un épisode de gel-dégel que sans évènement météorologique particulier, et 4,5 fois plus grande lors d'un épisode de pluie. De plus, elle devient 26 fois plus grande si l'intensité depuis le début de l'épisode est supérieure à 5 mm/h. A partir de ces résultats, une échelle de 4 niveaux d'aléa a pu être proposée pour la prévision de l'aléa. La base de données plus précise et la définition des épisodes de gel-dégel ont permis de distinguer différentes phases dans un épisode de gel-dégel : refroidissement à température négative, réchauffement à température négative, et dégel (à température positive). Il apparait que les éboulements rocheux se produisent plus fréquemment lors des périodes de réchauffement (à température négative) et de dégel que lors des périodes de refroidissement. Cela suggère que les éboulements sont causés par la dilatation thermique de la glace plutôt que par la dilatation due au changement de phase. Ils peuvent cependant ne se produire que lors du dégel, car la cohésion de l'interface roche-glace peut être suffisante pour tenir le bloc jusqu'à la fonte de la glace. Des expériences in situ et en laboratoire ont permis de mesurer la pression de glace dans une fissure avec écoulement d'eau. Elles montrent que la glace formée par accrétion (gel de gouttes ou films d'eau) n'exerce pas de pression sur les parois rocheuses. / Using laser scanner data, an exhaustive rockfall database (for volume larger than 0.1 m3) has been established for a rockwall located near the town of Grenoble (France). The study site is a long double cliff, on the eastern border of the Chartreuse Massif. The two cliffs consist respectively of thinly bedded and massive limestone, which show different structures, morphologies and rockfall activities.The 3D point clouds obtained by laser scanner allow to detect and model the fallen compartments in 3D. Information about cliff surface, and localization, dimensions, failure mechanism for each compartment were obtained and analyzed in order to characterize the morphological evolutions of the two cliffs. It appears that the morphology and the slope of the lower cliff is related to fracturing and torrential erosion which occurs in the marls below the cliff. The rockfall frequency for this lower cliff is 22 times higher than for the upper cliff. Moreover, in the lower cliff the erosion rate is at least 4 times higher for an elevation between 900 and 1000 m than between 1000 and 1100 m. These results show that the erosion process in the lower cliff is in a transient state, whereas it could be in a steady state in the upper cliff. The morphology and the slope of the upper cliff is more regular than for the lower cliff.Rockfalls have been dated by photographic surveys during 2.5 years. A near-continuous survey (1 photo each 10 mn) with a wide-angle lens have allowed dating 214 rockfalls larger than 0.1 m3, and a monthly survey with a telephoto lens, dating 854 rockfalls larger than 0.01 m3. The analysis of the two data bases shows that the rockfall frequency is 7 times higher during freeze-thaw episodes than without meteorological event, and 4.5 times higher during rainfall episodes. Moreover, it becomes 26 times higher when the mean rainfall intensity (since the beginning of the rainfall episode) is higher than 5 mm/h. Based on these results, a 4-level hazard scale has been proposed for hazard prediction. The more precise data base and freeze-thaw episode definition make it possible to distinguish different phases in freeze-thaw episodes: negative temperature cooling periods, negative temperature warming periods and thawing periods. It appears that rockfalls occur more frequently during warming and thawing periods than during cooling periods. It can be inferred that rockfalls are caused by thermal ice dilatation rather than by dilatation due to the phase transition. But they may occur only when the ice melt, because the cohesion of the ice-rock interface can be sufficient to hold the failed rock compartment until the ice melt. The formation of ice in rock cracks has been studied in the field and in laboratory to highlight its influence on rockfall triggering. It has been shown that ice forming by an accretion process (freezing of water drops) doesn't exert a pressure on the crack walls.

De nombreuses incertitudes liées aux mécanismes d'interaction entre les blocs rocheux et le versant naturel lors des chutes de blocs persistent ; la prévision de tels événements reste de ce fait encore incertaine. Néanmoins, les outils numériques et la puissance de calcul ne cessent d'évoluer. Si, auparavant, les calculs trajectographiques étaient restreints à des géométries simplifiées et à des mouvements balistiques en deux dimensions, il devient désormais possible d'y intégrer des raffinements tels que la forme complexe des blocs, des modèles numériques tridimensionnels de terrain d’une grande définition ou une prise en compte fine des mécanismes dissipatifs au niveau du point d'impact entre le bloc et le versant de propagation.L’objectif principal de la thèse est d’analyser, avec un code numérique discret en trois dimensions, l’influence des paramètres de forme et d’interaction sur la nature du rebond dans un contexte d’ingénierie. Nous présentons tout d’abord une méthodologie d'identification et d'étude de sensibilité des paramètres de contact, élaborée et validée à partir d’expérimentations de laboratoire. Cette méthodologie a été appliquée par la suite à deux expérimentations de chute de blocs menées sur sites réels à moyenne et à grande échelle.L’étude réalisée à moyenne échelle a permis de confronter le modèle numérique à des données obtenues lors d'une campagne expérimentale sur voies ferroviaires commanditée par la SNCF et menée en collaboration avec IRSTEA. Les analyses qui ont été réalisées ont porté sur les vitesses d’impact des blocs avec le ballast et les distances de propagation.L’étude menée à grande échelle s’appuie sur plusieurs séries de lâchés de blocs réalisées sur le site expérimental de la carrière d'Authume dans le cadre d’un Benchmark proposé dans le cadre du Projet National C2ROP. L’objectif principal du Benchmark est de tester et de comparer entre eux des logiciels trajectographiques, des codes de calculs numériques et les pratiques d’ingénierie pour en définir la pertinence et les domaines de validité. Dans le cadre de la thèse, ce travail a été conduit en plusieurs phases (à l’aveugle puis avec des données partielles mesurées lors de la campagne d'essais) et nous présentons l'évolution de ces analyses à l'issue de chacune des phases. L’étude a porté principalement sur les vitesses, les hauteurs et les énergies de passage des blocs en certains points du profil de propagation, ainsi que sur les positions d’arrêt des blocs. Une étude sur l'influence de la forme des blocs sur les distances de propagations est également présentée.Enfin, un Benchmark interne réalisé au sein de l'entreprise IMSRN montre l’importance, sur les analyses, de l'expertise de l'opérateur, et des conséquences de l'utilisation de différents outils trajectographiques (en 2D et en 3D). Ces travaux mettent en lumière les problématiques actuelles auxquelles sont souvent confrontés les bureaux d'études et les ingénieurs en charge des études de risques. / Numerous uncertainties related to the machanical interaction between rock boulders and the natural slope during block falls persist; and the forecast of such events is therefore still uncertain. Nevertheless, digital tools and computing power are constantly evolving. Previously, trajectory calculations were restricted to simplified geometries and two-dimensional ballistic movements, but it is now possible to incorporate refinements such as the complex shape of the blocks, three-dimensional numerical models of terrain of large sizes, as well as a better accounting of the dissipative mechanisms at the point of impact between the block and the run-out slope.The main objective of this work is to analyze, with a discrete elements code in three dimensions, the influence of the shape and interaction parameters on the nature of the rebound in an engineering context. We first present a methodology for identifying and studying the sensitivity of contact parameters, developed and validated from laboratory experiments. This methodology was subsequently applied to two block fall experiments conducted on medium and large real-scale scenarios.The study conducted on a medium scale allowed the numerical model to be compared with data obtained during an experimental rockfall tests campaign commissioned by the SNCF and conducted in collaboration with IRSTEA in a railway. The analyzes that were carried out mainly focused on the impact velocities of the blocks with ballast and propagation distances.The large-scale study is based on a series of block releases performed at the experimental site (Authume quarry, France) as part of a Benchmark proposed inside the National Project C2ROP. The mainly goal of this Benchmark is to access and compare trajectory softwares, numerical computation codes and engineering practices to define their relevance and validity domains. As part of the thesis, this work was conducted in several phases (blind phase and then conducted taking in account partial data measured during the experimental tests) and we present the evolution of these analyzes at the end of each one of these phases. The study focused on the velocities, heights and energies of the blocks at certain points of the propagation profile, as well as on the stopping positions of the blocks. The influence of block shapes on run-out distances is also presented.Lastly, an internal Benchmark performed within the IMSRN company shows the importance, on the analyzes, of the expertise of the operator, and the consequences derived from the application of different trajectography tools (in 2D and in 3D). This work highlights the current issues that are often faced by engineering offices and engineers in charge of risk quantification.

De nombreuses incertitudes interviennent lors de l’évaluation de l’aléa de départ des éboulements rocheux. Ce mémoire se propose de les étudier, en suivant la typologie de Baecher & Christian (2005) qui identifie (1) l’incertitude de décision, (2) l’incertitude de connaissance et (3) l’incertitude aléatoire. L’incertitude de décision intervient au travers de l’évaluation de l’aléa par expertise. Nous avons mis au point une expérimentation visant à étudier l’influence du niveau d’expertise de la personne en charge de l’étude et de la méthode utilisée (qualitative et quantitative) sur l’évaluation de l’aléa. Nous avons montré que l’utilisation d’une méthode qualitative permet de guider la pensée et de réduire la dispersion des niveaux d’aléas finaux, ce qui n’est pas le cas d’une méthode quantitative. Nous montrons également que dans des cas spectaculaires, l’expérience de l’ingénieur lui permet de ne pas surévaluer l’aléa, et que pour des cas plus classiques, il n’y a pas d’influence du niveau d’expertise dès lors que la méthode est suffisamment détaillée. L’incertitude aléatoire se manifeste par la variabilité temporelle des éboulements rocheux. Un travail sur des bases de données a été entrepris. Il a permis de mettre en évidence une corrélation statistique entre les éboulements rocheux et certains facteurs météorologiques (pluies et températures). Une méthode évaluant la probabilité d’occurrence des éboulements en fonction de l’intensité des pluies ou la valeur des températures a alors été développée. Elle peut être considérée comme un outil d’aide à la décision dans la gestion du risque. Un travail sur ces bases de données a également permis de montrer l’existence d’une composante de bruit de fond aléatoire, ne montrant aucune corrélation statistique, caractérisée par le fait que le temps qui sépare deux chutes de blocs suit une loi de Poisson. L’incertitude de connaissance résulte d’un manque d’informations concernant, en particulier, les facteurs préparatoires. Afin de mieux appréhender ces incertitudes, des modèles numériques (éléments distincts - UDEC) de glissement d’un bloc sur un plan incliné, dépendants de la présence et de la proportion de ponts rocheux, ont été réalisés. Ils ont permis de montrer l’existence de deux phases dans la rupture, contrôlées par les reports des contraintes le long du joint entre le bloc et le plan, à mesure que de la rupture des ponts rocheux se propage. Cependant, ces deux phases ne se retrouvent pas dans l’étude des déplacements du bloc. Nous montrons également l’influence de la position des ponts rocheux sur la vitesse de propagation de la rupture en fonction de la pente. Ces travaux fournissent un outil opérationnel et des indications sur le processus d’évaluation de l’aléa qui pourraient aider l’ingénieur en charge d’une étude d’aléa à justifier et à affiner son estimation / The present dissertation proposes to study uncertainties in rockfall hazard assessment process, on the basis of the Baecher & Christian typology (2005), which identifies (1) decision uncertainties, (2) knowledge uncertainties, and (3) natural variability. Decision uncertainties are due to the subjectivity of experts’ assessment. An experiment has been realized in order to evaluate the influence of the expertise level, and the chosen method (qualitative or quantitative) on the rockfall hazard. The rockfall hazard levels obtained by the qualitative method are quite uniform while the quantitative method produces more disparate results. We have shown that, classically, the expertise level has no influence on the assessment if the method is precisely detailed. In the case of spectacular sites, experts do not overestimate the hazard. The natural variability is associated to the temporal variability of the rockfalls. We have studied the statistical correlation between meteorological and rockfall databases. We have developed a method that takes into account the probability of occurrence of the studied triggering factor (rainfalls and temperatures). This new approach is easy to use, and also helps to determine the conditional probability of rockfall according to a given meteorological factor. This approach will help to optimize risk management in studied areas based on their meteorological conditions. A work on these databases also has allowed the unpredictability of the rockfalls to be highlighted. Indeed, the delay between two rockfalls follows a Poisson distribution. Knowledge uncertainties may concern the preparatory factors. We have studied them with numerical models of a rock block sliding along a planar joint, depending on the proportion and position of rock bridges (distinct elements - UDEC). Two stages of the failure process have been highlighted; they have been controlled by the stress redistribution induced by the failure, and the propagation of the rock bridges. However, these two phases are not identified when studying block displacements. We also have shown the influence of the position of rock bridges on the propagation of failure. This research provides an operational tool, and guidelines on assessment process that may help engineers in charge of a rockfall risk analysis to justify and refine the obtained risk estimations

Résumé La prédiction efficace des éboulements rocheux constitue un des moyens pour se prémunirvis-à-vis du risque naturel représenté par les chutes de blocs. Les joints rocheux jouant un rôle décisif dans le déclenchement de ces éboulements, il faut décrire du mieux possible le comportement mécanique de ceux-ci, et en particulier leur rupture. Une nouvelle loi de comportement de joint rocheux, incrémentalement non linéaire, est ainsi proposée, parallèlement à l'utilisation du critère du “travail du second ordre” pour détecter les conditions impliquant un éboulement. La définition de la loi se base sur un modèle numérique de joint rocheux, utilisant la méthode aux Éléments Discrets via le code Yade. Une fois le comportement des joints rocheux investigué grâce à ce modèle numérique (en lien avec des résultats expérimentaux), et la loi définie puis validée, cette-dernière est étudiée vis-à-vis de ce critère du travail du second ordre. L'existence de “directions instables” de sollicitation, susceptibles d'entraîner la rupture du joint rocheux avant le critère de Mohr-Coulomb, est ainsi mise en évidence. Ces directions instables dépendent tout particulièrement des couplages entre les directions normale et tangentielle du joint rocheux. Une falaise existante est enfin analysée : les “Gorges de Valabres”, situées dans les Alpes-Maritimes. L'analyse est effectuée en utilisant la nouvelle loi de comportement, et le critère du travail du second ordre, dans le cadre d'un modèle numérique discret de la falaise. On observe alors que des sollicitations du chargement simulé correspondent pour certains joints à des directions instables.