SEISMIC INVERSION OVERVIEW

Save this PDF as:
 WORD  PNG  TXT  JPG

Size: px
Start display at page:

Download "SEISMIC INVERSION OVERVIEW"

Transcription

1 DHI CONSORTIUM SEISMIC INVERSION OVERVIEW Rocky Roden September 2011

2 NOTE: Terminology for inversion varies, depending on the different contractors and service providers, emphasis on certain approaches, and the goals to be accomplished. There are numerous combinations, hybrids, and variations of these inversion methods. This overview attempts to summarize the most common terminology and inversion approaches described in the industry.

3 Seismic Inversion Outline Inversion Definition and Terms Required Data Types of Seismic Inversion Role of Inversion Realistic Expectation from Inversion Inversion and New SEC Regulations

4 What is Inversion? Inversion transforms seismic reflection data into rock and fluid properties. The objective of seismic inversion is to convert reflectivity data (interface properties) to layer properties. To determine elastic parameters, the reflectivity from AVO effects must be inverted. The most basic inversion calculates acoustic impedance (density X velocity) of layers from which predictions about lithology and porosity can be made. The more advanced inversion methods attempt to discriminate specifically between lithology, porosity, and fluid effects. Inversions can be grouped into categories: pre-stack vs. post-stack, deterministic vs. geostatistical, or relative vs. absolute. Inversion is the flip side of forward modeling.

5 Common Inversion Terms Relative Impedance Relative changes in inversion, not real values. Absolute Impedance Actual acoustic impedance, contains low frequency trend. Deterministic Inversion A single output value is determined from the input. Stochastic (Probabilistic) Inversion A range of equally probable outputs are derived from the input (under controlled number of inputs). Objective Function A quantitative measure of the misfit between the observed data and the data predicted by inversion (L1 and L2 Norm). L2 Norm Most common measure of data misfit, a least squares difference (sum of the square of the data residuals). L1 Norm Misfit measure that is the sum of the absolute data residuals. Global More than one trace is inverted at the same time within a common objective function. The number of traces inverted at once depends on the algorithm. Simulated Annealing A global optimization technique based on crystal growth in a cooling volcanic melt. Computes difference between seismic and convolutional model with applied wavelet. The model is perturbed and a new model is simulated, and the difference is again measured. Differences are compared with the smallest difference accepted and the other differences accepted with a probability (Metropolis Criterion). The process is repeated until and there is a very small residual or a threshold has been reached.

6 Forward Modeling takes well logs, combines with wavelet to produce synthetic seismic trace Inversion takes seismic trace, removes effects of estimated wavelet and creates acoustic impedance values Oilfield Review, 2008

7 Data Considerations for Inversion Pre-stack and/or post-stack time migrated seismic data Well data conditioning and editing is necessary Wavelet estimation required for all modern inversion approaches -Deterministic -Statistical Low Frequency Trend absolute acoustic impedance contains a low frequency trend that must be obtained usually from well control or stacking velocities NOTE: Without wavelet estimation and well calibration, the inversion solution is non-unique.

8 Data Considerations for Inversion Pillar, 2011

9 Data Considerations for Inversion Low Frequency Trend in Red Oilfield Review, 2008

10 Data Considerations for Inversion Missing low frequencies and their modeling from well log data cause bias in inversion, regardless of the methodology. Determining low frequency issues: Lithologies encountered at the wells thin or thicken away from the wells. In fact, some lens-like possible plays between wells are often difficult or impossible to model accurately. Lateral geological changes away from the wells have no resemblance to the geological layers at wells used fro low frequency model generation, e.g., sand injections. There are excessive lateral changes in rock and reservoir properties (e.g., porosity changes and shaling out). Ozdemir, 2009

11 Data Considerations for Inversion Rock physics approach to well log preparation for inversion: Continuous reconstruction of formation properties along the wellbore, both for reservoirs and non-reservoirs Characterization of rocks in terms of their elastic properties Establishment of causal relationships between elastic and petrophysical properties such as porosity, clay content and fluid saturation NOTE: Depends on inversion approach and stage of exploration through development cycle.

12 Types of Inversion Running Sum Inversion like (Runsum) Recursive Trace integration (relative acoustic impedance) Colored Inversion Sparse Spike (CSSI) Model-Based Inversion AVO Inversion Elastic Impedance Extended Elastic Impedance Simultaneous Inversion Stochastic Inversion Geostatistical Bayesian

13 Running Sum Running Sum is the integration or adding of seismic amplitudes on a seismic trace. It is the cumulative sum of amplitude values at any time sample starting from the top to the end of the trace. Running Sum provides an approximation to acoustic impedance (velocity X density), while lacking the low frequency trend of a full inversion this can be useful over short intervals. The less noise and the broader the bandwidth (narrower the wavelet), the better will be the approximation to impedance.

14 Running Sum Normal Amplitude Running Sum Courtesy SMT

15 Recursive Trace Integration Recursive/Trace Integration (RTI) employs the discrete recursive inversion formula (below) which indicates that the acoustic impedance of a particular layer and the reflection coefficient at its base can be used to calculate the acoustic impedance of the next layer. An estimate of the acoustic impedance of the first layer is required. Normal incidence reflectivity Recursive inversion formula r = reflection coefficient Z = impedance

16 Recursive Trace Integration Seismic section from Alberta highlighting an anomalous bright spot zone. Recursive inversion of the section above Russell and Hampson, 2006

17 Colored Inversion Colored Inversion transforms migrated seismic data into a bandlimited acoustic impedance volume by shaping the mean seismic spectrum into the impedance log spectrum. This approach only gives a relative acoustic impedance which does not contain any low frequency components. This inversion method can make it simpler to interpret the variations in thickness related to various lithologic packages, properties of thin beds, and changes in acoustic impedance associated with fluid effects and rock property variations. Veeken and Da Silva, 2004

18 Colored Inversion Colored Inversion was first introduced by Lancaster and Whitcombe (2000) and involves deriving an operator that transforms the seismic amplitude spectrum to the acoustic impedance spectrum from well logs with a -90 phase change. This typically requires crossplotting the log of the well acoustic impedances versus log of amplitudes to derive a best-fit line to the data. The slope of this line is used to develop an operator to apply to the seismic data to convert to acoustic impedance. The logic employed is that the gross spectral properties from the well acoustic impedances in any given field are reasonably constant. A phase rotation may be applied to the seismic data, but this approach assumes the seismic is zero phase. Colored inversion may establish a base case to compare against more sophisticated inversion techniques Courtesy SMT

19 Colored Inversion Conventional Seismic Line Colored Inversion Colored Inversion Colored Inversion Courtesy SMT

20 Sparse Spike Inversion Sparse Spike Inversion assumes that seismic reflectivity is a series of large spikes embedded in a background of small spikes (assumes only large spikes are meaningful). This approach seeks the simplest possible reflectivity model of spikes, that when convolved with the wavelet produces a synthetic that matches the real data. Non-linear Sparse Spike Inversion (Constrained Sparse Spike Inversion-CSSI) employs an optimization loop where the wavelet is updated (non-linearly) so that the mis-match between the synthetic and real seismic is minimized.

