scholarly journals Permeability Changes Resulting from Quartz Precipitation and Dissolution around Upper Crustal Intrusions

Geofluids ◽  
2018 ◽  
Vol 2018 ◽  
pp. 1-19 ◽  
Author(s):  
Samuel W. Scott ◽  
Thomas Driesner

It has long been recognized that quartz precipitation from circulating hydrothermal fluids may reduce porosity and permeability near intrusions. However, the magnitude of permeability changes and potential feedbacks between flow, heat transfer, and quartz precipitation/dissolution remain largely unquantified. Here, we present numerical simulations of fluid convection around upper crustal intrusions which explicitly incorporate the feedback between quartz solubility and rock permeability. As groundwater is heated to ~350°C, silica dissolves from the host rock, increasing porosity and permeability. Further heating to supercritical conditions leads to intensive quartz precipitation and consequent permeability reduction. The initial host rock permeability and porosity are found to be main controls on the magnitude and timescales of permeability changes. While the permeability changes induced by quartz precipitation are moderate in host rocks with a primary porosity ≥ 0.05, quartz precipitation may reduce rock permeability by more than an order of magnitude in host rocks with a primary porosity of 0.025. Zones of quartz precipitation transiently change locations as the intrusion cools, thereby limiting the clogging effect, except for host rocks with low initial porosity. This permeability reduction occurs in timescales of hundreds of years in host rocks with initial high permeability and thousands of years in host rocks with intermediate permeability.

Minerals ◽  
2019 ◽  
Vol 9 (3) ◽  
pp. 195 ◽  
Author(s):  
Wenheng Liu ◽  
Xiaodong Liu ◽  
Jiayong Pan ◽  
Kaixing Wang ◽  
Gang Wang ◽  
...  

The Qingshanbao complex, part of the uranium metallogenic belt of the Longshou-Qilian mountains, is located in the center of the Longshou Mountain next to the Jiling complex that hosts a number of U deposits. However, little research has been conducted in this area. In order to investigate the origin and formation of mafic enclaves observed in the Qingshanbao body and the implications for magmatic-tectonic dynamics, we systematically studied the mineralogy, petrography, and geochemistry of these enclaves. Our results showed that the enclaves contain plagioclase enwrapped by early dark minerals. These enclaves also showed round quartz crystals and acicular apatite in association with the plagioclase. Electron probe analyses showed that the plagioclase in the host rocks (such as K-feldspar granite, adamellite, granodiorite, etc.) show normal zoning, while the plagioclase in the mafic enclaves has a discontinuous rim composition and shows instances of reverse zoning. Major elemental geochemistry revealed that the mafic enclaves belong to the calc-alkaline rocks that are rich in titanium, iron, aluminum, and depleted in silica, while the host rocks are calc-alkaline to alkaline rocks with enrichment in silica. On Harker diagrams, SiO2 contents are negatively correlated with all major oxides but K2O. Both the mafic enclaves and host rock are rich in large ion lithophile elements such as Rb and K, as well as elements such as La, Nd, and Sm, and relatively poor in high field strength elements such as Nb, Ta, P, Ti, and U. Element ratios of Nb/La, Rb/Sr, and Nb/Ta indicate that the mafic enclaves were formed by the mixing of mafic and felsic magma. In terms of rare earth elements, both the mafic enclaves and the host rock show right-inclined trends with similar weak to medium degrees of negative Eu anomaly and with no obvious Ce anomaly. Zircon LA-ICP-MS (Laser ablation inductively coupled plasma mass spectrometry) U-Pb concordant ages of the mafic enclaves and host rock were determined to be 431.8 5.2 Ma (MSWD (mean standard weighted deviation)= 1.5, n = 14) and 432.8 4.2 Ma (MSWD = 1.7, n = 16), respectively, consistent with that for the zircon U-Pb ages of the granite and medium-coarse grained K-feldspar granites of the Qingshanbao complex. The estimated ages coincide with the timing of the late Caledonian collision of the Alashan Block. This comprehensive analysis allowed us to conclude that the mafic enclaves in the Qingshanbao complex were formed by the mixing of crust-mantle magma with mantle-derived magma due to underplating, which caused partial melting of the ancient basement crust during the collisional orogenesis between the Alashan Block and Qilian rock mass in the early Silurian Period.


