Laboratory experiments were performed on a polymethyl methacrylate (PMMA)-PMMA frictional interface in a direct shear apparatus in order to gain understanding of fault dynamics leading to gross rupture. Actual asperity sizes and locations along the interface were characterized using a pressure-sensitive film. Slow aseismic slip accumulated nonuniformly along the fault and showed dependency on the applied normal force-increased normal force resulted in higher slip gradients. The slow slip front propagated from the trailing (pushed) edge into a region of more densely distributed asperities at rates between 1 and 9.5 mm/s. Foreshocks were detected and displayed impulsive signals with source radii ranging between 0.21 and 1.09 mm; measurements made using the pressure-sensitive film were between 0.05 and 1.2 mm. The spatiotemporal clustering of foreshocks and their relation to the elastodynamic energy released was dependent on the normal force. In the region where foreshocks occurred, qualitative optical measurements of the asperities along the interface were used to visualize dynamic changes occurring during the slow slip phase. To better understand the nucleation process, a quasi-static asperity finite element (FE) model was developed and focused in the region where foreshocks clustered. The FE model consisted of 172 asperities, located and sized based on pressure-sensitive film measurements. The numerical model provides a plausible explanation as to why foreshocks cluster in space and observed a normal force dependency and lend credence to Ohnaka's nucleation model.

We investigate energy partitioning using seismological methods of seismic ruptures with estimated submicron levels of slip in the laboratory. Estimates inferred from recorded seismic waves are founded on microscale phenomenological friction experiments in the laboratory and appear to be constrained by inherent assumptions. In this concerted study, we build on the methods used to absolutely calibrate an array of piezoelectric transducers in a direct shear laboratory apparatus. We found that flat-broadband sensor behavior allowed us to study source extent parameters using spectral source models that are typically used to interpret small to moderate-sized earthquakes. We computed the corner frequencies, low-frequency plateaus, and high-frequency spectral falloff exponent using single-station assumptions. Moment magnitude ranged from −9

Abstract. We performed a series of 12 hydraulic stimulation experiments in a 20 × 20 × 20 m foliated, crystalline rock volume intersected by two distinct fault sets at the Grimsel Test Site, Switzerland. The goal of these experiments was to improve our understanding of stimulation processes associated with high-pressure fluid injection used for reservoir creation in enhanced or engineered geothermal systems. In the first six experiments, pre-existing fractures were stimulated to induce shear dilation and enhance permeability. Two types of shear zones were targeted for these hydroshearing experiments: i) ductile ones with intense foliation and ii) brittle-ductile ones associated with a fractured zone. The second series of six stimulations were performed in borehole intervals without natural fractures to initiate and propagate hydraulic fractures that connect the wellbore to the existing fracture network. The same injection protocol was used for all experiments within each stimulation series so that the differences observed will give insights into the effect of geology on the seismo-hydro-mechanical response rather than differences due to the injection protocols. Deformations and fluid pressure were monitored using a dense sensor network in boreholes surrounding the injection locations. Seismicity was recorded with sensitive in-situ acoustic emission sensors both in boreholes and at the tunnel walls. We observed high variability in the seismic response in terms of seismogenic indices, b-values, spatial and temporal evolution during both hydroshearing and hydrofracturing experiments, which we attribute to local geological heterogeneities. Seismicity was most pronounced for injections into the highly conductive brittle-ductile shear zones, while injectivity increase on these structures was only marginal. No significant differences between the seismic response of hydroshearing and hydrofracturing was identified, possibly because the hydrofractures interact with the same pre-existing fracture network that is reactivated during the hydroshearing experiments. Fault slip during the hydroshearing experiments was predominantly aseismic. The results of our hydraulic stimulations indicate that stimulation of short borehole intervals with limited fluid volumes (i.e., the concept of zonal insulation) may be an effective approach to limit induced seismic hazard if highly seismogenic structures can be avoided.

Abstract. We performed a series of 12 hydraulic stimulation experiments in a 20m×20m×20m foliated, crystalline rock volume intersected by two distinct fault sets at the Grimsel Test Site, Switzerland. The goal of these experiments was to improve our understanding of stimulation processes associated with high-pressure fluid injection used for reservoir creation in enhanced or engineered geothermal systems. In the first six experiments, pre-existing fractures were stimulated to induce shear dilation and enhance permeability. Two types of shear zones were targeted for these hydroshearing experiments: (i) ductile ones with intense foliation and (ii) brittle–ductile ones associated with a fractured zone. The second series of six stimulations were performed in borehole intervals without natural fractures to initiate and propagate hydraulic fractures that connect the wellbore to the existing fracture network. The same injection protocol was used for all experiments within each stimulation series so that the differences observed will give insights into the effect of geology on the seismo-hydromechanical response rather than differences due to the injection protocols. Deformations and fluid pressure were monitored using a dense sensor network in boreholes surrounding the injection locations. Seismicity was recorded with sensitive in situ acoustic emission sensors both in boreholes and at the tunnel walls. We observed high variability in the seismic response in terms of seismogenic indices, b values, and spatial and temporal evolution during both hydroshearing and hydrofracturing experiments, which we attribute to local geological heterogeneities. Seismicity was most pronounced for injections into the highly conductive brittle–ductile shear zones, while the injectivity increase on these structures was only marginal. No significant differences between the seismic response of hydroshearing and hydrofracturing was identified, possibly because the hydrofractures interact with the same pre-existing fracture network that is reactivated during the hydroshearing experiments. Fault slip during the hydroshearing experiments was predominantly aseismic. The results of our hydraulic stimulations indicate that stimulation of short borehole intervals with limited fluid volumes (i.e., the concept of zonal insulation) may be an effective approach to limit induced seismic hazard if highly seismogenic structures can be avoided.

