Show simple item record

dc.contributor.advisorTempone, Raul
dc.contributor.authorBen Issaid, Chaouki
dc.date.accessioned2015-05-13T14:09:23Z
dc.date.available2015-05-13T14:09:23Z
dc.date.issued2015-05-12
dc.identifier.doi10.25781/KAUST-U2H82
dc.identifier.urihttp://hdl.handle.net/10754/552705
dc.description.abstractExperimental design can be vital when experiments are resource-exhaustive and time-consuming. In this work, we carry out experimental design in the Bayesian framework. To measure the amount of information that can be extracted from the data in an experiment, we use the expected information gain as the utility function, which specifically is the expected logarithmic ratio between the posterior and prior distributions. Optimizing this utility function enables us to design experiments that yield the most informative data about the model parameters. One of the major difficulties in evaluating the expected information gain is that it naturally involves nested integration over a possibly high dimensional domain. We use the Multilevel Monte Carlo (MLMC) method to accelerate the computation of the nested high dimensional integral. The advantages are twofold. First, MLMC can significantly reduce the cost of the nested integral for a given tolerance, by using an optimal sample distribution among different sample averages of the inner integrals. Second, the MLMC method imposes fewer assumptions, such as the asymptotic concentration of posterior measures, required for instance by the Laplace approximation (LA). We test the MLMC method using two numerical examples. The first example is the design of sensor deployment for a Darcy flow problem governed by a one-dimensional Poisson equation. We place the sensors in the locations where the pressure is measured, and we model the conductivity field as a piecewise constant random vector with two parameters. The second one is chemical Enhanced Oil Recovery (EOR) core flooding experiment assuming homogeneous permeability. We measure the cumulative oil recovery, from a horizontal core flooded by water, surfactant and polymer, for different injection rates. The model parameters consist of the endpoint relative permeabilities, the residual saturations and the relative permeability exponents for the three phases: water, oil and microemulsions. We also compare the performance of the MLMC to the LA and the direct Double Loop Monte Carlo (DLMC). In fact, we show that, in the case of the aforementioned examples, MLMC combined with LA turns to be the best method in terms of computational cost.
dc.language.isoen
dc.subjectbayesian experimental design
dc.subjectMultilevel Monte Carlo
dc.subjectexpected information gain
dc.subjectlaplace approximation
dc.subjectsensor deployment
dc.subjectcore flooding
dc.titleBayesian Optimal Experimental Design Using Multilevel Monte Carlo
dc.typeThesis
dc.contributor.departmentComputer, Electrical and Mathematical Sciences and Engineering (CEMSE) Division
thesis.degree.grantorKing Abdullah University of Science and Technology
dc.contributor.committeememberBisetti, Fabrizio
dc.contributor.committeememberSun, Ying
dc.contributor.committeememberLong, Quan
thesis.degree.disciplineApplied Mathematics and Computational Science
thesis.degree.nameMaster of Science
refterms.dateFOA2016-05-12T00:00:00Z


Files in this item

Thumbnail
Name:
Thesis.pdf
Size:
1.182Mb
Format:
PDF

This item appears in the following Collection(s)

Show simple item record