A Bayesian model selection analysis of equilibrium and nonequilibrium models for multiphase flow in porous media. (March 2017)

• Record Type: Journal Article
• Authors: Ren, Guotong; Rafiee, Javad; Aryana, Saman A.; Younis, Rami M.
• Is Part Of: International journal of multiphase flow. Volume 89(2017)
• Page Start: 313
• Page End: 320
• Publication Date: 2017-03
• Subjects: Nonequilibrium two-phase flow -- Buckley–Leverett -- Model selection -- Core scale drainage experiments; Multiphase flow -- Periodicals; Écoulement polyphasique -- Périodiques; Multiphase flow; Periodicals; 620.1064
• DOI: 10.1016/j.ijmultiphaseflow.2016.11.006
• Languages: English
• ISSNs: 0301-9322

Abstract

Highlights: Assumption of instantaneous local equilibrium is inadequate for unstable processes. Candidate models are calibrated to core-scale drainage experimental data. The maximum likelihood estimate is calculated for candidate models. The classical formulation appears to be adequate for stable displacement processes. Nonequilibrium effects must be taken into account in unstable displacement processes. Abstract: The classical constitutive relations for multiphase flows in porous media assume instantaneous and local phase-equilibrium. Several alternative nonequilibrium/dynamic constitutive relations have been proposed in the literature including the works of Barenblatt, and Hassanizadeh and Gray. This work applies a Bayesian model selection framework in order to examine the relative efficacy of these three models to represent experimental observations. Experimental observations of multiphase displacement processes in natural porous media are often sparse and indirect, leading to considerable uncertainty in control conditions. Data from three core-scale drainage experiments are considered. Gaussian prior probability models are assumed for key multiphase flow parameters and measurements. Accurate numerical simulation approximations using the three constitutive relation models are implemented. The model selection analysis comprises a data-assimilation stage that calibrates the assumed model to the data while quantifying uncertainty. The second stage is the computation of the maximum likelihood estimate and its application to compute the Bayesian Information Criterion. It is observed that Barenblatt's nonequilibrium model is more likely to match data from unstable displacements that involve higher viscosity ratios of the invading phase to the resident fluid. At the lowest viscosity ratio, there is no delineation between the goodness of fit obtained using the classical model and the model proposed by Hassanizadeh and Gray, and both outperform Barenblatt's nonequilibrium model.