Monitoring and Characterization Based Upon Geophysical Onset Times*

Similar documents
ACOUSTIC LOGGING GUIDED WAVES IN TRANSVERSELY ISOTROPIC FORMATIONS

the results to larger systems due to prop'erties of the projection algorithm. First, the number of hidden nodes must

Main Menu. SEG Houston 2009 International Exposition and Annual Meeting. Summary

Kinetic Model Completeness

Fall 2013 Physics 172 Recitation 3 Momentum and Springs

Principles of Carbonate Reservoirs

A Few Basic Facts About Isothermal Mass Transfer in a Binary Mixture

UN Committee of Experts on Environmental Accounting New York, June Peter Cosier Wentworth Group of Concerned Scientists.

Petrel TIPS&TRICKS from SCM

CAUSAL INFERENCE. Technical Track Session I. Phillippe Leite. The World Bank

A Scalable Recurrent Neural Network Framework for Model-free

ChE 471: LECTURE 4 Fall 2003

WRITING THE REPORT. Organizing the report. Title Page. Table of Contents

ECEN 4872/5827 Lecture Notes

Bootstrap Method > # Purpose: understand how bootstrap method works > obs=c(11.96, 5.03, 67.40, 16.07, 31.50, 7.73, 11.10, 22.38) > n=length(obs) >

Math Foundations 20 Work Plan

ENSC Discrete Time Systems. Project Outline. Semester

History the Hood Way. Amy Shell-Gellasch Betty Mayfield Hood College. MD-DC-VA Section October 27, 2012

Rangely RE 4 Curriculum Development 5 th Grade Mathematics

5 th grade Common Core Standards

NUMBERS, MATHEMATICS AND EQUATIONS

Admissibility Conditions and Asymptotic Behavior of Strongly Regular Graphs

Differentiation Applications 1: Related Rates

Floating Point Method for Solving Transportation. Problems with Additional Constraints

A Simple Set of Test Matrices for Eigenvalue Programs*

Revision: August 19, E Main Suite D Pullman, WA (509) Voice and Fax

Study Group Report: Plate-fin Heat Exchangers: AEA Technology

Computational modeling techniques

Chapter 3: Cluster Analysis

Optimization Programming Problems For Control And Management Of Bacterial Disease With Two Stage Growth/Spread Among Plants

B. Definition of an exponential

5 th Grade Goal Sheet

1 The limitations of Hartree Fock approximation

CHAPTER 24: INFERENCE IN REGRESSION. Chapter 24: Make inferences about the population from which the sample data came.

NAME TEMPERATURE AND HUMIDITY. I. Introduction

Lecture 7: Damped and Driven Oscillations

February 28, 2013 COMMENTS ON DIFFUSION, DIFFUSIVITY AND DERIVATION OF HYPERBOLIC EQUATIONS DESCRIBING THE DIFFUSION PHENOMENA

CARBONATED WATER INJECTION FOR EOR IN ONE DIMENSIONAL FLOW WITH CONSTANT PRESSURE BOUNDARIES

Unit 2 Expressions, Equations, and Inequalities Math 7

Chapter 23 Electromagnetic Waves Lecture 14

Mathematics and Computer Sciences Department. o Work Experience, General. o Open Entry/Exit. Distance (Hybrid Online) for online supported courses

[COLLEGE ALGEBRA EXAM I REVIEW TOPICS] ( u s e t h i s t o m a k e s u r e y o u a r e r e a d y )

A NOTE ON THE EQUIVAImCE OF SOME TEST CRITERIA. v. P. Bhapkar. University of Horth Carolina. and

Department: MATHEMATICS

, which yields. where z1. and z2

Modeling the Nonlinear Rheological Behavior of Materials with a Hyper-Exponential Type Function

8 th Grade Math: Pre-Algebra

Foundations for Dynamic and Sensitive Equipment

