Influence of damage on pore size distribution and permeability of rocks

The model proposed in this article relates permeability to porosity measurements that can easily be performed in the laboratory. The pore size distribution (PSD) curve is updated with strains and damage. The updated volumetric fractions of natural pores and cracks are introduced in the expression of permeability. Contrary to classical permeability models based on PSD integrations, the model proposed in this article accounts for possible changes in the porosity modes: one mode for undamaged samples and two modes for cracked samples. The proposed approach also accounts for varying states of damage, as opposed to classical fracture network models, in which the cracks pattern is fixed. The only material parameters that are required to describe the microstructure are the lower and upper bounds of the pores size for both natural pores and cracks. All the other PSD parameters involved in the model are related to macroscopic parameters that can easily be determined in the laboratory, such as the initial void ratio. The framework proposed in this article can be used in any damage constitutive model to determine the permeability of a brittle porous medium. Drained triaxial compression tests have been simulated. Before cracks initiation, permeability decreases while the larger natural pores are getting squeezed. After the occurrence of damage, permeability grows due to the increase of cracks density. The model performs well to represent the influence of the confining pressure on damage evolution and permeability variations. Copyright (c) 2012 John Wiley & Sons, Ltd.