The inclusion record of fluid evolution,crack healing and trapping from a heterogeneous system during rapid cooling of pegmatitic veins (Dronning Maud Land; Antarctica)

Granitoid (pegmatite and aplite) veins in metamorphic rocks and intrusive syenites of central Dronning Maud Land, Antarctica, are flanked by conspicuous light‐coloured alteration halos, which represent the damage zone of fracture propagation. The damage zone is characterized by a high density of sealed or healed microcracks, about 1 order of magnitude above background. Fluid inclusions along healed microcracks in quartz of both pegmatite and alteration halos are inspected by optical and scanning electron microscopy, and their composition is analysed by microthermometry and quadrupole mass spectrometry. The similar inclusion record in the granitoid vein and in the damaged host rock indicates the derivation of the fluids from the hydrous melt phase. The aqueous inclusions bear abundant daughter crystals, mainly silicates, and may represent a hydrous melt. The volatile composition is variable in the system H₂O–CO₂, with mostly subordinate amounts of N₂. Phase separation with partitioning of CO₂ into the fluid phase coexisting with the hydrous melt, and possibly immiscibility in the subsolidus range, govern fluid evolution during cooling. The variable CO₂/N₂ ratio suggests mixing with fluids from an external source in the host rock and vigorous circulation at an early stage of high transient permeability. Experiments have shown that healing of microcracks at high temperatures is a matter of hours to weeks, hence similar in time scale to the cooling of the cm‐ to dm‐thick granitoid veins. In this case, rapid cooling and concomitant crack healing in a system undergoing phase separation causes a broad compositional variability of the inclusions due to necking down, and the underpressure developing in closed compartments precludes a meaningful thermobarometric interpretation.     

Hydrothermal dolomitization,mixing corrosion and deep burial porosity formation: numerical results from 1‐D reactive transport models

Major corrosion has been found at depth in carbonate hydrocarbon reservoirs from different geologic provinces. Fluid inclusion microthermometry and stable isotopic compositions of carbonate cements, predating major corrosion, constrain the interpretation of the evolution of parental fluids during progressive burial and prior to the major corrosion event. Post‐major corrosion mineral paragenesis includes pyrite (‐marcasite), anhydrite, kaolinite (dickite) and fluorite. Although the post‐corrosion mineral paragenesis represents minor volumes of rock, it may provide valuable insights into the post‐corrosion brine chemistry. Using reactive transport numerical models, the roles of cooling and/or mixing of brines on corrosion have been evaluated as controls for dolomitization, deep burial corrosion and precipitation of the post‐corrosion mineral paragenesis. Modelling results show that cooling of deep‐seated fluids moving upward along a fracture may cause minor calcite dissolution and porosity generation. Significant dolomitization along a fracture zone and nearby host‐rock only occurs when deep‐seated fluids have high salinities (4 mol Cl kg^?1 of solution) and Darcian flow rates are relatively high (1 m³ m^?2 year^?1). Only minor volumes of quartz and fluorite precipitate in the newly formed porosity. Moreover, modelling results cannot reproduce the authigenic precipitation of kaolinite (dickite at high temperatures) by cooling. As an alternative to cooling as a cause of corrosion, mixing between two brines of different compositions and salinities is represented by two main cases. One case consists of the flow up along a fracture of deep‐seated fluids with higher salinities than the fluid in the wall rock. Dolomite does not precipitate at a fracture zone. Nevertheless, minor volumes of dolomite are formed away from the fracture. The post‐corrosion mineral paragenesis can be partly reproduced, and the results are comparable to those obtained from cooling calculations. Minor volumes of quartz and fluorite are formed, and kaolinite‐dickite does not precipitate. The major outputs of this scenario are calcite dissolution and slight net increase in porosity. A second case corresponds to the mixing of low salinity deep‐seated fluids, flowing up along fractures, with high salinity brines within the wall rock. Calculations predict major dissolution of calcite and precipitation of dolomite. The post‐corrosion mineral paragenesis can be reproduced. High volumes of quartz, fluorite and kaolinite‐dickite precipitate and may even completely occlude newly formed porosity.     

Hybrid finite element–finite volume discretization of complex geologic structures and a new simulation workflow demonstrated on fractured rocks

