A tour of optimisation methods for facies estimation in AVO seismic inversion using Markov Random Fields
|
|
- Brook McKinney
- 5 years ago
- Views:
Transcription
1 A tour of optimisation methods for facies estimation in AVO seismic inversion using Markov Random Fields James Gunning, CSIRO ESRE Michel Kemper, Denis Saussus, Ikon Science Adrian Pelham, Elaine Fitzgerald, Tullow Oil January 16, 2013 Abstract We overview some joint lithology-fluid-class/rock-properties inversion techniques for AVO/quant- -itative-interpretation work from true-amplitude imaged seismic data, using a Bayesian hierarchical model. We use a Markov random-field (MRF) model for facies labels, and a facies-conditional Normal model for elastic rock properties. The MRF forces spatial smoothness in categories, and forbids non-hydrodynamical fluid-class transitions. The rock properties model captures loading effects and mutual correlation between elastic variables within a facies. We study the optimisation problem rather than sampling. The optimisation is a mixed integer-continuous problem, provably NP-hard, but susceptible to good heuristics in certain regimes we discuss. A natural approach to the problem is the EM algorithm, which provides fast local solutions using loopy-belief propagation for the E-step and conjugate-gradient least-squares for the M-step. This yields soft/marginalised estimates of facies labels. Industry-standard Bayesian classifications are related to single cycles of this EM algorithm, but without spatial coupling. A comparable method to the EM algorithm is relaxed joint inversion with continuous facies indicator variables confined to a simplex. The absence of middle frequencies and the multimodality of the objective make fast gradient methods, including EM and relaxation, vulnerable to local minima. Two globalising approaches, deterministic annealing, and analytical reduction to integer-optimisation, are discussed. Introduction Inversion of reflection seismic data for rock properties is a notoriously non unique problem. All acquisition frameworks have bandwidth and physics sensitivity limitations, which imply the data has insufficient information to allow unambiguous inference of absolute rock properties. Standard imaging/migration workflows acknowledge this problem by generating bandlimited images of reflection coefficient only, and, where seismic amplitudes are expected to be useful in reservoir characterisation, aim to produce true amplitude volumes where the imaged amplitudes are well modelled by specified incidence reflectivity equations. This passes much of the non uniqueness problem on to downstream workflows, where quantitative interpretation or AVO teams can work on a simplified inversion problem with migrated seismic as 1D bandlimited data, and the 3 dimensional aspects of the wave propagation largely removed. In AVO inversions, additional information brought in to address non uniqueness often consists of low frequency trends for rock properties, extrapolated from well data into space. These extrapolated trends are usually pooled/weighted across rock types, and their contribution to AVO inversion is usually a convex stabilising term added to an AVO amplitude mismatch objective function, taken in either a Bayesian or regularising sense. Wavelets and noise levels are assumed known from well ties. The overall AVO inversion problem is then convex and efficiently solvable using standard gradient based techniques, characteristics which we think very attractive. A significant problem with this approach is that, far from
2 wells, the background trend information yields only very low frequency drift. After inversion, the middle frequencies get lost, we do not recover blocky looking models, and direct facies classification becomes unreliable. (Broadband seismic is obviously a measurement approach to this issue). Most convex stabilising terms, especially those of l 2 form, are hostile to blocky model restorations, and the latter is what we commonly like to see when the subsurface geology is expected to have distinct and extensive facies. An approach to this problem is to drop the convex, facies pooled stabilising apparatus, replacing this with facies dependent rock property distributions in a Bayesian framework, and attempting joint inversion for facies types and rock properties. The inversion problem becomes a nonconvex, mixed continuous/discrete one, and thus computationally very challenging. The first difficulty is in formulating a discrete prior distribution for the facies labels in space. Many techniques for stochastic facies simulation, such as multipoint or indicator models, generate this distribution implicitly, and the consequent lack of closed forms make these choices problematic in either an optimisation or MCMC framework. In grid free or off lattice approaches, stochastic geometries represented by marked point models do have closed forms, and can be very useful in appraisal/production phases, since the facies sequence and positions may be nearly known. For models on a lattice, Markov random fields (MRFs) are a workable choice, and have the advantage that spatial interactions operate only locally. Recent treatments of joint facies and elastic inversion using MRFs are Larsen et al. (2006), Ulvmoen and Omre (2010). These approaches hinge on linearisations about a particular agnostic background model, from which the continuous rock properties can be approximately marginalised out of the problem. This leaves an inversion problem for the posterior distribution of the discrete facies labels only, and the authors demonstrate how to sample this using MCMC. In the example problems they discuss, the facies uncertainty is noticeably small in the posterior samples. We think this implies that optimisation for the maximum aposteriori (MAP) state would therefore be very useful. In this paper we focus on the problem of joint MAP estimation of facies labels and rock properties, using a Bayesian model with MRFs for the facies labels, and conditional normal distributions for the rock properties which capture rock physics effects and loading. The full optimisation problem is without doubt NP hard, so we present several approaches spanning a range of rigour and CPU cost. We sketch some local approaches based on the expectation maximisation (EM) algorithm, including graph cutting and loopy belief propagation. Closely related to the EM algorithm are gradient methods based a relaxation of the discrete labelling space to indicator variables. A further layer of globalisation can be attempted by embedding the algorithms above inside a deterministic annealing loop (Rangarajan, 2000). Finally, we discuss some more exact methods based on integer quadratic programming over moving templates. Theory In rough notation, the model to be estimated is the model vector of elastic properties m = {v p, v s, ρ} together with facies labels F (F i [1...N f ]) at all the lattice sample locations i = 1...N p in the volume (x, y, t). The data y are angle stacks of true amplitude imaged seismic reflectivities. The prior distribution of rock properties is a facies conditional joint distribution, constructed from rock physics regression equations developed from regional well log analysis, which we write as P(m F). The facies labels are assigned a discrete prior distribution P(F) over the volume. Given the forward model f(m), and a suitable likelihood of the data given the model L(y m) = L(y f(m), the joint posterior distribution of the model m and F can be written P(m, F y) L(y m)p(m F)P(F). (1) The MRF distribution for the facies labels, set on a 3D stratigraphic grid with pairwise neighbours denoted ( j i), is of form log(p(f)) = i (e Fi + j i 2β (D i j) I(F i F j )). Parameters e Fi are for
3 proportions and β (D i j) for connectivity of labels in direction D i j, presumed known in this note. This allows control of facies transitions in all directions. The conditional prior for the rock properties P(m F) is a product of independent Normals of form P(m i ={v p,i,v s,i,ρ i }) N( m i (F i ),C p,i (F i )). For simplicity, we focus on the simplest case of a fully linearised likelihood (Fatti model reflectivity) with N(0,C d ) Gaussian noise model, where the likelihood is L(y m)=n(y Xm,C d ), and X a multi stack block matrix representing the matrix product of vertical differentiation, Fatti AVO coefficients, and wavelet convolution. Finding the joint MAP estimate of facies and properties then amounts to minimising the -ve log posterior of (1), which looks like E(m,F) = (y Xm) T Cd 1 +log C p (F) + i (y Xm)+(m m(f)) T C 1 (F)(m m(f)) (e Fi + A hierarchy of increasingly difficult optimisation strategies follow j i p 2β (D i j) I(F i F j )). (2) 1. One shot vanilla Bayesian classification (closely related to various industry standards). Amounts to fixing m(f) at some background/hybrid mean, (say m ), minimising (2) over m (F bits neglected), then fixing the resulting m, & minimising over F (with all β = 0) for most likely facies. Many approaches depend on the accuracy of marginal distributions estimated/approximated from this inversion, e.g. Ulvmoen and Omre (2010), and some compensations can be made (Buland et al., 2008). Usually sensitive to the initial hybrid. 2. Alternating (EM like) methods. (i) K means/graph cut: minimise (2) alternately in spaces of m (large scale least squares via conjugate gradient) andf (using graph cutting, possible since the objective is submodular (Boykov et al., 2001)). The MRF cleans up short spatial scale noise in the model considerably. CPU costs are a small multiple of classical least squares inversion. Usually rather sensitive to choice of initial hybrid. (ii) Full EM algorithm: here the E step estimates soft facies memberships by loopy belief propagation with current model m fixed. The M step to update m is least squares but with a prior formed by a current membership weighted average of the facies prior means and covariances. With membership step limiting this helps mitigate the sensitivity to initial guesses appreciably. CPU costs are O(10-100) times costs of classical inversion. 3. Relaxed membership joint optimisation. Here, the discrete variables F are relaxed to continuous indicators γ il [0,1] (labels l, voxel i), and the MRF is recast to a quadratic form in the γs yielding the same energy transitions. The optimisation is constrained to the simplex l γ il = 1. The objective is then positive definite in m, but negative definite in the γs, which is a quadratic programming problem known to be NP hard. Local optima can be found efficiently using projected gradient descent, and again, step limiting in the membership (γ) space has proved helpful. This method is about comparably effective (and expensive) to the full EM algorithm, but bypasses belief propagation. The initial guess sensitivity seems similar. 4. Deterministic annealing (DA) variations. In either of (2) or (3) from above, the nonconvexity of the objective makes for initial guess dependence. DA methods wrap the optimisation in an outer cooling loop, with a temperature parameter gradually evolving the objective from some initially convex hybrid to its final multimodal form. The MAP point is tracked from its initial unique minima as the cooling proceeds. Optimisation is gradient based and efficient for each temperature. Our DA peeling scheme consists in pulling the facies normal distributions slowly away from a common starting hybrid which is surely convex. 5. Moving template integer quadratic programming (IQP). Facies clusters are approximated, sharing a common (abundance weighted average) rock properties covariance. The m optimisation of (2)
4 is done analytically, yielding an IQP problem for the facies labels. For realistic signal to noise ratios, the IQP matrix is not even close to submodular, so efficient graph cutting techniques are inapplicable. For smallish moving spatial templates, this can be converted to a binary linear programming problem and fed to a commodity IP solver. This has no startup dependence (template optima are global), but computational demands are very fierce. Two Examples Problems with 2 or 3 facies clusters which align well in P impedance space, in a dominant lithology background, yield to most of the faster methods above, with some tuning required. Example 1, NW Australian shelf. Here, Fig. 1 shows facies inversions for a turbidite field with strong Gassmann oil response. Peak S/N ratio is perhaps 5 here. The insets d,e,f, (detail at a trace) shows EM can be mildly robust to starting guesses cf a global optima. Example 2, North Sea: when the facies clusters are not t(s) t t (a) near-stack seismic + well data 2400 (b) facies, EM algorithm 2400 (c) facies, IQP algorithm (d) trace properties, EM run1 (e) trace properties, EM run2 (f) trace properties, IQP run Figure 1: (a-c) Seismic and MAP facies estimates from EM and IQP algorithms. Single trace detail: (d,e,f) show facies clusters (brown=shale, blue=brine, orange=oil), P impedance contours (thin red), (d,e) hollow squares=1st sweep inversion from 2 different facies-proportioned starting hybrids, filled squares=final model (f)=direct MAP properties from IQP optimum aligned, images are noisier etc, inferences can be much more fragile to the starting/background hybrid model and/or convexity approximations in the faster methods; see Fig 2. Here, the S/N ratio is much worse, and classification is much more reliant on the model and optimisation algorithm. Cheap methods can work very well with excellent S/N. Our Monte Carlo studies of noise effects show that facies misclassification rates, even with perfect optimization, grow sharply for S/N< 4, so the hard regime, where expensive optimization is warranted, may not be large. Conclusions Joint MAP inversion for rock properties and facies is a hard problem. Fast, local optimisation methods based on the EM and related algorithms can be remarkably effective, but good starting guesses and S/N
5 (a) Basic EM (b) EM within DA (c) IQP Figure 2: More difficult, 4 stack, noisy data, non aligned facies; behaviour at a single trace. (a) Basic EM algorithm (b) DA algorithm with inner EM optimisation (c) IQP. The DA algorithm, of controllable CPU budget (perhaps 10-fold the cost of basic EM ), is not too far from very expensive IQP. are needed. More ambiguous rock physics mixtures will require expensive globalising methods, and here, DA looks an affordable alternative to fearsome integer methods. References Boykov, Y., Veksler, O. and Zabih, R. [2001] Fast approximate energy minimization via graph cuts. IEEE Trans. Pattern Anal. Mach. Intell., 23, , ISSN , doi: / Buland, A., Kolbjornsen, O., Hauge, R., Skjaeveland, O. and Duffaut, K. [2008] Bayesian lithology and fluid prediction from seismic prestack data. Geophysics, 73, 13 21, doi: / Larsen, A.L., Ulvmoen, M., Omre, H. and Buland, A. [2006] Bayesian lithology/fluid prediction and simulation on the basis of a Markov-chain prior model. Geophysics, 71(5), R69 R78. Rangarajan, A. [2000] Self annealing and self annihilation: Unifying deterministic annealing and relaxation labeling. In Pattern Recognition, Ulvmoen, M. and Omre, H. [2010] Improved resolution in Bayesian lithology/fluid inversion from prestack seismic data and well observations: Part 1 methodology. Geophysics, 75(2), R21 R35.
23855 Rock Physics Constraints on Seismic Inversion
23855 Rock Physics Constraints on Seismic Inversion M. Sams* (Ikon Science Ltd) & D. Saussus (Ikon Science) SUMMARY Seismic data are bandlimited, offset limited and noisy. Consequently interpretation of
More informationDownloaded 10/02/18 to Redistribution subject to SEG license or copyright; see Terms of Use at
Multi-scenario, multi-realization seismic inversion for probabilistic seismic reservoir characterization Kester Waters* and Michael Kemper, Ikon Science Ltd. Summary We propose a two tiered inversion strategy
More informationBayesian lithology/fluid inversion comparison of two algorithms
Comput Geosci () 14:357 367 DOI.07/s596-009-9155-9 REVIEW PAPER Bayesian lithology/fluid inversion comparison of two algorithms Marit Ulvmoen Hugo Hammer Received: 2 April 09 / Accepted: 17 August 09 /
More informationBayesian Lithology-Fluid Prediction and Simulation based. on a Markov Chain Prior Model
Bayesian Lithology-Fluid Prediction and Simulation based on a Markov Chain Prior Model Anne Louise Larsen Formerly Norwegian University of Science and Technology, N-7491 Trondheim, Norway; presently Schlumberger
More informationStatistical Rock Physics
Statistical - Introduction Book review 3.1-3.3 Min Sun March. 13, 2009 Outline. What is Statistical. Why we need Statistical. How Statistical works Statistical Rock physics Information theory Statistics
More informationWe LHR3 06 Detecting Production Effects and By-passed Pay from 3D Seismic Data Using a Facies Based Bayesian Seismic Inversion
We LHR3 06 Detecting Production Effects and By-passed Pay from 3D Seismic Data Using a Facies Based Bayesian Seismic Inversion K.D. Waters* (Ikon Science Ltd), A.V. Somoza (Ikon Science Ltd), G. Byerley
More informationReservoir connectivity uncertainty from stochastic seismic inversion Rémi Moyen* and Philippe M. Doyen (CGGVeritas)
Rémi Moyen* and Philippe M. Doyen (CGGVeritas) Summary Static reservoir connectivity analysis is sometimes based on 3D facies or geobody models defined by combining well data and inverted seismic impedances.
More informationQUANTITATIVE 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 informationSeismic reservoir characterization in offshore Nile Delta.
Seismic reservoir characterization in offshore Nile Delta. Part II: Probabilistic petrophysical-seismic inversion M. Aleardi 1, F. Ciabarri 2, B. Garcea 2, A. Mazzotti 1 1 Earth Sciences Department, University
More informationNet-to-gross from Seismic P and S Impedances: Estimation and Uncertainty Analysis using Bayesian Statistics
Net-to-gross from Seismic P and S Impedances: Estimation and Uncertainty Analysis using Bayesian Statistics Summary Madhumita Sengupta*, Ran Bachrach, Niranjan Banik, esterngeco. Net-to-gross (N/G ) is
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 informationStochastic vs Deterministic Pre-stack Inversion Methods. Brian Russell
Stochastic vs Deterministic Pre-stack Inversion Methods Brian Russell Introduction Seismic reservoir analysis techniques utilize the fact that seismic amplitudes contain information about the geological
More 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 informationThe GIG consortium Geophysical Inversion to Geology Per Røe, Ragnar Hauge, Petter Abrahamsen FORCE, Stavanger
www.nr.no The GIG consortium Geophysical Inversion to Geology Per Røe, Ragnar Hauge, Petter Abrahamsen FORCE, Stavanger 17. November 2016 Consortium goals Better estimation of reservoir parameters from
More informationPattern Recognition and Machine Learning
Christopher M. Bishop Pattern Recognition and Machine Learning ÖSpri inger Contents Preface Mathematical notation Contents vii xi xiii 1 Introduction 1 1.1 Example: Polynomial Curve Fitting 4 1.2 Probability
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 informationNORGES TEKNISK-NATURVITENSKAPELIGE UNIVERSITET. Directional Metropolis Hastings updates for posteriors with nonlinear likelihoods
NORGES TEKNISK-NATURVITENSKAPELIGE UNIVERSITET Directional Metropolis Hastings updates for posteriors with nonlinear likelihoods by Håkon Tjelmeland and Jo Eidsvik PREPRINT STATISTICS NO. 5/2004 NORWEGIAN
More informationMultiple Scenario Inversion of Reflection Seismic Prestack Data
Downloaded from orbit.dtu.dk on: Nov 28, 2018 Multiple Scenario Inversion of Reflection Seismic Prestack Data Hansen, Thomas Mejer; Cordua, Knud Skou; Mosegaard, Klaus Publication date: 2013 Document Version
More informationQuantitative 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 informationSensitivity Analysis of Pre stack Seismic Inversion on Facies Classification using Statistical Rock Physics
Sensitivity Analysis of Pre stack Seismic Inversion on Facies Classification using Statistical Rock Physics Peipei Li 1 and Tapan Mukerji 1,2 1 Department of Energy Resources Engineering 2 Department of
More informationAn Integrated Workflow for Seismic Data Conditioning and Modern Prestack Inversion Applied to the Odin Field. P.E.Harris, O.I.Frette, W.T.
An Integrated Workflow for Seismic Data Conditioning and Modern Prestack Inversion Applied to the Odin Field P.E.Harris, O.I.Frette, W.T.Shea Talk Outline Introduction Motivation Introducing Pcube+ Gather
More informationGaussian Process Approximations of Stochastic Differential Equations
Gaussian Process Approximations of Stochastic Differential Equations Cédric Archambeau Centre for Computational Statistics and Machine Learning University College London c.archambeau@cs.ucl.ac.uk CSML
More informationHampsonRussell. A comprehensive suite of reservoir characterization tools. cgg.com/geosoftware
HampsonRussell A comprehensive suite of reservoir characterization tools cgg.com/geosoftware HampsonRussell Software World-class geophysical interpretation HampsonRussell Software is a comprehensive suite
More informationSeismic inversion for reservoir properties combining statistical rock physics and geostatistics: A review
GEOPHYSICS, VOL. 75, NO. 5 SEPTEMBER-OCTOBER 2010 ; P. 75A165 75A176, 8 FIGS. 10.1190/1.3478209 Seismic inversion for reservoir properties combining statistical rock physics and geostatistics: A review
More informationEmpirical comparison of two Bayesian lithology fluid prediction algorithms
Empirical comparison of two Bayesian lithology fluid prediction algorithms Hugo Hammer and Marit Ulvmoen Department of Mathematical Sciences Norwegian University of Science and Technology Trondheim, Norway
More informationUncertainty quantification for Wavefield Reconstruction Inversion
Uncertainty quantification for Wavefield Reconstruction Inversion Zhilong Fang *, Chia Ying Lee, Curt Da Silva *, Felix J. Herrmann *, and Rachel Kuske * Seismic Laboratory for Imaging and Modeling (SLIM),
More informationQuantitative interpretation using inverse rock-physics modeling on AVO data
Quantitative interpretation using inverse rock-physics modeling on AVO data Erling Hugo Jensen 1, Tor Arne Johansen 2, 3, 4, Per Avseth 5, 6, and Kenneth Bredesen 2, 7 Downloaded 08/16/16 to 129.177.32.62.
More informationIntroduction to Optimization
Introduction to Optimization Blackbox Optimization Marc Toussaint U Stuttgart Blackbox Optimization The term is not really well defined I use it to express that only f(x) can be evaluated f(x) or 2 f(x)
More informationApproximate Inference Part 1 of 2
Approximate Inference Part 1 of 2 Tom Minka Microsoft Research, Cambridge, UK Machine Learning Summer School 2009 http://mlg.eng.cam.ac.uk/mlss09/ Bayesian paradigm Consistent use of probability theory
More informationApproximate Inference Part 1 of 2
Approximate Inference Part 1 of 2 Tom Minka Microsoft Research, Cambridge, UK Machine Learning Summer School 2009 http://mlg.eng.cam.ac.uk/mlss09/ 1 Bayesian paradigm Consistent use of probability theory
More informationA Petroleum Geologist's Guide to Seismic Reflection
A Petroleum Geologist's Guide to Seismic Reflection William Ashcroft WILEY-BLACKWELL A John Wiley & Sons, Ltd., Publication Contents Preface Acknowledgements xi xiii Part I Basic topics and 2D interpretation
More informationComparative Study of AVO attributes for Reservoir Facies Discrimination and Porosity Prediction
5th Conference & Exposition on Petroleum Geophysics, Hyderabad-004, India PP 498-50 Comparative Study of AVO attributes for Reservoir Facies Discrimination and Porosity Prediction Y. Hanumantha Rao & A.K.
More informationMultiple horizons mapping: A better approach for maximizing the value of seismic data
Multiple horizons mapping: A better approach for maximizing the value of seismic data Das Ujjal Kumar *, SG(S) ONGC Ltd., New Delhi, Deputed in Ministry of Petroleum and Natural Gas, Govt. of India Email:
More informationAFI (AVO Fluid Inversion)
AFI (AVO Fluid Inversion) Uncertainty in AVO: How can we measure it? Dan Hampson, Brian Russell Hampson-Russell Software, Calgary Last Updated: April 2005 Authors: Dan Hampson, Brian Russell 1 Overview
More informationConsistent Downscaling of Seismic Inversions to Cornerpoint Flow Models SPE
Consistent Downscaling of Seismic Inversions to Cornerpoint Flow Models SPE 103268 Subhash Kalla LSU Christopher D. White LSU James S. Gunning CSIRO Michael E. Glinsky BHP-Billiton Contents Method overview
More informationP191 Bayesian Linearized AVAZ Inversion in HTI Fractured Media
P9 Bayesian Linearized AAZ Inversion in HI Fractured Media L. Zhao* (University of Houston), J. Geng (ongji University), D. Han (University of Houston) & M. Nasser (Maersk Oil Houston Inc.) SUMMARY A new
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 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 informationBayesian Networks: Construction, Inference, Learning and Causal Interpretation. Volker Tresp Summer 2016
Bayesian Networks: Construction, Inference, Learning and Causal Interpretation Volker Tresp Summer 2016 1 Introduction So far we were mostly concerned with supervised learning: we predicted one or several
More informationSTA 4273H: Sta-s-cal Machine Learning
STA 4273H: Sta-s-cal Machine Learning Russ Salakhutdinov Department of Computer Science! Department of Statistical Sciences! rsalakhu@cs.toronto.edu! h0p://www.cs.utoronto.ca/~rsalakhu/ Lecture 2 In our
More informationSTA414/2104 Statistical Methods for Machine Learning II
STA414/2104 Statistical Methods for Machine Learning II Murat A. Erdogdu & David Duvenaud Department of Computer Science Department of Statistical Sciences Lecture 3 Slide credits: Russ Salakhutdinov Announcements
More informationDownloaded 09/09/15 to Redistribution subject to SEG license or copyright; see Terms of Use at
Reservoir properties estimation from marine broadband seismic without a-priori well information: A powerful de-risking workflow Cyrille Reiser*, Matt Whaley and Tim Bird, PGS Reservoir Limited Summary
More informationMessage-Passing Algorithms for GMRFs and Non-Linear Optimization
Message-Passing Algorithms for GMRFs and Non-Linear Optimization Jason Johnson Joint Work with Dmitry Malioutov, Venkat Chandrasekaran and Alan Willsky Stochastic Systems Group, MIT NIPS Workshop: Approximate
More informationDevelopment of Stochastic Artificial Neural Networks for Hydrological Prediction
Development of Stochastic Artificial Neural Networks for Hydrological Prediction G. B. Kingston, M. F. Lambert and H. R. Maier Centre for Applied Modelling in Water Engineering, School of Civil and Environmental
More informationUncertainty analysis for the integration of seismic and CSEM data Myoung Jae Kwon & Roel Snieder, Center for Wave Phenomena, Colorado School of Mines
Myoung Jae Kwon & Roel Snieder, Center for Wave Phenomena, Colorado School of Mines Summary Geophysical inverse problems consist of three stages: the forward problem, optimization, and appraisal. We study
More informationReservoir properties inversion from AVO attributes
Reservoir properties inversion from AVO attributes Xin-gang Chi* and De-hua Han, University of Houston Summary A new rock physics model based inversion method is put forward where the shaly-sand mixture
More informationNew Frontier Advanced Multiclient Data Offshore Uruguay. Advanced data interpretation to empower your decision making in the upcoming bid round
New Frontier Advanced Multiclient Data Offshore Uruguay Advanced data interpretation to empower your decision making in the upcoming bid round Multiclient data interpretation provides key deliverables
More informationDelivery: an open-source Bayesian seismic inversion tool. James Gunning, CSIRO Petroleum Michael Glinsky,, BHP Billiton
Delivery: an open-source Bayesian seismic inversion tool James Gunning, CSIRO Petroleum Michael Glinsky,, BHP Billiton Application areas Development, appraisal and exploration situations Requisite information
More informationA Generative Perspective on MRFs in Low-Level Vision Supplemental Material
A Generative Perspective on MRFs in Low-Level Vision Supplemental Material Uwe Schmidt Qi Gao Stefan Roth Department of Computer Science, TU Darmstadt 1. Derivations 1.1. Sampling the Prior We first rewrite
More informationLecture 9: PGM Learning
13 Oct 2014 Intro. to Stats. Machine Learning COMP SCI 4401/7401 Table of Contents I Learning parameters in MRFs 1 Learning parameters in MRFs Inference and Learning Given parameters (of potentials) and
More informationIntegrating rock physics modeling, prestack inversion and Bayesian classification. Brian Russell
Integrating rock physics modeling, prestack inversion and Bayesian classification Brian Russell Introduction Today, most geoscientists have an array of tools available to perform seismic reservoir characterization.
More informationProbabilistic seismic inversion using pseudo-wells
Seismic Rock Physics Seminar Probabilistic seismic inversion using pseudo-wells Patrick Connolly*, PCA Ltd Patrick Connolly Associates Ltd. geophysics for integration Outline ODiSI: probabilistic inversion
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 Multi-scale geological analysis GEO4
More informationSTA 4273H: Statistical Machine Learning
STA 4273H: Statistical Machine Learning Russ Salakhutdinov Department of Statistics! rsalakhu@utstat.toronto.edu! http://www.utstat.utoronto.ca/~rsalakhu/ Sidney Smith Hall, Room 6002 Lecture 3 Linear
More informationProteomics and Variable Selection
Proteomics and Variable Selection p. 1/55 Proteomics and Variable Selection Alex Lewin With thanks to Paul Kirk for some graphs Department of Epidemiology and Biostatistics, School of Public Health, Imperial
More informationCheng Soon Ong & Christian Walder. Canberra February June 2018
Cheng Soon Ong & Christian Walder Research Group and College of Engineering and Computer Science Canberra February June 218 Outlines Overview Introduction Linear Algebra Probability Linear Regression 1
More informationDownloaded 09/16/16 to Redistribution subject to SEG license or copyright; see Terms of Use at
Ehsan Zabihi Naeini*, Ikon Science & Russell Exley, Summit Exploration & Production Ltd Summary Quantitative interpretation (QI) is an important part of successful Central North Sea exploration, appraisal
More informationWe LHR3 04 Realistic Uncertainty Quantification in Geostatistical Seismic Reservoir Characterization
We LHR3 04 Realistic Uncertainty Quantification in Geostatistical Seismic Reservoir Characterization A. Moradi Tehrani* (CGG), A. Stallone (Roma Tre University), R. Bornard (CGG) & S. Boudon (CGG) SUMMARY
More informationBayesian Regression Linear and Logistic Regression
When we want more than point estimates Bayesian Regression Linear and Logistic Regression Nicole Beckage Ordinary Least Squares Regression and Lasso Regression return only point estimates But what if we
More informationCSC 2541: Bayesian Methods for Machine Learning
CSC 2541: Bayesian Methods for Machine Learning Radford M. Neal, University of Toronto, 2011 Lecture 10 Alternatives to Monte Carlo Computation Since about 1990, Markov chain Monte Carlo has been the dominant
More informationLecture 16: Small Sample Size Problems (Covariance Estimation) Many thanks to Carlos Thomaz who authored the original version of these slides
Lecture 16: Small Sample Size Problems (Covariance Estimation) Many thanks to Carlos Thomaz who authored the original version of these slides Intelligent Data Analysis and Probabilistic Inference Lecture
More informationNonparametric Bayesian Methods (Gaussian Processes)
[70240413 Statistical Machine Learning, Spring, 2015] Nonparametric Bayesian Methods (Gaussian Processes) Jun Zhu dcszj@mail.tsinghua.edu.cn http://bigml.cs.tsinghua.edu.cn/~jun State Key Lab of Intelligent
More informationA6523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring
A6523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring 2015 http://www.astro.cornell.edu/~cordes/a6523 Lecture 23:! Nonlinear least squares!! Notes Modeling2015.pdf on course
More informationHeeyoul (Henry) Choi. Dept. of Computer Science Texas A&M University
Heeyoul (Henry) Choi Dept. of Computer Science Texas A&M University hchoi@cs.tamu.edu Introduction Speaker Adaptation Eigenvoice Comparison with others MAP, MLLR, EMAP, RMP, CAT, RSW Experiments Future
More information3D geologic modelling of channellized reservoirs: applications in seismic attribute facies classification
first break volume 23, December 2005 technology feature 3D geologic modelling of channellized reservoirs: applications in seismic attribute facies classification Renjun Wen, * president and CEO, Geomodeling
More informationElements of 3D Seismology Second Edition
Elements of 3D Seismology Second Edition Copyright c 1993-2003 All rights reserved Christopher L. Liner Department of Geosciences University of Tulsa August 14, 2003 For David and Samantha And to the memory
More informationPart 1: Expectation Propagation
Chalmers Machine Learning Summer School Approximate message passing and biomedicine Part 1: Expectation Propagation Tom Heskes Machine Learning Group, Institute for Computing and Information Sciences Radboud
More informationBayesian Networks: Construction, Inference, Learning and Causal Interpretation. Volker Tresp Summer 2014
Bayesian Networks: Construction, Inference, Learning and Causal Interpretation Volker Tresp Summer 2014 1 Introduction So far we were mostly concerned with supervised learning: we predicted one or several
More informationBayesian Machine Learning
Bayesian Machine Learning Andrew Gordon Wilson ORIE 6741 Lecture 2: Bayesian Basics https://people.orie.cornell.edu/andrew/orie6741 Cornell University August 25, 2016 1 / 17 Canonical Machine Learning
More informationInference in Bayesian Networks
Andrea Passerini passerini@disi.unitn.it Machine Learning Inference in graphical models Description Assume we have evidence e on the state of a subset of variables E in the model (i.e. Bayesian Network)
More informationCheng Soon Ong & Christian Walder. Canberra February June 2018
Cheng Soon Ong & Christian Walder Research Group and College of Engineering and Computer Science Canberra February June 2018 Outlines Overview Introduction Linear Algebra Probability Linear Regression
More informationMachine Learning Techniques for Computer Vision
Machine Learning Techniques for Computer Vision Part 2: Unsupervised Learning Microsoft Research Cambridge x 3 1 0.5 0.2 0 0.5 0.3 0 0.5 1 ECCV 2004, Prague x 2 x 1 Overview of Part 2 Mixture models EM
More informationBayesian Learning in Undirected Graphical Models
Bayesian Learning in Undirected Graphical Models Zoubin Ghahramani Gatsby Computational Neuroscience Unit University College London, UK http://www.gatsby.ucl.ac.uk/ Work with: Iain Murray and Hyun-Chul
More informationPhysics 403. Segev BenZvi. Numerical Methods, Maximum Likelihood, and Least Squares. Department of Physics and Astronomy University of Rochester
Physics 403 Numerical Methods, Maximum Likelihood, and Least Squares Segev BenZvi Department of Physics and Astronomy University of Rochester Table of Contents 1 Review of Last Class Quadratic Approximation
More informationLecture 6: Graphical Models: Learning
Lecture 6: Graphical Models: Learning 4F13: Machine Learning Zoubin Ghahramani and Carl Edward Rasmussen Department of Engineering, University of Cambridge February 3rd, 2010 Ghahramani & Rasmussen (CUED)
More informationC002 Petrophysical Seismic Inversion over an Offshore Carbonate Field
C002 Petrophysical Seismic Inversion over an Offshore Carbonate Field T. Coleou* (CGGVeritas), F. Allo (CGGVeritas), O. Colnard (CGGVeritas), I. Machecler (CGGVeritas), L. Dillon (Petrobras), G. Schwedersky
More informationBayesian Methods for Machine Learning
Bayesian Methods for Machine Learning CS 584: Big Data Analytics Material adapted from Radford Neal s tutorial (http://ftp.cs.utoronto.ca/pub/radford/bayes-tut.pdf), Zoubin Ghahramni (http://hunch.net/~coms-4771/zoubin_ghahramani_bayesian_learning.pdf),
More informationSimultaneous Inversion of Pre-Stack Seismic Data
6 th International Conference & Exposition on Petroleum Geophysics Kolkata 006 Summary Simultaneous Inversion of Pre-Stack Seismic Data Brian H. Russell, Daniel P. Hampson, Brad Bankhead Hampson-Russell
More informationIntroduction to Machine Learning
Introduction to Machine Learning Brown University CSCI 1950-F, Spring 2012 Prof. Erik Sudderth Lecture 25: Markov Chain Monte Carlo (MCMC) Course Review and Advanced Topics Many figures courtesy Kevin
More informationP235 Modelling Anisotropy for Improved Velocities, Synthetics and Well Ties
P235 Modelling Anisotropy for Improved Velocities, Synthetics and Well Ties P.W. Wild* (Ikon Science Ltd), M. Kemper (Ikon Science Ltd), L. Lu (Ikon Science Ltd) & C.D. MacBeth (Heriot Watt University)
More informationA graph contains a set of nodes (vertices) connected by links (edges or arcs)
BOLTZMANN MACHINES Generative Models Graphical Models A graph contains a set of nodes (vertices) connected by links (edges or arcs) In a probabilistic graphical model, each node represents a random variable,
More informationSemi-Markov/Graph Cuts
Semi-Markov/Graph Cuts Alireza Shafaei University of British Columbia August, 2015 1 / 30 A Quick Review For a general chain-structured UGM we have: n n p(x 1, x 2,..., x n ) φ i (x i ) φ i,i 1 (x i, x
More informationBayesian Inference: Principles and Practice 3. Sparse Bayesian Models and the Relevance Vector Machine
Bayesian Inference: Principles and Practice 3. Sparse Bayesian Models and the Relevance Vector Machine Mike Tipping Gaussian prior Marginal prior: single α Independent α Cambridge, UK Lecture 3: Overview
More informationOptimal Experimental Design (Survey)
Optimal Experimental Design (Survey) Sudelfeld, Germany, March 2018 Andrew Curtis University of Edinburgh Guest Professor ETH Zurich Survey and Experimental Design Methods Can not avoid: potentially one
More informationThe Origin of Deep Learning. Lili Mou Jan, 2015
The Origin of Deep Learning Lili Mou Jan, 2015 Acknowledgment Most of the materials come from G. E. Hinton s online course. Outline Introduction Preliminary Boltzmann Machines and RBMs Deep Belief Nets
More informationIJMGE Int. J. Min. & Geo-Eng. Vol.49, No.1, June 2015, pp
IJMGE Int. J. Min. & Geo-Eng. Vol.49, No.1, June 2015, pp.131-142 Joint Bayesian Stochastic Inversion of Well Logs and Seismic Data for Volumetric Uncertainty Analysis Moslem Moradi 1, Omid Asghari 1,
More informationBertrand Six, Olivier Colnard, Jean-Philippe Coulon and Yasmine Aziez CGGVeritas Frédéric Cailly, Total
4-D Seismic Inversion: A Case Study Offshore Congo Bertrand Six, Olivier Colnard, Jean-Philippe Coulon and Yasmine Aziez CGGVeritas Frédéric Cailly, Total Summary The first 4D seismic survey in Congo was
More informationOTC OTC PP. Abstract
OTC OTC-19977-PP Using Modern Geophysical Technology to Explore for Bypassed Opportunities in the Gulf of Mexico R.A. Young/eSeis; W.G. Holt, G. Klefstad/ Fairways Offshore Exploration Copyright 2009,
More informationMachine Learning for Signal Processing Bayes Classification and Regression
Machine Learning for Signal Processing Bayes Classification and Regression Instructor: Bhiksha Raj 11755/18797 1 Recap: KNN A very effective and simple way of performing classification Simple model: For
More informationMachine Learning Summer School
Machine Learning Summer School Lecture 3: Learning parameters and structure Zoubin Ghahramani zoubin@eng.cam.ac.uk http://learning.eng.cam.ac.uk/zoubin/ Department of Engineering University of Cambridge,
More informationReservoir Characterization using AVO and Seismic Inversion Techniques
P-205 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 informationBayesian closed skew Gaussian inversion of seismic AVO data into elastic material properties
Bayesian closed skew Gaussian inversion of seismic AVO data into elastic material properties Omid Karimi 1,2, Henning Omre 2 and Mohsen Mohammadzadeh 1 1 Tarbiat Modares University, Iran 2 Norwegian University
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 informationCondensed Table of Contents for Introduction to Stochastic Search and Optimization: Estimation, Simulation, and Control by J. C.
Condensed Table of Contents for Introduction to Stochastic Search and Optimization: Estimation, Simulation, and Control by J. C. Spall John Wiley and Sons, Inc., 2003 Preface... xiii 1. Stochastic Search
More informationBayesian Machine Learning - Lecture 7
Bayesian Machine Learning - Lecture 7 Guido Sanguinetti Institute for Adaptive and Neural Computation School of Informatics University of Edinburgh gsanguin@inf.ed.ac.uk March 4, 2015 Today s lecture 1
More informationHuman Pose Tracking I: Basics. David Fleet University of Toronto
Human Pose Tracking I: Basics David Fleet University of Toronto CIFAR Summer School, 2009 Looking at People Challenges: Complex pose / motion People have many degrees of freedom, comprising an articulated
More information4D stress sensitivity of dry rock frame moduli: constraints from geomechanical integration
Title 4D stress sensitivity of dry rock frame moduli: constraints from geomechanical integration Authors Bloomer, D., Ikon Science Asia Pacific Reynolds, S., Ikon Science Asia Pacific Pavlova, M., Origin
More informationProbabilistic Graphical Models
Probabilistic Graphical Models Brown University CSCI 2950-P, Spring 2013 Prof. Erik Sudderth Lecture 13: Learning in Gaussian Graphical Models, Non-Gaussian Inference, Monte Carlo Methods Some figures
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 informationDiscrete Inference and Learning Lecture 3
Discrete Inference and Learning Lecture 3 MVA 2017 2018 h
More information