Abstract:Limiting the injection rate to restrict the pressure below a threshold at a critical location can be an important goal of simulations that model the subsurface pressure between injection and extraction wells. The pressure is approximated by the solution of Darcy's partial differential equation (PDE) for a given permeability field. The subsurface permeability is modeled as a random field since it is known only up to statistical properties. This induces uncertainty in the computed pressure. Solving the PDE for an ensemble of random permeability simulations enables estimating a probability distribution for the pressure at the critical location. These simulations are computationally expensive, and practitioners often need rapid online guidance for real-time pressure management. An ensemble of numerical PDE solutions is used to construct a Gaussian process regression model that can quickly predict the pressure at the critical location as a function of the extraction rate and permeability realization. Our first novel contribution is to identify a sampling methodology for the random environment and matching kernel technology for which fitting the Gaussian process regression model scales as O(n log n) instead of the typical O(n^3) rate in the number of samples n used to fit the surrogate. The surrogate model allows almost instantaneous predictions for the pressure at the critical location as a function of the extraction rate and permeability realization. Our second contribution is a novel algorithm to calibrate the uncertainty in the surrogate model to the discrepancy between the true pressure solution of Darcy's equation and the numerical solution. Although our method is derived for building a surrogate for the solution of Darcy's equation with a random permeability field, the framework broadly applies to solutions of other PDE with random coefficients.

Computing first passage times for Markov-modulated fluid models using numerical PDE problem solvers

A Computational Approach to First Passage Problems of Reflected Hyper-Exponential Jump Diffusion Processes

Efficient computation of first passage times in Kou's jump-diffusion model

Data-driven Identification for Approximate Analytical Solution of First-Passage Problem

Markov chain Monte Carlo with Gaussian processes for fast parameter estimation and uncertainty quantification in a 1D fluid‐dynamics model of the pulmonary circulation

Computing probability density of the first passage time for state transition in stochastic dynamical systems driven by Brownian motions: A singular integral method

Analytic and Monte Carlo Approximations to the Distribution of the First Passage Time of the Drifted Diffusion with Stochastic Resetting and Mixed Boundary Conditions

Mean First Passage Times and Eyring-Kramers formula for Fluctuating Hydrodynamics

Scaling up Probabilistic PDE Simulators with Structured Volumetric Information

First-Passage Kinetic Monte Carlo method

A convergent discretisation method for transition path theory for diffusion processes

Transient Analysis of Fluid Flow Models via Matrix Decomposition

Studying First Passage Problems using Neural Networks: A Case Study in the Slit-Well Microfluidic Device

Numerically Modelling Stochastic Lie Transport in Fluid Dynamics

Benchmarking Autoregressive Conditional Diffusion Models for Turbulent Flow Simulation

Numerical Approximations and Convergence Analysis of Piecewise Diffusion Markov Processes, with Application to Glioma Cell Migration

The Computation of the Mean First Passage Times for Markov Chains

Numerically Modeling Stochastic Lie Transport in Fluid Dynamics

Computationally Efficient and Error Aware Surrogate Construction for Numerical Solutions of Subsurface Flow Through Porous Media

Analytical Solution of the Time-Dependent Microfluidic Poiseuille Flow in Rectangular Channel Cross-Sections and its Numerical Implementation in Microsoft Excel