2014 ◽  
Vol 41 (2) ◽  
pp. 399-406 ◽  
Author(s):  
Benoit Lamy-Chappuis ◽  
Doug Angus ◽  
Quentin Fisher ◽  
Carlos Grattoni ◽  
Bruce W. D. Yardley

2009 ◽  
Vol 73 (5) ◽  
pp. 777-796 ◽  
Author(s):  
R. Macdonald ◽  
H. E. Belkin ◽  
F. Wall ◽  
B. Baginski

AbstractElectron microprobe analyses are presented of chevkinite-group minerals from Canada, USA, Guatemala, Norway, Scotland, Italy and India. The host rocks are metacarbonates, alkaline and subalkaline granitoids, quartz-bearing pegmatites, carbonatite and an inferred K-rich tuff. The analyses extend slightly the range of compositions in the chevkinite group, e.g. the most MgO-rich phases yet recorded, and we report two further examples where La is the dominant cation in the A site. Patchily- zoned crystals from Virginia and Guatemala contain both perrierite and chevkinite compositions. The new and published analyses are used to review compositional variation in minerals of the perrierite subgroup, which can form in a wide range of host rock compositions and over a substantial pressure- temperature range. The dominant substitutions in the various cation sites and a generalized substitution scheme are described.


2018 ◽  
Vol 55 (2) ◽  
pp. 130-137
Author(s):  
David E. Newton ◽  
Amy G. Ryan ◽  
Luke J. Hilchie

We use analogue experimentation to test the hypothesis that host rock competence primarily determines the morphology of kimberlite pipes. Natural occurrences of kimberlite pipes are subdivided into three classes: class 1 pipes are steep-sided diatremes emplaced into crystalline rock; class 2 pipes have a wide, shallow crater emplaced into sedimentary rock overlain by unconsolidated sediments; class 3 pipes comprise a steep-sided diatreme with a shallow-angled crater emplaced into competent crystalline rock overlain by unconsolidated sediments. We use different configurations of three analogue materials with varying cohesions to model the contrasting geological settings observed in nature. Pulses of compressed air, representing the energy of the gas-rich head of a kimberlitic magma, are used to disrupt the experimental substrate. In our experiments, the competence and configuration of the analogue materials control the excavation processes as well as the final shape of the analogue pipes: eruption through competent analogue strata results in steep-sided analogue pipes; eruption through weak analogue strata results in wide, shallow analogue pipes; eruption through intermediate strength analogue strata results in analogue pipes with a shallow crater and a steep-sided diatreme. These experimental results correspond with the shapes of natural kimberlite pipes, and demonstrate that variations in the lithology of the host rock are sufficient to generate classic kimberlite pipe shapes. These findings are consistent with models that ascribe the pipe morphologies of natural kimberlites to the competence of the host rocks in which they are emplaced.


2019 ◽  
Vol 219 (3) ◽  
pp. 1786-1804 ◽  
Author(s):  
Rémi Vachon ◽  
Christoph F Hieronymus

