L'extraction de l'or se fait par cyanuration depuis plus d'un siècle. Ce procédé est donc bien connu. Le document qui suit présente la problématique reliée à l'effet de la granulométrie sur la cyanuration de l'or. Plus spécifiquement, l'effet de l'historique du broyage, soit l'histoire relié à la préparation d'un échantillon sur la mise en solution de différentes classes granulométriques, sera étudié dans ce document. Neuf tests ont été fait avec trois temps de broyage, 15, 35 et 65 minutes, et avec trois temps de cyanuration différents, 6, 12 et 24 heures. Les calculs de rendements sont basés sur trois échantillons broyés pendant 15, 35 et 65 minutes sans cyanuration. L'or a été analysé dans les solides et dans les solutions à différents moments en cours de réaction. L'argent, le cuivre et le soufre ont seulement été analysés dans les résidus solides. Les rendements sont calculés à partir des données brutes et des données réconciliées. Avec le traitement utilisé pour les essais, les rendements obtenus sont inférieurs à ceux observés en industrie. Certaines observations laissent penser qu'un pré-traitement d'aération et qu'un prétraitement au dinitrate de plomb auraient dû être effectué avant la cyanuration. Toutefois, l'analyse des résultats montre que la distribution granulométrique de l'or change avec la taille du minerai, l'or ayant tendance à demeurer grossier alors que la finesse du minerai augmente. Les résultats montrent aussi que le comportement de l'or dans une classe granulométrique semble être indépendant de l'histoire de la production de la classe granulométrique. Le traitement des données a permis de calibrer des modèles de dissolution de l'or qui ont servis à simuler l'effet de la granulométrie sur la cyanuration du minerai. / The gold extraction is made by cyanidation from more than one century. So it is a well known technique. The following document presents the problematic related to the effect of the particles' size on gold cyanidation. Specifically, the effect of the mill history, being the story related to the preparation of the ore sample, on the recovery in different particles' size is studied in this document. Nine tests have been made with three times of mill, 15, 35 and 65 minutes, and three different times of cyanidation, 6, 12 and 24 hours. The calculations of the recovery are based on three ore samples that have been milled for 15, 35 and 65 minutes and haven't been leached. The gold have been analysed in the ore and in the solutions at different moments trough the cyanidations. Silver, copper and sulphur have been analysed only in the residual ore. The recoveries are calculated with the raw and the reconciliated data. With the treatment used for the tests, the recoveries results obtained are smaller than the recovery observed in the industry. Some observations reveal that an aeration and lead dinitrate pre-treatment should have been done before the cyanidations. However, the analyses of the results show that the particles' size distribution changes with the size of the ore, since gold tends to stay in the bigger particles wile other particles get smaller. The results also show that the gold's behaviour seems to be independent of the history of the particles' size distribution production. The data treatment gives the possibility to calibrate dissolution models of gold that have been use to simulate the effect of particles' size on the gold cyanidation.

