Random-Walk Simulation of Transport in Heterogeneous Porous Media: Local Mass-Conservation Problem and Implementation Methods

The random-walk method for simulating solute transport in porous media is typically based on the assumption that the velocity and velocity-dependent dispersion tensor vary smoothly in space. However, in cases where sharp interfaces separate materials with contrasting hydraulic properties, these quantities may be discontinuous. Normally, velocities are interpolated to arbitrary particle locations when finite difference or finite element methods are used to solve the flow equation. The use of interpolation schemes that preserve discontinuities in velocity at material contacts can result in a random-walk model that does not locally conserve mass unless a correction is applied at these contacts. Test simulations of random-walk particle tracking with and without special treatment of material contacts demonstrate the problem. Techniques for resolving the problem, including interpolation schemes and a reflection principle, are reviewed and tested. Results from simulations of transport in porous media with discontinuities in the dispersion tensor show which methods satisfy continuity. Simulations of transport in two-dimensional heterogeneous porous media demonstrate the potentially significant effect of using a nonconservative model to compute spatial moments and breakthrough of a solute plume.

[1]  D. W. Pollock.  Semianalytical Computation of Path Lines for Finite‐Difference Models , 1988 .

[2]  L. Rosenhead.  Conduction of Heat in Solids , 1947, Nature.

[3]  Andrew F. B. Tompson,et al.  Numerical simulation of chemical migration in physically and chemically heterogeneous porous media , 1993 .

[4]  G. J. M. Uffink,et al.  A random walk method for the simulation of macrodispersion in a stratified aquifer , 1985 .

[5]  E. Sudicky,et al.  An Advection‐Diffusion Concept for Solute Transport in Heterogeneous Unconsolidated Geological Deposits , 1984 .

[6]  Steven G. Smith,et al.  Use of high performance computing to examine the effectiveness of aquifer remediation , 1994 .

[7]  T. A. Prickett,et al.  A "random-walk" solute transport model for selected groundwater quality evaluations , 1981 .

[8]  W. Kinzelbach,et al.  The Random Walk Method in Pollutant Transport Simulation , 1988 .

[9]  C. Axness,et al.  Three‐dimensional stochastic analysis of macrodispersion in aquifers , 1983 .

[10]  Daniel J. Goode,et al.  Particle velocity interpolation in block‐centered finite difference groundwater flow models , 1990 .

[11]  C. R. Cole,et al.  Multicomponent mass transport model: theory and numerical implementation (discrete-parcel-random-walk version) , 1977 .

[12]  K. Vahala.  Handbook of stochastic methods for physics, chemistry and the natural sciences , 1986, IEEE Journal of Quantum Electronics.

[13]  A.F.B. Tompson,et al.  Numerical simulation of solute transport in randomly heterogeneous porous media: Motivation, model development, and application , 1987 .

[14]  Arlen W. Harbaugh,et al.  A modular three-dimensional finite-difference ground-water flow model , 1984 .

[15]  G. Matheron,et al.  Is transport in porous media always diffusive? A counterexample , 1980 .

[16]  Peter E. Kloeden,et al.  A survey of numerical methods for stochastic differential equations , 1989 .

[17]  Amvrossios C. Bagtzoglou,et al.  Projection functions for particle‐grid methods , 1992 .

[18]  D. E. Dougherty,et al.  Particle-grid methods for reacting flows in porous media with application to Fisher's equation , 1992 .

[19]  Graham E. Fogg,et al.  STOCHASTIC SIMULATION OF INTERWELL-SCALE HETEROGENEITY FOR IMPROVED PREDICTION OF SWEEP EFFICIENCY IN A CARBONATE RESERVOIR , 1991 .

[20]  R. H. Brown,et al.  Theory of aquifer tests , 1962 .

[21]  J. Bredehoeft,et al.  Computer model of two-dimensional solute transport and dispersion in ground water , 1978 .

[22]  Andrew F. B. Tompson,et al.  Numerical simulation of solute transport in three-dimensional, randomly heterogeneous porous media , 1990 .

[23]  W. Kinzelbach,et al.  Continuous Groundwater Velocity Fields and Path Lines in Linear, Bilinear, and Trilinear Finite Elements , 1992 .

[24]  J. Westwater,et al.  The Mathematics of Diffusion. , 1957 .

[25]  John L. Wilson,et al.  Efficient and accurate front tracking for two‐dimensional groundwater flow models , 1991 .

[26]  Robert W. Gillham,et al.  Experimental Investigation of Solute Transport in Stratified Porous Media: 2. The Reactive Case , 1985 .

[27]  J. Bear.  Dynamics of Fluids in Porous Media , 1975 .

[28]  Gerard J. M. Uffink,et al.  Modeling of Solute Transport with the Random Walk Method , 1988 .

[29]  Feller William,et al.  An Introduction To Probability Theory And Its Applications , 1950 .

[30]  E. C. Childs.  Dynamics of fluids in Porous Media , 1973 .