Direct simulation of flow on microcomputed-tomography (micro-CT) images of rocks is widely used for the calculation of permeability. However, direct numerical methods are computationally demanding. A rapid and robust method… Click to show full abstract
Direct simulation of flow on microcomputed-tomography (micro-CT) images of rocks is widely used for the calculation of permeability. However, direct numerical methods are computationally demanding. A rapid and robust method is proposed to solve the elliptic flow equation. Segmented micro-CT images are used for the calculation of local conductivity in each voxel. The elliptic flow equation is then solved on the images using the finite-volume method. The numerical method is optimized in terms of memory usage using sparse matrix modules to eliminate memory overhead associated with both the inherent sparsity of the finite-volume two-point flux-approximation (TPFA) method, and the presence of zero conductivity for impermeable grain cells. The estimated permeabilities for a number of sandstone and carbonate micro-CT images are compared against estimation of other solvers, and results show a difference of approximately 11%. However, the computational time is 80% lower. Local conductivity can furthermore be assigned directly into matrix voxels without a loss in generality, hence allowing the pore-scale finite-volume solver (PFVS) to be able to solve for flow in a permeable matrix as well as open pore space. This has been developed to include the effect of microporosity in flow simulation.
               
Click one of the above tabs to view related content.