Model Inversion for Induced Seismicity
|
|
- Hortense Watson
- 5 years ago
- Views:
Transcription
1 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
2 Outline Slide 2 1. Introduction CCS Induced Seismicity 2. Methodology Forward Modeling Inverse Modeling 3. Results for CO 2 sequestration Model Model description Model uncertainties Bayesian solution for inverse problem 4. Conclusions
3 CO 2 sequestration and Induced Seismicity Slide 3 Carbon Capture and Storage (CCS) challenges: Economic (and political!) Engineering Risks and uncertainties Risks of Induced Seismicity: Uncertainty is always present (e.g., fault stability) Evaluating impact of CO 2 injection on fault stability is more difficult. Need to quantify the risk of induced seismicity due to CO 2 injection. Source: University of Edinburgh Proposed solution to evaluate seismic risk: To develop (forward) model that simulate the fundamental physics. To implement a (inverse) problem workflow to match observations To deploy model for designing operational constraints. Source: BGR
4 Forward Modeling: Coupled flow-geo simulation Slide 4 Existing computational solutions for coupling flow and deformation. Numerical simulation framework to couple a multiphase flow simulator with a mechanics simulator. Unconditionally stable fixed-stress scheme for sequential solution of two-way coupling. Rigorous formulation of nonlinear multiphase geomechanics to handle strong capillary effects. Mechanical stability of faults modeled from fluid pressures and multiphase flow behavior. Faults represented as surfaces embedded in a 3D medium by zero-thickness interf. elements (to accurately model fault slip under dynamically evolving fluid pressure and fault strength). Reference: Jha, B., and R. Juanes (2014), Coupled multiphase flow and poromechanics: A computational model of pore pressure effects on fault slip and earthquake triggering, Water Resour. Res., 50, doi: /2013wr Applications: Algorithms can be used for prediction fault slip, fault activation, subsidence, identification of unswept zones in both real and synthetic reservoirs, etc. Research: Improving forward-model Reservoir Simulation Analytical models Data- driven models
5 Inverse Modeling: Uncertain parameters Slide 5 Large-scale parameters Tectonic boundary conditions Size of the reservoir. Aquifer description: Size of the aquifer and aquifer/reservoir permeability ratio. Fault description: strike, dip, permeability PVT parameters Rock parameters Permeability, Porosity, NTG... Geomechanical properties Poisson ratio (u ), Static coefficient of fraction, dynamic friction coefficient, compressibility intercept (a), compressible wave velocity (vp), shear wave velocity (vs), etc.
6 Inverse Modeling: Observations Slide 6 Driving index Samaria (P50) 1 Depletion Drive Index (DDI) Connate water and rock expansion Index (CRI) Water Drive Index (WDI) Gas Injection Drive Index (GIDI) Reservoir flow data Well static bottom hole pressure (WBHP data) Well flow rates (production/injection) Mechanical data Surface deformation (measured from InSAR, GPS, etc) Reservoir compaction (measured using radioactive markers). Seismic data Seismic data reflects earthquake dynamic; it should brings valuable information to the inverse problem. Problem with seismic data: it s massive, noisy, difficult to interpret without specialized knowledge. We need an effective approach to incorporate seismic data into inverse problem workflow: earthquake attributes Driving Index Date SL No Stand Pipe Parasite String TD Ecd Ql Qg Max Cuttings ECD 5.ppg % % % % % % % % ECD 4.ppg % % % % % % ECD 3.ppg % % % % % Cas Shoe ECD 2009
7 Inverse Modeling: Bayesian Solution Slide 7 Observations (D) Model parameters (q ) q D) = Posterior Distribution Bayes Rule: q D) = D q ) q ) D) = ò D q ) q ) D q ) q ) dq D q ) q ) D) = Likelihood function = Prior distribution = Evidence (or marginal distribution) Why Bayesian? Challenges: Inference from noise and limited data. Complex structure of the posterior distributions Natural mechanism for regularization Expensive forward model and large dimensionality of Natural framework for uncertainty quantification and the inversion parameters VoI analysis.
8 Inverse Modeling: Gaussian Processes Slide 8 Gaussian Processes (GP) represents a powerful regression technique for I/O data. GP finds joint distribution of infinite-dimensional multivariate normally distributed random variables. GP is actually a distribution over functions! One needs to define Kernel matrix (i.e., covariance of multivariate normal distribution, e.g., square exponential), which controls the level of smoothness in your non-linear regression model. (It) can be used for multiple input/output data Multivariate Normal Distribution: æ N ç 1 d q ik - q K( q i, q j; a ) = s exp - ç å 2 è 2 k= 1 lk 2 2 jk ö ø
9 Recap: Modeling Induced Seismicity Slide 9 Source: BGR PHYSICAL MODEL Forward Model Computational Model Predictions Data Science Observations Inverse Problem 1. Updated geological model 2. Risk Analysis and Operational Constraints
10 Case Study: 3D MODEL FOR CO 2 SEQUESTRATION 2016 Annual Founding Members Meeting 10
11 CO 2 injection (2D model) Slide 11 CO 2 injection in a deep confined aquifer for the purpose of geologic carbon sequestration. CO 2 is injected at 1500 m Aquifer is bounded on the top and bottom by a low-permeability caprock, and the fault is (in principle) impermeable to flow. Aquifer is hydraulically compartmentalized with a sealing fault that cuts across it. 2D plane-strain model with fault under normal faulting conditions (vertical principal stress due to gravity is the largest among the three principal stresses). Storage capacity of the aquifer is limited by overpressurization and possible slip on the fault. -> Modeling needed (!)
12 CO 2 injection (3D model) Slide 12 CO 2 injection in a deep confined aquifer for the purpose of geologic carbon sequestration. 3D model with (4 x 4 x 2) km simulation domain. CO 2 is injected at 1500 m in 200m anticlinal aquifer. t Aquifer is hydraulically compartmentalized with a sealing fault that cuts across it. Aquifer is bounded on the top and bottom by a low-permeability caprock. failure stable s Normal faulting conditions (vertical principal stress due to gravity is the largest among the three principal stresses). Storage capacity of the aquifer is limited by overpressurization and possible slip on the fault. -> Modeling needed (!)
13 CO 2, Uncertain Model Parameters Slide 13 Static coefficient of friction (µ f ) : It s assumed that sliding begins when the ratio of shear to normal stress on contact surfaces reaches the value of static coefficient. Dynamic coefficient of friction(µ d ): after sliding occurs the coefficient of friction decreases and reaches the dynamic friction coefficient. Slip-weakening model (d 0 ): The linear slip-weakening friction model produces shear tractions equal to the cohesive stress plus a contribution proportional to the fault normal traction that decreases from a static value to a dynamic value as slip progresses Rock permeability (k) : Permeability affects how quickly pressure affects the fault stability. 4 permeability multipliers (for 4 reservoir layers) Rock porosity (f ): Rock porosity affects the total volume of fluid in the system. Fault permeability (k fx and k fz ): To model leakage across and along fault.
14 Sensitivity analysis on 2D Model Slide 14 Parameter (Units) Earthquake Time Earthquake Location Earthquake Magnitude µ s µ d d c (meters) k (mdarcy) k ft (mdarcy) k fl (mdarcy) 0 0 ~ 0 Fault Plane
15 Prior Distributions and True (synthetic) Model Slide 15 PRIOR DISTRIBUTIONS TRUE MODEL Parameter Distribution Mean Var Min Max Parameter True Value µ s TruncNormal µ d TruncNormal d c (meters) TruncNormal k 1 (mdarcy) TruncNormal k 2 (mdarcy) TruncNormal k 3 (mdarcy) TruncNormal k 4 (mdarcy) TruncNormal f TruncNormal k f,t (mdarcy) TruncNormal 5.00e e E E- 05 k f,l (mdarcy) TruncNormal 5.00e e E E- 05 µ s 0.6 µ d 0.15 d c (meters) 0.16 k 1 (mdarcy) 200 k 2 (mdarcy) 300 k 3 (mdarcy) 250 k 4 (mdarcy) 250 f 0.09 k f,t (mdarcy) 4.33E- 06 k f,l (mdarcy) 4.87E- 06
16 Observations/Model responses Slide 16 Earthquake Time: This time is defined as the first time when rupture is observed along the fault since the start of the simulation. Earthquake Magnitude: While inertial effects are ignored in our quasi-static simulation, a finite-element integration of the total amount of slip over the fault, and the energy liberation associated to that slip, can be used as an estimate of the earthquake magnitude. ò M = G dg 0 d G f 2 M w = log10 ( M o) where Mo = seismic moment Mw = earthquake magnitude G = shear modulus d = magnitude of slip vector Earthquake Location: Hypocenter is defined as the exact location (measured as distance from origin) where faultrupture is first observed (Fixed for this study!) 16
17 Gaussian Processes (GP) Total of 1000 simulations run to build GP model (80% training, 20% validation) Slide 17 Hyperparameter selection via maximum likelihood and cross-validation
18 Solution Inverse problem using MCMC + GP Slide 18 2D Model Solution q z) = z q ) q ) = z) ò z q ) q ) z q ) q ) dq Parameter Distribution Mean Std p µs Normal µd Normal dc (meters) Normal k (mdarcy) Normal kf,t (mdarcy) Normal 7e-6 1e-6 - kf,l (mdarcy) Normal 7e-6 1e-6 - Fault Plane Bernoulli
19 Solution Inverse problem using MCMC + GP Slide 19 3D Model Solution ò = = q q q q q q q q d p z p p z p z p p z p z p ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( µ s µ d d c f k fl k ft
20 Solution Inverse problem using MCMC + GP Slide 20 k 1 k 3 Converge plots (sampling history and posterior) k 2 k 4 Geweke plots (two-sample test of means)
21 Solution Inverse problem using MCMC + GP Slide 21 Earthquake Time Earthquake Magnitude
22 Posterior Predictive Earthquake Risk Assessment Slide 22 Earthquake Magnitude CO 2 injection rate (MSCF/D)
23 Conclusions Slide Uncertainty analysis seems critical for large-scale deployment of CCS 2. Data-fit models (e.g., GP) make inverse modeling (e.g., MCMC) attainable. 3. Earthquake attributes (time, location and magnitude) can be used to infer key geological properties 4. Posterior models can be used to quantify risk envelopes for key operational parameters.
24 Questions?
Reservoir 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 informationSimplified In-Situ Stress Properties in Fractured Reservoir Models. Tim Wynn AGR-TRACS
Simplified In-Situ Stress Properties in Fractured Reservoir Models Tim Wynn AGR-TRACS Before the What and the How is Why Potential decrease in fault seal capacity Potential increase in natural fracture
More information4D stress sensitivity of dry rock frame moduli: constraints from geomechanical integration
Title 4D stress sensitivity of dry rock frame moduli: constraints from geomechanical integration Authors Bloomer, D., Ikon Science Asia Pacific Reynolds, S., Ikon Science Asia Pacific Pavlova, M., Origin
More informationAn Open Air Museum. Success breeds Success. Depth Imaging; Microseismics; Dip analysis. The King of Giant Fields WESTERN NEWFOUNDLAND:
VOL. 7, NO. 4 2010 GEOSCIENCE & TECHNOLOGY EXPLAINED GEO EXPRO VOL. 7, NO. 4 2010 Success breeds Success geoexpro.com Country Profile: Senegal Ocean Bottom Node Seismic WESTERN NEWFOUNDLAND: An Open Air
More informationJ.V. Herwanger* (Ikon Science), A. Bottrill (Ikon Science) & P. Popov (Ikon Science)
29829. One 4D geomechanical model and its many applications J.V. Herwanger* (Ikon Science), A. Bottrill (Ikon Science) & P. Popov (Ikon Science) Main objectives (i) Field case study demonstrating application
More informationMicroseismic Reservoir Monitoring
Microseismic Reservoir Monitoring Ann-Sophie Boivineau Geosciences Domain Leader, SIS Paris 1 aboivineau@slb.com Presentation Outline What are microseismic events? Applications of microseismic monitoring
More informationStatistical Rock Physics
Statistical - Introduction Book review 3.1-3.3 Min Sun March. 13, 2009 Outline. What is Statistical. Why we need Statistical. How Statistical works Statistical Rock physics Information theory Statistics
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 informationReservoir characterization in an underground gas storage field using joint inversion of flow and geodetic data
INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS Int. J. Numer. Anal. Meth. Geomech. 2015; 39:1619 1638 Published online 26 August 2015 in Wiley Online Library (wileyonlinelibrary.com)..2427
More informationInternal microseismic event pattern revealed by waveform cross-correlation analysis
Internal microseismic event pattern revealed by waveform cross-correlation analysis Bettina P. Goertz-Allmann, Volker Oye, Steven J. Gibbons, Robert Bauer, and Robert Will MGSC annual meeting September
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 informationLab 5 - Aquifer Elasticity and Specific Storage. Due October 12 14, 2010
Lab 5 - Aquifer Elasticity and Specific Storage Due October 12 14, 2010 The goal of this experiment is to measure the specific storage S s of a balloon, which simulates aquifer elasticity. The experiment
More informationModeling pressure response into a fractured zone of Precambrian basement to understand deep induced-earthquake hypocenters from shallow injection
Modeling pressure response into a fractured zone of Precambrian basement to understand deep induced-earthquake hypocenters from shallow injection S. Raziperchikolaee 1 and J. F. Miller 1 Abstract Analysis
More informationStorage: Deep Monitoring and Verification
Storage: Deep Monitoring and Verification IEA GHG Summer School, Austin, TX David White 7 th to 11 th July, 2014 www.slb.com/carbonservices CS1406-063-DW Why Monitor? Manage Risk Risk = (Impact of Undesirable
More informationTensor character of pore pressure/stress coupling in reservoir depletion and injection
Tensor character of pore pressure/stress coupling in reservoir depletion and injection Müller, B., Altmann, J.B., Müller, T.M., Weißhardt, A., Shapiro, S., Schilling, F.R., Heidbach, O. Geophysical Institute
More informationGeomechanics for reservoir and beyond Examples of faults impact on fluid migration. Laurent Langhi Team Leader August 2014
Geomechanics for reservoir and beyond Examples of faults impact on fluid migration Laurent Langhi Team Leader August 2014 Reservoir Geomechanics It is critical to understand the mechanical behaviour of
More informationGAUSSIAN PROCESS REGRESSION
GAUSSIAN PROCESS REGRESSION CSE 515T Spring 2015 1. BACKGROUND The kernel trick again... The Kernel Trick Consider again the linear regression model: y(x) = φ(x) w + ε, with prior p(w) = N (w; 0, Σ). The
More informationCurrent challenges at CO 2 Sites
Current challenges at CO 2 Sites Ola Eiken et al., Statoil R&D Force seminar on injection safety 4 th December 2013 Offshore Sleipner Onshore In Salah Sub-sea Snøhvit 1 - Classification: External 2010-09-23
More informationJihoon Kim, George J. Moridis, John Edmiston, Evan S. Um, Ernest Majer. Earth Sciences Division, Lawrence Berkeley National Laboratory 24 Mar.
TOUGH+ROCMECH for the Analysis of coupled Flow, Thermal, Geomechanical and Geophysical Processes Code Description and Applications to Tight/Shale Gas Problems Jihoon Kim, George J. Moridis, John Edmiston,
More informationCoupled geomechanics and InSAR inversion of CO2 injection parameters
11512. Scrutinizing CO2 Sequestration - A Case Study Coupling InSAR and Geomechanical Modelling to Monitor Spatial and Temporal Characteristics of CO2 Injection at In Salah, Algeria N. Gourmelen (University
More informationFully Coupled Geomechanics, Multi-Phase, Thermal and. Equation of State Compositional Simulator
Fully Coupled Geomechanics, Multi-Phase, Thermal and Equation of State Compositional Simulator Jean H. Prévost, Lee Y. Chin*, Zhihua Weng e-mail: prevost@princeton.edu URL: http://www.princeton.edu/~prevost
More informationReservoir Modeling for Wabamun Area CO2 Sequestration Project (WASP) Davood Nowroozi
Reservoir Modeling for Wabamun Area CO2 Sequestration Project (WASP) Davood Nowroozi Don Lawton 1 Effect of production/injection on Geophysical parameters Production or injection makes change in fluid
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting 5 Sep. 2017 www.geosc.psu.edu/courses/geosc508 Work of deformation, shear and volume strain Importance of volume change and diltancy rate (rate of volume strain with
More informationGaussian Processes. Le Song. Machine Learning II: Advanced Topics CSE 8803ML, Spring 2012
Gaussian Processes Le Song Machine Learning II: Advanced Topics CSE 8803ML, Spring 01 Pictorial view of embedding distribution Transform the entire distribution to expected features Feature space Feature
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 informationWP 4.1. Site selection criteria and ranking methodology. Karen Kirk
WP 4.1 Site selection criteria and ranking methodology Karen Kirk 1 Basic site selection criteria Sufficient depth and storage capacity supercritical CO 2 below 700-800 m (rule of thumb) 2 Variation of
More informationKnow Before You Go. Steve O Connor. Bryony Youngs. How GeoPrediction helps with Production Optimization and Assessing the Size of the Prize
Know Before You Go How GeoPrediction helps with Production Optimization and Assessing the Size of the Prize Steve O Connor Global Technical Lead, Geopressure Bryony Youngs Reservoir Portfolio Development
More information5 IEAGHG CCS Summer School. Geological storage of carbon dioxide (a simple solution)
Storage 1- Reservoirs, Traps, Seals and Storage Capacity for Storage Geological storage of carbon dioxide (a simple solution) Professor John Kaldi Chief Scientist, CO2CRC Australian School of Petroleum,
More informationRisk Assessment for CO 2 Sequestration- Uncertainty Quantification based on Surrogate Models of Detailed Simulations
Risk Assessment for CO 2 Sequestration- Uncertainty Quantification based on Surrogate Models of Detailed Simulations Yan Zhang Advisor: Nick Sahinidis Department of Chemical Engineering Carnegie Mellon
More informationModeling 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,
More informationTraining Venue and Dates Ref # Reservoir Geophysics October, 2019 $ 6,500 London
Training Title RESERVOIR GEOPHYSICS Training Duration 5 days Training Venue and Dates Ref # Reservoir Geophysics DE035 5 07 11 October, 2019 $ 6,500 London In any of the 5 star hotels. The exact venue
More informationDynamic System Identification using HDMR-Bayesian Technique
Dynamic System Identification using HDMR-Bayesian Technique *Shereena O A 1) and Dr. B N Rao 2) 1), 2) Department of Civil Engineering, IIT Madras, Chennai 600036, Tamil Nadu, India 1) ce14d020@smail.iitm.ac.in
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 informationThe 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
More informationFRACTURE REORIENTATION IN HORIZONTAL WELL WITH MULTISTAGE HYDRAULIC FRACTURING
SPE Workshop OILFIELD GEOMECHANICS Slide 1 FRACTURE REORIENTATION IN HORIZONTAL WELL WITH MULTISTAGE HYDRAULIC FRACTURING A. Pimenov, R. Kanevskaya Ltd. BashNIPIneft March 27-28, 2017 Moscow, Russia Slide
More informationReliability Monitoring Using Log Gaussian Process Regression
COPYRIGHT 013, M. Modarres Reliability Monitoring Using Log Gaussian Process Regression Martin Wayne Mohammad Modarres PSA 013 Center for Risk and Reliability University of Maryland Department of Mechanical
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 informationCO 2 storage capacity and injectivity analysis through the integrated reservoir modelling
CO 2 storage capacity and injectivity analysis through the integrated reservoir modelling Dr. Liuqi Wang Geoscience Australia CO 2 Geological Storage and Technology Training School of CAGS Beijing, P.
More informationAssessing the Value of Information from Inverse Modelling for Optimising Long-Term Oil Reservoir Performance
Assessing the Value of Information from Inverse Modelling for Optimising Long-Term Oil Reservoir Performance Eduardo Barros, TU Delft Paul Van den Hof, TU Eindhoven Jan Dirk Jansen, TU Delft 1 Oil & gas
More informationNumerical Study of Relationship Between Landslide Geometry and Run-out Distance of Landslide Mass
Numerical Study of Relationship Between Landslide Geometry and Run-out Distance of Landslide Mass Muneyoshi Numada Research Associate, Institute of Industrial Science, The University of Tokyo, Japan Kazuo
More informationSynthetic Seismicity Models of Multiple Interacting Faults
Synthetic Seismicity Models of Multiple Interacting Faults Russell Robinson and Rafael Benites Institute of Geological & Nuclear Sciences, Box 30368, Lower Hutt, New Zealand (email: r.robinson@gns.cri.nz).
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 informationUse of Seismic and EM Data for Exploration, Appraisal and Reservoir Characterization
Use of Seismic and EM Data for Exploration, Appraisal and Reservoir Characterization Anton Ziolkowski and Folke Engelmark Petroleum Geo-Services CSEG, Calgary, 6 May 2009 Outline Exploration, appraisal,
More informationURTeC: Abstract
URTeC: 2902950 Can Seismic Inversion Be Used for Geomechanics? A Casing Deformation Example Jeremy J. Meyer 1*, Jeremy Gallop 1, Alvin Chen 1, Scott Reynolds 1, Scott Mildren 1 ; 1. Ikon Science Copyright
More informationImplementing 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;
More informationNonparametric Bayesian Methods (Gaussian Processes)
[70240413 Statistical Machine Learning, Spring, 2015] Nonparametric Bayesian Methods (Gaussian Processes) Jun Zhu dcszj@mail.tsinghua.edu.cn http://bigml.cs.tsinghua.edu.cn/~jun State Key Lab of Intelligent
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 informationIntegration of Geophysical and Geomechanical
Integration of Geophysical and Geomechanical Modeling to Monitor Integrity of Carbon Storage Birendra Jha Assistant Professor gemlab.usc.edu Department of Chemical engineering and Materials Science University
More informationLecture: Gaussian Process Regression. STAT 6474 Instructor: Hongxiao Zhu
Lecture: Gaussian Process Regression STAT 6474 Instructor: Hongxiao Zhu Motivation Reference: Marc Deisenroth s tutorial on Robot Learning. 2 Fast Learning for Autonomous Robots with Gaussian Processes
More informationWorkflows for Sweet Spots Identification in Shale Plays Using Seismic Inversion and Well Logs
Workflows for Sweet Spots Identification in Shale Plays Using Seismic Inversion and Well Logs Yexin Liu*, SoftMirrors Ltd., Calgary, Alberta, Canada yexinliu@softmirrors.com Summary Worldwide interest
More informationDevelopments on Microseismic Monitoring and Risk Assessment of Large-scale CO 2 Storage
Developments on Microseismic Monitoring and Risk Assessment of Large-scale CO 2 Storage Bettina Goertz-Allmann NORSAR CCS Technical Workshop, RITE, Tokyo, January 16, 2019 Outline Introduction: Induced
More informationPractical Geomechanics
www.bakerhughes.com Practical Geomechanics Baker Hughes - RDS Geomechanics Services 2015 Baker Hughes Incorporated. All rights reserved Copyright By accepting these materials you agree that all materials
More informationA PRESSURE VESSEL FOR TRUE-TRIAXIAL DEFORMATION & FLUID FLOW DURING FRICTIONAL SHEAR
A PRESSURE VESSEL FOR TRUE-TRIAXIAL DEFORMATION & FLUID FLOW DURING FRICTIONAL SHEAR Chris Marone, Brett Carperter, Derek Elsworth, Igor Faoro, Matt Ikari, Matt Knuth, André Niemeijer, Demian Saffer, and
More informationGain information on the development of Engineered Geothermal Systems (EGS) and enhance our understanding of long-term reservoir behaviour.
1 GEOTHERMAL ICL was involved in the first commercial Enhanced Geothermal System (EGS) project in the EU and for over 15 years has undertaken international research programs for the monitoring of rock
More informationApply 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
More informationExploration _Advanced geophysical methods. Research Challenges. Séverine Pannetier-Lescoffit and Ute Mann. SINTEF Petroleum Research
Exploration _Advanced geophysical methods * Research Challenges Séverine Pannetier-Lescoffit and Ute Mann SINTEF Petroleum Research 1 Exploration and Reservoir Characterization * Research Challenges 29%
More informationInterpreting lab creep experiments: from single contacts... to interseismic fault models
Interpreting lab creep experiments: from single contacts...... to interseismic fault models Delphine Fitzenz, EOST Strasbourg Steve Hickman, U.S.G.S. Menlo Park Yves Bernabé, EOST, now at MIT Motivation
More informationDynamic source inversion for physical parameters controlling the 2017 Lesvos earthquake
Dynamic source inversion for physical parameters controlling the 2017 Lesvos earthquake F. Kostka 1, E. Sokos 2, J. Zahradník 1, F. Gallovič 1 1 Charles University in Prague, Faculty of Mathematics and
More informationThe Interaction of Reservoir Engineering and Geomechanics (a story)
The Interaction of Reservoir Engineering and Geomechanics (a story) Brian G D Smart FREng, FRSE, FIMMM, CEng Petromall Ltd Why is the interaction a good thing? Assertion - Reservoir Geomechanics enables
More informationUncertainties in rock pore compressibility and effects on time lapse seismic modeling An application to Norne field
Uncertainties in rock pore compressibility and effects on time lapse seismic modeling An application to Norne field Amit Suman and Tapan Mukerji Department of Energy Resources Engineering Stanford University
More informationCharacterizing Seal Bypass Systems at the Rock Springs Uplift, Southwest Wyoming, Using Seismic Attribute Analysis*
Characterizing Seal Bypass Systems at the Rock Springs Uplift, Southwest Wyoming, Using Seismic Attribute Analysis* Yuri Ganshin 1 and Fred McLaughlin 1 Search and Discovery Article #80483 (2015)** Posted
More informationRegional-Scale Salt Tectonics Modelling: Bench-Scale Validation and Extension to Field-Scale Predictions
Software Limited Bench-Scale to Field-Scale Regional-Scale Salt Tectonics Modelling: Bench-Scale Validation and Extension to Field-Scale Predictions Thornton, Dean A., Roberts, Daniel T., Crook, Anthony
More informationDefmod, an earthquake simulator that adaptively switches between quasi-static and dynamic states
Defmod, an earthquake simulator that adaptively switches between quasi-static and dynamic states C. Meng (cmeng@mit.edu) B. Hager ERL, MIT May 16, 2016 Reservoir production/injection induced seismicity
More informationOpportunities in Oil and Gas Fields Questions TABLE OF CONTENTS
TABLE OF CONTENTS A. Asset... 3 1. What is the size of the opportunity (size the prize)?... 3 2. Volumetric Evaluation... 3 3. Probabilistic Volume Estimates... 3 4. Material Balance Application... 3 5.
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 informationStochastic Spectral Approaches to Bayesian Inference
Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to
More informationON NEAR-FIELD GROUND MOTIONS OF NORMAL AND REVERSE FAULTS FROM VIEWPOINT OF DYNAMIC RUPTURE MODEL
1 Best Practices in Physics-based Fault Rupture Models for Seismic Hazard Assessment of Nuclear ON NEAR-FIELD GROUND MOTIONS OF NORMAL AND REVERSE FAULTS FROM VIEWPOINT OF DYNAMIC RUPTURE MODEL Hideo AOCHI
More informationAn Investigation on the Effects of Different Stress Regimes on the Magnitude Distribution of Induced Seismic Events
An Investigation on the Effects of Different Stress Regimes on the Magnitude Distribution of Induced Seismic Events Afshin Amini, Erik Eberhardt Geological Engineering, University of British Columbia,
More informationSTA414/2104 Statistical Methods for Machine Learning II
STA414/2104 Statistical Methods for Machine Learning II Murat A. Erdogdu & David Duvenaud Department of Computer Science Department of Statistical Sciences Lecture 3 Slide credits: Russ Salakhutdinov Announcements
More informationInduced seismicity and GeoEnergies: lessons learned from coupled hydro-mechanical modeling
Induced seismicity and GeoEnergies: lessons learned from coupled hydro-mechanical modeling Antonio P. Rinaldi and many others. 4 th Cargèse Summer School: Flow and Transport in Porous and Fractured Media
More informationNeutron inverse kinetics via Gaussian Processes
Neutron inverse kinetics via Gaussian Processes P. Picca Politecnico di Torino, Torino, Italy R. Furfaro University of Arizona, Tucson, Arizona Outline Introduction Review of inverse kinetics techniques
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 informationPetrophysical Data Acquisition Basics. Coring Operations Basics
Petrophysical Data Acquisition Basics Coring Operations Basics Learning Objectives By the end of this lesson, you will be able to: Understand why cores are justified and who needs core data Identify which
More informationMaterial is perfectly elastic until it undergoes brittle fracture when applied stress reaches σ f
Material is perfectly elastic until it undergoes brittle fracture when applied stress reaches σ f Material undergoes plastic deformation when stress exceeds yield stress σ 0 Permanent strain results from
More informationSENSITIVITY ANALYSIS IN NUMERICAL SIMULATION OF MULTIPHASE FLOW FOR CO 2 STORAGE IN SALINE AQUIFERS USING THE PROBABILISTIC COLLOCATION APPROACH
XIX International Conference on Water Resources CMWR 2012 University of Illinois at Urbana-Champaign June 17-22,2012 SENSITIVITY ANALYSIS IN NUMERICAL SIMULATION OF MULTIPHASE FLOW FOR CO 2 STORAGE IN
More informationStorage 4 - Modeling for CO 2 Storage. Professor John Kaldi Chief Scientist, CO2CRC Australian School of Petroleum, University of Adelaide, Australia
Storage 4 - Modeling for CO 2 Storage Professor John Kaldi Chief Scientist, CO2CRC Australian School of Petroleum, University of Adelaide, Australia 1 Modelling 2 On Models. All models are wrong. some
More informationFully Bayesian Deep Gaussian Processes for Uncertainty Quantification
Fully Bayesian Deep Gaussian Processes for Uncertainty Quantification N. Zabaras 1 S. Atkinson 1 Center for Informatics and Computational Science Department of Aerospace and Mechanical Engineering University
More informationMUDLOGGING, CORING, AND CASED HOLE LOGGING BASICS COPYRIGHT. Coring Operations Basics. By the end of this lesson, you will be able to:
LEARNING OBJECTIVES MUDLOGGING, CORING, AND CASED HOLE LOGGING BASICS Coring Operations Basics By the end of this lesson, you will be able to: Understand why cores are justified and who needs core data
More informationJournal of Geophysical Research Letters Supporting Information for
Journal of Geophysical Research Letters Supporting Information for InSAR observations of strain accumulation and fault creep along the Chaman Fault system, Pakistan and Afghanistan H. Fattahi 1, F. Amelung
More informationNonparameteric Regression:
Nonparameteric Regression: Nadaraya-Watson Kernel Regression & Gaussian Process Regression Seungjin Choi Department of Computer Science and Engineering Pohang University of Science and Technology 77 Cheongam-ro,
More informationNet-to-gross from Seismic P and S Impedances: Estimation and Uncertainty Analysis using Bayesian Statistics
Net-to-gross from Seismic P and S Impedances: Estimation and Uncertainty Analysis using Bayesian Statistics Summary Madhumita Sengupta*, Ran Bachrach, Niranjan Banik, esterngeco. Net-to-gross (N/G ) is
More informationChecking up on the neighbors: Quantifying uncertainty in relative event location
Checking up on the neighbors: Quantifying uncertainty in relative event location The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation
More informationRecommendations for Injection and Storage Monitoring
Energy and Environmental Systems Group Institute for Sustainable Energy, Environment and Economy (ISEEE) Recommendations for Injection and Storage Monitoring WABAMUN AREA CO 2 SEQUESTRATION PROJECT (WASP)
More informationResearch in Induced Seismicity
Research in Induced Seismicity Grant Bromhal, NETL GWPC Annual UIC Meeting January 22-24. 2013 1 NAS Study on Induced Seismicity Three major findings emerged from the study: hydraulic fracturing does not
More informationSource parameters II. Stress drop determination Energy balance Seismic energy and seismic efficiency The heat flow paradox Apparent stress drop
Source parameters II Stress drop determination Energy balance Seismic energy and seismic efficiency The heat flow paradox Apparent stress drop Source parameters II: use of empirical Green function for
More informationMultivariate Bayesian Linear Regression MLAI Lecture 11
Multivariate Bayesian Linear Regression MLAI Lecture 11 Neil D. Lawrence Department of Computer Science Sheffield University 21st October 2012 Outline Univariate Bayesian Linear Regression Multivariate
More informationSTA414/2104. Lecture 11: Gaussian Processes. Department of Statistics
STA414/2104 Lecture 11: Gaussian Processes Department of Statistics www.utstat.utoronto.ca Delivered by Mark Ebden with thanks to Russ Salakhutdinov Outline Gaussian Processes Exam review Course evaluations
More informationADVANCED MACHINE LEARNING ADVANCED MACHINE LEARNING. Non-linear regression techniques Part - II
1 Non-linear regression techniques Part - II Regression Algorithms in this Course Support Vector Machine Relevance Vector Machine Support vector regression Boosting random projections Relevance vector
More informationMicroseismic Geomechanical Modelling of Asymmetric Upper Montney Hydraulic Fractures
Microseismic Geomechanical Modelling of Asymmetric Upper Montney Hydraulic Fractures Drew Chorney, Byungtark Lee, Shawn Maxwell (IMaGE) Summary Geomechanical modelling is a powerful tool to quantitatively
More informationModeling of Fault Reactivation and Induced Seismicity During Hydraulic Fracturing of Shale-Gas Reservoirs
Modeling of Fault Reactivation and Induced Seismicity During Hydraulic Fracturing of Shale-Gas Reservoirs Jonny Rutqvist *, Antonio P. Rinaldi, Frédéric Cappa, George J. Moridis Lawrence Berkeley National
More informationBRIEFING MEMO ON RESERVOIR TRIGGERED SEISMICITY (RTS)
BRIEFING MEMO ON RESERVOIR TRIGGERED SEISMICITY (RTS) 1. General. The issue of reservoir-triggered seismicity (RTS) has been controversial, and hotly debated, for many decades. There has been general recognition
More informationAbstracts ESG Solutions
Abstracts ESG Solutions 2015-2016 For more information, please contact Melissa Hoy, Technical Marketing Coordinator at melissa.hoy@esgsolutions.com Visit us online at www.esgsolutions.com Abstract #1 Fracture
More informationProbabilistic numerics for deep learning
Presenter: Shijia Wang Department of Engineering Science, University of Oxford rning (RLSS) Summer School, Montreal 2017 Outline 1 Introduction Probabilistic Numerics 2 Components Probabilistic modeling
More information(Refer Slide Time: 01:15)
Soil Mechanics Prof. B.V.S. Viswanathan Department of Civil Engineering Indian Institute of Technology, Bombay Lecture 56 Stability analysis of slopes II Welcome to lecture two on stability analysis of
More informationEstimating fault slip rates, locking distribution, elastic/viscous properites of lithosphere/asthenosphere. Kaj M. Johnson Indiana University
3D Viscoelastic Earthquake Cycle Models Estimating fault slip rates, locking distribution, elastic/viscous properites of lithosphere/asthenosphere Kaj M. Johnson Indiana University In collaboration with:
More informationRock and fluid thermodynamics control the dynamics of induced earthquakes
Rock and fluid thermodynamics control the dynamics of induced earthquakes M. Acosta, F. Passelègue, A. Schubnel et M. Violay (mateo.acosta@epfl.ch) Context Induced Seismicity in Enhanced Geothermal Systems
More informationTexNet and CISR: An Update on Monitoring and Understanding Seismicity in Texas
TexNet and CISR: An Update on Monitoring and Understanding Seismicity in Texas Michael H. Young a, Alexandros Savvaidis a, Peter Hennings a, Ellen Rathje b a Bureau of Economic Geology, Jackson School
More informationNew Insights into History Matching via Sequential Monte Carlo
New Insights into History Matching via Sequential Monte Carlo Associate Professor Chris Drovandi School of Mathematical Sciences ARC Centre of Excellence for Mathematical and Statistical Frontiers (ACEMS)
More informationResults of the Lacq pilot s monitoring Focus on microseismicity. Joëlle HY-BILLIOT 10 may - Venice
DGEP/SCR/RD Results of the Lacq pilot s monitoring Focus on microseismicity Joëlle HY-BILLIOT 10 may - Venice CCS pilot, Lacq, France Introduction Key areas Presentation of the Lacq pilot Monitoring program
More informationOne or two deficiencies of current fault seal analysis methods
Shell Exploration & Production One or two deficiencies of current fault seal analysis methods Scott J. Wilkins and Stephen J. Naruk Copyright 2003 SIEP B.V. Structure, Traps and Seals Team Bellaire Technology
More information