# 3D Finite Element Modeling of fault-slip triggering caused by porepressure

Save this PDF as:

Size: px
Start display at page:

## Transcription

1 3D Finite Element Modeling of fault-slip triggering caused by porepressure changes Arsalan Sattari and David W. Eaton Department of Geoscience, University of Calgary Suary We present a 3D model using a widely used finite-element package, ABAQUS, to simulate activation of slip on a pre-existing fault due to an imposed perturbation in pore pressure on a subregion of the fault. The model inputs are 3 initial principal background stresses, elastic properties of the medium, coefficient of friction and geometry of the fault. Initial stresses are assigned according to Mohr Coulomb theory based on the assumption that the crust is in a critically stressed state. The calculation is shown to be stable, robust and consistent with Okada's analytical displacement calculation. Our results indicate that the fault rupture area is generally larger than the initially perturbed area. Moreover, the rupture area, and thus the earthquake magnitude, increases with the size of the perturbed area, irrespective of the overall fault dimensions. For example, our model predicts that a 2km by 2km area perturbed by MPa pore pressure results in a magnitude 4.4 earthquake with 1 MPa (1 bar) of shear stress drop. For a given fault geometry and initial state, we interpret our calculations to represent an upper limit for earthquake magnitude, as the equivalent stress release could be achieved with a sequence of smaller events. Introduction Fluid injection can trigger earthquake activity on a nearby fault, either by increasing pore pressure or by poroelastic stress coupling between fractured lithologies and formation fluids (Ellsworth, 213; Davies et al., 213). A pore-pressure increase tends to reduce the effective normal stress on the fault plane and thus stimulates slip (Healy, 1968). Zoback (212) showed that the majority of faults are subject to failure in response to a small perturbation in stress, relative to ambient stress levels. Keranen et al (213) showed that pore pressure can travel tens of kilometers away from the injection site and can trigger slip on a fault with a pore-pressure perturbation as small as.7 MPa. Although the majority of injectioninduced earthquakes are small in magnitude, the relatively high rates of induced seismicity significantly adds to seismic hazard in some areas (Ellsworth, 213). Knowing the potential maximum magnitude of induced-earthquake is useful in seismic hazard assessment (McGarr, 214); however it is highly dependent on factors including fluid-injection parameters, rheology and size of the fault. This study is motivated by the need to understand the dependency of maximum magnitude of induced earthquake to these factors. Theory In this study we aim to investigate the maximum earthquake magnitude that can be triggered by an increase of pore pressure around a fault. The examples presented here are confined to a strike-slip regime, but our method is applicable to any stress regime and fault geometry. According to Mohr- Coulomb failure criterion, fault slip occurs when the Mohr circle intersects or crosses the failure line GeoConvention 21: New Horizons 1

2 Normal Stress (Figure 1), where the slope of failure line is equal to coefficient of friction. Individual points on the Mohr circle represent the shear and normal stress for different fault orientations. A pore pressure (P p ) increase shifts the Mohr circle to left, and at the point of failure in critically stressed regime the maximum principal stress (σ max ) is given by (Zoback, 21) σ max = [ μ μ] 2 σ min, where σ min is the minimum principal stress and μ is the coefficient of friction. For μ =.6 the orientation of critically stressed fault deviates by 3 degrees from the maximum stress direction. The vertical background stress S v is defined by the pressure by the rocks lying above a particular depth. It is described by (Twiss and Moores, 27) as S v = ρ Layer g where ρ Layer is density and g is gravitational acceleration. In this study, we use a stress state in which the vertical principal stress is the average of minimum and maximum principal horizontal stresses. Hence, Methodology S h = dz 2 1+[ μ 2 +1+μ] 2 S v. The finite element modeling in this study is calculated using ABAQUS. A 3D cubical model is used, covering a volume of 9km depth, 7km length and 9 km width. A 3km wide fault with a depth extent of 3 km is imposed in the middle of the model, starting at 3km depth. The coefficient of friction is set to.6. Young s modulus and Poisson s ratio are 8 MPa and.33, respectively. An average density of 26 kg m 3 is assumed for the entire model. All sides of the model except the surface have no degrees of freedom to move. At the top of the model, a free-surface boundary condition is applied. Hexahedral elements with full integration points (C3D8) are used. The dimensions of elements around the fault are about 2 m increasing to ~2 m at the sides. The average horizontal and vertical principal stresses are calculated on each element subject to the condition that the fault is critically stressed. Principal stresses are perturbed by a pore pressure increase in the target area. The first step in the procedure is to fix all the nodes and compute a reaction-force solution on each node, which equilibrates the initial stress condition. In the second step fault nodes are released to calculate the equilibrium point and slip on the fault. Convergence is validated by verifying a unique response subject to decreasing mesh size. Failure Line µ=.6 Figure 1 The blue Mohr circle is critically stressed at the tangent point with the failure line which makes a 6 degree angle with the shear stress axis. The black Mohr circle is shifted to the left to represent an increases in pore pressure. The part of black Mohr circle which is above the red line is in the unstable regime and is subject to failure. 6 Smin Pore Pressure Smax GeoConvention 21: New Horizons 2

