skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: A (Simplified) Biogeochemical Numerical Model to Predict Saturation, Porosity and Permeability During Microbially Induced Desaturation and Precipitation
Abstract Microbially Induced Desaturation and Precipitation (MIDP) through denitrification is an emerging ground improvement method in which indigenous nitrate reducing bacteria are stimulated to introduce biogas, biominerals and biomass in the soil matrix. In this study, a numerical model is developed to evaluate the effect of biogas, biominerals and biomass on the hydraulic properties of soils treated with MIDP. The proposed model couples the biochemical conversions to changes of porosity and water saturation and predicts changes in permeability through two separate power law equations. Experimental studies from the literature are used to calibrate the model. Comparing the results with other studies on bioclogging or biomineralization in porous media reveals that the combined production of biogas, biomass, and biominerals results in efficient clogging, in the sense that only a small amount of products leads to a substantial permeability reduction. Based on this comparison, the authors postulate that biogenic gas bubbles preferably form within the larger pore bodies. The presence of biogenic gas in the larger pore bodies forces calcium carbonate minerals and biomass to be formed mainly at the pore throats. The interaction between the different phases results in more efficient clogging than observed in other studies which focus on a single product only.  more » « less
Award ID(s):
2349122
PAR ID:
10493574
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
AGU
Date Published:
Journal Name:
Water Resources Research
Volume:
59
Issue:
1
ISSN:
0043-1397
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Microbially induced desaturation and precipitation (MIDP) via denitrification has the potential to reduce earthquake-induced liquefaction potential by two mechanisms: calcium carbonate precipitation to mechanically strengthen soil and biogenic gas production to desaturate and dampen pore pressure changes in soil. Lab-scale tests have demonstrated effective desaturation and improved mechanical strength by MIDP. However, in laboratory tests, gas pockets and lenses form causing upheaval as a result of low overburden pressures. The characteristics of biogenic gas formation, distribution, and retention need to be evaluated to gain comprehensive understanding of the effectiveness of this treatment at depth before and after an earthquake event. MIDP treatment during centrifuge loading conditions is being performed to simulate field stress conditions, prior to complete process scale-up for field application. A simplified numerical model was developed to evaluate the scaling effects on biogenic gas generation between the centrifuge model and prototype scale. The results indicate that diffusion of soluble N2 is negligible at both the model and prototype scales for the simulated reaction rate. However, the simplified model did not consider other pore-scale influences and mixing from liquid-gas transfer and transport. Future modeling work will need to add these features. 
    more » « less
  2. Hagfishes defend themselves from gill-breathing predators by producing large volumes of fibrous slime when attacked. The slime's effectiveness comes from its ability to clog predators' gills, but the mechanisms by which hagfish slime clogs are uncertain, especially given its remarkably dilute concentration of solids. We quantified the clogging performance of hagfish slime over a range of concentrations, measured the contributions of its mucous and thread components, and measured the effect of turbulent mixing on clogging. To assess the porous structure of hagfish slime, we used a custom device to measure its Darcy permeability. We show that hagfish slime clogs at extremely dilute concentrations like those found in native hagfish slime and displays clogging performance that is superior to three thickening agents. We report an extremely low Darcy permeability for hagfish slime, and an effective pore size of 10–300 nm. We also show that the mucous and thread components play distinct yet crucial roles, with mucus being responsible for effective clogging and low permeability and the threads imparting mechanical strength and retaining clogging function over time. Our results provide new insights into the mechanisms by which hagfish slime clogs gills and may inspire the development of ultra-soft materials with novel properties. 
    more » « less
  3. Abstract Geophysical fluid‐granular flows, such as pyroclastic currents and debris flows, owe much of their runout and hazard behavior to the occurrence and time‐variant decay of a flow‐internal fluid pore pressure. However, modeling the effects of fluid pore pressure to forecast hazards is challenging because a unified method in Earth Sciences to quantitatively determine the permeability of these natural mixtures is currently missing. Here we combine experiments on fluidization and defluidization of pyroclastic materials, eolian sediments, and glass beads mixtures with numerical multiphase simulations to compare previous attempts to compute the permeability of complex natural particle‐fluid mixtures. In analogy to particle‐engineering studies on simple gas‐particle mixtures, we demonstrate that the effective length‐scale in the characterization of the fluid‐particle interaction of complex natural mixtures is the product of the Sauter mean diameter and the particle sphericity. Its use in the Kozeny‐Carman equation allows accurate prediction of mixture permeability, and we suggest the routine calculation of the Sauter mean from grain size distributions of the deposits of geophysical mass flows in Earth Sciences. We also show, through defluidization experiments, that the duration of gas retention in natural mixtures is well described when using the Sauter mean as the effective particle size. Further, we show through multiphase simulations that initial bed expansion extends the pore pressure diffusion timescale up to nine times. These results can be applied to small‐to‐large volume dense pyroclastic currents where the ranges of Sauter mean diameter predict gas retention for long duration and to debris flows and snow avalanches. 
    more » « less
  4. Pore-scale modeling is essential in understanding and predicting flow and transport properties of rocks. Generally, pore-scale modeling is dependent on imaging technologies such as Micro Computed Tomography (micro-CT), which provides visual confirmation into the pore microstructures of rocks at a representative scale. However, this technique is limited in the ability to provide high resolution images showing the pore-throats connecting pore bodies. Pore scale simulations of flow and transport properties of rocks are generally done on a single 3D pore microstructure image. As such, the simulated properties are only representative of the simulated pore-scale rock volume. These are the technological and computational limitations which we address here by using a stochastic pore-scale simulation approach. This approach consists of constructing hundreds of 3D pore microstructures of the same pore size distribution and overall porosity but different pore connectivity. The construction of the 3D pore microstructures incorporates the use of Mercury Injection Capillary Pressure (MICP) data to account for pore throat size distribution, and micro-CT images to account for pore body size distribution. The approach requires a small micro-CT image volume (7–19 mm3) to reveal key pore microstructural features that control flow and transport properties of highly heterogeneous rocks at the core-scale. Four carbonate rock samples were used to test the proposed approach. Permeability calculations from the introduced approach were validated by comparing laboratory measured permeability of rock cores and permeability estimated using five well-known core-scale empirical model equations. The results show that accounting for the stochastic connectivity of pores results in a probabilistic distribution of flow properties which can be used to upscale pore-scale simulated flow properties to the core-scale. The use of the introduced stochastic pore-scale simulation approach is more beneficial when there is a higher degree of heterogeneity in pore size distribution. This is shown to be the case with permeability and hydraulic tortuosity which are key controls of flow and transport processes in rocks. 
    more » « less
  5. Abstract Fluids influence fault zone strength and the occurrence of earthquakes, slow slip events, and aseismic slip. We introduce an earthquake sequence model with fault zone fluid transport, accounting for elastic, viscous, and plastic porosity evolution, with permeability having a power‐law dependence on porosity. Fluids, sourced at a constant rate below the seismogenic zone, ascend along the fault. While the modeling is done for a vertical strike‐slip fault with 2D antiplane shear deformation, the general behavior and processes are anticipated to apply also to subduction zones. The model produces large earthquakes in the seismogenic zone, whose recurrence interval is controlled in part by compaction‐driven pressurization and weakening. The model also produces a complex sequence of slow slip events (SSEs) beneath the seismogenic zone. The SSEs are initiated by compaction‐driven pressurization and weakening and stalled by dilatant suctions. Modeled SSE sequences include long‐term events lasting from a few months to years and very rapid short‐term events lasting for only a few days; slip is ∼1–10 cm. Despite ∼1–10 MPa pore pressure changes, porosity and permeability changes are small and hence fluid flux is relatively constant except in the immediate vicinity of slip fronts. This contrasts with alternative fault valving models that feature much larger changes in permeability from the evolution of pore connectivity. Our model demonstrates the important role that compaction and dilatancy have on fluid pressure and fault slip, with possible relevance to slow slip events in subduction zones and elsewhere. 
    more » « less