Lattice Boltzmann modeling and evaluation of fluid flow in heterogeneous porous media involving multiple matrix constituents

Gao, Jing; Xing, Huilin; Tian, Zhilong; Mühlhaus, H. B.

doi:10.1016/j.cageo.2013.07.019

Cited by 41 publications

(17 citation statements)

References 22 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…In the simulation of gas flow through a porous system, the GLBM model can be employed by simply replacing the usual computational lattices with porous lattices in the domain occupied by the porous medium in both 2-D and 3-D situations [38]. In this study, the porosity remains the same in each component, and each porous lattice in the simulation domain is characterized by a pore radius based on the pore size distribution in OM and IOM.…”

Section: Generalized Lattice Boltzmann Modelmentioning

confidence: 99%

Apparent permeability prediction of organic shale with generalized lattice Boltzmann model considering surface diffusion effect

Wang

Chen

Kang

et al. 2016

Fuel

View full text Add to dashboard Cite

a b s t r a c tGas flow in shale is associated with both organic matter (OM) and inorganic matter (IOM) which contain nano-pores ranging in size from a few to hundreds of nano-meters. In addition to the non-continuum effect which leads to an apparent permeability of gas higher than the intrinsic permeability, the surface diffusion of adsorbed gas in organic pores also can influence the apparent permeability through its own transport mechanism. In this study, a generalized lattice Boltzmann model (GLBM) is employed for gas flow through the reconstructed shale matrix consisting of OM and IOM. The ExpectationMaximization (EM) algorithm is used to assign the pore size distribution to each component, and the dusty gas model (DGM) and generalized Maxwell-Stefan model (GMS) are adopted to calculate the apparent permeability accounting for multiple transport mechanisms including viscous flow, Knudsen diffusion and surface diffusion. Effects of pore radius and pressure on permeability of both IOM and OM as well as effects of Langmuir parameters on OM are investigated. The effect of total organic content and distribution on the apparent permeability of the reconstructed shale matrix at different surface diffusivity is also studied. It is found that the influence of pore size and pressure on the apparent permeability of organic matter is affected by the surface diffusion of adsorbed gas. Moreover, surface diffusion plays a significant role in determining apparent permeability and the velocity distribution of shale matrix.

show abstract

Section: Generalized Lattice Boltzmann Modelmentioning

confidence: 99%

Apparent permeability prediction of organic shale with generalized lattice Boltzmann model considering surface diffusion effect

Wang

Chen

Kang

et al. 2016

Fuel

View full text Add to dashboard Cite

show abstract

“…By simulating streaming and collision processes across a limited number of particles, the intrinsic particle interactions evince a microcosm of viscous flow behaviour applicable across the greater mass (e.g., Sukop and Thorne, 2007). Due to its particulate nature and local dynamics, LBM has several advantages over other conventional CFD methods, especially in dealing with complex boundaries such as multiphase fluid flow in complicated fractures/porous media with detailed microstructures, incorporating microscopic interactions as well as chemical reactions, and parallelization of the algorithm (Tian et al, 2014;Gao et al, 2013;Gao and Xing 2012a, b). (5) ESyS_Crustal is a finite element method based module developed for the interacting fault system simulation.…”

Section: An Integrated Geothermal Reservoir Simulatormentioning

confidence: 99%

Recent development in numerical simulation of enhanced geothermal reservoirs

et al. 2015

Self Cite

View full text Add to dashboard Cite

