Leakage of stored bulk phase CO₂ is one risk for sequestration in deep saline aquifers. As the less dense CO₂ migrates upward within a storage formation or in layers above the formation, the security of its storage depends upon the trapping mechanisms that counteract the migration. The trapping mechanism motivating this research is local capillary trapping (LCT), which occurs during buoyancy-driven migration of bulk phase CO₂ within a saline aquifer with spatially heterogeneous petrophysical properties. When a CO₂ plume rising by buoyancy encounters a region where capillary entry pressure is locally larger than average, CO₂ accumulates beneath the region. One benefit of LCT, applied specifically to CO₂ sequestration and storage, is that saturation of stored CO₂ phase is larger than the saturation for other permanent trapping mechanisms. Another potential benefit is security: CO₂ that occupies local capillary traps remains there, even if the overlying formation that provides primary containment were to be compromised and allow leakage. Most work on LCT has involved numerical simulation (Saadatpoor 2010, Ganesh 2012); the research work presented here is a step toward understanding local capillary trapping at the bench scale. An apparatus and set of fluids are described which allow examining the extent of local capillary trapping, i.e. buoyant nonwetting phase immobilization beneath small-scale capillary barriers, which can be expected in typical heterogeneous storage formation. The bench scale environment analogous to CO₂ and brine in a saline aquifer is created in a quasi-two dimensional experimental apparatus with dimension of 63 cm by 63 cm by 5 cm, which allows for observation of plume migration with physically representative properties but at experimentally convenient ambient conditions. A surrogate fluid pair is developed to mimic the density, viscosity and interfacial tension relationship found at pressure and temperature typical of storage aquifers. Porous media heterogeneity, pressure boundary conditions, migration modes of uprising nonwetting phase, and presence of fracture/breach in the capillary barrier are studied in series of experiments for their influences on LCT. A variety of heterogeneous porous media made of a range of sizes of loosely packed silica beads are used to validate and test the persistence of local capillary trapping mechanism. By adjusting the boundary conditions (fluid levels in reservoirs attached to top and to bottom ports of the apparatus), the capillary pressure gradient across the domain was manipulated. Experiments were conducted with and without the presence of fracture/potential leakage pathway in the capillary seal. The trapped buoyant phase remained secure beneath the local capillary barriers, as long as the effective capillary pressure exerted by the trapped phase (proportional to column height of the phase) is smaller than the capillary entry pressure of the barrier. The local capillary trapping mechanism remained persistent even under forced imbibition, in which a significantly higher hydraulic potential gradient, and therefore a larger gradient in capillary pressure, was applied to the system. The column height of buoyant fluid that remained beneath the local capillary barrier was smaller by a factor corresponding to the increase in capillary pressure gradient. Mimicking a breach of the caprock by opening valves at the top of the apparatus allowed buoyant mobile phase held beneath the valves to escape, but buoyant phase held in local traps at saturations above residual, and therefore potentially mobile, was undisturbed. This work provides systematic validation of a novel concept, namely the long-term security of CO₂ that fills local (small-scale) capillary traps in heterogeneous storage formations. Results from this work reveal the first ever unequivocal experimental evidence on persistence of local capillary trapping mechanism. Attempts to quantify the nonwetting phase saturation and extent of LCT persistence serve as the initial steps to potentially reduce the risks associated with long-term storage security. / text

This work examines the impact of a viscosity force parameter, fluid velocity, and a capillary force parameter, interfacial tension, on the saturation, morphology, and topology of NW fluid in Bentheimer sandstone after primary imbibition, drainage, and secondary imbibition. Brine and air (used as a proxy for supercritical CO₂) flow experiments were performed on 6 mm diameter Bentheimer cores and were quantified via imaging with x-ray computed microtomography (x-ray CMT), which allows for three dimensional, non-destructive, pore-scale analysis of the amount and distribution of NW phase fluid within the sandstone cores. It was found that trapped NW phase saturation decreases with increases in capillary number, average blob size decreases with increases in capillary number, and the number of NW blobs increases with increases in capillary number. In addition, it was found that NW phase trapping is dependent on initial NW phase connectivity within the porous medium; with more negative values of initial NW phase Euler number resulting in less trapping. We suggest that the Euler number-saturation and the capillary number-saturation relationships for a given medium should be taken into consideration when designing a CO₂ sequestration scenario. / Graduation date: 2013