The generation of computational meshes of complex geological objects is a challenge: their shape needs to be retained, resolution has to adapt to local detail, and variations of material properties in the objects have to be represented. Also mesh refinement and adaptation must be sufficient to resolve variations in the computed variable(s). Here, we present an unstructured hybrid finite element, node‐centred finite‐volume discretization suitable for solving fluid flow, reactive transport, and mechanical partial differential equations on a complex geometry with inhomogeneous material domains. We show that resulting meshes accurately capture free‐form material interfaces as defined by non‐uniform rational B‐spline curves and surfaces. The mesh discretization error is analysed for the elliptic pressure equation and an error metric is introduced to guide mesh refinement. Finite elements and finite volumes are represented in parametric space and integrations are conducted numerically. Subsequently, integral properties are mapped to physical space using Jacobian transformations. This method even retains its validity when the mesh is deformed. The resulting generic formulation is demonstrated for a transport calculation performed on a complex discrete fracture model.     

Role of compressive tectonics on gas charging into the Ordovician sandstone reservoirs in the Sbaa Basin,Algeria: constrained by fluid inclusions and mineralogical data

Structure‐ and tectonic‐related gas migration into Ordovician sandstone reservoirs and its impact on diagenesis history were reconstructed in two gas fields in the Sbaa Basin, in SW Algeria. This was accomplished by petrographical observations, fluid inclusion microthermometry and stable isotope geochemistry on quartz, dickite and carbonate cements and veins. Two successive phases of quartz cementation (CQ1 and CQ2) occurred in the reservoirs. Two phase aqueous inclusions show an increase in temperatures and salinities from the first CQ1 diagenetic phase toward CQ2 in both fields. Microthermometric data on gas inclusions in quartz veins reveal the presence of an average of 92 ± 5 mole% of CH₄ considering a CH₄‐CO₂ system, which is similar to the present‐day gas composition in the reservoirs. The presence of primary methane inclusions in early quartz overgrowths and in quartz and calcite veins suggests that hydrocarbon migration into the reservoir occurred synchronically with early quartz cementation in the sandstones located near the contact with the Silurian gas source rock at 100–140°C during the Late Carboniferous period and the late Hercynian episode fracturing at temperatures between 117 and 185°C, which increased in the NW‐direction of the basin. During the fracture filling, three main types of fluids were identified with different salinities and formation temperatures. A supplementary phase of higher fluid temperature (up to 226°C) recorded in late quartz, and calcite veins is related to a Jurassic thermal event. The occurrence of dickite cements close to the Silurian base near the main fault areas in both fields is mainly correlated with the sandstones where the early gas was charged. It implies that dickite precipitation is related to acidic influx. Late carbonate cements and veins (calcite – siderite – ankerite and strontianite) occurred at the same depths resulting from the same groundwater precipitation. The absence of methane inclusions in calcite cements result from methane flushing by saline waters.     

Carbon dioxide controlled earthquake distribution pattern in the NW Bohemian swarm earthquake region,western Eger Rift,Czech Republic – gas migration in the crystalline basement

The ascent of magmatic carbon dioxide in the western Eger (Oh?e) Rift is interlinked with the fault systems of the Variscian basement. In the Cheb Basin, the minimum CO₂ flux is about 160 m³ h^?1, with a diminishing trend towards the north and ceasing in the main epicentral area of the Northwest Bohemian swarm earthquakes. The ascending CO₂ forms Ca‐Mg‐HCO₃ type waters by leaching of cations from the fault planes and creates clay minerals, such as kaolinite, as alteration products on affected fault planes. These mineral reactions result in fault weakness and in hydraulically interconnected fault network. This leads to a decrease in the friction coefficient of the Coulomb failure stress (CFS) and to fault creep as stress build‐up cannot occur in the weak segments. At the transition zone in the north of the Cheb Basin, between areas of weak, fluid conductive faults and areas of locked faults with frictional strength, fluid pressure can increase resulting in stress build‐up. This can trigger strike‐slip swarm earthquakes. Fault creep or movements in weak segments may support a stress build‐up in the transition area by transmitting fluid pressure pulses. Additionally to fluid‐driven triggering models, it is important to consider that fluids ascending along faults are CO₂‐supersaturated thus intensifying the effect of fluid flow. The enforced flow of CO₂‐supersaturated fluids in the transitional zone from high to low permeability segments through narrowings triggers gas exsolution and may generate pressure fluctuations. Phase separation starts according to the phase behaviour of CO₂‐H₂O systems in the seismically active depths of NW Bohemia and may explain the vertical distribution of the seismicity. Changes in the size of the fluid transport channels in the fault systems caused, or superimposed, by fault movements, can produce fluid pressure increases or pulses, which are the precondition for triggering fluid‐induced swarm earthquakes.     

