Search

Geothermal energy is an infinite energy source for the present human society. Energy extraction from the deep subsurface requires engineering using a working fluid that circulates between well doublet. Due to its thermal properties, CO2 is an ideal option as a heat transfer fluid. By using CO2, working fluid loss is an advantage compared to other working fluids. This study developed a field-scale hydro-thermal model to examine the heat extraction potential from Soultz-sous-For\^ets with CO2 as the working fluid. Results are compared for the same scenario with water as the working fluid. A better understanding of the heat extraction mechanism is established by considering the reservoir response and the wellbore heat exchange. Sensitivity analyses are performed for different injection temperatures and flow rates for 50 years. Results show that the wellbore effect is multiple times higher than the reservoir response to the production temperature. Furthermore, lowering the injection temperature eventuates to a smaller temperature reduction at the subsurface, enhancing the overall heat extraction potential with a minor impact on thermal breakthrough. The cold region developed around the injection wellbore may affect the production fluid temperature due to its proximity to the production wellbore. To reach higher heat extraction efficiency, it is essential to use sufficient wellbore spacing. CO2 can be used as working fluid for over 50 years as it does not show significant thermal breakthrough and temperature plume evolution in the reservoir under studied conditions. CO2 shows lower temperature reduction for all injection rates and temperatures for 50 years of operation., Comment: 17 pages, 8 figures

Small-scale bed form topographies control hyporheic exchange and biogeochemical processes within aquatic sediments, which ultimately affect water quality and nutrient cycling at the watershed scale. The impact of three-dimensional (3D) and small-scale bed form topographies on hyporheic exchange and solute mixing is investigated in the present work. The effect of bed form morphologies on the development of zones of enhanced reaction rates (i.e., hot spots) is also studied. A computational fluid dynamics model to simulate river flow over bed forms is combined with a subsurface flow and multicomponent reactive solute transport model. A wide variety of bed form topographies are generated using geometric models by varying parameters controlling curvature as well as bed form wavelength and amplitude. Our results show that the pressure distribution at the sediment-water interface is strongly affected by the bed form geometry. Higher phase shifts in bed form shapes result in overall higher average velocity, larger zones of enhanced pressure and reaction rates, and higher amounts of solute exchange. Moreover, the bed form shapes control the reaction process for a wide range of sediment conductivities. This study advances the understanding of the effects of complex and small scale morphological features on hyporheic exchange processes.