Les modèles à milieu poreux équivalent ou à double continuum sont couramment utilisés lors de la modélisation du transport de masse à l'échelle régionale lorsqu'aucune fracture discrète n'est prise en compte. Par contre, la migration d'un contaminant lors de la modélisation du transport de masse de domaines régionaux où des plans de litage ou des failles majeures sont présents peut être importante. Ainsi, de tels chemins préférentiels doivent être explicitement pris en compte lors de modélisations. Les modèles régionaux devraient donc être autant capable de capter le mouvement du contaminant par advection dans les fractures que la diffusion matricielle dans le milieu poreux. Par contre, le ratio entre l'ouverture d'une fracture, qui est normalement de l'ordre du centimètre au maximum, par rapport à l'étendue latérale de celle-ci est extrêmement grand. L'ordre de grandeur des concentrations rencontrées dans le processus de diffusion matricielle par rapport au transport advection est aussi très grand étant donné le développement de gradient de concentration important à l'intérieur du milieu poreux. Ainsi, étant donné ces différences d'ordre de grandeur, le développement d'un modèle numérique robuste, précis et permettant la génération du maillage adéquat en présence de domaines régionaux fracturés est extrêmement intéressant. Présentement, les modèles permettant de simuler de tels domaines requièrent des maillages extrêmement denses et des temps de calculs dépassant les limites acceptables. L'objectif du projet est donc de développer une méthode permettant de simuler avec précision et efficacité le transport de masse en milieu fracturé à l'échelle régionale, c'est-à-dire sur des distances d'au moins 100 mètres et plus. Ainsi, un nouveau modèle numérique a été développé pour la simulation à l'échelle régionale du transport de masse dans un milieu poreux à fractures discrètes où la diffusion matricielle et l'advection dans les fractures sont considérées. Un code générique d'éléments finis, MEF++, a été utilisé pour résoudre les équations d'écoulement de l'eau et de transport de masse en employant des éléments de fractures 2D et des éléments 3D pour le milieu poreux. Le nouveau modèle résout séparément les équations pour les fractures et pour le milieu poreux qui sont considérés comme deux domaines distincts. Ensuite, ces deux domaines sont couplés par l'entremise d'un terme de transfert de fluide et de masse. Ce terme de transfert est calculé à partir des gradients de fluide et de concentration pouvant exister entre le milieu poreux et la fracture, et il est pondéré par un coefficient dépendant des propriétés du milieu poreux tel que la porosité et le coefficient de diffusion. Cette approche permet une représentation générale du transfert entre fractures et milieu poreux. De plus, pour combler les problèmes de discrétisation du milieu fracturé, MEF++ permet, à l'aide d'un algorithme d'adaptation de maillage utilisant un estimateur d'erreur a posteriori, d'augmenter efficacement la densité d'éléments aux endroits où les gradients de charge hydraulique ou de concentration sont élevés. À partir de ce nouveau modèle numérique, de concert avec l'utilisation de l'adaptation de maillage, des vérifications et des exemples illustratifs ont été effectués pour démontrer les capacités du modèle. Il a été montré que le modèle permet de reproduire, avec une grande précision, les résultats d'une solution analytique connue. De plus, des simulations de vérification ont permis de confirmer l'exactitude de l'adaptation de maillage. Enfin, un domaine ressemblant au site de Smithville, Ontario, contaminé par des BPC, a été utilisé afin d'effectuer l'analyse de sensibilité de certains paramètres sur le modèle numérique. / Regional scale modelling of solute transport in fractured porous media is commonly done with equivalent porous medium or dual-continuum models, where discrete fractures are not explicitly accounted for. However, for regional scale transport over distances of hundreds of meters or more, discrete features such as extensive bedding planes or regional scale faults could have a profound impact on solute migration. These discrete features represent preferential pathways for fluid flow and should ideally be explicitly accounted for in modelling studies instead of relying on a continuum approach. Models should be able to capture the rapid advection along these regional scale features as well as solute diffusion into the surrounding matrix. However, the aspect ratio of such conductive features is extremely high, with typical fracture aperture on the order of centimeters at the most, compared to possible lateral extent on the order of hundreds of meters. The length-scale of diffusion into the matrix is also much less than that of regional scale advective transport, with large concentration gradients developing over very short distances into the matrix. Because of these discrepancies in scales, the generation of a robust and adequate mesh for the numerical simulation of regional scale transport in discrete features is extremely challenging. Currently, models that simulate such Systems need extremely dense grids and computational time that often preclude their use. The objective of this research is to develop an efficient and accurate method for the simulation of regional scale (100 meters and greater) solute transport in the presence of discrete high-conductivity features and diffusion into the matrix. Therefore, a new numerical model has been developed for regional scale mass transport in discretely fractured porous media where matrix diffusion and fractures advection are accounted for. A generic finite elements code, MEF++, had been used to solve the flow and mass transport equations using 2D elements for the fracture features and 3D elements for the porous medium. In the model, fractures are represented as 2D features that are embedded into a 3D porous matrix. The model solves independently the equations for the fracture System and the porous medium System. Then, these two Systems and coupled together using mass and fluid transfer terms. The transfer term is calculated using the concentration or head gradients that exist at the fracture / porous medium interface and it is weighted using a coefficient that is dependant of the fracture and porous medium properties such as the porosity and diffusion coefficient. This approach allows a general representation of the mass and fluid transfer between the fractures features and the porous medium. Also, to prevent discretization problems of the porous medium, MEF++ allows, using an adaptive mesh refinement method that uses an a posteriori error estimator, to increase the mesh density where high gradients are occurring. Using this new numerical model, coupled with the adaptive mesh refinement technique, verifications and illustrative examples have been done to show the model performance and capacities. It is shown that the model allows to reproduce, with great accuracy, the results of a known and accepted analytical solution. Some examples are conducted to show the accuracy and the adaptive mesh refinement method as well as a sensitivity analysis on a hypothetical fractured rock aquifer that resembles the Smithville, Ontario, fractured carbonate aquifer where groundwater contamination by PCB has been recorded and where previous modelling effort were based on continuum approaches.

