scholarly journals MODELING OF MULTISCALE POROUS MEDIA

2011 ◽  
Vol 27 (1) ◽  
pp. 23 ◽  
Author(s):  
Bibhu Biswal ◽  
Pål-Eric Øren ◽  
Rudolf J Held ◽  
Stig Bakke ◽  
Rudolf Hilfer

A stochastic geometrical modeling method for reconstructing three dimensional pore scale microstructures of multiscale porous media is presented. In this method the porous medium is represented by a random but spatially correlated structure of objects placed in the continuum. The model exhibits correlations with the sedimentary textures, scale dependent intergranular porosity over many decades, vuggy or dissolution porosity, a percolating pore space, a fully connected matrix space, strong resolution dependence and wide variability in the permeabilities and other properties. The continuum representation allows discretization at arbitrary resolutions providing synthetic micro-computertomographic images for resolution dependent fluid flow simulation. Model implementations for two different carbonate rocks are presented. The method can be used to generate pore scale models of a wide class of multiscale porous media.

Georesursy ◽  
2020 ◽  
Vol 22 (1) ◽  
pp. 4-12
Author(s):  
Timur R. Zakirov ◽  
Maxim G. Khramchenkov

This paper presents the results of numerical simulations of two-phase flows in porous media under capillary forces dominance. For modeling of immiscible two-phase flow, the lattice Boltzmann equations with multi relaxation time operator were applied, and the interface phenomena was described with the color-gradient method. The objective of study is to establish direct links between quantitative characteristics of the flow and invasion events, using high temporal resolution when detecting simulation results. This is one of the few works where Haines jumps (rapid invasion event which occurs at meniscus displacing from narrow pore throat to its wide body) are considered in three-dimensional natural pore space, but the focus is also on the displacement mechanics after jumps. It was revealed the sequence of pore scale events which can be considered as a period of drainage process: rapid invasion event during Haines jump; finish of jump and continuation of uniform invasion in current pore; switching of mobile interfaces and displacement in new region. The detected interface change, along with Haines jump, is another distinctive feature of the capillary forces action. The change of the mobile interfaces is manifested in step-like behavior of the front movement. It was obtained that statistical distributions of pressure drops during Haines jumps obey lognormal law. When investigating the flow rate and surface tension effect on the pressure drop statistics it was revealed that these parameters practically don’t affect on the statistical distribution and influence only on the magnitude of pressure drops and the number of individual Haines jumps.


2014 ◽  
Vol 136 (4) ◽  
Author(s):  
Zhenglun Alan Wei ◽  
Zhongquan Charlie Zheng ◽  
Xiaofan Yang

A parallel implementation of an immersed-boundary (IB) method is presented for low Reynolds number flow simulations in a representative elementary volume (REV) of porous media that are composed of a periodic array of regularly arranged structures. The material of the structure in the REV can be solid (impermeable) or microporous (permeable). Flows both outside and inside the microporous media are computed simultaneously by using an IB method to solve a combination of the Navier–Stokes equation (outside the microporous medium) and the Zwikker–Kosten equation (inside the microporous medium). The numerical simulation is firstly validated using flow through the REVs of impermeable structures, including square rods, circular rods, cubes, and spheres. The resultant pressure gradient over the REVs is compared with analytical solutions of the Ergun equation or Darcy–Forchheimer law. The good agreements demonstrate the validity of the numerical method to simulate the macroscopic flow behavior in porous media. In addition, with the assistance of a scientific parallel computational library, PETSc, good parallel performances are achieved. Finally, the IB method is extended to simulate species transport by coupling with the REV flow simulation. The species sorption behaviors in an REV with impermeable/solid and permeable/microporous materials are then studied.


2021 ◽  
Author(s):  
Pietro de Anna ◽  
Amir A. Pahlavan ◽  
Yutaka Yawata ◽  
Roman Stocker ◽  
Ruben Juanes

<div> <div> <div> <p>Natural soils are host to a high density and diversity of microorganisms, and even deep-earth porous rocks provide a habitat for active microbial communities. In these environ- ments, microbial transport by disordered flows is relevant for a broad range of natural and engineered processes, from biochemical cycling to remineralization and bioremediation. Yet, how bacteria are transported and distributed in the sub- surface as a result of the disordered flow and the associ- ated chemical gradients characteristic of porous media has remained poorly understood, in part because studies have so far focused on steady, macroscale chemical gradients. Here, we use a microfluidic model system that captures flow disorder and chemical gradients at the pore scale to quantify the transport and dispersion of the soil-dwelling bacterium Bacillus subtilis in porous media. We observe that chemotaxis strongly modulates the persistence of bacteria in low-flow regions of the pore space, resulting in a 100% increase in their dispersion coefficient. This effect stems directly from the strong pore-scale gradients created by flow disorder and demonstrates that the microscale interplay between bacterial behaviour and pore-scale disorder can impact the macroscale dynamics of biota in the subsurface.</p> </div> </div> </div>


2021 ◽  
Author(s):  
Marco Dentz ◽  
Alexandre Puyguiraud ◽  
Philippe Gouze

