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

Title: Multiscale solute transport upscaling for a three-dimensional hierarchical porous medium

 [1];  [1]
  1. Department of Geology and Geophysics, University of Wyoming, Laramie Wyoming USA
Publication Date:
Sponsoring Org.:
USDOE Office of Fossil Energy (FE)
OSTI Identifier:
Grant/Contract Number:
Resource Type:
Journal Article: Publisher's Accepted Manuscript
Journal Name:
Water Resources Research
Additional Journal Information:
Journal Volume: 51; Journal Issue: 3; Related Information: CHORUS Timestamp: 2017-10-23 17:26:09; Journal ID: ISSN 0043-1397
Wiley Blackwell (John Wiley & Sons)
Country of Publication:
United States

Citation Formats

Zhang, Mingkan, and Zhang, Ye. Multiscale solute transport upscaling for a three-dimensional hierarchical porous medium. United States: N. p., 2015. Web. doi:10.1002/2014WR016202.
Zhang, Mingkan, & Zhang, Ye. Multiscale solute transport upscaling for a three-dimensional hierarchical porous medium. United States. doi:10.1002/2014WR016202.
Zhang, Mingkan, and Zhang, Ye. 2015. "Multiscale solute transport upscaling for a three-dimensional hierarchical porous medium". United States. doi:10.1002/2014WR016202.
title = {Multiscale solute transport upscaling for a three-dimensional hierarchical porous medium},
author = {Zhang, Mingkan and Zhang, Ye},
abstractNote = {},
doi = {10.1002/2014WR016202},
journal = {Water Resources Research},
number = 3,
volume = 51,
place = {United States},
year = 2015,
month = 3

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record at 10.1002/2014WR016202

Citation Metrics:
Cited by: 1work
Citation information provided by
Web of Science

Save / Share:
  • A field experiment is reported which monitored the three-dimensional movement of cubic solute plumes through an unsaturated, loamy sand soil. The plumes were created with one of two methods, a two-dimensional flux application and an initial resident distribution. Soil coring was used to sample resident concentrations for the three solutes studied. The data were analyzed using the method of moments. In addition to the solute transport experiments, a detailed set of physical properties of the field was obtained by excavating three pits to a depth of 5.0 m and also by taking soil cores throughout the study area. This papermore » explains the experimental methodology, summarizes the relevant site characteristics. Mass balance varied between 78 and 138%. The field-averaged gravimetric water content and dry bulk density were used to accurately predict the mean vertical plume displacements. The plumes spread relatively little in the horizontal direction.« less
  • Solute transport in fractured porous media is typically 'non-Fickian'; that is, it is characterized by early breakthrough and long tailing and by nonlinear growth of the Green function-centered second moment. This behavior is due to the effects of (1) multirate diffusion occurring between the highly permeable fracture network and the low-permeability rock matrix, (2) a wide range of advection rates in the fractures and, possibly, the matrix as well, and (3) a range of path lengths. As a consequence, prediction of solute transport processes at the macroscale represents a formidable challenge. Classical dual-porosity (or mobile-immobile) approaches in conjunction with anmore » advection-dispersion equation and macroscopic dispersivity commonly fail to predict breakthrough of fractured porous media accurately. It was recently demonstrated that the continuous time random walk (CTRW) method can be used as a generalized upscaling approach. Here we extend this work and use results from high-resolution finite element-finite volume-based simulations of solute transport in an outcrop analogue of a naturally fractured reservoir to calibrate the CTRW method by extracting a distribution of retention times. This procedure allows us to predict breakthrough at other model locations accurately and to gain significant insight into the nature of the fracture-matrix interaction in naturally fractured porous reservoirs with geologically realistic fracture geometries.« less
  • 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 three-dimensional solute transport model has been developed to study detailed contaminant movements through large heterogeneous flow systems in porous media. The model is based upon a random walk particle method (RWPM) which can readily treat multidimensional advection and dispersion processes in saturated or unsaturated media in a computationally efficient manner. The transport simulations are used to examine the large time and spatial effects of the variable flow field on developing solute plumes, and, in particular, to investigate the nature of the large-scale dispersive behavior. Numerical transport experiments were conducted using single realizations of random hydraulic conductivity fields with threemore » different degrees of heterogeneity. Experiments with different source locations were used to investigate preasymptotic and nonergodic effects that would appear as differences in plume evolution among the experiments. Analyses of the simulations indicate the spatial moments of the particle distributions in statistically isotropic saturated media compare favorably with stochastic theory predictions in terms of longitudinal advection and mixing, but differ markedly from predictions of transverse mixing. The simulations also demonstrate that significant nonergodic effects occur, as reflected in strong differences in the second moment evolution curves among the individual experiments and as predicted from the ensemble stochastic theory.« less
  • Solute plumes were created in an unsaturated field soil with either flux application or by leaching an initial resident distribution. The spatial variance of the plumes initially increased with time between the soil surface and a depth of 2.5 m, within which the soil was a nearly structureless loamy sand. Below this depth, the plumes were observed to compress in the vertical direction as they moved into, and through, a region of subangular blocky structure and loam texture (between 2.5 and 4.0 m depth). As the solute moved below the layer of fine texture, the plume variance again increased withmore » time. Using a transformed advection-dispersion equation description, two constant, field-averaged transport coefficients, V* and D*{sub zz}, were determined in a scaled coordinate system from the moment equations. These two constant parameters were then used to predict the observed local, or plot scale, transport. Results indicate that the two constant parameters describe transport reasonably well as each plot site and over all sampling depths.« less