Les gîtes du Dôme de Lemieux et de Sainte-Marguerite sont deux exemples de minéralisation épithermale et orogénique mise en place dans les roches sédimentaires et volcaniques siluro-dévoniennes de Gaspésie. Durant le siluro-dévonien, la Ceinture de Gaspé est interprétée comme un environnement d’arrière-arc, passant ensuite en régime de transpression lors de la déformation acadienne. Les indices du Dôme de Lemieux montrent des caractéristiques minéralogiques typiques des gisements épithermaux neutres à intermédiaires. Les veines minéralisées sont encaissées dans les roches sédimentaires siluro-dévoniennes mais leur mise en place est liée à la présence de volcanites bimodales dévoniennes. La minéralisation est probablement contemporaine du volcanisme et donc d’âge Dévonien inférieur à moyen, antérieur à la déformation acadienne. Le gîte de Sainte-Marguerite, associé à la faille acadienne de Sainte-Florence, est encaissé dans les roches volcaniques bimodales de Sainte-Marguerite. La minéralisation est représentée par cinq stades de minéralisation : (1) un stade à calcite-chlorite; (2) un stade à quartz-pyrite; (3) un stade à quartz-sphalérite-galène-chalcopyrite ± arsénopyrite-adulaire reprenant le stade 2; (4) un stade à quartz-arsénopyrite ± pyrite-sphalérite; et (5) un stade à quartz-carbonate. Le stade 2 à quartz-pyrite présente des sulfures recristallisés, indiquant que ce stade est précoce et a été affecté par le métamorphisme acadien. La pyrite et le quartz de stade 2 sont ensuite cataclasés puis cimentés par la sphalérite, la galène et la chalcopyrite de stade 3. L’attitude des veines de stade 3 et 4 est en accord avec la déformation régionale acadienne du Dévonien moyen à tardif (398 à 359 Ma) et suggère une affinité orogénique épizonale des stades 3 et 4, en accord avec les données isotopiques et les inclusions fluides. Les isotopes de l’oxygène et du carbone ainsi que la microthermométrie indiquent que les deux gîtes ont été formés à partir de saumures de bassins à des températures comprises entre 150° et 200°C et que de l’eau météorique a participé à la formation des veines du Dôme de Lemieux. Les isotopes du soufre et la composition isotopique du plomb suggèrent que le soufre et les métaux sont issus du lessivage des roches sédimentaires et volcaniques environnantes. La présence de ces gîtes dans la ceinture siluro-dévonienne de Gaspé démontre un potentiel minéral épithermal pré-orogenèse acadienne et orogénique épizonal syn-orogenèse acadienne dans les roches siluro-dévoniennes des Appalaches. / The Lemieux Dome and Sainte-Marguerite deposits are examples of epithermal and orogenic mineralization in Silurian-Devonian sedimentary and volcanic rocks of the Gaspé Peninsula. During the Silurian-Devonian, the Gaspé Belt is interpreted as a back-arc setting that eventually transformed into a transpressive setting during the Acadian orogeny. Lemieux Dome veins display mineralogical features typical of low- to intermediate-sulfidation epithermal deposits. Mineralized veins hosted by Silurian-Devonian sedimentary rocks are related to the emplacement of bimodal Devonian volcanic rocks. Mineralisation is probably coeval with volcanic rocks eruption (lower to middle Devonian), before the Acadian orogeny. The Sainte-Marguerite deposit, associated to the Acadian Sainte-Florence fault, is hosted by the Sainte-Marguerite bimodal volcanic rocks. Mineralization forms in 5 stages: (1) a calcite-chlorite stage; (2) a quartz-pyrite stage; (3) a quartz-sphalerite-galena-chalcopyrite ± arsenopyrite-adularia stage; (4) a quartz-arsenopyrite-pyrite-sphalerite stage and; (5) a quartz-carbonate stage. Stage 2 pyrite and quartz are recrystallized, indicating that sulfides formed early and were affected by Acadian metamorphism. Stage 2 pyrite and quartz were then cataclased and cemented by Stage 3 sphalerite, galena and chalcopyrite. Structural attitude of both Stage 3 and 4 veins is consistent with regional Acadian deformation during the middle to late Devonian (398 to 359 Ma), and suggests an epizonal orogenic affinity for Stage 3 and 4 mineralization, compatible with isotope and fluid inclusion data. δ18O, δ13C and microthermometric data in both deposits indicate that they were formed by infiltration of basinal brines at temperatures between 150 and 200°C. Lemieux Dome veins formed by mixing of the basinal brine with meteoric water. δ34S values and the lead isotope composition suggest that sulfur and metals were leached from surrounding sedimentary and volcanic rocks. Formation of both deposits in the Silurian-Devonian Gaspé Belt shows a mineral potential for pre-Acadian epithermal deposits, and for syn-Acadian epizonal orogenic deposits in other Silurian-Devonian rocks of the Appalachians.

