Search

A multiporosity extension of classical double and triple porosity fractured rock flow models for slightly compressible fluids is presented. The multiporosity model is an adaptation of the multirate solute transport model of Haggerty and Gorelick (1995) to viscous flow in fractured rock reservoirs. It is a generalization of both pseudo-steady-state and transient interporosity flow double porosity models. The model includes a fracture continuum and an overlapping distribution of multiple rock matrix continua, whose fracture-matrix exchange coefficients are specified through a discrete probability mass function. Semi-analytical cylindrically symmetric solutions to the multiporosity mathematical model are developed using the Laplace transform to illustrate its behavior. The multiporosity model presented here is conceptually simple, yet flexible enough to simulate common conceptualizations of double and triple porosity flow. This combination of generality and simplicity makes the multiporosity model a good choice for flow in low-permeability fractured rocks., 5 figures, 3 tables

A semianalytical solution is presented for transient streaming potentials associated with flow to a pumping well in an unconfined aquifer, taking into account the effect of flow in the unsaturated zone above the water table. Flow in the unsaturated zone is modeled with a linearized form of Richards' equation using an exponential model for soil moisture retention and unsaturated hydraulic conductivity. Archie's law is invoked for unsaturated electrical conductivity. The unsaturated electrokinetic coupling coefficient is modeled with a decaying exponential, where the maximum value is at and below the water table. The coupled flow and electrokinetic problem is solved using Laplace and Hankel transforms. The results of the model predicted behavior are presented and compared to that observed in laboratory simulations of pump­ ing tests. The early time polarity reversal predicted the model is observable in the experiments. Other non- monotonic streaming potential behaviors predicted by the model are also evident in experimental measurements. The model is used to estimate hydraulic parameters from SP data and these compare well to those obtained from drawdown data. For example, a hydraulic conductivity of 3.6 3 10 24 m/s is obtained from SP data compared to 3.4 3 10 24 m/s from drawdown data.

Model applicability to core-scale solute transport is evaluated using breakthrough data from column experiments conducted with conservative tracers tritium (H-3) and sodium-22, and the retarding solute uranium-232. The three models considered are single-porosity, double-porosity with single-rate mobile-immobile mass-exchange, and the multirate model, which is a deterministic model that admits the statistics of a random mobile-immobile mass-exchange rate coefficient. The experiments were conducted on intact Culebra Dolomite core samples. Previously, data were analyzed using single- and double-porosity models although the Culebra Dolomite is known to possess multiple types and scales of porosity, and to exhibit multirate mobile-immobile-domain mass transfer characteristics at field scale. The data are reanalyzed here and null-space Monte Carlo analysis is used to facilitate objective model selection. Prediction (or residual) bias is adopted as a measure of the model structural error. The analysis clearly shows single- and double-porosity models are structurally deficient, yielding late-time residual bias that grows with time. On the other hand, the multirate model yields unbiased predictions consistent with the late-time -5/2 slope diagnostic of multirate mass transfer. The analysis indicates the multirate model is better suited to describing core-scale solute breakthrough in the Culebra Dolomite than the other two models.

[1] Riparian zones in semiarid regions often exhibit high rates of evapotranspiration (ET) in spite of low-soil moisture content due to the presence of phreatophytic vegetation that is able to withdraw water from shallow aquifers. This work seeks to better define the relationship between ET, the saturated zone and the river boundary by comparing observed water table drawdown records to analytically modeled drawdown in fully penetrating wells of an unconfined aquifer in response to daily ET flux. ET at the Boise Hydrogeophysical Research Site (BHRS), a riparian zone in a temperate, semiarid environment, is calculated using a radiation-based method to provide ET values at four different wells with different vegetation densities. Analytically modeled drawdown response to ET forcing shows that drawdown magnitude increases with increasing distance from the river edge even as the surficial ET forcing remains constant. This behavior is also observed in well hydrographs and shows the buffering effect that flow from the river has on drawdown in fully penetrating riparian wells in highly permeable, unconfined aquifers. Relative contributions of river water to aquifer storage are calculated for ET-induced diurnal fluctuations of the water table at increasing distances from the river boundary. Failure to account for these spatial differences in drawdown related to the river source may explain some errors associated with estimating ET from well hydrographs alone.