SUMMARY The dynamics of dyke emplacement are typically modelled by assuming an elastic rheology for the host rock. However, the resulting stress field predicts significant shear failure in the region surrounding the dyke tip. Here, we model the dyking process in an elastic-perfectly plastic host rock in order to simulate distributed shear fracturing and subsequent frictional slip on the fracture surfaces. The fluid mechanical aspects of the magma are neglected as we are interested only in the fracture mechanics of the process. Magma overpressure in dykes is typically of the same order of magnitude as the yield stress of the host rock in shear, especially when the pressure effect of volatiles exsolving from the magma is taken into account. Under these conditions, the plastic deformation zone has spatial dimensions that approach the length of the dyke itself, and concepts based on linear elastic fracture mechanics (LEFM) no longer apply. As incremental plasticity is path dependent, we describe two geologically meaningful endmember cases, namely dyke propagation at constant driving pressure, and gradual inflation of a pre-existing crack. For both models, we find that plastic deformation surrounding the fracture tip enhances dyke opening, and thus increases the energy input into the system due to pressure work integrated over the fracture wall. At the same time, energy is dissipated by plastic deformation. Dissipation in the propagation model is greater by about an order of magnitude than it is in the inflation model because the propagating dyke tip leaves behind it a broad halo of deformation due to plastic bending and unbending in the relict process zone. The net effect is that plastic deformation impedes dyke growth in the propagation model, while it enhances dyke growth in the inflation model. The results show that, when the plastic failure zone is large, a single parameter such as fracture toughness is unable to capture the physics that underpin the resistance of a fracture or dyke against propagation. In these cases, plastic failure has to be modelled explicitly for the given conditions. We provide analytical approximations for the propagation forces and the maximum dyke aperture for the two endmember cases, that is, the propagating dyke and the dyke formed by inflation of a crack. Furthermore, we show that the effect of plasticity on dyke energetics, together with an overestimate of magma pressure when interpreting dyke aspect ratios using elastic host rock models, offers a possible explanation for the long-standing paradox that laboratory measurements of fracture toughness of rocks consistently indicate values about two orders of magnitude lower than those derived from dyke observations.


2015 ◽  
Vol 79 (6) ◽  
pp. 1353-1367 ◽  
Author(s):  
J. Quirke ◽  
C. M. B. Henderson ◽  
R. A. D. Pattrick ◽  
K. M. Rosso ◽  
A. Dent ◽  
...  

AbstractGeological disposal facilities (GDF) are intended to isolate and contain radioactive waste within multiple protective barriers, deep underground, to ensure that no harmful quantities of radioactivity reach the surface environment. The last line of defense in a multi-barrier GDF is the geosphere, where iron is present in the host rock mineralogy as either Fe(II) or Fe(III), and in groundwater as Fe(II) under reducing conditions. The mobility of risk-driving radionuclides, including uranium and technetium, in the environment is affected significantly by their valence state. Due to its low redox potential, Fe(II) can mediate reduction of these radionuclides from their oxidized, highly mobile, soluble state to their reduced, insoluble state, preventing them from reaching the biosphere. Here a study of five types of potential host rocks, two granitoids, an andesite, a mudstone and a clay-rich carbonate, is reported. The bulk rocks and their minerals were analysed for iron content, Fe(II/III) ratio, and for the speciation and fine-grained nature of alteration product minerals that might have important controls on groundwater interaction. Total iron content varies between 0.9% in clays to 5.6% in the andesite. X-ray absorption spectroscopy reveals that Fe in the granitoids and andesite is predominantly Fe(II), and in mudstones, argillaceous limestone and terrestrial sandstone is predominantly Fe(III). The redox reactivity of the potential host rocks both in the presence and absence of Fe(II)-containing 'model' groundwater was investigated using an azo dye as a probe molecule. Reduction rates as determined by reactivity with the azo dye were correlated with the ability of the rocks to uptake Fe(II) from groundwater rather than with initial Fe(II) content. Potential GDF host rocks must be characterized in terms of mineralogy, texture, grain size and bulk geochemistry to assess how they might interact with groundwater. This study highlights the importance of redox reactivity, not just total iron and Fe(II)/(III) ratio, when considering the host rock performance as a barrier material to limit transport of radionuclides from the GDF.


2020 ◽  
Author(s):  
Juliette Lamarche ◽  
Nicolas Espurt ◽  
Tassadit Kaci ◽  
Marié Lionel ◽  
Richard P. Pascal