After the injection of CO₂ into a subsurface formation, various storage mechanisms help immobilize the CO₂. Injection strategies that promote the buoyant movement of CO₂ during the post-injection period can increase immobilization by the mechanisms of dissolution and residual phase trapping. In this work, we argue that the heterogeneity intrinsic to sedimentary rocks gives rise to another category of trapping, which we call local capillary trapping. In a heterogeneous storage formation where capillary entry pressure of the rock is correlated with other petrophysical properties, numerous local capillary barriers exist and can trap rising CO₂ below them. The size of barriers depends on the correlation length, i.e., the characteristic size of regions having similar values of capillary entry pressure. This dissertation evaluates the dynamics of the local capillary trapping and its effectiveness to add an element of increased capacity and containment security in carbon storage in heterogeneous permeable media. The overall objective is to obtain the rigorous assessment of the amount and extent of local capillary trapping expected to occur in typical storage formations. A series of detailed numerical simulations are used to quantify the amount of local capillary trapping and to study the effect of local capillary barriers on CO₂ leakage from the storage formation. Also, a research code is developed for finding clusters of local capillary trapping from capillary entry pressure field based on the assumption that in post-injection period the viscous forces are negligible and the process is governed solely by capillary forces. The code is used to make a quantitative assessment of an upper bound for local capillary trapping capacity in heterogeneous domains using the geologic data, which is especially useful for field projects since it is very fast compared to flow simulation. The results show that capillary heterogeneity decreases the threshold capacity for non-leakable storage of CO₂. However, in cases where the injected volume is more than threshold capacity, capillary heterogeneity adds an element of security to the structural seal, regardless of how CO₂ is accumulated under the seal, either by injection or by buoyancy. In other words, ignoring heterogeneity gives the worst-case estimate of the risk. Nevertheless, during a potential leakage through failed seals, a range of CO₂ leakage amounts may occur depending on heterogeneity and the location of the leak. In geologic CO₂ storage in typical saline aquifers, the local capillary trapping can result in large volumes that are sufficiently trapped and immobilized. In fact, this behavior has significant implications for estimates of permanence of storage, for assessments of leakage rates, and for predicting ultimate consequences of leakage. / text

Geological storage of carbon dioxide (CO2) in deep saline aquifers is one approach to mitigate release from large point sources to the atmosphere. Understanding of in-situ processes providing trapping is important to the development of realistic models and the planning of future storage projects. This thesis covers both field- and pore-scale numerical modeling studies of such geohydrological processes, with focus on residual trapping. The setting is a CO2-injection experiment at the Heletz test site, conducted within the frame of the EU FP7 MUSTANG and TRUST projects. The objectives of the thesis are to develop and analyze alternative experimental characterization test sequences for determining in-situ residual CO2 saturation (Sgr), as well as to analyze the impact of the injection strategy on trapping, the effect of model assumptions (coupled wellbore-reservoir flow, geological heterogeneity, trapping model) on the predicted trapping, and to develop a pore-network model (PNM) for simulating and analyzing pore-scale mechanisms. The results include a comparison of alternative characterization test sequences for estimating Sgr. The estimates were retrieved through parameter estimation. The effect on the estimate of including various data sets was determined. A new method, using withdrawal and an indicator-tracer, for obtaining a residual zone in-situ was also introduced. Simulations were made of the CO2 partitioning between layers in a multi-layered formation, and parameters influencing this were identified. The results showed the importance of accounting for coupled wellbore-reservoir flow in simulations of such scenarios. Simulations also showed that adding chase-fluid stages after a conventional CO2 injection enhances the (residual and dissolution) trapping. Including geological heterogeneity generally decreased the estimated trapping. The choice of trapping model may largely effect the quantity of the predicted residual trapping (although most of them produced similar results). The use of an appropriate trapping model and description of geological heterogeneity for a site when simulating CO2 sequestration is vital, as different assumptions may give significant discrepancies in predicted trapping. The result also includes a PNM code, for multiphase quasi-static flow and trapping in porous materials. It was used to investigate trapping and obtain an estimated trapping (IR) curve for Heletz sandstone.