A.M.H. Pluymakers
Please Note
46 records found
1
In CCS wells, cyclic injection of cold CO2 into the hot subsurface may lead to debonding between sealant and steel casing. We test how thermal cycling affects the sealing ability of five different types of sealant (S1 to S5) surrounding a simulated steel wellbore. We use cylindrical sealant samples with a stainless steel (AISI 316 L) pipe in the centre, cured at 150°C and 30 MPa for 28 days. Using 3 bar N2 leak tests at room temperature, we test how much the sealant-steel interface leaks before and after thermal cycling under unconfined and confined (1.5 MPa) conditions. We also conduct push-off experiments using a 500 kN loading frame before and after. For the unconfined test, we place the sample on a custom-built jig, whereas for confined tests we have a similar assembly inside a conventional triaxial vessel. The samples are brought to 60°C. Subsequently, we inject 5°C water through the central pipe at 80 mL/min for 2 mins, and let the sample reheat for 12 mins. We repeat this 16 times. Afterwards, we allow the sample to cool to room temperature, and repeat the N2 leak test in-situ. The results show that under unconfined conditions, the interface leaks more for all sealant types except S3. The key parameter controlling performance is the linear thermal expansion coefficient, where an expansion coefficient closer to that of steel indicates better performance. Under confinement, all sealant types perform better post-thermal cycling, due to the prolonged exposure to confining pressure.
Coupling Neutron and X‐Ray Imaging of Fluid Mixing and Precipitation in Rocks
Challenges and Opportunities
Solute mixing in rocks plays a central role in a wide range of reactive processes. However, how the complex 3D pore structure of rocks governs mixing rates remains largely unknown. Moreover, some mixing-driven reactions—such as dissolution and precipitation—can modify the pore space, with poorly understood consequences for mixing itself. Recent advances in X-ray imaging techniques have significantly enhanced our ability to visualize the pore-scale rock architecture of rocks and a wide range of fluid processes. However, capturing solute mixing and its impact on chemical reactions—such as mineralization—remains a major challenge. Here, we investigated the potential of coupling time-lapse 3D neutron and X-ray imaging to characterize reactive fluid mixing and subsequent calcium carbonate mineralization in porous basalt. Two flow-through experiments were performed with co-injected CaCl2 and Na2CO3, leading to precipitation. Neutron imaging tracked fluid mixing, while X-ray imaging distinguished the solid matrix from pore space for fluid analysis. The first experiment showed steady transverse mixing, while a second experiment revealed temporal fluctuations due to trapped air, causing multiphase flow. Neutron images indicated significant fluid mixing driven by these fluctuations. A synchrotron X-ray image post-experiment indicated additional mineral precipitates from long-term diffusive mixing. Despite the promising results, several challenges remain, including resolution limits, temporal synchronization between modalities, and accurate fluid phase segmentation. Overall, our findings highlight both the potential and limitations of integrated neutron and X-ray imaging for studying pore-scale reactive transport and mineralization processes.
Plain Language Summary
Understanding fluid movement inside rocks is crucial for enhancing CO2 storage and other underground applications. This study used neutron and X-ray imaging to explore how reactive fluids mix in basalt rocks and how this process forms calcium carbonate, aiding long-term carbon storage. We conducted experiments by pumping fluids through basalt samples at varying flow rates and employed advanced imaging to observe the mixing and subsequent calcite precipitation. Neutron imaging tracked fluid movement, while X-ray imaging revealed the pores filled with fluids and calcite. Our findings underscore the challenges and opportunities of using these imaging techniques to study pore-scale mixing and reactive transport, providing a framework for future research on optimizing mineral precipitation processes.
...
Solute mixing in rocks plays a central role in a wide range of reactive processes. However, how the complex 3D pore structure of rocks governs mixing rates remains largely unknown. Moreover, some mixing-driven reactions—such as dissolution and precipitation—can modify the pore space, with poorly understood consequences for mixing itself. Recent advances in X-ray imaging techniques have significantly enhanced our ability to visualize the pore-scale rock architecture of rocks and a wide range of fluid processes. However, capturing solute mixing and its impact on chemical reactions—such as mineralization—remains a major challenge. Here, we investigated the potential of coupling time-lapse 3D neutron and X-ray imaging to characterize reactive fluid mixing and subsequent calcium carbonate mineralization in porous basalt. Two flow-through experiments were performed with co-injected CaCl2 and Na2CO3, leading to precipitation. Neutron imaging tracked fluid mixing, while X-ray imaging distinguished the solid matrix from pore space for fluid analysis. The first experiment showed steady transverse mixing, while a second experiment revealed temporal fluctuations due to trapped air, causing multiphase flow. Neutron images indicated significant fluid mixing driven by these fluctuations. A synchrotron X-ray image post-experiment indicated additional mineral precipitates from long-term diffusive mixing. Despite the promising results, several challenges remain, including resolution limits, temporal synchronization between modalities, and accurate fluid phase segmentation. Overall, our findings highlight both the potential and limitations of integrated neutron and X-ray imaging for studying pore-scale reactive transport and mineralization processes.
Plain Language Summary
Understanding fluid movement inside rocks is crucial for enhancing CO2 storage and other underground applications. This study used neutron and X-ray imaging to explore how reactive fluids mix in basalt rocks and how this process forms calcium carbonate, aiding long-term carbon storage. We conducted experiments by pumping fluids through basalt samples at varying flow rates and employed advanced imaging to observe the mixing and subsequent calcite precipitation. Neutron imaging tracked fluid movement, while X-ray imaging revealed the pores filled with fluids and calcite. Our findings underscore the challenges and opportunities of using these imaging techniques to study pore-scale mixing and reactive transport, providing a framework for future research on optimizing mineral precipitation processes.
This study explores two approaches to assess stress changes: a semi-analytical geomechanical proxy and a fully-coupled Thermo-Hydro-Mechanical (THM) model using open-DARTS. The THM model simulates coupled thermal, hydraulic, and mechanical processes in complex rock formations, while the proxy method approximates displacements and stress changes using reservoir simulation outputs and homogeneous geomechanical rock properties assumptions.
The proxy model has been applied to matrix- and fault-dominated systems, including the Brugge dataset. Results include pressure, temperature, displacements, stress changes predictions over 30 years. Fault stability is evaluated using Mohr-Coulomb criteria with a constant friction coefficient.
In fracture-dominated systems, faults often control flow but. Discrete Fracture Model (DFM) has been used for flow modelling.
Combining proxy and THM models can optimize the balance between accuracy and computational cost. The study emphasizes the differing impacts of pressure and temperature on fault stability during geothermal operations. ...
This study explores two approaches to assess stress changes: a semi-analytical geomechanical proxy and a fully-coupled Thermo-Hydro-Mechanical (THM) model using open-DARTS. The THM model simulates coupled thermal, hydraulic, and mechanical processes in complex rock formations, while the proxy method approximates displacements and stress changes using reservoir simulation outputs and homogeneous geomechanical rock properties assumptions.
The proxy model has been applied to matrix- and fault-dominated systems, including the Brugge dataset. Results include pressure, temperature, displacements, stress changes predictions over 30 years. Fault stability is evaluated using Mohr-Coulomb criteria with a constant friction coefficient.
In fracture-dominated systems, faults often control flow but. Discrete Fracture Model (DFM) has been used for flow modelling.
Combining proxy and THM models can optimize the balance between accuracy and computational cost. The study emphasizes the differing impacts of pressure and temperature on fault stability during geothermal operations.
This paper reports exposure of five different sealants to CO2-saturated water and wet supercritical CO2 at in-situ conditions (80 °C and 10 MPa). Three of the sealants investigated are based on Portland Cement, while the other two are based on Calcium Aluminate Cement, and a rock-based geopolymer specifically developed for Geological CO2 Storage (GCS). The five sealants were selected to represent different methods for improving wellbore seal integrity, such as restricting permeability (and porosity), or modifying how the material interacts with CO2-bearing fluids. Exposures were carried out in a purpose-built batch apparatus, enabling simultaneous exposure of up to 10 samples in total to CO2-saturated water and wet supercritical CO2.
After exposure, changes in the sealants’ microstructures and chemical and mineralogical compositions were assessed using scanning electron microscopy with energy-dispersive X-ray spectroscopy, computed tomography scanning, and fluid chemical analysis. The impact of exposure to CO2-bearing fluids was interpreted in terms of alteration and degradation of the materials, to compare how different sealant design modifications can be employed to enhance wellbore integrity. ...
This paper reports exposure of five different sealants to CO2-saturated water and wet supercritical CO2 at in-situ conditions (80 °C and 10 MPa). Three of the sealants investigated are based on Portland Cement, while the other two are based on Calcium Aluminate Cement, and a rock-based geopolymer specifically developed for Geological CO2 Storage (GCS). The five sealants were selected to represent different methods for improving wellbore seal integrity, such as restricting permeability (and porosity), or modifying how the material interacts with CO2-bearing fluids. Exposures were carried out in a purpose-built batch apparatus, enabling simultaneous exposure of up to 10 samples in total to CO2-saturated water and wet supercritical CO2.
After exposure, changes in the sealants’ microstructures and chemical and mineralogical compositions were assessed using scanning electron microscopy with energy-dispersive X-ray spectroscopy, computed tomography scanning, and fluid chemical analysis. The impact of exposure to CO2-bearing fluids was interpreted in terms of alteration and degradation of the materials, to compare how different sealant design modifications can be employed to enhance wellbore integrity.
Effectively mitigating induced seismicity in subsurface engineering operations within highly permeable, porous geo-energy reservoirs requires a clear understanding of how fluid injection parameters influence the seismic response. In this study, we performed injection-driven fault reactivation experiments on highly permeable saw-cut Red Felser sandstone to provide new insight into the effect of injection pattern and rate on fault slip behavior and seismicity evolution. Three different pressurization rates were applied: high, medium, and low rates of 2, 1, and 0.2 MPa/min, respectively. Three injection patterns were also used: cyclic recursive, monotonic, and stepwise injections. Our results reveal that a high pressurization rate leads to increased slip velocity, more microseismic events, higher total acoustic emission (AE) energy, and a lower b-value compared to tests with low pressurization rates. We postulate that a high pressurization rate enhances the likelihood of a sudden reduction in effective normal stress, leading to fault opening and the disruption of asperity contacts. Furthermore, results from samples subjected to various injection patterns demonstrate that the cyclic recursive pattern exhibits a higher maximum slip velocity, more episodes of slow slip, and greater radiated AE energy than a monotonic pattern. In the case of the cyclic recursive pattern, increasing the number of cycles increases shear stress drop, shear slip, and maximum slip velocity. Our findings suggest that using a monotonic injection pattern and low pressurization rate may mitigate seismicity on pre-existing faults in a highly permeable, porous reservoir.
Plain Language Summary
Human activities involving subsurface fluid injection projects, such as geothermal energy recovery and/or gas storage (CO2, H2 or methane), are widely acknowledged to cause earthquakes occasionally. This is a cause for public concern. Although several studies demonstrate that injection patterns and rates can play an essential role, the underlying physical mechanisms responsible for induced earthquakes still need to be better understood. Therefore, we performed laboratory tests on highly permeable Red Felser sandstone containing a simulated geological fault. We pumped water from the bottom of the sample using different pressurization rates and patterns while monitoring the effects on fault movement behavior. Our results showed that faster fluid injections tend to cause more rapid fault slips and generate more laboratory micro-earthquakes compared to slow injections. Among the injection patterns, the cyclic injection pattern resulted in the highest slip velocity and higher earthquake activity, indicating that the pattern of injection can impact fault movement. Our results can help improve the design of fluid injection projects to minimize the risk of inducing small earthquakes, especially in areas with pre-existing geological faults.
...
Effectively mitigating induced seismicity in subsurface engineering operations within highly permeable, porous geo-energy reservoirs requires a clear understanding of how fluid injection parameters influence the seismic response. In this study, we performed injection-driven fault reactivation experiments on highly permeable saw-cut Red Felser sandstone to provide new insight into the effect of injection pattern and rate on fault slip behavior and seismicity evolution. Three different pressurization rates were applied: high, medium, and low rates of 2, 1, and 0.2 MPa/min, respectively. Three injection patterns were also used: cyclic recursive, monotonic, and stepwise injections. Our results reveal that a high pressurization rate leads to increased slip velocity, more microseismic events, higher total acoustic emission (AE) energy, and a lower b-value compared to tests with low pressurization rates. We postulate that a high pressurization rate enhances the likelihood of a sudden reduction in effective normal stress, leading to fault opening and the disruption of asperity contacts. Furthermore, results from samples subjected to various injection patterns demonstrate that the cyclic recursive pattern exhibits a higher maximum slip velocity, more episodes of slow slip, and greater radiated AE energy than a monotonic pattern. In the case of the cyclic recursive pattern, increasing the number of cycles increases shear stress drop, shear slip, and maximum slip velocity. Our findings suggest that using a monotonic injection pattern and low pressurization rate may mitigate seismicity on pre-existing faults in a highly permeable, porous reservoir.
Plain Language Summary
Human activities involving subsurface fluid injection projects, such as geothermal energy recovery and/or gas storage (CO2, H2 or methane), are widely acknowledged to cause earthquakes occasionally. This is a cause for public concern. Although several studies demonstrate that injection patterns and rates can play an essential role, the underlying physical mechanisms responsible for induced earthquakes still need to be better understood. Therefore, we performed laboratory tests on highly permeable Red Felser sandstone containing a simulated geological fault. We pumped water from the bottom of the sample using different pressurization rates and patterns while monitoring the effects on fault movement behavior. Our results showed that faster fluid injections tend to cause more rapid fault slips and generate more laboratory micro-earthquakes compared to slow injections. Among the injection patterns, the cyclic injection pattern resulted in the highest slip velocity and higher earthquake activity, indicating that the pattern of injection can impact fault movement. Our results can help improve the design of fluid injection projects to minimize the risk of inducing small earthquakes, especially in areas with pre-existing geological faults.
Sealants that can guarantee long-term wellbore sealing integrity are of great significance to the safe and sustainable storage of CO2 in carbon capture and storage (CCS). In this study, we investigate how abrupt cyclic thermal shocks affect the integrity of four sealants of different compositions. These sealants include two reference OPC-based blends (S1 and S2), one newly-designed OPC-based blend that contains CO2-sequestering additives (S3), and one calcium aluminate cement (CAC)-based blend designed for CCS applications (S4). We have measured the thermal properties of these samples, followed by quenching and flow-through experiments to apply strong cyclic thermal shocks on samples of the four sealants, where we heated the samples to 120 °C, and quenched them in, or flowed through water of 20 °C. Using X-ray tomography (32 µm/voxel) before and after the experiment showed that both S1, S2 (reference OPC-based) and S4 (CAC-based) broke after thermal-shocking experiments. Cracks and new voids developed in the samples. Post-treatment strength testing shows that thermal shocks reduce the unconfined compressive strength of these three sealants. This implies that these compositions may not be optimal materials for long-term wellbore sealing during CO2 injection and storage afterward. For all these three sealant compositions, quenching resulted in a greater reduction in strength (by 53 % on average) than flow-through experiments (by 29 % on average). On the contrary, we have not observed any cracks after either quenching or flow-through experiments in S3 sealant (OPC with CO2-sequestering additives). We attribute the intactness of this sealant after thermal shocks to its higher thermal diffusivity than the other three sealants. Heat transfers more rapidly in this sealant and the associated thermal stresses are mild and insufficient to cause any damage to its integrity, which makes this sealant a good candidate for wellbore sealing material that can effectively withstand strong thermal shocks encountered during CCS, though further studies are required.
Induced earthquakes are still highly unpredictable, and often caused by variations in pore fluid pressure. Monitoring and understanding the mechanisms of fluid-induced fault slip is essential for seismic risk mitigation and seismicity forecasting. Fluid-induced slip experiments were performed on critically stressed faulted sandstone samples, and the evolution of the actively sent ultrasonic waves throughout the experiment was measured. Two different fault types were used: smooth saw-cut fault samples at a 35° angle, and a rough fault created by in situ faulting of the samples. Variations in the seismic slip velocity and friction along the fault plane were identified by the coda of the ultrasonic waves. Additionally, ultrasonic amplitudes show precursory signals to laboratory fault reactivation. Our results show that small and local variations in stress before fault failure can be inferred using coda wave interferometry for time-lapse monitoring, as coda waves are more sensitive to small perturbations in a medium than direct waves. Hence, these signals can be used as precursors to laboratory fault slip and to give insight into reactivation mechanisms. Our results show that time-lapse monitoring of coda waves can be used to monitor local stress changes associated with fault reactivation in this laboratory setting of fluid-induced fault reactivation. This is a critical first step toward a method for continuous monitoring of natural fault zones, contributing to seismic risk mitigation of induced and natural earthquakes.
In many geological systems, the porosity of rock or soil may evolve during mineral precipitation, a process that controls fluid transport properties. Here, we investigate the use of 4D neutron imaging to image flow and transport in Bentheim sandstone core samples before and after in-situ calcium carbonate precipitation. First, we demonstrate the applicability of neutron imaging to quantify the solute dispersion along the interface between heavy water and a cadmium aqueous solution. Then, we monitor the flow of heavy water within two Bentheim sandstone core samples before and after a step of in-situ mineral precipitation. The precipitation of calcium carbonate is induced by reactive mixing of two solutions containing CaCl2 and Na2CO3, either by injecting these two fluids one after each other (sequential experiment) or by injecting them in parallel (co-flow experiment). We use the contrast in neutron attenuation from time-resolved tomograms to derive three-dimensional fluid velocity field by using an inversion technique based on the advection-dispersion equation. Results show mineral precipitation induces a wider distribution of local flow velocities and leads to alterations in the main flow pathways. The flow distribution appears to be independent of the initial distribution in the sequential experiment, while in the co-flow experiment, we observed that higher initial local fluid velocities tended to increase slightly following precipitation. The outcome of this study contributes to progressing the knowledge in the domain of reactive solute and contaminant transport in the subsurface using the promising technique of neutron imaging.
This study presents a method to address the significant uncertainties in subsurface modeling that impact the efficiency of energy transition applications such as geothermal energy extraction and CO2 geological sequetsration. The approach combines a physics-based geomechanical proxy model with an ensemble smoother with multiple data assimilation (ES-MDA), aimed at enhancing uncertainty quantification through the integration of vertical displacement measurements from fluid production and injection. The data from wells is limited in spatial coverage, while these measurements offer extensive spatial information, improving the understanding of subsurface behavior by reflecting changes in reservoir pressure and temperature. The open-DARTS simulator for fluid flow and a geomechanical proxy are used to perform data assimilation with ES-MDA. By generating an ensemble of model realizations with varied permeability, calculating vertical displacements at the surface, and applying ES-MDA, we effectively identify the probability distribution of the vertical displacement of the model conditioned to observed subsidence data. Entropy is used as a statistical measure to quantify the reduction of uncertainty of subsurface models based on observations. Our approach was tested on a 2D conceptual and 3D realistic datasets, demonstrating its capability to provide data assimilation. This workflow represents an advancement in subsurface modeling, supporting informed decision-making in geothermal energy production and CO2 sequestration by offering an improved alternative for data assimilation and enhancing tools for uncertainty quantification.
Geothermal energy production often involves use of corrosion inhibitors. We performed rock mechanical experiments (room temperature; confining pressure of 10/20/30 MPa) on typical reservoir rocks (Bentheim sandstone and Treuchtlinger limestone) in contact with two different inhibitor solutions or with demineralized water. The sandstone experiments show no discernible difference in rock strength between inhibitors or water, attributed to low quartz reactivity. The limestone experiments show a significant difference in rock strength (and Mohr–Coulomb envelope), dependent on inhibitor type, attributed to high carbonate reactivity. This implies that, depending on the reactivity of the rocks and local stress conditions, inhibitor leakage may lead to unpredicted reservoir failure.
REFLECT DELIVERABLE D2.2
Report on gas solubility and degassing kinetic (type C)
Many rocks contain planar heterogeneities, in the form of open fractures, veins and/or stylolites, but scarce data exist on how strength and fracture pattern formation is affected by the presence of a singular planar heterogeneity in an otherwise uniform matrix. The mechanics of stylolite-bearing and/or fractured limestone is of interest to several engineering applications, from quarries to subsurface gas or geothermal reservoirs. We have performed Brazilian Disc tests on pre-fractured Indiana limestone samples and Treuchtlinger Marmor discs which contain cohesive stylolites, investigating Brazilian test Strength and the resulting fracture pattern. All experiments were filmed, and where possible analyzed with particle image velocimetry. When viewed in 2D, the planar discontinuity was set at different rotation angles compared to the principal loading direction, where perpendicular to the loading direction is defined as 0⁰. The results show that all samples are weaker than their intact counterparts. For the pre-fractured Indiana limestone, there is 10–75% angle-dependent weakening. However, in the samples with a stylolite, strength is weakened by 35–75%, independent of direction. Several new cracks appeared when fracturing a stylolite-sample, where the orientation is heavily influenced by the stylolite orientation. The fracture pattern and associated stress drops are more complex for high angles. In these samples always more than one fracture formed, whereas in pre-fractured samples usually only one new fracture formed. This suggests a potential for more permeability increase when hydrofracturing a stylolite-rich interval. Comparison with Finite Element Models indicates that this difference in fracture pattern is caused by the strength contrast between the anastomosing stylolite zone and the matrix material, leading to stress concentrations effects. This causes (micro-) fracture nucleation to occur locally, promotes fracture coalescence and fracture growth at lower overall sample-load conditions compared to intact samples.