skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A two-stage adaptive stochastic collocation method on nested sparse grids for multiphase flow in randomly heterogeneous porous media

Abstract

A new computational method is proposed for efficient uncertainty quantification of multiphase flow in porous media with stochastic permeability. For pressure estimation, it combines the dimension-adaptive stochastic collocation method on Smolyak sparse grids and the Kronrod–Patterson–Hermite nested quadrature formulas. For saturation estimation, an additional stage is developed, in which the pressure and velocity samples are first generated by the sparse grid interpolation and then substituted into the transport equation to solve for the saturation samples, to address the low regularity problem of the saturation. Numerical examples are presented for multiphase flow with stochastic permeability fields to demonstrate accuracy and efficiency of the proposed two-stage adaptive stochastic collocation method on nested sparse grids.

Authors:
 [1];  [1];  [2]
  1. ERE & BIC-ESAT, College of Engineering, Peking University, Beijing (China)
  2. Department of Energy Resources Engineering, Stanford University, Stanford, CA (United States)
Publication Date:
OSTI Identifier:
22622250
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Computational Physics; Journal Volume: 330; Other Information: Copyright (c) 2016 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND SUPERFLUIDITY; ACCURACY; COMPUTERIZED SIMULATION; EFFICIENCY; INTERPOLATION; MULTIPHASE FLOW; PERMEABILITY; POROUS MATERIALS; QUADRATURES; RANDOMNESS; SATURATION; STOCHASTIC PROCESSES; VELOCITY

Citation Formats

Liao, Qinzhuo, E-mail: liaoqz@pku.edu.cn, Zhang, Dongxiao, and Tchelepi, Hamdi. A two-stage adaptive stochastic collocation method on nested sparse grids for multiphase flow in randomly heterogeneous porous media. United States: N. p., 2017. Web. doi:10.1016/J.JCP.2016.10.061.
Liao, Qinzhuo, E-mail: liaoqz@pku.edu.cn, Zhang, Dongxiao, & Tchelepi, Hamdi. A two-stage adaptive stochastic collocation method on nested sparse grids for multiphase flow in randomly heterogeneous porous media. United States. doi:10.1016/J.JCP.2016.10.061.
Liao, Qinzhuo, E-mail: liaoqz@pku.edu.cn, Zhang, Dongxiao, and Tchelepi, Hamdi. Wed . "A two-stage adaptive stochastic collocation method on nested sparse grids for multiphase flow in randomly heterogeneous porous media". United States. doi:10.1016/J.JCP.2016.10.061.
@article{osti_22622250,
title = {A two-stage adaptive stochastic collocation method on nested sparse grids for multiphase flow in randomly heterogeneous porous media},
author = {Liao, Qinzhuo, E-mail: liaoqz@pku.edu.cn and Zhang, Dongxiao and Tchelepi, Hamdi},
abstractNote = {A new computational method is proposed for efficient uncertainty quantification of multiphase flow in porous media with stochastic permeability. For pressure estimation, it combines the dimension-adaptive stochastic collocation method on Smolyak sparse grids and the Kronrod–Patterson–Hermite nested quadrature formulas. For saturation estimation, an additional stage is developed, in which the pressure and velocity samples are first generated by the sparse grid interpolation and then substituted into the transport equation to solve for the saturation samples, to address the low regularity problem of the saturation. Numerical examples are presented for multiphase flow with stochastic permeability fields to demonstrate accuracy and efficiency of the proposed two-stage adaptive stochastic collocation method on nested sparse grids.},
doi = {10.1016/J.JCP.2016.10.061},
journal = {Journal of Computational Physics},
number = ,
volume = 330,
place = {United States},
year = {Wed Feb 01 00:00:00 EST 2017},
month = {Wed Feb 01 00:00:00 EST 2017}
}
  • In this study, a probabilistic collocation method (PCM) on sparse grids was used to solve stochastic equations describing flow and transport in three-dimensional in saturated, randomly heterogeneous porous media. Karhunen-Lo\`{e}ve (KL) decomposition was used to represent the three-dimensional log hydraulic conductivitymore » $$Y=\ln K_s$$. The hydraulic head $h$ and average pore-velocity $$\bf v$$ were obtained by solving the three-dimensional continuity equation coupled with Darcy's law with random hydraulic conductivity field. The concentration was computed by solving a three-dimensional stochastic advection-dispersion equation with stochastic average pore-velocity $$\bf v$$ computed from Darcy's law. PCM is an extension of the generalized polynomial chaos (gPC) that couples gPC with probabilistic collocation. By using the sparse grid points, PCM can handle a random process with large number of random dimensions, with relatively lower computational cost, compared to full tensor products. Monte Carlo (MC) simulations have also been conducted to verify accuracy of the PCM. By comparing the MC and PCM results for mean and standard deviation of concentration, it is evident that the PCM approach is computational more efficient than Monte Carlo simulations. Unlike the conventional moment-equation approach, there is no limitation on the amplitude of random perturbation in PCM. Furthermore, PCM on sparse grids can efficiently simulate solute transport in randomly heterogeneous porous media with large variances.« less
  • A computational methodology is developed to efficiently perform uncertainty quantification for fluid transport in porous media in the presence of both stochastic permeability and multiple scales. In order to capture the small scale heterogeneity, a new mixed multiscale finite element method is developed within the framework of the heterogeneous multiscale method (HMM) in the spatial domain. This new method ensures both local and global mass conservation. Starting from a specified covariance function, the stochastic log-permeability is discretized in the stochastic space using a truncated Karhunen-Loeve expansion with several random variables. Due to the small correlation length of the covariance function,more » this often results in a high stochastic dimensionality. Therefore, a newly developed adaptive high dimensional stochastic model representation technique (HDMR) is used in the stochastic space. This results in a set of low stochastic dimensional subproblems which are efficiently solved using the adaptive sparse grid collocation method (ASGC). Numerical examples are presented for both deterministic and stochastic permeability to show the accuracy and efficiency of the developed stochastic multiscale method.« less
  • In recent years, there has been intense interest in understanding various physical phenomena in random heterogeneous media. Any accurate description/simulation of a process in such media has to satisfactorily account for the twin issues of randomness as well as the multilength scale variations in the material properties. An accurate model of the material property variation in the system is an important prerequisite towards complete characterization of the system response. We propose a general methodology to construct a data-driven, reduced-order model to describe property variations in realistic heterogeneous media. This reduced-order model then serves as the input to the stochastic partialmore » differential equation describing thermal diffusion through random heterogeneous media. A decoupled scheme is used to tackle the problems of stochasticity and multilength scale variations in properties. A sparse-grid collocation strategy is utilized to reduce the solution of the stochastic partial differential equation to a set of deterministic problems. A variational multiscale method with explicit subgrid modeling is used to solve these deterministic problems. An illustrative example using experimental data is provided to showcase the effectiveness of the proposed methodology.« less