1 DHI CONSORTIUM SEISMIC INVERSION OVERVIEW Rocky Roden September 2011
2 NOTE: Terminology for inversion varies, depending on the different contractors and service providers, emphasis on certain approaches, and the goals to be accomplished. There are numerous combinations, hybrids, and variations of these inversion methods. This overview attempts to summarize the most common terminology and inversion approaches described in the industry.
3 Seismic Inversion Outline Inversion Definition and Terms Required Data Types of Seismic Inversion Role of Inversion Realistic Expectation from Inversion Inversion and New SEC Regulations
4 What is Inversion? Inversion transforms seismic reflection data into rock and fluid properties. The objective of seismic inversion is to convert reflectivity data (interface properties) to layer properties. To determine elastic parameters, the reflectivity from AVO effects must be inverted. The most basic inversion calculates acoustic impedance (density X velocity) of layers from which predictions about lithology and porosity can be made. The more advanced inversion methods attempt to discriminate specifically between lithology, porosity, and fluid effects. Inversions can be grouped into categories: pre-stack vs. post-stack, deterministic vs. geostatistical, or relative vs. absolute. Inversion is the flip side of forward modeling.
5 Common Inversion Terms Relative Impedance Relative changes in inversion, not real values. Absolute Impedance Actual acoustic impedance, contains low frequency trend. Deterministic Inversion A single output value is determined from the input. Stochastic (Probabilistic) Inversion A range of equally probable outputs are derived from the input (under controlled number of inputs). Objective Function A quantitative measure of the misfit between the observed data and the data predicted by inversion (L1 and L2 Norm). L2 Norm Most common measure of data misfit, a least squares difference (sum of the square of the data residuals). L1 Norm Misfit measure that is the sum of the absolute data residuals. Global More than one trace is inverted at the same time within a common objective function. The number of traces inverted at once depends on the algorithm. Simulated Annealing A global optimization technique based on crystal growth in a cooling volcanic melt. Computes difference between seismic and convolutional model with applied wavelet. The model is perturbed and a new model is simulated, and the difference is again measured. Differences are compared with the smallest difference accepted and the other differences accepted with a probability (Metropolis Criterion). The process is repeated until and there is a very small residual or a threshold has been reached.
6 Forward Modeling takes well logs, combines with wavelet to produce synthetic seismic trace Inversion takes seismic trace, removes effects of estimated wavelet and creates acoustic impedance values Oilfield Review, 2008
7 Data Considerations for Inversion Pre-stack and/or post-stack time migrated seismic data Well data conditioning and editing is necessary Wavelet estimation required for all modern inversion approaches -Deterministic -Statistical Low Frequency Trend absolute acoustic impedance contains a low frequency trend that must be obtained usually from well control or stacking velocities NOTE: Without wavelet estimation and well calibration, the inversion solution is non-unique.
8 Data Considerations for Inversion Pillar, 2011
9 Data Considerations for Inversion Low Frequency Trend in Red Oilfield Review, 2008
10 Data Considerations for Inversion Missing low frequencies and their modeling from well log data cause bias in inversion, regardless of the methodology. Determining low frequency issues: Lithologies encountered at the wells thin or thicken away from the wells. In fact, some lens-like possible plays between wells are often difficult or impossible to model accurately. Lateral geological changes away from the wells have no resemblance to the geological layers at wells used fro low frequency model generation, e.g., sand injections. There are excessive lateral changes in rock and reservoir properties (e.g., porosity changes and shaling out). Ozdemir, 2009
11 Data Considerations for Inversion Rock physics approach to well log preparation for inversion: Continuous reconstruction of formation properties along the wellbore, both for reservoirs and non-reservoirs Characterization of rocks in terms of their elastic properties Establishment of causal relationships between elastic and petrophysical properties such as porosity, clay content and fluid saturation NOTE: Depends on inversion approach and stage of exploration through development cycle.
13 Running Sum Running Sum is the integration or adding of seismic amplitudes on a seismic trace. It is the cumulative sum of amplitude values at any time sample starting from the top to the end of the trace. Running Sum provides an approximation to acoustic impedance (velocity X density), while lacking the low frequency trend of a full inversion this can be useful over short intervals. The less noise and the broader the bandwidth (narrower the wavelet), the better will be the approximation to impedance.
14 Running Sum Normal Amplitude Running Sum Courtesy SMT
15 Recursive Trace Integration Recursive/Trace Integration (RTI) employs the discrete recursive inversion formula (below) which indicates that the acoustic impedance of a particular layer and the reflection coefficient at its base can be used to calculate the acoustic impedance of the next layer. An estimate of the acoustic impedance of the first layer is required. Normal incidence reflectivity Recursive inversion formula r = reflection coefficient Z = impedance
16 Recursive Trace Integration Seismic section from Alberta highlighting an anomalous bright spot zone. Recursive inversion of the section above Russell and Hampson, 2006
17 Colored Inversion Colored Inversion transforms migrated seismic data into a bandlimited acoustic impedance volume by shaping the mean seismic spectrum into the impedance log spectrum. This approach only gives a relative acoustic impedance which does not contain any low frequency components. This inversion method can make it simpler to interpret the variations in thickness related to various lithologic packages, properties of thin beds, and changes in acoustic impedance associated with fluid effects and rock property variations. Veeken and Da Silva, 2004
18 Colored Inversion Colored Inversion was first introduced by Lancaster and Whitcombe (2000) and involves deriving an operator that transforms the seismic amplitude spectrum to the acoustic impedance spectrum from well logs with a -90 phase change. This typically requires crossplotting the log of the well acoustic impedances versus log of amplitudes to derive a best-fit line to the data. The slope of this line is used to develop an operator to apply to the seismic data to convert to acoustic impedance. The logic employed is that the gross spectral properties from the well acoustic impedances in any given field are reasonably constant. A phase rotation may be applied to the seismic data, but this approach assumes the seismic is zero phase. Colored inversion may establish a base case to compare against more sophisticated inversion techniques Courtesy SMT
20 Sparse Spike Inversion Sparse Spike Inversion assumes that seismic reflectivity is a series of large spikes embedded in a background of small spikes (assumes only large spikes are meaningful). This approach seeks the simplest possible reflectivity model of spikes, that when convolved with the wavelet produces a synthetic that matches the real data. Non-linear Sparse Spike Inversion (Constrained Sparse Spike Inversion-CSSI) employs an optimization loop where the wavelet is updated (non-linearly) so that the mis-match between the synthetic and real seismic is minimized.
21 Sparse Spike Inversion Saxena and Bhatnagar, 2008
22 Model-Based Inversion Model-Based Inversion typically starts with a low frequency model of the P-impedance and then the model is perturbed until a good fit is obtained between the seismic data and a synthetic trace using the recursive formula. The Generalized Linear Inversion (GLI) method perturbs an initial acoustic impedance estimate and determines the differences between the synthetic trace and the real data (simulated annealing and least square fits are often applied). Principal Component Analysis (PCA) method computes a standard response from which the input can be generated by applying specific weighting factors determined from well control (linear interpolation). The results of the convolution with the seismic wavelet are compared with the seismic traces and the velocity and density models are perturbed to reduce the discrepancy.
23 Model-Based Inversion Veeken and Da Silva, 2004
24 Model-Based Inversion Wavelet Initial Model Seismic Data Synthetic Data Misfit Function Stochastic Model Update No Optimisation by Simulated Annealing Convergent? Yes Model Based Inversion Utilizing Simulating Annealing Stop Courtesy Equipoise
25 Model-Based Inversion Seismic section from Alberta highlighting an anomalous bright spot zone. Model based inversion of the section above Russell and Hampson, 2006
26 AVO Elastic Impedance Connolly (1999) developed an angle dependent analogy to acoustic impedance to account for AVO effects. Elastic impedance is a function of P-wave velocity, S-wave velocity, density, and incidence angle. Input for elastic impedance is some form of angle stack (e.g., near, mid, and far). Wavelets are determined for each offset angle volume. This approach is accurate for small to moderate impedance changes.
27 AVO Elastic Impedance Connolly relates the Aki-Richards (Shuey) linear Zoeppritz approximation: to an elastic impedance relationship: Combining the two expressions gives elastic impedance: Variations of this approach include using the Shuey two-term approximation (good to ), assume Vp/Vs=2 or Vs calculation, and Extended Elastic Impedance (Whitcombe et al., 2002)
28 AVO Elastic Impedance Veeken and Da Silva, 2004
29 AVO Elastic Impedance Elastic impedance displays for different offset angles Veeken and Da Silva, 2004
30 AVO Extended Elastic Impedance Whitcombe et al (2002) describes an approach that extends beyond a typical incident angle range of 0 to 30 for seismic data to -90 to +90. The sin²ɵ term in the Shuey 2-term Zoeppritz linear approximation limits the range at which reflectivities can be defined. If this approximation is re-written as R = A + Btanx then scaled by cosx then a scaled reflectivity equation can be written: R = A cosx + B sinx. Extended Elastic Impedance logs can be generated from P-wave, S- wave, and density logs for each angle x. Once x is known then equivalent seismic sections can be generated.
31 AVO Extended Elastic Impedance EEI logs at different angles: Compressional modulus = +12 Lame s parameter = =20 Shear impedance = -50 Vp/Vs = +45 Acoustic impedance = 0 Francis and Hicks, 2006
32 AVO Extended Elastic Impedance Connolly, 2010
33 AVO Extended Elastic Impedance Connolly, 2010
34 AVO Extended Elastic Impedance Francis and Hicks, 2006
35 AVO Simultaneous Inversion Simultaneous inversion is a prestack inversion method that uses multiple offset or angle stacks. This method solves for S impedance, P impedance, and density, which are key for descriminating lithology, porosity, and fluid effects. For each input partial stack a wavelet is estimated. All models, partial stacks, and wavelets are input into a single inversion algorithm and solved simultaneously compensating for offset dependent phase, bandwidth, tuning and NMO stretch effects.
36 AVO Simultaneous Inversion NMO Gather (Full Offset) Wavelets Macro-models Simultaneous Inversion Constrained Model-driven Global (SA) P-wave Impedance Model S-wave Impedance Model Simultaneous Inversion Workflow Courtesy Equipoise
37 AVO Simultaneous Inversion Russell and Hampson, 2006 Simultaneous P-Impedance Inversion
38 AVO Simultaneous Inversion Simultaneous P-Impedance Inversion Vp/Vs ratio from dividing P and S-impedance sections from simultaneous inversion Russell and Hampson, 2006
39 Predicting Other Rock Properties from Pre-Stack Inversion Pre-stack Inversion Vp/Vs Poisson s Ratio Acoustic (Ip) Impedance Shear (Is) Impedance Elastic Impedance EI Lamé Parameter Lamé Parameter Lithology Fluid Content Porosity Pore Pressure Lamé Parameter Courtesy Equipoise
40 Pre-Stack Inversion allows discriminating lithologies and fluids Bunch and Dromgoole, 1995 Pillar, 2011
41 AVO Elastic Impedance Reservoir attributes from Elastic Impedance Inversion Veeken and Da Silva, 2004
42 Elastic Impedance Simultaneous Inversion Rasmussen et al., 2004
43 Elastic Impedance-Simultaneous Inversion Poisson s Ratio From Elastic Impedance Inversion From Simultaneous Inversion Saxena and Bhatnagar, 2008
44 Stochastic Inversion - Geostatistical Geostatistical or probabilistic inversion uses quantification of uncertainties attached to the inversion input data. Probability density functions (PDF) are defined and an earth model is simulated, which is perturbed to minimize the discrepancy between the modeled and measured seismic data (simulated annealing) producing multiple realizations. The PDF determination comes from well logs, spatial properties (variograms) and lithological distributions. It is critical that the interpreter quantify the uncertainties in a realistic way, especially in the absence or minimal well control in specific areas.
45 Stochastic Inversion - Geostatistical The workflow for performing a geostatistical inversion showing the input data and the outputted model realizations. McCrank et al., 2009
46 Stochastic Inversion - Geostatistical The acoustic impedance inversion results with the mean acoustic impedance results of the multiple realizations in the upper left. McCrank et al., 2009
47 Stochastic Inversion - Bayesian Seismic inversion that uses Baye s rule allows the information from all available measurements to be integrated into a consistent image of the reservoir and constrain these solutions based on a priori knowledge about the subsurface parameters. Prior knowledge about the model parameters usually is combined with a likelihood function, which depends on the misfit between the model response and the observed seismic data
48 Stochastic Inversion - Bayesian Pillar, 2011
49 Stochastic Inversion - Bayesian Pillar, 2011
50 Categories of Seismic Inversion Running Sum Well Control Post Stack Pre Stack Relative Imp. Absolute Imp. Deterministic Probabilistic Recursive Trace Integration Colored Inversion Sparse Spike (CSSI) Model Based Inversion Elastic Impedance Simultaneous Inversion Geostatistical Inversion Bayesian Inversion
51 Applicable Inversions for Exploration through Development Running Sum Exploration Exploitation Development Recursive Trace Integration Colored Inversion Sparse Spike (CSSI) Model Based Inversion Elastic Impedance Simultaneous Inversion Geostatistical Inversion Bayesian Inversion
52 A Few Additional Types of Inversion Simultaneous geostatistical partial stack inversion Spectral inversion Pre-stack waveform inversion Multi-component inversion Multi-azimuth inversion Inversion incorporating EM and FTG data 3D Full Waveform Inversion
53 Seismic Inversion Checklist (modified from Ikon website) 1. Check log data and edit accordingly. 2. Perform rock physics analysis to determine if inversion is useful. 3. Check seismic data for proper processing and conditioning, S/N, accurate partial stacks, etc. 4. Well ties and low frequency background model: -Do you get credible wavelets? -Are wavelets consistent across wells? -Is a broadband inversion required? 5. Be sure and employ appropriate inversion algorithm. 6. QC is paramount: -Check all wells and match of wavelets. -From impedances forward model to synthetic seismic; compare with actual seismic. -Does it make geological sense? -Is low frequency trend applied correctly? 7. Use impedance results intelligently they are not the end goal!
54 SEC Regulations and Inversion Pillar, 2011
55 SEC Regulations and Inversion Oil and gas companies may use any "reliable technology" to establish reserves volumes in addition to those established by production and flow test data. Oil and gas companies may classify proved undeveloped reserves ("PUDs") any distance from known proved reserves (rather than only in immediately offsetting units) based on a reasonable certainty standard. Use of reliable technology: A registrant is required to disclose, in general terms, the technologies used in ascertaining the reasonable certainty of the PUD locations producing hydrocarbons. Any technology stated must be field tested to demonstrate consistency and repeatability. In reviewing disclosures of the technologies used, some very simple statements that a combination of these technologies was used were noted: analogy, 2-D and 3-D seismic data, volumetric and material balance analysis, decline curves, petrophysics, and log analysis. A company does not need to disclose proprietary technologies or the mix of proprietary methods.
56 Conclusions 1) Seismic Inversion is not a unique process. 2) Seismic and well log information should be optimally conditioned for reliable inversion results. 3) There is a trade-off between work involved/cost/time and quality of the final inversion. 4) A rock physics study can help determine what the end product from inversion should be and whether this is attainable with the existing data. 5) Even when advanced inversion algorithms are chosen, a simpler deterministic inversion perhaps should be run for feasibility and as a yardstick of what can be resolved. 6) Most seismic inversion workflows are not linear and require reasonable inputs and usually numerous iterations. 7) Will the SEC accept inversion results for booking reserves?
57 Indeed, given that we can now directly invert to the acoustic properties using simultaneous inversion, some workers hold the view that conventional two-term AVO techniques are now passé. This is naïve given that deriving an accurate estimate of Poisson ratio from seismic is actually quite difficult, requiring a stringent set of data conditions. In many cases, bias of one form or another is introduced into the result, for example when merging the low frequency component. The interpreter needs to understand both conventional AVO approaches as well as the latest trends in inversion. Rob Simm What makes the wiggle waggle: a perspective on rock physics, First Break, June 2011
Stochastic vs Deterministic Pre-stack Inversion Methods Brian Russell Introduction Seismic reservoir analysis techniques utilize the fact that seismic amplitudes contain information about the geological
Statistical - Introduction Book review 3.1-3.3 Min Sun March. 13, 2009 Outline. What is Statistical. Why we need Statistical. How Statistical works Statistical Rock physics Information theory Statistics
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
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
Comput Geosci () 14:357 367 DOI.07/s596-009-9155-9 REVIEW PAPER Bayesian lithology/fluid inversion comparison of two algorithms Marit Ulvmoen Hugo Hammer Received: 2 April 09 / Accepted: 17 August 09 /
Quantifying Bypassed Pay Through 4-D Post-Stack Inversion* Robert Woock 1, Sean Boerner 2 and James Gamble 1 Search and Discovery Article #40799 (2011) Posted August 12, 2011 *Adapted from oral presentation
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,
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
NORGES TEKNISK-NATURVITENSKAPELIGE UNIVERSITET Directional Metropolis Hastings updates for posteriors with nonlinear likelihoods by Håkon Tjelmeland and Jo Eidsvik PREPRINT STATISTICS NO. 5/2004 NORWEGIAN
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
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
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
Seismic techniques for imaging fractures, cracks and faults in the Earth Michael Kendall Issues and Challanges Geometry (aspect ratio, size, orientation, density) Non-uniqueness (e.g., single set of aligned
Impedance inversion algorithm A simple algorithm for band-limited impedance inversion Robert J. Ferguson and Gary F. Margrave ABSTRACT This paper describes a seismic inversion method which has been cast
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
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,
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,
F-K Characteristics of the Seismic Response to a Set of Parallel Discrete Fractures Yang Zhang 1, Xander Campman 1, Samantha Grandi 1, Shihong Chi 1, M. Nafi Toksöz 1, Mark E. Willis 1, Daniel R. Burns
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
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
Interpretation of PP and PS seismic data from the Mackenzie Delta, N.W.T. Carlos E. Nieto* CREWES. The University of Calgary email@example.com and Robert R. Stewart CREWES. The University of Calgary
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
Science in China Ser. G Physics, Mechanics & Astronomy 2004 Vol.47 No.3 265 276 265 A new method for multi-exponential inversion of NMR relaxation measurements WANG Zhongdong 1, 2, XIAO Lizhi 1 & LIU Tangyan
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)
Th-01-08 Guided Waves - Inversion and Attenuation D. Boiero* (WesternGeco), C. Strobbia (WesternGeco), L. Velasco (WesternGeco) & P. Vermeer (WesternGeco) SUMMARY Guided waves contain significant information
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
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.
Pressure Prediction and Hazard Avoidance through Improved Seismic Imaging 12121-6002-02 Bill Barkhouse SEG Advanced Modeling Corporation (SEAM) RPSEA / SPE-GCS Ultra-Deepwater Technology Conference September
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,
Chapter 6 Conclusions 6.1 Conclusions and perspectives In this thesis an approach is presented for the in-situ characterization of rocks in terms of the distribution of hydraulic parameters (called SBRC
Well Evaluation An Integrated Well-Evaluation Process for Shale Gas Reservoirs APPLICATION Gas shale analysis for unconventional wells BENEFITS Optimized stimulation by thorough study and planning Increased
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.
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
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
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,
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
Available online at www.sciencedirect.com ScienceDirect Energy Procedia 114 (2017 ) 2772 2780 13th International Conference on Greenhouse Gas Control Technologies, GHGT-13, 14-18 November 2016, Lausanne,
Copyright by Barry Michael Borgman 2016 The Thesis Committee for Barry Michael Borgman Certifies that this is the approved version of the following thesis: Characterization of the Cana-Woodford Shale using
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
10-810: Advanced Algorithms and Models for Computational Biology Optimal leaf ordering and classification Hierarchical clustering As we mentioned, its one of the most popular methods for clustering gene
Integrating Geomechanics and Reservoir Characterization Examples from Canadian Shale Plays AAPG Geosciences Technology Workshops Geomechanics and Reservoir Characterization of Shale and Carbonates July
Frequently Asked Questions 1 Is THERMOVISION TOMOGRAPHY technology independent, or it requires previous seismic or other data? THERMOVISION TOMOGRAPHY technology is a stand-alone technology. Independent
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,
GEOPHYSICS, VOL. 58, NO. 8 (AUGUST 1993); P. 1183-1194, 12 FIGS. Determination of the signature of a dynamite source using source scaling, Part 2: Experiment Anton Ziolkowski* and Karel Bokhorst ABSTRACT
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,
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
UCRL-ID- 12873 1 Modified Kriging: Evaluation, Modification, Recommendations and C.A. Schultz S.C. Myers September 1997 This is an informal report intended primarily forinternal or limited external distribution.
GEOPHYSICS, VOL. 75, NO. 6 NOVEMBER-DECEMBER 21 ; P. WB265 WB274, 11 FIGS., 1 TABLE. 1.119/1.347539 The effect of interpolation on imaging and : A Viking case study Lee Hunt 1, Jon Downton 2, Scott Reynolds
34 Chapter 3 Isotropic scattering and crustal imaging with wide-angle prestack migration 3.1 Introduction Earthquake seismologists may be reluctant to believe that one can successfully image crustal faults
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,
PROCEEDINGS, Thirty-Eighth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 11-13, 2013 SGP-TR-198 FULL MOMENT TENSOR ANALYSIS USING FIRST MOTION DATA AT
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.
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,
Determination of Duvernay Formation Reservoir Properties through Probabilistic Petrophysical Analysis calibrated to Core Studies. Nasir Rahim, Neil Watson Canadian Discovery Ltd. Summary The petrophysical
Exploration Well Failures from the Moray Firth & Central North Sea (UK) 21 st Century Exploration Road Map Project Christian Mathieu 1. Project Objectives, Time-Line and Status 2. A few statistics 3. Reasons
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
Receiver Function Inversion Advanced Studies Institute on Seismological Research Kuwait City, Kuwait - January 19-22, 2013 Jordi Julià Universidade Federal do Rio Grande do Norte, Brasil Outline Introduction
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
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
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
GEOPHYSICS, VOL. 60, NO. 5 (SEPTEMBER-OCTOBER 1995); P. 1409-1425, 11 FIGS. Body-wave radiation patterns and AVO in transversely isotropic media llya Tsvankin* ABSTRACT The angular dependence of reflection
Bayesian Modeling and Classification of Neural Signals Michael S. Lewicki Computation and Neural Systems Program California Institute of Technology 216-76 Pasadena, CA 91125 lewickiocns.caltech.edu Abstract
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
umerical Modeling for Different Types of Fractures Xiaoqin Cui* CREWES Department of Geoscience University of Calgary Canada firstname.lastname@example.org and Laurence R. Lines Edward S. Krebes Department of Geoscience
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
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
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
URTeC: 2444366 Using Depletion-Zone Microseismicity to Understand Producing Volumes Jonathan P. McKenna*, Michael H. Grealy, Michael S. Blaz and Nathan M. Toohey, MicroSeismic, Inc. Copyright 2016, Unconventional
Pn Q UNDER TIBET AND TIENSHAN WITH PRACTICAL AND SCIENTIFIC IMPLICATIONS Jiakang Xie Lamont-Doherty Earth Observatory, Columbia University Sponsored by Defense Threat Reduction Agency Contract No. DTRA01-00-C-0048
Anisotropic velocity changes in seismic time-lapse data Jorg Herwanger*, Ed Palmer (WesternGeco), Christian Rau Schiøtt (Hess) Summary Reservoir production causes subsurface deformations and changes in
Chapter 2 Theory - SBRC fundamentals The approach proposed to estimate hydraulic properties of rocks using microseismicity is called Seismicity Based Reservoir Characterization (SBRC). It uses a spatio-temporal
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
Multi-offset vertical seismic profiling: an experiment to assess petrophysical-scale parameters at the JAPEX/JNOC/GSC et al. Mallik 5L-38 gas hydrate production research well B. Milkereit 1, E. Adam 2,
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
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
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