Prediction of water–rock interaction and porosity evolution in a granitoid-hosted enhanced geothermal system, using constraints from the 5 km Basel-1 well

Numerical simulations based on plans for a deep geothermal system in Basel, Switzerland are used here to understand chemical processes that occur in an initially dry granitoid reservoir during hydraulic stimulation and long-term water circulation to extract heat. An important question regarding the sustainability of such enhanced geothermal systems (EGS), is whether water–rock reactions will eventually lead to clogging of flow paths in the reservoir and thereby reduce or even completely block fluid throughput. A reactive transport model allows the main chemical reactions to be predicted and the resulting evolution of porosity to be tracked over the expected 30-year operational lifetime of the system. The simulations show that injection of surface water to stimulate fracture permeability in the monzogranite reservoir at 190 °C and 5000 m depth induces redox reactions between the oxidised surface water and the reduced wall rock. Although new calcite, chlorite, hematite and other minerals precipitate near the injection well, their volumes are low and more than compensated by those of the dissolving wall-rock minerals. Thus, during stimulation, reduction of injectivity by mineral precipitation is unlikely. During the simulated long-term operation of the system, the main mineral reactions are the hydration and albitization of plagioclase, the alteration of hornblende to an assemblage of smectites and chlorites and of primary K-feldspar to muscovite and microcline. Within a closed-system doublet, the composition of the circulated fluid changes only slightly during its repeated passage through the reservoir, as the wall rock essentially undergoes isochemical recrystallization. Even after 30 years of circulation, the calculations show that porosity is reduced by only ∼0.2%, well below the expected fracture porosity induced by stimulation. This result suggests that permeability reduction owing to water–rock interaction is unlikely to jeopardize the long-term operation of deep, granitoid-hosted EGS systems. A peculiarity at Basel is the presence of anhydrite as fracture coatings at ∼5000 m depth. Simulated exposure of the circulating fluid to anhydrite induces a stronger redox disequilibrium in the reservoir, driving dissolution of ferrous minerals and precipitation of ferric smectites, hematite and pyrite. However, even in this scenario the porosity reduction is at most 0.5%, a value which is unproblematic for sustainable fluid circulation through the reservoir.

Predictive modelling of mineral scaling, corrosion and the performance of solute geothermometers in a granitoid-hosted, enhanced geothermal system

Experience is lacking with mineral scaling and corrosion in enhanced geothermal systems (EGS) in which surface water is circulated through hydraulically stimulated crystalline rocks. As an aid in designing EGS projects we have conducted multicomponent reactive-transport simulations to predict the likely characteristics of scales and corrosion that may form when exploiting heat from granitoid reservoir rocks at ∼200 °C and 5 km depth. The specifications of an EGS project at Basel, Switzerland, are used to constrain the model. The main water–rock reactions in the reservoir during hydraulic stimulation and the subsequent doublet operation were identified in a separate paper (Alt-Epping et al., 2013b). Here we use the computed composition of the reservoir fluid to (1) predict mineral scaling in the injection and production wells, (2) evaluate methods of chemical geothermometry and (3) identify geochemical indicators of incipient corrosion. The envisaged heat extraction scheme ensures that even if the reservoir fluid is in equilibrium with quartz, cooling of the fluid will not induce saturation with respect to amorphous silica, thus eliminating the risk of silica scaling. However, the ascending fluid attains saturation with respect to crystalline aluminosilicates such as albite, microcline and chlorite, and possibly with respect to amorphous aluminosilicates. If no silica-bearing minerals precipitate upon ascent, reservoir temperatures can be predicted by classical formulations of silica geothermometry. In contrast, Na/K concentration ratios in the production fluid reflect steady-state conditions in the reservoir rather than albite–microcline equilibrium. Thus, even though igneous orthoclase is abundant in the reservoir and albite precipitates as a secondary phase, Na/K geothermometers fail to yield accurate temperatures. Anhydrite, which is present in fractures in the Basel reservoir, is predicted to dissolve during operation. This may lead to precipitation of pyrite and, at high exposure of anhydrite to the circulating fluid, of hematite scaling in the geothermal installation. In general, incipient corrosion of the casing can be detected at the production wellhead through an increase in H2(aq) and the enhanced precipitation of Fe-bearing aluminosilicates. The appearance of magnetite in scales indicates high corrosion rates.

