Advective Gas Transport Research Articles

Derivation of heterogeneous material distributions and their sensitivity to HM-coupled two-phase flow models exemplified with the LASGIT experiment

Advective gas transport in bentonite, a possible buffer material in repositories for radioactive materials, is difficult to simulate in numerical continuum models, partly due to the complicated microstructure of bentonite. To generate reliable models of repositories nevertheless, spatially distributed heterogeneous material properties can be used to allow localization of gas flow. In this study, a pore-size-dependent stochastic approach of the gas entry pressure is derived from Mercury Intrusion Porosimetry, which is used to replicate measurements from the LASGIT experiment. In addition, three benchmark tests are simulated to investigate the dependence of heterogeneous distributions of material properties on the mesh discretization, the temporal dependence, and the coupling between the processes influenced by the heterogeneous parameters. The numerical modeling results of the LASGIT experiment show that the onset of gas flow into the system and the subsequent increase in pressure and stress can be well reproduced using heterogeneous distributions. Compared to a model with homogeneous material properties, heterogeneous distributions may allow the generation of dilatancy-controlled microfractures—an important feature with regard to the advective gas flow in bentonites. However, it can be observed that the heterogeneous distributions in LASGIT are less significant, as technical gaps or differences in material types could have a greater impact.

Hydro-mechanical modelling of gas transport in clayey host rocks for nuclear waste repositories

Host rocks, as the final impediment to waste migration, play a significant role in the nuclear waste repositories. Modelling of gas migration in saturated host rocks as well as its coupled hydro-mechanical (HM) behaviors is of importance to the assessment of geological disposal facilities. A comprehensive literature review is carried out on the models for simulating gas migration in saturated rock materials. Several aspects have been provided and discussed, including the material properties and experimental interpretations, governing equations, constitutive models for hydraulic and mechanical processes, fracture propagation models. Specifically, these models are discussed in detail with respect to their performance in describing the recorded experimental behaviors. It is found that the visco-capillary two phase flow model with enriched intrinsic permeability is commonly used to describe the advective gas transport in saturated host rocks. The embedded fracture model (EFM) or enriched EFM seems to be the most favored model as it accounts for the fracturing mechanism which is more representative to implicitly simulate the preferential gas pathways. To describe the mechanical behavior of rocks during the gas migration process, linear/nonlinear elastic, elastoplastic and damage models have been included in the mechanical processes. The HM models incorporating plasticity or damage may not be applicable in most experimental studies, which are not competent to simulate the key experimental behavior associated with the development of gas preferential pathways. Current models that can explicitly describe the gas induced micro-fracturing are seldomly reported, the existed ones are not able to represent all the key experimental behaviors related to preferential gas flow. Advanced approaches, i.e., phase field (PF) method, extended finite element method (XFEM), discrete element method (DEM), hybrid finite discrete element method (FDEM) can be integrated with cohesive zone model (CZM) to provide some promising aspects in representing the development of preferential gas pathways. Lastly, the conclusions and recommendations for future modelling are given.

A fully implicit coupled pore-network/free-flow model for the pore-scale simulation of drying processes

We present a fully coupled pore-network/free-flow model providing pore-scale insight into drying processes. We solve the Navier–Stokes equations with component transport in the free-flow region, coupled to a dynamic two-phase, two-component pore-network model (PNM) in the porous domain. The dynamic multi-physics model allows to temporally resolve drying processes in-between capillary equilibrium states. All simulations are non-isothermal and use pressure- and temperature-depended fluid properties. Carefully chosen coupling conditions and a monolithic solver ensure local conservation of mass, momentum, and energy fluxes, in particular at the interface between both model domains. We solve for wetting and non-wetting fluid pressure fields and consider advective gas transport in the network. Numerical examples demonstrate that the coupled model is able to cover a wide range of physical processes relevant for drying and show the mutual interaction of the two subregions. The model is implemented in the modular open-source framework DuMu such that extensions are straight-forward.

Insights into gas migration behavior in saturated GMZ bentonite under flexible constraint conditions

