Magnetic resonance microscopy of biofouling induced scale dependent transport in porous media

Non-invasive magnetic resonance microscopy (MRM) methods are applied to study biofouling of a homogeneous model porous media. MRM of the biofilm biomass using magnetic relaxation weighting shows the heterogeneous nature of the spatial distribution of the biomass as a function of growth. Spatially resolved MRM velocity maps indicate a strong variation in the pore scale velocity as a function of biofilm growth. The hydrodynamic dispersion dynamics for flow through the porous media is quantitatively characterized using a pulsed gradient spin echo technique to measure the propagator of the motion. The propagator indicates a transition in transport dynamics from a Gaussian normal diffusion process following a normal advection diffusion equation to anomalous transport as a function of biofilm growth. Continuous time random walk models resulting in a time fractional advection diffusion equation are shown to model the transition from normal to anomalous transport in the context of a conceptual model for the biofouling. The initially homogeneous porous media is transformed into a more complex heterogeneous porous media by the biofilm growth.     

Does pore water velocity affect the reaction rates of adsorptive solute transport in soils? Demonstration with pore-scale modelling

Modelling adsorptive solute transport in soils needs a number of parameters to describe its reaction kinetics and the values of these parameters are usually determined from batch and displacement experiments. Some experimental results reveal that when describing the adsorption as first-order kinetics, its associated reaction rates are not constants but vary with pore water velocity. Explanation of this varies but an independent verification of each explanation is difficult because simultaneously measuring the spatiotemporal distributions of dissolved and adsorbed solutes in soils is formidable. Pore-scale modelling could play an important role to address this gap and has received increased attention over the past few years. This paper investigated the transport of adsorptive solute in a simple porous medium using pore-scale modelling. Fluid flow through the void space of the medium was assumed to be laminar and in saturated condition, and solute transport consisted of advection and molecular diffusion; the sorption and desorption occurring at the fluid–solid interface were modelled as linear first-order kinetics. Based on the simulated spatiotemporal distribution of dissolved and adsorbed solutes at pore scale, volumetric-average reaction kinetics at macroscopic scale and its associated reactive parameters were measured. Both homogeneous adsorption where the reaction rates at microscopic scale are constant, and heterogeneous adsorption where the reaction rates vary from site to site, were investigated. The results indicate that, in contrast to previously thought, the macroscopic reaction rates directly measured from the pore-scale simulations do not change with pore velocity under both homogeneous and heterogeneous adsorptions. In particular, we found that for the homogeneous adsorption, the macroscopic adsorption remains first-order kinetic and can be described by constant reaction rates, regardless of flow rate; whilst for the heterogeneous adsorption, the macroscopic adsorption kinetics continues not to be affected by flow rate but is no longer first-order kinetics that can be described by constant reaction rates. We discuss how these findings could help explain some contrary literature reports over the dependence of reaction rates on pore water velocity.     

Analysis of laboratory data on ultrasonic monitoring of permeability reduction due to biopolymer formation in unconsolidated granular media

We show how to estimate the fluid permeability changes due to accumulated biopolymer within the pore space of a granular material using laboratory measurements of overall permeability, together with various well‐known quantitative measures (e.g., porosity, specific surface area, and formation factor) of the granular medium microstructure. The main focus of the paper is on mutual validation of existing theory and a synthesis of new experimental results. We find that the theory and data are in good agreement within normal experimental uncertainties. We also establish quantitative empirical relationships between seismic and/or acoustic attenuation and overall permeability for these same systems.     

A multiple-point statistics algorithm for 3D pore space reconstruction from 2D images

