Microseismic Event Estimation Via Full Waveform Inversion
|
|
- Justina Williamson
- 5 years ago
- Views:
Transcription
1 Microseismic Event Estimation Via Full Waveform Inversion Susan E. Minkoff 1, Jordan Kaderli 1, Matt McChesney 2, and George McMechan 2 1 Department of Mathematical Sciences, University of Texas at Dallas 2 Department of Geosciences, University of Texas at Dallas Richardson, TX 758, USA IMA Hot Topics Workshop Hydraulic Fracturing: From Modeling and Simulation to Reconstruction and Characterization May 13, 215
2 Exploration Seismology Exploration is the search for commercial deposits of useful minerals, including hydrocarbons, geothermal resources, etc. Source generates a disturbance (wave). Wave reflects when it encounters a change in material properties. Returning wave is recorded using seismometers at the earth s surface. We can analyze the recorded data to image the earth s interior. Source can be intentional or active (e.g. explosion) or unintentional or passive (e.g. earthquake). Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 2 / 31
3 Hydraulic Fracturing (a) Hydraulic Fracturing Fracking is used to extract oil and gas from materials with low permeability such as shale. High pressure liquid is injected into the well to create fracture openings that allow oil and gas to flow more freely. Buildup of pressure and stress may result in a microseismic event (small earthquake), thus generating a passive source for imaging. Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 3 / 31
4 Hydraulic Fracturing Simulation and Elastic Response (the Big Picture) Processing Workflow: Using the Complex Fracture Research Code (CFRAC) (courtesy of Marc McClure, UT Austin, CPGE), model coupled flow and deformation to synthesize microseismic events produced by hydraulic fracturing. The code models hydraulic injection of a single phase fluid into an impermeable isotropic medium with discrete fractures. Deformation is modeled in a 2D plane assuming quasistatic equilibrium. Extract synthesized microseismic events from the flow/deformation code and inject the microseisms into an 3D (visco)elastic wave modeling code (courtesy of George McMechan, UTD, Geosciences). Model the elastic wave response resulting from these microseismic sources. Result could be used as forward model for inversion. Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 4 / 31
5 Numerical Example: Problem Setup 4 Fracture Geometry 8 Well Stimulation Treatment Curves Y (m) Injection Rate (kg/s) Injection Pressure (MPa) X (m) (b) Fracture Geometry: black line is well; blue line is natural fracture; red lines are (potential) hydraulic fractures Time (sec) (c) Well Treatment Curve Student Version of MATLAB Student Version of MATLAB Model contains a single natural fracture transecting an open wellbore. Strike of natural fracture is 6. Two potentially forming hydraulic fractures at ends of natural fracture (strike is 9 ). Injection pressure held constant at 6 MPa. Injection rate varies with time as failures occur along the fracture(s). Note: additional injection rate variation after about 5 s due to failure on hydraulic fracture. Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 5 / 31
6 Numerical Example: Microseismic Event Generation 4 Fracture Geometry with Microseismic Hypocenters 8 Well Stimulation Treatment Curves with Microseismic Events Y (m) Injection Rate (kg/s) Microseismic Moment Magnitude X (m) Time (sec) (d) Fracture Geometry with microseisms (e) Well Treatment Curve with microseisms Microseismic events defined to occur when sliding velocity reaches.5 m/s. Student Version of MATLAB One event has occurred Student along Version of MATLAB the potentially forming hydraulic fracture. Microseismic events correlate with changes in flow rate. Microseismic event at 525 s corresponds to failure on hydraulic fracture. Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 6 / 31
7 Microseism 13 Well Natural Fractures Hydraulic Fractures Microseismic Hypocenters Geophone Locations Numerical Example: Elastic Wave Modeling from Microseismic Event 13 4 v x snapshot 8 ms x y cordianate [m] Seismic Amplitude (m/s) x coordinate [m] (f) Fracture Geometry with microseismic hypocenters (g) X-component of particle velocity modeled using microseismic event 13 Student Version of MATLAB 4 v y snapshot 8 ms x v z snapshot 8 ms x y cordianate [m] Seismic Amplitude (m/s) y cordianate [m] Seismic Amplitude (m/s) x coordinate [m] x coordinate [m] (h) Particle velocity (y) (i) Particle velocity (z) Student Version of MATLAB Student Version of MATLAB Microseismic source has dip of 9, strike of 9, and rake of. Small microseismic event occurred along potentially forming hydraulic fracture. Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 7 / 31
8 Microseism 14 Well Natural Fractures Hydraulic Fractures Microseismic Hypocenters Geophone Locations Numerical Example: Elastic Wave Modeling from Microseismic Event 14 4 v x snapshot 8 ms y cordianate [m] Seismic Amplitude (m/s) x coordinate [m] (j) Fracture Geometry with microseismic hypocenters 4 35 v y snapshot 8 ms.1.8 (k) X-component of particle velocity modeled using microseismic event v z snapshot 8 ms Student Version of MATLAB y cordianate [m] Seismic Amplitude (m/s) y cordianate [m] Seismic Amplitude (m/s) x coordinate [m] x coordinate [m] (l) Particle velocity (y) (m) Particle velocity (z) Microseismic source has dip of 9, strike of 6, and rake of. Student Version of MATLAB Student Version of MATLAB Microseismic event occurred along natural fracture. Note the directional polarity of the emitted wavefield. Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 8 / 31
9 Numerical Example: Receiver Data for Microseismic Event 13 X-component of particle velocity as recorded at the receiver locations. Note the existence of both P and S-wave arrivals in the recorded seismograms. Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 9 / 31
10 Outline Goal is to estimate the microseismic source using full waveform inversion. I. Description of Problem and Approach II. Overview of Full Waveform Inversion III. Numerical Experiments for Source Location and Onset Time IV. Discussion of Future Work Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 1 / 31
11 Traditional Methods for Microseismic Event Detection Figure: Traditional Methods rely on traveltime picking of first arrivals. (a) Event estimation using clean data (b) Event estimation using noisy data Reference (courtesy of Pioneer Natural Resources): 1. Han L., J. Wong, and J. C. Bancroft, Hypocenter location using hodogram analysis of noisy 3C microseismograms, CREWES Research Report - Vol. 21, University of Calgary, (29), pp Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 11 / 31
12 Seismic Full Waveform Inversion 1 2 u c 2 t 2 2 u = s (c) Forward Problem (d) Velocity Inversion (e) Source Inversion Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 12 / 31
13 Assumptions We model wave propagation through a constant density medium using the acoustic wave equation: m(x, z) 2 u(x, z, t) 2 u(x, z, t) = f (x, z)w(t) t 2 where m(x, z) is the known squared slowness, i.e. the reciprocal of the sound velocity squared, f (x, z) is the spatial component of the source of acoustic energy, w(t) is a Ricker wavelet, and u(x, z, t) is the acoustic pressure. In the typical inverse problem, we know f (x, z)w(t), and we want to determine the sound velocity, i.e. we want to find m(x, z). In our source inversion problem, we know m(x, z), and we want to find f (x, z)w(t). Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 13 / 31
14 Overview of Full Waveform Inversion 1 1 Choose initial guess f. 2 Given the correct sound velocity and the current source estimate f n, solve the wave equation u n = F [f n]. 3 Compute the residual u n d obs, where d obs is the recorded data. 4 Solve the adjoint equation backwards in time using the residual as the source F (u n d obs ). 5 Compute the gradient to find a descent direction for the least squares objective function J(f ) = 1 2 F (f ) d obs 2. 6 Update the model f n+1 = f n α gradient. 7 Increment n and loop back to Step 2. 1 Note: all experiments shown were implemented using a modified version of PySIT, a seismic inversion toolbox developed originally by the Imaging and Computing Group at MIT. Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 14 / 31
15 Comparison of gradient for velocity inversion vs. source inversion In the typical velocity inversion problem, J(m) = 1 2 In our source inversion problem, J(f ) = 1 2 P In the typical velocity inversion problem, δ mj = P s s,r P s,r h R i T (Ss,r u ds,r )2 dt. h R i T (Ss,r u ds,r )2 dt. h R i T ( m 2 u, λ)dt. t 2 In our source inversion problem, δ f J = P s h R T (w, λ)dt i. We do not need to store the forward model data for all time steps. Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 15 / 31
16 Overview of Full Waveform Inversion Step 1: Given the correct sound velocity and an initial guess for the source f, solve the wave equation u = F [f ]. (f) Experimental Setup (g) Receiver Data for Initial Source Note that F is the forward operator, i.e. the wave equation operator, and maps from the x-z domain to the t-x domain. Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 16 / 31
17 Overview of Full Waveform Inversion Step 2: Compute the residual u n d obs. (h) Receiver Data for Initial Source (i) Recorded Data (j) Residual Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 17 / 31
18 Overview of Full Waveform Inversion Step 3: Solve the adjoint equation backwards in time using the residual as the source F (u n d obs ). (k) Residual (l) Adjoint Field at time t = 1.5 (m) Adjoint Field at time t = Note that F (the adjoint operator) maps from the t-x domain to the x-z domain. Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 18 / 31
19 Overview of Full Waveform Inversion Step 4: Compute the gradient to find a descent direction for the least squares objective function J(f ) = 1 2 F (f ) d obs 2. (n) Gradient Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 19 / 31
20 Overview of Full Waveform Inversion Step 5: Update the model f n+1 = f n α gradient. (o) f n (p) α gradient (q) f n+1 Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 2 / 31
21 Overview of Full Waveform Inversion Step 6: Increment n and loop back to Step 1: Given the correct sound velocity and the current source estimate f n, solve the wave equation u n = F [f n]. (r) Receiver Data after 1 Iteration (s) Recorded Data Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 21 / 31
22 Numerical Experiments: inverting for spatial component of the source (f ) Note: for all the experiments shown, Domain is size 181 x 141. Velocity is assumed constant and known. For this experiment Time-dependent component of source is 2 Hz Ricker wavelet (known). Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 22 / 31
23 Inverting for spatial component of the source using noisy data Initial source estimate is zero everywhere. Data has a.8 signal-to-noise ratio!! Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 23 / 31
24 Inverting for time-dependent part of the source (w) Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 24 / 31
25 Inverting for time-dependent part of the source using noisy data Data has.83 signal-to-noise ratio!! Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 25 / 31
26 Joint inversion for both the time and space components of the source Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 26 / 31
27 Joint inversion for both the time and space components of the source Problem: source is too spread out in space which has a corresponding impact on the time dependence. Need to focus the source. Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 27 / 31
28 Adding L 1 Regularization With regularization our new objective function for the spatial component of the source is J(f ) = 1 2 F (f ) d obs λ f 1 where λ is a penalty parameter. Anywhere in the domain where the source is zero (in our case, almost everywhere), the second term above is not differentiable. Several L 1 minimization algorithms attempt to handle this problem. We focus on Iteratively Reweighted Least Squares (IRLS 12 ): Away from f i =, gradient of the regularization term f 1 contains components f i / f i. Defining a diagonal matrix W with elements W i,i = 1/ f i we can handle the non-differentiability at f i = by j 1/ fi f i ɛ W i,i = 1/ɛ f i < ɛ or by where ɛ is a tolerance, and f 1 = Wf. q W i,i = 1/ (fi 2 + ɛ 2 ) 1 Daubechies, Devore, Fornasier, Gunturk, Iteratively reweighted least squares minimization for sparse recovery. Communications on Pure and Applied Mathematics, v.63, no.1., Aster, Borchers, and Thurber. Parameter Estimation and Inverse Problems. Elsevier, 25. Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 28 / 31
29 Numerical Experiment: inverting for spatial component of the source (f ) with L 1 Regularization. Domain is size 181 x 141. Velocity is assumed constant and known. Time-dependent component of source is 2 Hz Ricker wavelet (known). ɛ = Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 29 / 31
30 Numerical Experiment: inverting for spatial component of the source (f ) with L 1 Regularization. Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 3 / 31
31 Conclusions and Future Work Hydraulic fracturing of low permeability shale rock creates fractures. Buildup of pressure and mechanical stress may lead to small releases of elastic energy (earthquakes) which are called microseismic events. These passive sources can be used to further refine our images of the earth s subsurface (e.g., estimates of the velocity). To estimate these sources traditional methods have relied on picking first arrivals of seismic energy in the data. Unfortunately, first arrival picking is very difficult in the presence of noise. By using full waveform inversion, we have successfully found the spatial location of the microseismic event as well as the onset time even in the presence of substantial noise and with poor initial guesses. Joint inversion for both the space and time component of the source may require additional information from well logs, etc. Our future work will include Extending the work to more realistic problems, including more realistic representations of the source. Coupling the big picture with the source inversion via FWI. Adding uncertainty quantification in the context of joint inversion for velocity and the source. Susan E. Minkoff (UTD) Microseismic Source Estimation using FWI 5/13/15 31 / 31
Checking 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 informationLocation uncertainty for a microearhquake cluster
Analysis of location uncertainty for a microearhquake cluster: A case study Gabriela Melo, Alison Malcolm, Oleg Poliannikov, and Michael Fehler Earth Resources Laboratory - Earth, Atmospheric, and Planetary
More informationCharacterization of Induced Seismicity in a Petroleum Reservoir: A Case Study
Characterization of Induced Seismicity in a Petroleum Reservoir: A Case Study Edmond Sze, M. Nafi Toksöz, and Daniel R. Burns Earth Resources Laboratory Dept. of Earth, Atmospheric and Planetary Sciences
More informationRegistration-guided least-squares waveform inversion
Registration-guided least-squares waveform inversion Hyoungsu Baek 1, Henri Calandra, Laurent Demanet 1 1 MIT Mathematics department, TOTAL S.A. January 15 013 Abstract Full waveform inversion with frequency
More informationMicroseismicity applications in hydraulic fracturing monitoring
Available online atwww.scholarsresearchlibrary.com Archives of Applied Science Research, 2016, 8 (4):13-19 (http://scholarsresearchlibrary.com/archive.html) ISSN 0975-508X CODEN (USA) AASRC9 Microseismicity
More informationComparison between least-squares reverse time migration and full-waveform inversion
Comparison between least-squares reverse time migration and full-waveform inversion Lei Yang, Daniel O. Trad and Wenyong Pan Summary The inverse problem in exploration geophysics usually consists of two
More informationMonitoring induced microseismic events usually
SPECIAL M i c r SECTION: o s e i s m M i ci c r o s e i s m i c Beyond the dots in the box: microseismicity-constrained fracture models for reservoir simulation Leo Eisner, Sherilyn Williams-Stroud, Andrew
More informationTwo-Scale Wave Equation Modeling for Seismic Inversion
Two-Scale Wave Equation Modeling for Seismic Inversion Susan E. Minkoff Department of Mathematics and Statistics University of Maryland Baltimore County Baltimore, MD 21250, USA RICAM Workshop 3: Wave
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 informationPREDICTIVE MODELING OF INDUCED SEISMICITY: NUMERICAL APPROACHES, APPLICATIONS, AND CHALLENGES
PREDICTIVE MODELING OF INDUCED SEISMICITY: NUMERICAL APPROACHES, APPLICATIONS, AND CHALLENGES Mark McClure Assistant Professor Petroleum and Geosystems Engineering The University of Texas at Austin Overview
More informationVollständige Inversion seismischer Wellenfelder - Erderkundung im oberflächennahen Bereich
Seminar des FA Ultraschallprüfung Vortrag 1 More info about this article: http://www.ndt.net/?id=20944 Vollständige Inversion seismischer Wellenfelder - Erderkundung im oberflächennahen Bereich Thomas
More informationHijiori HDR Reservoir Evaluation by Micro-Earthquake Observation
GRC Transactions, Vol. 38, 2014 Hijiori HDR Reservoir Evaluation by Micro-Earthquake Observation Hideshi Kaieda Central Research Institute of Electric Power Industry, Abiko, Chiba, Japan Keywords HDR,
More informationObservation of shear-wave splitting from microseismicity induced by hydraulic fracturing: A non-vti story
Observation of shear-wave splitting from microseismicity induced by hydraulic fracturing: A non-vti story Petr Kolinsky 1, Leo Eisner 1, Vladimir Grechka 2, Dana Jurick 3, Peter Duncan 1 Summary Shear
More informationAnisotropic Seismic Imaging and Inversion for Subsurface Characterization at the Blue Mountain Geothermal Field in Nevada
PROCEEDINGS, 43rd Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 12-14, 2018 SGP-TR-213 Anisotropic Seismic Imaging and Inversion for Subsurface Characterization
More informationCalculation of Focal mechanism for Composite Microseismic Events
Calculation of Focal mechanism for Composite Microseismic Events Hongliang Zhang, David W. Eaton Department of Geoscience, University of Calgary Summary It is often difficult to obtain a reliable single-event
More informationAttenuation compensation in least-squares reverse time migration using the visco-acoustic wave equation
Attenuation compensation in least-squares reverse time migration using the visco-acoustic wave equation Gaurav Dutta, Kai Lu, Xin Wang and Gerard T. Schuster, King Abdullah University of Science and Technology
More informationWenyong Pan and Lianjie Huang. Los Alamos National Laboratory, Geophysics Group, MS D452, Los Alamos, NM 87545, USA
PROCEEDINGS, 44th Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 11-13, 019 SGP-TR-14 Adaptive Viscoelastic-Waveform Inversion Using the Local Wavelet
More informationFull Waveform Inversion via Matched Source Extension
Full Waveform Inversion via Matched Source Extension Guanghui Huang and William W. Symes TRIP, Department of Computational and Applied Mathematics May 1, 215, TRIP Annual Meeting G. Huang and W. W. Symes
More informationComparison of Microseismic Results in Complex Geologies Reveals the Effect of Local Stresses on Fracture Propagation
Comparison of Microseismic Results in Complex Geologies Reveals the Effect of Local Stresses on Fracture Propagation Mike Preiksaitis, Sheri Bowman, Ted Urbancic ESG Canada Inc. Summary Understanding the
More informationWaveform inversion and time-reversal imaging in attenuative TI media
Waveform inversion and time-reversal imaging in attenuative TI media Tong Bai 1, Tieyuan Zhu 2, Ilya Tsvankin 1, Xinming Wu 3 1. Colorado School of Mines 2. Penn State University 3. University of Texas
More informationLaboratory Shear Stimulation and Hydraulic Fracture Characterization Using Acoustic Emission
31 st Conference of the European Working Group on Acoustic Emission (EWGAE) Th.2.A.1 More Info at Open Access Database www.ndt.net/?id=17553 Laboratory Shear Stimulation and Hydraulic Fracture Characterization
More informationTowards Modelling Elastic and Viscoelastic Seismic Wave Propagation in Boreholes
Towards Modelling Elastic and Viscoelastic Seismic Wave Propagation in Boreholes NA WANG, DONG SHI, BERND MILKEREIT Department of Physics, University of Toronto, Toronto, Canada M5S 1A7 Summary We are
More informationSUMMARY INTRODUCTION. f ad j (t) = 2 Es,r. The kernel
The failure mode of correlation focusing for model velocity estimation Hyoungsu Baek 1(*), Henri Calandra 2, and Laurent Demanet 1 1 Dept. of Mathematics and Earth Resources Lab, Massachusetts Institute
More informationChapter 7. Seismic imaging. 7.1 Assumptions and vocabulary
Chapter 7 Seismic imaging Much of the imaging procedure was already described in the previous chapters. An image, or a gradient update, is formed from the imaging condition by means of the incident and
More informationMicroseismic Aids In Fracturing Shale By Adam Baig, Sheri Bowman and Katie Jeziorski
AUGUST 2014 The Better Business Publication Serving the Exploration / Drilling / Production Industry Microseismic Aids In Fracturing Shale By Adam Baig, Sheri Bowman and Katie Jeziorski KINGSTON, ONTARIO
More informationOil and natural gas production from shale formations
SPECIAL Passive SECTION: seismic Pand a s s ive microseismic Part and microseismic Part 2 2 Source characteristics of seismicity associated with underground wastewater disposal: A case study from the 2008
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 informationThe effect of location error on microseismic mechanism estimation: synthetic and real field data examples
The effect of location error on microseismic mechanism estimation: synthetic and real field data examples Michael Kratz 1 and Michael Thornton 1 present an issue that is of primary concern for all basins
More informationEOS 350 MIDTERM OCT 4, 2013 STUDENT NAME: TEAM #:
EOS 350 MIDTERM OCT 4, 2013 STUDENT NAME: TEAM #: Some equations which may, or may not, be useful: Distance from sensor to a dipole z ~ x ½, Distance to line of dipoles z ~ 0.75x ½ B = μh, M = κh Seismic
More information23855 Rock Physics Constraints on Seismic Inversion
23855 Rock Physics Constraints on Seismic Inversion M. Sams* (Ikon Science Ltd) & D. Saussus (Ikon Science) SUMMARY Seismic data are bandlimited, offset limited and noisy. Consequently interpretation of
More informationScienceDirect. Model-based assessment of seismic monitoring of CO 2 in a CCS project in Alberta, Canada, including a poroelastic approach
Available online at www.sciencedirect.com ScienceDirect Energy Procedia 63 (2014 ) 4305 4312 GHGT-12 Model-based assessment of seismic monitoring of CO 2 in a CCS project in Alberta, Canada, including
More informationChapter 6. Conclusions. 6.1 Conclusions and perspectives
Chapter 6 Conclusions 6.1 Conclusions and perspectives In this thesis an approach is presented for the in-situ characterization of rocks in terms of the distribution of hydraulic parameters (called SBRC
More informationSHALE GAS AND HYDRAULIC FRACTURING
SHALE GAS AND HYDRAULIC FRACTURING JAMES VERDON FRONTIERS IN EARTH SCIENCES (EASC M0016) WEDNESDAY 29.11.2011 Tight gas reservoirs have permeability of 0.1mD or less difficult to produce at economic rates.
More informationSeismic tomography with co-located soft data
Seismic tomography with co-located soft data Mohammad Maysami and Robert G. Clapp ABSTRACT There is a wide range of uncertainties present in seismic data. Limited subsurface illumination is also common,
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 informationUncertainty quantification for Wavefield Reconstruction Inversion
Uncertainty quantification for Wavefield Reconstruction Inversion Zhilong Fang *, Chia Ying Lee, Curt Da Silva *, Felix J. Herrmann *, and Rachel Kuske * Seismic Laboratory for Imaging and Modeling (SLIM),
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 informationPassive seismic monitoring in unconventional oil and gas
Passive seismic monitoring in unconventional oil and gas Michael Kendall, James Verdon, Alan Baird, Anna Stork and Philip Usher Bristol University Microseismicity Projects (BUMPS) Microseismicity and
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 informationFocal Mechanism Analysis of a Multi-lateral Completion in the Horn River Basin
Focal Mechanism Analysis of a Multi-lateral Completion in the Horn River Basin Paige Snelling*, Cameron Wilson, MicroSeismic Inc., Calgary, AB, Canada psnelling@microseismic.com Neil Taylor, Michael de
More informationInterpretation of baseline surface seismic data at the Violet Grove CO 2 injection site, Alberta
Violet Grove seismic interpretation Interpretation of baseline surface seismic data at the Violet Grove CO 2 injection site, Alberta Fuju Chen and Don Lawton ABSTRACT Time-lapse seismic technology has
More informationImaging complex structure with crosswell seismic in Jianghan oil field
INTERPRETER S CORNER Coordinated by Rebecca B. Latimer Imaging complex structure with crosswell seismic in Jianghan oil field QICHENG DONG and BRUCE MARION, Z-Seis, Houston, Texas, U.S. JEFF MEYER, Fusion
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 informationTheory. Summary. Introduction
Waveform similarity for quality control of event locations, time picking and moment tensor solutions Fernando Castellanos, University of Alberta. Edmonton, AB. Mirko van der Baan, University of Alberta.
More informationThis paper was prepared for presentation at the Unconventional Resources Technology Conference held in Denver, Colorado, USA, August 2014.
URTeC 1922263 Utilizing Ant-tracking to Identify Slowly Slipping Faults in the Barnett Shale Noha Sameh Farghal* and Mark D. Zoback, Stanford University, Stanford, CA, USA Copyright 2014, Unconventional
More informationP191 Bayesian Linearized AVAZ Inversion in HTI Fractured Media
P9 Bayesian Linearized AAZ Inversion in HI Fractured Media L. Zhao* (University of Houston), J. Geng (ongji University), D. Han (University of Houston) & M. Nasser (Maersk Oil Houston Inc.) SUMMARY A new
More informationMagnitude, scaling, and spectral signature of tensile microseisms
Magnitude, scaling, and spectral signature of tensile microseisms David W. Eaton Department of Geoscience, University of Calgary Summary The spatial dimensions and rupture characteristics of microseismic
More informationLocating Events using Borehole Microseismic Monitoring by Inclusion of Particle Motion Analysis
Locating Events using Borehole Microseismic Monitoring by Inclusion of Particle Motion Analysis Balikpapan, October 5th-8th, 2015 REXHA VERDHORA RY Geophysical Engineering Institut Teknologi Bandung ANDRI
More information2012 SEG SEG Las Vegas 2012 Annual Meeting Page 1
Wei Huang *, Kun Jiao, Denes Vigh, Jerry Kapoor, David Watts, Hongyan Li, David Derharoutian, Xin Cheng WesternGeco Summary Since the 1990s, subsalt imaging in the Gulf of Mexico (GOM) has been a major
More informationVelocity Measurements of Pore Fluids at Pressure and Temperature: Application to bitumen
Velocity Measurements of Pore Fluids at Pressure and Temperature: Application to bitumen Arif Rabbani 1*, Douglas R Schmitt 1, Jason Nycz 2, and Ken Gray 3 1 Institute for Geophysical Research, Department
More informationElements of 3D Seismology Second Edition
Elements of 3D Seismology Second Edition Copyright c 1993-2003 All rights reserved Christopher L. Liner Department of Geosciences University of Tulsa August 14, 2003 For David and Samantha And to the memory
More informationMicroseismic Monitoring Shale Gas Plays: Advances in the Understanding of Hydraulic Fracturing 20 MAR 16 HANNAH CHITTENDEN
Microseismic Monitoring Shale Gas Plays: Advances in the Understanding of Hydraulic Fracturing 20 MAR 16 HANNAH CHITTENDEN Introduction Early days: Microseismic monitoring has been around since the early
More informationTowards full waveform inversion: A torturous path
FWI Towards full waveform inversion: A torturous path J. Helen Isaac and Gary F. Margrave ABSTRACT Full waveform inversion (FWI) can be viewed as an iterative cycle involving forward modelling, pre-stack
More informationMEASUREMENT OF HYDRAULICALLY ACTIVATED SUBSURFACE FRACTURE SYSTEM IN GEOTHERMAL RESERVOIR BY USING ACOUSTIC EMISSION MULTIPLET-CLUSTERING ANALYSIS
MEASUREMENT OF HYDRAULICALLY ACTIVATED SUBSURFACE FRACTURE SYSTEM IN GEOTHERMAL RESERVOIR BY USING ACOUSTIC EMISSION MULTIPLET-CLUSTERING ANALYSIS HIROKAZU MORIYA 1, HIROAKI NIITSUMA 1 and ROY BARIA 2
More informationIwan Yandika Sihotang, Tommy Hendriansyah, Nanang Dwi Ardi
Proceedings Indonesia International Geothermal Convention & Exhibition 2014 Jakarta Convention Center, Indonesia 4-6 June 2014 FOCAL MECHANISM FOR DETERMINING FAULT PLANES ORIENTATION IN GAMMA GEOTHERMAL
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 informationNumerical Study: Time-Reversed Reciprocal Method and Damage Detection Method for Weld Fracture
Chapter 4 Numerical Study: Time-Reversed Reciprocal Method and Damage Detection Method for Weld Fracture A numerical study is performed to gain insight into applying the proposed method of detecting high-frequency
More informationGeothermal Reservoir Imaging Using 2016 Walkaway VSP Data from the Raft River Geothermal Field
PROCEEDINGS, 44th Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 11-13, 2019 SGP-TR-214 Geothermal Reservoir Imaging Using 2016 Walkaway VSP Data from
More informationF021 Detetection of Mechanical Failure During Hyraulic Fracturing Through Passive Seismic Microseismic Monitoring
F021 Detetection of Mechanical Failure During Hyraulic Fracturing Through Passive Seismic Microseismic Monitoring A. De La Pena* (Microseismic Inc.), L. Eisner (Microseismic Inc.), M.P. Thornton (Microseismic
More informationOverview Avoiding Cycle Skipping: Model Extension MSWI: Space-time Extension Numerical Examples
Overview Avoiding Cycle Skipping: Model Extension MSWI: Space-time Extension Numerical Examples Guanghui Huang Education University of Chinese Academy of Sciences, Beijing, China Ph.D. in Computational
More informationSeismic reservoir characterisation
Seismic reservoir characterisation Unconventional reservoir (shale gas) Robert Porjesz 1 2014 B A K E R H U G H E S I N C O R P O R A TED. A LL R I G H TS R E S E R V E D. TERMS A N D C O N D I TI O N
More informationElastic least-squares reverse time migration
CWP-865 Elastic least-squares reverse time migration Yuting Duan, Paul Sava, and Antoine Guitton Center for Wave Phenomena, Colorado School of Mines ABSTRACT Least-squares migration (LSM) can produce images
More informationResearch Project Report
Research Project Report Title: Prediction of pre-critical seismograms from post-critical traces Principal Investigator: Co-principal Investigators: Mrinal Sen Arthur Weglein and Paul Stoffa Final report
More informationWaveform inversion for attenuation estimation in anisotropic media Tong Bai & Ilya Tsvankin Center for Wave Phenomena, Colorado School of Mines
Waveform inversion for attenuation estimation in anisotropic media Tong Bai & Ilya Tsvankin Center for Wave Phenomena, Colorado School of Mines SUMMARY Robust estimation of attenuation coefficients remains
More informationOverview of moment-tensor inversion of microseismic events
Overview of moment-tensor inversion of microseismic events Thomas S. Eyre 1 and Mirko van der Baan 1 Downloaded 08/11/15 to 142.244.191.52. Redistribution subject to SEG license or copyright; see Terms
More informationFull Waveform Inversion (FWI) with wave-equation migration. Gary Margrave Rob Ferguson Chad Hogan Banff, 3 Dec. 2010
Full Waveform Inversion (FWI) with wave-equation migration (WEM) and well control Gary Margrave Rob Ferguson Chad Hogan Banff, 3 Dec. 2010 Outline The FWI cycle The fundamental theorem of FWI Understanding
More informationThomas Bayes versus the wedge model: An example inference using a geostatistical prior function
Thomas Bayes versus the wedge model: An example inference using a geostatistical prior function Jason M. McCrank, Gary F. Margrave, and Don C. Lawton ABSTRACT The Bayesian inference is used to estimate
More informationAnalytic Tools for Quantitative Interpretation. Mrinal K.Sen
Analytic Tools for Quantitative Interpretation Mrinal K.Sen Tools Modeling of fractures static, dynamic, numerical Modeling seismic wave propagation in anisotropic media Travel time an amplitude based
More informationOpportunities for Geothermal Development Created by New Technologies S2-1-2
Opportunities for Geothermal Development Created by New Technologies S2-1-2 Geothermal Academy Masami Nakagawa Director of Geothermal Academy Fulbright Energy Specialist CSM/NREL Joint Appointee Department
More informationDownloaded 11/08/14 to Redistribution subject to SEG license or copyright; see Terms of Use at
Velocity tomography based on turning waves and finite-frequency sensitivity kernels Xiao-Bi Xie*, University of California at Santa Cruz, Jan Pajchel, Statoil, Norway Summary We present a velocity tomography
More informationSeismic techniques for imaging fractures, cracks and faults in the Earth. Michael Kendall
Seismic techniques for imaging fractures, cracks and faults in the Earth Michael Kendall Issues and Challanges Geometry (aspect ratio, size, orientation, density) Non-uniqueness (e.g., single set of aligned
More informationA new methodology for estimating field horizontal stress from microseismic focal mechanisms
A new methodology for estimating field horizontal stress from microseismic focal mechanisms Alireza Agharazi 1*, Peter Duncan 1 and Michael Thornton 1 introduce a new methodology for the estimation of
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 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 informationResearch Themes in Stimulation Geomechanics. How do we optimize slickwater frac ing?
Research Themes in Stimulation Geomechanics How do stress, fractures and rock properties affect the success of stimulation? How do pressure and stress (and formation properties) evolve during stimulation?
More informationSimultaneous estimation of wavefields & medium parameters
Simultaneous estimation of wavefields & medium parameters reduced-space versus full-space waveform inversion Bas Peters, Felix J. Herrmann Workshop W- 2, The limit of FWI in subsurface parameter recovery.
More informationBandlimited impedance inversion: using well logs to fill low frequency information in a non-homogenous model
Bandlimited impedance inversion: using well logs to fill low frequency information in a non-homogenous model Heather J.E. Lloyd and Gary F. Margrave ABSTRACT An acoustic bandlimited impedance inversion
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 informationPengcheng Fu, Yue Hao, and Charles R. Carrigan
Pengcheng Fu, Yue Hao, and Charles R. Carrigan Math, Science and Computation of Hydraulic Fracturing Workshop Stanford, March 21, 2013 This work was performed under the auspices of the U.S. Department
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 informationSeismic methods in heavy-oil reservoir monitoring
Seismic methods in heavy-oil reservoir monitoring Duojun A. Zhang and Laurence R. Lines ABSTRACT Laboratory tests show that a significant decrease in acoustic velocity occurs as the result of heating rock
More informationCoupled seismoelectric wave propagation in porous media. Mehran Gharibi Robert R. Stewart Laurence R. Bentley
Coupled seismoelectric wave propagation in porous media Mehran Gharibi Robert R. Stewart Laurence R. Bentley a Introduction Seismic waves induce electric and magnetic fields. Conversion of acoustic energy
More informationIn situ stress estimation using acoustic televiewer data
Underground Mining Technology 2017 M Hudyma & Y Potvin (eds) 2017 Australian Centre for Geomechanics, Perth, ISBN 978-0-9924810-7-0 https://papers.acg.uwa.edu.au/p/1710_39_goodfellow/ SD Goodfellow KORE
More informationRadiation pattern in homogeneous and transversely isotropic attenuating media
Radiation pattern in homogeneous and transversely isotropic attenuating media Satish Sinha*, Sergey Abaseyev** and Evgeni Chesnokov** *Rajiv Gandhi Institute of Petroleum Technology, Rae Bareli, UP 229010
More informationMicrodeformation: combining direct fracture height measurement with microseismic response Natalia Verkhovtseva, Greg Stanley
Microdeformation: combining direct fracture height measurement with microseismic response Natalia Verkhovtseva, Greg Stanley Principal Geophysicist Pinnacle a Halliburton service Outline Introduction Microseismic
More informationAdvanced Workshop on Evaluating, Monitoring and Communicating Volcanic and Seismic Hazards in East Africa.
2053-11 Advanced Workshop on Evaluating, Monitoring and Communicating Volcanic and Seismic Hazards in East Africa 17-28 August 2009 Seismic monitoring on volcanoes in a multi-disciplinary context Jürgen
More informationExtending the magnitude range of seismic reservoir monitoring by Utilizing Hybrid Surface Downhole Seismic Networks
Extending the magnitude range of seismic reservoir monitoring by Utilizing Hybrid Surface Downhole Seismic Networks Gisela Viegas*, ESG, Kingston, Ontario, Canada Gisela.Fernandes@esgsolutions.com and
More informationA damage model for fracking
Article A damage model for fracking J Quinn Norris 1, Donald L Turcotte 2 and John B Rundle 1,2,3 arxiv:1503.01703v1 [physics.geo-ph] 5 Mar 2015 Abstract Injections of large volumes of water into tight
More informationUse of S-wave attenuation from perforation shots to map the growth of the stimulated reservoir volume in the Marcellus gas shale
Use of S-wave attenuation from perforation shots to map the growth of the stimulated reservoir volume in the Marcellus gas shale Yunhui Tan 1, Chengping Chai 1, and Terry Engelder 1 Downloaded 10/07/14
More informationSEISMIC TOMOGRAPHY MEASUREMENTS IN THE RADIOACTIVE WASTE DISPOSAL SITE IN BÁTAAPÁTI
Seismic tomography measurements in Bátaapáti SEISMIC TOMOGRAPHY MEASUREMENTS IN THE RADIOACTIVE WASTE DISPOSAL SITE IN BÁTAAPÁTI 1. Introduction PETER TILDY 1, ZSOLT PRONAY 1, ENDRE TÖRÖS 1 As a part of
More informationEstimating energy balance for hydraulic fracture stimulations: Lessons Learned from Basel
Estimating energy balance for hydraulic fracture stimulations: Lessons Learned from Basel David W. Eaton*, Department of Geoscience, University of Calgary, Calgary, Canada eatond@ucalgary.ca and Neda Boroumand,
More informationInduced seismicity in Enhanced Geothermal Systems (EGS) in Alsace, France. Jean Schmittbuhl 1
Induced seismicity in Enhanced Geothermal Systems (EGS) in Alsace, France Jean Schmittbuhl 1 Olivier Lengliné 1, Mohamed Boubacard 1, Lucie Lamourette 1, Lilas Vivin 1, Nicolas Cuenot 2 1 EOST-IPGS, Université
More informationDETAILED IMAGE OF FRACTURES ACTIVATED BY A FLUID INJECTION IN A PRODUCING INDONESIAN GEOTHERMAL FIELD
PROCEEDINGS, Thirty-Fourth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 9-11, 2009 SGP-TR-187 DETAILED IMAGE OF FRACTURES ACTIVATED BY A FLUID INJECTION
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 information7.2.1 Seismic waves. Waves in a mass- spring system
7..1 Seismic waves Waves in a mass- spring system Acoustic waves in a liquid or gas Seismic waves in a solid Surface waves Wavefronts, rays and geometrical attenuation Amplitude and energy Waves in a mass-
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 informationGeophysical model response in a shale gas
Geophysical model response in a shale gas Dhananjay Kumar and G. Michael Hoversten Chevron USA Inc. Abstract Shale gas is an important asset now. The production from unconventional reservoir like shale
More informationApproximate- vs. full-hessian in FWI: 1D analytical and numerical experiments
Approximate- vs. full-hessian in FWI: 1D analytical and numerical experiments Raul Cova and Kris Innanen ABSTRACT Feasibility of using Full Waveform Inversion (FWI) to build velocity models has been increasing
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 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 information