Dissolution Trapping of Carbon Dioxide in Heterogeneous Aquifers.

The geologic architecture in sedimentary reservoirs affects the behavior of density-driven flow and the dispersion of CO2-rich brine. The spatial organization and connectivity of facies types play an important role. Low-permeability facies may suppress fingering and reduce vertical spreading, but may also increase transverse mixing. This is more pronounced when geologic structures create preferential flow pathways through connected facies types. We perform high-resolution simulations of three-dimensional (3D) heterogeneous formations whose connectivity cannot be represented in two-dimensional models consistent with percolation theory. This work focuses on the importance of 3D facies-based heterogeneity and connectivity on advection-diffusion transport of dissolved CO2. Because the dissolution of CO2 and the subsequent density increase of brine are the driving force for gravitational instabilities, we model the phase behavior with the accurate cubic-plus-association equation-of-state, which accounts for the self-association of polar water molecules and the cross-association between CO2 and water. Our results elucidate how the spatial organization of facies affects the dynamics of CO2 convective mixing. Scaling relations for the evolution of a global dispersion-width provide insights that can be universally applied. The results suggest that the long-term evolution and scaling of dispersion are surprisingly similar for homogeneous and (binary and multiscale) heterogeneous porous media.

Geographic axes and the persistence of cultural diversity.

Jared Diamond's Guns, Germs, and Steel [Diamond J, (1997) Guns, Germs, and Steel (WW Norton, NY)] has provided a scientific foundation for answering basic questions, such as why Eurasians colonized the global South and not the other way around, and why there is so much variance in economic development across the globe. Diamond's explanatory variables are: (i) the susceptibility of local wild plants to be developed for self-sufficient agriculture; (ii) the domesticability of large wild animals for food, transport, and agricultural production; and (iii) the relative lengths of the axes of continents with implications for the spread of human populations and technologies. This third "continental axis" thesis is the most difficult of Diamond's several explanatory factors to test, given that the number of continents are too few for statistical analysis. This article provides a test of one observable implication of this thesis, namely that linguistic diversity should be more persistent to the degree that a geographic area is oriented more north-south than east-west. Using both modern states and artificial geographic entities as the units of analysis, the results provide significant confirmation of the relationship between geographic orientation and cultural homogenization. Beyond providing empirical support for one observable implication of the continental axis theory, these results have important implications for understanding the roots of cultural diversity, which is an important determinant of economic growth, public goods provision, local violence, and social trust.

Geometry-Driven Detection, Tracking and Visual Analysis of Viscous and Gravitational Fingers.

Viscous and gravitational flow instabilities cause a displacement front to break up into finger-like fluids. The detection and evolutionary analysis of these fingering instabilities are critical in multiple scientific disciplines such as fluid mechanics and hydrogeology. However, previous detection methods of the viscous and gravitational fingers are based on density thresholding, which provides limited geometric information of the fingers. The geometric structures of fingers and their evolution are important yet little studied in the literature. In this article, we explore the geometric detection and evolution of the fingers in detail to elucidate the dynamics of the instability. We propose a ridge voxel detection method to guide the extraction of finger cores from three-dimensional (3D) scalar fields. After skeletonizing finger cores into skeletons, we design a spanning tree based approach to capture how fingers branch spatially from the finger skeletons. Finally, we devise a novel geometric-glyph augmented tracking graph to study how the fingers and their branches grow, merge, and split over time. Feedback from earth scientists demonstrates the usefulness of our approach to performing spatio-temporal geometric analyses of fingers.

A higher-order finite element reactive transport model for unstructured and fractured grids.

This work presents a new reactive transport framework that combines a powerful geochemistry engine with advanced numerical methods for flow and transport in subsurface fractured porous media. Specifically, the PhreeqcRM interface (developed by the USGS) is used to take advantage of a large library of equilibrium and kinetic aqueous and fluid-rock reactions, which has been validated by numerous experiments and benchmark studies. Fluid flow is modeled by the Mixed Hybrid Finite Element (FE) method, which provides smooth velocity fields even in highly heterogenous formations with discrete fractures. A multilinear Discontinuous Galerkin FE method is used to solve the multicomponent transport problem. This method is locally mass conserving and its second order convergence significantly reduces numerical dispersion. In terms of thermodynamics, the aqueous phase is considered as a compressible fluid and its properties are derived from a Cubic Plus Association (CPA) equation of state. The new simulator is validated against several benchmark problems (involving, e.g., Fickian and Nernst-Planck diffusion, isotope fractionation, advection-dispersion transport, and rock-fluid reactions) before demonstrating the expanded capabilities offered by the underlying FE foundation, such as high computational efficiency, parallelizability, low numerical dispersion, unstructured 3D gridding, and discrete fraction modeling.

Numerical Modeling of Methane Leakage from a Faulty Natural Gas Well into Fractured Tight Formations.

