High-resolution mapping of yield curve shape and evolution for high porosity sandstones

The onset of permanent inelastic deformation for porous rock is typically defined by a yield curve plotted in P-Q space, where P is the effective mean stress and Q is the differential stress. Sandstones usually have broadly elliptical shaped yield curves, with the low pressure side of the ellipse associated with localized brittle faulting (dilation) and the high pressure side with distributed ductile deformation (compaction). However recent works have shown that these curves might not be perfectly elliptical and that significant evolution in shape occurs with continued deformation. We therefore use a novel stress-probing methodology to map in high-resolution the yield curve shape for Boise and Idaho Gray sandstones (36-38% porosity) and also investigate curve evolution with increasing deformation. The data reveal yield curves with a much flatter geometry than previously recorded for porous sandstone and that the compactive side of the curve is partly comprised of a near vertical limb. The yield curve evolution is found to be strongly dependent on the nature of inelastic strain. Samples that were compacted under a deviatoric load, with a component of inelastic shear strain, were found to have yield curves with peaks that are approximately 50% higher than similar porosity samples that were hydrostatically compacted (i.e. purely volumetric strain). The difference in yield curve evolution along the different loading paths is attributed to mechanical anisotropy that develops during deviatoric loading by the closure of preferentially orientated fractures. Increased shear strain also leads to the formation of a plateau at the peak of the yield curve as samples deform along the deviatoric loading path. These results have important implications for understanding how the strength of porous rock evolves along different stress paths, including during fluid extraction from hydrocarbon reservoirs where the stress state is rarely isotropic.