<p>Transport of dissolved substances through porous media is determined by the complexity of the pore space and diffusive mass transfer within and between pores. The interplay of diffusive pore-scale mixing and spatial flow variability are key for the understanding of transport and reaction phenomena in porous media. We study the interplay of pore-scale mixing and network-scale advection through heterogeneous porous media, and its role for the evolution and asymptotic behavior of hydrodynamic dispersion. In a Lagrangian framework, we identify three fundamental mechanisms of pore-scale mixing that determine large scale particle motion: (i) The smoothing of intra-pore velocity contrasts, (ii) the increase of the tortuosity of particle paths, and (iii) the setting of a maximum time for particle transitions. Based on these mechanisms, we derive an upscaled approach that predicts anomalous and normal hydrodynamic dispersion based on the characteristic pore length, Eulerian velocity distribution and Péclet number. The theoretical developments are supported and validated by direct numerical flow and transport simulations in a three-dimensional digitized Berea sandstone sample obtained using X-Ray microtomography. Solute breakthrough curves, are characterized by an intermediate power-law behavior and exponential cut-off, which reflect pore-scale velocity variability and intra-pore solute mixing. Similarly, dispersion evolves from molecular diffusion at early times to asymptotic hydrodynamics dispersion via an intermediate superdiffusive regime. The theory captures the full evolution form anomalous to normal transport behavior at different Péclet numbers as well as the Péclet-dependence of asymptotic dispersion. It sheds light on hydrodynamic dispersion behaviors as a consequence of the interaction between pore-scale mixing and Eulerian flow variability. </p>


2019 ◽  
Vol 51 (1) ◽  
pp. 429-449 ◽  
Author(s):  
Kamaljit Singh ◽  
Michael Jung ◽  
Martin Brinkmann ◽  
Ralf Seemann

Liquid invasion into a porous medium is a phenomenon of great importance in both nature and technology. Despite its enormous importance, there is a surprisingly sparse understanding of the processes occurring on the scale of individual pores and of how these processes determine the global invasion pattern. In particular, the exact influence of the wettability remains unclear besides the limiting cases of very small or very large contact angles of the invading fluid. Most quantitative pore-scale experiments and theoretical considerations have been conducted in effectively two-dimensional (2D) micromodels and Hele–Shaw geometries. Although these pioneering works helped to unravel some of the physical aspects of the displacement processes, the relevance of 2D models has not always been appreciated for natural porous media. With the availability of X-ray microtomography, 3D imaging has become a standard for exploring pore-scale processes in porous media. Applying advanced postprocessing routines and synchrotron microtomography, researchers can image even slow flow processes in real time and extract relevant material parameters like the contact angle from the interfaces in the pore space. These advances are expected to boost both theoretical and experimental understanding of pore-scale processes in natural porous media.


2017 ◽  
Vol 5 (1) ◽  
pp. SB57-SB67 ◽  
Author(s):  
Nattavadee Srisutthiyakorn ◽  
Gerald M. Mavko

Hydraulic tortuosity is an important parameter in characterizing fluid-flow heterogeneity in porous media. The most basic definition of tortuosity is the ratio of the average flow path length to the sample length. Although this definition seems straightforward, the lack of understanding and the lack of proper ways to measure tortuosity make it one of the most abused parameters in rock physics. Hydraulic tortuosity is often treated merely as a fitting factor, or worse, it is neglected by being combined with a geometric factor in the Kozeny-Carman (KC) equation. Often, the tortuosity is obtained from laboratory measurements of porosity, permeability, and specific surface area by inverting the KC equation. This approach has a major pitfall because it treats tortuosity as a fitting factor, and the inverted tortuosity is often unphysically high. In contrast, we obtained the tortuosity from 3D segmented binary images of porous media using streamlines extracted from a local flux, the output from the lattice Boltzmann method (LBM) flow simulation. After obtaining streamlines from each sample, we calculated the distribution of tortuosities and flux-weighted average tortuosity. With the tortuosity measurement from streamlines, every parameter in the KC equation can be measured accurately from 3D segmented binary images. We found, however, that the KC equation is still missing some important geometric information needed to predict permeability. With known parameters and without a fitting factor, the KC equation predicts permeability higher by one to two orders of magnitude than that predicted by the LBM. We searched for a missing parameter by exploring various concepts such as connected pore space and pore throat distribution. We found that the connected pore space does not contribute to the difference between the KC permeability and LBM permeability, whereas, as we learn with sinusoidal pipe examples, the pore throat distribution captures what is missing from the KC equation.


1998 ◽  
Vol 09 (06) ◽  
pp. 837-849 ◽  
Author(s):  
A. M. Vidales ◽  
J. L. Riccardo ◽  
G. Zgrablich

Immiscible displacement at pore level on a three-dimensional correlated porous network is simulated allowing flow of the wetting phase along crevices of the pore walls (possibility of snap-off in throats) and advance through the centers of the pore space with different pore and throat filling conditions, leading to a cooperative filling. When these two mechanisms compete, different patterns arise. We study the effect of the correlation strength on the onset of each pattern. We do not take buoyancy forces into account.


Sign in / Sign up

Export Citation Format

Share Document