Fluid flow behavior in a porous medium is a function of the geometry and topology of its pore space. The construction of a three dimensional pore space model of a porous medium is therefore an important first step in characterizing the medium and predicting its flow properties. A stochastic technique for reconstruction of the 3D pore structure of unstructured random porous media from a 2D thin section training image is presented. The proposed technique relies on successive 2D multiple point statistics simulations coupled to a multi-scale conditioning data extraction procedure. The Single Normal Equation Simulation Algorithm (SNESIM), originally developed as a tool for reproduction of long-range, curvilinear features of geological structures, serves as the simulation engine. Various validating criteria such as marginal distributions of pore and grain, directional variograms, multiple-point connectivity curves, single phase effective permeability and two phase relative permeability calculations are used to analyze the results. The method is tested on a sample of Berea sandstone for which a 3D micro-CT scanning image is available. The results confirm that the equi-probable 3D realizations obtained preserve the typical patterns of the pore space that exist in thin sections, reproduce the long-range connectivities, capture the characteristics of anisotropy in both horizontal and vertical directions and have single and two phase flow characteristics consistent with those of the measured 3D micro-CT image.     

A transport phase diagram for pore-level correlated porous media

Transport in porous media is often characterized by the advection–dispersion equation, with the dispersion coefficient as the most important parameter that links the hydrodynamics to the transport processes. Morphological properties of any porous medium, such as pore size distribution, network topology, and correlation length control transport. In this study we explore the impact of correlation length on transport regime using pore-network modelling. Earlier direct simulation studies of dispersion in carbonate and sandstone rocks showed larger dispersion compared to granular homogenous sandpacks. However, in these studies, isolation of the impact of correlation length on transport regime was not possible due to the fundamentally different pore morphologies and pore-size distributions. Against this limitation, we simulate advection–dispersion transport for a wide range of Péclet numbers in unstructured irregular networks with “different” correlation lengths but “identical” pore size distributions and pore morphologies. Our simulation results show an increase in the magnitudes of the estimated dispersion coefficients in correlated networks compared to uncorrelated ones in the advection-controlled regime. The range of the Péclet numbers which dictate mixed advection–diffusion regime considerably reduces in the correlated networks. The findings emphasize the critical role of correlation length which is depicted in a conceptual transport phase diagram and the importance of accounting for the micro-scale correlation lengths into predictive stochastic pore-scale modelling.     

Permeability-porosity Relationships in Rocks Subjected to Various Evolution Processes

— It is well known that there is no “universal” permeability-porosity relationship valid in all porous media. However, the evolution of permeability and porosity in rocks can be constrained provided that the processes changing the pore space are known. In this paper, we review observations of the relationship between permeability and porosity during rock evolution and interpret them in terms of creation/destruction of effectively and non-effectively conducting pore space. We focus on laboratory processes, namely, plastic compaction of aggregates, elastic-brittle deformation of granular rocks, dilatant and thermal microcracking of dense rocks, chemically driven processes, as a way to approach naturally occurring geological processes. In particular, the chemically driven processes and their corresponding evolution permeability-porosity relationships are discussed in relation to sedimentary rocks diagenesis.     

Seeing through porous media: An experimental study for unveiling interstitial flows

We describe a novel inexpensive method, utilizing particle image velocimetry (PIV) and refractive index‐matching (RIM) for visualizing and quantifying the flow field within bio‐amended porous media. To date, this technique has been limited to idealized particles, whose refractive index does not match that of fresh water and thus requires specialized and often toxic or hazardous fluids. Here, we use irregularly shaped grains made of hydrogel as the solid matrix and water as the fluid. The advantage of using water is that it provides, for the first time, the opportunity to study both hydraulic and biological processes, which typically occur in soils and streambeds. By using RIM coupled with PIV (RIM‐PIV), we measured the interstitial flow field within a cell packed with granular material consisting of hydrogel grains in a size range of 1–8 mm, both in the presence and in the absence of Sinorhizobium meliloti bacteria (strain Rm8530). We also performed experiments with fluorescent tracer (fluorescein) and fluorescent microbes (Shewanella GPF MR‐1) to test the capability of visualizing solute transport and microbial movements. Results showed that the RIM‐PIV can measure the flow field for both biofilm‐free and biofilm‐covered hydrogel grains. The fluorescent tracer injection showed the ability to visualize both physical (concave surfaces and eddies) and biological (biofilms) transient storage zones, whereas the fluorescent microbe treatment showed the ability to track microbial movements within fluids. We conclude that the proposed methodology is a promising tool to visualize and quantify biofilm attachment, growth, and detachment in a system closer to natural conditions than a 2D flow cell experiment.     

