La performance de l'effet de barrière hydraulique d'un noyau de barrage en remblai est hautement dépendante des conditions au compactage observées au cours de la construction, et conjointement aux propriétés géotechniques des matériaux, peuvent varier significativement dans l'espace. Des études récentes (Qin, 2015; Venkovic et al., 2013; Smith et Konrad, 2011) ont montré que les géostatistiques peuvent être utilisées pour modéliser la variabilité spatiale des propriétés géotechniques des barrages, tel que la conductivité hydraulique, à l'aide de données géoréférencées issus des essais de contrôle réalisés au cours de la construction. Ces études ont permis une meilleure compréhension de divers problèmes liées à nature stratifiée de la structure des barrages et ont élargi le champ d'application des géostatistiques dans le cadre de l'analyse des barrages en remblai. Ce mémoire porte sur l'amélioration de l'approche géostatistique utilisée pour modéliser la conductivité hydraulique du noyau en till compacté d'un barrage en remblai. Le site étudié est le barrage Sainte-Marguerite-3 (SM-3), érigé sur la rivière Sainte-Marguerite au nord-ouest de la ville de Sept-Îles au Québec. Une démarche de calcul faisant appel au taux d'agrégation des matériaux a d'abord été développée afinn d'obtenir une évaluation plus représentativede la conductivité hydraulique. Une structure d'interpolation spatiale quasi-3D a été utilisée dans le cadre de l'analyse géostatistique afin d'estimer la conductivité hydraulique du noyau en tout point de la structure. La priorité d'analyse des paramètres géotechniques utilisés dans le calcul de la conductivité hydraulique a ensuite été étudiée. Pour cet aspect du projet, les géostatistiques ont été utilisées pour analyser chaque paramètre individuellement, puis les estimations obtenues ont servi à calculer la conductivité hydraulique. Finalement, une analyse géostatistique par domaine a été réalisée dans le but de modéliser plus fidèlement l'importante variabilité spatiale des propriétés géotechniques, observées dans certaines zones du noyau, engendrée par la variation de l'origine des matériaux. Les résultats montrent que l'utilisation d'une structure d'interpolation spatiale quasi-3D est une alternative efficace à une analyse en 1D utilisée dans le cadre de travaux précédents et que ce type de structure d'interpolation respecte la continuité spatiale des matériaux. Les résultats montrent aussi qu'une analyse multi-paramétrique par domaine, motivée par la réalité du terrain, offre une meilleure estimation, particulièrement en ce qui concerne l'évaluation de l'influence des conditions au compactage. / The hydraulic barrier performance of the compacted till core of embankment dams is highly dependent on compaction conditions during construction, which alongside with basic geotechnical properties, may significantly vary spatially. Recent studies showed that geostatistics can be used to represent the spatial variability of the geotechnical properties of dam materials, such as hydraulic conductivity, using georeferenced construction control data. These studies led to advances in the understanding of various issues related to the stratified nature of dams' structure and broadened perspectives for geostatistical analyses of large embankment dams. This research focuses on the improvement of the geostatistical approach used to model the hydraulic conductivity of an embankment dam's compacted till core. The studied site is the Sainte-Marguerite-3 (SM-3) dam, built on the Sainte-Marguerite river, northwest to the city of Sept-Îles in Québec. A mathematical approach using an aggregation rate was first developed as a way to obtain a more representative assessment of the dam's hydraulic conductivity. A quasi-3D spatial interpolation structure was used for the geostatistical analysis to estimate the hydraulic conductivity for every location in the core. A priority analysis of the multiple geotechnical parameters used to compute hydraulic conductivity was then performed. For this aspect of the project, geostatistics were used to individually analyze the spatial variability of the parameters measured during construction, which were then used to compute the hydraulic conductivity. Finally, in order to provide a solution to a specific problem observed in the SM-3 dam, a geostatistical analysis using multiple domains was performed. This approach was used to accurately model the important spatial variability of the geotechnical parameters, observed in some areas of the dam, caused by variations in the origins of the materials. Results showed that the quasi-3D spatial interpolation structure used for this project is an effective alternative to a 1D analysis used in past projects and that this type of interpolation structure respects the materials' spatial continuity. The results also showed that a multi-parametric approach using domains, motivated by the terrain's reality, offers a better estimation, especially when it comes to assessing the influence of compaction conditions.

