skip to main content


Title: Method of Distributions for Two‐Phase Flow in Heterogeneous Porous Media
Abstract

Multiscale heterogeneity and insufficient characterization data for a specific subsurface formation of interest render predictions of multi‐phase fluid flow in geologic formations highly uncertain. Quantification of the uncertainty propagation from the geomodel to the fluid‐flow response is typically done within a probabilistic framework. This task is computationally demanding due to, for example, the slow convergence of Monte Carlo simulations (MCS), especially when computing the tails of a distribution that are necessary for risk assessment and decision‐making under uncertainty. The frozen streamlines method (FROST) accelerates probabilistic predictions of immiscible two‐phase fluid flow problems; however, FROST relies on MCS to compute the travel‐time distribution, which is then used to perform the transport (phase saturation) computations. To alleviate this computational bottleneck, we replace MCS with a deterministic equation for the cumulative distribution function (CDF) of travel time. The resulting CDF‐FROST approach yields the CDF of the saturation field without resorting to sampling‐based strategies. Our numerical experiments demonstrate the high accuracy of CDF‐FROST in computing the CDFs of both saturation and travel time. For the same accuracy, it is about 5 and 10 times faster than FROST and MCS, respectively.

 
more » « less
PAR ID:
10444059
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Water Resources Research
Volume:
58
Issue:
12
ISSN:
0043-1397
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Two-phase porous media flow is important in many applications from drug delivery to groundwater diffusion and oil recovery and is of particular interest to biomedical diagnostic test developers using cellulose and nitrocellulose membranes with limited fluid sample volumes. This work presents a new two-phase porous media flow model based on the incompressible Navier−Stokes equation. The model aims to address the limitations of existing methods by incorporating a partial saturation distribution in porous media to account for limited fluid volumes. The basic parameters of the model are the pore size distribution and the contact angle. To validate the model, we solved five analytical solutions and compared them to corresponding experimental data. The experimentally measured penetration length data agreed with the model predictions, demonstrating model accuracy. Our findings suggest that this new two-phase porous media flow model can provide a valuable tool for researchers developing fluidic assays in paper and other porous media. 
    more » « less
  2. Abstract

    Flow through partially frozen pores in granular media containing ice or gas hydrate plays an essential role in diverse phenomena including methane migration and frost heave. As freezing progresses, the frozen phase grows in the pore space and constricts flow paths so that the permeability decreases. Previous works have measured the relationship between permeability and volumetric fraction of the frozen phase, and various correlations have been proposed to predict permeability change in hydrology and the oil industry. However, predictions from different formulae can differ by orders of magnitude, causing great uncertainty in modeling results. We present a floating random walk method to approximate the porous flow field and estimate the effective permeability in isotropic granular media with specified particle size distributions, without solving for the entire flow field in the pore space. In packed spherical particles, the method compares favorably with the Kozeny‐Carman formula. We further extend this method with a probabilistic interpretation of the volumetric fraction of the frozen phase, simulate the effect of freezing in irregular pores, and predict the evolution of permeability. Employing no adjustable parameters, our results can provide insight into the coupling between phase transitions and permeability change, which plays important roles in hydrate formation and dissociation, as well as in the thawing and freezing of permafrost and ice‐bed coupling beneath glaciers.

     
    more » « less
  3. The remarkable success of the Transformer model in Natural Language Processing (NLP) is increasingly capturing the attention of vision researchers in contemporary times. The Vision Transformer (ViT) model effectively models long-range dependencies while utilizing a self-attention mechanism by converting image information into meaningful representations. Moreover, the parallelism property of ViT ensures better scalability and model generalization compared to Recurrent Neural Networks (RNN). However, developing robust ViT models for high-risk vision applications, such as self-driving cars, is critical. Deterministic ViT models are susceptible to noise and adversarial attacks and incapable of yielding a level of confidence in output predictions. Quantifying the confidence (or uncertainty) level in the decision is highly important in such real-world applications. In this work, we introduce a probabilistic framework for ViT to quantify the level of uncertainty in the model's decision. We approximate the posterior distribution of network parameters using variational inference. While progressing through non-linear layers, the first-order Taylor approximation was deployed. The developed framework propagates the mean and covariance of the posterior distribution through layers of the probabilistic ViT model and quantifies uncertainty at the output predictions. Quantifying uncertainty aids in providing warning signals to real-world applications in case of noisy situations. Experimental results from extensive simulation conducted on numerous benchmark datasets (e.g., MNIST and Fashion-MNIST) for image classification tasks exhibit 1) higher accuracy of proposed probabilistic ViT under noise or adversarial attacks compared to the deterministic ViT. 2) Self-evaluation through uncertainty becomes notably pronounced as noise levels escalate. Simulations were conducted at the Texas Advanced Computing Center (TACC) on the Lonestar6 supercomputer node. With the help of this vital resource, we completed all the experiments within a reasonable period. 
    more » « less
  4. Abstract

    We study the dependence of the transverse mass distribution of charged leptons and the missing energy on parton distributions (PDFs) adapted toWboson mass measurements at the CDF and ATLAS experiments. We compare the shape variations of the distribution induced by different PDFs and find that the spread of predictions from different PDF sets can be significantly larger than the PDF uncertainty predicted by a specific PDF set. We suggest analyzing the experimental data using up-to-date PDFs to gain a better understanding of the PDF uncertainties inWboson mass measurements. We also perform a series of Lagrange multiplier scans to identify the constraints on the transverse mass distribution imposed by individual data sets in the CT18 global analysis. In the case of the CDF measurement, the distribution is mostly sensitive tod-quark PDFs in the intermediatexregion, which are largely constrained by DIS and Drell-Yan data on deuteron targets and Tevatron lepton charge asymmetry data.

     
    more » « less
  5. Learning multi-agent dynamics is a core AI problem with broad applications in robotics and autonomous driving. While most existing works focus on deterministic prediction, producing probabilistic forecasts to quantify uncertainty and assess risks is critical for downstream decision-making tasks such as motion planning and collision avoidance. Multi-agent dynamics often contains internal symmetry. By leveraging symmetry, specifically rotation equivariance, we can improve not only the prediction accuracy but also uncertainty calibration. We introduce Energy Score, a proper scoring rule, to evaluate probabilistic predictions. We propose a novel deep dynamics model, Probabilistic Equivariant Continuous COnvolution (PECCO) for probabilistic prediction of multi-agent trajectories. PECCO extends equivariant continuous convolution to model the joint velocity distribution of multiple agents. It uses dynamics integration to propagate the uncertainty from velocity to position. On both synthetic and real-world datasets, PECCO shows significant improvements in accuracy and calibration compared to non-equivariant baselines. 
    more » « less