Image-based modelling of rock non-linear deformation under low-stress levels
Abstract
Rock geophysical properties are widely reported to exhibit non-linear behaviours under low-stress conditions (below 10-20 MPa) before transitioning to the linear elastic stage, primarily due to the closure of microcracks and grain interfaces. Image-based modelling of rock deformation struggles to effectively characterise the microcrack closure effect because of the partial-volume effect, where image voxels are larger than microcracks and contain both pore and solid phases. This study presents a novel method to simulate non-linear rock deformation under elevated stress conditions. The method reconstructs digital rock models by treating partial-volume voxels as transitional phases that incorporate microcracks. By assigning intermediate elastic moduli and assuming that the pore portion within each partial-volume voxel deforms before the remaining solid content, the method employs the finite element method to simulate rock deformation and calculate the porosity of the deformed model. The method is tested on two Bentheimer sandstone models, and the results demonstrate its ability to predict the non-linear changes in porosity and elastic properties as the effective stress increases. This work provides a new pathway for image-based modelling of non-linear rock deformation considering the microcrack closure effect, offering valuable insights into the complex mechanical behaviour of rocks under confinement.