An isotopic analysis of geothermal brine and calcite scaling from the Blue Mountain geothermal field, Winnemucca, Nevada

Understanding, and controlling, the conditions under which calcite precipitates within geothermal energy production systems is a key step in maintaining production efficiency. In this study, I apply methods of bulk and clumped isotope thermometry to an operating geothermal energy facility in northern Nevada to see how those methods can better inform the facility owner, AltaRock Energy, Inc., about the occurrence of calcite scale in their power plant. I have taken water samples from five production wells, the combined generator effluent, shallow cold-water wells, monitoring wells, and surface water. I also collected calcite scale samples from within the production system. Water samples were analyzed for stable oxygen isotope composition (d18O). Calcite samples were analyzed for stable oxygen and carbon (d13C) composition, and clumped isotope composition (D47). With two exceptions, the water compositions are very similar, likely indicating common origin and a well-mixed hydrothermal system. The calcite samples are likewise similar to one another. Apparent temperatures calculated from d18O values of water and calcite are lower than those recorded for the system. Apparent temperatures calculated from D47 are several degrees higher than the recorded well temperatures. The lower temperatures from the bulk isotope data are consistent with temperatures that could be expected during a de-pressurization of the production system, which would cause boiling in the pipes, a reduction in system temperature, and rapid precipitation of calcite scale. However, the high apparent temperature indicated by the D47 data suggests that the calcite is depleted in clumped isotopes given the known temperature of the system, which is inconsistent with this hypothesis. This depletion could instead result from disequilibrium isotopic fractionation during the aforementioned boil events, which would make both the apparent d18O-based and D47-based temperatures unrepresentative of the actual water temperature. This research can help improve our understanding of how isotopic analyses can better inform us about the movement of water through geothermal systems of the past and how it now moves through modern systems. Increased understanding of water movement in these systems could potentially allow for more efficient utilization of geothermal energy as a renewable resource.

Analyses of permeability and porosity of sedimentary rocks in terms of unconventional geothermal resource explorations in Poland

Petrophysical investigations are fundamental to natural resource exploration. In order to recognise the geothermal potential of sedimentary rocks in central Poland, 259 samples were collected from prospective deep-lying geothermal reservoirs. Parameters measured include bulk density, skeletal density, effective porosity, permeability, average pore diameter and specific surface. Results indicate that at great depths (mostly > 3,000 m below surface) sedimentary rocks show low values of porosity (mainly less than 5%) and permeability (only sporadically in excess of 1 md). These values call for a petrothermal use of reservoirs, for which an Enhanced Geothermal System (EGS) was developed. Reser- voirs suited for the EGS are Carboniferous and Lower Triassic sandstones in the central part of Poland (Mogilno-Łódź Trough region and a small part of the Kujawy Swell and Fore-Sudetic regions). In addition, Carboniferous limestones in this area are potentially prospective.

Different mineralization styles in a volcanic-hosted ore deposit: the fluid and isotopic signatures of the Mt Morgan Au-Cu deposit, Australia

