Computational Modelling of Multi-scale Solute Dispersion in Porous Media - An Approach Based on Stochastic Calculus

10.5772/3189 ◽  
2011 ◽  
2021 ◽  
Author(s):  
Andres Gonzalez ◽  
Zoya Heidari ◽  
Olivier Lopez

Abstract Depositional mechanisms of sediments and post-depositional process often cause spatial variation and heterogeneity in rock fabric, which can impact the directional dependency of petrophysical, electrical, and mechanical properties. Quantification of the directional dependency of the aforementioned properties is fundamental for the appropriate characterization of hydrocarbon-bearing reservoirs. Anisotropy quantification can be accomplished through numerical simulations of physical phenomena such as fluid flow, gas diffusion, and electric current conduction in porous media using multi-scale image data. Typically, the outcome of these simulations is a transport property (e.g., permeability). However, it is also possible to quantify the tortuosity of the media used as simulation domain, which is a fundamental descriptor of the microstructure of the rock. The objectives of this paper are (a) to quantify tortuosity anisotropy of porous media using multi-scale image data (i.e., whole-core CT-scan and micro-CT-scan image stacks) through simulation of electrical potential distribution, diffusion, and fluid flow, and (b) to compare electrical, diffusional, and hydraulic tortuosity. First, we pre-process the images (i.e., CT-scan images) to remove non-rock material visual elements (e.g., core barrel). Then, we perform image analysis to identify different phases in the raw images. Then, we proceed with the numerical simulations of electric potential distribution. The simulation results are utilized as inputs for a streamline algorithm and subsequent direction-dependent electrical tortuosity estimation. Next, we conduct numerical simulation of diffusion using a random walk algorithm. The distance covered by each walker in each cartesian direction is used to compute the direction-dependent diffusional tortuosity. Finally, we conduct fluid-flow simulations to obtain the velocity distribution and compute the direction-dependent hydraulic tortuosity. The simulations are conducted in the most continuous phase of the segmented whole-core CT-scan image stacks and in the segmented pore-space of the micro-CT-scan image stacks. Finally, the direction-dependent tortuosity values obtained with each technique are employed to assess the anisotropy of the evaluated samples. We tested the introduced workflow on dual energy whole-core CT-scan images and on smaller scale micro-CT-scan images. The whole-core CT-scan images were obtained from a siliciclastic depth interval, composed mainly by spiculites. Micro-CT-scan images we obtained from Berea Sandstone and Austin Chalk formations. We observed numerical differences in the estimates of direction-dependent electrical, diffusional, and hydraulic tortuosity for both types of image data employed. The highest numerical differences were observed when comparing electrical and hydraulic tortuosity with diffusional tortuosity. The observed differences were significant specially in anisotropic samples. The documented comparison provides useful insight in the selection process of techniques for estimation of tortuosity. The use of core-scale image data in the proposed workflow provides semi-continuous estimates of tortuosity and tortuosity anisotropy which is typically not attainable when using pore-scale images. Additionally, the semi-continuous nature of the tortuosity and tortuosity anisotropy estimates in whole-core CT-scan image data provides an excellent tool for the selection of core plugs coring locations.


2021 ◽  
Vol MA2021-01 (27) ◽  
pp. 961-961
Author(s):  
Marc Secanell ◽  
Seongyeop Jung ◽  
Alexandre Jarauta-Arabi ◽  
Fei Wei ◽  
Mayank Sabharwal ◽  
...  

Symmetry ◽  
2020 ◽  
Vol 12 (3) ◽  
pp. 364 ◽  
Author(s):  
Huimin Wang ◽  
Jianguo Wang ◽  
Xiaolin Wang ◽  
Andrew Chan

Low-permeability porous medium usually has asymmetric distributions of pore sizes and pore-throat tortuosity, thus has a non-linear flow behavior with an initial pressure gradient observed in experiments. A threshold pressure gradient (TPG) has been proposed as a crucial parameter to describe this non-linear flow behavior. However, the determination of this TPG is still unclear. This study provides multi-scale insights on the TPG in low-permeability porous media. First, a semi-empirical formula of TPG was proposed based on a macroscopic relationship with permeability, water saturation, and pore pressure, and verified by three sets of experimental data. Second, a fractal model of capillary tubes was developed to link this TPG formula with structural parameters of porous media (pore-size distribution fractal dimension and tortuosity fractal dimension), residual water saturation, and capillary pressure. The effect of pore structure complexity on the TPG is explicitly derived. It is found that the effects of water saturation and pore pressure on the TPG follow an exponential function and the TPG is a linear function of yield stress. These effects are also spatially asymmetric. Complex pore structures significantly affect the TPG only in the range of low porosity, but water saturation and yield stress have effects on a wider range of porosity. These results are meaningful to the understanding of non-linear flow mechanism in low-permeability reservoirs.


Sign in / Sign up

Export Citation Format

Share Document