La télédétection des conditions physiques du pergélisol dans trois buttes de pergélisol près d'Umiujaq au Nunavik (Québec) a été réalisée à l'aide de la tomographie électromagnétique forage-forage et surface-forage. Les objectifs de ces travaux géophysiques étaient d'étudier la propagation des ondes électromagnétiques dans un milieu physique partiellement gelé et de déterminer la cryostratigraphie des buttes de pergélisol. Des antennes à géométrie cylindrique de géoradar à ondes courtes de 50, 100 et 200 MHz ont été glissées dans des regards géophysiques secs mis en place préalablement à l'étude ou utilisées en surface. La cryostratigraphie générale des buttes de pergélisol étudiées définie à partir des tomogrammes de vitesse de propagation des ondes électromagnétiques correspond à celle obtenue de forages et de l'échantillonnage de ces buttes. La résolution des tomogrammes est meilleure aux plus hautes fréquences.

Une géomorphologie complexe avec d'importantes structures de mouvements de masse sous-marins a été mise en évidence lors de levés géophysiques dans l'estuaire du St-Laurent entre Betsiamites et Rimouski, Québec. Cette étude présente une description morpho-sédimentologique du secteur, une analyse spatio-temporelle des mouvements de masse et finalement une chronologie proposant l'âge pour certaines des ruptures de pente répertoriées. Les structures observées sont associées à au moins quatre épisodes de mouvement de masse. Un premier mouvement de masse paraglaciaire daterait de 9 ka BP. Un second aurait causé de larges couloirs de glissement sur le plateau ainsi qu'un cône de sédiments dans le Chenal laurentien et daterait de 7.25 ka cal BP. Un troisième événement plus récent aurait été causé par le séisme de 1663 AD (M~7) et serait associé au glissement subaérien identifié sur la côte. Finalement une dernière coulée daterait de l'an 1860 AD ou 1870 AD.