Quantitative laser ablation (LA)-ICP-MS analyses of fluid inclusions, trace element chemistry of sulfides, stable isotope (S), and Pb isotopes have been used to discriminate the formation of two contrasting mineralization styles and to evaluate the origin of the Cu and Au at Mt Morgan. The Mt Morgan Au-Cu deposit is hosted by Devonian felsic volcanic rocks that have been intruded by multiple phases of the Mt Morgan Tonalite, a low-K, low-Al2O3 tonalite-trondhjemite-dacite (TTD) complex. An early, barren massive sulfide mineralization with stringer veins is conforming to VHMS sub-seafloor replacement processes, whereas the high-grade Au-Cu. ore is associated with a later quartz-chalcopyrite-pyrite stock work mineralization that is related to intrusive phases of the Tonalite complex. LA-ICP-MS fluid inclusion analyses reveal high As (avg. 8850 ppm) and Sb (avg. 140 ppm) for the Au-Cu mineralization and 5 to 10 times higher Cu concentration than in the fluids associated with the massive pyrite mineralization. Overall, the hydrothermal system of Mt Morgan is characterized by low average fluid salinities in both mineralization styles (45-80% seawater salinity) and temperatures of 210 to 270 degreesC estimated from fluid inclusions. Laser Raman Spectroscopic analysis indicates a consistent and uniform array Of CO2-bearing fluids. Comparison with active submarine hydrothermal vents shows an enrichment of the Mt Morgan fluids in base metals. Therefore, a seawater-dominated fluid is assumed for the barren massive sulfide mineralization, whereas magmatic volatile contributions are implied for the intrusive related mineralization. Condensation of magmatic vapor into a seawater-dominated environment explains the CO2 occurrence, the low salinities, and the enriched base and precious metal fluid composition that is associated with the Au-Cu. mineralization. The sulfur isotope signature of pyrite and chalcopyrite is composed of fractionated Devonian seawater and oxidized magmatic fluids or remobilized sulfur from existing sulfides. Pb isotopes indicate that Au and Cu. originated from the Mt Morgan intrusions and a particular volcanic strata that shows elevated Cu background. (C) 2002 Elsevier Science B.V. All rights reserved.

Time-lapse and probabilistic inversion strategies for plane-wave electromagnetic data.

