Hydraulic fracturing has been inferred to trigger the majority of injection-induced earthquakes in western Canada, in contrast to the Midwestern United States, where massive saltwater disposal is the dominant triggering mechanism. A template-based earthquake catalog from a seismically active Canadian shale play, combined with comprehensive injection data during a 4-month interval, shows that earthquakes are tightly clustered in space and time near hydraulic fracturing sites. The largest event [moment magnitude (MW) 3.9] occurred several weeks after injection along a fault that appears to extend from the injection zone into crystalline basement. Patterns of seismicity indicate that stress changes during operations can activate fault slip to an offset distance of >1 km, whereas pressurization by hydraulic fracturing into a fault yields episodic seismicity that can persist for months.

https://science.sciencemag.org/content/sci/354/6318/1406.full.pdf

Fault activation by hydraulic fracturing in western Canada.

3D computer simulations and experiments are employed to study random packings of compressible spherical grains under external confining stress. In the rigid ball limit, we find a continuous transition in which the stress vanishes as (straight phi-straight phi(c))(beta), where straight phi is the (solid phase) volume density. The value of straight phi(c) depends on whether the grains interact via only normal forces (giving rise to random close packings) or by a combination of normal and friction generated transverse forces (producing random loose packings). In both cases, near the transition, the system's response is controlled by localized force chains.

Packings of Compressible Granular Materials

A recently developed unconventional fracture model (UFM) is able to simulate complex fracture network propagation in a formation with pre-existing natural fractures. A method for computing the stress shadow from fracture branches in a complex hydraulic fracture network (HFN) based on an enhanced 2D displacement discontinuity method with correction for finite fracture height is implemented in UFM and is presented in detail including approach validation and examples. The influence of stress shadow effect from the HFN generated at previous treatment stage on the HFN propagation and shape at new stage is also discussed.

Numerical Modeling of Hydraulic Fractures Interaction In Complex Naturally Fractured Formations

Numerical Modeling of Hydraulic Fractures Interaction in Complex Naturally Fractured Formations

Multistage fracturing of the horizontal well is recognized as the main stimulation technology for shale gas development. The hydraulic fracture geometry and stimulated reservoir volume (SRV) is interpreted by using the microseismic mapping technology. In this paper, we used a computerized tomography (CT) scanning technique to reveal the fracture geometry created in natural bedding-developed shale (cubic block of 30 cm × 30 cm × 30 cm) by laboratory fracturing. Experimental results show that partially opened bedding planes are helpful in increasing fracture complexity in shale. However, they tend to dominate fracture patterns for vertical stress difference Δσ
 v
  ≤ 6 MPa, which decreases the vertical fracture number, resulting in the minimum SRV. A uniformly distributed complex fracture network requires the induced hydraulic fractures that can connect the pre-existing fractures as well as pulverize the continuum rock mass. In typical shale with a narrow (<0.05 mm) and closed natural fracture system, it is likely to create complex fracture for horizontal stress difference Δσ
 h
  ≤ 6 MPa and simple transverse fracture for Δσ
 h
  ≥ 9 MPa. However, high naturally fractured shale with a wide open natural fracture system (>0.1 mm) does not agree with the rule that low Δσ
 h
 is favorable for uniformly creating a complex fracture network in zone. In such case, a moderate Δσ
 h
 from 3 to 6 MPa is favorable for both the growth of new hydraulic fractures and the activation of a natural fracture system. Shale bedding, natural fracture, and geostress are objective formation conditions that we cannot change; we can only maximize the fracture complexity by controlling the engineering design for fluid viscosity, flow rate, and well completion type. Variable flow rate fracturing with low-viscosity slickwater fluid of 2.5 mPa s was proved to be an effective treatment to improve the connectivity of induced hydraulic fracture with pre-existing fractures. Moreover, the simultaneous fracturing can effectively reduce the stress difference and increase the fracture number, making it possible to generate a large-scale complex fracture network, even for high Δσ
 h
 from 6 MPa to 12 MPa.

Ripudaman Manchanda

Papers

Diagnosing Multi-Cluster Fracture Propagation Using Dynamic Poroelastic Pressure Transient Analysis

Fully 3-D Simulation of Diagnostic Fracture Injection Tests with Application in Depleted Reservoirs

An extended finite volume model for implicit cohesive zone fracture propagation in a poroelastic medium

Analyzing Pressure Interference Between Horizontal Wells During Fracturing

Interpreting Inter-Well Poroelastic Pressure Transient Data: An Analytical Approach Validated with Field Case Studies