21 Sparse Spike Inversion Saxena and Bhatnagar, 2008

22 Model-Based Inversion Model-Based Inversion typically starts with a low frequency model of the P-impedance and then the model is perturbed until a good fit is obtained between the seismic data and a synthetic trace using the recursive formula. The Generalized Linear Inversion (GLI) method perturbs an initial acoustic impedance estimate and determines the differences between the synthetic trace and the real data (simulated annealing and least square fits are often applied). Principal Component Analysis (PCA) method computes a standard response from which the input can be generated by applying specific weighting factors determined from well control (linear interpolation). The results of the convolution with the seismic wavelet are compared with the seismic traces and the velocity and density models are perturbed to reduce the discrepancy.

23 Model-Based Inversion Veeken and Da Silva, 2004

24 Model-Based Inversion Wavelet Initial Model Seismic Data Synthetic Data Misfit Function Stochastic Model Update No Optimisation by Simulated Annealing Convergent? Yes Model Based Inversion Utilizing Simulating Annealing Stop Courtesy Equipoise

25 Model-Based Inversion Seismic section from Alberta highlighting an anomalous bright spot zone. Model based inversion of the section above Russell and Hampson, 2006

26 AVO Elastic Impedance Connolly (1999) developed an angle dependent analogy to acoustic impedance to account for AVO effects. Elastic impedance is a function of P-wave velocity, S-wave velocity, density, and incidence angle. Input for elastic impedance is some form of angle stack (e.g., near, mid, and far). Wavelets are determined for each offset angle volume. This approach is accurate for small to moderate impedance changes.

27 AVO Elastic Impedance Connolly relates the Aki-Richards (Shuey) linear Zoeppritz approximation: to an elastic impedance relationship: Combining the two expressions gives elastic impedance: Variations of this approach include using the Shuey two-term approximation (good to ), assume Vp/Vs=2 or Vs calculation, and Extended Elastic Impedance (Whitcombe et al., 2002)

28 AVO Elastic Impedance Veeken and Da Silva, 2004

29 AVO Elastic Impedance Elastic impedance displays for different offset angles Veeken and Da Silva, 2004

30 AVO Extended Elastic Impedance Whitcombe et al (2002) describes an approach that extends beyond a typical incident angle range of 0 to 30 for seismic data to -90 to +90. The sin²ɵ term in the Shuey 2-term Zoeppritz linear approximation limits the range at which reflectivities can be defined. If this approximation is re-written as R = A + Btanx then scaled by cosx then a scaled reflectivity equation can be written: R = A cosx + B sinx. Extended Elastic Impedance logs can be generated from P-wave, S- wave, and density logs for each angle x. Once x is known then equivalent seismic sections can be generated.

31 AVO Extended Elastic Impedance EEI logs at different angles: Compressional modulus = +12 Lame s parameter = =20 Shear impedance = -50 Vp/Vs = +45 Acoustic impedance = 0 Francis and Hicks, 2006

32 AVO Extended Elastic Impedance Connolly, 2010

33 AVO Extended Elastic Impedance Connolly, 2010

34 AVO Extended Elastic Impedance Francis and Hicks, 2006

35 AVO Simultaneous Inversion Simultaneous inversion is a prestack inversion method that uses multiple offset or angle stacks. This method solves for S impedance, P impedance, and density, which are key for descriminating lithology, porosity, and fluid effects. For each input partial stack a wavelet is estimated. All models, partial stacks, and wavelets are input into a single inversion algorithm and solved simultaneously compensating for offset dependent phase, bandwidth, tuning and NMO stretch effects.

36 AVO Simultaneous Inversion NMO Gather (Full Offset) Wavelets Macro-models Simultaneous Inversion Constrained Model-driven Global (SA) P-wave Impedance Model S-wave Impedance Model Simultaneous Inversion Workflow Courtesy Equipoise

37 AVO Simultaneous Inversion Russell and Hampson, 2006 Simultaneous P-Impedance Inversion

38 AVO Simultaneous Inversion Simultaneous P-Impedance Inversion Vp/Vs ratio from dividing P and S-impedance sections from simultaneous inversion Russell and Hampson, 2006

39 Predicting Other Rock Properties from Pre-Stack Inversion Pre-stack Inversion Vp/Vs Poisson s Ratio Acoustic (Ip) Impedance Shear (Is) Impedance Elastic Impedance EI Lamé Parameter Lamé Parameter Lithology Fluid Content Porosity Pore Pressure Lamé Parameter Courtesy Equipoise

40 Pre-Stack Inversion allows discriminating lithologies and fluids Bunch and Dromgoole, 1995 Pillar, 2011

41 AVO Elastic Impedance Reservoir attributes from Elastic Impedance Inversion Veeken and Da Silva, 2004

42 Elastic Impedance Simultaneous Inversion Rasmussen et al., 2004

43 Elastic Impedance-Simultaneous Inversion Poisson s Ratio From Elastic Impedance Inversion From Simultaneous Inversion Saxena and Bhatnagar, 2008

44 Stochastic Inversion - Geostatistical Geostatistical or probabilistic inversion uses quantification of uncertainties attached to the inversion input data. Probability density functions (PDF) are defined and an earth model is simulated, which is perturbed to minimize the discrepancy between the modeled and measured seismic data (simulated annealing) producing multiple realizations. The PDF determination comes from well logs, spatial properties (variograms) and lithological distributions. It is critical that the interpreter quantify the uncertainties in a realistic way, especially in the absence or minimal well control in specific areas.

45 Stochastic Inversion - Geostatistical The workflow for performing a geostatistical inversion showing the input data and the outputted model realizations. McCrank et al., 2009

46 Stochastic Inversion - Geostatistical The acoustic impedance inversion results with the mean acoustic impedance results of the multiple realizations in the upper left. McCrank et al., 2009

47 Stochastic Inversion - Bayesian Seismic inversion that uses Baye s rule allows the information from all available measurements to be integrated into a consistent image of the reservoir and constrain these solutions based on a priori knowledge about the subsurface parameters. Prior knowledge about the model parameters usually is combined with a likelihood function, which depends on the misfit between the model response and the observed seismic data

48 Stochastic Inversion - Bayesian Pillar, 2011

49 Stochastic Inversion - Bayesian Pillar, 2011

50 Categories of Seismic Inversion Running Sum Well Control Post Stack Pre Stack Relative Imp. Absolute Imp. Deterministic Probabilistic Recursive Trace Integration Colored Inversion Sparse Spike (CSSI) Model Based Inversion Elastic Impedance Simultaneous Inversion Geostatistical Inversion Bayesian Inversion

51 Applicable Inversions for Exploration through Development Running Sum Exploration Exploitation Development Recursive Trace Integration Colored Inversion Sparse Spike (CSSI) Model Based Inversion Elastic Impedance Simultaneous Inversion Geostatistical Inversion Bayesian Inversion