3 Slip Slip Slip a b c d e Figure 2 a), b) and c) Calculated slip contours in on a 3 km by 3 km fault plane. The initially perturbed region is indicated by the blue square. d), e) and f) Shear stress contours. Faults are located from 3km to 6km deep. f Results Figure 2 shows slip pattern and shear stress contours for varying size of initially perturbed area (Table 1) on a 3km by 3km fault plane. Pore pressure perturbed volume is cubical and surrounds the fault plane syetrically. The average slip increases proportionally to the perturbed volume. The slip is maximum in the middle of the pressure perturbed area and shrinks to zero as it propagates away. The rupture area is not bounded by the initially perturbed part of the fault and can be slightly larger. Asyetry is displayed in slip and displacement pattern, which is due to changes in initial stress distribution applied by ABAQUS to reach equilibrium in the first step of the calculation. Asyetry may also arise from the use of elements with different sizes around the fault. This issue can be resolved by using a smaller size of elements at the expense of a longer simulation time. A constant shear stress drop of ~1 MPa (1 bar) occurs for all three earthquakes, which is consistent with observed stress drops for induced earthquakes within our modelled depth range (e.g., Hough, 214). Okada (198) calculated displacement inside a semi-infinite medium due to slip on a finite rectangular fault. The inputs of the model are two elastic parameters of medium, geometry and orientation of fault, strike-slip, dip-slip and tensile components of dislocation. Figure 3 shows a comparison between the finite element Fault Size Pore Pressure Change Perturbed Area Average Slip Moment Magnitude Drop 3km * 3km MPa.km *.km 12 (fig.1.a) MPa (fig.1.d) 3km * 3km MPa 1km * 1km 19 (fig.1.b) MPa (fig.1.e) 3km * 3km MPa 2km *2km 2 (fig.1.c) MPa (fig.1.f) Table 1 Model parameters and calculated results. GeoConvention 21: New Horizons 3

4 Displacement Magnitude 2 Displacement Magnitude Finite Element Solution Solution by Okada Figure 3 Both picures show a map veiw of 9km by 9km horizontal surface at 4. km depth. The black line presents a fault with 3km width and 3km length, starts at 3km depth. Left) Displacement magnitude of nodes resulted from an average of 4 cm slip according to finite element simulation. Right) Displacement magnitude resulted from 4cm uniform slip according to analythical solution by Okada (198). modelling result and Okada s model, verifying that our finite element model estimation is consistent with the theoretical solution. Conclusions The 3D finite element model that we proposed here is capable of simulating maximum slip on a fault plane in any faulting regime, which is useful in seismic hazard modeling. Also it can simulate the maximum magnitude of an earthquake due to pore-pressure increase on a fault. The model results are in agreement with Okada s theoretical solution of displacements. We find that the rupture area can be larger than the initially perturbed region. Moreover, slip increases as the perturbed area becomes larger. A 2km by 2km area, disturbed by MPa pore pressure can trigger a 2 average slip on the fault. Seismic moment (M ) in dyne.cm is calculated by multiplying shear modulus (3 GPa), area of slip (~ 2.km by 2.km) and average slip. Therefore, moment magnitude 4.4 is calculated by using the formula given by (Hanks & Kanamori, 1979) Acknowledgments M = 2 3 logm 1.7. This work was supported by a Collaborative Research and Development Grant from the Natural Sciences and Engineering Research Council of Canada. TransAlta Utilities and Nanometrics Inc. are also thanked for their support of this initiative. References Davies, R., G. Fouler, A. Bindley and P. Styles (213), Induced seismicity and hydraulic fracturing for the recovery of hydrocarbons. Mar Petrol. Geol., 4, Ellsworth, W. L. (213), Injection-induced earthquakes. Science, 341 Frohlich, C., and M. Brunt (213), Two-year survey of earthquakes and injection/production wells in the Eagle Ford Shale, Texas, prior to the M October 211 earthquake. Earth Planet. Sci. Lett., 379, Hanks, T. C., and H. Kanamori (1979), A moment-magnitude scale. J. Geophys. Res, 84, GeoConvention 21: New Horizons 4

5 Healy, J. H., W. W. Rubey, D. T. Griggs, and C. B. Raleigh (1968), The Denver earthquakes. Science, 161, Hough, S. (214). Shaking from injection-induced earthquakes in the central and eastern United States. Bull. Seismol. Soc. Am., 14(), Keranen, K., H. Savage, G. Abers, and E. S. Cochran (213), Potentially induced earthquake in Oklahoma, USA: links between wastewater injection and the 211 M.7 earthquake sequence. Geology, 41, McGarr, A. (214) Maximum magnitude earthquake induced by fluid injection. J. Geophys. Res., 119, Okada, Y. (198) Surface deformation due to shear and tensile faults in a half-space. Bull. Seism. Soc. Am., 7, Twiss, R., and E. Moores (27), Structural geology. 2 ed.: W.H. Freeman and Company. Zoback, M.D and Gorelick, S.M, (212), Earthquake triggering and large-scale geologic storage of carbon dioxide. Proc. Natl. Acad. Sci. U.S.A., 19, Zoback, M. D., (21), Reservoir Geomechanics. Cambridge University Press. GeoConvention 21: New Horizons

