[CODE] A MATLAB code for digitized reconstruction of the rock porosity distribution from Computed Tomography (CT) images
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 acc...
| Autores: | , , , , |
|---|---|
| Tipo de recurso: | conjunto de datos |
| Fecha de publicación: | 2022 |
| País: | España |
| Institución: | Consejo Superior de Investigaciones Científicas (CSIC) |
| Repositorio: | DIGITAL.CSIC. Repositorio Institucional del CSIC |
| OAI Identifier: | oai:digital.csic.es:10261/284837 |
| Acceso en línea: | http://hdl.handle.net/10261/284837 |
| Access Level: | acceso abierto |
| Palabra clave: | Rock porosity distribution Computed Tomography (CT) images |
| Sumario: | 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. |
|---|