Optimisation under Uncertainty with Stochastic PDEs for the History Matching Problem in Reservoir Engineering

Similar documents
Parametric Problems, Stochastics, and Identification

Numerical Approximation of Stochastic Elliptic Partial Differential Equations

Hierarchical Parallel Solution of Stochastic Systems

Quantifying Uncertainty: Modern Computational Representation of Probability and Applications

Partial Differential Equations with Stochastic Coefficients

Stochastic Collocation Methods for Polynomial Chaos: Analysis and Applications

At A Glance. UQ16 Mobile App.

At-A-Glance SIAM Events Mobile App

Sampling and low-rank tensor approximation of the response surface

Monte Carlo Methods for Uncertainty Quantification

A Spectral Approach to Linear Bayesian Updating

Introduction to Computational Stochastic Differential Equations

NON-LINEAR APPROXIMATION OF BAYESIAN UPDATE

Organization. I MCMC discussion. I project talks. I Lecture.

MULTISCALE FINITE ELEMENT METHODS FOR STOCHASTIC POROUS MEDIA FLOW EQUATIONS AND APPLICATION TO UNCERTAINTY QUANTIFICATION

Stochastic Finite Elements: Computational Approaches to SPDEs

Gaussian Process Approximations of Stochastic Differential Equations

On a Data Assimilation Method coupling Kalman Filtering, MCRE Concept and PGD Model Reduction for Real-Time Updating of Structural Mechanics Model

Implementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs

Sampling and Low-Rank Tensor Approximations

PRECONDITIONING MARKOV CHAIN MONTE CARLO SIMULATIONS USING COARSE-SCALE MODELS

Stochastic Spectral Approaches to Bayesian Inference

Galerkin Methods for Linear and Nonlinear Elliptic Stochastic Partial Differential Equations

DATA ASSIMILATION FOR FLOOD FORECASTING

Numerical Analysis of Elliptic PDEs with Random Coefficients

NONLOCALITY AND STOCHASTICITY TWO EMERGENT DIRECTIONS FOR APPLIED MATHEMATICS. Max Gunzburger

Solving the Stochastic Steady-State Diffusion Problem Using Multigrid

Utilizing Adjoint-Based Techniques to Improve the Accuracy and Reliability in Uncertainty Quantification

Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen

Overview. Bayesian assimilation of experimental data into simulation (for Goland wing flutter) Why not uncertainty quantification?

INTERPOLATION AND UPDATE IN DYNAMIC DATA-DRIVEN APPLICATION SIMULATIONS

Uncertainty updating in the description of heterogeneous materials

Multilevel Monte Carlo for two phase flow and transport in random heterogeneous porous media

Dynamic Data Driven Simulations in Stochastic Environments

Robust Ensemble Filtering With Improved Storm Surge Forecasting

Schwarz Preconditioner for the Stochastic Finite Element Method

Macroscopic Failure Analysis Based on a Random Field Representations Generated from Material Microstructures

A Posteriori Adaptive Low-Rank Approximation of Probabilistic Models

Lecture 1: Center for Uncertainty Quantification. Alexander Litvinenko. Computation of Karhunen-Loeve Expansion:

Bayesian parameter estimation in predictive engineering

Reactive Transport in Porous Media

Gaussian Process Approximations of Stochastic Differential Equations

Sequential Monte Carlo Samplers for Applications in High Dimensions

Adaptive Posterior Approximation within MCMC

Probabilistic Collocation Method for Uncertainty Analysis of Soil Infiltration in Flood Modelling

Uncertainty quantification of engineering systems using the multilevel Monte Carlo method

Bayesian Inverse problem, Data assimilation and Localization

Characterization of reservoir simulation models using a polynomial chaos-based ensemble Kalman filter

Ergodicity in data assimilation methods

Lagrangian Data Assimilation and Manifold Detection for a Point-Vortex Model. David Darmon, AMSC Kayo Ide, AOSC, IPST, CSCAMM, ESSIC

Homogenization Theory

Assessing the Value of Information from Inverse Modelling for Optimising Long-Term Oil Reservoir Performance

Measure-Theoretic parameter estimation and prediction for contaminant transport and coastal ocean modeling