Gas migration behavior of bentonite-based buffer/backfill materials is one of the most important issues for safety assessment of deep geological repository for disposal of high-level radioactive wastes. In this paper, a series of gas injection tests were conducted on water-saturated bentonite specimens under stepwise increasing gas injection pressure and a flexible boundary condition. Prior to the gas injection tests, the steady state water permeability (stst kw) was determined by the water injection tests with deionized water and the stst kw values obtained were in an order of magnitude of 10-19-10-20 m2, which were slightly lower than the non-steady state ones calculated from the results obtained at the initial stage of gas injection tests. The effective gas permeability measured before and after gas breakthrough on the initially water-saturated bentonite specimens with different dry densities ranges between 5.45 × 10-25 m2 to 2.76 × 10-22 m2. The Peclet number (Pe) was adopted to determine the relative importance of advective gas transport (pressure-driven gas flow) to diffusive gas transport through saturated bentonite specimens. Calculation shows that the Pe values are significantly lower than 1, indicating the dominant role of gas diffusion across the specimens before and after gas breakthrough. As dry density of the specimen increases from 1.3 g/cm3 to 1.7 g/cm3, the gas breakthrough pressure increases from 4.92 MPa to 8.68 MPa. Meanwhile, the residual pressure differences measured after gas breakthrough ranges from 4.04 MPa to 7.88 MPa for specimens with different dry densities. A slight decrease of water content of the specimens experienced gas injection tests and a stable increase of specimen volume at the breakthrough indicate that both the displacement of water from interconnected pathways (capillary effect) and the pathway dilation (mechanical effect) play important roles in gas breakthrough for specimen tested under flexible conditions.

Influence of total organic carbon deposition on the inventory of gas hydrate in the Indian continental margins

Total organic carbon (TOC) content of marine sediments represents residual carbon, originally derived from terrestrial and marine sources, which has survived seafloor and shallow subseafloor diagenesis. Ultimately, its preservation below the sulfate reduction zone in marine sediments drives methanogenesis. Within the gas hydrate stability zone (GHSZ), methane production along continental margins can supersaturate pore fluids and lead to the formation of gas hydrate. In this paper we examine the inventory and sources of TOC in sediments collected from four regions within the GHSZ along the Indian continental margins. The recovered sediments vary in age from Oligocene to recent. Mean TOC abundance is greatest in the Krishna–Godavari (K–G) Basin and decreases progressively to the Mahanadi basin, Andaman wedge, and Kerala–Konkan (K–K) Basin. This decrease in TOC is matched by a progressive increase in biogenic CaCO3 and increasing distance from terrestrial sources of organic matter and lithogenic materials. Organic carbon sources inferred from C/N and δ13CTOC range from terrestrial (K–G Basin) to mixed marine and terrestrial (Mahanadi Basin), to marine dominant (Andaman wedge and K–K Basin). In the K–G Basin, variation in the bulk δ13CTOC is consistent with changes in C3 and C4 vegetation driven by monsoon variability on glacial-interglacial timescales, whereas in the Mahanadi Basin a shift in the δ13CTOC likely reflects the onset of C4 plant deposition in the Late Miocene. A large shift the δ13CTOC in the K–K basin is consistent with a change from C3 to C4 dominated plants during the middle Miocene. We observe a close relationship between TOC content and gas hydrate saturation, but consider the role of sedimentation rates on the preservation of TOC in the zone of methanogenesis and advective flow of methane from depth. Although TOC contents are sufficient for in situ methanogenesis at all the sites where gas hydrates were observed or inferred from proxy data, seismic, borehole log, pressure core, and gas composition data coupled with relatively high observed gas hydrate saturations suggest that advective gas transport may also play a role in the saturation of methane and the formation of gas hydrates in these regions. Although TOC content may be a first order indicator for gas hydrate potential, the structural and stratigraphic geologic environment along a margin will most likely dictate where the greatest gas hydrate saturations will occur.

Numerical Model Investigation for Potential Methane Explosion and Benzene Vapor Intrusion Associated with High-Ethanol Blend Releases