52 A Few Additional Types of Inversion Simultaneous geostatistical partial stack inversion Spectral inversion Pre-stack waveform inversion Multi-component inversion Multi-azimuth inversion Inversion incorporating EM and FTG data 3D Full Waveform Inversion

53 Seismic Inversion Checklist (modified from Ikon website) 1. Check log data and edit accordingly. 2. Perform rock physics analysis to determine if inversion is useful. 3. Check seismic data for proper processing and conditioning, S/N, accurate partial stacks, etc. 4. Well ties and low frequency background model: -Do you get credible wavelets? -Are wavelets consistent across wells? -Is a broadband inversion required? 5. Be sure and employ appropriate inversion algorithm. 6. QC is paramount: -Check all wells and match of wavelets. -From impedances forward model to synthetic seismic; compare with actual seismic. -Does it make geological sense? -Is low frequency trend applied correctly? 7. Use impedance results intelligently they are not the end goal!

54 SEC Regulations and Inversion Pillar, 2011

55 SEC Regulations and Inversion Oil and gas companies may use any "reliable technology" to establish reserves volumes in addition to those established by production and flow test data. Oil and gas companies may classify proved undeveloped reserves ("PUDs") any distance from known proved reserves (rather than only in immediately offsetting units) based on a reasonable certainty standard. Use of reliable technology: A registrant is required to disclose, in general terms, the technologies used in ascertaining the reasonable certainty of the PUD locations producing hydrocarbons. Any technology stated must be field tested to demonstrate consistency and repeatability. In reviewing disclosures of the technologies used, some very simple statements that a combination of these technologies was used were noted: analogy, 2-D and 3-D seismic data, volumetric and material balance analysis, decline curves, petrophysics, and log analysis. A company does not need to disclose proprietary technologies or the mix of proprietary methods.

56 Conclusions 1) Seismic Inversion is not a unique process. 2) Seismic and well log information should be optimally conditioned for reliable inversion results. 3) There is a trade-off between work involved/cost/time and quality of the final inversion. 4) A rock physics study can help determine what the end product from inversion should be and whether this is attainable with the existing data. 5) Even when advanced inversion algorithms are chosen, a simpler deterministic inversion perhaps should be run for feasibility and as a yardstick of what can be resolved. 6) Most seismic inversion workflows are not linear and require reasonable inputs and usually numerous iterations. 7) Will the SEC accept inversion results for booking reserves?

57 Indeed, given that we can now directly invert to the acoustic properties using simultaneous inversion, some workers hold the view that conventional two-term AVO techniques are now passé. This is naïve given that deriving an accurate estimate of Poisson ratio from seismic is actually quite difficult, requiring a stringent set of data conditions. In many cases, bias of one form or another is introduced into the result, for example when merging the low frequency component. The interpreter needs to understand both conventional AVO approaches as well as the latest trends in inversion. Rob Simm What makes the wiggle waggle: a perspective on rock physics, First Break, June 2011

58 DHI CONSORTIUM Thank You

Stochastic vs Deterministic Pre-stack Inversion Methods. Brian Russell

Stochastic vs Deterministic Pre-stack Inversion Methods. Brian Russell Stochastic vs Deterministic Pre-stack Inversion Methods Brian Russell Introduction Seismic reservoir analysis techniques utilize the fact that seismic amplitudes contain information about the geological

More information

Statistical Rock Physics

Statistical 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 information

A021 Petrophysical Seismic Inversion for Porosity and 4D Calibration on the Troll Field

A021 Petrophysical Seismic Inversion for Porosity and 4D Calibration on the Troll Field A021 Petrophysical Seismic Inversion for Porosity and 4D Calibration on the Troll Field T. Coleou* (CGG), A.J. van Wijngaarden (Hydro), A. Norenes Haaland (Hydro), P. Moliere (Hydro), R. Ona (Hydro) &

More information

URTeC: Summary

URTeC: Summary URTeC: 2665754 Using Seismic Inversion to Predict Geomechanical Well Behavior: a Case Study From the Permian Basin Simon S. Payne*, Ikon Science; Jeremy Meyer*, Ikon Science Copyright 2017, Unconventional

More information

Thin Sweet Spots Identification in the Duvernay Formation of North Central Alberta*

Thin Sweet Spots Identification in the Duvernay Formation of North Central Alberta* Thin Sweet Spots Identification in the Duvernay Formation of North Central Alberta* Ritesh K. Sharma 1 and Satinder Chopra 1 Search and Discovery Article #10902 (2017)** Posted January 16, 2017 *Adapted

More information

Rock Physics & Formation Evaluation. Special Topic

Rock Physics & Formation Evaluation. Special Topic Volume 30 Issue 5 May 2012 Special Topic Technical Articles Dual representation of multiscale fracture network modelling for UAE carbonate field AVO and spectral decomposition for derisking Palaeogene

More information

F003 Geomodel Update Using 4-D Petrophysical Seismic Inversion on the Troll West Field

F003 Geomodel Update Using 4-D Petrophysical Seismic Inversion on the Troll West Field F003 Geomodel Update Using 4-D Petrophysical Seismic Inversion on the Troll West Field K. Gjerding* (Statoil), N. Skjei (Statoil), A. Norenes Haaland (Statoil), I. Machecler (CGGVeritas Services) & T.

More information

Bayesian lithology/fluid inversion comparison of two algorithms

Bayesian lithology/fluid inversion comparison of two algorithms Comput Geosci () 14:357 367 DOI.07/s596-009-9155-9 REVIEW PAPER Bayesian lithology/fluid inversion comparison of two algorithms Marit Ulvmoen Hugo Hammer Received: 2 April 09 / Accepted: 17 August 09 /

More information

Geological Classification of Seismic-Inversion Data in the Doba Basin of Chad*

Geological Classification of Seismic-Inversion Data in the Doba Basin of Chad* Geological Classification of Seismic-Inversion Data in the Doba Basin of Chad* Carl Reine 1, Chris Szelewski 2, and Chaminda Sandanayake 3 Search and Discovery Article #41899 (2016)** Posted September

More information

DHI Analysis Using Seismic Frequency Attribute On Field-AN Niger Delta, Nigeria

DHI Analysis Using Seismic Frequency Attribute On Field-AN Niger Delta, Nigeria IOSR Journal of Applied Geology and Geophysics (IOSR-JAGG) e-issn: 2321 0990, p-issn: 2321 0982.Volume 1, Issue 1 (May. Jun. 2013), PP 05-10 DHI Analysis Using Seismic Frequency Attribute On Field-AN Niger

More information

Quantifying Bypassed Pay Through 4-D Post-Stack Inversion*

Quantifying Bypassed Pay Through 4-D Post-Stack Inversion* Quantifying Bypassed Pay Through 4-D Post-Stack Inversion* Robert Woock 1, Sean Boerner 2 and James Gamble 1 Search and Discovery Article #40799 (2011) Posted August 12, 2011 *Adapted from oral presentation

More information

Geophysical and geomechanical rock property templates for source rocks Malleswar Yenugu, Ikon Science Americas, USA