Ensemble Kalman filter assimilation of transient groundwater flow data via stochastic moment equations

A Stochastic Collocation based. for Data Assimilation

Estimating functional uncertainty using polynomial chaos and adjoint equations

The Ensemble Kalman Filter:

The Ensemble Kalman Filter:

Introduction to Aspects of Multiscale Modeling as Applied to Porous Media

Dinesh Kumar, Mehrdad Raisee and Chris Lacor

Error Budgets: A Path from Uncertainty Quantification to Model Validation

A regularized least-squares method for sparse low-rank approximation of multivariate functions

Development of Stochastic Artificial Neural Networks for Hydrological Prediction

A Note on the Particle Filter with Posterior Gaussian Resampling

Consistent Downscaling of Seismic Inversions to Cornerpoint Flow Models SPE

Uncertainty Quantification and related areas - An overview

Dynamic System Identification using HDMR-Bayesian Technique

Uncertainty Quantification and hypocoercivity based sensitivity analysis for multiscale kinetic equations with random inputs.

Scalable Algorithms for Optimal Control of Systems Governed by PDEs Under Uncertainty

Probabilistic collocation and Lagrangian sampling for tracer transport in randomly heterogeneous porous media

UNCERTAINTY QUANTIFICATION IN LIQUID COMPOSITE MOULDING PROCESSES

Polynomial Chaos and Karhunen-Loeve Expansion

Kernel-based Approximation. Methods using MATLAB. Gregory Fasshauer. Interdisciplinary Mathematical Sciences. Michael McCourt.

Seminar: Data Assimilation

Uncertainty Quantification for multiscale kinetic equations with random inputs. Shi Jin. University of Wisconsin-Madison, USA

Winter 2019 Math 106 Topics in Applied Mathematics. Lecture 1: Introduction

Sparse polynomial chaos expansions in engineering applications

Employing Model Reduction for Uncertainty Visualization in the Context of CO 2 Storage Simulation

Non-Intrusive Solution of Stochastic and Parametric Equations

Stochastic Elastic-Plastic Finite Element Method for Performance Risk Simulations

STOCHASTIC SAMPLING METHODS

Multilevel Sequential 2 Monte Carlo for Bayesian Inverse Problems

Gaussian Filtering Strategies for Nonlinear Systems

Analysis and Computation of Hyperbolic PDEs with Random Data

Uncertainty Quantification for multiscale kinetic equations with high dimensional random inputs with sparse grids

Lagrangian Data Assimilation and Its Application to Geophysical Fluid Flows

Uncertainty analysis of large-scale systems using domain decomposition

PATTERN RECOGNITION AND MACHINE LEARNING CHAPTER 13: SEQUENTIAL DATA

A new Hierarchical Bayes approach to ensemble-variational data assimilation

Ensemble Kalman Filter

Parameter Selection Techniques and Surrogate Models

Autonomous Navigation for Flying Robots

arxiv: v2 [math.na] 8 Sep 2017

Fundamentals of Data Assimila1on

NEW ALGORITHMS FOR UNCERTAINTY QUANTIFICATION AND NONLINEAR ESTIMATION OF STOCHASTIC DYNAMICAL SYSTEMS. A Dissertation PARIKSHIT DUTTA

Stochastic collocation and mixed finite elements for flow in porous media

UNCERTAINTY ASSESSMENT USING STOCHASTIC REDUCED BASIS METHOD FOR FLOW IN POROUS MEDIA

Research Article Multiresolution Analysis for Stochastic Finite Element Problems with Wavelet-Based Karhunen-Loève Expansion

Machine Learning. Probability Basics. Marc Toussaint University of Stuttgart Summer 2014

Modelling Under Risk and Uncertainty

Transcription:

Optimisation under Uncertainty with Stochastic PDEs for the History Matching Problem in Reservoir Engineering Hermann G. Matthies Technische Universität Braunschweig wire@tu-bs.de http://www.wire.tu-bs.de

Overall Problem 2 Motivation: Computational tools for uncertainty quantification Situation: Geostatistical model with stochastic uncertainty and physical/ mathematical flow model Identification: History matching of well production data Goal: Predict well production and uncertainty Reservoir Long Term Goal: Optimal control of well production under uncertainty