Ethanol-blended fuel releases usually stimulate methanogenesis in the subsurface, which could pose an explosion risk if methane accumulates in a confined space above the ground where ignitable conditions exist. Ethanol-derived methane may also increase the vapor intrusion potential of toxic fuel hydrocarbons by stimulating the depletion of oxygen by methanotrophs, and thus inhibiting aerobic biodegradation of hydrocarbon vapors. To assess these processes, a three-dimensional numerical vapor intrusion model was used to simulate the degradation, migration, and intrusion pathway of methane and benzene under different site conditions. Simulations show that methane is unlikely to build up to pose an explosion hazard (5% v/v) if diffusion is the only mass transport mechanism through the deeper vadose zone. However, if methanogenic activity near the source zone is sufficiently high to cause advective gas transport, then the methane indoor concentration may exceed the flammable threshold under simulated conditions. During subsurface migration, methane biodegradation could consume soil oxygen that would otherwise be available to support hydrocarbon degradation, and increase the vapor intrusion potential for benzene. Vapor intrusion would also be exacerbated if methanogenic activity results in sufficiently high pressure to cause advective gas transport in the unsaturated zone. Overall, our simulations show that current approaches to manage the vapor intrusion risk for conventional fuel released might need to be modified when dealing with some high ethanol blend fuel (i.e., E20 up to E95) releases.

Evidence for gas-induced pathways in clay using a nanoparticle injection technique

AbstractCorrosion, water radiolysis and microbial degradation will result in the generation of gas within repositories designed for the geological disposal of high-level radioactive waste. It is therefore crucial in the design of such facilities that the relevant mechanisms allowing gas migration through repository materials, both engineered barriers and clay-based candidate host rocks, are correctly identified. In Belgium, the Boom Clay represents a candidate host material for which the advective gas breakthrough characteristics and transport properties have been extensively tested and are well defined by numerous studies. The Boom Clay displays a significant capacity for self-sealing and both laboratory and field tests indicate that advective gas transport occurs not by visco-capillary flow, but instead through the formation of pressure-induced dilatant pathways. In this study, we present results from a gas injection test designed to demonstrate the presence of these features by injecting nanoparticulate tracers with helium gas into a sample of Boom Clay. The results provide conclusive evidence for the formation of transient, dilatant gas pathways within a candidate clay-based host rock. This technique provides a novel diagnostic tool for the identification of processes governing multi-phase flow, supporting robust long-term assessments of repository performance.

Surface heat flow and CO2 emissions within the Ohaaki hydrothermal field, Taupo Volcanic Zone, New Zealand