Fluid effect on hydraulic fracture propagation behavior: a comparison between water and supercritical CO2‐like fluid

The initiation of hydraulic fractures during fluid injection in deep formations can be either engineered or induced unintentionally. Upon injection of CO₂, the pore fluids in deep formations can be changed from oil/saline water to CO₂ or CO₂ dominated. The type of fluid is important not only because the fluid must fracture the rock, but also because rocks saturated with different pore fluids behave differently. We investigated the influence of fluid properties on fracture propagation behavior by using the cohesive zone model in conjunction with a poroelasticity model. Simulation results indicate that the pore pressure fields are very different for different pore fluids even when the initial field conditions and injection schemes (rate and time) are kept the same. Low viscosity fluids with properties of supercritical CO₂ will create relatively thin and much shorter fractures in comparison with fluids exhibiting properties of water under similar injection schemes. Two significant times are recognized during fracture propagation: the time at which a crack ceases opening and the later time point at which a crack ceases propagating. These times are very different for different fluids. Both fluid compressibility and viscosity influence fracture propagation, with viscosity being the more important property. Viscosity can greatly affect hydraulic conductivity and the leak‐off coefficient. This analysis assumes the in‐situ pore fluid and injected fluid are the same and the pore space is 100% saturated by that fluid at the beginning of the simulation.     

Precipitation of lead–zinc ores in the Mississippi Valley‐type deposit at Trèves,Cévennes region of southern France

The Trèves zinc–lead deposit is one of several Mississippi Valley‐type (MVT) deposits in the Cévennes region of southern France. Fluid inclusion studies show that the ore was deposited at temperatures between approximately 80 and 150°C from a brine that derived its salinity mainly from the evaporation of seawater past halite saturation. Lead isotope studies suggest that the metals were extracted from local basement rocks. Sulfur isotope data and studies of organic matter indicate that the reduced sulfur in the ores was derived from the reduction of Mesozoic marine sulfate by thermochemical sulfate reduction or bacterially mediated processes at a different time or place from ore deposition. The large range of δ³⁴S values determined for the minerals in the deposit (12.2–19.2‰ for barite, 3.8–13.8‰ for sphalerite and galena, and 8.7 to ?21.2‰ for pyrite), are best explained by the mixing of fluids containing different sources of sulfur. Geochemical reaction path calculations, based on quantitative fluid inclusion data and constrained by field observations, were used to evaluate possible precipitation mechanisms. The most important precipitation mechanism was probably the mixing of fluids containing different metal and reduced sulfur contents. Cooling, dilution, and changes in pH of the ore fluid probably played a minor role in the precipitation of ores. The optimum results that produced the most metal sulfide deposition with the least amount of fluid was the mixing of a fluid containing low amounts of reduced sulfur with a sulfur‐rich, metal poor fluid. In this scenario, large amounts of sphalerite and galena are precipitated, together with smaller quantities of pyrite precipitated and dolomite dissolved. The relative amounts of metal precipitated and dolomite dissolved in this scenario agree with field observations that show only minor dolomite dissolution during ore deposition. The modeling results demonstrate the important control of the reduced sulfur concentration on the Zn and Pb transport capacity of the ore fluid and the volumes of fluid required to form the deposit. The studies of the Trèves ores provide insights into the ore‐forming processes of a typical MVT deposit in the Cévennes region. However, the extent to which these processes can be extrapolated to other MVT deposits in the Cévennes region is problematic. Nevertheless, the evidence for the extensive migration of fluids in the basement and sedimentary cover rocks in the Cévennes region suggests that the ore forming processes for the Trèves deposit must be considered equally viable possibilities for the numerous fault‐controlled and mineralogically similar MVT deposits in the Cévennes region.     

Fluid flow in shear zones: insights from the geometry and evolution of ore bodies at Renco gold mine,Zimbabwe