Horizontal drilling and hydraulic fracturing have enabled hydrocarbon recovery from unconventional reservoirs, but led to natural gas contamination of shallow groundwaters. We describe and apply numerical models of gas-phase migration associated with leaking natural gas wells. Three leakage scenarios are simulated: (1) high-pressure natural gas pulse released into a fractured aquifer; (2) continuous slow leakage into a tilted fractured formation; and (3) continuous slow leakage into an unfractured aquifer with fluvial channels, to facilitate a generalized evaluation of natural gas transport from faulty natural gas wells. High-pressure pulses of gas leakage into sparsely fractured media are needed to produce the extensive and rapid lateral spreading of free gas previously observed in field studies. Transport in fractures explains how methane can travel vastly different distances and directions laterally away from a leaking well, which leads to variable levels of methane contamination in nearby groundwater wells. Lower rates of methane leakage (≤1 Mcf/day) produce shorter length scales of gas transport than determined by the high-pressure scenario or field studies, unless aquifers have low vertical permeabilities (≤1 millidarcy) and fractures and bedding planes have sufficient tilt (∼10°) to allow a lateral buoyancy component. Similarly, in fractured rock aquifers or where permeability is controlled by channelized fluvial deposits, lateral flow is not sufficiently developed to explain fast-developing gas contamination (0-3 months) or large length scales (∼1 km) documented in field studies. Thus, current efforts to evaluate the frequency, mechanism, and impacts of natural gas leakage from faulty natural gas wells likely underestimate contributions from small-volume, low-pressure leakage events.

Impacts of Methane on Carbon Dioxide Storage in Brine Formations.

In the context of geological carbon sequestration (GCS), carbon dioxide (CO2 ) is often injected into deep formations saturated with a brine that may contain dissolved light hydrocarbons, such as methane (CH4 ). In this multicomponent multiphase displacement process, CO2 competes with CH4 in terms of dissolution, and CH4 tends to exsolve from the aqueous into a gaseous phase. Because CH4 has a lower viscosity than injected CO2 , CH4 is swept up into a 'bank' of CH4 -rich gas ahead of the CO2 displacement front. On the one hand, this may provide a useful tracer signal of an approaching CO2 front. On the other hand, the emergence of gaseous CH4 is undesirable because it poses a leakage risk of a far more potent greenhouse gas than CO2 if the cap rock is compromised. Open fractures or faults and wells could result in CH4 contamination of overlying groundwater aquifers as well as surface emissions. We investigate this process through detailed numerical simulations for a large-scale GCS pilot project (near Cranfield, Mississippi) for which a rich set of field data is available. An accurate cubic-plus-association equation-of-state is used to describe the non-linear phase behavior of multiphase brine-CH4 -CO2 mixtures, and breakthrough curves in two observation wells are used to constrain transport processes. Both field data and simulations indeed show the development of an extensive plume of CH4 -rich (up to 90 mol%) gas as a consequence of CO2 injection, with important implications for the risk assessment of future GCS projects.

Hydrocarbon-Rich Groundwater above Shale-Gas Formations: A Karoo Basin Case Study.

Horizontal drilling and hydraulic fracturing have enhanced unconventional hydrocarbon recovery but raised environmental concerns related to water quality. Because most basins targeted for shale-gas development in the USA have histories of both active and legacy petroleum extraction, confusion about the hydrogeological context of naturally occurring methane in shallow aquifers overlying shales remains. The Karoo Basin, located in South Africa, provides a near-pristine setting to evaluate these processes, without a history of conventional or unconventional energy extraction. We conducted a comprehensive pre-industrial evaluation of water quality and gas geochemistry in 22 groundwater samples across the Karoo Basin, including dissolved ions, water isotopes, hydrocarbon molecular and isotopic composition, and noble gases. Methane-rich samples were associated with high-salinity, NaCl-type groundwater and elevated levels of ethane, 4 He, and other noble gases produced by radioactive decay. This endmember displayed less negative δ13 C-CH4 and evidence of mixing between thermogenic natural gases and hydrogenotrophic methane. Atmospheric noble gases in the methane-rich samples record a history of fractionation during gas-phase migration from source rocks to shallow aquifers. Conversely, methane-poor samples have a paucity of ethane and 4 He, near saturation levels of atmospheric noble gases, and more negative δ13 C-CH4 ; methane in these samples is biogenic and produced by a mixture of hydrogenotrophic and acetoclastic sources. These geochemical observations are consistent with other basins targeted for unconventional energy extraction in the USA and contribute to a growing data base of naturally occurring methane in shallow aquifers globally, which provide a framework for evaluating environmental concerns related to unconventional energy development (e.g., stray gas).

Critical Dynamics of Gravito-Convective Mixing in Geological Carbon Sequestration.

When CO2 is injected in saline aquifers, dissolution causes a local increase in brine density that can cause Rayleigh-Taylor-type gravitational instabilities. Depending on the Rayleigh number, density-driven flow may mix dissolved CO2 throughout the aquifer at fast advective time-scales through convective mixing. Heterogeneity can impact density-driven flow to different degrees. Zones with low effective vertical permeability may suppress fingering and reduce vertical spreading, while potentially increasing transverse mixing. In more complex heterogeneity, arising from the spatial organization of sedimentary facies, finger propagation is reduced in low permeability facies, but may be enhanced through more permeable facies. The connectivity of facies is critical in determining the large-scale transport of CO2-rich brine. We perform high-resolution finite element simulations of advection-diffusion transport of CO2 with a focus on facies-based bimodal heterogeneity. Permeability fields are generated by a Markov Chain approach, which represent facies architecture by commonly observed characteristics such as volume fractions. CO2 dissolution and phase behavior are modeled with the cubic-plus-association equation-of-state. Our results show that the organization of high-permeability facies and their connectivity control the dynamics of gravitationally unstable flow. We discover new flow regimes in both homogeneous and heterogeneous media and present quantitative scaling relations for their temporal evolution.