L'utilisation efficace des systèmes géothermaux, la séquestration du CO2 pour limiter le changement climatique et la prévention de l'intrusion d'eau salée dans les aquifères costaux ne sont que quelques exemples qui démontrent notre besoin en technologies nouvelles pour suivre l'évolution des processus souterrains à partir de la surface. Un défi majeur est d'assurer la caractérisation et l'optimisation des performances de ces technologies à différentes échelles spatiales et temporelles. Les méthodes électromagnétiques (EM) d'ondes planes sont sensibles à la conductivité électrique du sous-sol et, par conséquent, à la conductivité électrique des fluides saturant la roche, à la présence de fractures connectées, à la température et aux matériaux géologiques. Ces méthodes sont régies par des équations valides sur de larges gammes de fréquences, permettant détudier de manières analogues des processus allant de quelques mètres sous la surface jusqu'à plusieurs kilomètres de profondeur. Néanmoins, ces méthodes sont soumises à une perte de résolution avec la profondeur à cause des propriétés diffusives du champ électromagnétique. Pour cette raison, l'estimation des modèles du sous-sol par ces méthodes doit prendre en compte des informations a priori afin de contraindre les modèles autant que possible et de permettre la quantification des incertitudes de ces modèles de façon appropriée. Dans la présente thèse, je développe des approches permettant la caractérisation statique et dynamique du sous-sol à l'aide d'ondes EM planes. Dans une première partie, je présente une approche déterministe permettant de réaliser des inversions répétées dans le temps (time-lapse) de données d'ondes EM planes en deux dimensions. Cette stratégie est basée sur l'incorporation dans l'algorithme d'informations a priori en fonction des changements du modèle de conductivité électrique attendus. Ceci est réalisé en intégrant une régularisation stochastique et des contraintes flexibles par rapport à la gamme des changements attendus en utilisant les multiplicateurs de Lagrange. J'utilise des normes différentes de la norme l2 pour contraindre la structure du modèle et obtenir des transitions abruptes entre les régions du model qui subissent des changements dans le temps et celles qui n'en subissent pas. Aussi, j'incorpore une stratégie afin d'éliminer les erreurs systématiques de données time-lapse. Ce travail a mis en évidence l'amélioration de la caractérisation des changements temporels par rapport aux approches classiques qui réalisent des inversions indépendantes à chaque pas de temps et comparent les modèles. Dans la seconde partie de cette thèse, j'adopte un formalisme bayésien et je teste la possibilité de quantifier les incertitudes sur les paramètres du modèle dans l'inversion d'ondes EM planes. Pour ce faire, je présente une stratégie d'inversion probabiliste basée sur des pixels à deux dimensions pour des inversions de données d'ondes EM planes et de tomographies de résistivité électrique (ERT) séparées et jointes. Je compare les incertitudes des paramètres du modèle en considérant différents types d'information a priori sur la structure du modèle et différentes fonctions de vraisemblance pour décrire les erreurs sur les données. Les résultats indiquent que la régularisation du modèle est nécessaire lorsqu'on a à faire à un large nombre de paramètres car cela permet d'accélérer la convergence des chaînes et d'obtenir des modèles plus réalistes. Cependent, ces contraintes mènent à des incertitudes d'estimations plus faibles, ce qui implique des distributions a posteriori qui ne contiennent pas le vrai modèledans les régions ou` la méthode présente une sensibilité limitée. Cette situation peut être améliorée en combinant des méthodes d'ondes EM planes avec d'autres méthodes complémentaires telles que l'ERT. De plus, je montre que le poids de régularisation des paramètres et l'écart-type des erreurs sur les données peuvent être retrouvés par une inversion probabiliste. Finalement, j'évalue la possibilité de caractériser une distribution tridimensionnelle d'un panache de traceur salin injecté dans le sous-sol en réalisant une inversion probabiliste time-lapse tridimensionnelle d'ondes EM planes. Etant donné que les inversions probabilistes sont très coûteuses en temps de calcul lorsque l'espace des paramètres présente une grande dimension, je propose une stratégie de réduction du modèle ou` les coefficients de décomposition des moments de Legendre du panache de traceur injecté ainsi que sa position sont estimés. Pour ce faire, un modèle de résistivité de base est nécessaire. Il peut être obtenu avant l'expérience time-lapse. Un test synthétique montre que la méthodologie marche bien quand le modèle de résistivité de base est caractérisé correctement. Cette méthodologie est aussi appliquée à un test de trac¸age par injection d'une solution saline et d'acides réalisé dans un système géothermal en Australie, puis comparée à une inversion time-lapse tridimensionnelle réalisée selon une approche déterministe. L'inversion probabiliste permet de mieux contraindre le panache du traceur salin gr^ace à la grande quantité d'informations a priori incluse dans l'algorithme. Néanmoins, les changements de conductivités nécessaires pour expliquer les changements observés dans les données sont plus grands que ce qu'expliquent notre connaissance actuelle des phénomenès physiques. Ce problème peut être lié à la qualité limitée du modèle de résistivité de base utilisé, indiquant ainsi que des efforts plus grands devront être fournis dans le futur pour obtenir des modèles de base de bonne qualité avant de réaliser des expériences dynamiques. Les études décrites dans cette thèse montrent que les méthodes d'ondes EM planes sont très utiles pour caractériser et suivre les variations temporelles du sous-sol sur de larges échelles. Les présentes approches améliorent l'évaluation des modèles obtenus, autant en termes d'incorporation d'informations a priori, qu'en termes de quantification d'incertitudes a posteriori. De plus, les stratégies développées peuvent être appliquées à d'autres méthodes géophysiques, et offrent une grande flexibilité pour l'incorporation d'informations additionnelles lorsqu'elles sont disponibles. -- The efficient use of geothermal systems, the sequestration of CO2 to mitigate climate change, and the prevention of seawater intrusion in coastal aquifers are only some examples that demonstrate the need for novel technologies to monitor subsurface processes from the surface. A main challenge is to assure optimal performance of such technologies at different temporal and spatial scales. Plane-wave electromagnetic (EM) methods are sensitive to subsurface electrical conductivity and consequently to fluid conductivity, fracture connectivity, temperature, and rock mineralogy. These methods have governing equations that are the same over a large range of frequencies, thus allowing to study in an analogous manner processes on scales ranging from few meters close to the surface down to several hundreds of kilometers depth. Unfortunately, they suffer from a significant resolution loss with depth due to the diffusive nature of the electromagnetic fields. Therefore, estimations of subsurface models that use these methods should incorporate a priori information to better constrain the models, and provide appropriate measures of model uncertainty. During my thesis, I have developed approaches to improve the static and dynamic characterization of the subsurface with plane-wave EM methods. In the first part of this thesis, I present a two-dimensional deterministic approach to perform time-lapse inversion of plane-wave EM data. The strategy is based on the incorporation of prior information into the inversion algorithm regarding the expected temporal changes in electrical conductivity. This is done by incorporating a flexible stochastic regularization and constraints regarding the expected ranges of the changes by using Lagrange multipliers. I use non-l2 norms to penalize the model update in order to obtain sharp transitions between regions that experience temporal changes and regions that do not. I also incorporate a time-lapse differencing strategy to remove systematic errors in the time-lapse inversion. This work presents improvements in the characterization of temporal changes with respect to the classical approach of performing separate inversions and computing differences between the models. In the second part of this thesis, I adopt a Bayesian framework and use Markov chain Monte Carlo (MCMC) simulations to quantify model parameter uncertainty in plane-wave EM inversion. For this purpose, I present a two-dimensional pixel-based probabilistic inversion strategy for separate and joint inversions of plane-wave EM and electrical resistivity tomography (ERT) data. I compare the uncertainties of the model parameters when considering different types of prior information on the model structure and different likelihood functions to describe the data errors. The results indicate that model regularization is necessary when dealing with a large number of model parameters because it helps to accelerate the convergence of the chains and leads to more realistic models. These constraints also lead to smaller uncertainty estimates, which imply posterior distributions that do not include the true underlying model in regions where the method has limited sensitivity. This situation can be improved by combining planewave EM methods with complimentary geophysical methods such as ERT. In addition, I show that an appropriate regularization weight and the standard deviation of the data errors can be retrieved by the MCMC inversion. Finally, I evaluate the possibility of characterizing the three-dimensional distribution of an injected water plume by performing three-dimensional time-lapse MCMC inversion of planewave EM data. Since MCMC inversion involves a significant computational burden in high parameter dimensions, I propose a model reduction strategy where the coefficients of a Legendre moment decomposition of the injected water plume and its location are estimated. For this purpose, a base resistivity model is needed which is obtained prior to the time-lapse experiment. A synthetic test shows that the methodology works well when the base resistivity model is correctly characterized. The methodology is also applied to an injection experiment performed in a geothermal system in Australia, and compared to a three-dimensional time-lapse inversion performed within a deterministic framework. The MCMC inversion better constrains the water plumes due to the larger amount of prior information that is included in the algorithm. The conductivity changes needed to explain the time-lapse data are much larger than what is physically possible based on present day understandings. This issue may be related to the base resistivity model used, therefore indicating that more efforts should be given to obtain high-quality base models prior to dynamic experiments. The studies described herein give clear evidence that plane-wave EM methods are useful to characterize and monitor the subsurface at a wide range of scales. The presented approaches contribute to an improved appraisal of the obtained models, both in terms of the incorporation of prior information in the algorithms and the posterior uncertainty quantification. In addition, the developed strategies can be applied to other geophysical methods, and offer great flexibility to incorporate additional information when available.