GEOL 2920C The Sedimentary Rock Cycle of Mars & Earth Eberswalde fan deposits: deltaic or alluvial? Deltas at Aeolis Dorsa(?

WELL LOGS AND ROCK PHYSICS IN SEISMIC RESERVOIR CHARACTERIZATION

1. INTRODUCTION. In many polymer processing operations, molten polymers. emerge from dies into a stress field which deforms the

A.H. Helou Ph.D.~P.E.

The calculation method of small-scale water injection multiple in water drive reservoirs

Function notation & composite functions Factoring Dividing polynomials Remainder theorem & factor property

Revisiting the Socrates Example

EASTERN ARIZONA COLLEGE Precalculus Trigonometry

Associated Students Flacks Internship

Introductory Thoughts

Negative Propagating Light Daniel Bauer University of Arizona

Analysis on the Stability of Reservoir Soil Slope Based on Fuzzy Artificial Neural Network

This section is primarily focused on tools to aid us in finding roots/zeros/ -intercepts of polynomials. Essentially, our focus turns to solving.

Compressibility Effects

Building to Transformations on Coordinate Axis Grade 5: Geometry Graph points on the coordinate plane to solve real-world and mathematical problems.

Nye County Groundwater Evaluation Drilling Program. Levi Kryder Nye County NWRPO May 5, 2011

Find this material useful? You can help our team to keep this site up and bring you even more content consider donating via the link on our site.

5 th Grade Goal Sheet

Theoretical study of third virial coefficient with Kihara potential

NOTE ON THE ANALYSIS OF A RANDOMIZED BLOCK DESIGN. Junjiro Ogawa University of North Carolina

A new Type of Fuzzy Functions in Fuzzy Topological Spaces

Lifting a Lion: Using Proportions

Systems and Modeling Introduction

Free Vibrations of Catenary Risers with Internal Fluid

Thermodynamics EAS 204 Spring 2004 Class Month Day Chapter Topic Reading Due 1 January 12 M Introduction 2 14 W Chapter 1 Concepts Chapter 1 19 M MLK

Numerical Simulation of the Thermal Resposne Test Within the Comsol Multiphysics Environment

Design and Simulation of Dc-Dc Voltage Converters Using Matlab/Simulink

On Boussinesq's problem

Modelling of Clock Behaviour. Don Percival. Applied Physics Laboratory University of Washington Seattle, Washington, USA

Competency Statements for Wm. E. Hay Mathematics for grades 7 through 12:

Displacement and Deflection Sensitivity of Gas-coupled Laser Acoustic. Detector

Dataflow Analysis and Abstract Interpretation

ES201 - Examination 2 Winter Adams and Richards NAME BOX NUMBER

A Transition to Advanced Mathematics. Mathematics and Computer Sciences Department. o Work Experience, General. o Open Entry/Exit

Emphases in Common Core Standards for Mathematical Content Kindergarten High School

5.4 Measurement Sampling Rates for Daily Maximum and Minimum Temperatures

Resampling Methods. Cross-validation, Bootstrapping. Marek Petrik 2/21/2017

ECE 2100 Circuit Analysis

Thermodynamics Partial Outline of Topics

Statistics, Numerical Models and Ensembles

NGSS High School Physics Domain Model

The standards are taught in the following sequence.

A TIME-DOMAIN FINITE-DIFFERENCE METHOD WITH ATTENUATION BY A RECURSIVE ALGORITHM

Biocomputers. [edit]scientific Background

Figure 1a. A planar mechanism.

Chem 115 POGIL Worksheet - Week 8 Thermochemistry (Continued), Electromagnetic Radiation, and Line Spectra

4F-5 : Performance of an Ideal Gas Cycle 10 pts

Intelligent Pharma- Chemical and Oil & Gas Division Page 1 of 7. Global Business Centre Ave SE, Calgary, AB T2G 0K6, AB.

Lead/Lag Compensator Frequency Domain Properties and Design Methods

7 TH GRADE MATH STANDARDS

Thermodynamics and Equilibrium

Web-based GIS Systems for Radionuclides Monitoring. Dr. Todd Pierce Locus Technologies

Transcription:

Mnitring and haracterizatin Based Upn Gephysical Onset Times* Dn Vasc 1 Search and Discvery Article #4237 (217)** Psted March 2, 217 *Adapted frm ral presentatin given at AAPG Internatinal nference & xhibitin, Barcelna Spain, April 3-6, 216 **Datapages 217 Serial rights given by authr. Fr all ther rights cntact authr directly. 1 Lawrence Berkeley Natinal Labratry (rdica_negulescu@htmail.cm) Abstract Gephysical methds are increasingly called upn t mnitr fluid injectin and prductin. Gephysical mnitring f flw related prcesses is hampered by difficulties in relating changes in fluid saturatin and pressure t bservable quantities such as seismic amplitude changes, variatins in electric fields, and surface defrmatin. Prblems are particularly acute when relating changes in the amplitudes f bservable gephysical quantities t the magnitudes f fluid related changes in a reservir. In such cases the relatinship depends strngly upn the rck physics mdel that is chsen. In this presentatin I present an alternative apprach t reservir mnitring and characterizatin that is based n the ntin f an nset time. An nset time is the calendar time at which a measured quantity begins t deviate frm its backgrund r initial value. In many cases invlving injectin and prductin, the nset time is related t the arrival f a fluid frnt r a rapid increase in pressure. Under the relatively general cnditin that a change in fluid saturatin and/r pressure leads t a change in a gephysical attribute, it can be shwn that nset times are sensitive t flw prperties, such as permeability, and nt sensitivity f the details f the rck physics mdel. Several examples f the use f nset times fr mnitring and characterizatin will be given, invlving bth gedetic and seismic bservatins. One example invlves the gedetic mnitring f grund defrmatin assciated with the gelgical strage f injected carbn dixide at In Salah, Algeria. Seismic timelapse mnitring f injected carbn dixide at the Fri site in Texas prvides yet anther example. References ited Vasc, D.W., A. Bakulin, H. Baek, and L.R. Jhnsn, 215, Reservir characterizatin based upn the nset f time-lapse amplitude changes: Gephysics, v. 8, January-February. Vasc, D.W., T.M. Daley, and A. Bakulin, 214, Utilizing the nset f time-lapse changes: a rbust basis fr reservir mnitring and characterizatin: Gephysical Jurnal Internatinal.

Presenter s ntes: Our gal is a very cmplete reservir characterizatin/making use f all majr data types. mphasis n practicality f the algrithms. It is imprtant t include all available reservir engineering data as this is directly sensitive t reservir flw prperties. Field testing is imprtant t ensure that ur methds are truly effective. Finally, we need t get these tls int the hands f industry.

Issues: Wuld like t estimate prperties such as permeability Gephysical data indirectly related t flw prperties ven time-lapse bservatins depend upn a rck physics mdel AAPG/SG Internatinal nference and xhibitin 216, 3-6 April

Apprach: Avid using the magnitude f gephysical changes due t variatins in the state [saturatin, pressure, stress] f the reservir Use the nset time f changes in gephysical bservatins AAPG/SG Internatinal nference and xhibitin 216, 3-6 April

Outline: Onset times Illustratin based upn time-lapse seismic amplitude changes Inversin fr reservir permeability nclusins AAPG/SG Internatinal nference and xhibitin 216, 3-6 April

Definitin: Onset time - The time at which a seismic bservatin begins t change frm its backgrund value. AAPG/SG Internatinal nference and xhibitin 216, 3-6 April

Seismic Time-Lapse Reflectin Amplitudes AAPG/SG Internatinal nference and xhibitin 216, 3-6 April

Apprach: Use reservir mdel t cmpute seismic amplitude changes due t O 2 injectin mpare inversin based upn seismic amplitudes with ne based upn nset times AAPG/SG Internatinal nference and xhibitin 216, 3-6 April

AAPG/SG Internatinal nference and xhibitin 216, 3-6 April

1.RO- 1.83 Km..c 1: z <D U L 2 '" 75. 288 56.5 BIV3 1123 r:tt AJ\B F.AS"'" OJ} 2 24O PI"':S.5UII: (M -)A) ItH, 1.83-1.86 Km 1.86-1.89 Km 185.<: 1: Z 13 " U l "' '"".<: Z <D U -;;; l 2 2 75 13 185 24 2 75 1:S 185 2 1 Distance ast (m). Arrival -ime (cays) 2O. AAPG/SG Internatinal nference and xhibitin 216, 3-6 April

. IX=3, IY=4.1 IX=4, IY=4 5. -.---... 1. - "----' I f D- w 15. - I I 2. i=======;::=:-l== L I I (j) I- Z => U w.5 U.u => t: -,5 --' L :::;; <:(.1 U :::;; {I) -.15 W (f) -.2 +------.-----,------.-----,----1 1. 1.1 1. 2 1.3 1.4,.5G TIM (5) 25. -------'-I------'I--------r- I ------,LI---- 2. 3. 4. 5. 6. 7. OMPRSS IONA L V LOITY (KMjS) AAPG/SG Internatinal nference and xhibitin 216, 3-6 April

Issues: The relatinship between the reservir state and seismic velcities is cmplicated Thught t depend upn messcpic [cm t m] variatins in saturatin, pressure, and mechanical prperties AAPG/SG Internatinal nference and xhibitin 216, 3-6 April

Seismic velcities depend upn hw fluids are distributed and averaged by the prpagating wave. K Vigt = 1 = K Reuss 3 ås K i i i=1 3 å S i K i=1 i Vigt average Reuss average K = 1 ( Hill 2 K + K ) Hill average Vigt Reuss AAPG/SG Internatinal nference and xhibitin 216, 3-6 April

3.4 (j)..y. '--" >- -+- U Q) > 3.3 3.2 Reuss (j) (j) (]) I..... U 3. 1..2.4.6.8 1. Sa tura ti n (f rlct i n ) AAPG/SG Internatinal nference and xhibitin 216, 3-6 April

Hi li 24-.-----------------------------------, 185 '--./ -. +- '-- c 13 Q) U +- (/) 75 2--------._--------._--------._------ 2 75 13 185 24 Dista nce e a s t (m ) Amplitude c ha n ge (fractin) 1. AAPG/SG Internatinal nference and xhibitin 216, 3-6 April

