Multiscale Lattice Boltzmann Method for flow simulations in highly heterogenous porous media
Donald L. Brown,Jun Li,V. Calo,Y. Efendiev,O. Iliev
DOI: https://doi.org/10.2118/165985-MS
2013-09-16
Abstract:A lattice Boltzmann method (LBM) for flow simulations in highly heterogeneous porous media at both pore and Darcy scales is proposed in the paper. In the pore scale simulations, flow of two phases (e.g., oil and gas) or two immiscible fluids (e.g., water and oil) are modeled using cohesive or repulsive forces, respectively. The relative permeability can be computed using pore-scale simulations and seamlessly applied for intermediate and Darcy-scale simulations. A multiscale LBM that can reduce the computational complexity of existing LBM and transfer the information between different scales is implemented. The results of coarse-grid, reduced-order, simulations agree very well with the averaged results obtained using fine grid. Introduction Detailed flow simulations in porous media are modeled using the Darcy or Brinkman approximations. In these models, the effective parameters such as absolute and relative permeabilities depend on the pore-scale geometry. To compute these effective parameters, pore-scale simulations accounting for relevant geometric features in a Representative Elementary Volume (REV) are commonly used as in (Khan et al. (2012)). The lattice Boltzmann method (LBM), see (McNamara and Zanetti (1988); Chen et al. (1992); Qian et al. (1992)), is well developed for pore scale flow simulations and extended to model two-phase systems or two immiscible fluids (Shan and Chen (1993); Shan and Chen (1994)). After computing the effective parameters, we can perform Darcy-scale simulations using traditional finite volume or element method. However, these computations are limited to small REVs (compared to the computational grid) and rely on two distinct idealized scale concepts. In previous findings, simulation techniques using the modified LBM algorithm is proposed (Kang et al. (2002); Guo and Zhao (2002)), where LBM is considered as a unified framework for simulations at all scales. However, these simulations require significant computational resources to converge since the permeability distribution usually has drastic changes in space, which requires very fine grid for the high spatial resolution. To overcome this difficulty, we propose a multiscale LBM algorithm that is applicable at the pore and coarser, e.g., Darcy, scales. The proposed algorithm is capable of simulating single and multiphase flows and appropriately compute coarse-grid flow properties. The proposed multiscale LBM uses coarser grid to represent the effects of the fine-grid (pore-scale) simulations. For the multiscale LBM, the effective properties and reduced-order models are proposed as we coarsen the grid. To show the convenience of LBM in modeling complex physical processes in porous media, a single component system is modeled first to show phase separation process. The location of each phase inside the geometry depends on the contact angle of the materials simulated. Then, the flow of two immiscible fluids inside a cross channel is modeled, which shows the applicability of different boundary conditions in LBM simulations. For the coarse scale simulation, we reduce the memory usage and improve the efficiency using coarse grid (see Efendiev and Hou (2009) for general overview of multiscale methods). The agreement between the coarse and fine grids simulations demonstrates the validity of the multiscale technique. Lattice Boltzmann method for pore scale simulations based on Shan-Chen model General algorithm The grid (lattice points) are uniformly distributed inside the computational domain and all computational quantities are defined at those discrete points. We use ∆x and ∆t for the grid distance and time step, respectively.
Environmental Science,Physics,Engineering