1
|
Bultreys T, Ellman S, Schlepütz CM, Boone MN, Pakkaner GK, Wang S, Borji M, Van Offenwert S, Moazami Goudarzi N, Goethals W, Winardhi CW, Cnudde V. 4D microvelocimetry reveals multiphase flow field perturbations in porous media. Proc Natl Acad Sci U S A 2024; 121:e2316723121. [PMID: 38478686 PMCID: PMC10962996 DOI: 10.1073/pnas.2316723121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2023] [Accepted: 02/04/2024] [Indexed: 03/27/2024] Open
Abstract
Many environmental and industrial processes depend on how fluids displace each other in porous materials. However, the flow dynamics that govern this process are still poorly understood, hampered by the lack of methods to measure flows in optically opaque, microscopic geometries. We introduce a 4D microvelocimetry method based on high-resolution X-ray computed tomography with fast imaging rates (up to 4 Hz). We use this to measure flow fields during unsteady-state drainage, injecting a viscous fluid into rock and filter samples. This provides experimental insight into the nonequilibrium energy dynamics of this process. We show that fluid displacements convert surface energy into kinetic energy. The latter corresponds to velocity perturbations in the pore-scale flow field behind the invading fluid front, reaching local velocities more than 40 times faster than the constant pump rate. The characteristic length scale of these perturbations exceeds the characteristic pore size by more than an order of magnitude. These flow field observations suggest that nonlocal dynamic effects may be long-ranged even at low capillary numbers, impacting the local viscous-capillary force balance and the representative elementary volume. Furthermore, the velocity perturbations can enhance unsaturated dispersive mixing and colloid transport and yet, are not accounted for in current models. Overall, this work shows that 4D X-ray velocimetry opens the way to solve long-standing fundamental questions regarding flow and transport in porous materials, underlying models of, e.g., groundwater pollution remediation and subsurface storage of CO2 and hydrogen.
Collapse
Affiliation(s)
- Tom Bultreys
- Ghent University Centre for X-ray Tomography (UGCT), Ghent University, Ghent9000, Belgium
- Department of Geology, Ghent University, Ghent9000, Belgium
| | - Sharon Ellman
- Ghent University Centre for X-ray Tomography (UGCT), Ghent University, Ghent9000, Belgium
- Department of Geology, Ghent University, Ghent9000, Belgium
| | | | - Matthieu N. Boone
- Ghent University Centre for X-ray Tomography (UGCT), Ghent University, Ghent9000, Belgium
- Department of Physics and Astronomy, Ghent University, Ghent9000, Belgium
| | - Gülce Kalyoncu Pakkaner
- Ghent University Centre for X-ray Tomography (UGCT), Ghent University, Ghent9000, Belgium
- Department of Geology, Ghent University, Ghent9000, Belgium
| | - Shan Wang
- Ghent University Centre for X-ray Tomography (UGCT), Ghent University, Ghent9000, Belgium
- Department of Geology, Ghent University, Ghent9000, Belgium
| | - Mostafa Borji
- Ghent University Centre for X-ray Tomography (UGCT), Ghent University, Ghent9000, Belgium
- Department of Geology, Ghent University, Ghent9000, Belgium
| | - Stefanie Van Offenwert
- Ghent University Centre for X-ray Tomography (UGCT), Ghent University, Ghent9000, Belgium
- Department of Geology, Ghent University, Ghent9000, Belgium
| | - Niloofar Moazami Goudarzi
- Ghent University Centre for X-ray Tomography (UGCT), Ghent University, Ghent9000, Belgium
- Department of Physics and Astronomy, Ghent University, Ghent9000, Belgium
| | - Wannes Goethals
- Ghent University Centre for X-ray Tomography (UGCT), Ghent University, Ghent9000, Belgium
- Department of Physics and Astronomy, Ghent University, Ghent9000, Belgium
| | - Chandra Widyananda Winardhi
- Ghent University Centre for X-ray Tomography (UGCT), Ghent University, Ghent9000, Belgium
- Department of Geology, Ghent University, Ghent9000, Belgium
| | - Veerle Cnudde
- Ghent University Centre for X-ray Tomography (UGCT), Ghent University, Ghent9000, Belgium
- Department of Geology, Ghent University, Ghent9000, Belgium
- Department of Earth Sciences, Utrecht University, CB Utrecht3584, The Netherlands
| |
Collapse
|
2
|
Integrating Pore-Scale Flow MRI and X-ray μCT for Validation of Numerical Flow Simulations in Porous Sedimentary Rocks. Transp Porous Media 2022. [DOI: 10.1007/s11242-022-01770-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Abstract
AbstractSingle-phase fluid flow velocity maps in Ketton and Estaillades carbonate rock core plugs are computed at a pore scale, using the lattice Boltzmann method (LBM) simulations performed directly on three-dimensional (3D) X-ray micro-computed tomography (µCT) images (≤ 7 µm spatial resolution) of the core plugs. The simulations are then benchmarked on a voxel-by-voxel and pore-by-pore basis to quantitative, 3D spatially resolved magnetic resonance imaging (MRI) flow velocity maps, acquired at 35 µm isotropic spatial resolution for flow of water through the same rock samples. Co-registration of the 3D experimental and simulated velocity maps and coarse-graining of the simulation to the same resolution as the experimental data allowed the data to be directly compared. First, the results are demonstrated for Ketton limestone rock, for which good qualitative and quantitative agreement was found between the simulated and experimental velocity maps. The flow-carrying microstructural features in Ketton rock are mostly larger than the spatial resolution of the µCT images, so that the segmented images are an adequate representation of the pore space. Second, the flow data are presented for Estaillades limestone, which presents a more heterogeneous case with microstructural features below the spatial resolution of the µCT images. Still, many of the complex flow patterns were qualitatively reproduced by the LBM simulation in this rock, although in some pores, noticeable differences between the LBM and MRI velocity maps were observed. It was shown that 80% of the flow (fractional summed z-velocities within pores) in the Estaillades rock sample is carried by just 10% of the number of macropores, which is an indication of the high structural heterogeneity of the rock; in the more homogeneous Ketton rock, 50% of the flow is carried by 10% of the macropores. By analysing the 3D MRI velocity map, it was found that approximately one-third of the total flow rate through the Estaillades rock is carried by microporosity—a porosity that is not captured at the spatial resolution of the µCT image.
Collapse
|