Crack-initiation and fracture propagation in reservoir rocks

Similar documents
arxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 12 Jan 2007

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

Integrating Lab and Numerical Experiments to Investigate Fractured Rock

Anisotropic permeabilities evolution of reservoir rocks under pressure:

Laboratory Shear Stimulation and Hydraulic Fracture Characterization Using Acoustic Emission

Downloaded 12/20/17 to Redistribution subject to SEG license or copyright; see Terms of Use at

Rock Mechanics Laboratory Tests for Petroleum Applications. Rob Marsden Reservoir Geomechanics Advisor Gatwick

Acquisition of Complete Acoustic Emission Amplitude Records during Rock Fracture Experiments

Laboratory Triaxial and Permeability Tests on Tournemire Shale and Cobourg Limestone (Contract : R413.8) FINAL REPORT

Temperature Dependent Mechanical Properties of Reservoir s Overburden Rocks During SAGD Process

Fracture of andesite in the brittle and brittleductile transition regimes

3D HM-DEM model for Hydro-Fracturing

Wellbore stability analysis in porous carbonate rocks using cap models

Integration of Geophysical and Geomechanical

RockLab Details. Rock Physics Testing

Experimental study of mechanical and thermal damage in crystalline hard rock

Geomechanical controls on fault and fracture distribution with application to structural permeability and hydraulic stimulation

Computers and Geotechnics

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

Simulation of the cutting action of a single PDC cutter using DEM

Activated creeping shale to remove the open annulus

J.V. Herwanger* (Ikon Science), A. Bottrill (Ikon Science) & P. Popov (Ikon Science)

Chapter 6. Conclusions. 6.1 Conclusions and perspectives

FRICTIONAL HEATING DURING AN EARTHQUAKE. Kyle Withers Qian Yao

Module 5: Failure Criteria of Rock and Rock masses. Contents Hydrostatic compression Deviatoric compression

The Frictional Regime

Module 9 : Foundation on rocks. Content

3D simulations of an injection test done into an unsaturated porous and fractured limestone

Rock and fluid thermodynamics control the dynamics of induced earthquakes

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

Rock visualization using micro-ct scanner and X-ray transparent triaxial apparatus

PETROPHYSICAL STUDIES OF SANDSTONES UNDER TRUE-TRIAXIAL STRESS CONDITIONS

Hijiori HDR Reservoir Evaluation by Micro-Earthquake Observation

20. Rheology & Linear Elasticity

Crosswell tomography imaging of the permeability structure within a sandstone oil field.

Discrete Element Modelling of a Reinforced Concrete Structure

EDEM DISCRETIZATION (Phase II) Normal Direction Structure Idealization Tangential Direction Pore spring Contact spring SPRING TYPES Inner edge Inner d

MULTISTAGE TRIAXIAL TESTING OF ACTUAL RESERVOIR CORES UNDER SIMULATED RESERVOIR CONDITIONS

Modeling pressure response into a fractured zone of Precambrian basement to understand deep induced-earthquake hypocenters from shallow injection

6298 Stress induced azimuthally anisotropic reservoir - AVO modeling

Gas Shale Hydraulic Fracturing, Enhancement. Ahmad Ghassemi

Current challenges at CO 2 Sites

Some New Developments in Rock Mechanics Research and Application

DEVELOPMENT OF AUTOMATIC CONTROL OF MULTI-STAGE TRIAXIAL TESTS AT THE UNIVERSITY OF MISKOLC

Estimating Permeability from Acoustic Velocity and Formation Resistivity Factor

Reservoir Geomechanics and Faults

Numerical modeling of standard rock mechanics laboratory tests using a finite/discrete element approach

Permeability Variations During Crack Damage Evolution in Rocks

TITLE Brazilian tensile strength test: post-failure behavior of Jurassic and Cretaceous shales from Svalbard

Measurement of acousto-elastic effect by using surface acoustic wave

Understanding hydraulic fracture variability through a penny shaped crack model for pre-rupture faults

Yifei Cui. Doctor of Philosophy. in Geotechnical Engineering. Department of Civil and Environmental Engineering University of Alberta

VMS-GeoMil. Background

Modeling of stress-dependent static and dynamic moduli of weak sandstones

Modeling of Propagation of Interacting Cracks Under Hydraulic Pressure Gradient. Hai Huang, Earl Mattson and Rob Podgorney

Gain information on the development of Engineered Geothermal Systems (EGS) and enhance our understanding of long-term reservoir behaviour.

Jihoon Kim, George J. Moridis, John Edmiston, Evan S. Um, Ernest Majer. Earth Sciences Division, Lawrence Berkeley National Laboratory 24 Mar.

DAMAGE EVALUATION OF CRACKED CONCRETE BY DECAT