La caratterizzazione probabilistica del sottosuolo come strumento per l’ottimizzazione della progettazione integrata dei sistemi geotermici

L’utilizzo del reservoir geotermico superficiale a scopi termici / frigoriferi è una tecnica consolidata che permette di sfruttare, tramite appositi “geoscambiatori”, un’energia presente ovunque ed inesauribile, ad un ridotto prezzo in termini di emissioni climalteranti. Pertanto, il pieno sfruttamento di questa risorsa è in linea con gli obiettivi del Protocollo di Kyoto ed è descritto nella Direttiva Europea 2009/28/CE (Comunemente detta: Direttiva Rinnovabili). Considerato il notevole potenziale a fronte di costi sostenibili di installazione ed esercizio, la geotermia superficiale è stata sfruttata già dalla metà del ventesimo secolo in diversi contesti (geografici, geologici e climatici) e per diverse applicazioni (residenziali, commerciali, industriali, infrastrutturali). Ciononostante, solo a partire dagli anni 2000 la comunità scientifica e il mercato si sono realmente interessati ed affacciati all’argomento, a seguito di sopraggiunte condizioni economiche e tecniche. Una semplice ed immediata dimostrazione di ciò si ritrova nel fatto che al 2012 non esiste ancora un chiaro riferimento tecnico condiviso a livello internazionale, né per la progettazione, né per l’installazione, né per il testing delle diverse applicazioni della geotermia superficiale, questo a fronte di una moltitudine di articoli scientifici pubblicati, impianti realizzati ed associazioni di categoria coinvolte nel primo decennio del ventunesimo secolo. Il presente lavoro di ricerca si colloca all’interno di questo quadro. In particolare verranno mostrati i progressi della ricerca svolta all’interno del Dipartimento di Ingegneria Civile, Ambientale e dei Materiali nei settori della progettazione e del testing dei sistemi geotermici, nonché verranno descritte alcune tipologie di geoscambiatori innovative studiate, analizzate e testate nel periodo di ricerca.