Physical / Mathematical / Stochastic Model 3 Physical: multiphase flow in porous media (Darcy, Buckley-Leverett,...) Mathematical: coupled PDEs, conservation and transport/constitutive laws C t w + A(w) = g Stochastic: Prior geostatistic model of data and stochastic sources/sinks uncertainty. Simple Example: Darcy flow with stochastic coefficients c t u (κ u) = f Goal is the prediction of functional of future state u and its uncertainty. Explicit modelling and quantificaton of uncertainty. System identification under prior uncertainty and observed data. Predict future state and quantify uncertainty given posteriori model. Long term goal: Optimal Control of well yields.

Partly Open Problems 4 General Formulation: Coherent stochastic forward / inverse problem formulation, adequate type of identfication, geostatistical (non-gaussian) models. Forward Problem: Fast SPDE solvers based on stochastic Galerkin methods, much faster than Monte Carlo. Inverse Problem: Efficient solvers for nonlinear stochastic problem. Error Control: combined ( spatial & stochastic) error estimates and adaptivity. Uncertainty: Efficient quantification, high-dimensional integration. Software: Modular, independent, software management conformity. Long Term: Optimal learning / adaptive control.

Novel SPDE Techniques 5 General: Functional analytic formulation of stochastic problem, white noise analysis, Malliavin calculus. Solver: Stochastic Galerkin via polynomial chaos, tensor product solver, combination with fast deterministic solvers, low-rank approximation, Karhunen- Loève (POD, SVD), model reduction, efficient parallelisation. Error Control: Goal oriented, dual weighted residual (DWR), adjoint methods. Uncertainty: High-dimensional integration, sparse grids Smolyak methods. Other: Standard state of the art techniques, FEM/FVM & FDM spatial/ temporal discretisation, fast (multilevel) deterministic solvers. 4 3 2 1 1.5.5 x Mean.4.3.2.1 1.5.5 x 1.5.5 y 1 1 1.5.5 y 1 1 Standard Deviation.8.6.4.2 1.5.5 x 1.5.5 y 1 1 Pr{u(x) > 8}

System Identification Approaches 6 Approaches for system identification Maximum Likelihood: find most likely point. Bayesian like update of prior information based on observation stochastic posterior. Least Squares in equilibrium equation. Constitutive relation error CRE Separate certain (conservation: c u t + q = f) from uncertain assumptions (CR: q = κ u). Realisation Kalman filter EnKF Ensemble (via Monte Carlo) Kalman filtering.

Problem Sequence 7 Sequence of problems with increasing difficulty: linear diffusion: Darcy flow with stochastic conductivity stochastic non-gaussian velocity field. diffusion / advection: Solute transport (single phase) with stochastic velocity and (anomalous) dispersion. non-linear multi-phase flow, water air. Simulation NAPL multi-phase flow oil gas water...

Software Framework 8 Loose software coupling: Modular, independent parts. Software engineering: Components communicate via middleware. Tight numerical coupling: Integration of components into overall numerical scheme. Middleware: Component Template Library (CTL), allows distributed / coarse grain parallelisation. Optimiser Platon FEM ParaFep Eigenvalue (P)Arpack Communication Middleware CTL Component Template Library (MPI, PVM, TCP IP, etc...) Software KLE PCE StoFEL KLE StoFEL PCE Tensor Product Solver StoFEL PS

Related Collaboration / Work 9 Numerical / stochastic upscaling: Helmig IWS Uni Stuttgart; Cirpka ETH Zürich; Kulasiri Lincoln University, Christchurch. Stochastic material properties, error estimates: Ibrahimbegović, Ladevèze LMT ENS Cachan/ Paris; Villon UTC Compiegne; Wriggers IBNM Uni Hannover; Carstensen HU Berlin. Numerics for SPDEs: Schwab ETH Zürich; Tempone Florida State University; Babuška U of Texas, Austin; Burrage U of Queensland, Brisbane. System Identification: Faßbender. Coupled Problems: Ohayon CNAM, Paris. Software framework for history matching: Scandpower, Hamburg. Please visit poster early, I have to teach tomorrow morning.