Geophysical and geomechanical rock property templates for source rocks Malleswar Yenugu, Ikon Science Americas, USA Geophysical and geomechanical rock property templates for source rocks Malleswar Yenugu, Ikon Science Americas, USA Summary Sweet spot identification for source rocks involve detection of organic rich,

More information

Neural Networks Find Meaning In Data

Neural Networks Find Meaning In Data t Neural Networks Find Meaning n Data 1,_ By Brian Russell CALGARYThe ultimate goal of any seismic analysis is to extend reservoir knowledge beyond the wells. Neural networks represent a class of statistical

More information

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

Downloaded 11/22/17 to Redistribution subject to SEG license or copyright; see Terms of Use at t Special section: Thin beds Downloaded 11/22/17 to 37.44.193.54. Redistribution subject to SEG license or copyright; see Terms of Use at http://library.seg.org/ Tuning of flat spots with overlying bright

More information

NORGES TEKNISK-NATURVITENSKAPELIGE UNIVERSITET. Directional Metropolis Hastings updates for posteriors with nonlinear likelihoods

NORGES TEKNISK-NATURVITENSKAPELIGE UNIVERSITET. Directional Metropolis Hastings updates for posteriors with nonlinear likelihoods NORGES TEKNISK-NATURVITENSKAPELIGE UNIVERSITET Directional Metropolis Hastings updates for posteriors with nonlinear likelihoods by Håkon Tjelmeland and Jo Eidsvik PREPRINT STATISTICS NO. 5/2004 NORWEGIAN

More information

V003 How Reliable Is Statistical Wavelet Estimation?

V003 How Reliable Is Statistical Wavelet Estimation? V003 How Reliable Is Statistical Wavelet Estimation? J.A. Edgar* (BG Group) & M. van der Baan (University of Alberta) SUMMARY Well logs are often used for the estimation of seismic wavelets. The phase

More information

Nonlinear Bayesian joint inversion of seismic reflection

Nonlinear Bayesian joint inversion of seismic reflection Geophys. J. Int. (2) 142, Nonlinear Bayesian joint inversion of seismic reflection coefficients Tor Erik Rabben 1,, Håkon Tjelmeland 2, and Bjørn Ursin 1 1 Department of Petroleum Engineering and Applied

More information

EOS 350 MIDTERM OCT 4, 2013 STUDENT NAME: TEAM #:

EOS 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 information

Special section: Seismic attributes. Relative acoustic impedance from wavelet transform

Special section: Seismic attributes. Relative acoustic impedance from wavelet transform t Special section: Seismic attributes Relative acoustic impedance from wavelet transform Marcílio Castro de Matos 1, Rodrigo Penna 2, Paulo Johann 2, and Kurt Marfurt 3 Abstract Most deconvolution algorithms

More information

Seismic 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 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 information

A simple algorithm for band-limited impedance inversion

A simple algorithm for band-limited impedance inversion Impedance inversion algorithm A simple algorithm for band-limited impedance inversion Robert J. Ferguson and Gary F. Margrave ABSTRACT This paper describes a seismic inversion method which has been cast

More information

Research - Recent work:

Research - Recent work: Research - Recent work: Henning Omre - www.math.ntnu.no/~omre Professor in Statistics Department of Mathematical Sciences Norwegian University of Science & Technology Trondheim Norway December 2016 Current

More information

Seismic Attributes and Their Applications in Seismic Geomorphology

Seismic Attributes and Their Applications in Seismic Geomorphology Academic article Seismic Attributes and Their Applications in Seismic Geomorphology Sanhasuk Koson, Piyaphong Chenrai* and Montri Choowong Department of Geology, Faculty of Science, Chulalongkorn University,

More information

Relinquishment Report. Licence P2016 Block 205/4c

Relinquishment Report. Licence P2016 Block 205/4c Relinquishment Report Licence P2016 Block 205/4c September 2016 JX Nippon Exploration and Production (U.K.) Limited Idemitsu Petroleum UK Ltd. Table of Contents 1. Licence Information... 1 2. Licence Synopsis...

More information

Parameter Estimation and Sensitivity Analysis in Clastic Sedimentation Modeling

Parameter Estimation and Sensitivity Analysis in Clastic Sedimentation Modeling Parameter Estimation and Sensitivity Analysis in Clastic Sedimentation Modeling A. Acevedo 1, A. Khramtsov 2, H. A. Madhoo 3, L. Noomee 4, and D. Tetzlaff 5 1 Schlumberger Information Solutions,Gatwick,

More information

F-K Characteristics of the Seismic Response to a Set of Parallel Discrete Fractures

F-K Characteristics of the Seismic Response to a Set of Parallel Discrete Fractures F-K Characteristics of the Seismic Response to a Set of Parallel Discrete Fractures Yang Zhang 1, Xander Campman 1, Samantha Grandi 1, Shihong Chi 1, M. Nafi Toksöz 1, Mark E. Willis 1, Daniel R. Burns

More information

Main Challenges and Uncertainties for Oil Production from Turbidite Reservoirs in Deep Water Campos Basin, Brazil*

Main Challenges and Uncertainties for Oil Production from Turbidite Reservoirs in Deep Water Campos Basin, Brazil* Main Challenges and Uncertainties for Oil Production from Turbidite Reservoirs in Deep Water Campos Basin, Brazil* Carlos H. Bruhn 1, Antonio Pinto 1, and Paulo R. Johann 1 Search and Discovery Article

More information

ROCK PHYSICS MODELING FOR LITHOLOGY PREDICTION USING HERTZ- MINDLIN THEORY

ROCK PHYSICS MODELING FOR LITHOLOGY PREDICTION USING HERTZ- MINDLIN THEORY ROCK PHYSICS MODELING FOR LITHOLOGY PREDICTION USING HERTZ- MINDLIN THEORY Ida Ayu PURNAMASARI*, Hilfan KHAIRY, Abdelazis Lotfy ABDELDAYEM Geoscience and Petroleum Engineering Department Universiti Teknologi

More information

P-wave reflection coefficients for transversely isotropic models with vertical and horizontal axis of symmetry

P-wave reflection coefficients for transversely isotropic models with vertical and horizontal axis of symmetry GEOPHYSICS, VOL. 62, NO. 3 (MAY-JUNE 1997); P. 713 722, 7 FIGS., 1 TABLE. P-wave reflection coefficients for transversely isotropic models with vertical and horizontal axis of symmetry Andreas Rüger ABSTRACT

More information

Interpretation of PP and PS seismic data from the Mackenzie Delta, N.W.T.

Interpretation of PP and PS seismic data from the Mackenzie Delta, N.W.T. Interpretation of PP and PS seismic data from the Mackenzie Delta, N.W.T. Carlos E. Nieto* CREWES. The University of Calgary nieto@geo.ucalgary.ca and Robert R. Stewart CREWES. The University of Calgary

More information

A Modular NMF Matching Algorithm for Radiation Spectra