Modelización y simulación numérica de sistemas geotérmicos planos

En el presente proyecto se estudiará la conveniencia de utilizar simulación numérica para evaluar, diseñar y mejorar el diseño de instalaciones geotérmicas planas. Para ello, se modelizará el conjunto de terreno y sistema geotérmico mediante el acoplamiento de un modelo tridimensional (terreno) y un modelo unidimensional, formado por la red de tubería geotérmica, que estará dispuesta en un plano cuya orientación dependerá de las aplicaciones. Para la simulación numérica se utilizará un software de análisis y resolución por elementos finitos tipo comercial. Finalmente, se definirán unos parámetros geométricos de diseño y, teniendo en cuenta las propiedades térmicas de cada material, se realizarán estudios paramétricos para determinar la influencia de cada grupo de parámetros. Ello requerirá la realización de cientos de simulaciones numéricas que permitirán optimizar el rendimiento de sistemas geotérmicos planos tanto horizontales (soleras), como verticales (muros pantalla). ABSTRACT In this project it is going to be studied the convenience of numerical simulation as to evaluate, design and improve the design of geothermal planar installations. With this purpose, the terrain and the geothermal system will be modeled by coupling a three-dimensional model (terrain) and a one-dimensional model formed by the geothermal piping network disposed in a plane with an orientation which will depend on the applications. For the numerical simulation commercial software of analyzing and resolving by finite elements will be used. Finally, geometric design parameters will be defined, and taking into account the thermal properties of each material, parametric studies will be performed to determine the influence of each group of parameters. This will require the completion of hundreds of numerical simulations which will enhance the performance of both horizontal (slabs) and vertical (diaphragm walls) geothermal systems.

Rare earth element fractionation in authigenic illite–smectite from Late Permian clastic rocks, Bowen Basin, Australia: implications for physico-chemical environments of fluids during illitization

