Porosity evolution and fluid flow in the basalts of the Skaergaard magma-hydrothermal system, East Greenland

1991 ◽  
Vol 291 (3) ◽  
pp. 201-257 ◽  
Author(s):  
C. E. Manning ◽  
D. K. Bird
2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Anna Jentsch ◽  
Walter Duesing ◽  
Egbert Jolie ◽  
Martin Zimmer

AbstractCarbon dioxide is the most abundant, non-condensable gas in volcanic systems, released into the atmosphere through either diffuse or advective fluid flow. The emission of substantial amounts of CO2 at Earth’s surface is not only controlled by volcanic plumes during periods of eruptive activity or fumaroles, but also by soil degassing along permeable structures in the subsurface. Monitoring of these processes is of utmost importance for volcanic hazard analyses, and is also relevant for managing geothermal resources. Fluid-bearing faults are key elements of economic value for geothermal power generation. Here, we describe for the first time how sensitively and quickly natural gas emissions react to changes within a deep hydrothermal system due to geothermal fluid reinjection. For this purpose, we deployed an automated, multi-chamber CO2 flux monitoring system within the damage zone of a deep-rooted major normal fault in the Los Humeros Volcanic Complex (LHVC) in Mexico and recorded data over a period of five months. After removing the atmospheric effects on variations in CO2 flux, we calculated correlation coefficients between residual CO2 emissions and reinjection rates, identifying an inverse correlation of ρ = − 0.51 to − 0.66. Our results indicate that gas emissions respond to changes in reinjection rates within 24 h, proving an active hydraulic communication between the hydrothermal system and Earth’s surface. This finding is a promising indication not only for geothermal reservoir monitoring but also for advanced long-term volcanic risk analysis. Response times allow for estimation of fluid migration velocities, which is a key constraint for conceptual and numerical modelling of fluid flow in fracture-dominated systems.


2021 ◽  
Vol 9 (1) ◽  
Author(s):  
Drew L. Siler ◽  
Jeff D. Pepin ◽  
Velimir V. Vesselinov ◽  
Maruti K. Mudunuru ◽  
Bulbul Ahmmed

AbstractIn this paper, we present an analysis using unsupervised machine learning (ML) to identify the key geologic factors that contribute to the geothermal production in Brady geothermal field. Brady is a hydrothermal system in northwestern Nevada that supports both electricity production and direct use of hydrothermal fluids. Transmissive fluid-flow pathways are relatively rare in the subsurface, but are critical components of hydrothermal systems like Brady and many other types of fluid-flow systems in fractured rock. Here, we analyze geologic data with ML methods to unravel the local geologic controls on these pathways. The ML method, non-negative matrix factorization with k-means clustering (NMFk), is applied to a library of 14 3D geologic characteristics hypothesized to control hydrothermal circulation in the Brady geothermal field. Our results indicate that macro-scale faults and a local step-over in the fault system preferentially occur along production wells when compared to injection wells and non-productive wells. We infer that these are the key geologic characteristics that control the through-going hydrothermal transmission pathways at Brady. Our results demonstrate: (1) the specific geologic controls on the Brady hydrothermal system and (2) the efficacy of pairing ML techniques with 3D geologic characterization to enhance the understanding of subsurface processes.


2001 ◽  
Vol 172 (4) ◽  
pp. 427-436 ◽  
Author(s):  
Philippe Gouze ◽  
Riad Hassani ◽  
Dominique Bernard ◽  
Anne Coudrain-Ribstein