Re uss 24-.--------------======----------, V igt 24-.--------------------------------, r 185 185 -"..c L L 13 13O OJ OJ U U III 75 III 75!J O D D DD!: 2--------------_,r_-------------- 2 --------------_,r_------------ 2 75 13 185 24 2 75 13 185 24 Distance east ( m ) Distance east ( m ) Amplitude change (fractin) 1. Amplitude change (fractin) 1. AAPG/SG Internatinal nference and xhibitin 216, 3-6 April

5 5,.-." 4 - U1 V 3 - '-' <l> - -+- 2 - - «1 - A, nset ', I \ I I \ I.\..... '- \ -.-...., -...- "'- ++,-,-,,, r, I I I --... i'2!._, \"iill r " Reuss 4-3 - 2-1 -,, '. nset.,,.,",, I I I I,\, " +... Vigt, ' \ '...,.' -\ + +'t! Reuss 18 36 54 72 9 a 18 36 54 72 9 S u rvey t im e ( d ays ) B, S u rvey time (days) AAPG/SG Internatinal nference and xhibitin 216, 3-6 April

Re u ss 24-r----------------- V igt 24,------------------ '--../ 185 '--../ 185 L: L: "-- Z 13 "-- Z 1.3 Q) U (!) 75 Q) U (!) 75 2-r------------_,--------r_---- 2 +-------r--------,---------,,--------i 2 75 Dista n ce 13 a st 185 (m) 24 2 75 Distance 13 ast 185 (m) 24 Onset time (days) 9 Onset ti m e (days) 9 AAPG/SG Internatinal nference and xhibitin 216, 3-6 April