A Modular NMF Matching Algorithm for Radiation Spectra A Modular NMF Matching Algorithm for Radiation Spectra Melissa L. Koudelka Sensor Exploitation Applications Sandia National Laboratories mlkoude@sandia.gov Daniel J. Dorsey Systems Technologies Sandia

More information

Chapter 7: Reflection Seismology Homework Solutions (Jan. 2010)

Chapter 7: Reflection Seismology Homework Solutions (Jan. 2010) Chapter 7: eflection Seismology Homework Solutions (Jan. 200). Why do marine seismic reflection surveys not record (a) S waves? (b) refracted rays? 2 μ a) For ideal fluid, μ=0, thus, v s = = 0 ρ b) eflection

More information

A new method for multi-exponential inversion of NMR relaxation measurements

A new method for multi-exponential inversion of NMR relaxation measurements Science in China Ser. G Physics, Mechanics & Astronomy 2004 Vol.47 No.3 265 276 265 A new method for multi-exponential inversion of NMR relaxation measurements WANG Zhongdong 1, 2, XIAO Lizhi 1 & LIU Tangyan

More information

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

Crosswell tomography imaging of the permeability structure within a sandstone oil field. Crosswell tomography imaging of the permeability structure within a sandstone oil field. Tokuo Yamamoto (1), and Junichi Sakakibara (2) (1) University of Miami and Yamamoto Engineering Corporation, (2)

More information

Th Guided Waves - Inversion and Attenuation

Th Guided Waves - Inversion and Attenuation Th-01-08 Guided Waves - Inversion and Attenuation D. Boiero* (WesternGeco), C. Strobbia (WesternGeco), L. Velasco (WesternGeco) & P. Vermeer (WesternGeco) SUMMARY Guided waves contain significant information

More information

Acoustic Anisotropy Measurements and Interpretation in Deviated Wells

Acoustic Anisotropy Measurements and Interpretation in Deviated Wells Acoustic Anisotropy Measurements and Interpretation in Deviated Wells X. M. Tang, and D. Patterson, Houston Technology Center, Baker Atlas, Houston, Texas, USA ABSTRACT Many acoustic anisotropy measurements

More information

APPLICATION OF OPTIMAL BASIS FUNCTIONS IN FULL WAVEFORM INVERSION

APPLICATION OF OPTIMAL BASIS FUNCTIONS IN FULL WAVEFORM INVERSION METHODS AND APPLICATIONS OF ANALYSIS. c 2004 International Press Vol. 11, No. 3, pp. 345 352, September 2004 004 APPLICATION OF OPTIMAL BASIS FUNCTIONS IN FULL WAVEFORM INVERSION PING SHENG, GANG SUN,

More information

IPTC Abstract. Introduction

IPTC Abstract. Introduction IPTC 14817 Application of PSDM Imaging for Reservoir Characterisation in the Northern Malay Basin: A Case Study Loic Michel and Chee Hau Hoo, Hampson Russell, and Mark.J. Flynn, Carigali Hess Sdn. Bhd.

More information

Pressure Prediction and Hazard Avoidance through Improved Seismic Imaging

Pressure Prediction and Hazard Avoidance through Improved Seismic Imaging Pressure Prediction and Hazard Avoidance through Improved Seismic Imaging 12121-6002-02 Bill Barkhouse SEG Advanced Modeling Corporation (SEAM) RPSEA / SPE-GCS Ultra-Deepwater Technology Conference September

More information

Impact of Rock Physics Depth Trends and Markov Random. Fields on Hierarchical Bayesian Lithology/Fluid Prediction.

Impact of Rock Physics Depth Trends and Markov Random. Fields on Hierarchical Bayesian Lithology/Fluid Prediction. Impact of Rock Physics Depth Trends and Markov Random Fields on Hierarchical Bayesian Lithology/Fluid Prediction. Kjartan Rimstad 1 & Henning Omre 1 1 Norwegian University of Science and Technology, Trondheim,

More information

Chapter 6. Conclusions. 6.1 Conclusions and perspectives

Chapter 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 information

Well Evaluation. An Integrated Well-Evaluation Process for Shale Gas Reservoirs

Well Evaluation. An Integrated Well-Evaluation Process for Shale Gas Reservoirs Well Evaluation An Integrated Well-Evaluation Process for Shale Gas Reservoirs APPLICATION Gas shale analysis for unconventional wells BENEFITS Optimized stimulation by thorough study and planning Increased

More information

TUTORIAL WORKFLOWS. What s the added value to my interpretation project?

TUTORIAL WORKFLOWS. What s the added value to my interpretation project? Rock Physics Models TUTORIAL WORKFLOWS This section contains: Rock Physics Models Tutorial workflow names What is a workflow? Symbols used in rock physics formulas Tutorial workflows overview Starting

More information

GRAIN SORTING, POROSITY, AND ELASTICITY. Jack Dvorkin and Mario A. Gutierrez Geophysics Department, Stanford University ABSTRACT

GRAIN SORTING, POROSITY, AND ELASTICITY. Jack Dvorkin and Mario A. Gutierrez Geophysics Department, Stanford University ABSTRACT GRAIN SORTING, POROSITY, AND ELASTICITY Jack Dvorkin and Mario A. Gutierrez Geophysics Department, Stanford University July 24, 2001 ABSTRACT Grain size distribution (sorting) is determined by deposition.

More information

Geologic and Reservoir Characterization and Modeling

Geologic and Reservoir Characterization and Modeling Geologic and Reservoir Characterization and Modeling Scott M. Frailey and James Damico Illinois State Geological Survey Midwest Geologic Sequestration Science Conference November 8 th, 2011 Acknowledgements

More information

Quantitative Use of Seismic Attributes for Reservoir Characterization Introduction Primitive Hybrid Taner

Quantitative Use of Seismic Attributes for Reservoir Characterization Introduction Primitive Hybrid Taner Quantitative Use of Seismic Attributes for Reservoir Characterization Richard L. Chambers* and Jeffrey M. Yarus** Quantitative Geosciences, Inc. Broken Arrow, Ok, U.S.A* Houston, Texas, U.S.A** Introduction

More information

3D Seismic Reservoir Characterization and Delineation in Carbonate Reservoir*

3D Seismic Reservoir Characterization and Delineation in Carbonate Reservoir* 3D Seismic Reservoir Characterization and Delineation in Carbonate Reservoir* M. N. Alamsyah 1, Bambang W. Handono 1, and Andri Syafriya 1 Search and Discovery Article #41760 (2016) Posted January 25,

More information

AN INTEGRATED STUDY OF THE GRAYBURG/SAN ANDRES RESERVOIR, FOSTER AND SOUTH COWDEN FIELDS, ECTOR COUNTY TEXAS. Annual Technical Report

AN INTEGRATED STUDY OF THE GRAYBURG/SAN ANDRES RESERVOIR, FOSTER AND SOUTH COWDEN FIELDS, ECTOR COUNTY TEXAS. Annual Technical Report AN INTEGRATED STUDY OF THE GRAYBURG/SAN ANDRES RESERVOIR, FOSTER AND SOUTH COWDEN FIELDS, ECTOR COUNTY TEXAS Annual Technical Report Reporting Period: 08/08/1996 to 07/31/1997 Authors: Robert C. Trentham,