Titre de l'écran-titre (visionné le 4 décembre 2023) / À l'échelle mondiale, les zones de cisaillement sont les plus importants hôtes de minéralisation aurifère de type « or orogénique ». Ces failles fragiles-ductiles sont la conséquence de processus physiques liés aux cycles orogéniques et se forment dans des conditions métamorphiques postérieurs au pic. Les périodes de raccourcissement crustale génèrent des zones de cisaillement inverses à l'échelle régionale qui servent de conduits principaux pour les fluides hydrothermaux ascendants, formant ainsi des veines mésothermales de quartz-carbonate. L'origine des fluides hydrothermaux, les processus physico-chimiques et la chronologie qui mènent à la précipitation de l'or dans ces veines sont l'objet de recherches de longue date, en particulier à l'échelle des gisements et des corps minéralisés. Néanmoins, à l'échelle des districts aurifères orogéniques, comprenant généralement des champs filoniens régionaux de plusieurs kilomètres de large, la nature de ces paramètres demeure énigmatique. Dans cette thèse, on explore, grâce à une combinaison de techniques micro-analytiques chimiques et isotopiques, les contrôles spatiaux et temporels de plusieurs gisements aurifères appartenant au district de classe mondiale de Malartic-Val-d'Or, au sud de la sous-province néoarchéenne de l'Abitibi (Canada). Dans ce district, les gîtes d'or sont concentrés à proximité de la zone de faille Larder Lake-Cadillac orientée est-ouest et sont encaissés dans des structures d'ordres inférieurs recoupant les roches volcaniques et les intrusions de la ceinture de roches vertes. Dans ces structures, l'or est généralement associé à deux ensembles de veines qui se forment avant ou après la fabrique pénétrative d'étendue régionale (S2) développée au cours d'un raccourcissement N-S majeur au pic du métamorphisme. Il s'agit notamment (1) de veines de quartz-carbonate fortement déformées intégrées dans la fabrique S2 et (2) de veines fragiles-ductiles laminées de quartz-tourmaline-carbonate à cisaillement inverse ou à extension sub-horizontale qui sont subparallèles ou obliques à la fabrique S2. Ces dernières font partie du Champ filonien de Val-d'Or (CVO), qui contient la majeure partie de l'or orogénique du district de Malartic-Val-d'Or. Dans cette thèse, j'ai étudié neuf corps minéralisés aurifères (Kiena-S50, Canadian Malartic East, Kiena-Deep, Goldex, Triangle, Plug #4, Beaufor, Pascalis Gold Trend et Akasaba West) afin d'aborder trois sujets clés. 1) La chronologie à l'échelle régionale des événements hydrothermaux aurifères dans la région a été résolue par géochronologie U-Pb in situ. 2) L'identification des processus physico-chimiques contrôlant la saturation en or dans les fluides hydrothermaux aurifères a été effectuée à l'aide d'analyses in situ des éléments traces et de la compositions isotopique du soufre dans les sulfures. 3) La remobilisation hydrothermale de l'or contenu dans la pyrite a été étudiée à l'échelle nanométrique par microscopie électronique à transmission. [...] Les données montrent que la minéralisation aurifère orogénique à l'échelle d'un district est le résultat de processus multiples qui opèrent à différentes échelles, dans l'espace et dans le temps, tout au long d'un cycle orogénique. Des événements hydrothermaux de courtes durées, discrets au niveau régional, au cours du trajet métamorphique rétrograde d'un cycle orogénique, constituent une condition préalable essentielle à la formation d'un district aurifère. Au cours d'un tel épisode de minéralisation aurifère, la majeure partie de la CVO est formée par un mécanisme principal, à savoir la sulfuration de la roche encaissante dans les roches méta-volcaniques et intrusives, conduisant à la précipitation d'inclusions d'Au-Ag-Te-Bi dans la pyrite. Les événements hydrothermaux ultérieurs peuvent façonner et modifier la minéralisation aurifère orogénique à l'échelle du district, comme le montrent les assemblages Au-Ag-Te-Bi remobilisés dans les fractures de pyrites. Malgré un cadre structural bien compris et des processus physiques communs conduisant à la formation de minéralisations aurifères orogéniques, les techniques de micro-analyse in situ offrent l'opportunité unique de relier les observations structurales aux données géochimiques détaillées à l'échelle du grain minéral. Cette approche améliorera considérablement la chronologie de la minéralisation aurifère et la reconnaissance des processus physico-chimiques clés agissant à l'échelle du disctrict, et présentera différentes perspectives sur les processus physico-chimiques associés à la précipitation de l'or à faible et moyenne teneur. / Brittle-ductile reverse shear zones are the most important host for orogenic gold mineralization globally. These reverse shear zones are known to form by common physical processes associated with orogenic cycles, at post-peak metamorphic conditions. Crustal shortening periods induce regional-scale reverse shear zones that act as main conduits for ascending hydrothermal fluids, forming mesothermal quartz-carbonate veins. Fluid source reservoir(s), physico-chemical processe(s) and the timing that lead to gold precipitation within these veins have been the subject of long-standing research activities, particularly focusing at the deposit and ore body-scale. Nevertheless, the nature of these parameters remains elusive at the orogenic gold district-scale, commonly encompassing regional, km wide vein fields. In this thesis, I explore through a combination of micro-analytical chemical and isotopic techniques, spatial and temporal controls on several gold ore bodies contained in the world-class Malartic-Val-d'Or district of the southern Neoarchean Abitibi subprovince (Canada). Structurally-controlled gold deposit clusters in this district are distributed along the E-W striking Larder Lake-Cadillac-fault zone and hosted in lower-order structures that cut volcanic rocks and intrusions of the greenstone belt. In these structures, gold is commonly associated with two vein sets that form pre- or syn- the regional-penetrative fabric (S2) developed during major N-S shortening and peak-metamorphism. This include: (1) highly-deformed quartz-carbonate veins entrained in the S2-fabric and (2) laminated brittle-ductile quartz-tourmaline-carbonate reverse shear and sub-horizontal extension-veins that are sub-parallel to oblique to the S2-fabric, the latter form part of the Val-d'Or Vein Field (VVF), which contains the bulk of orogenic gold in the Malartic-Val-d'Or district. In this thesis, I studied nine gold ore bodies (Kiena-S50, Canadian Malartic East, Kiena-Deep, Goldex, Triangle, Plug #4, Beaufor, Pascalis Gold Trend and Akasaba West) to tackle three key topics. 1) The regional scale timing of auriferous hydrothermal events in the area was resolved by in-situ U-Pb geochronology. 2) Physico-chemical processes controlling gold saturation in auriferous hydrothermal fluids by in-situ sulfide multiple sulfur isotope and trace element compositions. 3) Hydrothermal remobilization of gold hosted in pyrite were studied at the nano-scale by transmission electron microscopy. [...] Acquired data show that district-scale orogenic gold mineralization is the result of multiple processes that operate at various scales in space and time throughout an orogenic cycle. Regionally discrete, short-lived hydrothermal fluid events during the retrograde path of an orogenic cycle are a main pre-condition to form district scale gold mineralization. During such an auriferous fluid event, bulk of the VVF is formed through a main fluid mechanism, wallrock sulfidation in meta-volcanic and intrusive rocks, which controls precipitation of Au-Ag-Te-Bi inclusions in pyrite. Subsequent hydrothermal events have the potential to shape and modify district-scaleorogenic gold mineralization, as recognized in remobilized Au-Ag-Te-Bi assemblages within pyrite fractures. Despite a well understood structural framework and common physical processes that lead to the formation of orogenic gold mineralization, in-situ micro-analytical techniques offer the unique opportunity to link structural observations with detailed geochemical data at the mineral grain scale. This approach will significantly improve the timing of gold mineralization and recognition of key physico-chemical processes acting at the disctrict-scale, as well as present different perspectives on physico-chemical processes associated with low- and medium-grade gold precipitation.

Les travaux présentés sont relatifs à l'analyse des mouvements de versant et au dimensionnement de talus artificiels de mines et carrières à ciel ouvert de grande profondeur. Une démarche méthodologique pluridisciplinaire est développée, dans laquelle l'apport des sciences de la Terre est souligné, de même que l'apport des sciences mécaniques et des sciences hydrologiques. La reconnaissance détaillée des structures et des matériaux géologiques ainsi que l'analyse et la modélisation des processus physiques contrôlant la déformation et la rupture des versants ou talus instables constituent les axes principaux de cette recherche. Une juste appréciation des apports et limites de la modélisation reste cependant une préoccupation majeure, considérant la variabilité des objets géologiques et la complexité des processus naturels étudiés. La recherche s'est orientée, dans une première étape, vers le développement d'outils logiciels permettant une représentation géométrique des milieux rocheux discontinus (logiciel SIMBLOC) ou la modélisation géomécanique des milieux rocheux discontinus (logiciel DEGRES, logiciel BRIG3D). L'analyse géologique et géomécanique de terrain, confrontée à la simulation numérique des processus de déformation et de rupture des grands talus de mines à ciel ouvert ou du génie civil, prenant appui sur ces logiciels et d'autres outils (logiciel Bishop-Carter du CGI, logiciel Flac d'Itasca) a permis de mettre en évidence la multiplicité de configurations géomécaniques toujours prédéterminées par la structure géologique du massif considéré. Ainsi nous avons pu montrer l'importance d'une connaissance détaillée du système de fracturation affectant un massif rocheux pour analyser les conditions de déformabilité et de rupture. L'attention doit porter sur les familles directionnelles de discontinuités, leur hiérarchisation et leur degré de connectivité, leurs paramètres géométriques et géotechniques. Des applications nombreuses ont été expérimentées en situation de grands talus de mines à ciel ouvert (mines en France, au Maroc, en Iran, au Brésil, en Russie, ...) ou de talus du génie civil (écluse à bateaux du barrage des Trois-Gorges en Chine, ...). La réalisation de simulations numériques a permis de mettre en évidence différents aspects de la déformation et de la rupture progressive d'un versant ainsi que le rôle des pressions interstitielles dans ces processus. Des effets de seuil, résultant de l'évolution lente d'un facteur déclenchant, tel qu'un niveau piézométrique, mais se traduisant dans le massif par une évolution brutale des conditions d'équilibre mécanique ont pu être mis en évidence. L'identification de telles situations reste primordiale dans une perspective d'alarme et de gestion des risques. Des processus fondamentaux associés aux conditions de la déformation ont pu être analysés tels que la dilatance d'une masse rocheuse en cours de déformation, contrôlée par les directions d'anisotropie du massif, celles-ci résultant de la structure géologique. Le suivi numérique des processus de localisation des déformations suivant une surface ou une zone de cisaillement et l'analyse du rôle de la structure géologique (favorisant ou contrariant ce processus) font actuellement l'objet de nouvelles investigations. L'objectif est de mettre en évidence des conditions d'initiation, de progression, d'arrêt ou de bifurcation dans le développement de ces processus. L'utilisation des outils numériques permet également de mieux comprendre certains scénarios types tels que celui du déclenchement ou de la réactivation de mouvements de versant en réponse à la montée d'un plan d'eau d'une retenue de barrage ou en réponse à des fluctuations de ce plan d'eau. Nous avons analysé ce sujet sur le plan géomécanique et hydrogéologique, nous référant en particulier à la catastrophe du Vaïont (Italie). Puis, dans une application relative à la retenue du futur barrage des Trois-Gorges (Chine), nous avons confirmé et illustré ces résultats avec l'utilisation du logiciel Flac, et apporté des conclusions opérationnelles. Dans le domaine des risques naturels, de nombreux travaux de terrain, complétés par des caractérisations de matériaux au laboratoire et des simulations numériques, ont permis de préciser les conditions de déclenchement de glissements de terrain et de génération des laves torrentielles. Des conclusions importantes ont été apportées sur les conditions de déclenchement de laves torrentielles, les processus d'initiation et de mobilisation des matériaux de zones sources, la caractérisation géotechnique de ces matériaux, l'analyse de scénarios d'événement utilisant des techniques d'arbres d'événement et d'arbres de défaut. Ces travaux se poursuivent avec l'objectif de mettre en évidence l'influence des actions anthropiques sur l'évolution de la torrentialité. De nouveaux travaux sur l'analyse des grands mouvements de versant : écroulements rocheux, avalanches de roches, sturzstroms sont en cours de développement. Ces phénomènes de grande ou très grande ampleur (glissement rocheux - éboulement ou sturzstrom de Flims, en Suisse : 12.109 m3) font l'objet d'investigations détaillées de terrain et d'analyses géomécaniques. Les conditions du déclenchement sont analysées ainsi que les processus de propagation des masses rocheuses. Dans certains des cas étudiés, ces processus de propagation résultent de l'apparition de quelques surfaces ou zones de cisaillement majeures au sein desquelles une très forte réduction granulométrique et un remaniement intense se réalisent, alors que la structure originelle de la masse rocheuse est préservée au sein des très grands volumes déplacés sur des distances pluri-kilométriques. Le thème délicat de la prévision de rupture pour différentes catégories de mouvements de versant fera bientôt l'objet de nouveaux travaux, prenant en compte les divers contextes géologiques, géomorphologiques et géomécaniques, s'appuyant sur une meilleure compréhension des processus de déformation et de rupture, dans une perspective de prévention des risques. D'une façon plus générale, considérant que tout projet de génie civil ou d'aménagement de l'espace se doit d'intégrer la part d'incertitude inhérente à toute construction intellectuelle basée sur des éléments naturels qui ne sont connus qu'avec une précision limitée ou qui présentent une variabilité impossible à quantifier avec précision, j'envisage d'approfondir le sujet suivant : comment apprendre à gérer les incertitudes, plus spécialement dans un contexte d'analyse de mouvements de versants, éventuellement en situation de crise. Un certain nombre d'outils de diverses natures sont présentés dans ce mémoire, mais c'est toute une démarche d'apprentissage face à une situation de risque qui mérite d'être développée, intégrant l'identification de scénarios d'évolution, la caractérisation des paramètres d'intensité de ces scénarios jugés réalistes et précisant pour un scénario d'intensité donnée un délai d'occurrence. Cette problématique fera l'objet de travaux dans les années à venir.

Ce travail de recherche est axé sur l'étude des instabilités gravitaires dans les Pyrénées et sur le rôle des facteurs géologiques, hydrogéologiques et mécaniques dans le déclenchement des glissements de terrain. Cette réflexion est menée à travers l'étude de deux types différents de glissements de terrain : l'un en milieu rocheux, l'autre en milieu granulaire. Dans le premier cas il s'agit du glissement rocheux historique de Pleysse datant de 1982 et situé en Vallée d'Ossau, dans l'autre cas il s'agit du glissement de Gouetsoule activé en 1997 et localisé en Vallée d'Aspe. Au terme de ce travail, les différents facteurs géologiques, hydrogéologiques et mécaniques ont été définis pour chaque glissement de terrain afin de déterminer les facteurs déclenchants. Leurs études ont aussi permis d'améliorer la compréhension d'une instabilité gravitaire tant en milieu granulaire que rocheux.

Les études de terrain indiquent qu'une faille n'est pas une simple surface de discontinuité où se concentre le glissement mais correspond à une bande endommagée où se distribue la déformation, la zone de faille (ZF). Caractériser la structure de la ZF permet de mieux comprendre les processus de croissance d'une faille i.e. les modes de localisation et de propagation de la rupture. Le couplage de l'analyse structurale fine avec l'analyse statistique des marqueurs de la déformation permet de décrire à plusieurs échelles la structure interne de ZF affectant des grès. La typologie, les modes d'association, et la distribution des marqueurs dans une ZF indiquent que les paramètres intrinsèques du grès (porosité, ciment) et les conditions physico-chimiques environnementales (P,T) contrôlent les mécanismes de déformation. Ces mécanismes déterminent la géométrie et les modes de développement de la ZF, et donc le comportement de la faille par rapport à la rupture et aux circulations de fluides.

Ce travail est une contribution au lever géologique à l'échelle du 1/100000 du CAMEROUN.

Cette étude concerne les régions des Monts de Judée, d'Artuf, vallée d'Ayalon; Gezer; plaine maritime. La stratigraphie, et tectonique sont abpordées.

L'aménagement du Nahr Beyrouth ou fleuve de Beyrouth entre dans le cadre de la planification des eaux libanaises. Le barrage de QARAOUN , pour l'aménagement du Nahr el Litani qui est en cours de réalisation dans la Bekaa Sud, représente le projet principal de cette campagne. Etant donné le profit appréciable que l' économie nationale peut en tirer, tant sur l e plan industriel pour la production de l'énergie que sur le plan agricole (irrigation et mise en valeur des terrains en friche), le gouvernement libanais a décidé de projeter une série de barrages pour l'utilisation systématique de certains cours d ' eau permanents ou périodiques de nature torrentielle. L'aménagement du fleuve de Beyrouth vient s'intégrer dans ce projet

Une brèche est un ensemble de blocs anguleux noyés dans un ciment de nature variable. Les brèches hydrothermales sont générées par un processus de fracturation, de dissolution des fragments, ainsi que des changements de composition causés par des eaux souterraines sous pression à haute température. La nature de la majorité des processus impliqués dans la formation des brèches hydrothermales est bien comprise d'un point de vue géochimique et plusieurs modèles basés sur cette perspective existent. Par contre, il n'existe pas de modèles approchant ces processus d'un point de vue géométrique. Dans ce mémoire, nous proposons un modèle basé sur les automates cellulaires, capable de simuler les principaux processus qui interviennent dans la formation des brèches. Un automate cellulaire est un modèle discret qui consiste en une grille de cellules pouvant chacune prendre à un instant donné un nombre fini d'états. Le temps est également discret et l'état d'une cellule au temps t est fonction de l'état au temps t -1 d'un nombre fini de cellules appelé son voisinage. À chaque nouvelle unité de temps, les mêmes règles sont appliquées pour toutes les cellules de la grille, produisant une nouvelle génération de cellules dépendant entièrement de la génération précédente. Cette approche est compatible avec l'aspect discret de la dissolution des minéraux et permet l'étude de l'évolution géométrique de fragments de roche virtuelle. Plus spécifiquement, on veut mesurer la complexité morphologique des fragments par leur dimension fractale de bordure, une méthode de mesure utilisée sur des échantillons réels et permettant de valider notre modèle avec des données analogiques. Un simulateur a été conçu pour mettre en oeuvre un tel modèle. Celui-ci est codé en Java et l'interface graphique est en HTML. Des expériences sur le simulateur ont mis en évidence deux régimes de dissolution: l'un limité par la diffusion (Diffusion Limited Regime -DLR), l'autre cinétique. Le premier régime dépend de la surface exposée et on y observe l'arrondissement et le lissage progressif des fragments. Le second régime est indépendant de la surface et on observe la formation de cavités dendritiques et une augmentation progressive de la complexité morphologique. D'un point de vue géochimique, le régime DLR est dit «contrôlé par la surface» alors que le régime cinétique est dit «contrôlé par le transport». Les extensions possibles au modèle sont variées et nombreuses. ______________________________________________________________________________ MOTS-CLÉS DE L’AUTEUR : Brèche hydrothermale, Automate cellulaire, Modélisation, Dissolution, Dimension fractale.

La mine Marban est située dans le sud de la ceinture de roches vertes de l'Abitibi, dans le secteur minier de Val d'Or-Malartic. Les roches encaissantes de la mine sont une alternance d'unités mafiques et ultramafiques de la Formation de Jacola. La mine Marban a produit 1 983 112 tonnes à 5.27 g/t d'or (Au) pour un total de 330 027 onces entre 1961 et 1974. Le minerai aurifère de la mine Marban est encaissé par une unité mafique, principalement dans sa partie cisaillée et plissée. Le gisement est situé dans la zone de cisaillement Marbenite, orientée est-sud-est -ouest-nord-ouest. La géologie locale montre une alternance d'unités mafiques et ultramafiques cisaillées et plissées, recoupées par des dykes dioritiques et gabbroïques ainsi que par des intrusions granodioritiques dans la partie nord du secteur. De récents forages carottés ont permis d'établir la géométrie 3D des unités présentes. Les unités mafiques et ultramafiques forment quatre plis isoclinaux serrés dont le plan axial est généralement orienté est-sud est/ouest-nord ouest, et incliné vers le nord de 45° à 60°. Il semble donc que la minéralisation exploitée à la mine Marban soit logée dans les charnières des plis antiformes développés dans une unité de roches mafiques, à proximité du contact avec une unité de roches ultramafiques, ou le long du plan axial de ces plis. L'étude des figures de déformation des roches de Marban montrent que le secteur a subit des déformations cassantes et ductiles, provoquées par l'orogénie Kénoréenne. La forme des minéraux métalliques dont celle des grains d'or suggère que la minéralisation a subit au moins une partie de ces déformations. Il semble donc que le gisement de Marban soit précoce par rapport aux évènements associés à l'orogénie Kénoréenne. D'après ces éléments, l'âge de précipitation de l'or de Marban est pré-à syn-tectonique. Une étude métallogénique et minéralogique des carottes de sondage a permis de déterminer qu'une grande partie de la minéralisation aurifère récemment découverte dans les forages se concentre le long du contact entre les unités de roches mafiques et ultramafiques. On y a par exemple obtenu des teneurs de 6.08 g/t Au sur 6.5 m, et de 2,44 g/t Au sur 9 m sur les flancs du pli antiforme contenant le gisement de la mine Marban. D'autre part, des mesures de la susceptibilité magnétique des lithologies recoupées par les sondages font apparaître une diminution de la teneur en magnétite des roches ultramafiques à l'approche du contact avec les roches mafiques adjacentes. L'hypothèse de la mise en place de l'or est que des fluides enrichis en CO₂ et en complexes bisulfurés porteurs d'or sont associés à une altération des différentes lithologies de Marban, forçant ainsi des échanges métasomatiques. L'altération des roches ultramafiques (komatiites) en roches à talc-carbonate aurait provoqué la transformation de la magnétite en pyrite par interaction avec des fluides riches en CO₂ et en soufre réduit. Des cristaux de magnétite corrodés par de la pyrite ont notamment été observés. Dans ce type d'environnement, la cristallisation de la pyrite entraîne indirectement la précipitation de l'or transporté par les fluides. La minéralisation aurifère apparaît donc comme le produit d'une altération métasomatique entre des unités de composition contrastée, mafique-ultramafique. Cette altération se manifeste souvent sous forme d'auréoles aurifères enrichies en chlorite, carbonates, et pyrite le long des contacts entre les différentes lithologies. ______________________________________________________________________________ MOTS-CLÉS DE L’AUTEUR : Abitibi, Formation de Jacola, Mine Marban, Minéralisation aurifère, Roches basaltiques, Roches ultramafiques, Altération métasomatique, Complexes bisulfurés, Fluides, Pyrite, Magnétite, Chlorites.