Preferential Paths of Air-water Two-phase Flow in Porous Structures with Special Consideration of Channel Thickness Effects.

Accurate understanding and predicting the flow paths of immiscible two-phase flow in rocky porous structures are of critical importance for the evaluation of oil or gas recovery and prediction of rock slides caused by gas-liquid flow. A 2D phase field model was established for compressible air-water two-phase flow in heterogenous porous structures. The dynamic characteristics of air-water two-phase interface and preferential paths in porous structures were simulated. The factors affecting the path selection of two-phase flow in porous structures were analyzed. Transparent physical models of complex porous structures were prepared using 3D printing technology. Tracer dye was used to visually observe the flow characteristics and path selection in air-water two-phase displacement experiments. The experimental observations agree with the numerical results used to validate the accuracy of phase field model. The effects of channel thickness on the air-water two-phase flow behavior and paths in porous structures were also analyzed. The results indicate that thick channels can induce secondary air flow paths due to the increase in flow resistance; consequently, the flow distribution is different from that in narrow channels. This study provides a new reference for quantitatively analyzing multi-phase flow and predicting the preferential paths of immiscible fluids in porous structures

Convex hull approach for determining rock representative elementary volume for multiple petrophysical parameters using pore-scale imaging and Lattice-Boltzmann modelling

In the last decade, the study of fluid flow in porous media has developed considerably due to the combination of X-ray Micro Computed Tomography (micro-CT) and advances in computational methods for solving complex fluid flow equations directly or indirectly on reconstructed three-dimensional pore space images. In this study, we calculate porosity and single phase permeability using micro-CT imaging and Lattice Boltzmann (LB) simulations for 8 different porous media: beadpacks (with bead sizes 50 µm and 350 µm), sandpacks (LV60 and HST95), sandstones (Berea, Clashach and Doddington) and a carbonate (Ketton). Combining the observed porosity and calculated single phase permeability, we shed new light on the existence and size of the Representative Element of Volume (REV) capturing the different scales of heterogeneity from the pore-scale imaging. Our study applies the concept of the ‘Convex Hull’ to calculate the REV by considering the two main macroscopic petrophysical parameters, porosity and single phase permeability, simultaneously. The shape of the hull can be used to identify strong correlation between the parameters or greatly differing convergence rates. To further enhance computational efficiency we note that the area of the convex hull (for well-chosen parameters such as the log of the permeability and the porosity) decays exponentially with sub-sample size so that only a few small simulations are needed to determine the system size needed to calculate the parameters to high accuracy (small convex hull area). Finally we propose using a characteristic length such as the pore size to choose an efficient absolute voxel size for the numerical rock

Multiscale prediction of acoustic properties for glass wools: Computational study and experimental validation

International audienceThis work is concerned with the multiscale prediction of the transport and sound absorption propertiesassociated with industrial glass wool samples. In the first step, an experimental characterizationis performed on various products using optical granulometry and porosity measurements. A morphological analysis, based on scanning electron imaging, is further conducted to identify the probability density functions associated with the fiber angular orientation. The key morphological characterization parameters of the microstructure, which serve as input parameters of the model, include the porosity, the weighted volume diameter accounting for both lengths and diameters of the analyzed fibers (and therefore the specific surface area of the random fibrous material), and the preferred out-of-plane fiber orientation generated by the manufacturing process. A computational framework is subsequently proposed and allows for the reconstruction of an equivalent fibrous network. A fully stochastic microstructural model, parameterized by the probability laws inferred from the database, is also proposed herein. Multiscale simulations are carried out to estimate transport properties and sound absorption. With no adjustable parameter, the results accounting for ten different samples obtained with various processing parameters are finally compared with the experimental data and used to assess the relevance of the reconstruction procedures and the multiscale computations