Refine
Has Fulltext
- no (11)
Document Type
- Article (11) (remove)
Language
- English (11)
Is part of the Bibliography
- yes (11)
Keywords
- PHREEQC (2)
- digital rock physics (2)
- elastic properties (2)
- micro-CT (2)
- numerical simulation (2)
- Bentheim sandstone (1)
- Cantera (1)
- Code_Aster (1)
- Mallik (1)
- Pitzer (1)
Institute
Geochemical processes change the microstructure of rocks and thereby affect their physical behaviour at the macro scale. A micro-computer tomography (micro-CT) scan of a typical reservoir sandstone is used to numerically examine the impact of three spatial alteration patterns on pore morphology, permeability and elastic moduli by correlating precipitation with the local flow velocity magnitude. The results demonstrate that the location of mineral growth strongly affects the permeability decrease with variations by up to four orders in magnitude. Precipitation in regions of high flow velocities is characterised by a predominant clogging of pore throats and a drastic permeability reduction, which can be roughly described by the power law relation with an exponent of 20. A continuous alteration of the pore structure by uniform mineral growth reduces the permeability comparable to the power law with an exponent of four or the Kozeny-Carman relation. Preferential precipitation in regions of low flow velocities predominantly affects smaller throats and pores with a minor impact on the flow regime, where the permeability decrease is considerably below that calculated by the power law with an exponent of two. Despite their complete distinctive impact on hydraulics, the spatial precipitation patterns only slightly affect the increase in elastic rock properties with differences by up to 6.3% between the investigated scenarios. Hence, an adequate characterisation of the spatial precipitation pattern is crucial to quantify changes in hydraulic rock properties, whereas the present study shows that its impact on elastic rock parameters is limited. The calculated relations between porosity and permeability, as well as elastic moduli can be applied for upscaling micro-scale findings to reservoir-scale models to improve their predictive capabilities, what is of paramount importance for a sustainable utilisation of the geological subsurface.
Quantifying rock weakening due to decreasing calcite mineral content by numerical simulations
(2018)
The quantification of changes in geomechanical properties due to chemical reactions is of paramount importance for geological subsurface utilisation, since mineral dissolution generally reduces rock stiffness. In the present study, the effective elastic moduli of two digital rock samples, the Fontainebleau and Bentheim sandstones, are numerically determined based on micro-CT images. Reduction in rock stiffness due to the dissolution of 10% calcite cement by volume out of the pore network is quantified for three synthetic spatial calcite distributions (coating, partial filling and random) using representative sub-cubes derived from the digital rock samples. Due to the reduced calcite content, bulk and shear moduli decrease by 34% and 38% in maximum, respectively. Total porosity is clearly the dominant parameter, while spatial calcite distribution has a minor impact, except for a randomly chosen cement distribution within the pore network. Moreover, applying an initial stiffness reduced by 47% for the calcite cement results only in a slightly weaker mechanical behaviour. Using the quantitative approach introduced here substantially improves the accuracy of predictions in elastic rock properties compared to general analytical methods, and further enables quantification of uncertainties related to spatial variations in porosity and mineral distribution.
Diagenetic trends of synthetic reservoir sandstone properties assessed by digital rock physics
(2021)
Quantifying interactions and dependencies among geometric, hydraulic and mechanical properties of reservoir sandstones is of particular importance for the exploration and utilisation of the geological subsurface and can be assessed by synthetic sandstones comprising the microstructural complexity of natural rocks. In the present study, three highly resolved samples of the Fontainebleau, Berea and Bentheim sandstones are generated by means of a process-based approach, which combines the gravity-driven deposition of irregularly shaped grains and their diagenetic cementation by three different schemes. The resulting evolution in porosity, permeability and rock stiffness is examined and compared to the respective micro-computer tomographic (micro-CT) scans. The grain contact-preferential scheme implies a progressive clogging of small throats and consequently produces considerably less connected and stiffer samples than the two other schemes. By contrast, uniform quartz overgrowth continuously alters the pore space and leads to the lowest elastic properties. The proposed stress-dependent cementation scheme combines both approaches of contact-cement and quartz overgrowth, resulting in granulometric, hydraulic and elastic properties equivalent to those of the respective micro-CT scans, where bulk moduli slightly deviate by 0.8%, 4.9% and 2.5% for the Fontainebleau, Berea and Bentheim sandstone, respectively. The synthetic samples can be further altered to examine the impact of mineral dissolution or precipitation as well as fracturing on various petrophysical correlations, which is of particular relevance for numerous aspects of a sustainable subsurface utilisation.
Injection of fluids into deep saline aquifers causes a pore pressure increase in the storage formation, and thus displacement of resident brine. Via hydraulically conductive faults, brine may migrate upwards into shallower aquifers and lead to unwanted salinisation of potable groundwater resources. In the present study, we investigated different scenarios for a potential storage site in the Northeast German Basin using a three-dimensional (3-D) regional-scale model that includes four major fault zones. The focus was on assessing the impact of fault length and the effect of a secondary reservoir above the storage formation, as well as model boundary conditions and initial salinity distribution on the potential salinisation of shallow groundwater resources. We employed numerical simulations of brine injection as a representative fluid. Our simulation results demonstrate that the lateral model boundary settings and the effective fault damage zone volume have the greatest influence on pressure build-up and development within the reservoir, and thus intensity and duration of fluid flow through the faults. Higher vertical pressure gradients for short fault segments or a small effective fault damage zone volume result in the highest salinisation potential due to a larger vertical fault height affected by fluid displacement. Consequently, it has a strong impact on the degree of shallow aquifer salinisation, whether a gradient in salinity exists or the saltwater-freshwater interface lies below the fluid displacement depth in the faults. A small effective fault damage zone volume or low fault permeability further extend the duration of fluid flow, which can persist for several tens to hundreds of years, if the reservoir is laterally confined. Laterally open reservoir boundaries, large effective fault damage zone volumes and intermediate reservoirs significantly reduce vertical brine migration and the potential of freshwater salinisation because the origin depth of displaced brine is located only a few decametres below the shallow aquifer in maximum. The present study demonstrates that the existence of hydraulically conductive faults is not necessarily an exclusion criterion for potential injection sites, because salinisation of shallower aquifers strongly depends on initial salinity distribution, location of hydraulically conductive faults and their effective damage zone volumes as well as geological boundary conditions.
Leaching zones within potash seams generally represent a significant risk to subsurface mining operations and the construction of technical caverns in salt rocks, but their temporal and spatial formation has been investigated only rudimentarily to date. To the knowledge of the authors, current reactive transport simulation implementations are not capable to address hydraulic-chemical interactions within potash salt. For this reason, a reactive transport model has been developed and complemented by an innovative approach to calculate the interchange of minerals and solution at the water-rock interface. Using this model, a scenario analysis was carried out based on a carnallite-bearing potash seam. The results show that the evolution of leaching zones depends on the mineral composition and dissolution rate of the original salt rock, and that the formation can be classified by the dimensionless parameters of Peclet (Pe) and Damkohler (Da). For Pe > 2 and Da > 1, a funnel-shaped leaching zone is formed, otherwise the dissolution front is planar. Additionally, Da > 1 results in the formation of a sylvinitic zone and a flow barrier. Most scenarios represent hybrid forms of these cases. The simulated shapes and mineralogies are confirmed by literature data and can be used to assess the hazard potential.
How insoluble inclusions and intersecting layers affect the leaching process within potash seams
(2021)
Potash seams are a valuable resource containing several economically interesting, but also highly soluble minerals. In the presence of water, uncontrolled leaching can occur, endangering subsurface mining operations. In the present study, the influence of insoluble inclusions and intersecting layers on leaching zone evolution was examined by means of a reactive transport model. For that purpose, a scenario analysis was carried out, considering different rock distributions within a carnallite-bearing potash seam. The results show that reaction-dominated systems are not affected by heterogeneities at all, whereas transport-dominated systems exhibit a faster advance in homogeneous rock compositions. In return, the ratio of permeated rock in vertical direction is higher in heterogeneous systems. Literature data indicate that most natural potash systems are transport-dominated. Accordingly, insoluble inclusions and intersecting layers can usually be seen as beneficial with regard to reducing hazard potential as long as the mechanical stability of leaching zones is maintained. Thereby, the distribution of insoluble areas is of minor impact unless an inclined, intersecting layer occurs that accelerates leaching zone growth in one direction. Moreover, it is found that the saturation dependency of dissolution rates increases the growth rate in the long term, and therefore must be considered in risk assessments.
Underground coal gasification (UCG) is an in situ conversion technique that enables the production of high-calorific synthesis gas from resources that are economically not minable by conventional methods. A broad range of end-use options is available for the synthesis gas, including fuels and chemical feedstock production. Furthermore, UCG also offers a high potential for integration with Carbon Capture and Storage (CCS) to mitigate greenhouse gas emissions. In the present study, a stoichiometric equilibrium model, based on minimization of the Gibbs function has been used to estimate the equilibrium composition of the synthesis gas. Thereto, we further developed and applied a proven thermodynamic equilibrium model to simulate the relevant thermochemical coal conversion processes (pyrolysis and gasification). Our modeling approach has been validated against thermodynamic models, laboratory gasification experiments and UCG field trial data reported in the literature. The synthesis gas compositions have been found to be in good agreement under a wide range of different operating conditions. Consequently, the presented modeling approach enables an efficient quantification of synthesis gas quality resulting from UCG, considering varying coal and oxidizer compositions at deposit-specific pressures and temperatures.
Underground coal gasification (UCG) enables utilization of coal reserves, currently not economically exploitable due to complex geological boundary conditions. Hereby, UCG produces a high-calorific synthesis gas that can be used for generation of electricity, fuels, and chemical feedstock. The present study aims to identify economically-competitive, site-specific end-use options for onshore- and offshore-produced UCG synthesis gas, taking into account the capture and storage (CCS) and/or utilization (CCU) of produced CO2. Modeling results show that boundary conditions favoring electricity, methanol, and ammonia production expose low costs for air separation, low compression power requirements, and appropriate shares of H-2/N-2. Hereby, a gasification agent ratio of more than 30% oxygen by volume is not favorable from the economic and CO2 mitigation viewpoints. Compared to the costs of an offshore platform with its technical equipment, offshore drilling costs are marginal. Thus, uncertainties related to parameters influenced by drilling costs are negligible. In summary, techno-economic process modeling results reveal that air-blown gasification scenarios are the most cost-effective ones, while offshore UCG-CCS/CCU scenarios are up to 1.7 times more expensive than the related onshore processes. Hereby, all investigated onshore scenarios except from ammonia production under the assumed worst-case conditions are competitive on the European market.
Power-to-gas (PtG) stores chemical energy by converting excess electrical energy from renewable sources into an energy-dense gas. Due to its higher available capacity compared to surface-based storage technologies, subsurface storage in geological systems is the most promising approach for efficient and economic realization of the PtG system’s storage component. For this purpose, methane (CH4) produced by methanation by means of hydrogen (H2) and carbon dioxide (CO2) is stored in a geological reservoir until required for further use. In this context, CO2 is used as the cushion gas to maintain reservoir pressure and limiting working gas, i.e., (CH4) losses during withdrawal periods. Consequently, mixing of both gases in the reservoir is inevitable. Therefore, it is necessary to minimize the gas mixing region to optimize the efficiency of the PtG system’s storage component. In the present study, the physical properties of CH4, CO2 and their mixtures are reviewed. Then, a multicomponent flow model is implemented and validated against published data. Next, a hydromechanically coupled model is established, considering fluid flow through porous media and effective stresses to investigate the mixing behavior of both gases and the mechanical reservoir stability. The simulation results show that, with increasing reservoir thickness and dip angle, the mixing region is reduced during gas injection if CO2 is employed as the cushion gas. In addition, the degree of mixing is lower at higher temperatures. Feasible injection rates and injection schedules can be derived from the integrated reservoir stability analysis. The methodology developed in the present study allows the determination of optimum strategies for storage reservoir selection and gas injection scheduling by minimizing the gas mixing region.
The LArge-scale Reservoir Simulator (LARS) has been previously developed to study hydrate dissociation in hydrate-bearing systems under in-situ conditions. In the present study, a numerical framework of equations of state describing hydrate formation at equilibrium conditions has been elaborated and integrated with a numerical flow and transport simulator to investigate a multi-stage hydrate formation experiment undertaken in LARS. A verification of the implemented modeling framework has been carried out by benchmarking it against another established numerical code. Three-dimensional (3D) model calibration has been performed based on laboratory data available from temperature sensors, fluid sampling, and electrical resistivity tomography. The simulation results demonstrate that temperature profiles, spatial hydrate distribution, and bulk hydrate saturation are consistent with the observations. Furthermore, our numerical framework can be applied to calibrate geophysical measurements, optimize post-processing workflows for monitoring data, improve the design of hydrate formation experiments, and investigate the temporal evolution of sub-permafrost methane hydrate reservoirs.