Abstract Carbon dioxide emissions and heat flow have been determined from the Ohaaki hydrothermal field, Taupo Volcanic Zone (TVZ), New Zealand following 20 a of production (116 MW e ). Soil CO 2 degassing was quantified with 2663 CO 2 flux measurements using the accumulation chamber method, and 2563 soil temperatures were measured and converted to equivalent heat flow (W m −2 ) using published soil temperature heat flow functions. Both CO 2 flux and heat flow were analysed statistically and then modelled using 500 sequential Gaussian simulations. Forty subsoil CO 2 gas samples were also analysed for stable C isotopes. Following 20 a of production, current CO 2 emissions equated to 111 ± 6.7 T/d. Observed heat flow was 70 ± 6.4 MW, compared with a pre-production value of 122 MW. This 52 MW reduction in surface heat flow is due to production-induced drying up of all alkali–Cl outflows (61.5 MW) and steam-heated pools (8.6 MW) within the Ohaaki West thermal area (OHW). The drying up of all alkali–Cl outflows at Ohaaki means that the soil zone is now the major natural pathway of heat release from the high-temperature reservoir. On the other hand, a net gain in thermal ground heat flow of 18 MW (from 25 MW to 43.3 ± 5 MW) at OHW is associated with permeability increases resulting from surface unit fracturing by production-induced ground subsidence. The Ohaaki East (OHE) thermal area showed no change in distribution of shallow and deep soil temperature contours despite 20 a of production, with an observed heat flow of 26.7 ± 3 MW and a CO 2 emission rate of 39 ± 3 T/d. The negligible change in the thermal status of the OHE thermal area is attributed to the low permeability of the reservoir beneath this area, which has limited production (mass extraction) and sheltered the area from the pressure decline within the main reservoir. Chemistry suggests that although alkali–Cl outflows once contributed significantly to the natural surface heat flow (∼50%) they contributed little ( 2 emissions due to the loss of >99% of the original CO 2 content due to depressurisation and boiling as the fluids ascended to the surface. Consequently, the soil has persisted as the major (99%) pathway of CO 2 release to the atmosphere from the high temperature reservoir at Ohaaki. The CO 2 flux and heat flow surveys indicate that despite 20 a of production the variability in location, spatial extent and magnitude of CO 2 flux remains consistent with established geochemical and geophysical models of the Ohaaki Field. At both OHW and OHE carbon isotopic analyses of soil gas indicate a two-stage fractionation process for moderate-flux (>60 g m −2 d −1 ) sites; boiling during fluid ascent within the underlying reservoir and isotopic enrichment as CO 2 diffuses through porous media of the soil zone. For high-flux sites (>300 g m −2 d −1 ), the δ 13 CO 2 signature (−7.4 ± 0.3‰ OHW and −6.5 ± 0.6‰ OHE) is unaffected by near-surface (soil zone) fractionation processes and reflects the composition of the boiled magmatic CO 2 source for each respective upflow. Flux thresholds of −2 d −1 for purely diffusive gas transport, between 30 and 300 g m −2 d −1 for combined diffusive–advective transport, and ⩾300 g m −2 d −1 for purely advective gas transport at Ohaaki were assigned. δ 13 CO 2 values and cumulative probability plots of CO 2 flux data both identified a threshold of ∼15 g m −2 d −1 by which background (atmospheric and soil respired) CO 2 may be differentiated from hydrothermal CO 2 .

Temporal variability of soil gas composition in landfill covers

In order to assess the temporal variability of the conditions for the microbial oxidation of methane in landfill cover soils and their driving variables, gas composition at non-emissive and strongly emissive locations (hotspots) was monitored on a seasonal, daily and hourly time scale on an old, unlined landfill in northern Germany. Our study showed that the impact of the various environmental factors varied with the mode of gas transport and with the time scale considered. At non-emissive sites, governed by diffusive gas transport, soil gas composition was subject to a pronounced seasonal variation. A high extent of aeration, low methane concentrations and a high ratio of CO 2 to CH 4 were found across the entire depth of the soil cover during the warm and dry period, whereas in the cool and moist period aeration was less and landfill gas migrated further upward. Statistically, variation in soil gas composition was best explained by the variation in soil temperature. At locations dominated by advective gas transport and showing considerable emissions of methane, this pattern was far less pronounced with only little increase in the extent of aeration during drier periods. Here, the change of barometric pressure was found to impact soil gas composition. On a daily scale under constant conditions of temperature, gas transport at both types of locations was strongly impacted by the change in soil moisture. On an hourly scale, under constant conditions of temperature and moisture, gas migration was impacted most by the change in barometric pressure. It was shown that at diffusion-dominated sites complete methane oxidation was achieved even under adverse wintry conditions, whereas at hotspots, even under favorable dry and warm conditions, aerobic biological activity can be limited to the upper crust of the soil.

Characterizing the impact of diffusive and advective soil gas transport on the measurement and interpretation of the isotopic signal of soil respiration

By measuring the isotopic signature of soil respiration, we seek to learn the isotopic composition of the carbon respired in the soil ( δ 13C R-s) so that we may draw inferences about ecosystem processes. Requisite to this goal is the need to understand how δ 13C R-s is affected by both contributions of multiple carbon sources to respiration and fractionation due to soil gas transport. In this study, we measured potential isotopic sources to determine their contributions to δ 13C R-s and we performed a series of experiments to investigate the impact of soil gas transport on δ 13C R-s estimates. The objectives of these experiments were to: i) compare estimates of δ 13C R-s derived from aboveground and belowground techniques, ii) evaluate the roles of diffusion and advection in a forest soil on the estimates of δ 13C R-s, and iii) determine the contribution of new and old carbon sources to δ 13C R-s for a Douglas-fir stand in the Pacific Northwest during our measurement period. We found a maximum difference of −2.36‰ between estimates of δ 13C R-s based on aboveground vs. belowground measurements; the aboveground estimate was enriched relative to the belowground estimate. Soil gas transport during the experiment was primarily by diffusion and the average belowground estimate of δ 13C R-s was enriched by 3.8–4.0‰ with respect to the source estimates from steady-state transport models. The affect of natural fluctuations in advective soil gas transport was little to non-existent; however, an advection–diffusion model was more accurate than a model based solely on diffusion in predicting the isotopic samples near the soil surface. Thus, estimates made from belowground gas samples will improve with an increase in samples near the soil surface. We measured a −1‰ difference in δ 13C R-s as a result of an experiment where advection was induced, a value which may represent an upper limit in fractionation due to advective gas transport in forest ecosystems. We found that aboveground measurements of δ 13C R-s may be particularly susceptible to atmospheric incursion, which may produce estimates that are enriched in 13C. The partitioning results attributed 69–98% of soil respiration to a source with a highly depleted isotopic signature similar to that of water-soluble carbon from foliage measured at our site.