More information

Nuclear Magnetic Resonance Log

Nuclear Magnetic Resonance Log Objective The development of the nuclear magnetic resonance (NMR) log was fueled by the desire to obtain an estimate of permeability from a continuous measurement. Previous work had relied on empirical

More information

Available online at ScienceDirect. Energy Procedia 114 (2017 )

Available online at  ScienceDirect. Energy Procedia 114 (2017 ) Available online at www.sciencedirect.com ScienceDirect Energy Procedia 114 (2017 ) 2772 2780 13th International Conference on Greenhouse Gas Control Technologies, GHGT-13, 14-18 November 2016, Lausanne,

More information

Copyright. Barry Michael Borgman

Copyright. Barry Michael Borgman Copyright by Barry Michael Borgman 2016 The Thesis Committee for Barry Michael Borgman Certifies that this is the approved version of the following thesis: Characterization of the Cana-Woodford Shale using

More information

OGIP Evaluation of Shale Gas and CBM with Basin Modeling and Gas Isotopes Interpretation*

OGIP Evaluation of Shale Gas and CBM with Basin Modeling and Gas Isotopes Interpretation* OGIP Evaluation of Shale Gas and CBM with Basin Modeling and Gas Isotopes Interpretation* Daniel Xia 1,2 and Yongchun Tang 1 Search and Discovery Article #40842 (2011) Posted November 30, 2011 *Adapted

More information

10-810: Advanced Algorithms and Models for Computational Biology. Optimal leaf ordering and classification

10-810: Advanced Algorithms and Models for Computational Biology. Optimal leaf ordering and classification 10-810: Advanced Algorithms and Models for Computational Biology Optimal leaf ordering and classification Hierarchical clustering As we mentioned, its one of the most popular methods for clustering gene

More information

Hague and London Oil Plc

Hague and London Oil Plc Hague and London Oil Plc Exploiting Shallow Gas Offshore Netherlands De-Gassing Shallow Hazards or De-Risking Shallow Gas November 2016 The Shallow Gas Play can be Summarized as Schematic evolution of

More information

Integrating Geomechanics and Reservoir Characterization Examples from Canadian Shale Plays

Integrating Geomechanics and Reservoir Characterization Examples from Canadian Shale Plays Integrating Geomechanics and Reservoir Characterization Examples from Canadian Shale Plays AAPG Geosciences Technology Workshops Geomechanics and Reservoir Characterization of Shale and Carbonates July

More information

Frequently Asked Questions

Frequently Asked Questions Frequently Asked Questions 1 Is THERMOVISION TOMOGRAPHY technology independent, or it requires previous seismic or other data? THERMOVISION TOMOGRAPHY technology is a stand-alone technology. Independent

More information

IDENTIFICATION FOR CONTROL

IDENTIFICATION FOR CONTROL IDENTIFICATION FOR CONTROL Raymond A. de Callafon, University of California San Diego, USA Paul M.J. Van den Hof, Delft University of Technology, the Netherlands Keywords: Controller, Closed loop model,

More information

Determination of the signature of a dynamite source using source scaling, Part 2: Experiment

Determination of the signature of a dynamite source using source scaling, Part 2: Experiment GEOPHYSICS, VOL. 58, NO. 8 (AUGUST 1993); P. 1183-1194, 12 FIGS. Determination of the signature of a dynamite source using source scaling, Part 2: Experiment Anton Ziolkowski* and Karel Bokhorst ABSTRACT

More information

Determination of Incompressibility, Elasticity and the Rigidity of Surface Soils and Shallow Sediments from Seismic Wave Velocities

Determination of Incompressibility, Elasticity and the Rigidity of Surface Soils and Shallow Sediments from Seismic Wave Velocities Journal of Earth Sciences and Geotechnical Engineering, vol. 6, no.1, 2016, 99-111 ISSN: 1792-9040 (print), 1792-9660 (online) Scienpress Ltd, 2016 Determination of Incompressibility, Elasticity and the

More information

Delivery: an open source model based Bayesian seismic inversion program

Delivery: an open source model based Bayesian seismic inversion program Delivery: an open source model based Bayesian seismic inversion program James Gunning a,, Michael E. Glinsky b a Australian Petroleum Cooperative Research Centre, CSIRO Division of Petroleum Resources,

More information

LAMINATION THEORY FOR THE STRENGTH OF FIBER COMPOSITE MATERIALS

LAMINATION THEORY FOR THE STRENGTH OF FIBER COMPOSITE MATERIALS XXII. LAMINATION THEORY FOR THE STRENGTH OF FIBER COMPOSITE MATERIALS Introduction The lamination theory for the elastic stiffness of fiber composite materials is the backbone of the entire field, it holds

More information

Modified Kriging: Evaluation, Modification, Recommendations

Modified Kriging: Evaluation, Modification, Recommendations UCRL-ID- 12873 1 Modified Kriging: Evaluation, Modification, Recommendations and C.A. Schultz S.C. Myers September 1997 This is an informal report intended primarily forinternal or limited external distribution.

More information

The effect of interpolation on imaging and AVO: A Viking case study

The effect of interpolation on imaging and AVO: A Viking case study GEOPHYSICS, VOL. 75, NO. 6 NOVEMBER-DECEMBER 21 ; P. WB265 WB274, 11 FIGS., 1 TABLE. 1.119/1.347539 The effect of interpolation on imaging and : A Viking case study Lee Hunt 1, Jon Downton 2, Scott Reynolds

More information

Study on Rayleigh Wave Inversion for Estimating Shear-wave Velocity Profile

Study on Rayleigh Wave Inversion for Estimating Shear-wave Velocity Profile PROC. ITB Eng. Science Vol. 35 B, No. 1, 2003, 11-22 11 Study on Rayleigh Wave Inversion for Estimating Shear-wave Velocity Profile T.A. Sanny*, Yoes Avianto* *Geophysical Engineering Department, Faculty

More information

with wide-angle prestack migration image crustal faults using only P, P, isotropic scattering characteristics and

with wide-angle prestack migration image crustal faults using only P, P, isotropic scattering characteristics and 34 Chapter 3 Isotropic scattering and crustal imaging with wide-angle prestack migration 3.1 Introduction Earthquake seismologists may be reluctant to believe that one can successfully image crustal faults

More information

Importance of regional geological and petrophysical analysis: Bakken/Three Forks formations of the Williston basin North Dakota and Montana

Importance of regional geological and petrophysical analysis: Bakken/Three Forks formations of the Williston basin North Dakota and Montana Importance of regional geological and petrophysical analysis: Bakken/Three Forks formations of the Williston basin North Dakota and Montana Presenter: Stefani Brakenhoff The Discovery Group Team: Bob Cluff,

More information

FULL MOMENT TENSOR ANALYSIS USING FIRST MOTION DATA AT THE GEYSERS GEOTHERMAL FIELD

