skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Speeding-up Simulation of Multiphase Flow in Digital Images of Heterogeneous Porous Media by Curvelet Transformation
Multiphase fluid flow in porous media is important to a wide variety of processes of fundamental scientific and practical importance. Developing a model for the pore space of porous media represents the first step for simulating such flows. With rapid increase in the computation power and advances in instrumentation and imaging processes, it has become feasible to carry out simulation of multiphase flow in two- and three-dimensional images of porous media, hence dispensing with development of models of pore space that are based on approximating their morphology. Image-based simulations are, however, very time consuming. We describe an approach for speeding-up image-based simulation of multiphase flow in porous media based on curvelet transformations, which are specifically designed for processing of images that contain complex curved surfaces. Most porous media contain correlations in their morphology and, therefore, their images carry redundant information that, in the curvelet transform space, can be removed efficiently and accurately in order to obtain a coarser image with which the computations are far less intensive. We utilize the methodology to simulate two-phase flow of oil and water in two-dimensional digital images of sandstone and carbonate samples, and demonstrate that while the results with the curvelet-processed images are as accurate as those with the original ones, the computations are speeded up by a factor of 110–150. Thus, the methodology opens the way toward achieving the ultimate goal of simulation of multiphase flow in porous media, namely, making image-based computations a standard practice.  more » « less
Award ID(s):
2000968
PAR ID:
10625566
Author(s) / Creator(s):
;
Editor(s):
Blunt, MJ
Publisher / Repository:
Springer
Date Published:
Journal Name:
Transport in Porous Media
Volume:
137
Issue:
1
ISSN:
0169-3913
Page Range / eLocation ID:
215 to 232
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; (, Advances in Water Resources)
    Advances in instrumentation have made it possible to obtain high-resolution images of heterogeneous porous media. Such advances and the rapid increase in computational power mean that direct numerical simulation of multiphase flow in two- and three-dimensional (3D) images of porous media is feasible and, therefore, models of pore space that represent simplification and approximation of the actual morphology are dispensable. The bottleneck for image-based simulation is its long computation time. We propose an approach for speeding-up simulation of multiphase flow in 3D images of porous media that utilizes curvelet transformations (CTs). The CTs are designed for denoising of images that contain complex curved surfaces, such as those of heterogeneous porous media. This is possible because the morphology of most porous media contain extended correlations, implying that their images carry redundant information that can be eliminated by a suitable CT. As a result, simulation of multiphase flow in the coarser images are far less computationally intensive. The new approach is used to simulate two-phase flow of CO2 and brine in a 3D image of Berea sandstone. We demonstrate that the results with the CT-processed image are as accurate as those with the original one, while computations are significantly faster by a speedup factor between one and more than two orders of magnitude 
    more » « less
  2. ; (, Physics of Fluids)
    Blunt, MJ (Ed.)
    Multiphase fluid flow in porous media is relevant to many fundamental scientific problems as well as numerous practical applications. With advances in instrumentations, it has become possible to obtain high-resolution three-dimensional (3D) images of complex porous media and use them directly in the simulation of multiphase flows. A prime method for carrying out such simulations is the color-fluid lattice Boltzmann method with multi-relaxation time (CFLB-MRT) collision operator. The simulations are, however, time consuming and intensive. We propose a method to accelerate image-based computations with the CFLB-MRT method, in which the 3D image is preprocessed by curvelet transforming it and eliminating those details that do not contribute significantly to multiphase flow. The coarsening is done by thresholding the image. After inverting the coarser image back to the real space, it is utilized in the simulation of multiphase flow by the CFLB-MRT approach. As the test of the method, we carry out simulation of a two-phase flow problem in which the porous media are initially saturated by brine or water, which is then displaced by CO2 or oil, injected into the pore space. The simulations are carried out with two types of sandstone. We show that the method accelerates the computations significantly by a factor of up to 35. 
    more » « less
  3. ; ; ; (, Physics of Fluids)
    This study aims to bridge length scales in immiscible multiphase flow simulation by connecting two published governing equations at the pore-scale and continuum-scale through a novel validation framework. We employ Niessner and Hassnaizadeh's [“A model for two-phase flow in porous media including fluid-fluid interfacial area,” Water Resour. Res. 44(8), W08439 (2008)] continuum-scale model for multiphase flow in porous media, combined with the geometric equation of state of McClure et al. [“Modeling geometric state for fluids in porous media: Evolution of the Euler characteristic,” Transp. Porous Med. 133(2), 229–250 (2020)]. Pore-scale fluid configurations simulated with the lattice-Boltzmann method are used to validate the continuum-scale results. We propose a mapping from the continuum-scale to pore-scale utilizing a generalized additive model to predict non-wetting phase Euler characteristics during imbibition, effectively bridging the continuum-to-pore length scale gap. Continuum-scale simulated measures of specific interfacial area, saturation, and capillary pressure are directly compared to up-scaled pore-scale simulation results. This research develops a numerical framework capable of capturing multiscale flow equations establishing a connection between pore-scale and continuum-scale simulations. 
    more » « less
  4. ; ; ; ; (, Transport in Porous Media)
    Abstract Wettability is one of the critical parameters affecting multiphase flow in porous media. The wettability is determined by the affinity of fluids to the rock surface, which varies due to factors such as mineral heterogeneity, roughness, ageing, and pore-space geometry. It is well known that wettability varies spatially in natural rocks, and it is still generally considered a constant parameter in pore-scale simulation studies. The accuracy of pore-scale simulation of multiphase flow in porous media is undermined by such inadequate wettability models. The advent of in situ visualization techniques, e.g. X-ray imaging and microtomography, enables us to characterize the spatial distribution of wetting more accurately. There are several approaches for such characterization. Most include the construction of a meshed surface of the interface surfaces in a segmented X-ray image and are known to have significant errors arising from insufficient resolution and surface-smoothing algorithms. This work presents a novel approach for spatial determination of wetting properties using local lattice-Boltzmann simulations. The scheme is computationally efficient as the segmented X-ray image is divided into subdomains before conducting the lattice-Boltzmann simulations, enabling fast simulations. To test the proposed method, it was applied to two synthetic cases with known wettability and three datasets of imaged fluid distributions. The wettability map was obtained for all samples using local lattice-Boltzmann calculations on trapped ganglia and optimization on surface affinity parameters. The results were quantitatively compared with a previously developed geometrical contact angle determination method. The two synthetic cases were used to validate the results of the developed workflow, as well as to compare the wettability results with the geometrical analysis method. It is shown that the developed workflow accurately characterizes the wetting state in the synthetic porous media with an acceptable uncertainty and is better to capture extreme wetting conditions. For the three datasets of imaged fluid distributions, our results show that the obtained contact angle distributions are consistent with the geometrical method. However, the obtained contact angle distributions tend to have a narrower span and are considered more realistic compared to the geometrical method. Finally, our results show the potential of the proposed scheme to efficiently obtain wettability maps of porous media using X-ray images of multiphase fluid distributions. The developed workflow can help for more accurate characterization of the wettability map in the porous media using limited experimental data, and hence more accurate digital rock analysis of multiphase flow in porous media. 
    more » « less
  5. ; (, Journal of Fluid Mechanics)
    null (Ed.)
    Numerical modelling of deformation in hydromechanical systems can be time-consuming using fully coupled classical numerical methods for large representative porous media samples. In this paper, we present a new two-way coupled partitioned fluid–solid system. The coupled system is applied for simulating geomechanical processes at the pore-scale. We track the deformation of the solid resulting from the drainage of resident fluids in the pores, as well as the evolution of fluid properties from dynamic loading. The finite element method is responsible for capturing the structural deformation in the coupled system while the dynamic pore network is used for modelling multiphase flow in the fluid subsystem. A fictitious fluid–solid interface is created at each pore network-finite element node junction via convex hulling, followed by data exchange using linear interpolation. The results show good agreement with a pre-existing coupled finite volume model and the computations are completed in much less time. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email