We present the first experiments of dissolution-driven convection of carbon dioxide (CO$_2$) in a confined brine-saturated porous medium at high pressures. We designed a novel Hele-Shaw cell that allows for both visual and quantitative analyses, and address the effects of free-phase CO$_2$ and brine composition on convective dissolution. The visual examination of the gas volume combined with the measurement of pressure, which both evolve with dissolution, enable us to yield insights into the dynamics of convection in conditions that more closely reflect the geologic conditions. We find and analyze different dissolution events, including diffusive, early and late convection, and shut-down regimes. Our experiments reveal that in intermediate regime, a so-called `quasi-steady' state actually never happens. Dissolution flux continuously decreases in this regime, which is due to a negative feedback loop: the rapid reduction of pressure following convective dissolution, in turn, decreases the solubility of CO$_2$ at the gas-brine interface and thus the instability strength. We introduce a new scaling factor that not only compensates the flux reduction but also the nonlinearities that arise from different salt types. We present robust scaling relations for the compensated flux and for the transition times between consecutive regimes in systems with NaCl (Ra $\sim$ 3271--4841) and NaCl+CaCl$_2$ mixtures (Ra $\sim$ 2919--4283). We also find that NaCl+CaCl$_2$ mixtures enjoy a longer intermediate period before the shut-down of dissolution, but with a lower dissolution flux, as compared to NaCl brines...

We numerically examine solutal convection in porous media, driven by the dissolution of carbon dioxide (CO2) into water---an effective mechanism for CO2 storage in saline aquifers. Dissolution is associated with slow diffusion of free-phase CO2 into the underlying aqueous phase followed by density-driven convective mixing of CO2 throughout the water-saturated layer. We study the fluid dynamics of CO2 convection in the single aqueous-phase region. A comparison is made between two different boundary conditions in the top of the formation: (i) a constant, maximum aqueous-phase concentration of CO2, and (ii) a constant, low injection-rate of CO2, such that all CO2 dissolves instantly and the system remains in single phase. The latter model is found to involve a nonlinear evolution of CO2 composition and associated aqueous-phase density, which depend on the formation permeability. We model the full nonlinear phase behavior of water-CO2 mixtures in a confined domain, consider dissolution and fluid compressibility, and relax the common Boussinesq approximation. We discover new flow regimes and present quantitative scaling relations for global characteristics of spreading, mixing, and a dissolution flux in two- and three-dimensional media for both boundary conditions. We also revisit the scaling behavior of Sherwood number (Sh) with Rayleigh number (Ra), which has been under debate for porous-media convection. Our measurements from the solutal convection in the range 1, 500

We present a new implicit higher-order finite element (FE) approach to efficiently model compressible multicomponent fluid flow on unstructured grids and in fractured porous subsurface formations. The scheme is sequential implicit: pressures and fluxes are updated with an implicit Mixed Hybrid Finite Element (MHFE) method, and the transport of each species is approximated with an implicit second-order Discontinuous Galerkin (DG) FE method. Discrete fractures are incorporated with a cross-flow equilibrium approach. This is the first investigation of all-implicit higher-order MHFE-DG for unstructured triangular, quadrilateral (2D), and hexahedral (3D) grids and discrete fractures. A lowest-order implicit finite volume (FV) transport update is also developed for the same grid types. The implicit methods are compared to an Implicit-Pressure-Explicit-Composition (IMPEC) scheme. For fractured domains, the unconditionally stable implicit transport update is shown to increase computational efficiency by orders of magnitude as compared to IMPEC, which has a time-step constraint proportional to the pore volume of discrete fracture grid cells. However, when lowest-order Euler time-discretizations are used, numerical errors increase linearly with the larger implicit time-steps, resulting in high numerical dispersion. Second-order Crank-Nicolson implicit MHFE-DG and MHFE-FV are therefore presented as well. Convergence analyses show twice the convergence rate for the DG methods as compared to FV, resulting in two to three orders of magnitude higher computational efficiency. Numerical experiments demonstrate the efficiency and robustness in modeling compressible multicomponent flow on irregular and fractured 2D and 3D grids, even in the presence of fingering instabilities., Comment: 24 pages, 3 tables, 10 figures

We studied the effects on oil sweep efficiency of the proportion, hierarchical organization, and connectivity of high-permeability open-framework conglomerate (OFC) cross-sets within the multi-scale stratal architecture found in fluvial deposits. Utilizing numerical simulations and the RVA/Paraview open-source visualization package, we analyzed oil production rate, water breakthrough time, and spatial and temporal distribution of residual oil saturation. The effective permeability of the reservoir exhibits large-scale anisotropy created by the organization of OFC cross-sets within unit bars, and the organization of unit bars within compound bars. As a result oil sweep efficiency critically depends on the direction of the pressure gradient. When pressure gradient is oriented normal to paleoflow direction, the total oil production and the water breakthrough time are larger, and remaining oil saturation is smaller. This result is found regardless of the proportion or connectivity of the OFC cross-sets, within the ranges examined. Contrary to expectations, the total amount of trapped oil due to the effect of capillary trapping does not depend on the pressure gradient within the examined range. Hence the pressure difference between production and injection wells does not affect sweep efficiency, although the spatial distribution of oil remaining in the reservoir depends on this value. Whether or not clusters of connected OFC span the domain does not affect sweep efficiency, only the absolute rate of oil production. The RVA/Paraview application allowed us to visualize and examine these non-intuitive results., Comment: 31 pages, 17 figures. arXiv admin note: text overlap with arXiv:1405.3998

Physical and chemical heterogeneities have a large impact on reactive transport in porous media. Examples of heterogeneous attributes affecting reactive mass transport are the hydraulic conductivity (K), and the equilibrium sorption distribution coefficient (Kd). This paper uses the Deng et al. (2013) conceptual model for multimodal reactive mineral facies and a Lagrangian-based stochastic theory in order to analyze the reactive solute dispersion in three-dimensional anisotropic heterogeneous porous media with hierarchical organization of reactive minerals. An example based on real field data is used to illustrate the time evolution trends of reactive solute dispersion. The results show that the correlation between the hydraulic conductivity and the equilibrium sorption distribution coefficient does have a significant effect on reactive solute dispersion. The anisotropy ratio does not have a significant effect on reactive solute dispersion. Furthermore, through a sensitivity analysis we investigate the impact of changing the mean, variance, and integral scale of K and Kd on reactive solute dispersion.