SEISMIC INVERSION OVERVIEW

 Arthur Potter
 9 months ago
 Views:
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: prestack vs. poststack, 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 Prestack and/or poststack 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 nonunique.
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 lenslike 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 nonreservoirs 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) ModelBased 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 bestfit 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. Nonlinear Sparse Spike Inversion (Constrained Sparse Spike InversionCSSI) employs an optimization loop where the wavelet is updated (nonlinearly) so that the mismatch between the synthetic and real seismic is minimized.
21 Sparse Spike Inversion Saxena and Bhatnagar, 2008
22 ModelBased Inversion ModelBased Inversion typically starts with a low frequency model of the Pimpedance 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 ModelBased Inversion Veeken and Da Silva, 2004
24 ModelBased 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 ModelBased 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 Pwave velocity, Swave 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 AkiRichards (Shuey) linear Zoeppritz approximation: to an elastic impedance relationship: Combining the two expressions gives elastic impedance: Variations of this approach include using the Shuey twoterm 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 2term Zoeppritz linear approximation limits the range at which reflectivities can be defined. If this approximation is rewritten 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 Pwave, 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 Macromodels Simultaneous Inversion Constrained Modeldriven Global (SA) Pwave Impedance Model Swave Impedance Model Simultaneous Inversion Workflow Courtesy Equipoise
37 AVO Simultaneous Inversion Russell and Hampson, 2006 Simultaneous PImpedance Inversion
38 AVO Simultaneous Inversion Simultaneous PImpedance Inversion Vp/Vs ratio from dividing P and Simpedance sections from simultaneous inversion Russell and Hampson, 2006
39 Predicting Other Rock Properties from PreStack Inversion Prestack 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 PreStack 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 ImpedanceSimultaneous 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 Prestack waveform inversion Multicomponent inversion Multiazimuth 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, 2D and 3D 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 tradeoff 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 twoterm 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 Prestack Inversion Methods. Brian Russell
Stochastic vs Deterministic Prestack Inversion Methods Brian Russell Introduction Seismic reservoir analysis techniques utilize the fact that seismic amplitudes contain information about the geological
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 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 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 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 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 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 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 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 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 informationQuantifying Bypassed Pay Through 4D PostStack Inversion*
Quantifying Bypassed Pay Through 4D PostStack 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 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 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 informationDownloaded 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 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 informationV003 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 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 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 informationSpecial 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 informationSeismic techniques for imaging fractures, cracks and faults in the Earth. Michael Kendall
Seismic techniques for imaging fractures, cracks and faults in the Earth Michael Kendall Issues and Challanges Geometry (aspect ratio, size, orientation, density) Nonuniqueness (e.g., single set of aligned
More informationA simple algorithm for bandlimited impedance inversion
Impedance inversion algorithm A simple algorithm for bandlimited impedance inversion Robert J. Ferguson and Gary F. Margrave ABSTRACT This paper describes a seismic inversion method which has been cast
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 informationSeismic 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 informationRelinquishment 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 informationParameter 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 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 informationMain 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 informationROCK 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 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 informationInterpretation 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 informationA 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 informationChapter 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 informationA new method for multiexponential 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 multiexponential inversion of NMR relaxation measurements WANG Zhongdong 1, 2, XIAO Lizhi 1 & LIU Tangyan
More informationCrosswell 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 informationTh Guided Waves  Inversion and Attenuation
Th0108 Guided Waves  Inversion and Attenuation D. Boiero* (WesternGeco), C. Strobbia (WesternGeco), L. Velasco (WesternGeco) & P. Vermeer (WesternGeco) SUMMARY Guided waves contain significant information
More informationAcoustic 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 informationAPPLICATION 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 informationIPTC 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 informationPressure Prediction and Hazard Avoidance through Improved Seismic Imaging
Pressure Prediction and Hazard Avoidance through Improved Seismic Imaging 12121600202 Bill Barkhouse SEG Advanced Modeling Corporation (SEAM) RPSEA / SPEGCS UltraDeepwater Technology Conference September
More informationImpact 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 informationChapter 6. Conclusions. 6.1 Conclusions and perspectives
Chapter 6 Conclusions 6.1 Conclusions and perspectives In this thesis an approach is presented for the insitu characterization of rocks in terms of the distribution of hydraulic parameters (called SBRC
More informationWell Evaluation. An Integrated WellEvaluation Process for Shale Gas Reservoirs
Well Evaluation An Integrated WellEvaluation Process for Shale Gas Reservoirs APPLICATION Gas shale analysis for unconventional wells BENEFITS Optimized stimulation by thorough study and planning Increased
More informationTUTORIAL 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 informationGRAIN 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 informationGeologic 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 informationQuantitative 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 information3D 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 informationAN 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 informationNuclear 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 informationAvailable 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, GHGT13, 1418 November 2016, Lausanne,
More informationCopyright. 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 CanaWoodford Shale using
More informationOGIP 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 information10810: Advanced Algorithms and Models for Computational Biology. Optimal leaf ordering and classification
10810: 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 informationHague and London Oil Plc
Hague and London Oil Plc Exploiting Shallow Gas Offshore Netherlands DeGassing Shallow Hazards or DeRisking Shallow Gas November 2016 The Shallow Gas Play can be Summarized as Schematic evolution of
More informationIntegrating 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 informationFrequently Asked Questions
Frequently Asked Questions 1 Is THERMOVISION TOMOGRAPHY technology independent, or it requires previous seismic or other data? THERMOVISION TOMOGRAPHY technology is a standalone technology. Independent
More informationIDENTIFICATION 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 informationDetermination of the signature of a dynamite source using source scaling, Part 2: Experiment
GEOPHYSICS, VOL. 58, NO. 8 (AUGUST 1993); P. 11831194, 12 FIGS. Determination of the signature of a dynamite source using source scaling, Part 2: Experiment Anton Ziolkowski* and Karel Bokhorst ABSTRACT
More informationDetermination 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, 99111 ISSN: 17929040 (print), 17929660 (online) Scienpress Ltd, 2016 Determination of Incompressibility, Elasticity and the
More informationDelivery: 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 informationLAMINATION 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 informationModified Kriging: Evaluation, Modification, Recommendations
UCRLID 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 informationThe effect of interpolation on imaging and AVO: A Viking case study
GEOPHYSICS, VOL. 75, NO. 6 NOVEMBERDECEMBER 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 informationStudy on Rayleigh Wave Inversion for Estimating Shearwave Velocity Profile
PROC. ITB Eng. Science Vol. 35 B, No. 1, 2003, 1122 11 Study on Rayleigh Wave Inversion for Estimating Shearwave Velocity Profile T.A. Sanny*, Yoes Avianto* *Geophysical Engineering Department, Faculty
More informationwith wideangle prestack migration image crustal faults using only P, P, isotropic scattering characteristics and
34 Chapter 3 Isotropic scattering and crustal imaging with wideangle prestack migration 3.1 Introduction Earthquake seismologists may be reluctant to believe that one can successfully image crustal faults
More informationImportance 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 informationFULL MOMENT TENSOR ANALYSIS USING FIRST MOTION DATA AT THE GEYSERS GEOTHERMAL FIELD
PROCEEDINGS, ThirtyEighth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 1113, 2013 SGPTR198 FULL MOMENT TENSOR ANALYSIS USING FIRST MOTION DATA AT
More informationMicroseismic 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 informationEnhanced 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 informationSearch and Discovery Article #51409 (2017)** Posted August 7, Abstract. Selected References
Saturations of Migrating Buoyant Fluids from Invasion Percolation Flow Simulation Using SmallScale, High Resolution Geologic Models With Realistic Heterogeneity* Timothy A. Meckel 1, Luca Trevisan 2,
More informationDetermination 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 informationSearch 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 PreTertiary Fractured Granite Basement Reservoir in Cuu Long
More informationand 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 informationExploration 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, TimeLine and Status 2. A few statistics 3. Reasons
More informationRelinquishment 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 informationHeterogeneity 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 informationReceiver Function Inversion
Receiver Function Inversion Advanced Studies Institute on Seismological Research Kuwait City, Kuwait  January 1922, 2013 Jordi Julià Universidade Federal do Rio Grande do Norte, Brasil Outline Introduction
More informationApplication 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 informationIMPEDANCE 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 informationBorehole Geophysics. Acoustic logging measurements
Acoustic logging measurements  Review of basic physics background  Concept of P and Swave measurements and logging tools  Tube waves  Seismic imaging  Synthetic seismograms  Field application examples
More informationASSESSMENT 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 informationEXTENDED 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 informationBodywave radiation patterns and AVO in transversely isotropic media
GEOPHYSICS, VOL. 60, NO. 5 (SEPTEMBEROCTOBER 1995); P. 14091425, 11 FIGS. Bodywave radiation patterns and AVO in transversely isotropic media llya Tsvankin* ABSTRACT The angular dependence of reflection
More informationBayesian 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 21676 Pasadena, CA 91125 lewickiocns.caltech.edu Abstract
More informationIBIS  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 informationNumerical 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 informationRelinquishment 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 informationBayesian 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 informationECE 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 informationP11 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 informationThis paper was prepared for presentation at the Unconventional Resources Technology Conference held in San Antonio, Texas, USA, 13 August 2016.
URTeC: 2444366 Using DepletionZone Microseismicity to Understand Producing Volumes Jonathan P. McKenna*, Michael H. Grealy, Michael S. Blaz and Nathan M. Toohey, MicroSeismic, Inc. Copyright 2016, Unconventional
More information25th Seismic Research Review  Nuclear Explosion Monitoring: Building the Knowledge Base
Pn Q UNDER TIBET AND TIENSHAN WITH PRACTICAL AND SCIENTIFIC IMPLICATIONS Jiakang Xie LamontDoherty Earth Observatory, Columbia University Sponsored by Defense Threat Reduction Agency Contract No. DTRA0100C0048
More informationAnisotropic velocity changes in seismic timelapse data Jorg Herwanger*, Ed Palmer (WesternGeco), Christian Rau Schiøtt (Hess)
Anisotropic velocity changes in seismic timelapse data Jorg Herwanger*, Ed Palmer (WesternGeco), Christian Rau Schiøtt (Hess) Summary Reservoir production causes subsurface deformations and changes in
More informationTheory  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 spatiotemporal
More informationBayesian 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 informationSEAPEX 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 informationB. Milkereit 1, E. Adam 2, Z. Li 3, Wei Qian 1, T. Bohlen 4, D. Banerjee 1, and D.R. Schmitt 5
Multioffset vertical seismic profiling: an experiment to assess petrophysicalscale parameters at the JAPEX/JNOC/GSC et al. Mallik 5L38 gas hydrate production research well B. Milkereit 1, E. Adam 2,
More informationSimulated 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 informationL11: 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 informationADI 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