REE analyses were performed on authigenic illitic clay. minerals from Late Permian mudrocks, sandstones and bentonites from the Bowen Basin (Australia). The mixed-layer illite-smectite exhibit REE patterns with an obvious fractionation of the HREE from the LREE and MREE, which is an apparent function of degree of illitization reaction. The highly illitic (R greater than or equal to 3) illite-smectite from the northern Bowen Basin show a depletion of LREE relative to the less illitic (R=0 and 1) clays. In contrast, an enrichment of HREE for the illite-rich clays relative to less. illitic clays is evident for the southern Bowen Basin samples. The North American Shale Composite-normalized (La/Lu)(sn) ratios show negative correlations with the illite content in illite-smectite and positive correlations with the delta(18)O values of the clays for both the northern and southern Bowen Basin samples. These correlations indicate that the increasing depletion of LREE in hydrothermal fluids is a function of increasing water/rock ratios in the northern Bowen Basin. Good negative correlations between (La/Lu)(sn) ratios and illite content in illite-smectite from the southern Bowen Basin suggest the involvement of fluids with higher alkalinity and higher pH in low water/ rock ratio conditions. Increasing HREE enrichment with delta(18)O decrease indicates the effect of increasing temperature at low water/rock ratios in the southern Bowen Basin. Results of the present study confirm the conclusions of some earlier studies suggesting that REE in illitic clay minerals are mobile and fractionated during illitization and that this fact should be considered in studies of sedimentary processes and in identifying provenance. Moreover, our results show that REE systematic of illitic clay minerals can be applied as an useful technique to gain information about physico-chemical conditions during thermal and fluid flow events in certain sedimentary basins. (C) 2003 Elsevier Science B.V. All rights reserved.

Bajo de la Alumbrera copper-gold deposit: Stable isotope evidence for a porphyry-related hydrothermal system dominated by magmatic aqueous fluids

Alteration zones at the gold-rich Bajo de la Alumbrera porphyry copper deposit in northwestern Argentina are centered on several porphyritic intrusions. They are zoned from a central copper-iron sulfide and gold-mineralized potassic (biotite-K-feldspar +/- quartz) core outward to propylitic (chlorite-illite-epidote-calcite) assemblages. A mineralized intermediate argillic alteration assemblage (chlorite-illite +/- pyrite) has overprinted the potassic alteration zone across the top and sides of the deposit and is itself zoned outward into phyllic (quartzinuscovite-illite +/- pyrite) alteration. This study contributes new data to previously reported delta(18)O and delta D compositions of fluids responsible for the alteration at Bajo de la Alumbrera, and the data are used to infer likely ore-forming processes. Measured and calculated delta(18)O and delta D values of fluids (+8.3 to +10.2 and -33 to -81 parts per thousand, respectively) confirm a primary magmatic origin for the earliest potassic alteration phase. Lower temperature potassic alteration formed from magmatic fluids with lower delta D values (down to -123 parts per thousand). These depleted compositions are distinct from meteoric water and consistent with degassing and volatile exsolution of magmatic fluids derived from an underlying magma. Variability in the calculated composition of fluid associated with potassic alteration is explained in terms of phase separation (or boiling). if copper-iron sulfide deposition occurred during cooling (as proposed elsewhere), this cooling was largely a result of phase separation. Magmatic water was directly involved in the formation of overprinting intermediate argillic alteration assemblages at Bajo de la Alumbrera. Calculated delta(18)O and delta D values of fluids associated with this alteration range from +4.8 to +8.1 and -31 to -71 per mil, respectively Compositions determined for fluids associated with phyllic alteration (-0.8 to +10.2 and -31 to -119 parts per thousand) overlap with the values determined for the intermediate argillic alteration. We infer that phyllic alteration assemblages developed during two stages; the first was a high-temperature (400 degrees-300 degrees C) stage with D-depleted water (delta D = -66 to -119 parts per thousand). This compositional range may have resulted from magma degassing and/or the injection of new magmatic water into a compositionally evolved hydrothermal system. The isotopic variations also can be explained by increased fluid-rock interaction. The second stage of phyllic alteration occurred at a lower temperature (similar to 200 degrees C), and variations in the modeled isotopic compositions imply mixing of magmatic and meteoric waters. Ore deposition that occurred late in the evolution of the hydrothermal system was probably associated with further cooling of the magmatic fluid, in part caused by fluid-rock interaction and phase separation. Changing pH and/or oxygen fuoracity may have caused additional ore deposition. The ingress of meteoric water appears to postdate the bulk of mineralization and occurred as the system at Bajo de la Alumbrera waned.

