Anisotropy in Fracking: A Percolation Model for Observed Microseismicity

Hydraulic fracturing (fracking) using high pressures and a low viscosity fluid allow the extraction of large quantiles of oil and gas from very low permeability shale formations. The initial production of oil and gas at depth leads to high pressures and an extensive distribution of natural fractures which reduce the pressures. With time these fractures heal, sealing the remaining oil and gas in place. High volume fracking opens the healed fractures allowing the oil and gas to flow the horizontal productions wells. We model the injection process using invasion percolation. We utilize a 2D square lattice of bonds to model the sealed natural fractures. The bonds are assigned random strengths and the fluid, injected at a point, opens the weakest bond adjacent to the growing cluster of opened bonds. Our model exhibits burst dynamics in which the clusters extends rapidly into regions with weak bonds. We associate these bursts with the microseismic activity generated by fracking injections. A principal object of this paper is to study the role of anisotropic stress distributions. Bonds in the $y$-direction are assigned higher random strengths than bonds in the $x$-direction. We illustrate the spatial distribution of clusters and the spatial distribution of bursts (small earthquakes) for several degrees of anisotropy. The results are compared with observed distributions of microseismicity in a fracking injection. Both our bursts and the observed microseismicity satisfy Gutenberg-Richter frequency-size statistics.Comment: 14 pages, 10 figure

Natural Time, Nowcasting and the Physics of Earthquakes: Estimation of Seismic Risk to Global Megacities

This paper describes the use of the idea of natural time to propose a new method for characterizing the seismic risk to the world's major cities at risk of earthquakes. Rather than focus on forecasting, which is the computation of probabilities of future events, we define the term seismic nowcasting, which is the computation of the current state of seismic hazard in a defined geographic region.Comment: 9 Figures, 4 Table

Near mean-field behavior in the generalized Burridge-Knopoff earthquake model with variable range stress transfer

Simple models of earthquake faults are important for understanding the mechanisms for their observed behavior in nature, such as Gutenberg-Richter scaling. Because of the importance of long-range interactions in an elastic medium, we generalize the Burridge-Knopoff slider-block model to include variable range stress transfer. We find that the Burridge-Knopoff model with long-range stress transfer exhibits qualitatively different behavior than the corresponding long-range cellular automata models and the usual Burridge-Knopoff model with nearest-neighbor stress transfer, depending on how quickly the friction force weakens with increasing velocity. Extensive simulations of quasiperiodic characteristic events, mode-switching phenomena, ergodicity, and waiting-time distributions are also discussed. Our results are consistent with the existence of a mean-field critical point and have important implications for our understanding of earthquakes and other driven dissipative systems.Comment: 24 pages 12 figures, revised version for Phys. Rev.

A damage model based on failure threshold weakening

A variety of studies have modeled the physics of material deformation and damage as examples of generalized phase transitions, involving either critical phenomena or spinodal nucleation. Here we study a model for frictional sliding with long range interactions and recurrent damage that is parameterized by a process of damage and partial healing during sliding. We introduce a failure threshold weakening parameter into the cellular-automaton slider-block model which allows blocks to fail at a reduced failure threshold for all subsequent failures during an event. We show that a critical point is reached beyond which the probability of a system-wide event scales with this weakening parameter. We provide a mapping to the percolation transition, and show that the values of the scaling exponents approach the values for mean-field percolation (spinodal nucleation) as lattice size $L$ is increased for fixed $R$. We also examine the effect of the weakening parameter on the frequency-magnitude scaling relationship and the ergodic behavior of the model