Laboratory study of microseismicity spreading due to pore pressure change

Microseismicity variations in space and time due to pore pressure changes in an inhomogeneous permeable sample are considered. The analysis is based on laboratory experiments performed to study the relation between acoustic emission (which corresponds to microseismic events in real scale) and pore pressure changes due to water injection into, or release from, a porous sample under load. The study showed not only some characteristics of microseismicity variations during non-stationary fluid flow but also the possibility to solve the inverse problem of estimating local permeability from variations of microseismic activity in a particular volume of the porous medium.     

The effect of microporosity on transport properties in porous media

Sizeable amounts of connected microporosity with various origins can have a profound effect on important petrophysical properties of a porous medium such as (absolute/relative) permeability and capillary pressure relationships. We construct pore-throat networks that incorporate both intergranular porosity and microporosity. The latter originates from two separate mechanisms: partial dissolution of grains and pore fillings (e.g. clay). We then use the reconstructed network models to estimate the medium flow properties. In this work, we develop unique network construction algorithms and simulate capillary pressure–saturation and relative permeability–saturation curves for cases with inhomogeneous distributions of pores and micropores. Furthermore, we provide a modeling framework for variable amounts of cement and connectivity of the intergranular porosity and quantifying the conditions under which microporosity dominates transport properties. In the extreme case of a disconnected inter-granular network due to cementation a range of saturations within which neither fluid phase is capable of flowing emerges. To our knowledge, this is the first flexible pore scale model, from first principles, to successfully approach this behavior observed in tight reservoirs.     

Biofilm growth in porous media: Experiments,computational modeling at the porescale,and upscaling

Biofilm growth changes many physical properties of porous media such as porosity, permeability and mass transport parameters. The growth depends on various environmental conditions, and in particular, on flow rates. Modeling the evolution of such properties is difficult both at the porescale where the phase morphology can be distinguished, as well as during upscaling to the corescale effective properties. Experimental data on biofilm growth is also limited because its collection can interfere with the growth, while imaging itself presents challenges.In this paper we combine insight from imaging, experiments, and numerical simulations and visualization. The experimental dataset is based on glass beads domain inoculated by biomass which is subjected to various flow conditions promoting the growth of biomass and the appearance of a biofilm phase. The domain is imaged and the imaging data is used directly by a computational model for flow and transport. The results of the computational flow model are upscaled to produce conductivities which compare well with the experimentally obtained hydraulic properties of the medium. The flow model is also coupled to a newly developed biomass–nutrient growth model, and the model reproduces morphologies qualitatively similar to those observed in the experiment.     

Stochastic modeling of the permeability of randomly generated porous media

Permeability of porous media in subsurface environments is subject to potentially large uncertainties due to the heterogeneity of natural systems. In this study, a first-order reliability method (FORM) is combined with a lattice Boltzmann method (LBM) to estimate the permeability of randomly generated porous media. The proposed procedure provides an increased ease of addressing complex pore structures by employing LBM to model fluid flow, while inheriting the computational efficiency from FORM. Macroscale-equivalent permeability can thus be estimated with significantly reduced computational efforts, while maintaining a connection to the complex microscale fluid dynamics within a pore structure environment. Implemented on several randomly generated porous media domains, the proposed method provides 13–120 times the efficiency compared to Monte Carlo methods.     

Back front of seismicity induced by non‐linear pore pressure diffusion

