A lattice Boltzmann study on the impact of the geometrical properties of porous media on the steady state relative permeabilities on two-phase immiscible flows

2016 ◽  
Vol 95 ◽  
pp. 61-79 ◽  
Author(s):  
Duo Zhang ◽  
K. Papadikis ◽  
Sai Gu
2020 ◽  
Vol 146 ◽  
pp. 03002
Author(s):  
Marios S. Valavanides ◽  
Matthieu Mascle ◽  
Souhail Youssef ◽  
Olga Vizika

The phenomenology of steady-state two-phase flow in porous media is recorded in SCAL relative permeability diagrams. Conventionally, relative permeabilities are considered to be functions of saturation. Yet, this has been put into challenge by theoretical, numerical and laboratory studies that have revealed a significant dependency on the flow rates. These studies suggest that relative permeability models should include the functional dependence on flow intensities. Just recently a general form of dependence has been inferred, based on extensive simulations with the DeProF model for steady-state two-phase flows in pore networks. The simulations revealed a systematic dependence of the relative permeabilities on the local flow rate intensities that can be described analytically by a universal scaling functional form of the actual independent variables of the process, namely, the capillary number, Ca, and the flow rate ratio, r. In this work, we present the preliminary results of a systematic laboratory study using a high throughput core-flood experimentation setup, whereby SCAL measurements have been taken on a sandstone core across different flow conditions -spanning 6 orders of magnitude on Ca and r. The scope is to provide a preliminary proof-of-concept, to assess the applicability of the model and validate its specificity. The proposed scaling opens new possibilities in improving SCAL protocols and other important applications, e.g. field scale simulators.


SPE Journal ◽  
2010 ◽  
Vol 15 (04) ◽  
pp. 917-927 ◽  
Author(s):  
Thomas Ramstad ◽  
Pål-Eric Øren ◽  
Stig Bakke

Summary We present results from simulations of two-phase flow directly on digitized rock-microstructure images of porous media using a lattice Boltzmann (LB) method. The implemented method is performed on a D3Q19 lattice with fluid/fluid and fluid/solid interaction rules to handle interfacial tension and wetting properties. We demonstrate that the model accurately reproduces capillary and wetting effects in pores with a noncircular shape. The model is applied to study viscous coupling effects for two-phase concurrent annular flow in circular tubes. Simulated relative permeabilities for this case agree with analytical predictions and show that the nonwetting-phase relative permeability might greatly exceed unity when the wetting phase is less viscous than the nonwetting phase. Two-phase LB simulations are performed on microstructure images derived from X-ray microtomography and process-based reconstructions of Bentheimer sandstone. By imposing a flow regulator to control the capillary number of the flow, the LB model can closely mimic typical experimental setups, such as centrifuge capillary pressure and unsteady- and steady-state relative permeability measurements. Computed drainage capillary pressure curves are found to be in excellent agreement with experimental data. Simulated steady-state relative permeabilities at typical capillary numbers in the vicinity of 10−5 are in fair agreement with measured data. The simulations accurately reproduce the wetting-phase relative permeability but tend to underpredict the nonwetting-phase relative permeability at high wetting-phase saturations. We explain this by pointing to percolation threshold effects of the samples. For higher capillary numbers, we correctly observe increased relative permeability for the nonwetting phase caused by mobilization and flow of trapped fluid. It is concluded that the LB model is a powerful and promising tool for deriving physically meaningful constitutive relations directly from rock-microstructure images.


Author(s):  
Wojciech Sobieski

AbstractThe paper describes the so-called Waterfall Algorithm, which may be used to calculate a set of parameters characterising the spatial structure of granular porous media, such as shift ratio, collision density ratio, consolidation ratio, path length and minimum tortuosity. The study is performed for 1800 different two-dimensional random pore structures. In each geometry, 100 individual paths are calculated. The impact of porosity and the particle size on the above-mentioned parameters is investigated. It was stated in the paper, that the minimum tortuosity calculated by the Waterfall Algorithm cannot be used directly as a representative tortuosity of pore channels in the Kozeny or the Carman meaning. However, it may be used indirect by making the assumption that a unambiguous relationship between the representative tortuosity and the minimum tortuosity exists. It was also stated, that the new parameters defined in the present study are sensitive on the porosity and the particle size and may be therefore applied as indicators of the geometry structure of granular media. The Waterfall Algorithm is compared with other methods of determining the tortuosity: A-Star Algorithm, Path Searching Algorithm, Random Walk technique, Path Tracking Method and the methodology of calculating the hydraulic tortuosity based on the Lattice Boltzmann Method. A very short calculation time is the main advantage of the Waterfall Algorithm, what meant, that it may be applied in a very large granular porous media.


Author(s):  
Andreas G. Yiotis ◽  
John Psihogios ◽  
Michael E. Kainourgiakis ◽  
Aggelos Papaioannou ◽  
Athanassios K. Stubos

Sign in / Sign up

Export Citation Format

Share Document