Microstructure-based multi-scale evaluation of fluid flow in an anthracite coal sample with partially-percolating voxels

摘要

Understanding fluid flow behavior in coal is of great significance for coal-bed methane exploration. X-ray CT and image segmentation have been widely used to extract pore network and generate flow field grids for flow simulation in coal samples. However, these techniques have fundamental limitations for the multi-scale characterization of coal samples, where the sub-voxel scale details could not be resolved for millimeter scale macroscopic samples. This makes it difficult to simulate the multi-scale flow behavior of fluid transport in coal sample with varying pore scales. The primary challenge is to make connection between simulation results of different scales. In the present work, multi-scale fluid flow in an anthracite coal sample was simulated by incorporating the data-constrained modeling (DCM), molecular dynamics (MD) method and partially-percolating lattice Boltzmann method (PP-LBM). In this multi-scale simulation method, three-dimensional (3D) flow field containing multi-scale structural information of the coal sample was generated by combining DCM with multi-energy synchrotron radiation CT. Multi-scale fluid flow was simulated by PP-LBM. In PP-LBM, an effective percolation fraction parameter which represents the effective volume fraction of the fluid that contributed to the flow for the voxel was used as a bridge to connect the fluid flow pattern of sub-voxel scales and voxel scales. The effective percolation fraction of a voxel versus its porosity was derived by MD simulations at the sub-voxel size level. The 3D distribution of fluid speed in the coal sample and its permeability were obtained by this multi-scale method. The numerical results are consistent with published laboratory measurements. Our proposed approach incorporated multi-scale effects and offered a more realistic fluid transport simulation method for a coal sample with varying pore size scales from the microscopic to macroscopic level. The method would be applicable for fluid transport simul