Coupling between geochemical reactions and multicomponent gas and solute transport in unsaturated media: A reactive transport modeling study

The two‐way coupling that exists between biogeochemical reactions and vadose zone transport processes, in particular gas phase transport, determines the composition of soil gas. To explore these feedback processes quantitatively, multicomponent gas diffusion and advection are implemented into an existing reactive transport model that includes a full suite of geochemical reactions. Multicomponent gas diffusion is described on the basis of the dusty gas model, which accounts for all relevant gas diffusion mechanisms. The simulation of gas attenuation in partially saturated landfill soil covers, methane production, and oxidation in aquifers contaminated by organic compounds (e.g., an oil spill site) and pyrite oxidation in mine tailings demonstrate that both diffusive and advective gas transport can be affected by geochemical reactions. Methane oxidation in landfill covers reduces the existing upward pressure gradient, thereby decreasing the contribution of advective methane emissions to the atmosphere and enhancing the net flux of atmospheric oxygen into the soil column. At an oil spill site, methane oxidation causes a reversal in the direction of gas advection, which results in advective transport toward the zone of oxidation both from the ground surface and the deeper zone of methane production. Both diffusion and advection contribute to supply atmospheric oxygen into the subsurface, and methane emissions to the atmosphere are averted. During pyrite oxidation in mine tailings, pressure reduction in the reaction zone drives advective gas flow into the sediment column, enhancing the oxidation process. In carbonate‐rich mine tailings, calcite dissolution releases carbon dioxide, which partly offsets the pressure reduction caused by O2 consumption.

Shallow Subsurface Characterization of Gas Transport in a Playa Wetland

Controls on vadose-zone gas transport beneath and adjacent to a southern High Plains ephemeral lake (playa) were investigated. Under dry conditions, vertical gas permeability and tortuosity were enhanced by cracks and root tubules in the upper 2.5 m. Below this depth, the horizontal components of both permeability and tortuosity tensors were dominant. Both atmospheric pumping and pneumatic tests were used to determine the gas permeability tensor; whereas gas tracer tests were used to estimate the tortuosity tensor. The field data were elevated in a multilayer numerical model. The results suggest that gas movement can be dominated by both advection and diffusion, with vertical movement dominating in the shallow zone under dry conditions. After a large precipitation event, vertical gas permeability was reduced and large pressure differentials (as high as 1.5 kPa) were produced, creating driving forces for advective gas transport.

Estimation of cleanup time for continuous and pulsed soil vapor extraction

Soil vapor extraction (SVE) is a standard remediation technique for removing volatile organic compounds from contaminated soils. An important aspect in planning of SVE systems is the estimation of the time frame in which remediation can be achieved. In this study, analytical solutions are developed which can be used on the one hand for estimation of the required cleanup time and on the other hand for gaining insight into the processes which take place during the soil remediation. The analytical solutions are based on the well-mixed reservoir model as well as on the pure advective gas transport model, and their reliability is assessed by comparing their results with those of numerical simulations. Concerning mass transfer between the water phase and the solid phase, equilibrium as well as rate-limited sorption are considered. For the case of rate-limited sorption two extraction strategies (continuous extraction and pulsed extraction) are investigated.