Defmod, an earthquake simulator that adaptively switches between quasi-static and dynamic states
|
|
- Joel Carter
- 5 years ago
- Views:
Transcription
1 Defmod, an earthquake simulator that adaptively switches between quasi-static and dynamic states C. Meng B. Hager ERL, MIT May 16, 2016
2 Reservoir production/injection induced seismicity Groningen earthquakes, gas production induced. Oklahoma earthquakes, waste water disposal induced.
3 ground overburden reservoir producer reservoir underburden
4 Defmod code is supposed to Capture poro-visco-elasticity processes due to fluid injection/production, viscoelastic flow and external loadings; When failure criterion is met, allow the fault to have frictional slip, accompanied by seismic radiation; Exchange the fault slip and stress perturbation between the quasi-static and dynamic solvers, forming a hybrid solver. History match the earthquake event occurrence and waveforms.
5 Hybrid solver flowchart The hybrid model will return to the quasi-static (parent) loop once the dynamic run is over.
6 Dynamic fault slip schematic Nucleation starts at a fault patch. n f 0 -f 0 t
7 Dynamic fault slip schematic Slip and stress concentration. n f 1 -f 1 t
8 Dynamic fault slip schematic Rupture propagation. n f 2 -f 2 t
9 Dynamic fault slip schematic More slip. n f 3 -f 3 t
10 Dynamic fault slip schematic Shear stress drop and rupture arrest. n f 4 -f 4 f 4 ' -f 4 ' t
11 Fault constraint Coinciding nodes on the fault belong to different elements.
12 For locked and permeable fault, n x n z 0 n x n z 0 t x t z 0 t x t z For slipping and permeable fault, [ ] nx n z 0 n x n z u (1) x u (1) z p (1) u (3) x u (3) z p (3) u (1) x u (1) z p (1) u (3) x u (3) z p (3) = 0, = 0.
13 For locked and impermeable fault, [ ] nx n z n x n z t x t z t x t z For slipping and impermeable fault, [ ] nx n z n x n z u (1) x u (1) z u (3) x u (3) z u (1) x u (1) z u (3) x u (3) z = 0, = 0, A compact linear constraint equation can be written as GU = I, where, G is constraint matrix, I is constraint function, nonzero for prescribed slip (pressure jump).
14 where, governing equation { KU = F Mü + C u + Ku = f (quasi-)static, dynamic, K, M and C are stiffness, mass and damping matrices; [ ] u U is the solutions vector, e.g. U = displacement and p pressure for poroelastic problem; [ ] f F is exterior load, e.g. F =, force and flow rate. q
15 Write the time dependent quasi-static equation in the compact form: KU n = F n. The constrained equation can be written as [ ] [ ] [ ] K G T Un Fn =, G 0 λ n I n n is the time step index, and λ n is Lagrange Multiplier, i.e. the nodal force/flux needed to let the solution honor the constraint function I n. Solution is solved by inverting the global matrix.
16 With the unconstrained solution obtained from the Newmark method, ( ) u n = M 1 ( t 2 f n Ku n 1 ) tc (u n 1 u n 2 ) +2u n 1 u n 2. The Forward Increment Lagrange Multiplier method [Carpenter et al. 1991] is applied to impose the constraint, ( ) λ n = t 2 GM 1 G T 1 (Gun I n ), u n =u n t 2 M 1 G T λ n.
17 Function and benchmark functionality method benchmark implicit, poroelasticity pore pressure Mandel stabilization viscoelastic implicit Abaqus power law (quasi)static Lagrange Mohr-Coulomb constraint Multiplier elastodynamic forward increment constraint Lagrange Multiplier SCEC absorbing viscous damping fault/faulting implicit/explicit
18 SCEC benchmark SCEC benchmark problem TPV5: Strike-slip rupture on a heterogeneously stressed fault patch. z km 3 4 3km x y
19 Fault slip rate magnitude [m/s]
20 Rupture front comparison against EqSim (Pylith) and SGFD Defmod EqSim SGFD z [km] y [km]
21 Waveform comparison for station 1. v [m/s] Defmod 0.4 EqSim SGFD t [s]
22 3D quasi-static loading and fault rupture Hybrid output: seismic radiation, velocity magnitude [m/s].
23 Hybrid output: quasi-static state after failure, x displacement [m] discontinuous across fault.
24 2D production induced earthquake Quasi-static output: fault traction and pressure. Production/shutin is indicated by pressure perturbation. Contact force/traction always appears in symmetric pairs.
25 Induced earthquake: dynamic output for a small event. 2D plot of the velocity magnitude. 1D (fault profile) plot of the slip rate.
26 Induced earthquake: dynamic output for a large event. 2D plot of the velocity magnitude. 1D (fault profile) plot of the slip rate.
27 Summary The model adaptively switches between (quasi-)static and dynamic states describing the induced earthquake cycles. The model functionalities are well benchmarked against established results. Future work The model will be used for history matching the induced earthquake occurrence and waveforms in Groningen. Once a reasonable history match is achieved, the model is capable of predicting induced seismic risks for given scenarios of production/injection.
28 Code availability GNU General Public License, for bug report and contribution: Stable code: Developer version: https: //bitbucket.org/chunfangmeng/defmod-dev
29 Poroelastic equation The production/injection is typically of the time scale much longer than the seismic events. Such process can be considered quasi-static. The incremental loading scheme for poroelasticity [Smith and Griffiths, 2004] with the source and storage terms for the fluid: [ ] [ ] [ ] Ke H un f H T = n tk c + S p n q n tk c p n 1 Where, K e and K c are the solid and fluid stiffness matrices; H is the coupling matrix; S is the storage matrix, f and q are the nodal force increment and the fluid flux during one time step t. Absolute solution at step n: [ ] [ ] un 1 un U n = + p n p n 1 where, u n 1 and p n 1 are the pressure of previous step.
30 Stabilizing fluid pressure Unstable pressure is caused by using linear element known as the Ladyzenskaja-Babuska-Brezzi restrictions. Pore pressure changes, 2D triangular element domain, following co-seismic slip on a thrust fault, with (left) and without (right) stabilization.
31 Local pressure projection scheme [Bochev and Dohrmann, 2006] is implemented to stabilize the pore pressure, [White and Borja, 2008] for quad/hex element. [ ] 0 F n+1 =F n+1, H s p n [ ] 0 0 K n+1 =K n+1 +, 0 H s, where, H s = Ω (N (1/n e )I N ) T (N (1/n e )I N )/(2G)dΩ, where, N is the shape function; ( ) is a function averaged over all nodal points of each element; G is the shear modulus; τ scale is a scale constant.
32 Poroelastic benchmark: Mandel solution [M. Kurashige et al. 2004] Mandel Cartesian Mandel cylindrical
33
34 1.2 normalized pressure at different time 1.4 normalized pressure t = t = D Mandel 2D Defmod x/ a 0.2 3D Mandel 3D Defmod r/ a
35 Pressure at origin, 2D vs 3D
36 Static fault model benchmark Anisotropic loading model: Incrementally load the sample at the rate of σ H = σ V = 0.2 MPa/yr. At year 50, stop the loading increment in the horizontal direction, and keep the increment in the vertical direction until year 150.
37 Stress ratio τ/σ n calculated by Defmod at different depth and time against the analytical values yr yr τ/σ n yr LM analytical 1.5 t=0 2.0 z [km]
38 Viscoelastic problem The stiffness matrix and RHS vector of a viscoelastic media have [Melosh and Raefsky 1980] K n+1 = B T (D 1 + α tβ n) 1 BdΩ Ω, F n+1 = B T (D 1 + α tβ n) 1 ( tβ n )dω + F n+1 Ω where, B is displacement to strain matrix depending on the element geometry, D is the element stiffness matrix depending on the elastic constants, β(σ) = σe 1 4η C c : σ, e 1, σ = (σ ii σ jj ) 2 /(2d) + σ 2 ij, d = 2 or 3, i j
39 Viscoelastics: β and β β = β(σ) = σe 1 4η C c : σ C c = /3 2/3 2/3 2/3 4/3 2/3 2/3 2/3 4/3 0, 2D, , 3D. σ e 1 c 1 c 1 c 3 c 1 c 1 c 3 4η c 3 c 3 4c 2 4/3 + S 2 x 2/3 + S x S y 2/3 + S x S z S x T 1 S x T 2 S x T 3 2/3 + S x S y 4/3 + S 2 σ e 1 y 2/3 + S y S z S y T 1 S y T 2 S y T 3 2/3 + S x S z 2/3 + S y S z 4/3 + S 2 z S z T 1 S z T 2 S z T 3 4η S x T 1 S y T 1 S z T T 2 1 T 1 T 2 T 1 T 3 S x T 2 S y T 2 S z T 2 T 2 T T 2 2 T 2 T 3 S x T 3 S y T 3 S z T 3 T 3 T 1 T 3 T T 2 3, 2D,, 3D.
40 c 1 =1 + (e 1)((σ xx σ yy )/(2σ)) 2 c 2 =1 + (e 1)(σ xy /σ) 2 c 3 =(e 1)(σ xx σ yy σ yy σ xy )/σ 2 S x =c(2σ xx σ yy σ zz )/(3σ) where c = e 1 S y =c(2σ yy σ zz σ xx )/(3σ) S z =c(2σ zz σ xx σ yy )/(3σ) T 1 =2cσ xy /σ, T 2 = 2cσ yz /σ, T 3 = 2cσ xz /σ
41 When e = 1, β = 1 4η C c, the matrix K n is then independent of σ ij and n, as long as the step length t is a constant. In this case we only need to assemble the viscoelastic stiffness matrix K once for the first step, and keep updating the RHS F n. When e > 1 however, both the stiffness matrix and RHS need to be reassembled for every step. Since the scale factor α and the effective viscosity η always appear as α 4η, they are treated as one parameter. Therefore, in addition to the elastic constants, two parameters, e and η, need to be specified.
42 Viscoelastic benchmark Slip on a strike-slip fault: Elastic crust over viscoelastic mantle. Only a part of model domain is shown.
43 Comparison against Abaqus at t = 0 and t = 10 years. The displacement is plotted along a trajectory perpendicular to the fault plane through the viscoelastic layer. Displacement (m) Abaqus Defmod T=10 years 0.0 T=0 years Distance (km)
44 Explicit solver for the seismic radiation Mü + C u + Ku = f, where, M is mass matrix; C = αm + βk is damping matrix; α and β are Rayleigh damping coefficients. Newmark explicit scheme has, ( ) u n = M 1 ( t 2 f n Ku n 1 ) tc (u n 1 u n 2 ) +2u n 1 u n 2 The incremental form, u n =M 1 ( ( t 2 f n K u n 1 ) tc ( u n 1 u n 2 )) + 2 u n 1 u n 2
45 Absorbing boundary [Lysmer and Kuhlemeyer, 1969] propose that the absorbing boundary for the elastodynamic model can be achieved by adding additional terms to the damping matrix as, { c ii + c ii = Γ ρv pdγ, p wave i th axis, c ii + Γ ρv sdγ, p wave i th for all i-axes. axis, Assuming small incident angles (θ < 30 ), that the p-wave can be considered roughly perpendicular to the absorbing boundary. Therefore, this addition is irrespective to the coming wave directions.
46 TVP 5 wave forms Waveform comparison for station 2. v [m/s] Defmod EqSim SGFD t [s]
47 TVP 5 wave forms Waveform comparison for station 3. v [m/s] Defmod 0.2 EqSim SGFD t [s]
48 TVP 5 wave forms Waveform comparison for station 4. v [m/s] Defmod EqSim SGFD t [s]
49 Performance Subduction, 8 m elements, implicitly solved, under 120 secs. Thrust fault, 6 m elements, 2000 explicit steps, under 60 secs.
50 Winkler Foundation Winkler Foundation is to consider the gravity caused anisotropy by adding linear springs in the gravity direction. The element stiffness matrix is therefore modified by k (p) ii = k (p) ii + 1 n p Γ g ρgdγ g, for p Γ g, i th axis g, Γ g g, where, n p is the number of Gauss points on the facet Γ g.
Multiphysics and multiscale earthqauke modeling
Multiphysics and multiscale earthqauke modeling C. Meng (cmeng@mit.edu) B. Hager ERL, MIT May 29, 2017 Multiphysics Poro(visco-)elasticity, (quasi-)static and dynamic processes. ground overburden reservoir
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Overview Milestones in continuum mechanics Concepts of modulus and stiffness. Stress-strain relations Elasticity Surface and body
More informationData Repository Hampel et al., page 1/5
GSA DATA REPOSITORY 2138 Data Repositor Hampel et al., page 1/5 SETUP OF THE FINITE-ELEMENT MODEL The finite-element models were created with the software ABAQUS and consist of a 1-km-thick lithosphere,
More informationStress equilibrium in southern California from Maxwell stress function models fit to both earthquake data and a quasi-static dynamic simulation
Stress equilibrium in southern California from Maxwell stress function models fit to both earthquake data and a quasi-static dynamic simulation Peter Bird Dept. of Earth, Planetary, and Space Sciences
More informationA new hybrid numerical scheme for modelling elastodynamics in unbounded media with near-source heterogeneities
A new hybrid numerical scheme for modelling elastodynamics in unbounded media with near-source heterogeneities Setare Hajarolasvadi Ahmed E. Elbanna https://www.youtube.com/watch?v=bltx92tuwha MOTIVATION
More information1 Exercise: Linear, incompressible Stokes flow with FE
Figure 1: Pressure and velocity solution for a sinking, fluid slab impinging on viscosity contrast problem. 1 Exercise: Linear, incompressible Stokes flow with FE Reading Hughes (2000), sec. 4.2-4.4 Dabrowski
More informationElizabeth H. Hearn modified from W. Behr
Reconciling postseismic and interseismic surface deformation around strike-slip faults: Earthquake-cycle models with finite ruptures and viscous shear zones Elizabeth H. Hearn hearn.liz@gmail.com modified
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationFault Representation Methods for Spontaneous Dynamic Rupture Simulation. Luis A. Dalguer
Fault Representation Methods for Spontaneous Dynamic Rupture Simulation Luis A. Dalguer Computational Seismology Group Swiss Seismological Service (SED) ETH-Zurich July 12-18, 2011 2 st QUEST Workshop,
More information3D Finite Element Modeling of fault-slip triggering caused by porepressure
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
More informationInstabilities and Dynamic Rupture in a Frictional Interface
Instabilities and Dynamic Rupture in a Frictional Interface Laurent BAILLET LGIT (Laboratoire de Géophysique Interne et Tectonophysique) Grenoble France laurent.baillet@ujf-grenoble.fr http://www-lgit.obs.ujf-grenoble.fr/users/lbaillet/
More informationMathematical Modelling of a Fault Slip Induced by Water Injection
Mathematical Modelling of a Fault Slip Induced by Water Injection T. S. Nguyen, 1 J.Rutqvist 2 and Y. Gugliemi 2 1 Canadian Nuclear Safety Commission 2 Lawrence Berkeley National Laboratory ComGeo IV Symposium
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Standard Solids and Fracture Fluids: Mechanical, Chemical Effects Effective Stress Dilatancy Hardening and Stability Mead, 1925
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Surface and body forces Tensors, Mohr circles. Theoretical strength of materials Defects Stress concentrations Griffith failure
More informationSpectral Element simulation of rupture dynamics
Spectral Element simulation of rupture dynamics J.-P. Vilotte & G. Festa Department of Seismology, Institut de Physique du Globe de Paris, 75252 France ABSTRACT Numerical modeling is an important tool,
More informationOverview of the PyLith Finite Element Code
Overview of the PyLith Finite Element Code www.geodynamics.org Charles Williams (GNS Science) Brad Aagaard (USGS) Matt Knepley (University of Chicago) CIG Software PyLith is a Tool for Crustal Deformation
More informationSurface changes caused by erosion and sedimentation were treated by solving: (2)
GSA DATA REPOSITORY 214279 GUY SIMPSON Model with dynamic faulting and surface processes The model used for the simulations reported in Figures 1-3 of the main text is based on two dimensional (plane strain)
More informationPEAT SEISMOLOGY Lecture 12: Earthquake source mechanisms and radiation patterns II
PEAT8002 - SEISMOLOGY Lecture 12: Earthquake source mechanisms and radiation patterns II Nick Rawlinson Research School of Earth Sciences Australian National University Waveform modelling P-wave first-motions
More informationEarthquake and Volcano Deformation
Earthquake and Volcano Deformation Paul Segall Stanford University Draft Copy September, 2005 Last Updated Sept, 2008 COPYRIGHT NOTICE: To be published by Princeton University Press and copyrighted, c
More informationNumerical modeling of dynamic rupture propagation
IX. Slovenská geofyzikálna konferencia, 22. - 23. jún 2011 Fakulta matematiky, fyziky a informatiky UK, Bratislava Numerical modeling of dynamic rupture propagation Martin Galis Peter Moczo Jozef Kristek
More informationReview of Strain Energy Methods and Introduction to Stiffness Matrix Methods of Structural Analysis
uke University epartment of Civil and Environmental Engineering CEE 42L. Matrix Structural Analysis Henri P. Gavin Fall, 22 Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods
More informationDynamic Analysis Contents - 1
Dynamic Analysis Contents - 1 TABLE OF CONTENTS 1 DYNAMIC ANALYSIS 1.1 Overview... 1-1 1.2 Relation to Equivalent-Linear Methods... 1-2 1.2.1 Characteristics of the Equivalent-Linear Method... 1-2 1.2.2
More informationEDEM DISCRETIZATION (Phase II) Normal Direction Structure Idealization Tangential Direction Pore spring Contact spring SPRING TYPES Inner edge Inner d
Institute of Industrial Science, University of Tokyo Bulletin of ERS, No. 48 (5) A TWO-PHASE SIMPLIFIED COLLAPSE ANALYSIS OF RC BUILDINGS PHASE : SPRING NETWORK PHASE Shanthanu RAJASEKHARAN, Muneyoshi
More informationSeismic and aseismic processes in elastodynamic simulations of spontaneous fault slip
Seismic and aseismic processes in elastodynamic simulations of spontaneous fault slip Most earthquake simulations study either one large seismic event with full inertial effects or long-term slip history
More informationModeling seismic wave propagation during fluid injection in a fractured network: Effects of pore fluid pressure on time-lapse seismic signatures
Modeling seismic wave propagation during fluid injection in a fractured network: Effects of pore fluid pressure on time-lapse seismic signatures ENRU LIU, SERAFEIM VLASTOS, and XIANG-YANG LI, Edinburgh
More informationFinite element modelling of fault stress triggering due to hydraulic fracturing
Finite element modelling of fault stress triggering due to hydraulic fracturing Arsalan, Sattari and David, Eaton University of Calgary, Geoscience Department Summary In this study we aim to model fault
More informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
More informationBrittle 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
More informationDETAILS ABOUT THE TECHNIQUE. We use a global mantle convection model (Bunge et al., 1997) in conjunction with a
DETAILS ABOUT THE TECHNIQUE We use a global mantle convection model (Bunge et al., 1997) in conjunction with a global model of the lithosphere (Kong and Bird, 1995) to compute plate motions consistent
More informationIntroduction to Seismology Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 12.510 Introduction to Seismology Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Stress and Strain
More informationDevelopment of a Predictive Simulation System for Crustal Activities in and around Japan - II
Development of a Predictive Simulation System for Crustal Activities in and around Japan - II Project Representative Mitsuhiro Matsu'ura Graduate School of Science, The University of Tokyo Authors Mitsuhiro
More informationExercise: concepts from chapter 8
Reading: Fundamentals of Structural Geology, Ch 8 1) The following exercises explore elementary concepts associated with a linear elastic material that is isotropic and homogeneous with respect to elastic
More informationEarthquake response analysis of rock-fall models by discontinuous deformation analysis
c Earthquake response analysis of rock-fall models by discontinuous deformation analysis T. Sasaki, I. Hagiwara & K. Sasaki Rock Engineering Lab., Suncoh Consultants Co. Ltd., Tokyo, Japan R. Yoshinaka
More informationWeight-adjusted DG methods for elastic wave propagation in arbitrary heterogeneous media
Weight-adjusted DG methods for elastic wave propagation in arbitrary heterogeneous media Jesse Chan Department of Computational and Applied Math, Rice University ICOSAHOM 2018 July 12, 2018 Chan (CAAM)
More informationA Study on Numerical Solution to the Incompressible Navier-Stokes Equation
A Study on Numerical Solution to the Incompressible Navier-Stokes Equation Zipeng Zhao May 2014 1 Introduction 1.1 Motivation One of the most important applications of finite differences lies in the field
More informationDynamic Earthquake Triggering Due to Stress from Surface Wave Particle Displacement
Emily Morton Geop 523: Theoretical Seismology 29 April 2011 Introduction Dynamic Earthquake Triggering Due to Stress from Surface Wave Particle Displacement Earthquakes can be triggered by other earthquakes.
More informationThe Finite Element Method for Mechonics of Solids with ANSYS Applicotions
The Finite Element Method for Mechonics of Solids with ANSYS Applicotions ELLIS H. DILL 0~~F~~~~"P Boca Raton London New Vork CRC Press is an imprint 01 the Taylor & Francis Group, an Informa business
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting Lecture 20, 30 Nov. 2017 www.geosc.psu.edu/courses/geosc508 Seismic Spectra & Earthquake Scaling laws. Seismic Spectra & Earthquake Scaling laws. Aki, Scaling law
More informationGeology for Engineers Rock Mechanics and Deformation of Earth Materials
89.325 Geology for Engineers Rock Mechanics and Deformation of Earth Materials Why do rocks break? Rock mechanics experiments a first order understanding. Faults and Fractures Triaxial load machine. a)
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting Lectures & 3, 9/31 Aug 017 www.geosc.psu.edu/courses/geosc508 Discussion of Handin, JGR, 1969 and Chapter 1 Scholz, 00. Stress analysis and Mohr Circles Coulomb Failure
More informationStatic & Dynamic. Analysis of Structures. Edward L.Wilson. University of California, Berkeley. Fourth Edition. Professor Emeritus of Civil Engineering
Static & Dynamic Analysis of Structures A Physical Approach With Emphasis on Earthquake Engineering Edward LWilson Professor Emeritus of Civil Engineering University of California, Berkeley Fourth Edition
More informationApplication of pseudo-symmetric technique in dynamic analysis of concrete gravity dams
Application of pseudo-symmetric technique in dynamic analysis of concrete gravity dams V. Lotfi Department of Civil and Environmental Engineering, Amirkabir University, Iran Abstract A new approach is
More informationPractical methodology for inclusion of uplift and pore pressures in analysis of concrete dams
Practical methodology for inclusion of uplift and pore pressures in analysis of concrete dams Michael McKay 1 and Francisco Lopez 2 1 Dams Engineer, GHD Pty 2 Principal Dams/Structural Engineer, GHD Pty
More informationContinuum mechanism: Stress and strain
Continuum mechanics deals with the relation between forces (stress, σ) and deformation (strain, ε), or deformation rate (strain rate, ε). Solid materials, rigid, usually deform elastically, that is the
More informationFRICTIONAL HEATING DURING AN EARTHQUAKE. Kyle Withers Qian Yao
FRICTIONAL HEATING DURING AN EARTHQUAKE Kyle Withers Qian Yao Temperature Change Along Fault Mode II (plain strain) crack rupturing bilaterally at a constant speed v r Idealize earthquake ruptures as shear
More informationMacroscopic theory Rock as 'elastic continuum'
Elasticity and Seismic Waves Macroscopic theory Rock as 'elastic continuum' Elastic body is deformed in response to stress Two types of deformation: Change in volume and shape Equations of motion Wave
More informationRheology III. Ideal materials Laboratory tests Power-law creep The strength of the lithosphere The role of micromechanical defects in power-law creep
Rheology III Ideal materials Laboratory tests Power-law creep The strength of the lithosphere The role of micromechanical defects in power-law creep Ideal materials fall into one of the following categories:
More informationAbstract. 1 Introduction
Contact analysis for the modelling of anchors in concrete structures H. Walter*, L. Baillet** & M. Brunet* *Laboratoire de Mecanique des Solides **Laboratoire de Mecanique des Contacts-CNRS UMR 5514 Institut
More informationHydraulic 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
More informationNUMERICAL ANALYSIS OF A PILE SUBJECTED TO LATERAL LOADS
IGC 009, Guntur, INDIA NUMERICAL ANALYSIS OF A PILE SUBJECTED TO LATERAL LOADS Mohammed Younus Ahmed Graduate Student, Earthquake Engineering Research Center, IIIT Hyderabad, Gachibowli, Hyderabad 3, India.
More informationSummary so far. Geological structures Earthquakes and their mechanisms Continuous versus block-like behavior Link with dynamics?
Summary so far Geodetic measurements velocities velocity gradient tensor (spatial derivatives of velocity) Velocity gradient tensor = strain rate (sym.) + rotation rate (antisym.) Strain rate tensor can
More informationI hereby declare that, except where specifically indicated, the work submitted herein is my own original work.
Finite Element Studies on the Mechanical Stability of Arbitrarily Oriented and Inclined Wellbores Using a Vertical 3D Wellbore Model by Di Zien Low (F) Fourth-Year Undergraduate Project in Group D, 2010/2011
More informationIntegrating Lab and Numerical Experiments to Investigate Fractured Rock
Integrating Lab and Numerical Experiments to Investigate Fractured Rock Bradford H. Hager Director, Earth Resources Laboratory and Cecil and Ida Green Professor Department of Earth, Atmospheric and Planetary
More informationResolving Stress Components and Earthquake Triggering
Resolving Stress Components and Earthquake Triggering Earthquake Triggering Do certain events make an earthquake more likely to occur? Earthquakes Slow Slip Wastewater Fluids Dams The focus of this presentation
More informationSpectral Elements. Introduction recalling the elastic wave equation
Spectral Elements Introduction recalling the elastic wave equation The spectral-element method: General concept domain mapping from space-continuous to space-discrete time extrapolation Gauss-Lobatto-Legendre
More informationWave Propagation in Fractured Poroelastic Media
Wave Propagation in Fractured Poroelastic Media WCCM, MS170: Advanced Computational Techniques in Geophysical Sciences, Barcelona, Spain, July 2014 Juan E. Santos Instituto del Gas y del Petróleo (IGPUBA),
More informationComparison between the visco-elastic dampers And Magnetorheological dampers and study the Effect of temperature on the damping properties
Comparison between the visco-elastic dampers And Magnetorheological dampers and study the Effect of temperature on the damping properties A.Q. Bhatti National University of Sciences and Technology (NUST),
More informationAnalytical Mechanics: Elastic Deformation
Analytical Mechanics: Elastic Deformation Shinichi Hirai Dept. Robotics, Ritsumeikan Univ. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 1 / 60 Agenda Agenda
More informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More informationSEISMOLOGY I. Laurea Magistralis in Physics of the Earth and of the Environment. Elasticity. Fabio ROMANELLI
SEISMOLOGY I Laurea Magistralis in Physics of the Earth and of the Environment Elasticity Fabio ROMANELLI Dept. Earth Sciences Università degli studi di Trieste romanel@dst.units.it 1 Elasticity and Seismic
More informationContinuum Mechanics. Continuum Mechanics and Constitutive Equations
Continuum Mechanics Continuum Mechanics and Constitutive Equations Continuum mechanics pertains to the description of mechanical behavior of materials under the assumption that the material is a uniform
More informationThe 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
More informationLawrence Berkeley National Laboratory
Lawrence Berkeley National Laboratory Peer Reviewed Title: Fracture permeability and seismic wave scattering--poroelastic linear-slip interface model for heterogeneous fractures Author: Nakagawa, S. Publication
More informationA Hybrid Method for the Wave Equation. beilina
A Hybrid Method for the Wave Equation http://www.math.unibas.ch/ beilina 1 The mathematical model The model problem is the wave equation 2 u t 2 = (a 2 u) + f, x Ω R 3, t > 0, (1) u(x, 0) = 0, x Ω, (2)
More informationSensitivity Study of Physical Limits on Ground Motion at Yucca Mountain
Bulletin of the Seismological Society of America, Vol. 100, No. 6, pp. 2996 3019, December 2010, doi: 10.1785/0120090372 Sensitivity Study of Physical Limits on Ground Motion at Yucca Mountain by Benchun
More informationA review of friction laws and their application for simulation of microseismicity prior to hydraulic fracturing
A review of friction laws and their application for simulation of microseismicity prior to hydraulic fracturing Jiyang Ye, Mirko Van Der Baan (Email: jiyang1@ualberta.ca, Mirko.VanderBaan@ualberta.ca)
More informationModel Inversion for Induced Seismicity
Model Inversion for Induced Seismicity David Castiñeira Research Associate Department of Civil and Environmental Engineering In collaboration with Ruben Juanes (MIT) and Birendra Jha (USC) May 30th, 2017
More informationComputational Seismology: An Introduction
Computational Seismology: An Introduction Aim of lecture: Understand why we need numerical methods to understand our world Learn about various numerical methods (finite differences, pseudospectal methods,
More informationINVERSE ANALYSIS METHODS OF IDENTIFYING CRUSTAL CHARACTERISTICS USING GPS ARRYA DATA
Problems in Solid Mechanics A Symposium in Honor of H.D. Bui Symi, Greece, July 3-8, 6 INVERSE ANALYSIS METHODS OF IDENTIFYING CRUSTAL CHARACTERISTICS USING GPS ARRYA DATA M. HORI (Earthquake Research
More informationYusuke Mukuhira. Integration of Induced Seismicity and Geomechanics For Better Understanding of Reservoir Physics
Integration of Induced Seismicity and Geomechanics For Better Understanding of Reservoir Physics Yusuke Mukuhira Postdoctoral Fellow (JSPS research fellow) Department of Earth, Atmospheric, and Planetary
More informationNumerical Modelling of Dynamic Earth Force Transmission to Underground Structures
Numerical Modelling of Dynamic Earth Force Transmission to Underground Structures N. Kodama Waseda Institute for Advanced Study, Waseda University, Japan K. Komiya Chiba Institute of Technology, Japan
More informationMaterials and Methods The deformation within the process zone of a propagating fault can be modeled using an elastic approximation.
Materials and Methods The deformation within the process zone of a propagating fault can be modeled using an elastic approximation. In the process zone, stress amplitudes are poorly determined and much
More informationA FINITE-VOLUME DISCRETIZATION FOR DEFORMATION OF FRACTURED MEDIA
A FINITE-VOLUME DISCRETIZATION FOR DEFORMATION OF FRACTURED MEDIA Eren Ucar a, Eirik Keilegavlen a, Inga Berre a,b, Jan Martin Nordbotten a,c a Department of Mathematics, University of Bergen, Bergen,
More informationUsing transient stresses to monitor poroelastic and stress conditions in CO 2 reservoirs
Using transient stresses to monitor poroelastic and stress conditions in CO 2 reservoirs Andrew A. Delorey and Paul A. Johnson July 6, 2016 1 Stress, Pore Pressure, and Poroelastic Behavior Induced seismicity
More informationModeling Approaches That Reproduce a Range of Fault Slip Behaviors: What We Have and What We Need Nadia Lapusta. California Institute of Technology
Modeling Approaches That Reproduce a Range of Fault Slip Behaviors: What We Have and What We Need Nadia Lapusta California Institute of Technology Modeling Approaches That Reproduce a Range of Fault Slip
More informationRUPTURE OF FRICTIONALLY HELD INCOHERENT INTERFACES UNDER DYNAMIC SHEAR LOADING
RUPTURE OF FRICTIONALLY HELD INCOHERENT INTERFACES UNDER DYNAMIC SHEAR LOADING G. Lykotrafitis and A.J. Rosakis Graduate Aeronautical Laboratories, Mail Stop 105-50, California Institute of Technology,
More informationDraft Dated: Monday, August 25, 2003 Contents of Chapter 23
1. Introduction 2. Fluid-Structure Interaction Draft Dated: Monday, August 25, 2003 Contents of Chapter 23 3. Finite Element Model Of Dam-Foundation Interface 4. Loading Due To Uplift And Pore Water Pressure
More informationNumerical Methods in Geophysics. Introduction
: Why numerical methods? simple geometries analytical solutions complex geometries numerical solutions Applications in geophysics seismology geodynamics electromagnetism... in all domains History of computers
More informationChapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature
Chapter 1 Continuum mechanics review We will assume some familiarity with continuum mechanics as discussed in the context of an introductory geodynamics course; a good reference for such problems is Turcotte
More informationfriction 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
More informationApplicability of Multi-spring Model Based on Finite Strain Theory to Seismic Behavior of Embankment on Liquefiable Sand Deposit
Applicability of Multi-spring Model Based on Finite Strain Theory to Seismic Behavior of Embankment on Liquefiable Sand Deposit Kyohei Ueda Railway Technical Research Institute, Kokubunji, Tokyo, Japan
More informationReservoir Geomechanics and Faults
Reservoir Geomechanics and Faults Dr David McNamara National University of Ireland, Galway david.d.mcnamara@nuigalway.ie @mcnamadd What is a Geological Structure? Geological structures include fractures
More informationALASKA ENERGY AUTHORITY AEA ENGINEERING FEASIBILITY REPORT. Appendix B8. Finite Element Analysis
ALASKA ENERGY AUTHORITY AEA11-022 ENGINEERING FEASIBILITY REPORT Appendix B8 Finite Element Analysis Susitna-Watana Hydroelectric Project Alaska Energy Authority FERC Project No. 14241 December 2014 Seismic
More informationIn all of the following equations, is the coefficient of permeability in the x direction, and is the hydraulic head.
Groundwater Seepage 1 Groundwater Seepage Simplified Steady State Fluid Flow The finite element method can be used to model both steady state and transient groundwater flow, and it has been used to incorporate
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
More informationFig. 1. Circular fiber and interphase between the fiber and the matrix.
Finite element unit cell model based on ABAQUS for fiber reinforced composites Tian Tang Composites Manufacturing & Simulation Center, Purdue University West Lafayette, IN 47906 1. Problem Statement In
More informationImprovement in the Fault Boundary Conditions for a Staggered Grid Finite-difference Method
Pure appl. geophys. 63 (6) 977 99 33 553/6/9977 DOI.7/s-6-8- Ó Birkhäuser Verlag, Basel, 6 Pure and Applied Geophysics Improvement in the Fault Boundary Conditions for a Staggered Grid Finite-difference
More informationA.0 IDUKKI ARCH DAM - ANALYSIS FOR VARIOUS LOAD CASES AND DISCRETISATIONS
Annexure A.0 IDUKKI ARCH DAM - ANALYSIS FOR VARIOUS LOAD CASES AND DISCRETISATIONS In this Annexure, Idukki arch dam chosen for case study is analyzed with the developed program for various load conditions
More informationNumerical Modeling of Interface Between Soil and Pile to Account for Loss of Contact during Seismic Excitation
Numerical Modeling of Interface Between Soil and Pile to Account for Loss of Contact during Seismic Excitation P. Sushma Ph D Scholar, Earthquake Engineering Research Center, IIIT Hyderabad, Gachbowli,
More informationTwo ways to think about the dynamics of earthquake ruptures
Two ways to think about the dynamics of earthquake ruptures (1) In terms of friction (2) In terms of fracture mechanics Scholz describes conditions for rupture propagation (i.e. instability) via energy
More informationLecture 8. Stress Strain in Multi-dimension
Lecture 8. Stress Strain in Multi-dimension Module. General Field Equations General Field Equations [] Equilibrium Equations in Elastic bodies xx x y z yx zx f x 0, etc [2] Kinematics xx u x x,etc. [3]
More informationMHA042 - Material mechanics: Duggafrågor
MHA042 - Material mechanics: Duggafrågor 1) For a static uniaxial bar problem at isothermal (Θ const.) conditions, state principle of energy conservation (first law of thermodynamics). On the basis of
More informationIntroduction: Advancing Simulations of Sequences of Earthquakes and Aseismic Slip (SEAS)
Introduction: Advancing Simulations of Sequences of Earthquakes and Aseismic Slip (SEAS) Brittany Erickson (Portland State University) Junle Jiang (University of California, San Diego) SCEC DR-SEAS Workshop,
More informationElastoplastic Deformation in a Wedge-Shaped Plate Caused By a Subducting Seamount
Elastoplastic Deformation in a Wedge-Shaped Plate Caused By a Subducting Seamount Min Ding* 1, and Jian Lin 2 1 MIT/WHOI Joint Program, 2 Woods Hole Oceanographic Institution *Woods Hole Oceanographic
More informationProduction-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
More informationMPM Research Community. Anura3D MPM Software. Verification Manual
MPM Research Community Anura3D MPM Software Verification Manual Version: 2017.1 12 January 2017 Anura3D MPM Software, Verification Manual Edited by: Miriam Mieremet (Deltares Delft, The Netherlands) With
More informationCritical Borehole Orientations Rock Mechanics Aspects
Critical Borehole Orientations Rock Mechanics Aspects By R. BRAUN* Abstract This article discusses rock mechanics aspects of the relationship between borehole stability and borehole orientation. Two kinds
More informationChapter 4 Analysis of a cantilever
Chapter 4 Analysis of a cantilever Before a complex structure is studied performing a seismic analysis, the behaviour of simpler ones should be fully understood. To achieve this knowledge we will start
More information3.22 Mechanical Properties of Materials Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 3.22 Mechanical Properties of Materials Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Quiz #1 Example
More informationUnderstanding hydraulic fracture variability through a penny shaped crack model for pre-rupture faults
Penny shaped crack model for pre-rupture faults Understanding hydraulic fracture variability through a penny shaped crack model for pre-rupture faults David Cho, Gary F. Margrave, Shawn Maxwell and Mark
More information