FULL MOMENT TENSOR ANALYSIS USING FIRST MOTION DATA AT THE GEYSERS GEOTHERMAL FIELD PROCEEDINGS, Thirty-Eighth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 11-13, 2013 SGP-TR-198 FULL MOMENT TENSOR ANALYSIS USING FIRST MOTION DATA AT

More information

Microseismic Monitoring: Insights from Moment Tensor Inversion

Microseismic Monitoring: Insights from Moment Tensor Inversion Microseismic Monitoring: Insights from Moment Tensor Inversion Farshid Forouhideh and David W. Eaton ABSTRACT This paper reviews the mathematical tools used for describing microseismic source mechanisms.

More information

Enhanced Formation Evaluation of Shales Using NMR Secular Relaxation*

Enhanced Formation Evaluation of Shales Using NMR Secular Relaxation* Enhanced Formation Evaluation of Shales Using NMR Secular Relaxation* Hugh Daigle 1, Andrew Johnson 1, Jameson P. Gips 1, and Mukul Sharma 1 Search and Discovery Article #41411 (2014)** Posted August 11,

More information

Search and Discovery Article #51409 (2017)** Posted August 7, Abstract. Selected References

Search and Discovery Article #51409 (2017)** Posted August 7, Abstract. Selected References Saturations of Migrating Buoyant Fluids from Invasion Percolation Flow Simulation Using Small-Scale, High- Resolution Geologic Models With Realistic Heterogeneity* Timothy A. Meckel 1, Luca Trevisan 2,

More information

Determination of Duvernay Formation Reservoir Properties through Probabilistic Petrophysical Analysis calibrated to Core Studies.

Determination of Duvernay Formation Reservoir Properties through Probabilistic Petrophysical Analysis calibrated to Core Studies. Determination of Duvernay Formation Reservoir Properties through Probabilistic Petrophysical Analysis calibrated to Core Studies. Nasir Rahim, Neil Watson Canadian Discovery Ltd. Summary The petrophysical

More information

Search and Discovery Article #40507 (2010) Posted February 12, Abstract

Search and Discovery Article #40507 (2010) Posted February 12, Abstract Click to view slides prepared for presentation. Application of From Seismic Interpretation to Tectonic Reconstruction Methodology to Study Pre-Tertiary Fractured Granite Basement Reservoir in Cuu Long

More information

and a contribution from Offshore Europe

and a contribution from Offshore Europe Primary funding is provided by The SPE Foundation through member donations and a contribution from Offshore Europe The Society is grateful to those companies that allow their professionals to serve as

More information

Exploration Well Failures from the Moray Firth & Central North Sea (UK) 21 st Century Exploration Road Map Project Christian Mathieu

Exploration Well Failures from the Moray Firth & Central North Sea (UK) 21 st Century Exploration Road Map Project Christian Mathieu Exploration Well Failures from the Moray Firth & Central North Sea (UK) 21 st Century Exploration Road Map Project Christian Mathieu 1. Project Objectives, Time-Line and Status 2. A few statistics 3. Reasons

More information

Relinquishment Report for Licence Number P1471 Block 16/8f March 2009

Relinquishment Report for Licence Number P1471 Block 16/8f March 2009 Relinquishment Report for Licence Number P1471 Block 16/8f March 2009 Endeavour Energy UK Ltd Relinquishment Report for Licence Number P1471 Block 16/8f 1 Table of Contents 1.0 Header...3 2.0 Synopsis...3

More information

Heterogeneity Type Porosity. Connected Conductive Spot. Fracture Connected. Conductive Spot. Isolated Conductive Spot. Matrix.

Heterogeneity Type Porosity. Connected Conductive Spot. Fracture Connected. Conductive Spot. Isolated Conductive Spot. Matrix. Porosity Histogram Porosity Contribution 1.3.3.3 Connected 9.8 ohm.m U R D 9 18 7. 5.25 4.38 3.5 2.63 1.75 48 Heterogeneity Distribution Image Orientation, L U 27 36.4.3 X,X72.5 Depth, ft.3 1 Isolated.3

More information

Receiver Function Inversion

Receiver Function Inversion Receiver Function Inversion Advanced Studies Institute on Seismological Research Kuwait City, Kuwait - January 19-22, 2013 Jordi Julià Universidade Federal do Rio Grande do Norte, Brasil Outline Introduction

More information

Application of the Ensemble Kalman Filter to History Matching

Application of the Ensemble Kalman Filter to History Matching Application of the Ensemble Kalman Filter to History Matching Presented at Texas A&M, November 16,2010 Outline Philosophy EnKF for Data Assimilation Field History Match Using EnKF with Covariance Localization

More information

IMPEDANCE MICROSTRUCTURE OF KEROGEN SHALES

IMPEDANCE MICROSTRUCTURE OF KEROGEN SHALES IMPEDANCE MICROSTRUCTURE OF KEROGEN SHALES Manika Prasad Rock Abuse Laboratory Colorado School of Mines Presented at the 27th Oil Shale Symposium at CSM on October 17, 2007 OUTLINE Motivation Working Principles

More information

Borehole Geophysics. Acoustic logging measurements

Borehole Geophysics. Acoustic logging measurements Acoustic logging measurements - Review of basic physics background - Concept of P- and S-wave measurements and logging tools - Tube waves - Seismic imaging - Synthetic seismograms - Field application examples

More information

ASSESSMENT OF LIQUEFACTION RISK IN THE HARBOUR AREA OF GIOIA TAURO

ASSESSMENT OF LIQUEFACTION RISK IN THE HARBOUR AREA OF GIOIA TAURO ASSESSMENT OF LIQUEFACTION RISK IN THE HARBOUR AREA OF GIOIA TAURO International Doctoral Course on Risks management on the built environment Florence, 26 nov 2002 Johann Facciorusso In a strict sense

More information

EXTENDED ABSTRACT EVALUATING THE SHALY SAND OIL RESERVOIRS OF EL TORDILLO FIELD, ARGENTINA, USING MAGNETIC RESONANCE LOGS

EXTENDED ABSTRACT EVALUATING THE SHALY SAND OIL RESERVOIRS OF EL TORDILLO FIELD, ARGENTINA, USING MAGNETIC RESONANCE LOGS EXTENDED ABSTRACT EVALUATING THE SHALY SAND OIL RESERVOIRS OF EL TORDILLO FIELD, ARGENTINA, USING MAGNETIC RESONANCE LOGS Maged Fam, Halliburton Energy Services, Luis P. Stinco, and Julio. A. Vieiro, Tecpetrol

More information

Body-wave radiation patterns and AVO in transversely isotropic media

Body-wave radiation patterns and AVO in transversely isotropic media GEOPHYSICS, VOL. 60, NO. 5 (SEPTEMBER-OCTOBER 1995); P. 1409-1425, 11 FIGS. Body-wave radiation patterns and AVO in transversely isotropic media llya Tsvankin* ABSTRACT The angular dependence of reflection

More information

Bayesian Modeling and Classification of Neural Signals

