Digital.CSIC. Repositorio Institucional del CSIC

oai:digital.csic.es:10261/284837

The code contains a MATLAB m.file, named “GridBasedPorosityReconstruction.m”, to reconstruct a 3D grid-based porosity map for centimeter-scale rock cores. The code reads binarized cross-section CT images collected along the core length (examples from a limestone core (Fig. 1) are provided in the accompanying file “IMG.zip”) as well as the number and size of grid blocks for domain discretization as inputs. The code first correlates the pixel domain with the arbitrarily-sized grid blocks in each image and then calculates the porosity for all cells. It enables shifting the average rock porosity by a certain value (e.g. the difference between average CT and effective rock porosities) by randomly distributing pores on different CT images. An explicit permeability-porosity relationship is used to build a grid-based permeability map from that of porosity across the core. The code plots 2D and 3D distribution maps of rock porosity and exports the derived data in Excel files. The data can be used in reactive transport modeling (e.g. CrunchFlow code (Steefel et al., 2015; Comput. Geosci. 19, 445–478)). Detailed explanations of the performed calculations are provided as stepwise comments in the m.file to allow the user to set modifications so as to obtain, for example, the distribution maps of any other rock characteristics (e.g., elastic, thermal, and electromagnetic constants) related to porosity., Peer reviewed