MODELING GEOMATERIALS ACROSS SCALES JOSÉ E. ANDRADE DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING EPS SEMINAR SERIES MARCH 2008

Role of lithological layering on spatial variation of natural and induced fractures in hydraulic fracture stimulation

Estimation of Seismic Damage for Optimum Management of Irrigation Infrastructures in Service Conditions using Elastic Waves Information

Critical Borehole Orientations Rock Mechanics Aspects

Finite element modelling of fault stress triggering due to hydraulic fracturing

EFFECT OF BEDDING PLANES ON ROCK MECHANICAL PROPERTIES ANISOTROPY OF SANDSTONE FOR GEOMECHANICAL MODELING

D047 Change of Static and Dynamic Elastic Properties due to CO2 Injection in North Sea Chalk

Lecture 20. Measuring Pressure and Temperature (Chapter 9) Measuring Pressure Measuring Temperature MECH 373. Instrumentation and Measurements

Deliverable D4.3 - Final Report

ROTATING RING. Volume of small element = Rdθbt if weight density of ring = ρ weight of small element = ρrbtdθ. Figure 1 Rotating ring

Acoustic emission analysis for failure identification in composite materials

Long-Term Strength Degradation Testing of DGR-3 and DGR-4 Core

SEG/San Antonio 2007 Annual Meeting SUMMARY DATA

CHARACTERISTICS OF LIQUEFIED SILTY SANDS FROM MEIZOSEISMAL REGION OF SHILLONG PLATEAU, ASSAM AND BHUJ IN INDIA

John E. Gale 1 and Eunjeong Seok 2

Title. Author(s)Akm Badrul, Alam; Fujii, Yoshiaki; Aramaki, Noritaka. CitationProceedings of Japan Symposium on Rock Mechanics, 14

TRUST. Figure 1: Location of the shallow seismic wells, the injection well (H18A) and the monitoring well (H18B).

AFES Masterclass Overburden Correction. Evidence from Log Data

Main Means of Rock Stress Measurement

Geomechanics for reservoir and beyond Examples of faults impact on fluid migration. Laurent Langhi Team Leader August 2014

EXPERIMENTAL STUDY OF SHEAR AND COMPACTION BAND FORMATION IN BEREA SANDSTONE

ENERGY EMISSIONS FROM FAILURE PHENOMENA: ACOUSTIC, ELECTROMAGNETIC, PIEZONUCLEAR

Seismic Sources. Seismic sources. Requirements; Principles; Onshore, offshore. Reading: Telford et al., Section 4.5 Sheriff and Geldart, Chapter 7

An Investigation on the Effects of Different Stress Regimes on the Magnitude Distribution of Induced Seismic Events

The Mine Geostress Testing Methods and Design

Supplementary material: The relation between viscosity and acoustic emissions as a laboratory analogue for volcano seismology

CHAPTER 3 EXPERIMENTAL STUDY

Competing Effect of Pore Fluid and Texture -- Case Study

Mathematical Modelling of a Fault Slip Induced by Water Injection

DESIGN OF AN ULTRA-SPEED LAB-SCALE DRILLING RIG FOR SIMULATION OF HIGH SPEED DRILLING OPERATIONS IN HARD ROCKS. *V. Rasouli, B.

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

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

EXPERIMENTAL INVESTIGATION OF INDUCED SEISMICITY IN GRANITIC ROCK ON CENTIMETER SCALE. Linus Villiger

Research Article. Experimental Analysis of Laser Drilling Impacts on Rock Properties

5 ADVANCED FRACTURE MODELS

Insights into subdecimeter fracturing processes during the hydraulic fracture experiment in Äspö hard rock laboratory, Sweden

RESEARCH PROPOSAL. Effects of scales and extracting methods on quantifying quality factor Q. Yi Shen

scaling parameters of laboratory modelling of

Stress Damage in Borehole and Rock Cores; Developing New Tools to Update the Stress Map of Alberta

MMJ1133 FATIGUE AND FRACTURE MECHANICS A - INTRODUCTION INTRODUCTION

APPLICATION OF THE KAISER EFFECT TO THE MEASUREMENT OF IN-SITU STRESSES IN ARABIAN DEVONIAN SANDSTONE

Reservoir Geomechanics with ABAQUS

Transcription:

Crack-initiation and fracture propagation in reservoir rocks Srutarshi Pradhan Fracture & Flow meeting, January 2013, IMSc, Chennai, India 1

Dynamics of fracture Crack-growth in nano-materials onset of earthquake 2

Fracturing during fluid injection 3

Methods: 1) Fracturing test on core samples 2) PFC at pore scale 3) Beam-lattice model at mesoscopic scale 4) MDEM at macro scale (even resorvoir scale) 5) FEM using Abaqus Application: 1) Planning safe & efficient drilling (geothermal, shale-gas) 2) Reservoir characterisation (CO 2 storage) 3) Prediction of well collapse & leakage 4) Enhance production by increasing permeability. 4