<p>Fractures in rocks are sensitive cursers that may enhance porosity and permeability. This is particularly true in carbonates because background fractures might be ubiquitous after embrittlement at early burial (Lavenu & Lamarche, 2018). Barren fractures at depth are susceptible to chemical reactions with underground fluids and cementation that might totally or partially reduce porosity and permeability (Laubach et al., 2019; Aubert et al., 2019). Hence, early background fractures with long lasting tectonic history and structural diagenesis, in addition to fractures neo-formed at any time during burial, tectonic inversion and folding join the game of matrix/fracture permeability and porosity modification. To predict the fractures contribution to flow in Naturally Fractured Reservoirs, it is fundamental to know the fracture sequence and geometry resulting from the geological history in folded carbonates, from the host-rock embrittlement to the present-day situation. At any step, we intent quantifying the fracture geometry and estimating their contribution to the host reservoir properties.</p><p> </p><p>The study is performed in Upper Jurassic to Lower Cretaceous carbonates (Oxfordian, Tithonian, Berriasian) formed in the South-Provençal Basin. From deposition to present-day, the platform carbonates underwent alternating subsidence, uplift, erosion and folding. We sampled a scan-line along a horizontal path across both flanks of the Mirabeau Anticline (SE France). We measured all tectonic and stratigraphic features crossed by the line, checked their nature and position. We deciphered their chronological relationships with respect to each other and to the bed tilting. We compiled all cross-cutting relationships into a coherent sequence of deformation of pre-, syn- and post-fold structures and correlated it to burial, uplift and folding of the host rock. At each brittle stage, the fracture pattern was characterized in terms of architecture, mechanical stratigraphy and reservoir properties in order to draw a time-path in a matrix versus fracture permeability and porosity table (Nelson Reservoir types) during 150My. After embrittlement, the host-rocks bear fractures, pressure-solution, faulting, folding and erosion. If it was a reservoir, its Nelson type would have evolved from IV to III during the burial and initial brittle deformation. The tectonic inversion and onset of multiple-scale brittle structures would have increased and decreased the fracture and matrix contribution respectively and the reservoir evolved to types II and I. During the 150My history, fracture porosity and permeability depends on their geometry (veins versus tension gashes) and cementation. This results in several switches from type II to I as a function of the fracture timing, geometry, connectivity and diagenesis.</p><p> </p><p>Aubert I. et al. (2019). Imbricated structure and hydraulic path induced by strike-slip reactivation of a normal fault in carbonates. Fifth International Conference on Fault and Top Seals, 8-12 September 2019, Palermo, Italy.</p><p>Bestani L.et al. (2016) Reconstruction of the Provence Chain evolution, southeastern France., Tectonics Vol: 35, p.1506–1525</p><p>Laubach, S. E. et al. (2019) The role of chemistry in fracture pattern development and opportunities to advance interpretations of geological materials. Reviews Geophysics, 57.</p><p>Lavenu A.P.C., Lamarche J. (2018) What controls diffuse fractures in platform carbonates? Insights from Provence (France) and Apulia (Italy), JSG 108, p. 94-107</p>


2020 ◽  
Author(s):  
Dániel Kiss ◽  
Evangelos Moulas ◽  
Lisa Rummel ◽  
Boris Kaus

<p>A recent focus of studies in geodynamic modeling and magmatic petrology is to understand the coupled behavior between deformation and magmatic processes. Here, we present a 2D numerical model of an upper crustal magma (or mush) chamber in a visco-elastic host rock, with coupled thermal, mechanical and chemical (TMC) processes. The magma chamber is isolated from deeper sources of magma and it is cooling, and thus shrinking. We quantify the mechanical interaction between the shrinking magma chamber and the surrounding host rock, using a compressible visco-elastic formulation, considering several geometries of the magma chamber.</p><p>We present a self-consistent system of the conservation equations for coupled TMC processes, under the assumptions of slow (negligible inertial forces), visco-elastic deformation and constant chemical bulk composition. The thermodynamic melting/crystallization model is based on a pelitic melting model calculated with Perple_X, assuming a granitic composition and is incorporated as a look-up table. We will discuss the numerical implementation, show the results of systematic numerical simulations, and illustrate the effect of volume changes due to crystallization on stresses in the host rocks.</p>


Sign in / Sign up

Export Citation Format

Share Document