This paper briefly introduces the current state in computer modelling of geothermal reservoir system and then focuses on our research efforts in high performance simulation of enhanced geothermal reservoir system. A novel supercomputer simulation tool has been developing towards simulating the highly non-linear coupled geomechanical-fluid flow-thermal systems involving heterogeneously fractured geomaterials at different spatial and temporal scales. It is applied here to simulate and visualise the enhanced geothermal system (EGS), such as (1) visualisation of the microseismic events to monitor and determine where/how the underground rupture proceeds during a hydraulic stimulation, to generate the mesh using the recorded data for determining the domain of the ruptured zone and to evaluate the material parameters (i.e., the permeability) for the further numerical analysis and evaluation of the enhanced geothermal reservoir; (2) converting the available fractured rock image/fracture data as well as the reservoir geological geometry to suitable meshes/grids and further simulating the fluid flow in the complicated fractures involving the detailed description of fracture dimension and geometry by the lattice Boltzmann method and/or finite element method; (3) interacting fault system simulation to determine the relevant complicated rupture process for evaluating the geological setting and the in-situ reservoir properties; (4) coupled thermo-fluid flow analysis of a geothermal reservoir system for an optimised geothermal reservoir design and management. A few of application examples are presented to show its usefulness in simulating the enhanced geothermal reservoir system. KEY WORDS: numerical simulation, geothermal, EGS, microseismicity, finite element method, lattice Boltzmann method. INTRODUCTIONOver the past 30 years, a large amount of research and field testing on engineered geothermal system (EGS), also known as hot dry rock (HDR) and hot fractured rock (HFR) geothermal reservoirs, have been accomplished worldwide which include the reservoir construction, fluid circulation and heat extraction (e.g., Tester et al., 2006). A successful EGS reservoir depends on thermal-fluid flow at any given time which primarily determined by its mean temperature and pressure, the nature of the interconnected network of hydraulic stimulated joints and open fractures (including both stimulated and natural), the cumulative amount of fluid circulation (reservoir cooling) that has occurred and water loss (Rybach, 2010; Brown et al., 1999). The reservoir characteristics are complicated and functions of the applied reservoir pressure/stress that are controlling the nature and degree of interconnection within the network of fractures, therefore it is crucial to have good measures and understanding of such reservoir characteristics (i.e., permeability)

show abstract

“…Hitherto, these methods are limited to a single phase fluid. There are two main classes of lattice Boltzmann solvers that incorporate effective media: force-adjusted and partial bounce back methods [16][17][18][19][20][21][22][23][24][25][26][27][28][29][30]. We chose to follow the partial bounce back method of Walsh et al [26] that combines computational simplicity with an ability to accurately link the partial bounce back fraction to the permeability.…”

Section: Single Phase Fluids and Effective Mediamentioning

confidence: 99%

“…The first lattice models incorporating effective media properties were developed as early as the 1990s [16][17][18][19][20][21] but they started gaining traction only recently [22][23][24][25][26][27][28][29][30]. In these models, in addition to lattice nodes assigned to be either fluid space or solid, some nodes are assigned effective media properties.…”

Section: Introductionmentioning

confidence: 99%

Gray free-energy multiphase lattice Boltzmann model with effective transport and wetting properties

et al. 2016

View full text Add to dashboard Cite

The paper shows that it is possible to combine the free energy lattice Boltzmann approach to multi-phase modelling of fluids involving both liquid and vapour with the partial bounce back lattice Boltzmann approach to modelling effective media. Effective media models are designed to mimic the properties of porous materials with porosity much finer than the scale of the simulation lattice. In the partial bounce back approach, an effective media parameter or bounce back fraction controls fluid transport. In the combined model, a wetting potential is additionally introduced that controls the wetting properties of the fluid with respect to interfaces between free space (white nodes), effective media (grey nodes) and solids (black nodes). The use of the wetting potential combined with the bounce back parameter gives the model the ability to simulate transport and sorption of a wide range of fluid / material systems. Results for phase separation, permeability, contact angle and wicking in grey media are shown. Sorption is explored in small sections of model multi-scale porous systems to demonstrate two-step desorption, sorption hysteresis and the ink-bottle effect.

show abstract

Lattice Boltzmann modeling and evaluation of fluid flow in heterogeneous porous media involving multiple matrix constituents

Cited by 41 publications

References 22 publications

Apparent permeability prediction of organic shale with generalized lattice Boltzmann model considering surface diffusion effect

Apparent permeability prediction of organic shale with generalized lattice Boltzmann model considering surface diffusion effect

Recent development in numerical simulation of enhanced geothermal reservoirs

Gray free-energy multiphase lattice Boltzmann model with effective transport and wetting properties

Contact Info

Product

Resources

About