Iterative inversin AAPG/SG Internatinal nference and xhibitin 216, 3-6 April

Presenter s ntes: (1) In ur apprach we intrduce a new methd fr the matching f reservir engineering data, based upn streamline simulatin. (2) With this new methd we can develp efficient ways t integrate prductin and seismic data. (3) Field applicatins are an imprtant part f ur apprach as they indicate what will be feasible.

Iterative Updating Algrithm Iterative mdificatin f the reservir mdel t better match nset times: K m = K m-1 +DK m That is, t better match bservatins in a least squares sense: Misfit = ( dt - AdK ) T ( dt - AdK ) m m quivalent t slving the linear system: A T AdK m = A T dt AAPG/SG Internatinal nference and xhibitin 216, 3-6 April

Influence f the rck physics mdel AAPG/SG Internatinal nference and xhibitin 216, 3-6 April

Time-Lapse Amplitudes d F ( K) Hill Hill dreuss FRe uss ( K) d F ( K) Vigt Vigt Hill Inversin rrect Sensitivities T Reuss Inversin Vigt Inversin A A K A d Hill Hill m Hill Hill T T Hill Hill m Hill Reuss A A K A d T A A K A d Hill Hill m Hill Vigt T T AAPG/SG Internatinal nference and xhibitin 216, 3-6 April