Borehole fluid injections are accompanied by microseismic activity not only during but also after termination of the fluid injection. Previously, this phenomenon has been analysed, assuming that the main triggering mechanism is governed by a linear pressure diffusion in a hydraulically isotropic medium. In this context the so‐called back front of seismicity has been introduced, which allows to characterize the hydraulic transport from the spatiotemporal distribution of post‐injection induced events. However, rocks are generally anisotropic, and in addition, fluid injections can strongly enhance permeability. In this case, permeability becomes a function of pressure. For such situations, we carry out a comprehensive study about the behaviour and parametrization of the back front. Based on a model of a factorized anisotropic pressure dependence of permeability, we present an approach to reconstruct the principal components of the diffusivity tensor. We apply this approach to real microseismic data and show that the back front characterizes the least hydraulic transport. To investigate the back front of non‐linear pore‐fluid pressure diffusion, we numerically consider a power‐law and an exponential‐dependent diffusivity. To account for a post‐injection enhanced hydraulic state of the rock, we introduce a model of a frozen (i.e., nearly unchanged after the stimulation) medium diffusivity and generate synthetic seismicity. We find that, for a weak non‐linearity and 3D exponential diffusion, the linear diffusion back front is still applicable. This finding is in agreement with microseismic data from Ogachi and Fenton Hill. However, for a strong non‐linear fluid–rock interaction such as hydraulic fracturing, the back front can significantly deviate from a time dependence of a linear diffusion back front. This is demonstrated for a data set from the Horn River Basin. Hence, the behaviour of the back front is a strong indicator of a non‐linear fluid–rock interaction.     

Application of GA in optimization of pore network models generated by multi-cellular growth algorithms

In pore network modeling, the void space of a rock sample is represented at the microscopic scale by a network of pores connected by throats. Construction of a reasonable representation of the geometry and topology of the pore space will lead to a reliable prediction of the properties of porous media. Recently, the theory of multi-cellular growth (or L-systems) has been used as a flexible tool for generation of pore network models which do not require any special information such as 2D SEM or 3D pore space images. In general, the networks generated by this method are irregular pore network models which are inherently closer to the complicated nature of the porous media rather than regular lattice networks. In this approach, the construction process is controlled only by the production rules that govern the development process of the network. In this study, genetic algorithm has been used to obtain the optimum values of the uncertain parameters of these production rules to build an appropriate irregular lattice network capable of the prediction of both static and hydraulic information of the target porous medium.     

裂缝-孔隙介质储层渗透率表征及其影响因素分析

含裂缝多孔介质渗透率预测是非常规油气资源勘探开发的一个紧迫问题.现有多孔介质岩石物理模型通常利用圆形孔管模拟宏观岩石孔隙空间,难以定量描述软孔隙/裂缝在压力作用下的闭合情况,缺乏裂缝/孔隙间流量交换的连通机制.本文提出含三维裂缝/软孔隙网络多孔介质模型,将储层岩石裂缝/软孔隙表示为椭圆截面微管,建立了周期性压力作用下微观裂缝流量表达式,通过网络模型和流量守恒条件,得到含有三维裂缝/软孔隙网络的多孔介质渗透率计算方法.数值算例表明,预测结果与实验数据分布范围吻合很好,能够给出不同类型岩心对应孔隙纵横比的分布图.三维裂缝/软孔隙网络模型建立了宏观可观测量与裂缝参数之间关系,能够定量分析岩石渗透率随裂缝体密度、纵横比、孔隙流体类型和围压等因素的变化规律,为复杂条件下储层渗透率预测提供了一种有效方法.

     

Permeability changes in layered sediments: impact of particle release