SINTEF Laboratory MessTek system TerraTek system 5

Triaxial cell instrumentation Pore pressure (fluid flow) Steel piston Acoustic transducer Sintered plate Directions of radial strain measurements Sample Acoustic transducer Sleeve Axial strain LVDT Pore pressure (fluid flow) 6

AE set-up AE transducer Vallen system Max capacity => 16 channels Currently used => 9 channels Diameter = 3.5 mm Centre frequency = 1.3 MHz

Samples Sandstones Chalk Berea ρ 2150 kg/m 3 Vp ~ 2090 m/s Castlegate ρ 1920 kg/m 3 Vp ~ 1830 m/s Red Wildmoor ρ 1930 kg/m 3 Vp ~ 1590 m/s Mons chalk ρ 1520 kg/m 3 Vp ~ 2140 m/s Size of samples: Length ~ 135 mm & Diameter ~ 51 mm/10.50 mm

Stress conditions Initial load 2MPa confining pressure (Pc) & 1MPa borehole pressure (Pp) for AE calibration Main test @ const. confinement (Pc = 5 MPa) & stepwise increasing borehole pressure (Pp increase until fracture)

Berea sandstones after test Top Bottom

Berea sandstones AE results Filter 60 db Position of acoustic events (channel threshold 19.3 db)

Castlegate sandstones after test Top Bottom

Castlegate sandstones AE results Filter 60 db Position of acoustic events (channel threshold 19.3 db)

Castlegate sandstones AE results Location of acoustic events indicates two symmetric fractures localized between AE sensors 6 and 8. This fracture is visible as a clear core damage. Two additional (symmetric) fractures can be located in direction of AE sensor 5 and 7, but these are not visible with the bare eye.

Red Wildmoor sandstones after test Top Bottom

Red Wildmoor sandstones AE results Filter 60 db Position of acoustic events (channel threshold 19.3 db)

Mons chalk after test Top Bottom

Mons chalk AE results Filter 40 db Position of acoustic events (channel threshold 19.3 db)

AE events during the test sandstone limestone 19

Pore scale modeling: using PFC PFC (Particle Flow Code) is a code based on the Discrete Element Method (DEM). PFC solves the equations of motion directly. In each time step, the movements of all the particles are calculated according to the motion law, and the forces at all the contacts are calculated according to a contact law. ITASCA, USA 20

PFC: Particle-fluid interaction 21

Hydraulic fracture: PFC 2D 22

Fracture modeling at mesoscopic scale Bjørn, NTNU Elastic beam lattice model Fracture pattern and pressure distribution 23

Macro-scale modeling: MDEM Haitham, SINTEF 24

Finite element modeling (ABAQUS) Alexandre, SINTEF Fracture pattern for Gaussian dist. of tensile strength with different heterogeneity level 1 % 20 % ABAQUS model of a hollow cylinder specimen with 18 built-in cohesive interfaces 50 % 100 % 25

AE bursts Theory More fracturing activities near breakdown point Simulation Experiment 26

FBM: Avalanche statistics Burst or avalanche D( ) 5/ 2 Hemmer & Hansen, 1992 F = ( N + 1 k) k x k N = = 6 10 ; avg 20000 27

Prediction of failure point Crossover behavior D( ) (1 e 5/ 2 / c 1 c = 8( x x ) c 0 2 ) (Pradhan, Hemmer & Hansen; PRL 2005) N = = 6 10 ; avg 50000 Single sample 7 N =10 28

Crossover behavior: In other systems Seismic data prior to a mainshock (Kawamura et al, 2006) Burst distribution in Fuse model 29

Failure dynamics: Two-sided divergence Above critical stress t f (σ ) is the failure time t f ( σ ) π 1/ 2 = ( σ σ ) 2 c Below critical stress t f (σ ) is the time to reach the fixed point t f ( σ ) = ( σ σ ) ln( N ) 4 c 1/ 2 (Pradhan & Hemmer; PRE 2007) 30

The Fiber Bundle Modeling Failure in Materials A. Hansen, P. C. Hemmer and S. Pradhan In the book series Fracture, Breakdown & Earthquake Edited by B. K. Chakrabarti and P. Ray INDNOR project activity 31

Collaborators Jørn, Anna, Hans, Eyvind, Haitham, Alexandre, Erling SINTEF Alex Hansen, Per C. Hemmer, Bjørn Skjetne, NTNU Bikas K. Chakrabarti, P. Bhattacharyya, SINP FUNDING: RESEARCH COUNCIL OF NORWAY (NFR) THROUGH CLIMIT AND INDNOR (RENERGI) 32