(Supplement 1) Selected chemical composition of Icelandic river waters and bedrock age

In this study we map the spatial distribution of selected dissolved constituents in Icelandic river waters using GIS methods to study and interpret the connection between river chemistry, bedrock, hydrology, vegetation and aquatic ecology. Five parameters were selected: alkalinity, SiO2, Mo, F and the dissolved inorganic nitrogen and dissolved inorganic phosphorus mole ratio (DIN/DIP). The highest concentrations were found in rivers draining young rocks within the volcanic rift zone and especially those draining active central volcanoes. However, several catchments on the margins of the rift zone also had high values for these parameters, due to geothermal influence or wetlands within their catchment area. The DIN/DIP mole ratio was higher than 16 in rivers draining old rocks, but lowest in rivers within the volcanic rift zone. Thus primary production in the rivers is limited by fixed dissolved nitrogen within the rift zone, but dissolved phosphorus in the old Tertiary catchments. Nitrogen fixation within the rift zone can be enhanced by high dissolved molybdenum concentrations in the vicinity of volcanoes. The river catchments in this study were subdivided into several hydrological categories. Importantly, the variation in the hydrology of the catchments cannot alone explain the variation in dissolved constituents. The presence or absence of central volcanoes, young reactive rocks, geothermal systems and wetlands is important for the chemistry of the river waters. We used too many categories within several of the river catchments to be able to determine a statistically significant connection between the chem¬istry of the river waters and the hydrological categories. More data are needed from rivers draining one single hydrological category. The spatial dissolved constituent distribution clearly revealed the difference between the two extremes, the young rocks of the volcanic rift zone and the old Tertiary terrain.

Coupled wellbore-reservoir simulation of thermal effects during multiphase CO2 injection

T2Well-ECO2M is a coupled wellbore reservoir simulator still under development at Lawrence Berkeley National Laboratory (USA) with the ability to deal with a mixture of H2O-CO2-NaCl and includes the simulation of CO2 phase transition and multiphase flow. The code was originally developed for the simulation of CO2 injection into deep saline aquifers and the modelling of enhanced geothermal systems; however, the focus of this research was to modify and test T2Well-ECO2M to simulate CO2 injection into depleted gas reservoirs. To this end, the original code was properly changed in a few parts and a dedicated injection case was developed to study CO2 phase transition inside of a wellbore and the corresponding thermal effects. In the first scenario, the injection case was run applying the fully numerical approach of wellbore to formation heat exchange calculation. Results were analysed in terms of wellbore pressure and temperature vertical profiles, wellhead and bottomhole conditions, and characteristic reservoir displacement fronts. Special attention was given to the thorough analysis of bottomhole temperature as the critical parameter for hydrate formation. Besides the expected direct effect of wellbore temperature changes on reservoir conditions, the simulation results indicated also the effect of CO2 phase change in the near wellbore zone on BH pressure distribution. To test the implemented software changes, in a second scenario, the same injection case was reproduced using the improved semi-analytical time-convolution approach for wellbore to formation heat exchange calculation. The comparison of the two scenarios showed that the simulation of wellbore and reservoir parameters after one year of continuous CO2 injection are in good agreement with the computation time to solve the time-convolution semi-analytical reduced. The new updated T2Well-ECO2M version has shown to be a robust and performing wellbore-reservoir simulator that can be also used to simulate the CO2 injection into depleted gas reservoirs.