Reconciliation of contrasting theories for fracture spacing in layered rocks

Natural and man-made brittle layers embedded in a weaker matrix and subjected to layer-parallel extension typically develop an array of opening-mode fractures with a remarkably regular spacing. This spacing often scales with layer thickness, and it decreases as extension increases until fracture saturation is reached. Existing analytical one-dimensional (1-D) ‘full-slip’ models, which assume that interfacial slip occurs over the entire length of the fracture-bound blocks, predict that the ratio of fracture spacing to layer thickness at saturation is proportional to the ratio of layer tensile to interface shear strength (T/τ). Using 2-D discontinuum mechanical models run for conditions appropriate to layered rocks, we show that fracture spacing at saturation decreases linearly with decreasing T/τ ratio, as predicted by 1-D models. At low T/τ ratios (ca. <3.0), however, interfacial slip is suppressed and the heterogeneous 2-D stress distribution within fracture-bound blocks controls further fracture nucleation, as predicted by an existing 2-D ‘fracture infill criterion’. The applicability of the two theories is hence T/τ ratio dependent. Our models illustrate that fracture spacing in systems permitting interfacial slip is not necessarily an indicator of fracture system maturity. Fracture spacing is expected to decrease with increasing overburden pressure and decreasing layer tensile strength.