Reactivity of carbon dioxide (CO2), rocks and brine is important in a number of practical situations in carbon dioxide sequestration. Injectivity of CO2 will be affected by near wellbore dissolution or precipitation. Natural fractures or faults containing specific minerals may reactivate leading to induced seismicity. In this project, we first examined if the reactions between CO2, brine and rocks affect the nature of the porous medium and properties including petrophysical properties in the timeframe of the injection operations. This was done by carrying out experiments at sequestration conditions (2000 psi for corefloods and 2400 psi for batch experiments, and 600°C) with three different types of rocks – sandstone, limestone and dolomite. Experiments were performed in batch mode and corefloods were conducted over a two-week period. Batch experiments were performed with samples of differing surface area to understand the impact of surface area on overall reaction rates. Toughreact, a reactive transport model was used to interpret and understand the experimental results. The role of iron in dissolution and precipitation reactions was observed to be significant. Iron containing minerals – siderite and ankerite dissolved resulting in changes in porosity and permeability. Corefloods and batch experiments revealed similar patterns. With the right cationic balance, there is a possibility of precipitation of iron bearing carbonates. The results indicate that during injection operations mineralogical changes may lead to injectivity enhancements near the wellbore and petrophysical changes elsewhere in the system. Limestone and dolomite cores showed consistent dissolution at the entrance of the core. The dissolution led to formation of wormholes and interconnected dissolution zones. Results indicate that near wellbore dissolution in these rock-types may lead to rock failure. Micro-CT images of the cores before and after the experiments revealed that an initial high-permeability pathway facilitated the formation of wormholes. The peak cation concentrations and general trends were matched using Toughreact. Batch reactor modeling showed that the geometric factors obtained using powder data that related effective surface area to the BET surface area had to be reduced for fractured samples and cores. This indicates that the available surface area in consolidated samples is lower than that deduced from powder experiments. Field-scale modeling of reactive transport and geomechanics was developed in parallel at Idaho National Laboratory. The model is able to take into account complex chemistry, and consider interactions of natural fractures and faults. Poroelastic geomechanical considerations are also included in the model.
Abstract This paper is the second part of a two part sequence on multiphysics algorithms and software. The first [1] focused on the algorithms; this part treats the multiphysics software framework and applications based on it. Tight coupling is typically designed into the analysis application at inception, as such an application is strongly tied to a composite nonlinear solver that arrives at the final solution by treating all equations simultaneously. The application must also take care to minimize both time and space error between the physics, particularly if more than one mesh representation is needed in the solution process. This paper presents an application framework that was specifically designed to support tightly coupled multiphysics analysis. The Multiphysics Object Oriented Simulation Environment (MOOSE) is based on the Jacobian-free Newton-Krylov (JFNK) method combined with physics-based preconditioning to provide the underlying mathematical structure for applications. The report concludes with the presentation of a host of nuclear, energy and environmental applications that demonstrate the efficacy of the approach and the utility of a well-designed multiphysics framework.
Subsurface radionuclide and trace metal contaminants throughout the U.S. Department of Energy (DOE) complex pose one of DOE’s greatest challenges for long-term stewardship. One promising in situ immobilization approach of these contaminants is engineered mineral (co)precipitation of calcite driven by urea hydrolysis that is catalyzed by enzyme urease. The tight nonlinear coupling among flow, transport, reaction and reaction-induced property changes of media of this approach was studied by reactive transport simulations with systematically increasing level of complexities of reaction network and physical/chemical heterogeneities using a numerical simulator named STOMP. Sensitivity studies on the reaction rates of both urea hydrolysis and calcite precipitation are performed via controlling urease enzyme concentration and precipitation rate constant according to the rate models employed. We have found that the rate of ureolysis is a dominating factor in the amount of precipitated mineral; however, the spatial distribution of the precipitates depends on both rates of ureolysis and calcite precipitation. A maximum 5% reduction in the porosity was observed within the simulation time period of 6 pore volumes in our 1-dimensional (1D) column simulations. When a low permeability inclusion is considered in the 2D simulations, the altered flow fields redistribute mineral forming constituents, leading to a distortedmore » precipitation reaction front. The simulations also indicate that mineral precipitation occurs along the boundary of the low permeability zone, which implies that contaminants in the low permeability zone could be encapsulated and isolated from the flow paths.« less
Ureolytically driven calcite precipitation is a promising approach for inducing subsurface mineral precipitation, but engineered application requires the ability to control and predict precipitate distribution. To study the coupling between reactant transport and precipitate distribution, columns with defined zones of immobilized urease were used to examine the distribution of calcium carbonate precipitation along the flow path, at two different initial flow rates. As expected, with slower flow precipitate was concentrated toward the upstream end of the enzyme zone and with higher flow the solid was more uniformly distributed over the enzyme zone. Under constant hydraulic head conditions the flow rate decreased as precipitates decreased porosity and permeability. The hydrolysis/precipitation zone was expected to become compressed in the upstream direction. However, apparent reductions in the urea hydrolysis rate and changes in the distribution of enzyme activity, possibly due to CaCO3 precipitate hindering urea transport to the enzyme, or enzyme mobilization, mitigated reaction zone compression. Co-injected strontium was expected to be sequestered by coprecipitation with CaCO3, but the results suggested that coprecipitation was not an effective sequestration mechanism in this system. In addition, spectral induced polarization (SIP) was used to monitor the spatial and temporal evolution of the reaction zone.
'/%3e%3cuse xlink:href='%23a' transform='rotate(23.036 .093 25.536)'/%3e%3cuse xlink:href='%23a' transform='rotate(45.87 1.273 16.18)'/%3e%3cuse xlink:href='%23a' transform='rotate(69.945 .996 12.078)'/%3e%3cuse xlink:href='%23a' transform='rotate(20.66 -19.689 31.932)'/%3e%3c/svg%3e)