### The Frictional Regime

The Frictional Regime Processes in Structural Geology & Tectonics Ben van der Pluijm WW Norton+Authors, unless noted otherwise 1/25/2016 10:08 AM We Discuss The Frictional Regime Processes of Brittle Deformation

### Effect Of The In-Situ Stress Field On Casing Failure *

Effect Of The In-Situ Stress Field On Casing Failure * Tang Bo Southwest Petroleum Institute, People's Republic of China Lian Zhanghua Southwest Petroleum Institute, People's Republic of China Abstract

### Brittle Deformation. Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm

Lecture 6 Brittle Deformation Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm WW Norton, unless noted otherwise Brittle deformation EarthStructure (2 nd

### Earthquake and Volcano Clustering at Mono Basin (California)

Excerpt from the Proceedings of the COMSOL Conference 2010 Paris Earthquake and Volcano Clustering at Mono Basin (California) D. La Marra *,1, A. Manconi 2,3 and M. Battaglia 1 1 Dept of Earth Sciences,

### The Mechanics of Earthquakes and Faulting

The Mechanics of Earthquakes and Faulting Christopher H. Scholz Lamont-Doherty Geological Observatory and Department of Earth and Environmental Sciences, Columbia University 2nd edition CAMBRIDGE UNIVERSITY

### This paper was prepared for presentation at the Unconventional Resources Technology Conference held in San Antonio, Texas, USA, 1-3 August 2016.

URTeC: 2461621 Determining Maximum Horizontal Stress With Microseismic Focal Mechanisms Case Studies in the Marcellus, Eagle Ford, Wolfcamp Alireza Agharazi*, MicroSeismic Inc. Copyright 2016, Unconventional

### Geology 229 Engineering Geology. Lecture 5. Engineering Properties of Rocks (West, Ch. 6)

Geology 229 Engineering Geology Lecture 5 Engineering Properties of Rocks (West, Ch. 6) Common mechanic properties: Density; Elastic properties: - elastic modulii Outline of this Lecture 1. Uniaxial rock

### Rate and State-Dependent Friction in Earthquake Simulation

Rate and State-Dependent Friction in Earthquake Simulation Zac Meadows UC Davis - Department of Physics Summer 2012 REU September 3, 2012 Abstract To better understand the spatial and temporal complexity

### Shear Stresses and Displacement for Strike-slip Dislocation in an Orthotropic Elastic Half-space with Rigid Surface

International Journal of Applied Science-Research and Review (IJAS) www.ijas.org.uk Original Article Shear Stresses and Displacement for Strike-slip Dislocation in an Orthotropic Elastic Half-space with

### 5. What is an earthquake 6. Indicate the approximate radius of the earth, inner core, and outer core.

Tutorial Problems 1. Where Do Earthquakes Happen? 2. Where do over 90% of earthquakes occur? 3. Why Do Earthquakes Happen? 4. What are the formulae for P and S velocity 5. What is an earthquake 6. Indicate

### Production-induced stress change in and above a reservoir pierced by two salt domes: A geomechanical model and its applications

Production-induced stress change in and above a reservoir pierced by two salt domes: A geomechanical model and its applications Peter Schutjens, Jeroen Snippe, Hassan Mahani, Jane Turner, Joel Ita and

### Daniel O Connell and Robert Creed. Fugro Consultants, Inc Cole Blvd., Suite 230, Lakewood CO

PROCEEDINGS, Thirty-Eighth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 11-13, 2013 SGP-TR-198 DETERMINING INJECTION-INDUCED EARTHQUAKE MAGNITUDE AND

### Originally published as:

Originally published as: Lorenzo Martín, F., Wang, R., Roth, F. (2002): The effect of input parameters on visco-elastic models of crustal deformation. - Física de la Tierra, 14, 33-54 The effect of input

### 21. Earthquakes I (p ; 306)

21. Earthquakes I (p. 296-303; 306) How many people have been killed by earthquakes in the last 4,000 years? How many people have been killed by earthquakes in the past century? What two recent earthquakes

### Modelling Progressive Failure with MPM

Modelling Progressive Failure with MPM A. Yerro, E. Alonso & N. Pinyol Department of Geotechnical Engineering and Geosciences, UPC, Barcelona, Spain ABSTRACT: In this work, the progressive failure phenomenon

### SHEAR-SLIP ANALYSIS IN MULTIPHASE FLUID-FLOW RESERVOIR ENGINEERING APPLICATIONS USING TOUGH-FLAC

PROCEEDINGS, TOUGH Symposium 6 Lawrence Berkeley National Laboratory, Berkeley, California, May 5 7, 6 SHEAR-SLIP ANALYSIS IN MULTIPHASE FLUID-FLOW RESERVOIR ENGINEERING APPLICATIONS USING TOUGH-FLAC Jonny

### Reservoir Geomechanics with ABAQUS

Reservoir Geomechanics with ABAQUS B. Bostrøm and E. Skomedal Statoil ASA, Norway Abstract: The coupled hydro-mechanical behavior of two North Sea high-pressure/hightemperature gas-condensate fields during

### SPE DISTINGUISHED LECTURER SERIES is funded principally through a grant of the SPE FOUNDATION

SPE DISTINGUISHED LECTURER SERIES is funded principally through a grant of the SPE FOUNDATION The Society gratefully acknowledges those companies that support the program by allowing their professionals

### Distinguished Lecturer David A. Ferrill. Sponsored by AAPG Foundation

Distinguished Lecturer David A. Ferrill Sponsored by AAPG Foundation AAPG Membership Information Benefits: Monthly Receipt of AAPG Explorer & AAPG Bulletin Discounts on Books and Other Publications Registration

### Analysis of Fracture Network Response to Fluid Injection

PROCEEDINGS, Fourtieth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 26-28, 2015 SGP-TR-204 Analysis of Fracture Network Response to Fluid Injection Moien

### Tectonic Seismogenic Index of Geothermal Reservoirs

Tectonic Seismogenic Index of Geothermal Reservoirs C. Dinske 1, F. Wenzel 2 and S.A. Shapiro 1 1 Freie Universität Berlin 2 KIT Karlsruhe November 27, 2012 Introduction M max reservoir location -2.0 Barnett

### A new methodology for estimating field horizontal stress from microseismic focal mechanisms

A new methodology for estimating field horizontal stress from microseismic focal mechanisms Alireza Agharazi 1*, Peter Duncan 1 and Michael Thornton 1 introduce a new methodology for the estimation of

Supplement A. Effective Coulomb Stress Analysis We assessed changes in Coulomb stress beneath reservoirs by the general theory of decoupled linear poroelasticity [Biot, 1956; Roeloffs, 1988; Wang 2000;

### High-resolution analysis of microseismicity related to hydraulic stimulations in the Berlín Geothermal Field, El Salvador

High-resolution analysis of microseismicity related to hydraulic stimulations in the Berlín Geothermal Field, El Salvador Grzegorz Kwiatek 1, Fatih Bulut 1 Marco Bohnhoff 1, Georg Dresen 1 1) GFZ Potsdam,

### INVESTIGATION OF GEOLOGICAL FAULT REACTIVATION AND OPENING

Mecánica Computacional Vol XXIX, págs. 8687-8697 (artículo completo) Eduardo Dvorkin, Marcela Goldschmit, Mario Storti (Eds.) Buenos Aires, Argentina, 15-18 Noviembre 2010 INVESTIGATION OF GEOLOGICAL FAULT

### friction friction a-b slow fast increases during sliding

µ increases during sliding faster sliding --> stronger fault --> slows sliding leads to stable slip: no earthquakes can start velocity-strengthening friction slow fast µ velocity-strengthening friction

### Methods of Interpreting Ground Stress Based on Underground Stress Measurements and Numerical Modelling

University of Wollongong Research Online Coal Operators' Conference Faculty of Engineering and Information Sciences 2006 Methods of Interpreting Ground Stress Based on Underground Stress Measurements and

### Chapter 6. Conclusions. 6.1 Conclusions and perspectives

Chapter 6 Conclusions 6.1 Conclusions and perspectives In this thesis an approach is presented for the in-situ characterization of rocks in terms of the distribution of hydraulic parameters (called SBRC

### Force and Stress. Processes in Structural Geology & Tectonics. Ben van der Pluijm. WW Norton+Authors, unless noted otherwise 1/9/ :35 PM

Force and Stress Processes in Structural Geology & Tectonics Ben van der Pluijm WW Norton+Authors, unless noted otherwise 1/9/2017 12:35 PM We Discuss Force and Stress Force and Units (Trigonometry) Newtonian

### Hitoshi Hirose (1), and Kazuro Hirahara (2) Abstract. Introduction

Three dimensional simulation for the earthquake cycle at a subduction zone based on a rate- and state-dependent friction law: Insight into a finiteness and a variety of dip-slip earthquakes Hitoshi Hirose

### Investigation of Injection-Triggered Slip on Basement Faults: Role of Fluid Leakoff on Post Shut-In Seismicity

Proceedings World Geothermal Congress 2015 Melbourne, Australia, 19-25 April 2015 Investigation of Injection-Triggered Slip on Basement Faults: Role of Fluid Leakoff on Post Shut-In Seismicity Jack Norbeck

### Search and Discovery Article # (2015) Posted April 20, 2015

Considering the Vertical Variation in Rock Mechanical Properties of a Lithologic Zone Using Laboratory Derived Data Is it Time for Geomechanical Stratigraphy?* Douglas E. Wyatt 1, Jesse Hampton 1, Dandan

### Estimating energy balance for hydraulic fracture stimulations: Lessons Learned from Basel

Estimating energy balance for hydraulic fracture stimulations: Lessons Learned from Basel David W. Eaton*, Department of Geoscience, University of Calgary, Calgary, Canada eatond@ucalgary.ca and Neda Boroumand,

### Targeting of Potential Geothermal Resources in the Great Basin from Regional Relationships between Geodetic Strain and Geological Structures

Targeting of Potential Geothermal Resources in the Great Basin from Regional Relationships between Geodetic Strain and Geological Structures Geoffrey Blewitt and Mark Coolbaugh Great Basin Center for Geothermal

### Elastic rebound theory

Elastic rebound theory Focus epicenter - wave propagation Dip-Slip Fault - Normal Normal Fault vertical motion due to tensional stress Hanging wall moves down, relative to the footwall Opal Mountain, Mojave

### Seismic stability safety evaluation of gravity dam with shear strength reduction method

Water Science and Engineering, 2009, 2(2): 52-60 doi:10.3882/j.issn.1674-2370.2009.02.006 http://kkb.hhu.edu.cn e-mail: wse@hhu.edu.cn Seismic stability safety evaluation of gravity dam with shear strength

### San Andreas Movie Can It Happen?

San Andreas Movie Can It Happen? Learning Objectives (LO) Lecture 14: Faults and Quakes Read: Chapter 10 and 11 Homework #12 due Thursday 12pm What we ll learn today:! 1. Compare strike-slip to dip-slip

### Earthquake Measurements

Volume II, Chapter 4: Appendix 4.A Appendix 4.A Earthquake Measurements 4.A.1. Earthquake Recording and Analysis Seismic waves radiated by earthquakes are recorded by networks of seismometers placed on

### The role of Coulomb stress changes for injection induced seismicity: The Basel enhanced geothermal system

GEOPHYSICAL RESEARCH LETTERS, VOL. 4, 72 77, doi:1.129/212gl54147, 213 The role of Coulomb stress changes for injection induced seismicity: The Basel enhanced geothermal system Flaminia Catalli, 1 Men

### Lecture 2: Deformation in the crust and the mantle. Read KK&V chapter 2.10

Lecture 2: Deformation in the crust and the mantle Read KK&V chapter 2.10 Tectonic plates What are the structure and composi1on of tectonic plates? Crust, mantle, and lithosphere Crust relatively light

### Implementing and Monitoring a Geomechanical Model for Wellbore Stability in an Area of High Geological Complexity

AADE-10-DF-HO-15 Implementing and Monitoring a Geomechanical Model for Wellbore Stability in an Area of High Geological Complexity J. G. Vargas, Halliburton; J. Mateus, Halliburton; A. Jaramillo, Halliburton;

### A comparison of seismicity rates and fluid injection operations in Oklahoma and California: Implications for crustal stresses

A comparison of seismicity rates and fluid injection operations in Oklahoma and California: Implications for crustal stresses Thomas Goebel California Institute of Technology Seismological Laboratory,

### FULL MOMENT TENSOR ANALYSIS USING FIRST MOTION DATA AT THE GEYSERS GEOTHERMAL FIELD

PROCEEDINGS, Thirty-Eighth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 11-13, 2013 SGP-TR-198 FULL MOMENT TENSOR ANALYSIS USING FIRST MOTION DATA AT

### Automatic Moment Tensor Analyses, In-Situ Stress Estimation and Temporal Stress Changes at The Geysers EGS Demonstration Project

PROCEEDINGS, 42nd Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 13-15, 2017 SGP-TR-212 Automatic Moment Tensor Analyses, In-Situ Stress Estimation and

### Earth Tides Can Trigger Shallow Thrust Fault Earthquakes

Earth Tides Can Trigger Shallow Thrust Fault Earthquakes Elizabeth S. Cochran, 1 * John E. Vidale, 1 Sachiko Tanaka 2 1 Department of Earth and Space Sciences and Institute of Geophysics and Planetary

### DYNAMIC ANALYSIS OF PILES IN SAND BASED ON SOIL-PILE INTERACTION

October 1-17,, Beijing, China DYNAMIC ANALYSIS OF PILES IN SAND BASED ON SOIL-PILE INTERACTION Mohammad M. Ahmadi 1 and Mahdi Ehsani 1 Assistant Professor, Dept. of Civil Engineering, Geotechnical Group,

### Modeling Fault Reactivation, Induced Seismicity, and Leakage during Underground CO2 Injection

Modeling Fault Reactivation, Induced Seismicity, and Leakage during Underground CO2 Injection Jonny Rutqvist A.P. Rinaldi 1, F. Cappa 1,2 1 Lawrence Berkeley National Laboratory, California, USA 2 Geoazur,

### The geomechanical significance of clay in geothermal reservoirs

The geomechanical significance of clay in geothermal reservoirs Carola Meller, Emmanuel Gaucher, Thomas Kohl KIT Institute of Applied Geosciences, Div. Geothermal Research, Adenauerring 20b, 76131 Karlsruhe,

### Mining. Slope stability analysis at highway BR-153 using numerical models. Mineração. Abstract. 1. Introduction

Mining Mineração http://dx.doi.org/10.1590/0370-44672015690040 Ricardo Hundelshaussen Rubio Engenheiro Industrial / Doutorando Universidade Federal do Rio Grande do Sul - UFRS Departamento de Engenharia

### Maximize the potential of seismic data in shale exploration and production Examples from the Barnett shale and the Eagle Ford shale

Maximize the potential of seismic data in shale exploration and production Examples from the Barnett shale and the Eagle Ford shale Joanne Wang, Paradigm Duane Dopkin, Paradigm Summary To improve the success

### Recent Research on EPS Geofoam Seismic Buffers. Richard J. Bathurst and Saman Zarnani GeoEngineering Centre at Queen s-rmc Canada

Recent Research on EPS Geofoam Seismic Buffers Richard J. Bathurst and Saman Zarnani GeoEngineering Centre at Queen s-rmc Canada What is a wall (SEISMIC) buffer? A compressible inclusion placed between

### Evaluation of Stress and Strain Induced by the Rock Compaction on a Hydrocarbon Well Completion Using Contact Interfaces with Abaqus

Evaluation of Stress and Strain Induced by the Rock Compaction on a Hydrocarbon Well Completion Using Contact Interfaces with Abaqus G. Capasso and G. Musso * eni E&P, S. Donato Milanese, MI, Italy * now

### GSA DATA REPOSITORY

GSA DATA REPOSITORY 2016167 SUPPLEMENTAL INFORMATION APPENDIX DR1 Barenblatt (1962) shows that the stress intensity factor K 0 for a fluid-filled crack of length l with a uniform fluid pressure P f,0 within

### EAS MIDTERM EXAM

Ave = 98/150, s.d. = 21 EAS 326-03 MIDTERM EXAM This exam is closed book and closed notes. It is worth 150 points; the value of each question is shown at the end of each question. At the end of the exam,

### Example-3. Title. Description. Cylindrical Hole in an Infinite Mohr-Coulomb Medium

Example-3 Title Cylindrical Hole in an Infinite Mohr-Coulomb Medium Description The problem concerns the determination of stresses and displacements for the case of a cylindrical hole in an infinite elasto-plastic

### GEOL5690: Fold and Thrust Belts and Orogenic Wedges ! = S 0

GEOL569: Fold and Thrust Belts and Orogenic Wedges One of the earlier mysteries in geology was the discovery in the 19 th century of large scale overthrusts in the Alps. Sheets of rock were found to have

### BE Semester- I ( ) Question Bank (MECHANICS OF SOLIDS)

BE Semester- I ( ) Question Bank (MECHANICS OF SOLIDS) All questions carry equal marks(10 marks) Q.1 (a) Write the SI units of following quantities and also mention whether it is scalar or vector: (i)

### Using information about wave amplitudes to learn about the earthquake size.

Earthquake Magnitudes and Moments Using information about wave amplitudes to learn about the earthquake size. Need to correct for decrease with distance M = log(a/t) + F(h,Δ) + C A is the amplitude of

### Hydraulic fracturing in unconventional shale gas reservoirs. Dr.-Ing. Johannes Will, Dynardo GmbH, Weimar, Germany

Hydraulic fracturing in unconventional shale gas reservoirs Dr.-Ing. Johannes Will, Dynardo GmbH, Weimar, Germany Founded: 2001 (Will, Bucher, CADFEM International) More than 35 employees, offices at Weimar

### Deformation of a layered half-space due to a very long tensile fault

Deformation of a layered half-space due to a very long tensile fault Sarva Jit Singh and Mahabir Singh Department of Mathematics, Maharshi Dayanand University, Rohtak 124 1, India. e-mail: s j singh@yahoo.com

### WELLBORE STABILITY ANALYSIS IN CHEMICALLY ACTIVE SHALE FORMATIONS

S911 WELLBORE SABILIY ANALYSIS IN CHEMICALLY ACIVE SHALE FORMAIONS by Xiang-Chao SHI *, Xu YANG, Ying-Feng MENG, and Gao LI State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Southwest

### on the earthquake's strength. The Richter scale is a rating of an earthquake s magnitude based on the size of the

Earthquakes and Seismic Waves An earthquake is the shaking and trembling that results from the movement of rock beneath Earth's surface. The point beneath Earth s surface where rock under stress breaks

### 3D ANALYSIS OF STRESSES AROUND AN UNLINED TUNNEL IN ROCK SUBJECTED TO HIGH HORIZONTAL STRESSES

3D ANALYSIS OF STRESSES AROUND AN UNLINED TUNNEL IN ROCK SUBJECTED TO HIGH HORIZONTAL STRESSES Abdel Meguid, M. Graduate Student, Department of Civil Engineering, University of Western Ontario, London,

### Summary. Introduction

1540160 Using Microseismicity to Understand Subsurface Fracture Systems and Increase the Effectiveness of Completions: Eagle Ford Shale, TX John P. Detring 1 and Sherilyn Williams-Stroud 2 1. Analysis,

### THE EFFECT OF DIRECTIVITY ON THE STRESS PARAMETER DETERMINED FROM GROUND MOTION OBSERVATIONS

Bulletin of the Seismological Society of America, Vol. 79, No. 6, pp. 1984-1988, December 1989 THE EFFECT OF DIRECTIVITY ON THE STRESS PARAMETER DETERMINED FROM GROUND MOTION OBSERVATIONS BY DAVID M. BOORE

### Elastic Rebound Theory

Earthquakes Elastic Rebound Theory Earthquakes occur when strain exceeds the strength of the rock and the rock fractures. The arrival of earthquakes waves is recorded by a seismograph. The amplitude of

### Implementation of geomechanical model to identify total loss zones

SPE Workshop OILFIELD GEOMECHANICS Slide 1 Implementation of geomechanical model to identify total loss zones Session 3. Geomechanics for drilling and completion Kalinin Oleg Co-authors: Zhigulsky S.V.;

### Characterizing and Responding to Seismic Risk Associated with Earthquakes Potentially Triggered by Saltwater Disposal and Hydraulic Fracturing

Characterizing and Responding to Seismic Risk Associated with Earthquakes Potentially Triggered by Saltwater Disposal and Hydraulic Fracturing R. J. Walters 1, M. D. Zoback 1, J. W. Baker 2 and G. C. Beroza

### Anisotropy of Shale Properties: A Multi-Scale and Multi-Physics Characterization

Observation Scale Wavelength 10 0 10 4 10 6 10 8 10 12 10 16 10 18 10 20 Frequency (Hz) Anisotropy of Shale Properties: A Multi-Scale and Multi-Physics Characterization Elastic, Mechanical, Petrophysical

### In-Situ Stress and Active Faulting. In Oklahoma. Richard C. Alt II and Mark D. Zoback. Department of Geophysics. Stanford University

1 2 In-Situ Stress and Active Faulting In Oklahoma 3 4 5 6 Richard C. Alt II and Mark D. Zoback Department of Geophysics Stanford University 7 8 9 10 11 12 13 14 15 16 Corresponding Author Richard C. Alt

### Section 19.1: Forces Within Earth Section 19.2: Seismic Waves and Earth s Interior Section 19.3: Measuring and Locating.

CH Earthquakes Section 19.1: Forces Within Earth Section 19.2: Seismic Waves and Earth s Interior Section 19.3: Measuring and Locating Earthquakes Section 19.4: Earthquakes and Society Section 19.1 Forces

### Lateral Earth Pressure

1 of 11 6/2/2012 4:28 AM Lateral Earth Pressure The magnitude of lateral earth pressure depends on: 1. Shear strength characteristics of soil 2. Lateral strain condition 3. Pore water pressure 4. State

### Caprock Integrity Analysis in Thermal Operations: An Integrated Geomechanics Approach

WHOC11-609 Caprock Integrity Analysis in Thermal Operations: An Integrated Geomechanics Approach Safdar Khan 1, Hongxue Han, Sajjad Ansari, Morteza Vishteh and Nader Khosravi Data and Consulting Services,

### Asish Karmakar 1, Sanjay Sen 2 1 (Corresponding author, Assistant Teacher, Udairampur Pallisree Sikshayatan (H.S.), Udairampur, P.O.

IOSR Journal of Applied Geology and Geophysics (IOSR-JAGG) e-issn: 3 99, p-issn: 3 98.Volume 4, Issue 5 Ver. III (Sep. - Oct. 6), PP 39-58 www.iosrjournals.org A Sudden Movement across an Inclined Surface

### Predicting Subseismic Fracture Density and Orientation in the Gorm Field, Danish North Sea: An Approach Using Elastic Dislocation Models*

Predicting Subseismic Fracture Density and Orientation in the Gorm Field, Danish North Sea: An Approach Using Elastic Dislocation Models* Brett Freeman 1, David Quinn 1, Cathal G. Dillon 1, Michael Arnhild

### Wastewater disposal and earthquake swarm activity at the southern end of the Central Valley, California

GEOPHYSICAL RESEARCH LETTERS, VOL.???, XXXX, DOI:10.1029/, Wastewater disposal and earthquake swarm activity at the southern end of the Central Valley, California T. H. W. Goebel, 1* S. M. Hosseini, 2

### Identification of natural fractures and in situ stress at Rantau Dedap geothermal field

IOP Conference Series: Earth and Environmental Science PAPER OPEN ACCESS Identification of natural fractures and in situ stress at Rantau Dedap geothermal field To cite this article: Andika Artyanto et

### Earthquakes and Seismotectonics Chapter 5

Earthquakes and Seismotectonics Chapter 5 What Creates Earthquakes? The term Earthquake is ambiguous: Applies to general shaking of the ground and to the source of the shaking We will talk about both,

### Influence of material contrast on fault branching behavior

GEOPHYSICAL RESEARCH LETTERS, VOL. 38,, doi:10.1029/2011gl047849, 2011 Influence of material contrast on fault branching behavior Nora DeDontney, 1 James R. Rice, 1,2 and Renata Dmowska 2 Received 20 April

### An embedded fracture modeling framework for simulation of hydraulic fracturing and shear stimulation

Comput Geosci (2016) 20:1 18 DOI 10.1007/s10596-015-9543-2 ORIGINAL PAPER An embedded fracture modeling framework for simulation of hydraulic fracturing and shear stimulation Jack H. Norbeck 1 Mark W.

### Rock Failure. Topics. Compressive Strength Rock Strength from Logs Polyaxial Strength Criteria Anisotropic Rock Strength Tensile Strength

Rock Failure Topics Compressive Strength Rock Strength from Logs Polyaxial Strength Criteria Anisotropic Rock Strength Tensile Strength Key Points 1. When rock fails in compression, the compressive stress

### ANALYSIS OF VERTICAL GROWTH OF FRACTURES IN FRAC PACK OPERATIONS IN RESERVOIR ROCKS

ANALYSIS OF VERTICAL GROWTH OF FRACTURES IN FRAC PACK OPERATIONS IN RESERVOIR ROCKS Paulo Dore Fernandes PETROBRAS S. A., Centro de Pesquisas - CENPES Rio de Janeiro, RJ Paulo Roberto Ribeiro Universidade

### Activity Submitted by Tim Schroeder, Bennington College,

Structural Analysis of a Hot Dry Rock Geothermal Energy System Activity Submitted by Tim Schroeder, Bennington College, tschroeder@bennington.edu Description: This project applies basic geologic skills

### Earthquakes Chapter 19

Earthquakes Chapter 19 Does not contain complete lecture notes. What is an earthquake An earthquake is the vibration of Earth produced by the rapid release of energy Energy released radiates in all directions

### Stress-Strain Behavior

Stress-Strain Behavior 6.3 A specimen of aluminum having a rectangular cross section 10 mm 1.7 mm (0.4 in. 0.5 in.) is pulled in tension with 35,500 N (8000 lb f ) force, producing only elastic deformation.

### Fracture and vein patterns as indicators of deformation history: a numerical study

Geological Society, London, Special Publications Fracture and vein patterns as indicators of deformation history: a numerical study Daniel Koehn, Jochen Arnold and Cees W. Passchier Geological Society,

### Dangerous tsunami threat off U.S. West Coast

Earthquakes Ch. 12 Dangerous tsunami threat off U.S. West Coast Earthquakes What is an Earthquake? It s the shaking and trembling of the Earth s crust due to plate movement. The plates move, rocks along

### Rheology. What is rheology? From the root work rheo- Current: flow. Greek: rhein, to flow (river) Like rheostat flow of current

Rheology What is rheology? From the root work rheo- Current: flow Greek: rhein, to flow (river) Like rheostat flow of current Rheology What physical properties control deformation? - Rock type - Temperature

### Global geophysics and wave propagation

Global geophysics and wave propagation Reading: Fowler p76 83 Remote sensing Geophysical methods Seismology Gravity and bathymetry Magnetics Heat flow Seismology: Directly samples the physical properties

### Lecture # 6. Geological Structures

1 Lecture # 6 Geological Structures ( Folds, Faults and Joints) Instructor: Dr. Attaullah Shah Department of Civil Engineering Swedish College of Engineering and Technology-Wah Cantt. 2 The wavy undulations

### Apply Rock Mechanics in Reservoir Characterization Msc Julio W. Poquioma

Apply Rock Mechanics in Reservoir Characterization Msc Julio W. Poquioma Chief Reservoir Engineer PetroSA Mobil: +27 79 492.02.52 Email julio.poquioma@petrosa.co.za Aberdeen, April 18 20 / 2012 Content

### STABILITY CHECK AND SUPPORT DESIGNING FOR THE GR-2011 EXPLORATION DRIFT

UNDERGROUND MINING ENGINEERING 19 (2011) 83-91 UDK 62 FACULTY OF MINING AND GEOLOGY, BELGRADE YU ISSN 03542904 Professional paper STABILITY CHECK AND SUPPORT DESIGNING FOR THE GR-2011 EXPLORATION DRIFT

### IZMIT BAY BRIDGE SOUTH APPROACH VIADUCT: SEISMIC DESIGN NEXT TO THE NORTH ANATOLIAN FAULT

Istanbul Bridge Conference August 11-13, 2014 Istanbul, Turkey IZMIT BAY BRIDGE SOUTH APPROACH VIADUCT: SEISMIC DESIGN NEXT TO THE NORTH ANATOLIAN FAULT A. Giannakou 1, J. Chacko 2 and W. Chen 3 ABSTRACT

### Terminology for fluid-injection induced seismicity in oil and gas operations

FOCUS ARTICLE Coordinated by Mostafa Naghizadeh / Fereidoon Vasheghani Terminology for fluid-injection induced seismicity in oil and gas operations David W. Eaton1, Mirko van der Baan2 and Allan Ingelson3

### Main Means of Rock Stress Measurement

Real Stress Distributions through Sedimentary Strata and Implications for Reservoir Development and Potential Gas and Coal Development Strategies Ian Gray Sigra Pty Ltd 93 Colebard St West, Acacia Ridge,