Abstract

This paper investigates analytical solutions of stochastic Darcy flow in randomly heterogeneous porous media. We focus on infinite series solutions of the steady-state equations in the case of continuous porous media whose saturated log-conductivity (lnK) is a gaussian random field. The standard deviation of lnK is denoted ‘σ’. The solution method is based on a Taylor series expansion in terms of parameter σ, around the value σ=0, of the hydraulic head (H) and gradient (J). The head solution H is expressed, for any spatial dimension, as an infinite hierarchy of Green's function integrals, and the hydraulic gradient J is given by a linear first-order recursion involving a stochastic integral operator. The convergence of the ‘σ-expansion’ solution is not guaranteed a priori. In one dimension, however, we prove convergence by solving explicitly the hierarchical sequence of equations to all orders. An ‘infinite-order stochastic solution is obtained in the form of a σ-power series that converges for any finite value of σ. It is pointed out that other expansion methods based on K rather than lnK yield divergent series. The infinite-order solution depends on the integration method and the boundary conditions imposed on individual order equations. The most flexible and general method is that based on Laplacian Green's functions and boundary integrals. Imposing zero head conditions for all orders greater than one yields meaningful far-field gradient conditions. The whole approach can serve as a basis for treatment of higher-dimensional problems.

Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call