A008 THE PROBABILITY PERTURBATION METHOD AN ALTERNATIVE TO A TRADITIONAL BAYESIAN APPROACH FOR SOLVING INVERSE PROBLEMS
|
|
- Damon Gilbert
- 6 years ago
- Views:
Transcription
1 A008 THE PROAILITY PERTURATION METHOD AN ALTERNATIVE TO A TRADITIONAL AYESIAN APPROAH FOR SOLVING INVERSE PROLEMS Jef AERS Stanford University, Petroleum Engineering, Stanford A USA Abstract ayesian inverse theory provides a framework for solving inverse problems that are non-unique, e.g. history matching. The ayesian approach relies on the fact that the conditional probability of the model parameters given the data (the posterior) is proportional to the likelihood of observing the data and a prior belief expressed as a prior distribution of the model parameters. In case the prior is not Gaussian and the relation between data and parameters (forward model) is strongly non-linear, one has to resort to iterative samplers, often Markov chain Monte arlo methods, for generating samples that fit the data and reflect the prior model statistics. While theoretically sound, such methods can be slow to converge, and are often impractical when the forward model is PU demanding. In this paper we propose a new sampling method that allows to sample from a variety of priors and condition model parameters to a variety of data types. The method does not rely on the traditional ayesian decomposition of posterior into likelihood and prior, instead it uses socalled pre-posterior distributions, i.e. the probability of the model parameters given some subset of the data. The use of pre-posterior allows to decompose the data into so-called, easy data (or linear data) and difficult data (or non-linear data). The method relies on fast non-iterative sequential simulation to generate model realizations, instead of an iterative sampling. The difficult data is matched by a perturbing an initial realization using a perturbation mechanism termed probability perturbation. A simple example is used to illustrate the theoretical development, while a reservoir study is provided in an accompanying paper. Introduction Inverse problems are ubiquitous in the Petroleum Sciences. Sets of measurements d are used to determine the spatial distribution of a physical attribute, described mathematically by a model with a set of parameters m. In most applications, the measurement d is too sparse to determine uniquely the underlying subsurface phenomenon or the model m. For example, a well-test pressure (d) does not uniquely determine subsurface permeability at every location (m) unless the reservoir is purely homogeneous (single m). A set of geophysical measurements such as seismic data does not uniquely determine the subsurface porosity distribution. Most inverse problems are therefore fundamentally underdetermined (ill-posed), with many alternative solutions fitting the same data d. For this reason, a stochastic approach to solving the inverse problem is taken. Instead of determining a single solution, one generates a set of solutions distributed according to a probability distribution function. More specifically, in a ayesian context one aims at sampling such solutions from the posterior density distribution of the model parameters m given the data d, 9th European onference on the Mathematics of Oil Recovery annes, France, 30 August - 2 September 2004
2 2 ( ) ( ) ( ) f d m f f m m d () f ( d) The prior density f(m) describes the dependency between the model parameters and therefore constrains the set of possible inverse solutions. In a spatial context such dependency refers to the spatial structure of m, be it a covariance, oolean or training image-defined dependency. The likelihood density f(d m) models the stochastic relationship between the observed data d and each particular model m retained. This likelihood would account for model and measurement errors. Void of any such errors, the data and model parameters are related through a forward model g, d g( m ) (2) The density f(d) depends on the density f(m) and the forward model g, but its specification can often be avoided in most samplers of f(m d). The aim of ayesian inverse methods is to obtain samples of the posterior density f(m d). Other than the case where g is linear, such sampling will need to be iterative. Popular sampling methods are reection sampling and the Metropolis sampler (Omre and Telmeland, 996; Hegstad and Omre, 996). oth are Markov chain Monte arlo (McM) samplers that avoid the specification of f(d) but are iterative in nature in order to obtain a single sample m (l) of f(m d). Generating multiple (conditioned to d) samples m (l), l,,l in this manner quantifies the uncertainty modeled in f(m d). The problem with such approaches is that they either rely on simple probability models such as multi-gaussian distributions (Tarantola, 989) or require a very large number of forward model evaluations, hence are impractical when g is a flow simulation model. In this paper, we also treat the more general case where two types of dataset d and d 2 are available. In such case the ayesian decomposition of () is written as f( d, d2 m) f( m) f ( m d, d2) f ( d, d2) (3) f( d m) f( d2 m) f( m) ; f ( d, d2) where the latter expression relies on a conditional independence hypothesis between the two data types. We will consider the data d as the easy data, i.e. data than can be conditioned to within a more classical geostatistical context. d could consist of a set of point measurements (termed hard data ) combined with secondary information (termed soft data ) whose relationship to the model m is linear or pseudo-linear in nature. The data d 2 is termed the difficult data, such as production data, exhibiting a complex multi-point and non-linear relationship with the model m. In this paper we present a practical method, termed probability perturbation, that addresses the above mentioned limitations of traditional ayesian approaches. The method relies on: () The use of sequential simulation (Deutsch and Journel, 998) to sample the prior and posterior distributions. Sequential simulation is not iterative, hence more PU efficient than iterative McM samplers. (2) The use of so-called pre-posterior distributions f( m d) and f( m d 2), instead of likelihoods. We will show that this leads to an efficient sampling method, primarily by avoiding the calculation of f ( d, d 2). The purpose of this paper is to present the theory behind the probability perturbation method. For a practical reservoir example, we refer to an accompanying paper by Hoffman and aers (2004)
3 3 Methodology Sampling the prior Practical inverse methods should be able to incorporate a diversity of prior models, not limited to Gaussian ones. To emphasize the non-gaussianity of the proposed method, we consider a binary random function, modeling the absence or presence of a geological event, if the "event" occurs at u I( u ) 0 else where event could represent any spatially distributed, whether continuous or categorical phenomenon. The random function is discretized on a grid composed of a finite set of N grid node locations u i (x i,y i,z i ). The parameters of the inverse problem are given by the unknown indicator variables at each grid node location: m {( I u), I( u2), K, I( u N )} The grid need not be rectangular. The prior distribution is then simply the oint distribution of all indicator random variables at all grid node locations f( m) Prob{ I( u), I( u2), K, I( u N ) } A fast and general method for sampling a prior distribution is sequential simulation. Unlike McM methods, sequential simulation is non-iterative and is completed after a single pass over all grid node locations. Sequential simulation relies on the following decomposition of the prior: f( m) Prob{ I( u) } Prob{ I( u2) i( u)} K (4) Prob{ I( u ) i( u ), K, i( u )} N N This decomposition has led to a series of practical algorithms that can sample a wide variety of prior models (e.g. Deutsch and Journel, 998; Strebelle, 2002) ayesian inversion Priors need to be conditioned to data d and d 2, resulting in a posterior distribution of the model parameters given the data, as formulated in ayes rule () or (3). In determining the posterior, we will follow the same philosophy as in determining the prior, i.e. rely on a sequential decomposition into conditionals. We consider the case where the data d constitutes direct observations (hard data) of the model parameters at a set of n < N spatial locations (hard data), in the binary random function case, d {( i u α), α, K, n} d 2 contains any type of data that has a non-linear relationship with the model parameters d2 g( m) g( I( u), I( u2), K, I( u N )) (5) with g is the forward model. We will assume for now that there are no data and model errors. One is interested in generating samples from a posterior distribution: f( m d, d ) Prob{ I( u ), I( u ), K, I( u ) { i( u ), α, K, n}, d } 2 2 N α 2 To sample such posterior, we do not follow the traditional ayesian route of likelihood and prior. Instead, we rely on the a sequential decomposition similar as in the prior case: f( m d, d2) Prob{ I( u) { i( uα), α, K, n}, d2} Prob{ I( u2) i( u),{ i( uα), α, K, n}, d2} K (6) Prob{ I( u ) i( u ), K, i( u ),{ i( u ), α, K, n}, d } N N α 2 9th European onference on the Mathematics of Oil Recovery annes, France, 30 August - 2 September 2004
4 4 Sampling from the complex (oint) posterior f is equivalent to sequential sampling from a series of univariate conditionals of the type: Prob{ I( u ) i( u), K, i( u ),{ i( uα), α, K, n}, d2} Prob( A, ) with A { I( u ) } { i( u), K, i( u ),{ i( uα), α, K, n}} d2 where we have introduced a simpler notation A (unknown model parameter), data (easy data) and data (difficult data) to make further development clear. Since it is difficult to state the conditionals Prob( A, ) explicitly, we decompose it further into two pre-posterior distributions, Prob( A ) (pre-posterior to knowing ) and Prob( A ) (pre-posterior to knowing ), using Journel s (Journel, 2002) combination method: τ τ2 x b c Prob( A, ) with where: + x a a a (7) -Prob( A ) -Prob( A ) -Prob( A) b, c, a Prob( A ) Prob( A ) Prob( A ) In case, τ τ 2, Journel shows that Eq. (6) is equivalent (up to a simple standardization) to the hypothesis of conditional independence of Eq. (2). Assuming the prior Prob( A ) is known, the problem of stating the conditional Prob( A, ) is now decomposed into a problem of stating the pre-posteriors Prob( A ) and Prob( A ). Working with pre-posteriors will lead to an approach that is different from a classical ayesian inversion which would involve the likelihoods Prob( A) and Prob( A ). This difference may seem subtle but has important practical consequences. Stating pre-posteriors, instead of likelihoods, allows using (noniterative) sequential simulation via Eqns. (6) and (7), instead of (iterative) McM. The τ-values in Eq. (7) allow the user to model explicitly the dependency between the -data and -data. The τ-values can be interpreted as weights given to each data type (see Journel, 2002). Assuming a conditional independence (τ τ 2 ) results in a very particular dependency model that may often not reflect the actual dependency between and data. In the context of sequential simulation, the pre-posterior Prob( A ) is simply the conditional distribution of any previously simulated nodes i( u), K, i( u ) and the hard data i( u α), α, K, n. The remaining pre-posterior Prob( A ) cannot be directly estimated due to the complex nonlinear relationship (forward model) between A and. Instead, in the next section, we propose an iterative sampling method for determining this probability, termed probability perturbation.
5 5 Sampling Prob( A, ) : the probability perturbation method Given a random seed s, and given the pre-posterior Prob( A ) at all grid nodes u, a ( ) ( ) ( ) ( ) realization i s { s ( ), s ( 2), K, s i u i u i ( u N)} can be drawn by ways of sequential simulation. The subscript in i ( s) emphasizes that i ( s) is conditioned to data d (-data) only, not yet to data d 2 (-data), this would require somehow the use of the other pre-posterior, Prob( A ). To generate a realization i (l) matching both d and d 2, each grid node should be sequentially sampled from the oint distribution Prob( A, ). However, the oint distribution is unknown, since Prob( A ) is not known. The probability perturbation method performs a stochastic search for those probabilities Prob( A ) (at all N grid nodes) that achieve a match to the d 2 data. The key idea is to search for all N probabilities Prob( A ), such that, after combining Prob( A ) with Prob( A ) into Prob( A, ) using Eq. (7), a realization i (l) sequentially drawn from Prob( A, ),, K, N, matches the data d 2. Searching for all N probabilities Prob( A ) directly is impossible, particularly since N can be large. Therefore we rely on a perturbation parameterization of all N Prob( A ) using a single parameter as follows: ( ) Prob( A ) Prob( ( ) ) ( ) s I u r i ( u ) + rp( A),, K, N (8) where r is a parameter between [0,], not dependent on u. This independence on u is will be relaxed later. Given Eq. (5), Prob( A ) can be calculated for a given value of r and for a given initial realization ( s i ) ( u ). The importance of Eq. (8) is that it translates the search for N probabilities Prob( A ) into an optimization problem of a single parameter r as follows: hoose a random seed s Generate an initial realization i ( s) using the data d and random seed s Iterate: Until the data d 2 is matched, do the following:. hoose another random seed s 2. Determine a realization i ( s ) that matched better data d 2 as follows: From Prob( A ) in Eq. (8), Prob( A, ) can be determined using Eq. (7). This allows generating a realization ( s ) ( s ) ( s ) ( s ) i { i ( u ), i ( u ), K, i ( u )} drawn by sequentially sampling from r 2 N Prob( A, ),, K, N. This realization is dependent on the value of r, as well as a random seed s. To find an optimal r, hence per Eq. (8) best Prob( A ), the following obective function is formulated ( s ) Or (, s) d2 g( i r ) (9) which measures the distance between the data d 2 and the forward model evaluated on the realization ( s i ). The obective function (9) depends on r r 9th European onference on the Mathematics of Oil Recovery annes, France, 30 August - 2 September 2004
6 6 and a fixed random seed s. A simple one-dimensional optimization can be performed to find the value of r that best matches the data d 2. Once an optimal value r opt is found, generate a realization ( s i ) opt to be used r in Eq. (8) during the next iteration step. The two limit values r 0 and r clarify the choice of the parameterization in Eq. (8) ( ) () In case r 0, then Prob( A ) s i ( u ), hence per Eq (7) ( ) Prob( A, ) i s ( u ). Regardless of the random seeds s and s, any realization ( s ) ( s) r 0 i i. In other words, r 0 entails no perturbation of i ( s). (2) In case r, then Prob( A ) P( A), a simple calculation using Eq. (7) shows that Prob( A, ) PA ( ). Since the seed s is different from the seed s, the ( s ) ( ) realization s ir i is equiprobable with i ( s). In other words, r entails a maximum perturbation within the prior model constraints. A value r in between (0,) will therefore generate a perturbation ( s i ) r between some initial realization i ( s) and another equiprobable realization i ( s ) both conditioned to the data d and each following the prior model statistics. The optimization of r in (9), has to be repeated for multiple random seeds s since a single opt optimization of Or (, s ) with fixed random seed is likely not to reach satisfactory match to the data d 2. A simple illustrative example A simple example is presented to clarify the approach and illustrate various properties of the method in finding inverse solutions. The model consists of a grid with three nodes, u, u 2 and u 3. Each node can be either black, I(u) or white, I(u)0. The model m is therefore simply m {( I u), I( u2), I( u 3)} The spatial dependency of this simple D model is described by a training image shown in Figure. One can extract, by scanning the training image with a 3 template, the prior distribution (f(m)) of the model m, as shown in Figure. To test the probability perturbation method we consider two data: the first data is a point measurement (-data, or easy data ) namely, i(u 2 ), the second one is I( u) + I( u2) + I( u 3) 2 (-data or difficult data ). The problem posed is: What is Pr( I( u) i( u2), I( u) + I( u2) + I( u 3) 2)? Or in simpler notation, what is Pr( A, )? We consider three ways to solve this problem:. alculate the true posterior Pr( A, ) directly by simple elimination of those prior models that do not match the and data. 2. alculate Pr( A, ) by first calculating Pr( A ) and Pr( A ) based on the prior, then using Journel s Eq. (7) under conditional independence, namely τ, τ alculate Pr( A, ) by running a Monte arlo experiment using the probability perturbation method (PPM). On each realization drawn from the prior model, the PPM using the algorithm described above is applied. From the set of matched realizations, Pr( A, ) is calculated by checking how many times i(u ).
7 7 The solutions are as follows:. Directly: Pr( A, ) 3 2. Using Journel s Equation Pr( A) a, Pr( A) b, Pr( A ) c Applying Eq. (7) with τ τ 2: x Pr( A, ) Using Monte arlo on the PPM method, Pr( A, ) 0.39 The following observations can be made: omparing the results of [.] and [2.], it is clear that the assumption of conditional independence is not valid for this case. A simple calculation shows that using τ 2.59 and τ 2 in Eq. (7) would provide an exact approximation of the true posterior. This indicates that the single hard conditioning data should receive substantially more weight than assumed under conditional independence. The relative differences, between P(A), P(A ), P(A ) and P(A,) indicate that P(A ) carries more relevance in determining Prob(A,) then Prob(A ). In other words, the single hard () data carries more relevant information about the unknown A than the sum-equals-two () information. The difference of results between methods [2.] and [3.] can be explained as follows: in the PPM method, perturbations are created by means of sequential simulation using Eqns. (6), (7) and (8). In sequential simulation, the -data (hard data at u 2 and any previously simulated values) changes at every node visited and depend on the order of visited nodes (u first, then u 3 or vice versa). In method [2.] the -data is fixed and consists of the single conditioning value at u 2. While the PPM relies on the same assumption of conditional independence, the result is much closer to the true posterior probability. This result was confirmed by applying a variety of other training images, i.e. other prior distributions. Method [2.] produces consistently a considerable overestimation of Prob(A,), the PPM is off by a few percentages only. The reason for the latter observation can be explained by means of Eq (8). In this equation, the pre-posterior Prob( I( u ) ) is a function of the data through the parameter r, and, a function of an initial realization {i(u ),i(u 2 ),i(u 3 )}. This initial realization depends itself on the pre-posterior Prob( I( u ) ) with depending on the random path. Hence, Eq. (8) forces an explicit dependency between the Prob( I( u ) ) and Prob( I( u ) ) prior to combining both into Prob( I( u ), ) using a conditional independence hypothesis (Eq. (7) with τ, τ 2 ). At least from this simple example, one can conclude that the sequential decomposition of the posterior into pre-posteriors has robustified the estimate of the true posterior under the conditional independence hypothesis. onclusion 9th European onference on the Mathematics of Oil Recovery annes, France, 30 August - 2 September 2004
8 8 The probability perturbation method shares two important properties of other inverse algorithms: () Inverse solutions generated using the probability perturbation method honor the prior statistics. Whatever value of r or random seed s, all model realizations are drawn using a sequential sampling method that depends on Prob( A ) (depending on prior statistics). The probability perturbation is essentially a search method for finding those prior model realizations that honor the data d 2. (2) The probability perturbation method can cover the entire space of the prior model realizations. At each step of the outer loop, the random seed is changed, hence, if infinite computing time were available, the method would cover the entire space of possible realizations, ust like a reection sampler would. References aers, J., (2003). History matching under a training image-based geological model constraint. SPE Journal, SPE # 7476, p Deutsch,.V. and Journel, A.G. (998). GSLI: the geostatistical software library. Oxford University Press. Hegstad,.K. and Omre, H. (996). Uncertainty assessment in history matching and forecasting. In Proceeding of the Fifth International Geostatistics ongress, ed. E.Y aafi and N.A. Schofield, Wollongong Australia, v., p Hoffman,.T. and aers, J., History matching with the regional probability perturbation method applications to a North Sea reservoir In: Proceedings to the EMOR IX, annes, Aug 29 - Sept 2, Journel, A.G. (2002). ombining knowledge from diverse data sources: an alternative to traditional data independence hypothesis. Math. Geol., v Omre, H. and Telmeland, H. Petroleum Geostatistics. In Proceeding of the Fifth International Geostatistics ongress, ed. E.Y aafi and N.A. Schofield, Wollongong Australia, v., p Strebelle, S., (2002). onditional simulation of complex geological structures using multiple-point geostatistics. Math. Geol., v34, p-22. Tarantola, A. (989). Inverse problem theory. Kluwer Academic Publisher, Dordrecht Training image Prior probabilities of the following realizations: Figure : prior probabilities extracted from the training image
Comparing the gradual deformation with the probability perturbation method
Comparing the gradual deformation with the probability perturbation method Jef Caers Petroleum Engineering Department Stanford University, California USA 1. Introduction This paper is intended as a comparison
More informationMultiple-Point Geostatistics: from Theory to Practice Sebastien Strebelle 1
Multiple-Point Geostatistics: from Theory to Practice Sebastien Strebelle 1 Abstract The limitations of variogram-based simulation programs to model complex, yet fairly common, geological elements, e.g.
More informationLinear inverse Gaussian theory and geostatistics a tomography example København Ø,
Linear inverse Gaussian theory and geostatistics a tomography example Thomas Mejer Hansen 1, ndre Journel 2, lbert Tarantola 3 and Klaus Mosegaard 1 1 Niels Bohr Institute, University of Copenhagen, Juliane
More informationA Short Note on the Proportional Effect and Direct Sequential Simulation
A Short Note on the Proportional Effect and Direct Sequential Simulation Abstract B. Oz (boz@ualberta.ca) and C. V. Deutsch (cdeutsch@ualberta.ca) University of Alberta, Edmonton, Alberta, CANADA Direct
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 informationEntropy of Gaussian Random Functions and Consequences in Geostatistics
Entropy of Gaussian Random Functions and Consequences in Geostatistics Paula Larrondo (larrondo@ualberta.ca) Department of Civil & Environmental Engineering University of Alberta Abstract Sequential Gaussian
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 informationGeostatistical History Matching coupled with Adaptive Stochastic Sampling: A zonation-based approach using Direct Sequential Simulation
Geostatistical History Matching coupled with Adaptive Stochastic Sampling: A zonation-based approach using Direct Sequential Simulation Eduardo Barrela* Instituto Superior Técnico, Av. Rovisco Pais 1,
More informationLarge Scale Modeling by Bayesian Updating Techniques
Large Scale Modeling by Bayesian Updating Techniques Weishan Ren Centre for Computational Geostatistics Department of Civil and Environmental Engineering University of Alberta Large scale models are useful
More informationUSING GEOSTATISTICS TO DESCRIBE COMPLEX A PRIORI INFORMATION FOR INVERSE PROBLEMS THOMAS M. HANSEN 1,2, KLAUS MOSEGAARD 2 and KNUD S.
USING GEOSTATISTICS TO DESCRIBE COMPLEX A PRIORI INFORMATION FOR INVERSE PROBLEMS THOMAS M. HANSEN 1,2, KLAUS MOSEGAARD 2 and KNUD S. CORDUA 1 1 Institute of Geography & Geology, University of Copenhagen,
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 informationB008 COMPARISON OF METHODS FOR DOWNSCALING OF COARSE SCALE PERMEABILITY ESTIMATES
1 B8 COMPARISON OF METHODS FOR DOWNSCALING OF COARSE SCALE PERMEABILITY ESTIMATES Alv-Arne Grimstad 1 and Trond Mannseth 1,2 1 RF-Rogaland Research 2 Now with CIPR - Centre for Integrated Petroleum Research,
More information4th HR-HU and 15th HU geomathematical congress Geomathematics as Geoscience Reliability enhancement of groundwater estimations
Reliability enhancement of groundwater estimations Zoltán Zsolt Fehér 1,2, János Rakonczai 1, 1 Institute of Geoscience, University of Szeged, H-6722 Szeged, Hungary, 2 e-mail: zzfeher@geo.u-szeged.hu
More informationOptimizing Thresholds in Truncated Pluri-Gaussian Simulation
Optimizing Thresholds in Truncated Pluri-Gaussian Simulation Samaneh Sadeghi and Jeff B. Boisvert Truncated pluri-gaussian simulation (TPGS) is an extension of truncated Gaussian simulation. This method
More informationA013 HISTORY MATCHING WITH RESPECT TO RESERVOIR STRUCTURE
A3 HISTORY MATCHING WITH RESPECT TO RESERVOIR STRUCTURE SIGURD IVAR AANONSEN ; ODDVAR LIA ; AND OLE JAKOB ARNTZEN Centre for Integrated Research, University of Bergen, Allégt. 4, N-7 Bergen, Norway Statoil
More informationInverting hydraulic heads in an alluvial aquifer constrained with ERT data through MPS and PPM: a case study
Inverting hydraulic heads in an alluvial aquifer constrained with ERT data through MPS and PPM: a case study Hermans T. 1, Scheidt C. 2, Caers J. 2, Nguyen F. 1 1 University of Liege, Applied Geophysics
More informationA033 PRACTICAL METHODS FOR UNCERTAINTY ASSESSMENT
A33 PRACTICAL METHODS FOR UNCERTAINTY ASSESSMENT OF FLOW PREDICTIONS FOR RESERVOIRS WITH SIGNIFICANT HISTORY AFIELD CASE STUDY ALEXANDRE CASTELLINl, JORGE L. LANDA, JITENDRA KIKANI 2 () ChevronTexaco,
More informationReservoir Uncertainty Calculation by Large Scale Modeling
Reservoir Uncertainty Calculation by Large Scale Modeling Naeem Alshehri and Clayton V. Deutsch It is important to have a good estimate of the amount of oil or gas in a reservoir. The uncertainty in reserve
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 informationSequential Simulations of Mixed Discrete-Continuous Properties: Sequential Gaussian Mixture Simulation
Sequential Simulations of Mixed Discrete-Continuous Properties: Sequential Gaussian Mixture Simulation Dario Grana, Tapan Mukerji, Laura Dovera, and Ernesto Della Rossa Abstract We present here a method
More informationAdvances in Locally Varying Anisotropy With MDS
Paper 102, CCG Annual Report 11, 2009 ( 2009) Advances in Locally Varying Anisotropy With MDS J.B. Boisvert and C. V. Deutsch Often, geology displays non-linear features such as veins, channels or folds/faults
More informationCombining geological surface data and geostatistical model for Enhanced Subsurface geological model
Combining geological surface data and geostatistical model for Enhanced Subsurface geological model M. Kurniawan Alfadli, Nanda Natasia, Iyan Haryanto Faculty of Geological Engineering Jalan Raya Bandung
More informationBuilding Blocks for Direct Sequential Simulation on Unstructured Grids
Building Blocks for Direct Sequential Simulation on Unstructured Grids Abstract M. J. Pyrcz (mpyrcz@ualberta.ca) and C. V. Deutsch (cdeutsch@ualberta.ca) University of Alberta, Edmonton, Alberta, CANADA
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 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 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 informationA Statistical Input Pruning Method for Artificial Neural Networks Used in Environmental Modelling
A Statistical Input Pruning Method for Artificial Neural Networks Used in Environmental Modelling G. B. Kingston, H. R. Maier and M. F. Lambert Centre for Applied Modelling in Water Engineering, School
More informationMultiple realizations using standard inversion techniques a
Multiple realizations using standard inversion techniques a a Published in SEP report, 105, 67-78, (2000) Robert G Clapp 1 INTRODUCTION When solving a missing data problem, geophysicists and geostatisticians
More informationCharacterization of Geoobjects Continuity using Moments of Inertia
Characterization of Geoobjects Continuity using Moments of Inertia Saina Lajevardi, Olena Babak, and Clayton V. Deutsch Well-placement is one of the main challenges in reservoir engineering. The connectivity
More informationRC 3.3. Summary. Introduction
Application of multipoint geostatistics to honor multiple attribute constraints applied to a deepwater outcrop analog, Tanqua Karoo Basin, South Africa Daniel Tetzlaff, Roy Davies, David McCormick, Claude
More informationA MultiGaussian Approach to Assess Block Grade Uncertainty
A MultiGaussian Approach to Assess Block Grade Uncertainty Julián M. Ortiz 1, Oy Leuangthong 2, and Clayton V. Deutsch 2 1 Department of Mining Engineering, University of Chile 2 Department of Civil &
More informationConditional Simulation of Random Fields by Successive Residuals 1
Mathematical Geology, Vol 34, No 5, July 2002 ( C 2002) Conditional Simulation of Random Fields by Successive Residuals 1 J A Vargas-Guzmán 2,3 and R Dimitrakopoulos 2 This paper presents a new approach
More informationPETROLEUM GEOSTATISTICS HENNING OMRE AND HAKON TJELMELAND. Department of Mathematical Sciences. Norwegian University of Science and Technology
PETROLEUM GEOSTATISTICS HENNING OMRE AND HAKON TJELMELAND Department of Mathematical Sciences Norwegian University of Science and Technology Trondheim, Norway Abstract. For reservoir evaluation, general
More informationAdvanced analysis and modelling tools for spatial environmental data. Case study: indoor radon data in Switzerland
EnviroInfo 2004 (Geneva) Sh@ring EnviroInfo 2004 Advanced analysis and modelling tools for spatial environmental data. Case study: indoor radon data in Switzerland Mikhail Kanevski 1, Michel Maignan 1
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 informationBootstrap confidence intervals for reservoir model selection techniques
Bootstrap confidence intervals for reservoir model selection techniques Céline Scheidt and Jef Caers Department of Energy Resources Engineering Stanford University Abstract Stochastic spatial simulation
More informationDirect forecasting without full model inversion Jef Caers
Direct forecasting without full model inversion Jef Caers Professor of Geological Sciences Stanford University, USA Why moving to Geological Sciences? To attract a greater diversity of students To be able
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 informationMonte Carlo Simulation. CWR 6536 Stochastic Subsurface Hydrology
Monte Carlo Simulation CWR 6536 Stochastic Subsurface Hydrology Steps in Monte Carlo Simulation Create input sample space with known distribution, e.g. ensemble of all possible combinations of v, D, q,
More informationStanford Exploration Project, Report 105, September 5, 2000, pages 41 53
Stanford Exploration Project, Report 105, September 5, 2000, pages 41 53 40 Stanford Exploration Project, Report 105, September 5, 2000, pages 41 53 Short Note Multiple realizations using standard inversion
More informationSTA 4273H: Statistical Machine Learning
STA 4273H: Statistical Machine Learning Russ Salakhutdinov Department of Computer Science! Department of Statistical Sciences! rsalakhu@cs.toronto.edu! h0p://www.cs.utoronto.ca/~rsalakhu/ Lecture 7 Approximate
More informationRecent developments in object modelling opens new era for characterization of fluvial reservoirs
Recent developments in object modelling opens new era for characterization of fluvial reservoirs Markus L. Vevle 1*, Arne Skorstad 1 and Julie Vonnet 1 present and discuss different techniques applied
More informationModeling of Atmospheric Effects on InSAR Measurements With the Method of Stochastic Simulation
Modeling of Atmospheric Effects on InSAR Measurements With the Method of Stochastic Simulation Z. W. LI, X. L. DING Department of Land Surveying and Geo-Informatics, Hong Kong Polytechnic University, Hung
More informationLatin Hypercube Sampling with Multidimensional Uniformity
Latin Hypercube Sampling with Multidimensional Uniformity Jared L. Deutsch and Clayton V. Deutsch Complex geostatistical models can only be realized a limited number of times due to large computational
More informationComputer Vision Group Prof. Daniel Cremers. 14. Sampling Methods
Prof. Daniel Cremers 14. Sampling Methods Sampling Methods Sampling Methods are widely used in Computer Science as an approximation of a deterministic algorithm to represent uncertainty without a parametric
More informationCOLLOCATED CO-SIMULATION USING PROBABILITY AGGREGATION
COLLOCATED CO-SIMULATION USING PROBABILITY AGGREGATION G. MARIETHOZ, PH. RENARD, R. FROIDEVAUX 2. CHYN, University of Neuchâtel, rue Emile Argand, CH - 2009 Neuchâtel, Switzerland 2 FSS Consultants, 9,
More informationBest Practice Reservoir Characterization for the Alberta Oil Sands
Best Practice Reservoir Characterization for the Alberta Oil Sands Jason A. McLennan and Clayton V. Deutsch Centre for Computational Geostatistics (CCG) Department of Civil and Environmental Engineering
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 informationPorosity prediction using cokriging with multiple secondary datasets
Cokriging with Multiple Attributes Porosity prediction using cokriging with multiple secondary datasets Hong Xu, Jian Sun, Brian Russell, Kris Innanen ABSTRACT The prediction of porosity is essential for
More informationComputational Challenges in Reservoir Modeling. Sanjay Srinivasan The Pennsylvania State University
Computational Challenges in Reservoir Modeling Sanjay Srinivasan The Pennsylvania State University Well Data 3D view of well paths Inspired by an offshore development 4 platforms 2 vertical wells 2 deviated
More informationEarth models for early exploration stages
ANNUAL MEETING MASTER OF PETROLEUM ENGINEERING Earth models for early exploration stages Ângela Pereira PhD student angela.pereira@tecnico.ulisboa.pt 3/May/2016 Instituto Superior Técnico 1 Outline Motivation
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 informationApplications of Randomized Methods for Decomposing and Simulating from Large Covariance Matrices
Applications of Randomized Methods for Decomposing and Simulating from Large Covariance Matrices Vahid Dehdari and Clayton V. Deutsch Geostatistical modeling involves many variables and many locations.
More informationEffect of velocity uncertainty on amplitude information
Stanford Exploration Project, Report 111, June 9, 2002, pages 253 267 Short Note Effect of velocity uncertainty on amplitude information Robert G. Clapp 1 INTRODUCTION Risk assessment is a key component
More informationGeostatistical History Matching coupled with Adaptive Stochastic Sampling
Geostatistical History Matching coupled with Adaptive Stochastic Sampling A Geologically consistent approach using Stochastic Sequential Simulation Eduardo Barrela Nº 79909 Project Thesis Presentation
More informationThe Tau model for Data Redundancy: Part 2
The Tau model for Data Redundancy: Part 2 Sunderrajan Krishnan April 24, 25 Abstract Knowledge about earth properties arrive from diverse sources of information. Connectivity of permeability is one such
More informationLatent Variable Models
Latent Variable Models Stefano Ermon, Aditya Grover Stanford University Lecture 5 Stefano Ermon, Aditya Grover (AI Lab) Deep Generative Models Lecture 5 1 / 31 Recap of last lecture 1 Autoregressive models:
More informationReservoir Modeling with GSLIB. Overview
Reservoir Modeling with GSLIB Overview Objectives of the Course What is Geostatistics? Why Geostatistics / 3-D Modeling? Uncertainty Quantification and Decision Making Heterogeneous Reservoir Modeling
More informationDynamic Data Driven Simulations in Stochastic Environments
Computing 77, 321 333 (26) Digital Object Identifier (DOI) 1.17/s67-6-165-3 Dynamic Data Driven Simulations in Stochastic Environments C. Douglas, Lexington, Y. Efendiev, R. Ewing, College Station, V.
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 informationComputer Vision Group Prof. Daniel Cremers. 11. Sampling Methods
Prof. Daniel Cremers 11. Sampling Methods Sampling Methods Sampling Methods are widely used in Computer Science as an approximation of a deterministic algorithm to represent uncertainty without a parametric
More informationStepwise Conditional Transformation for Simulation of Multiple Variables 1
Mathematical Geology, Vol. 35, No. 2, February 2003 ( C 2003) Stepwise Conditional Transformation for Simulation of Multiple Variables 1 Oy Leuangthong 2 and Clayton V. Deutsch 2 Most geostatistical studies
More informationB005 A NEW FAST FOURIER TRANSFORM ALGORITHM FOR FLUID FLOW SIMULATION
1 B5 A NEW FAST FOURIER TRANSFORM ALGORITHM FOR FLUID FLOW SIMULATION LUDOVIC RICARD, MICAËLE LE RAVALEC-DUPIN, BENOÎT NOETINGER AND YVES GUÉGUEN Institut Français du Pétrole, 1& 4 avenue Bois Préau, 92852
More informationGeostatistical Determination of Production Uncertainty: Application to Firebag Project
Geostatistical Determination of Production Uncertainty: Application to Firebag Project Abstract C. V. Deutsch, University of Alberta (cdeutsch@civil.ualberta.ca) E. Dembicki and K.C. Yeung, Suncor Energy
More informationBackground. Developing a FracMan DFN Model. Fractures, FracMan and Fragmentation Applications of DFN Models to Block & Panel Caving
Background Golder Associates are one of the pioneering groups in the use of the Discrete Fracture Network (DFN) approach. DFN models seek to describe the heterogeneous nature of fractured rock masses by
More informationCOLLOCATED CO-SIMULATION USING PROBABILITY AGGREGATION
COLLOCATED CO-SIMULATION USING PROBABILITY AGGREGATION G. MARIETHOZ, PH. RENARD, R. FROIDEVAUX 2. CHYN, University of Neuchâtel, rue Emile Argand, CH - 2009 Neuchâtel, Switzerland 2 FSS Consultants, 9,
More informationCorrecting Variogram Reproduction of P-Field Simulation
Correcting Variogram Reproduction of P-Field Simulation Julián M. Ortiz (jmo1@ualberta.ca) Department of Civil & Environmental Engineering University of Alberta Abstract Probability field simulation is
More informationFacies Modeling in Presence of High Resolution Surface-based Reservoir Models
Facies Modeling in Presence of High Resolution Surface-based Reservoir Models Kevin Zhang Centre for Computational Geostatistics Department of Civil and Environmental Engineering University of Alberta
More informationEnsemble Kalman filter for automatic history matching of geologic facies
Journal of Petroleum Science and Engineering 47 (2005) 147 161 www.elsevier.com/locate/petrol Ensemble Kalman filter for automatic history matching of geologic facies Ning LiuT, Dean S. Oliver School of
More informationA Program for Data Transformations and Kernel Density Estimation
A Program for Data Transformations and Kernel Density Estimation John G. Manchuk and Clayton V. Deutsch Modeling applications in geostatistics often involve multiple variables that are not multivariate
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 information3D stochastic inversion of borehole and surface gravity data using Geostatistics
3D stochastic inversion of borehole and surface gravity data using Geostatistics P. Shamsipour ( 1 ), M. Chouteau ( 1 ), D. Marcotte( 1 ), P. Keating( 2 ), ( 1 ) Departement des Genies Civil, Geologique
More informationDownloaded 10/25/16 to Redistribution subject to SEG license or copyright; see Terms of Use at
Facies modeling in unconventional reservoirs using seismic derived facies probabilities Reinaldo J. Michelena*, Omar G. Angola, and Kevin S. Godbey, ireservoir.com, Inc. Summary We present in this paper
More informationSequential Monte Carlo Algorithms
ayesian Phylogenetic Inference using Sequential Monte arlo lgorithms lexandre ouchard-ôté *, Sriram Sankararaman *, and Michael I. Jordan *, * omputer Science ivision, University of alifornia erkeley epartment
More informationBasic Sampling Methods
Basic Sampling Methods Sargur Srihari srihari@cedar.buffalo.edu 1 1. Motivation Topics Intractability in ML How sampling can help 2. Ancestral Sampling Using BNs 3. Transforming a Uniform Distribution
More informationA Reservoir Sampling Algorithm with Adaptive Estimation of Conditional Expectation
A Reservoir Sampling Algorithm with Adaptive Estimation of Conditional Expectation Vu Malbasa and Slobodan Vucetic Abstract Resource-constrained data mining introduces many constraints when learning from
More informationUniversity of Alberta
University of Alberta Library Release Form Name of Author: Linan Zhang Title of Thesis: Production Data Integration in Geostatistical Reservoir Modeling Degree: Master of Science Year this Degree Granted:
More informationReliability of Seismic Data for Hydrocarbon Reservoir Characterization
Reliability of Seismic Data for Hydrocarbon Reservoir Characterization Geetartha Dutta (gdutta@stanford.edu) December 10, 2015 Abstract Seismic data helps in better characterization of hydrocarbon reservoirs.
More informationPETROLEUM RESERVOIR CHARACTERIZATION
PETROLEUM RESERVOIR CHARACTERIZATION Eliandro Admir da Rocha Marques ABSTRACT The reservoir simulation in the oil industry has become the standard for the solution of engineering problems of reservoirs.
More informationA010 MULTISCALE RESERVOIR CHARACTERIZATION USING
1 A010 MULTISCALE RESERVOIR CHARACTERIZATION USING RODUCTION AND TIME LASE SEISMIC DATA Mokhles MEZGHANI, Alexandre FORNEL, Valérie LANGLAIS, Nathalie LUCET IF, 1 & 4 av de Bois réau, 92852 RUEIL-MALMAISON
More informationEstimating the hydrocarbon volume from elastic and resistivity data: A concept
INTERPRETER S CORNER Coordinated by Rebecca B. Latimer Estimating the hydrocarbon volume from elastic and resistivity data: A concept CARMEN T. GOMEZ, JACK DVORKIN, and GARY MAVKO, Stanford University,
More information(5) Multi-parameter models - Gibbs sampling. ST440/540: Applied Bayesian Analysis
Summarizing a posterior Given the data and prior the posterior is determined Summarizing the posterior gives parameter estimates, intervals, and hypothesis tests Most of these computations are integrals
More informationMCMC Sampling for Bayesian Inference using L1-type Priors
MÜNSTER MCMC Sampling for Bayesian Inference using L1-type Priors (what I do whenever the ill-posedness of EEG/MEG is just not frustrating enough!) AG Imaging Seminar Felix Lucka 26.06.2012 , MÜNSTER Sampling
More informationMonte Carlo full-waveform inversion of crosshole GPR data using multiple-point geostatistical a priori information
Downloaded from orbit.dtu.dk on: Dec, 0 Monte Carlo full-waveform inversion of crosshole GPR data using multiple-point geostatistical a priori information Cordua, Knud Skou; Hansen, Thomas Mejer; Mosegaard,
More informationEnKF Applications to Thermal Petroleum Reservoir Characterization
EnKF Applications to Thermal Petroleum Reservoir Characterization Yevgeniy Zagayevskiy and Clayton V. Deutsch This paper presents a summary of the application of the ensemble Kalman filter (EnKF) to petroleum
More informationDownscaling Seismic Data to the Meter Scale: Sampling and Marginalization. Subhash Kalla LSU Christopher D. White LSU James S.
Downscaling Seismic Data to the Meter Scale: Sampling and Marginalization Subhash Kalla LSU Christopher D. White LSU James S. Gunning CSIRO Contents Context of this research Background Data integration
More information27 : Distributed Monte Carlo Markov Chain. 1 Recap of MCMC and Naive Parallel Gibbs Sampling
10-708: Probabilistic Graphical Models 10-708, Spring 2014 27 : Distributed Monte Carlo Markov Chain Lecturer: Eric P. Xing Scribes: Pengtao Xie, Khoa Luu In this scribe, we are going to review the Parallel
More informationGEOINFORMATICS Vol. II - Stochastic Modelling of Spatio-Temporal Phenomena in Earth Sciences - Soares, A.
STOCHASTIC MODELLING OF SPATIOTEMPORAL PHENOMENA IN EARTH SCIENCES Soares, A. CMRP Instituto Superior Técnico, University of Lisbon. Portugal Keywords Spacetime models, geostatistics, stochastic simulation
More informationAssessment of Hydraulic Conductivity Upscaling Techniques and. Associated Uncertainty
CMWRXVI Assessment of Hydraulic Conductivity Upscaling Techniques and Associated Uncertainty FARAG BOTROS,, 4, AHMED HASSAN 3, 4, AND GREG POHLL Division of Hydrologic Sciences, University of Nevada, Reno
More informationMULTISCALE FINITE ELEMENT METHODS FOR STOCHASTIC POROUS MEDIA FLOW EQUATIONS AND APPLICATION TO UNCERTAINTY QUANTIFICATION
MULTISCALE FINITE ELEMENT METHODS FOR STOCHASTIC POROUS MEDIA FLOW EQUATIONS AND APPLICATION TO UNCERTAINTY QUANTIFICATION P. DOSTERT, Y. EFENDIEV, AND T.Y. HOU Abstract. In this paper, we study multiscale
More informationGSLIB Geostatistical Software Library and User's Guide
GSLIB Geostatistical Software Library and User's Guide Second Edition Clayton V. Deutsch Department of Petroleum Engineering Stanford University Andre G. Journel Department of Geological and Environmental
More informationINTEGRATION OF BOREHOLE INFORMATION AND RESISTIVITY DATA FOR AQUIFER VULNERABILITY
INTEGRATION OF BOREHOLE INFORMATION AND RESISTIVITY DATA FOR AQUIFER VULNERABILITY ANDERS V. CHRISTIANSEN 1, ESBEN AUKEN 1 AND KURT SØRENSEN 1 1 HydroGeophysics Group, Aarhus University, Finlandsgade 8,
More informationMultiple realizations: Model variance and data uncertainty
Stanford Exploration Project, Report 108, April 29, 2001, pages 1?? Multiple realizations: Model variance and data uncertainty Robert G. Clapp 1 ABSTRACT Geophysicists typically produce a single model,
More informationConditional Distribution Fitting of High Dimensional Stationary Data
Conditional Distribution Fitting of High Dimensional Stationary Data Miguel Cuba and Oy Leuangthong The second order stationary assumption implies the spatial variability defined by the variogram is constant
More informationSPATIAL-TEMPORAL TECHNIQUES FOR PREDICTION AND COMPRESSION OF SOIL FERTILITY DATA
SPATIAL-TEMPORAL TECHNIQUES FOR PREDICTION AND COMPRESSION OF SOIL FERTILITY DATA D. Pokrajac Center for Information Science and Technology Temple University Philadelphia, Pennsylvania A. Lazarevic Computer
More informationEstimation of Cole-Cole parameters from time-domain electromagnetic data Laurens Beran and Douglas Oldenburg, University of British Columbia.
Estimation of Cole-Cole parameters from time-domain electromagnetic data Laurens Beran and Douglas Oldenburg, University of British Columbia. SUMMARY We present algorithms for the inversion of time-domain
More informationApplied Geostatisitcs Analysis for Reservoir Characterization Based on the SGeMS (Stanford
Open Journal of Yangtze Gas and Oil, 2017, 2, 45-66 http://www.scirp.org/journal/ojogas ISSN Online: 2473-1900 ISSN Print: 2473-1889 Applied Geostatisitcs Analysis for Reservoir Characterization Based
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 informationChecking up on the neighbors: Quantifying uncertainty in relative event location
Checking up on the neighbors: Quantifying uncertainty in relative event location The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation
More informationThe Particle Filter. PD Dr. Rudolph Triebel Computer Vision Group. Machine Learning for Computer Vision
The Particle Filter Non-parametric implementation of Bayes filter Represents the belief (posterior) random state samples. by a set of This representation is approximate. Can represent distributions that
More information