The geometry of mineral deposits can give insights into fluid flow in shear zones. Lode gold ore bodies at Renco Mine, in the Limpopo Belt, Zimbabwe, occur as siliceous breccias and mylonites within amphibolite facies shear zones that dip either gently or steeply. The two sets of ore bodies formed synchronously from hydrothermal fluids. The ore bodies are oblate, but have well‐defined long axes. Larger ore bodies are more oblate. High‐grade gold ore shoots have long axes that plunge down dip; this direction is perpendicular to the long axes of the low‐grade ore bodies. The centres of the high‐grade ore bodies align within the low‐grade ore bodies along strike in both gently and steeply dipping groups. The range of sizes and shapes of the ore bodies are interpreted as a growth sequence. Geometrical models are proposed for the gently and steeply dipping ore bodies, in which individual ore bodies grow with long axes plunging down dip, and merge to form larger, more oblate ore bodies. The models show that when three or more ore bodies coalesce, the long axis of the merged ore body is perpendicular to the component ore bodies, and that ore bodies in the deposit may have a range of shapes due to both growth of individual ore bodies, and their coalescence. The long axes of the high‐grade ore bodies are parallel to the shear directions of both the gently and steeply dipping dip slip shear zones, which were the directions of greatest permeability and fluid flow. The larger, lower grade bodies, which may have formed by coalescence, are elongate perpendicular to these directions.     

Thermal convection in faulted extensional sedimentary basins: theoretical results from finite-element modeling

Thermal convection has the potential to be a significant and widespread mechanism of fluid flow, mass transport, and heat transport in rift and other extensional basins. Based on numerical simulation results, large‐scale convection can occur on the scale of the basin thickness, depending on the Rayleigh number for the basin. Our analysis indicates that for syn‐rift and early post‐rift settings with a basin thickness of 5 km, thermal convection can occur for basal heat flows ranging from 80 to 150 mW m^?2, when the vertical hydraulic conductivity is on the order of 1.5 m year^?1 and lower. The convection cells have characteristic wavelengths and flow patterns depending on the thermal and hydraulic boundary conditions. Steeply dipping extensional faults can provide pathways for vertical fluid flow across large thicknesses of basin sediments and can modify the dynamics of thermal convection. The presence of faults perturbs the thermal convective flow pattern and can constrain the size and locations of convection cells. Depending on the spacing of the faults and the hydraulic properties of the faults and basin sediments, the convection cells can be spatially organized to align with adjacent faults. A fault‐bounded cell occurs when one convection cell is constrained to occupy a fault block so that the up‐flow zone converges into one fault zone and the down‐flow zone is centred on the adjacent fault. A fault‐bounded cell pair occurs when two convection cells occupy a fault block with the up‐flow zone located between the faults and the down‐flow zones centred on the adjacent faults or with the reverse pattern of flow. Fault‐bounded cells and cell pairs can be referred to collectively as fault‐bounded convective flow. The flow paths in fault‐bounded convective flow can be lengthened significantly with respect to those of convection cells unperturbed by the presence of faults. The cell pattern and sense of circulation depend on the fault spacing, sediment and fault permeabilities, lithologic heterogeneity, and the basal heat flow. The presence of fault zones also extends the range of conditions for which thermal convection can occur to basin settings with Rayleigh numbers below the critical value for large‐scale convection to occur in a basin without faults. The widespread potential for the occurrence of thermal convection suggests that it may play a role in controlling geological processes in rift basins including the acquisition and deposition of metals by basin fluids, the distribution of diagenetic processes, the temperature field and heat flow, petroleum generation and migration, and the geochemical evolution of basin fluids. Fault‐bounded cells and cell pairs can focus mass and heat transport from longer flow paths into fault zones, and their discharge zones are a particularly favourable setting for the formation of sediment‐hosted ore deposits near the sea floor.     

The application of structured‐light illumination microscopy to hydrocarbon‐bearing fluid inclusions

Structured‐light illumination (SLI)‐based microscopy offers geologists a new perspective for screening of hydrocarbon‐bearing (HCFI) and small aqueous fluid inclusion (AFI) assemblages. This optical‐sectioning technique provides rapid, confocal‐like imaging, using relatively simple and inexpensive instrumentation. The 3D fluorescent images of HCFI planes and large single HCFIs permit the visualization of the relationships between HCFI assemblages, examination of HCFI morphology, and volume estimates of the fluorescent components within HCFIs. By the use of normal white light illumination, SLI image capture, and varying acquisition time it is also possible to image AFI because of the random movements of vapour bubbles within the inclusions. This allows the near‐simultaneous visualization of hydrocarbon and AFI which is of significant importance for the study of sedimentary basins and petroleum reservoirs. SLI is a unique and accessible 3D petrographic tool, with practical advantages over conventional epifluorescence and confocal laser scanning microscopy.