Abstract We propose a model for simulating the changes in porosity and permeability caused by hydrothermal diagenesis in sedimentary aquifer where salinity, temperature and fluid flow vary in space and time. Such modifications of the hydrodynamic properties of the medium are bounded to geochemical reactions and groundwater flow. Fluid velocity is particularly low in deep reservoirs (typically less than 1 m/year). Then, the local equilibrium simplification, which is justified by a set of world-wide data of the chemical composition of groundwater, can be implemented toward straightforward transient calculations. In the model presented here, the coupled processes of fluid flow, temperature and chemical species transport are solved using well established methods. The originality of the model is the development carried on to predict the permeability evolution controlled by the mineral dissolution and precipitation. Usually to simulate permeability changes modelers use the classical porosity-permeability model based on statistical analyses of in situ or laboratory measurements. However, hydraulic conductivity changes are not controlled solely by porosity changes, but also depend on pore-scale structure transformations. Depending on the mineral type, the precipitation or dissolution of the same quantity of volumetric quantity will induce very different changes in the hydraulic conductivity. Principally clay minerals depict a wide range of atypical organisations of different microstructural characteristics of the porous media. The spatial distribution of these characteristics cannot be modelled at basin scale. Away from both too complicated and too unrealistically simplified approach, the model presented here is based on the calculation of the permeability evolution from the change in the mineral fraction due to mineral precipitation and dissolution. To simplify, the minerals are divided into two groups: clay minerals and non-clay minerals. The specific contribution of clay minerals is controlled by a single weighting coefficient. This coefficient is associated to the proportion of poorly connected porosity that characterize clay structure, albeit it is presently impossible to propose any quantitative relationship between the value of this parameter and the microstructural characteristics of the diagenetic clays. The model is tested here to simulate the evolution of the porosity and the permeability in a peculiar zone of the Paris Basin. The study area of several hundred meters large is inside the Dogger aquifer, close to the Bray fault zone where invasion of saline water from Triassic formation takes place. This zone is characterised by high thermal and salinity gradient as well as by the superposition of sub-horizontal regional flow and ascendant fault-controlled flow: it is an ideal case study for examining the importance of taking into account the specific contribution by clay minerals when computing permeability evolution. This study is proposed as a parameter sensibility analysis: - to compare the relative influence of the clay weighting coefficient, the temperature, the salinity, and the cementation exponent on the computed evolution of the permeability, - to discuss the consequences of the introduction of the clay weighting coefficient in comparison to the classical porosity - permeability evolution model, - to simulate various evolution scenarios of past and future thermal and geochemical constraints and their consequences on the evolution of the permeability changes in the Bray fault zone taking into account uncertainties on the value of the clay weighting coefficient and on the cementation exponent. Forty-one simulations of one million years were necessary to cover a large spectrum of the expected variations of each parameter. The results show that: - the local variation of the permeability depends on the time evolution of temperature and of salinity, and on the values of the cementation exponent of the porosity-permeability law and of the clay weighting coefficient. Within reasonable ranges of these four parameters, their influence on the permeability changes is of the same order of magnitude, - the influence of the clay weighting coefficient on the porosity evolution is negligible. Feedback effects of permeability evolution on the porosity evolution, through the change in the flow regime, is minor, - by the use of a classical model without a clay weighting coefficient, permeability and porosity present the same pattern of evolution: they both increase or decrease. By the use of the clay weighting coefficient, in some places the permeability and porosity can show opposite evolution. One increases when the other decreases even for low values of the coefficient, - in the vicinity of the fault, the model predict an increase of permeability independently of potential temperature and salinity modifications and whatever the clay mineral weighting coefficient is: Bray fault sealing is unlikely as long as head gradient is maintained in the fracture zone.


2007 ◽  
Vol 79 (3) ◽  
pp. 503-517 ◽  
Author(s):  
Everton M. Bongiolo ◽  
Daniela E. Bongiolo ◽  
Paul Sardini ◽  
André S. Mexias ◽  
Marja Siitari-Kauppi ◽  
...  

This work is an application of the 14C-Polymethylmethacrylate method to compare the porosity evolution between unaltered and propylitic-altered granites, using samples from Lavras do Sul region, Brazil. This method, when coupled with optical and electronic petrography has the advantage over other methods to provide the quantification and identification of total and local porosity of rocks. From petrographic observations, different kinds of porous zones were identified and quantified (microfractures, grain boundaries, alteration of minerals, etc). Results show that unaltered granites have 0.5 to 0.6% porosity and propylitic-altered ones have 1.7 to 1.8% porosity, even between samples with different textures. Porosity of altered rocks increases mainly due to higher porosity of neoformed chlorite, calcite, sericite and microfractures. Field observations show that later phyllic alteration halos are wider in equigranular than in porphyritic granites, which could not be explained by different original porosity between those rocks. The observed differences of phyllic halos diffusion were controlled by structural and fluid/rock ratio variations between the equigranular and porphyritic granitic facies during the later hydrothermal stage.


Sign in / Sign up

Export Citation Format

Share Document