Photoelastic and numerical investigation of stress distributions around fault models under biaxial compressive loading conditions

This paper presents the results of a photoelastic and numerical study of stress distributions (contours and trajectories) around fault models of various geometries, submitted to a biaxial compressive load. It aims to describe typical biaxial stress behaviours and emphasize the existing differences with the well-known uniaxial compressive load case. Stress trajectories are sometimes shown by joint sets acting as markers of a paleostress field, and they can be interpreted by particular shallow tension–compression situations. At depth fractures can be reactivated, or can dilate under conditions of triaxial compression, and behaviour is essentially controlled by a high stress ratio (high σ₃/σ₁ ratio). In spite of the potential importance of such stress states on fracture permeability and fluid flow, and although they are frequently found at depth in a reservoir context, such stress conditions are poorly investigated, particularly in terms of stress perturbations.The presented analogue experiments consisted of compression tests done on polymethylmethacrylate (PMMA) plates; this material has mechanical properties comparable to those of brittle rocks in the upper crust, and presents birefringence. The samples contained open defects acting as faults, and the stress trajectories around these faults were investigated using a photoelastic device. Comparable numerical experiments were realised with a finite-element code (Franc 2D), using frictionless fault models.First, the effect of an increasing biaxial compressive load ratio σ₃/σ₁ on stress trajectories around an isolated open defect was explored. It was shown that the stress trajectories were drastically modified when σ₃/σ₁ increased from 0.2 to 0.4, this result being consistent with previous studies. In particular, when σ₃/σ₁ was superior or equal to 0.4, external isotropic points around where trajectories diverged (called repulsive isotropic points) appeared near the tips, but away from the defects. They tended to move away from the defects towards the main load direction when σ₃/σ₁ increased. The described isotropic points were points of stress trajectory divergence, i.e. points where stress decreased, implying that zones around them were strongly unfavourable to shear reactivation.Second, stress trajectories around fault models of various geometries (oblique isolated defect, dilational jog, compressive jog, and complex patterns) were studied, the applied biaxial compressive load ratio being 0.7. These biaxial stress trajectories were compared with similar uniaxial stress trajectories in order to explore the existing differences between the two regimes. It enabled new stress trajectory geometries to be described and interpreted. In particular, numerous external repulsive isotropic points were observed, and defect tips were shown to be zones of high convergence of stress trajectories. Furthermore, in contrast to the uniaxial compressive load case, stress trajectories and stress contours were geometrically similar within dilational and compressive jogs under biaxial compressive load. For both jog types, the centre of the overlapping zone, and the areas along the wall of a jog-defining fault and facing the overlapping tip of the other fault, proved to be zones of low mean stress, implying that fluids may migrate towards them from the tips of the faults, in response to mean stress gradients. Furthermore, the centre of the jogs exhibiting a high differential stress was a favourable area in terms of fracture reactivation, which may facilitate fluids transfer and storage.It was also observed that for both the uniaxial and biaxial compressive loads, isotropic zones were localised at the acute angle between branching defects.