Numerical investigation of pore-scale heat conduction in rock matrix blocks subject to time-dependent fracture thermal conditions
- 1Department of Petroleum and Geo-Energy Engineering, Amirkabir University of Technology
- 2Chemical Engineering Department, Tarbiat Modarres University
Abstract
Geothermal energy extraction from fractured reservoirs is strongly controlled by heat transfer between the rock matrix and fractures. Although fractures act as the primary flow pathways, the low-permeability matrix often dominates the thermal response of the system. In this study, transient heat conduction within heterogeneous rock matrix blocks is investigated at the pore scale using binarized micro-CT images of sandstone and carbonate rocks. Fractures are not explicitly resolved in the computational domain; instead, they are represented by imposing time-dependent thermal boundary conditions on the external surfaces of the matrix blocks. The heat conduction equation is solved using an explicit finite difference method on the voxel grid, with temperature-dependent thermal properties assigned to solid grains and water-filled pores. The effects of mineralogy, porosity, and the spatial distribution of pores and grains on heat diffusion and the matrix–fracture heat transfer shape factor are examined. In addition, the influence of block size and proximity to the injection well, represented by constant and linearly varying fracture temperature boundary conditions, is analyzed. The results demonstrate that pore-scale heterogeneity of the rock matrix plays a significant role in controlling heat transfer in fractured geothermal reservoirs.
Related
- Experimental and numerical investigation of voltage-dependent active vibration control on a notch-cracked composite mono leaf spring
- Numerical investigation of a finite element model to inform the design of a synthetic abdominal wall
- Numerical investigation of radiative chemically reactive induced-mhd burgers’ nanofluid flow over a nonlinear stretching sheet
- Experimental and numerical investigation of the synergistic rock-breaking mechanisms of double disc cutters under multi-mode loading
- A reliability analysis method for Christmas tree in operating conditions based on FEMCA, FTA and numerical simulation
BibTeX
@article{abbasi2026numerical,
title = {{Numerical investigation of pore-scale heat conduction in rock matrix blocks subject to time-dependent fracture thermal conditions}},
author = {Mahdi Abbasi and Arman Aryanzadeh},
journal = {Scientific Reports},
year = {2026},
month = jul,
doi = {10.1038/s41598-026-63044-7},
url = {https://doi.org/10.1038/s41598-026-63044-7},
note = {CC BY},
}
RIS
TY - JOUR
TI - Numerical investigation of pore-scale heat conduction in rock matrix blocks subject to time-dependent fracture thermal conditions
AU - Mahdi Abbasi
AU - Arman Aryanzadeh
PY - 2026
DA - 2026/07/01
JO - Scientific Reports
DO - 10.1038/s41598-026-63044-7
UR - https://doi.org/10.1038/s41598-026-63044-7
AB - Geothermal energy extraction from fractured reservoirs is strongly controlled by heat transfer between the rock matrix and fractures. Although fractures act as the primary flow pathways, the low-permeability matrix often dominates the thermal response of the system. In this study, transient heat conduction within heterogeneous rock matrix blocks is investigated at the pore scale using binarized micro-CT images of sandstone and carbonate rocks. Fractures are not explicitly resolved in the computational domain; instead, they are represented by imposing time-dependent thermal boundary conditions on the external surfaces of the matrix blocks. The heat conduction equation is solved using an explicit finite difference method on the voxel grid, with temperature-dependent thermal properties assigned to solid grains and water-filled pores. The effects of mineralogy, porosity, and the spatial distribution of pores and grains on heat diffusion and the matrix–fracture heat transfer shape factor are examined. In addition, the influence of block size and proximity to the injection well, represented by constant and linearly varying fracture temperature boundary conditions, is analyzed. The results demonstrate that pore-scale heterogeneity of the rock matrix plays a significant role in controlling heat transfer in fractured geothermal reservoirs.
SN - 2045-2322
ER -