One of the mechanisms of sudden particle release from grain surfaces in natural porous media is a decrease in salt concentration of the permeating fluid to below the critical salt concentration. Particle release can cause a change in hydraulic conductivity of the matrix, either by washing out the fines and thus increasing the pore sizes or by the plugging of pore constrictions. The phenomenon of permeability changes as a result of particle detachment was investigated in a series of column experiments. Coarse and fine sediments from the Hanford Formation in southeast Washington were tested. Columns were subject to a pulse of highly saline solution (NaNO3) followed by a fresh water shock causing particle release. Outflow rates and changes in hydraulic head as well as electric conductivity and pH were monitored over time. No permeability decrease occurred within the coarse matrix alone. However, when a thin layer of fine sediment was embedded within the coarse material (mimicking field conditions at the Hanford site), permeability irreversibly decreased to 10% to 20% of the initial value. Evidence suggests that most of this permeability decrease was a result of particles detached within the fine layer and its subsequent clogging. An additional observation was a sudden increase in pH in the outflow solution, generated in situ during the fresh water shock. Because layered systems are common in natural settings, our results suggest that alteration between sodium solution and fresh water can lead to particle release and subsequently reduce the overall permeability of the matrix.     

The Mechanical Coupling of Fluid-Filled Granular Material Under Shear

The coupled mechanics of fluid-filled granular media controls the physics of many Earth systems, for example saturated soils, fault gouge, and landslide shear zones. It is well established that when the pore fluid pressure rises, the shear resistance of fluid-filled granular systems decreases, and, as a result, catastrophic events such as soil liquefaction, earthquakes, and accelerating landslides may be triggered. Alternatively, when the pore pressure drops, the shear resistance of these geosystems increases. Despite the great importance of the coupled mechanics of grain–fluid systems, the basic physics that controls this coupling is far from understood. Fundamental questions that must be addressed include: what are the processes that control pore fluid pressurization and depressurization in response to deformation of the granular skeleton? and how do variations of pore pressure affect the mechanical strength of the grains skeleton? To answer these questions, a formulation for the pore fluid pressure and flow has been developed from mass and momentum conservation, and is coupled with a granular dynamics algorithm that solves the grain dynamics, to form a fully coupled model. The pore fluid formulation reveals that the evolution of pore pressure obeys viscoelastic rheology in response to pore space variations. Under undrained conditions elastic-like behavior dominates and leads to a linear relationship between pore pressure and overall volumetric strain. Viscous-like behavior dominates under well-drained conditions and leads to a linear relationship between pore pressure and volumetric strain rate. Numerical simulations reveal the possibility of liquefaction under drained and initially over-compacted conditions, which were often believed to be resistant to liquefaction. Under such conditions liquefaction occurs during short compactive phases that punctuate the overall dilative trend. In addition, the previously recognized generation of elevated pore pressure under undrained compactive conditions is observed. Simulations also show that during liquefaction events stress chains are detached, the external load becomes completely supported by the pressurized pore fluid, and shear resistance vanishes.     

Efficient flow and transport simulations in reconstructed 3D pore geometries

Upscaling pore-scale processes into macroscopic quantities such as hydrodynamic dispersion is still not a straightforward matter for porous media with complex pore space geometries. Recently it has become possible to obtain very realistic 3D geometries for the pore system of real rocks using either numerical reconstruction or micro-CT measurements. In this work, we present a finite element–finite volume simulation method for modeling single-phase fluid flow and solute transport in experimentally obtained 3D pore geometries. Algebraic multigrid techniques and parallelization allow us to solve the Stokes and advection–diffusion equations on large meshes with several millions of elements. We apply this method in a proof-of-concept study of a digitized Fontainebleau sandstone sample. We use the calculated velocity to simulate pore-scale solute transport and diffusion. From this, we are able to calculate the a priori emergent macroscopic hydrodynamic dispersion coefficient of the porous medium for a given molecular diffusion D_m of the solute species. By performing this calculation at a range of flow rates, we can correctly predict all of the observed flow regimes from diffusion dominated to convection dominated.     

Prediction of diffusion coefficients in porous media using tortuosity factors based on interfacial areas