Earthquake faults, and all frictional surfaces, establish contact through asperities. A detailed knowledge of how asperities form will enable a better understanding of the manner in which they communicate during foreshock failure sequences that are observed, leading to the larger main shock. We present results of experiments where a pressure sensitive film was used to map, size and measure the magnitudes of the normal stresses at asperities along a seismogenic section of a laboratory simulated fault. We measured seismicity acoustically and foreshocks were found to be the result of localized asperity failure during the nucleation phase of gross fault rupture. Since surface roughness plays an important role in how asperities are formed, two Hurst exponents were measured to characterize a highly worn interface using roughness profiles: (i) long wavelength estimates (H ∼ 0.45) and (ii) short wavelength estimates (H ∼ 0.8-1.2). The short wavelength roughness estimates were computed at the scale of single asperity junction points. Macroscopically, the number of asperities and real contact area increased with additional application of normal force while the mean normal stress remained constant.The ratio of real to nominal contact area was low -ranging from 0.02 < Ar/A 0 < 0.05-predicting that the asperities should be elastically independent of each other. Results from the pressure sensitive film showed that asperities were closely spaced and could not be treated as mechanically independent. Larger asperities carried both higher levels of average normal stress and higher levels of normal stress heterogeneity than smaller ones. Using linear stability theorem, the critical slip distance on foreshocking asperities was estimated to be d 0 ∼ 0.65-3 μm. The critical slip distance d 0 was ∼1.8-11.5 per cent of the premonitory slip needed to initiate gross fault rupture of the interface (20-40 μm) and the overall slip necessary to initiate gross fault rupture was on the order of the average asperity diameter (52 μm). Foreshocks may be due to a change in the critical slip distance, at localized sections of the fault, caused by the two distinct roughness profiles measured at short and long length scales.

Reservoir core measurements can help guide seismic monitoring of fluid-induced pressure variations in tight fractured reservoirs, including those targeted for supercritical [Formula: see text] injection. We have developed the first seismic-frequency “room-dry” measurements of fracture-specific shear stiffness, using artificially fractured standard granite samples with different degrees of mating, a well-mated tensile fracture from a dolomite reservoir core, as well as simple roughened polymethyl methacrylate (PMMA) surfaces. We have adapted a low-frequency (0.01–100 Hz) shear modulus and attenuation apparatus to explore the seismic signature of fractures and understand the mechanics of asperity contacts under a range of normal stress conditions. Our instrument is unique in its ability to measure at low-normal stresses (0.5–20 MPa), simulating “open” fractures in shallow or high-fluid-pressure reservoirs. The accuracy of our instrument is demonstrated by calibration and comparison with ultrasonic measurements and low-frequency direct shear measurements of intact samples from the literature. Pressure-sensitive film was used to measure real contact area of the fracture surfaces. The fractured shear modulus for most of the samples shows an exponential dependence on the real contact area. A simple numerical model, with one bonded circular asperity, predicts this behavior and matches the data for the simple PMMA surfaces. The rock surfaces reach their intact moduli at lower contact area than the model predicts, likely due to more complex geometry. Finally, we apply our results to a linear-slip interface model to estimate reflection coefficients and calculate S-wave time delays due to the lower-wave velocities through the fractured zone. We find that cross-well surveys could detect even well-mated hard-rock fractures, assuming the availability of high-repeatability acquisition systems.

Slow slip events (SSEs) represent a slow faulting process leading to aseismic strain release often accompanied by seismic tremor or earthquake swarms. The larger SSEs last longer and are often associated with intense and energetic tremor activity, suggesting that aseismic slip controls tremor genesis. A similar pattern has been observed for SSEs that trigger earthquake swarms, although no comparative studies exist on the source parameters of SSEs and tremor or earthquake swarms. We analyze the source scaling of SSEs and associated tremor- or swarm-like seismicity through our newly compiled dataset. We find a correlation between the aseismic and seismic moment release indicating that the shallower SSEs produce larger seismic moment release than deeper SSEs. The scaling may arise from the heterogeneous frictional and rheological properties of faults prone to SSEs and is mainly controlled by temperature. Our results indicate that similar physical phenomena govern tremor and earthquake swarms during SSEs.

A pressure-sensitive film was used to characterize the asperity contacts along a polymethyl methacrylate (PMMA) interface in the laboratory. The film has structural health monitoring (SHM) applications for flanges and other precision fittings and train rail condition monitoring. To calibrate the film, simple spherical indentation tests were performed and validated against a finite element model (FEM) to compare normal stress profiles. Experimental measurements of the normal stress profiles were within −7.7% to 6.6% of the numerical calculations between 12 and 50 MPa asperity normal stress. The film also possessed the capability of quantifying surface roughness, an important parameter when examining wear and attrition in SHM applications. A high definition video camera supplied data for photometric analysis (i.e., the measure of visible light) of asperities along the PMMA-PMMA interface in a direct shear configuration, taking advantage of the transparent nature of the sample material. Normal stress over individual asperities, calculated with the pressure-sensitive film, was compared to the light intensity transmitted through the interface. We found that the luminous intensity transmitted through individual asperities linearly increased 0.05643 ± 0.0012 candelas for an increase of 1 MPa in normal stress between normal stresses ranging from 23 to 33 MPa.