  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.