Determination of aqueous phase diffusion coefficients of solutes through porous media is essential for understanding and modeling contaminant transport. Prediction of diffusion coefficients in both saturated and unsaturated zones requires knowledge of tortuosity and constrictivity factors. No methods are available for the direct measurement of these factors, which are empirical in their definition. In this paper, a new definition for the tortuosity factor is proposed, as the real to ideal interfacial area ratio. We define the tortuosity factor for saturated porous media (tau5) as the ratio S/S(o) (specific surface of real porous medium to that of an idealized capillary bundle). For unsaturated media, tortuosity factor (tau(a)) is defined as a(aw)/a(aw),o (ratio of the specific air-water interfacial area of real and the corresponding idealized porous medium). This tortuosity factor is suitably measured using sorptive tracers (e.g., nitrogen adsorption method) for saturated media and interfacial tracers for unsaturated media. A model based on this new definition of tortuosity factors, termed the interfacial area ratio (IAR) model, is presented for the prediction of diffusion coefficients as a function of the degree of water saturation. Diffusion coefficients and diffusive resistances measured in a number of saturated and unsaturated granular porous media, for solutes in dilute aqueous solutions, agree well with the predictions of the IAR model. A comparison of permeability of saturated sands estimated based on tau(s) and the same based on the Kozeny-Carman equation confirm the usefulness of the tau(s) parameter as a measure of tortuosity.     

Visualization experiments of biodegradation in porous media and calculation of the biodegradation rate

Biodegradation in porous media is studied with carefully controlled and well-characterized experiments in model porous media constructed of etched glass. Porous media of this type allow visual observation of the phenomena that take place at pore scale. An aqueous solution of five organic pollutants (toluene, phenol, o-cresol, naphthalene and 1,2,3-trimethylbenzene) was used as a model NAPL (representing creosote). The bacteria used were Pseudomonas fluorescens, which are indigenous (even predominant) in many contaminated soils. The maximum aqueous concentrations of the specific organic substances, below which biodegradation becomes possible, were determined as a function of temperature from toxicity experiments. Visualization experiments were made under various flow velocities and organic loadings to study the morphology and thickness of the biofilm as a function of the pore size and the distance from the entrance, and the efficiency of biodegradation. The efficiency of biodegradation decreased as the aqueous concentration of NAPL at the inlet increased and/or as the flow velocity increased. The thickness of biofilm decreased as the distance from the inlet increased and/or the pore diameter decreased. A quasi-steady-state theoretical model of biodegradation was used to calculate the values of the mesoscopic biochemical rates and to predict the profile of NAPL concentration in the porous medium and the thickness of biofilm in pores. The agreement between experimental data and model predictions is quite satisfactory.     

Flow and entrapment of dense nonaqueous phase liquids in physically and chemically heterogeneous aquifer formations

The migration and entrapment of dense nonaqueous phase liquids (DNAPLs) in aquifer formations is typically believed to be controlled by physical heterogeneities. This belief is based upon the assumption that permeability and capillary properties are determined by the soil texture. Capillarity and relative permeability, however, will also depend on porous medium wettability characteristics. This wettability may vary spatially in a formation due to variations in aqueous phase chemistry, contaminant aging, and/or variations in mineralogy and organic matter distributions. In this work, a two-dimensional multiphase flow simulator is modified to simulate coupled physical and chemical formation heterogeneity. To model physical heterogeneity, a spatially correlated permeability field is generated, and then related to the capillary pressure-saturation function according to Leverett scaling. Spatial variability of porous medium wettability is assumed to be correlated with the natural logarithm of the intrinsic permeability. The influence of wettability on the hysteretic hydraulic property relations is also modeled. The simulator is then employed to investigate the potential influence of coupled physical and chemical heterogeneity on DNAPL flow and entrapment. For reasonable ranges of wettability characteristics, simulations demonstrate that spatial variations in wettability can have a dramatic impact on DNAPL distributions. Higher organic saturations, increased lateral spreading, and decreased depth of infiltration were predicted when the contact angle was varied spatially. When chemical heterogeneity was defined by spatial variation of organic-wet solid fractions (fractional wettability porous media), however, the resultant organic saturation distributions were more similar to those for perfectly water-wet media, due to saturation dependent wettability effects on the hydraulic property relations.