Bayesian Modeling and Classification of Neural Signals Bayesian Modeling and Classification of Neural Signals Michael S. Lewicki Computation and Neural Systems Program California Institute of Technology 216-76 Pasadena, CA 91125 lewickiocns.caltech.edu Abstract

More information

IBIS - IN FLIGHT CALIBRATION ACTIVITIES

IBIS - IN FLIGHT CALIBRATION ACTIVITIES IBIS - IN FLIGHT CALIBRATION ACTIVITIES Natalucci on behalf of the IBIS Team INTEGRAL/IBIS Calibration Update THE IBIS INSTRUMENT coded mask, placed 3.2m above detector (1m2) Collection area ~ 3000 cm2

More information

Numerical Modeling for Different Types of Fractures

Numerical Modeling for Different Types of Fractures umerical Modeling for Different Types of Fractures Xiaoqin Cui* CREWES Department of Geoscience University of Calgary Canada xicui@ucalgary.ca and Laurence R. Lines Edward S. Krebes Department of Geoscience

More information

Relinquishment report P.1190 & P Blocks: 204/13, 204/14b

Relinquishment report P.1190 & P Blocks: 204/13, 204/14b Relinquishment report P.1190 & P.1262 Blocks: 204/13, 204/14b CONTENTS 1: LICENCE DETAILS... 2 2: LICENCE SYNOPSIS... 2 3: WORK PROGRAMME SUMMARY... 3 4. DATABASE... 5 5: PROSPECTIVITY UPDATE... 6 6. RESERVE

More information

Bayesian Inference in Astronomy & Astrophysics A Short Course

Bayesian Inference in Astronomy & Astrophysics A Short Course Bayesian Inference in Astronomy & Astrophysics A Short Course Tom Loredo Dept. of Astronomy, Cornell University p.1/37 Five Lectures Overview of Bayesian Inference From Gaussians to Periodograms Learning

More information

ECE 3800 Probabilistic Methods of Signal and System Analysis

ECE 3800 Probabilistic Methods of Signal and System Analysis ECE 3800 Probabilistic Methods of Signal and System Analysis Dr. Bradley J. Bazuin Western Michigan University College of Engineering and Applied Sciences Department of Electrical and Computer Engineering

More information

P11 A Methodology to Rescale the Vertical Coordinate into Geological Time Using a Geostatistical Analysis of Log Data

P11 A Methodology to Rescale the Vertical Coordinate into Geological Time Using a Geostatistical Analysis of Log Data P A Methodology to Rescale the Vertical Coordinate into Geological Time Using a Geostatistical Analysis of Log Data M. Lefranc* (Andra / EMP), B. Beaudoin (Ecole des Mines de Paris), J.P. Chilès (Ecole

More information

This paper was prepared for presentation at the Unconventional Resources Technology Conference held in San Antonio, Texas, USA, 1-3 August 2016.

This paper was prepared for presentation at the Unconventional Resources Technology Conference held in San Antonio, Texas, USA, 1-3 August 2016. URTeC: 2444366 Using Depletion-Zone Microseismicity to Understand Producing Volumes Jonathan P. McKenna*, Michael H. Grealy, Michael S. Blaz and Nathan M. Toohey, MicroSeismic, Inc. Copyright 2016, Unconventional

More information

25th Seismic Research Review - Nuclear Explosion Monitoring: Building the Knowledge Base

25th Seismic Research Review - Nuclear Explosion Monitoring: Building the Knowledge Base Pn Q UNDER TIBET AND TIENSHAN WITH PRACTICAL AND SCIENTIFIC IMPLICATIONS Jiakang Xie Lamont-Doherty Earth Observatory, Columbia University Sponsored by Defense Threat Reduction Agency Contract No. DTRA01-00-C-0048

More information

Anisotropic velocity changes in seismic time-lapse data Jorg Herwanger*, Ed Palmer (WesternGeco), Christian Rau Schiøtt (Hess)

Anisotropic velocity changes in seismic time-lapse data Jorg Herwanger*, Ed Palmer (WesternGeco), Christian Rau Schiøtt (Hess) Anisotropic velocity changes in seismic time-lapse data Jorg Herwanger*, Ed Palmer (WesternGeco), Christian Rau Schiøtt (Hess) Summary Reservoir production causes subsurface deformations and changes in

More information

Theory - SBRC fundamentals

Theory - SBRC fundamentals Chapter 2 Theory - SBRC fundamentals The approach proposed to estimate hydraulic properties of rocks using microseismicity is called Seismicity Based Reservoir Characterization (SBRC). It uses a spatio-temporal

More information

Bayesian Paradigm. Maximum A Posteriori Estimation

Bayesian Paradigm. Maximum A Posteriori Estimation Bayesian Paradigm Maximum A Posteriori Estimation Simple acquisition model noise + degradation Constraint minimization or Equivalent formulation Constraint minimization Lagrangian (unconstraint minimization)

More information

SEAPEX 2017 Central Luconia Exploration Update: Are we good or did we just get lucky? James Clark Paul Owen Sean O Brien Barry Dawe

SEAPEX 2017 Central Luconia Exploration Update: Are we good or did we just get lucky? James Clark Paul Owen Sean O Brien Barry Dawe SEAPEX 2017 Central Luconia Exploration Update: Are we good or did we just get lucky? James Clark Paul Owen Sean O Brien Barry Dawe Presentation Outline Introduction Exploration Recap Phase 2 results Impact

More information

B. Milkereit 1, E. Adam 2, Z. Li 3, Wei Qian 1, T. Bohlen 4, D. Banerjee 1, and D.R. Schmitt 5

B. Milkereit 1, E. Adam 2, Z. Li 3, Wei Qian 1, T. Bohlen 4, D. Banerjee 1, and D.R. Schmitt 5 Multi-offset vertical seismic profiling: an experiment to assess petrophysical-scale parameters at the JAPEX/JNOC/GSC et al. Mallik 5L-38 gas hydrate production research well B. Milkereit 1, E. Adam 2,

More information

Simulated and Observed Scaling in Earthquakes Kasey Schultz Physics 219B Final Project December 6, 2013

Simulated and Observed Scaling in Earthquakes Kasey Schultz Physics 219B Final Project December 6, 2013 Simulated and Observed Scaling in Earthquakes Kasey Schultz Physics 219B Final Project December 6, 2013 Abstract Earthquakes do not fit into the class of models we discussed in Physics 219B. Earthquakes

More information

L11: Pattern recognition principles

L11: Pattern recognition principles L11: Pattern recognition principles Bayesian decision theory Statistical classifiers Dimensionality reduction Clustering This lecture is partly based on [Huang, Acero and Hon, 2001, ch. 4] Introduction

More information

ADI iterations for. general elliptic problems. John Strain Mathematics Department UC Berkeley July 2013

ADI iterations for. general elliptic problems. John Strain Mathematics Department UC Berkeley July 2013 ADI iterations for general elliptic problems John Strain Mathematics Department UC Berkeley July 2013 1 OVERVIEW Classical alternating direction implicit (ADI) iteration Essentially optimal in simple domains

More information