Stochastic vs Deterministic Prestack Inversion Methods. Brian Russell

 Aubrey Davidson
 8 months ago
 Views:
Transcription
1 Stochastic vs Deterministic Prestack Inversion Methods Brian Russell
2 Introduction Seismic reservoir analysis techniques utilize the fact that seismic amplitudes contain information about the geological properties of the reservoir. The mathematics behind this observation was developed in the early 1900s, but its application to exploration seismic data did not start until the 1970s. We classify these methods into two categories: methods that analyze only the amplitudes, and methods that invert the amplitudes to reservoir properties. Newer methods analyze prestack data, where the analysis of the amplitudes without inversion is called Amplitude versus Offset, or AVO. Prestack inversion has many forms, where the major division is between deterministic and stochastic, or geostatistical, methods. In this talk I will discuss these methods and look at their assumptions and limitations.
3 A suggested workflow Well Log Data Seismic Data V P, ρ V P, V S, ρ Modeling for V S Build rock physics model Poststack only Poststack inversion Gathers, only offsets AVO & Prestack inversion Gathers with azimuths AVAz / Fracture Identification Integrate using multivariate or Bayesian statistics
4 Seismic Inversion Methods Inversion methods Poststack Prestack Model Based Recursive Sparse spike Colored Elastic Impedance LMR Simultaneous Joint PP/PS Inversion 4D Inversion Azimuthal Inversion Stochastic / Geostatistical Inversion
5 The basic model for inversion The zero offset, or stacked, seismic trace can be modeled as the convolution of the acoustic impedance (AI) reflectivity with the wavelet. As shown in the next slide, this is the basis for poststack inversion. Acoustic Impedance AI = ρvp Reflectivity AI R AI = 2AI Wavelet W Seismic W S = W * R AI
6 Poststack inversion Poststack seismic Inversion, developed in the 1970s, reverses the forward modeling procedure, allowing us to derive the impedance from the reflectivity: Inverse Wavelet 6 Impedance Reflectivity Seismic
7 Qualitative AVO In the 1980s, geophysicists observed that the amplitudes in a seismic gather could be written in linearized form using the amplitude versus offset (AVO) equation, a reformulation of the AkiRichards linearized solution to the Zoeppritz equations: Note that this has added two extra terms to the zerooffset case, a gradient term G and a curvature term C, often referred to as A, B and C, where the term A is called the intercept. This formed the basis to what I refer to as qualitative AVO. where :, sin tan sin ) ( θ θ θ θ C G R R AI P + + =. 2 and, 2 4 2, P P P S S S P S P P P P AI V V C V V V V V V V V G V V R = = + = ρ ρ ρ ρ
8 Intercept and gradient analysis The amplitudes are extracted at all times, two of which are shown here: Offset or Angle θ +R AI +G sin 2 θ Time The AVO equation predicts a linear relationship between these amplitudes and sin 2 θ. Regression curves are calculated to give R AI and G values for each time sample. R AI  G
9 Using the angle gathers for inversion Fatti et al. (1994) reformulated this equation to show that the prestack seismic data is a function of the acoustic impedance reflectivity (R AI ), shear impedance reflectivity (R SI ) and density reflectivity (R D ) term: R P( θ ) = arai + brsi + crd, where R AI AI VS ρ SI =, RSI = + =, SI = ρvs, RD 2AI 2V 2ρ 2SI S = ρ, 2ρ a = 1+ tan 2 θ, b V = 8 V S P 2 sin 2 θ, and c = V 4 V S P 2 sin 2 θ tan 2 θ.
10 Independent prestack inversion Independent prestack inversion is implemented by first extracting the reflectivity components, and then inverting them separately. To estimate the reflectivities, the amplitudes at each time t in an N trace angle gather are picked as shown here, to give R P (θ 1 ) R P (θ N ): We can then solve for the reflectivities at each time sample using leastsquares inversion. Finally, these estimates are inverted using a poststack type scheme. Time (ms) t R R R AI SI D Reflectivities Angle 1 N Generalized inverse weight = matrix 1 R R P P ( θ 1 ) ( θ ) N Observations
11 Simultaneous Prestack Seismic inversion Prestack inversion is also based on an extension of the Fatti formulation of the AkiRichards equation: S( θ ) = c1w ( θ ) DLP + c2w ( θ ) DLS + c3w ( θ ) DLD, where : S(θ) = seismic trace at angleθ, L, L, L = logarithms of Z, Z W(θ) = P S, and ρ, the extracted wavelet at angleθ, and D is the derivative operation. D P P As in our discussion of AVO and independent inversion, this can again be set up as a leastsquares problem: model parameters = generalized inverse x observations As we discussed earlier, there are two main types of prestack inversion, deterministic and stochastic.
12 Deterministic vs Stochastic Inversion First of all, let us define the fundamental difference between deterministic and stochastic inversion: In deterministic inversion we produce what we consider to be a single best solution. In stochastic inversion we produce many possible solutions, all plausible, which average to the deterministic solution. The advantage of deterministic inversion is that we get the best leastsquares solution to our problem. The advantages of stochastic inversion are its higher frequency nature and the calculation of uncertainty.
13 Deterministic prestack inversion example On the next slide, I will show an example of deterministic prestack inversion. A Gulf Coast dataset (shown on the left of the slide) was inverted for P impedance, Simpedance and density (which are shown on the right). The inverted volumes were transformed to Vshale, porosity and Sw (also shown on the right of the slide). Our assumption is that each inverted or transformed result is the correct answer. However, this will not allow us to obtain uncertainty estimates from of the rock properties.
14 Deterministic prestack inversion example Seismic Amplitude Map Inverted Acoustic Impedance
15 Deterministic prestack inversion example Seismic Amplitude Map Inverted Inverted Acoustic Shear Impedance Impedance
16 Deterministic prestack inversion example Seismic Amplitude Map Inverted Inverted Acoustic Inverted Shear Impedance Density Impedance
17 Deterministic prestack inversion example Seismic Amplitude Map Derived V shale Map Inverted Inverted Acoustic Inverted Derived Shear Impedance Density Vshale Impedance
18 Deterministic prestack inversion example Seismic Amplitude Map Derived V shale Map Inverted Inverted Acoustic Inverted Derived Shear Impedance Density Porosity Vshale Impedance
19 Deterministic prestack inversion example Seismic Amplitude Map Derived V shale Map Inverted Inverted Acoustic Inverted Derived Shear Impedance Density Porosity Vshale Sw Impedance
20 Stochastic inversion In stochastic inversion, the leastsquares inversion method is extended by formulating the problem using a Gaussian or Log Gaussian posterior probability density function, or pdf (Tarantola, 1987). This allows us to sample various scenarios from the pdf using the Monte Carlo (MC) or Markov Chain Monte Carlo (MCMC) approach. The earliest approach to stochastic inversion was by Haas and Dubrule, 1994, in which Sequential Gaussian Simulation (SGS) is used. Buland and Omre (2003) developed a fast approach to stochastic linearized inversion which utilized a Gaussian pdf. The GeoSI method that I will discuss today combines both a Gaussian pdf and the SGS approach (Doyen, Williamson et al., 2007) My colleague Ali Tehrani discussed the Jason StatMod approach yesterday.
21 Geostatistical inversion (Haas and Dubrulle) Actual seismic trace Adapted from Dubrule, 2003 (x,y) Best simulated synthetic trace wavelet AI simulations Populate model with AI data at wells Define a random path through all (x,y) trace locations At each trace location perform a local optimization Generate a large number of trial AI sequences using SGS with spatial and vertical variograms. Compute reflectivity series and convolve with extracted wavelet. Compute misfit against observed seismic. Retain best matching AI (ρ >0.8). Go to next trace
22 Variogram models Vertical (temporal) variogram Horizontal variogram map showing anisotropy Anisotropic variograms in principal directions Here are the variograms computed by Haas and Dubrulle (1994), showing the vertical, or temporal change, and the horizontal change including anisotropy.
23 Bayesian stochastic inversion Although geostatistical stochastic inversion produces reasonable results, it has two limitations: It is quite slow. It has difficulty in converging to an answer. Buland and Omre (2003) introduced a new type of stochastic inversion which was based on multivariate Bayesian statistics. To illustrate the concept of Bayesian statistics, I will first consider the standard leastsquares regression problem. We will then look at the general theory proposed by Buland and Omre. We will then extend this method by combining it with SGS.
24 Leastsquares regression Consider a regression fit to 16 measured porosity values (φ i ) plotted against seismic impedance (z i ), shown by the red line in the plot. This can be written: φi = a + bz i The regression line is the leastsquares fit between porosity and impedance and is considered the right answer, even though only one point falls on it.
25 Statistical interpretation In the statistical interpretation of this crossplot, each variable (porosity and impedance) is given as a Gaussian probability distribution function (pdf) defined by its mean (µ) and variance (σ). p(φ) µ φ σ φ p(z) σ z Joint pdf p(φ,z) µ z The joint pdf p(φ,z) is the probability of φ and z occurring, and is defined by the variances and means, as well as the covariance between φ and z.
26 Bayesian regression Bayesian statistics tells us that the conditional probability of φ given z, or the posterior, equals the joint probability divided by the probability of z, or the prior. p( φ z) = p( φ, z) p( z) σ φ µ φ z conditional pdf p(φ z) σ φ z The conditional mean µ φ z is the leastsquares fit, and the conditional variance σ φ z gives us the scatter in this fit. Note it is narrower than p(φ).
27 Bayesian stochastic inversion Generalizing the previous example to inversion, Buland and Omre (2003) showed that: where: µ C m d m and : µ = model T 1 1 = C ( G C d + C µ ) m d m d d m m = conditional mean, C C m d = ( G covariance, T C 1 d G + C 1 m = conditional This equation reduces to the leastsquares solution if we assume that µ m = 0, and C d = σ d2 I: T 2 1 mˆ = µ m d = ( G G + σ dcm ) d µ ) = data m 1 = model covariance mean (prior), 1 G covariance. T, d
28 GeoSI The GeoSI method, as implemented by CGG and ported to the Hampson Russell suite of software, involves the following steps: Build a stratigraphic grid using horizons, well logs and layerbased kriging. Bring in partial angle stacks and wavelets. Compute the Bayesian posterior distribution by combining the model, seismic data and well logs. Create multiple P and Simpedance realizations using the SGS technique. Compute the mean and standard deviations from the impedance realizations. These steps are shown diagrammatically on the next two slides.
29 Building the stratigraphic grid Stratigraphic grid Ip Horizons in time Layerbased Kriging Is Stratigraphic layering style Lowpass filtering Well logs in time (Vp, Vs, Density) Lowfrequency prior model in stratigraphic grid R. Moyen and J. Frelet
30 Stochastic Inversion Workflow Partial angle stack seismic cubes n IpIs realisations Well logs (Vp, Vs, Density) Well uncertainty Bayesian stochastic inversion time AI IpIs prior mean & standard deviation in stratigraphic grid Horizontal & vertical variograms R. Moyen and J. Frelet Posterior mean & standard deviation IpIs
31 Bandwidth components For all inversion methods, the prior model is constructed by interpolating filtered logs, and controls low frequencies. For both deterministic and stochastic inversion, the seismic amplitudes control intermediate frequencies within the seismic bandwidth. In stochastic inversion, the vertical variogram model controls the high frequencies. Prior model Seismic Variogram model Power Spectrum Frequency (Hz) Adapted from Moyen and Frelet
32 Offshore West Africa example Elastic inversion (IpIs) 3 seismic angle stacks ,000 traces Time window of 200 ms 132 layers in grid 500 realisations (59 Gb total) 3 wells with Vp, Vs and density logs Computations on standard workstation Courtesy of R. Moyen and J. Frelet
33 IpIs Prior Model P Impedance 200 ms m/s x g/cm 3 S Impedance Courtesy of R. Moyen and J. Frelet m/s x g/cm 3
34 IpIs Realisations P Impedance 200 ms 5200 m/s x g/cm S Impedance Courtesy of R. Moyen and J. Frelet 2200 m/s x g/cm3 4200
35 Ip Posterior Mean and Standard Deviation P Impedance mean 200 ms m/s x g/cm 3 P Impedance std. dev. Courtesy of R. Moyen and J. Frelet σ Ip (km/s. g/cm 3 )
36 Vp/Vs Mean vs Realisations 1.5 Vp/Vs 2.5 Mean Realisation Posterior mean Sand/shale Cutoff Realisation
37 Inversion Results Vp/Vs Ratio Mean of 500 realisations Vp/Vs Courtesy of R. Moyen and J. Frelet
38 Using the realizations One of the key questions about stochastic inversion is: what do we do with all the realizations? In other words, wouldn t a single answer (i.e. the deterministic solution) be better? The answer is that with multiple realizations we can generate a number of new results, such as: Seismic lithology prediction. Facies classification. Volumetric uncertainty analysis. Petrophysical property analysis. These concepts are illustrated in the next few slides.
39 Stochastic Lithology Prediction N sand / shale simulations Is Ip Sand probability cube N realisations of Ip, Is Histogram of sand volume Courtesy of R. Moyen and J. Frelet
40 Facies Discrimination 0.5 Poisson s ratio Ip/Is Courtesy of R. Moyen and J. Frelet VSH P Impedance
41 Individual Realisations Courtesy of R. Moyen and J. Frelet
42 Histogram of sand volume For each realization, we can compute the sand volume from the number of cells with sand. This can then be arranged in histogram format, and the probability percentiles can be computed. A percentile is computed from the total area under the probability curve. Note that the percentile maps do not indicate a higher probability of sand, only where the map falls within the distribution. As shown by an earlier slide, the percentile values will in general be larger than mean computation. These concepts are illustrated in the next few slides.
43 Histogram of Sand Volume from Realizations 40 P50 Number of realisations P10 P90 0 Sand volume P10 P50 P90 Ranked lithology simulations
44 Sand Volume from Realisations and Mean 40 P50 Number of realisations Sand volume from inversion mean P10 P90 0 Sand volume
45 Connected Sand Geobodies Geobodies connected to at least one well Courtesy of R. Moyen and J. Frelet Colorcode: geobody volume (only largest are displayed)
46 Facies Probability from Stochastic Inversion Sand probability Large volume but small probability Courtesy of R. Moyen and J. Frelet Smaller volume but high probability
47 Stochastic Petrophysical Modelling Multiple I p & I s models Multiple V sh and Φ models Φ I p Statistical petroelastic calibration Courtesy of R. Moyen and J. Frelet
48 Geostatistical inversion vs other modelling techniques Geostatistical reservoir modeling Interpolate between the wells Plausible details Accurate near wells Not elsewhere Deterministic seismic inversion Optimize PImpedance to minimize synthetictoseismic misfit Accurate within seismic bandwidth Unrealistically smooth Only one possibility StatMod/GeoSI geostatistical seismic inversion Subsumes geostatistical modeling and deterministic inversion Does both, simultaneously and in a statistically rigorous way Multiple plausible realizations at high detail (e.g. 1ms 25m) Yet also coherent interpretations of the seismic up to the km scale
49 Conclusions Stochastic inversion is a natural extension of deterministic inversion (mean of realizations deterministic inversion) But it can provide extra information, such as: Lithology probability Facies distribution Volumetrics Petrophysical parameters Our case study focussed on a channel sand play from West Africa.
50 Thank You
Quantitative Interpretation
Quantitative Interpretation The aim of quantitative interpretation (QI) is, through the use of amplitude analysis, to predict lithology and fluid content away from the well bore. This process should make
More informationStatistical Rock Physics
Statistical  Introduction Book review 3.13.3 Min Sun March. 13, 2009 Outline. What is Statistical. Why we need Statistical. How Statistical works Statistical Rock physics Information theory Statistics
More informationSEISMIC INVERSION OVERVIEW
DHI CONSORTIUM SEISMIC INVERSION OVERVIEW Rocky Roden September 2011 NOTE: Terminology for inversion varies, depending on the different contractors and service providers, emphasis on certain approaches,
More informationReservoir Characterization using AVO and Seismic Inversion Techniques
P205 Reservoir Characterization using AVO and Summary *Abhinav Kumar Dubey, IIT Kharagpur Reservoir characterization is one of the most important components of seismic data interpretation. Conventional
More informationAFI (AVO Fluid Inversion)
AFI (AVO Fluid Inversion) Uncertainty in AVO: How can we measure it? Dan Hampson, Brian Russell HampsonRussell Software, Calgary Last Updated: April 2005 Authors: Dan Hampson, Brian Russell 1 Overview
More informationComparative Study of AVO attributes for Reservoir Facies Discrimination and Porosity Prediction
5th Conference & Exposition on Petroleum Geophysics, Hyderabad004, India PP 49850 Comparative Study of AVO attributes for Reservoir Facies Discrimination and Porosity Prediction Y. Hanumantha Rao & A.K.
More informationReliability of Seismic Data for Hydrocarbon Reservoir Characterization
Reliability of Seismic Data for Hydrocarbon Reservoir Characterization Geetartha Dutta (gdutta@stanford.edu) December 10, 2015 Abstract Seismic data helps in better characterization of hydrocarbon reservoirs.
More informationURTeC: 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 informationQuantitative Seismic Interpretation An Earth Modeling Perspective
Quantitative Seismic Interpretation An Earth Modeling Perspective Damien Thenin*, RPS, Calgary, AB, Canada TheninD@rpsgroup.com Ron Larson, RPS, Calgary, AB, Canada LarsonR@rpsgroup.com Summary Earth models
More informationLinearized AVO and Poroelasticity for HRS9. Brian Russell, Dan Hampson and David Gray 2011
Linearized AO and oroelasticity for HR9 Brian Russell, Dan Hampson and David Gray 0 Introduction In this talk, we combine the linearized Amplitude ariations with Offset (AO) technique with the BiotGassmann
More informationReservoir properties inversion from AVO attributes
Reservoir properties inversion from AVO attributes Xingang Chi* and Dehua Han, University of Houston Summary A new rock physics model based inversion method is put forward where the shalysand mixture
More informationAVAZ and VVAZ practical analysis to estimate anisotropic properties
AVAZ and VVAZ practical analysis to estimate anisotropic properties Yexin Liu*, SoftMirrors Ltd., Calgary, Alberta, Canada yexinliu@softmirrors.com Summary Seismic anisotropic properties, such as orientation
More informationIntegration of Rock Physics Models in a Geostatistical Seismic Inversion for Reservoir Rock Properties
Integration of Rock Physics Models in a Geostatistical Seismic Inversion for Reservoir Rock Properties Amaro C. 1 Abstract: The main goal of reservoir modeling and characterization is the inference of
More informationPorosity prediction using cokriging with multiple secondary datasets
Cokriging with Multiple Attributes Porosity prediction using cokriging with multiple secondary datasets Hong Xu, Jian Sun, Brian Russell, Kris Innanen ABSTRACT The prediction of porosity is essential for
More informationPreStack Seismic Inversion and Amplitude Versus Angle Modeling Reduces the Risk in Hydrocarbon Prospect Evaluation
Advances in Petroleum Exploration and Development Vol. 7, No. 2, 2014, pp. 3039 DOI:10.3968/5170 ISSN 1925542X [Print] ISSN 19255438 [Online] www.cscanada.net www.cscanada.org PreStack Seismic Inversion
More informationPorosity. Downloaded 09/22/16 to Redistribution subject to SEG license or copyright; see Terms of Use at
Geostatistical Reservoir Characterization of Deepwater Channel, Offshore Malaysia Trisakti Kurniawan* and Jahan Zeb, Petronas Carigali Sdn Bhd, Jimmy Ting and Lee Chung Shen, CGG Summary A quantitative
More informationBayesian lithology/fluid inversion comparison of two algorithms
Comput Geosci () 14:357 367 DOI.07/s59600991559 REVIEW PAPER Bayesian lithology/fluid inversion comparison of two algorithms Marit Ulvmoen Hugo Hammer Received: 2 April 09 / Accepted: 17 August 09 /
More informationRock Physics and Quantitative Wavelet Estimation. for Seismic Interpretation: Tertiary North Sea. R.W.Simm 1, S.Xu 2 and R.E.
Rock Physics and Quantitative Wavelet Estimation for Seismic Interpretation: Tertiary North Sea R.W.Simm 1, S.Xu 2 and R.E.White 2 1. Enterprise Oil plc, Grand Buildings, Trafalgar Square, London WC2N
More informationPetrophysical Study of Shale Properties in Alaska North Slope
Petrophysical Study of Shale Properties in Alaska North Slope Minh Tran Tapan Mukerji Energy Resources Engineering Department Stanford University, CA, USA Region of Interest 1.5 miles 20 miles Stratigraphic
More informationAn empirical study of hydrocarbon indicators
An empirical study of hydrocarbon indicators Brian Russell 1, Hong Feng, and John Bancroft An empirical study of hydrocarbon indicators 1 HampsonRussell, A CGGVeritas Company, Calgary, Alberta, brian.russell@cggveritas.com
More informationA021 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 informationDeterministic and stochastic inversion techniques used to predict porosity: A case study from F3Block
Michigan Technological University Digital Commons @ Michigan Tech Dissertations, Master's Theses and Master's Reports 2015 Deterministic and stochastic inversion techniques used to predict porosity: A
More informationRock physics and AVO analysis for lithofacies and pore fluid prediction in a North Sea oil field
Rock physics and AVO analysis for lithofacies and pore fluid prediction in a North Sea oil field Downloaded 09/12/14 to 84.215.159.82. Redistribution subject to SEG license or copyright; see Terms of Use
More informationTowards Interactive QI Workflows Laurie Weston Bellman*
Laurie Weston Bellman* Summary Quantitative interpretation (QI) is an analysis approach that is widely applied (Aki and Richards, 1980, Verm and Hilterman, 1995, Avseth et al., 2005, Weston Bellman and
More informationMaximize the potential of seismic data in shale exploration and production Examples from the Barnett shale and the Eagle Ford shale
Maximize the potential of seismic data in shale exploration and production Examples from the Barnett shale and the Eagle Ford shale Joanne Wang, Paradigm Duane Dopkin, Paradigm Summary To improve the success
More informationPrinciples of 3D Seismic Interpretation and Applications
Principles of 3D Seismic Interpretation and Applications Instructor: Dominique AMILHON Duration: 5 days Level: IntermediateAdvanced Course Description This course delivers techniques related to practical
More informationTechnical Note: Coloured, Deterministic & Stochastic Inversion
Technical Note: Coloured, Deterministic & Stochastic Inversion http://www.sorviodvnvm.co.uk 1 Earthworks Ltd  Company Overview Earthworks Environment & Resources Ltd. was formed in June 2001 to provide
More informationSPECIAL Latin SECTION: America L a t i n A m e r i c a
SPECIAL Latin SECTION: America L a t i n A m e r i c a Seismic inversion and AVO analysis applied to predictivemodeling gascondensate sands for exploration and early production in the Lower Magdalena
More informationDATA ANALYSIS AND INTERPRETATION
III. DATA ANALYSIS AND INTERPRETATION 3.1. Rift Geometry Identification Based on recent analysis of modern and ancient rifts, many previous workers concluded that the basic structural unit of continental
More informationElastic impedance inversion from robust regression method
Elastic impedance inversion from robust regression method Charles Prisca Samba 1, Liu Jiangping 1 1 Institute of Geophysics and Geomatics,China University of Geosciences, Wuhan, 430074 Hubei, PR China
More informationA E. SEG/San Antonio 2007 Annual Meeting. exp. a V. a V. Summary
Timelapse simulatortoseismic study  Forties field, North Sea. Christophe Ribeiro *, Cyrille Reiser, Philippe Doyen, CGGeritas, London, UK August Lau, Apache Corp., Houston, US, Steve Adiletta, Apache
More informationTHE USE OF SEISMIC ATTRIBUTES AND SPECTRAL DECOMPOSITION TO SUPPORT THE DRILLING PLAN OF THE URACOABOMBAL FIELDS
THE USE OF SEISMIC ATTRIBUTES AND SPECTRAL DECOMPOSITION TO SUPPORT THE DRILLING PLAN OF THE URACOABOMBAL FIELDS Cuesta, Julián* 1, Pérez, Richard 1 ; Hernández, Freddy 1 ; Carrasquel, Williams 1 ; Cabrera,
More informationNonlinear 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 informationAVAZ inversion for fracture orientation and intensity: a physical modeling study
AVAZ inversion for fracture orientation and intensity: a physical modeling study Faranak Mahmoudian*, Gary F. Margrave, and Joe Wong, University of Calgary. CREWES fmahmoud@ucalgary.ca Summary We present
More informationAmplitude variation with offset AVO. and. Direct Hydrocarbon Indicators DHI. Reflection at vertical incidence. Reflection at oblique incidence
Amplitude variation with offset AVO and Direct Hydrocarbon Indicators DHI Reflection at vertical incidence Reflection coefficient R(p) c α 1 S wavespeed β 1 density ρ 1 α 2 S wavespeed β 2 density ρ 2
More informationImproving Resolution with Spectral Balancing A Case study
P299 Improving Resolution with Spectral Balancing A Case study M Fatima, Lavendra Kumar, RK Bhattacharjee, PH Rao, DP Sinha Western Offshore Basin, ONGC, Panvel, Mumbai Summary: The resolution limit
More informationQuantitative interpretation using faciesbased seismic inversion
t Special section: Seismic inversion Quantitative interpretation using faciesbased seismic inversion Ehsan Zabihi Naeini 1 and Russell Exley 2 Abstract Quantitative interpretation (QI) is an important
More informationPractical aspects of AVO modeling
Practical aspects of AVO modeling YONGYI LI, Paradigm Geophysical, Calgary, Canada JONATHAN DOWNTON, Veritas, Calgary, Canada, YONG XU, Arcis Corporation, Calgary, Canada AVO (amplitude variation with
More informationFluidproperty discrimination with AVO: A BiotGassmann perspective
Fluidproperty discrimination with AVO: A BiotGassmann perspective Brian H. Russell, Ken Hedlin 1, Fred J. Hilterman, and Laurence R. Lines ABSTRACT This paper draws together basic rock physics, AVO,
More informationNORGES TEKNISKNATURVITENSKAPELIGE UNIVERSITET. Directional Metropolis Hastings updates for posteriors with nonlinear likelihoods
NORGES TEKNISKNATURVITENSKAPELIGE UNIVERSITET Directional Metropolis Hastings updates for posteriors with nonlinear likelihoods by Håkon Tjelmeland and Jo Eidsvik PREPRINT STATISTICS NO. 5/2004 NORWEGIAN
More informationNeural 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 information3D and 4D Seismic Data Integration for Geomodel Infilling: A Deep Offshore Turbiditic Field Case Study.
IPTC18306MS 3D and 4D Seismic Data Integration for Geomodel Infilling: A Deep Offshore Turbiditic Field Case Study. V. Silva, T. Cadoret, L.Bergamo, and R.Brahmantio, TOTAL E&P France Copyright 2015,
More informationDerived Rock Attributes Analysis for Enhanced Reservoir Fluid and Lithology Discrimination
IOSR Journal of Applied Geology and Geophysics (IOSRJAGG) eissn: 2321 0990, pissn: 2321 0982.Volume 5, Issue 2 Ver. I (Mar.  Apr. 2017), PP 95105 www.iosrjournals.org Derived Rock Attributes Analysis
More informationBertrand Six, Olivier Colnard, JeanPhilippe Coulon and Yasmine Aziez CGGVeritas Frédéric Cailly, Total
4D Seismic Inversion: A Case Study Offshore Congo Bertrand Six, Olivier Colnard, JeanPhilippe Coulon and Yasmine Aziez CGGVeritas Frédéric Cailly, Total Summary The first 4D seismic survey in Congo was
More informationReservoir characterization
1/15 Reservoir characterization This paper gives an overview of the activities in geostatistics for the Petroleum industry in the domain of reservoir characterization. This description has been simplified
More informationImpact of Phase Variations on Quantitative AVO Analysis
Impact of Phase Variations on Quantitative AVO Analysis Summary A.K. Srivastava, V. Singh*, D.N. Tiwary and V. Rangachari GEOPIC, ONGC, Dehradun Email: ak_sri3@rediffmail.com Effectiveness of AVO techniques
More informationPETROLEUM GEOSCIENCES GEOLOGY OR GEOPHYSICS MAJOR
PETROLEUM GEOSCIENCES GEOLOGY OR GEOPHYSICS MAJOR APPLIED GRADUATE STUDIES Geology Geophysics GEO1 Introduction to the petroleum geosciences GEO2 Seismic methods GEO3 Multiscale geological analysis GEO4
More informationRESERVOIR SEISMIC CHARACTERISATION OF THIN SANDS IN WEST SYBERIA
www.senergyltd.com RESERVOIR SEISMIC CHARACTERISATION OF THIN SANDS IN WEST SYBERIA Erick Alvarez, Jaume Hernandez, Bolkhotivin E.A., Belov A.V., Hakima Ben Meradi,Jonathan Hall, Olivier Siccardi, Phil
More informationEstimation of Elastic Parameters Using TwoTerm Fatti Elastic Impedance Inversion
Journal of Earth Science, Vol. 6, No. 4, p. 556 566, August 15 ISSN 1674487X Printed in China DOI:.7/s158155645 Estimation of Elastic Parameters Using TwoTerm Fatti Elastic Impedance Inversion Jin
More informationPROCEEDINGS, INDONESIAN PETROLEUM ASSOCIATION ThirtyEighth Annual Convention & Exhibition, May 2014
IPA14G227 PROCEEDINGS, INDONESIAN PETROLEUM ASSOCIATION ThirtyEighth Annual Convention & Exhibition, May 2014 A COMPARISON OF DEPTH CONVERSION METHODS IN BUNTAL GAS FIELD, BLOCK B, NATUNA SEA, INDONESIA
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 informationUsing multicomponent seismic for reservoir characterization in Venezuela
Using multicomponent seismic for reservoir characterization in Venezuela REINALDO J. MICHELENA, MARÍA S. DONATI, ALEJANDRO A. VALENCIANO, and CLAUDIO D AGOSTO, Petróleos de Venezuela (Pdvsa) Intevep, Caracas,
More informationRock Physics of Organic Shale and Its Implication
Rock Physics of Organic Shale and Its Implication Lev Vernik, Marathon Oil Corporation, Houston, USA lvernik@marathonoil.com Yulia Khadeeva, Marathon Oil Corporation, Houston, USA Cris Tuttle, Marathon
More informationSadewa Field is in Kutei Basin in the Makassar Strait between
SPECIAL Asia SECTION: Pacific A s i a P acific Distinguishing gas sand from shale/brine sand using elastic impedance data and the determination of the lateral extent of channel reservoirs using amplitude
More informationThin 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 informationStochastic inversion of seismic PP and PS data for reservoir parameter. estimation
Stochastic inversion of seismic PP and PS data for reservoir parameter estimation Shortened Title: Stochastic inversion of PP and PS data Jinsong Chen (jchen@lbl.gov) Lawrence Berkeley National Laboratory,
More informationShaly Sand Rock Physics Analysis and Seismic Inversion Implication
Shaly Sand Rock Physics Analysis and Seismic Inversion Implication Adi Widyantoro (IkonScience), Matthew Saul (IkonScience/UWA) Rock physics analysis of reservoir elastic properties often assumes homogeneity
More informationUnconventional reservoir characterization using conventional tools
Ritesh K. Sharma* and Satinder Chopra Arcis Seismic Solutions, TGS, Calgary, Canada Summary Shale resources characterization has gained attention in the last decade or so, after the Mississippian Barnett
More information2011 SEG SEG San Antonio 2011 Annual Meeting 771. Summary. Method
Geological Parameters Effecting ControlledSource Electromagnetic Feasibility: A North Sea Sand Reservoir Example Michelle Ellis and Robert Keirstead, RSI Summary Seismic and electromagnetic data measure
More informationGeological Classification of SeismicInversion Data in the Doba Basin of Chad*
Geological Classification of SeismicInversion 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 informationThe prediction of reservoir
Risk Reduction in Gas Reservoir Exploration Using Joint SeismicEM Inversion GAS EXPLORATION By Yoram Rubin, G. Michael Hoversten, Zhangshuan Hou and Jinsong Chen, University of California, Berkeley A
More informationSearch and Discovery Article #42083 (2017)** Posted August 21, 2017
Geophysical Characterization of Carbonate Reservoirs Using Advanced Interpretation Techniques: Applications to Abenaki Formation, Penobscot Block, Nova Scotia* D. Coronel 1, E. Villamizar 1, E. Illidge
More informationMultiple realizations: Model variance and data uncertainty
Stanford Exploration Project, Report 108, April 29, 2001, pages 1?? Multiple realizations: Model variance and data uncertainty Robert G. Clapp 1 ABSTRACT Geophysicists typically produce a single model,
More informationChurning seismic attributes with principal component analysis
Satinder Chopra + * and Kurt J. Marfurt + Arcis Seismic Solutions, Calgary; The University of Oklahoma, Norman Summary Seismic attributes are an invaluable aid in the interpretation of seismic data. Different
More informationGeophysical 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 informationDownloaded 09/16/16 to Redistribution subject to SEG license or copyright; see Terms of Use at
Data Using a Facies Based Bayesian Seismic Inversion, Forties Field, UKCS Kester Waters* (Ikon Science Ltd), Ana Somoza (Ikon Science Ltd), Grant Byerley (Apache Corp), Phil Rose (Apache UK) Summary The
More informationConstraining Uncertainty in Static Reservoir Modeling: A Case Study from Namorado Field, Brazil*
Constraining Uncertainty in Static Reservoir Modeling: A Case Study from Namorado Field, Brazil* Juliana F. Bueno 1, Rodrigo D. Drummond 1, Alexandre C. Vidal 1, Emilson P. Leite 1, and Sérgio S. Sancevero
More informationHeriotWatt University
HeriotWatt University HeriotWatt University Research Gateway 4D seismic feasibility study for enhanced oil recovery (EOR) with CO2 injection in a mature North Sea field Amini, Hamed; Alvarez, Erick Raciel;
More informationBandlimited impedance inversion: using well logs to fill low frequency information in a nonhomogenous model
Bandlimited impedance inversion: using well logs to fill low frequency information in a nonhomogenous model Heather J.E. Lloyd and Gary F. Margrave ABSTRACT An acoustic bandlimited impedance inversion
More informationF003 Geomodel Update Using 4D Petrophysical Seismic Inversion on the Troll West Field
F003 Geomodel Update Using 4D Petrophysical Seismic Inversion on the Troll West Field K. Gjerding* (Statoil), N. Skjei (Statoil), A. Norenes Haaland (Statoil), I. Machecler (CGGVeritas Services) & T.
More informationEffects of Fracture Parameters in an Anisotropy Model on PWave Azimuthal Amplitude Responses
PROC. ITB Eng. Science Vol. 38 B, No. 2, 2006, 159170 159 Effects of Fracture Parameters in an Anisotropy Model on PWave Azimuthal Amplitude Responses Fatkhan Program Studi Teknik Geofisika FIKTMITB
More informationAVO inversion and Lateral prediction of reservoir properties of Amangi hydrocarbon field of the Niger Delta area of Nigeria
IOSR Journal of Applied Geology and Geophysics (IOSRJAGG) eissn: 2321 0990, pissn: 2321 0982.Volume 2, Issue 2 Ver. II. (MarApr. 2014), PP 0817 AVO inversion and Lateral prediction of reservoir properties
More informationDynamic GeoScience Martyn Millwood Hargrave Chief Executive OPTIMISE SUCCESS THROUGH SCIENCE
Dynamic GeoScience Martyn Millwood Hargrave Chief Executive OPTIMISE SUCCESS THROUGH SCIENCE Agenda 1. Ikon Science Where we are now 2. Geoscience 2012 A motion picture 3. Rock physics, AVO and Inversion
More informationEvaluation of Rock Properties from Logs Affected by Deep Invasion A Case Study
Evaluation of Rock Properties from Logs Affected by Deep Invasion A Case Study Jahan Zeb a, Reece Murrell b a CGG Australia, 1 Ord Street, West Perth, WA 6005 Contact email: Jahan.Zeb@cgg.com b Esso Australia,
More informationA Case Study on Simulation of Seismic Reflections for 4C Ocean Bottom Seismometer Data in Anisotropic Media Using Gas Hydrate Model
A Case Study on Simulation of Seismic Reflections for 4C Ocean Bottom Seismometer Data in Anisotropic Media Using Gas Hydrate Model Summary P. Prasada Rao*, N. K. Thakur 1, Sanjeev Rajput 2 National Geophysical
More informationStochastic rock physics modeling for seismic anisotropy
Stochastic rock physics modeling for seismic anisotropy Yunyue (Elita) Li, Biondo Biondi, Dave Nichols, Gary Mavko, and Robert Clapp ABSTRACT Anisotropic model building using surface seismic data is a
More informationAVO attribute inversion for finely layered reservoirs
GEOPHYSICS, VOL. 71, NO. 3 MAYJUNE 006 ; P. C5 C36, 16 FIGS., TABLES. 10.1190/1.197487 Downloaded 09/1/14 to 84.15.159.8. Redistribution subject to SEG license or copyright; see Terms of Use at http://library.seg.org/
More informationResearch  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 informationReservoir Characterization for Shales: a Barnett Shale Case Study
GUSS1424 Reservoir Characterization for Shales: a Barnett Shale Case Study JOHN PENDREL CGG GeoSoftware ALFONSO IUNIO MARINI Eni E&P This paper has been selected for presentation for the 2014 Gussow Geosciences
More informationInverting hydraulic heads in an alluvial aquifer constrained with ERT data through MPS and PPM: a case study
Inverting hydraulic heads in an alluvial aquifer constrained with ERT data through MPS and PPM: a case study Hermans T. 1, Scheidt C. 2, Caers J. 2, Nguyen F. 1 1 University of Liege, Applied Geophysics
More informationQUANTITATIVE ANALYSIS OF SEISMIC RESPONSE TO TOTALORGANICCONTENT AND THERMAL MATURITY IN SHALE GAS PLAYS
E: infoikonscience.com W: www.ikonscience.com QUANTITATIVE ANALYSIS OF SEISMIC RESPONSE TO TOTALORGANICCONTENT AND THERMAL MATURITY IN SHALE GAS PLAYS Ebrahim Zadeh 12, Reza Rezaee 1, Michel Kemper 2
More informationCorrelation of brittleness index with fractures and microstructure in the Barnett Shale
Correlation of brittleness index with fractures and microstructure in the Barnett Shale Z. Guo (British Geological Survey), M. Chapman (University of Edinburgh), X.Y. Li (British Geological Survey) SUMMARY
More informationRock 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 informationBest practices predicting unconventional reservoir quality
Introduction Best practices predicting unconventional reservoir quality Cristian Malaver, Michel Kemper, and Jorg Herwanger 1 Unconventional reservoirs have proven challenging for quantitative interpretation
More informationKYLE T. SPIKES Research Group
KYLE T. SPIKES Research Group 1 WELCOME TO THE MEETING Welcome everyone! Currently we have 8 PhD and 3 MS students in total Presentations today cover a variety of topics in various stages of development
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 informationDelivery: an opensource Bayesian seismic inversion tool. James Gunning, CSIRO Petroleum Michael Glinsky,, BHP Billiton
Delivery: an opensource Bayesian seismic inversion tool James Gunning, CSIRO Petroleum Michael Glinsky,, BHP Billiton Application areas Development, appraisal and exploration situations Requisite information
More informationCase study: AVO analysis in a highimpedance Atoka Sandstone (Pennsylvanian), North Arkoma Basin, McIntosh County, Oklahoma
Case study: AVO analysis in a highimpedance Atoka Sandstone (Pennsylvanian), North Arkoma Basin, McIntosh County, Oklahoma MOHAMED A. EISSA, Tanta University, Egypt, and University of Oklahoma, Norman,
More informationRock physics and AVO applications in gas hydrate exploration
Rock physics and AVO applications in gas hydrate exploration ABSTRACT Yong Xu*, Satinder Chopra Core Lab Reservoir Technologies Division, 301,4003rd Ave SW, Calgary, AB, T2P 4H2 yxu@corelab.ca Summary
More informationFK Characteristics of the Seismic Response to a Set of Parallel Discrete Fractures
FK 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 informationAcoustic impedance inversion analysis: Croatia offshore and onshore case studies
Acoustic impedance inversion analysis: Croatia offshore and onshore case studies Domagoj Vukadin, mag. ing. geol. Stipica Brnada, mag. geol. SPE Conference Hungarian Section Visegrád, November, 19. 2015
More information3D Geological Modeling and Uncertainty Analysis of Pilot Pad in the Long Lake Field with Lean Zone and Shale Layer
Datapages/Search and Discovery Article #9224 GeoConvention 214, FOCUS  Adapt, Refine, Sustain Calgary, Alberta, Canada, May 1216, 214 3D Geological Modeling and Uncertainty Analysis of Pilot Pad in the
More informationAVO Attributes of a Deep Coal Seam
AVO Attributes of a Deep Coal Seam Jinfeng Ma* State Key Laboratory of Continental Dynamics, Northwest University, China jinfengma@sohu.com Igor Morozov University of Saskatchewan, Saskatoon, SK, Canada
More informationAcoustic impedance inversion and CO 2 flood detection at the Alder Flats ECBM project
ECBM CO 2 flood detection Acoustic impedance inversion and CO 2 flood detection at the Alder Flats ECBM project Jason M. McCrank and Don C. Lawton ABSTRACT The 3D poststack vertical component data from
More informationEffects of fluid changes on seismic reflections: Predicting amplitudes at gas reservoir directly from amplitudes at wet reservoir
GEOPHYSICS, VOL. 77, NO. 4 (JULYAUGUST 2012); P. D129 D140, 15 FIGS., 2 TABLES. 10.1190/GEO20110331.1 Effects of fluid changes on seismic reflections: Predicting amplitudes at gas reservoir directly
More informationFacies Classifications for Seismic Inversion
Facies Classifications for Seismic Inversion Jeremy Gallop Ikon Science Canada Summary Identifying facies for classification for a seismic inversion project is an important step where one balances computational
More informationPwave reflection coefficients for transversely isotropic models with vertical and horizontal axis of symmetry
GEOPHYSICS, VOL. 62, NO. 3 (MAYJUNE 1997); P. 713 722, 7 FIGS., 1 TABLE. Pwave reflection coefficients for transversely isotropic models with vertical and horizontal axis of symmetry Andreas Rüger ABSTRACT
More informationMultiattribute seismic analysis on AVOderived parameters A case study
Multiattriute seismic analysis on AVOderived parameters A case study SATINDER CHOPRA, Core La Reservoir Technologies, Calgary, Canada DOUG PRUDEN, GEDCO, Calgary, Canada Prospecting for reservoir zones
More informationFeasibility and design study of a multicomponent seismic survey: Upper Assam Basin
P276 Summary Feasibility and design study of a multicomponent seismic survey: Upper Assam Basin K.L.Mandal*, R.K.Srivastava, S.Saha, Oil India Limited M.K.Sukla, Indian Institute of Technology, Kharagpur
More informationDHI Analysis Using Seismic Frequency Attribute On FieldAN Niger Delta, Nigeria
IOSR Journal of Applied Geology and Geophysics (IOSRJAGG) eissn: 2321 0990, pissn: 2321 0982.Volume 1, Issue 1 (May. Jun. 2013), PP 0510 DHI Analysis Using Seismic Frequency Attribute On FieldAN Niger
More information