1 Amplitude Inversin.-+- -.- 8 6 4 2... "- '", Reuss... -..., "a. -. -..., Vigt Hill 'a._ Hill Reuss _ Vigt + --... ----- - - --.--- +"t I +-+++-t--t--t--+-+-i-+ 5 1 Iterat i n 15 AAPG/SG Internatinal nference and xhibitin 216, 3-6 April

1.8-1.83 Km Hil l 24 24 185 ::: lase.r:: 1: s: Z 13e- Z 13 " c " c Q> " Vl 7Stl 75 '" 2 2 2:) 75 13 185 24 2 75 13 185 24O Reuss Vi t 2>1 24 O. 185 ISSG.r:: -<= 1: L Z 13 Z Ql Q) U " Ui i 75 Vl 13 75 7 ------,-------,-------,------- 2 ------.. ------.. ------.. ------ 2 75 13 185 24 75 13 le Distance a st (m) 24-1. Lg 1K Mu ti plier 1. AAPG/SG Internatinal nference and xhibitin 216, 3-6 April

Onset times d F ( K) Hill Hill dreuss FRe uss ( K) d F ( K) Vigt Vigt Hill Inversin rrect Sensitivities T Reuss Inversin Vigt Inversin A A K A d Hill Hill m Hill Hill T T Hill Hill m Hill Reuss A A K A d T A A K A d Hill Hill m Hill Vigt T T AAPG/SG Internatinal nference and xhibitin 216, 3-6 April

Onset Time Inversin 2 -;;- + Q) I- ::J 16 - - U) Hill Reuss _ Vigt + Vl > 12-8 - x -+ '+ Vl 4 - I 2,, I 4 6 Itera tin I.L.L I 8.L 1 AAPG/SG Internatinal nference and xhibitin 216, 3-6 April

1.8-1.83 Km Hi ll 24) 24 185J '85.c 1:: '- '- z 13:> Z ' 3 " OJ <.l <J Vl 7S Vl '" 75:) lbse." 2 2J 2 75 13 185 24 2 75 13GO 185 2/ Reu ss V igt 24 24.. '\ c 15 '- r.c -;: L Z 13 Z 13 (1J w U U Vl 7.') Vl 75 2 -------,-------,-------,------ 2 75 13 18 24 2 1-------,-------r------,--------l :.1 u lj.ou -1. Lg ' U K \lultiplle' 1. AAPG/SG Internatinal nference and xhibitin 216, 3-6 April

Variatin in the time inverval between seismic surveys AAPG/SG Internatinal nference and xhibitin 216, 3-6 April

211 15 Da y Interval 6 Day Interva l 185 '" L Z 13 w U 75 2':) +----,-----,----,------1 2+----,-----,----,------1 2ec' 75 13 185 24 2 75 13 185 24' 24 12 Day Interva l 3 Day Inter va l ll:io Ie L Z 13 r c -;;; 75 2 +----,,----,----,------1 2+----,-----,----,------1 2e 7 5 13 1.35 24 2 75 13 185 2 4 [,::1 Distance ast (m ) Dista nce ast (m) On s et time (days ) 9 Onset time (days) 9 AAPG/SG Internatinal nference and xhibitin 216, 3-6 April

1 5 Day I nt e r va I 6 Day Int e rval 24 24 185 185 "" t Z ID 13 13 U c 75 75 2 2 2 75 13 185 24 2 75 13 185 24 1 2 Day In t e rvl I 3 Day I nt e rva I?.ir, 24, I SSG 1 gs "" Z 13 13 v u c. w 75 75 I I 2 2 2 7.'iO I.'O 185 24 2 75 13 185 24 Distance ast (m) Distance ast (m) -1.:) Lg 1K Multiplier 1. -1. Lg 1K Multiplier 1, AAPG/SG Internatinal nference and xhibitin 216, 3-6 April

nclusins Relatinship between changes in the state f a reservir and gephysical bservatins can be cmplicated and cnfunded by prly knwn equatin-f-state parameters Inversin based upn nset times appear stable with respect t variatins in the saturatin/velcity relatinship AAPG/SG Internatinal nference and xhibitin 216, 3-6 April

Acknwlegments This material is based upn wrk supprted by the U. S. Department f nergy, Office f Science, Office f Basic nergy Sciences, hemical Sciences, Gesciences, and Bisciences Divisin under cntract number D-A2-5-H11231. Additinal supprt was prvided by Aramc Research. AAPG/SG Internatinal nference and xhibitin 216, 3-6 April