Publications

The authors consider the problem of defining the fracture permeability tensor for each grid lock in a rock mass from maps of natural fractures. For this purpose they implement a statistical model of cracked rock due to M. Oda [1985], where the permeability tensor is related to the crack geometry via a volume average of the contribution from each crack in the population. In this model tectonic stress is implicitly coupled to fluid flow through an assumed relationship between crack aperture and normal stress across the crack. The authors have included the following enhancements to the basic model: (1) a realistic model of crack closure under stress has been added along with the provision to apply tectonic stresses to the fracture system in any orientation, the application of stress results in fracture closure and consequently a reduction in permeability; (2) the fracture permeability can be superimposed onto an arbitrary anisotropic matrix permeability; (3) the fracture surfaces are allowed to slide under the application of shear stress, causing fractures to dilate and result in a permeability increase. Through an example, the authors demonstrate that significant changes in permeability magnitudes and orientations are possible when tectonic stress is applied to a fracture system.

An experimental study of electrokinetic effects (streaming potential) in earth materials was undertaken. The objective was to evaluate the measurement of electrokinetic effects as a method of monitoring and predicting the movement of groundwater, contaminant plumes, and other fluids in the subsurface. The laboratory experiments verified that the electrokinetic effects in earth materials are prominent, repeatable, and can be described well to first order by a pair of coupled differential equations.

The evaluation of the stability of the openings for the Exploratory Studies Facility and a potential repository for high-level nuclear waste at Yucca Mountain, Nevada will require computer codes capable of predicting slip on rock joints resulting from changes in thermal stresses. The geometrical method of analysis of moire fringe analysis was used to evaluate the magnitude and extent of frictional sliding in a layered polycarbonate rock mass model containing a circular hole. Slips were observed in confined zones around the hole and micron resolutions were obtained. Unpredicted and uncontrolled uniform slip of several interfaces in the model were observed giving considerable uncertainty in the boundary conditions of the model, perhaps making detailed comparison with numerical models impossible.

Prediction of the deformation behavior of large engineering structures in jointed rock under a specified loading history requires the extensive use of numerical simulation. For example, the evaluation of the stability of the openings for the Exploratory Studies Facility and a potential repository for high-level nuclear waste at Yucca Mountain, Nevada will require computer codes capable of predicting slip on rock joints resulting from changes in thermal stresses. The testing and ultimate validation of these complex finite element computer codes is an important step in their development before their use as a design tool for an engineering structure or for the study of some other practical problem. While field tests may be ultimately necessary, the authors propose a different and more thorough approach where early tests are done on a bench scale with easily characterized materials and geometries. For these bench-scale tests, the basic approach is to construct a laboratory specimen with a known geometry from an easily characterized material. Digital video imaging combined with the geometric moire fringe method of strain analysis is used to measure and derive the displacements on the sample under load. Here the authors present the method of acquiring and analyzing the moire data and give an analysis of its problems and benefits.