ph: 805 886 9856
john
1990-97: Research Associate, Department of Physics, University of Colorado, Boulder, on DOD, USAF and State of Nevada contracts for mathematical and numerical modeling of physical processes with emphasis on multi-component fluid flow in complex environments. Required project management, proposal preparation, results presentation at meetings with submission of regular reports and publication of papers in refereed journals. Researched, published and won awards for papers that include modeling of nonlinear processes such as flows through mineralized and fractured materials, seismic noise and prediction of subsurface fluid motions and atmospheric and ionospheric disturbances due to impulsive forces. Advised and worked closely with graduate students.
1991-95: Consultant on State of Nevada contract through TRAC Inc. for mathematical and numerical modeling of geophysical processes relevant to the safety of the proposed nuclear waste repository at Yucca Mountain, Nevada, for earthquake and flooding hazard. Researched, published and presented papers on the modeling of earthquakes and the interaction of subsurface fluid motions with stress and strain changes.
Davies, J.B. and C.B. Archambeau, 1997: "Geohydrological Models and Earthquake Effects at Yucca Mountain, Nevada", Environmental Geology, 32, 23-35; and presented at IUGG 1995.
ABSTRACT
Yucca Mountain, the proposed site for the high-level nuclear waste repository, is located just south of where the present water table begins a sharp rise in elevation. This large hydraulic gradient is a regional feature that extends for over 100 km. Yucca Mountain and its vicinity are underlain by faulted and fractured tuffs with hydraulic conductivities controlled by flow through the fractures. Close to and parallel with the region of large hydraulic gradient, and surrounding the core of the Timber Mountain Caldera, there is a 10- to 20-km-wide zone containing few faults and thus, most likely, few open fractures. Consequently, this zone should have a relatively low hydraulic conductivity, and this inference is supported by the available conductivity measurements in wells near the large hydraulic gradient. Also, slug injection tests indicate significantly higher pressures for fracture opening in wells located near the large hydraulic gradient compared to the opening pressures in wells further to the south, hence implying that lower extensional stresses prevail to the north with consequently fewer open fractures there. Analytical and numerical modeling shows that such a boundary between media of high and low conductivity can produce the observed, large hydraulic gradient, with the high conductivity medium having a lower elevation of the water table. Further, as fractures can close due to tectonic activity, the conductivity of the Yucca Mountain tuffs can be reduced to a value near that for the hydraulic barrier due to strain release by a moderate earthquake. Under these conditions, simulations show that the elevation of the steady-state water table could rise between 150 and 250 m at the repository site. This elevation rise is due to the projected shift in the location of the large hydraulic gradient to the south in response to a moderate earthquake, near magnitude 6, along one of the major normal faults adjacent to Yucca Mountain. As the proposed repository would only be 200–400 m above the present water table, this predicted rise in the water table indicates a potential hazard involving water intrusion.
Davies, J.B. and C.B. Archambeau, 1997: "Analysis of High-Pressure Fluid Flow in Fractures, with application to Yucca Mountain, Nevada, slug test data", Tectonophysics, 277, 83-98; invited paper after presentation at IUGG 1995.
ABSTRACT
Increases in water pressure, if large enough, can be sufficient to open existing fractures in a fractured rock-water system. Such large water pressure increases are expected during earthquake activity and are also produced artificially in high-pressure hydraulic fracturing and slug injection tests. In an in-situ slug injection test, the slug of water is initially released into a packed-off interval of the borehole where the contained fluid is at ambient pressures. We recognize three separate stages in the behaviour of the fracture as the water pressure decreases due to drainage, namely opening, closing and, finally, completely closed. In the first stage, the water will initially drain out along conduits between the faces of the opening fracture. As the fracture opens and grows, it can become long enough to intersect any pre-existing network of open joints and fractures. If this connection occurs, the opened fracture will stop growing and there will be further rapid flow of fluid out along these open conduits. Even if connection to an open conduit does not occur, as the pressure drops the fracture will stop growing and eventually halt. At the pressure where the fracture tip begins to close, drainage into any open conduits will cease and drainage during fracture closing will thus be by flow into the porous walls of the fracture and borehole. This effect of the conduit being shut off may manifest as a kink in the observed height decay of the slug at this pressure due to the effect of the different drainage terms. As the water pressure decreases, the fracture will continue to close until the water pressure drops below the critical pressure and the fracture becomes completely closed. This total closure occurs at the same pressure as that necessary to re-open the fracture, and has been taken equal to the minimum principal stress. This total closure of the fracture also may manifest as a kink in the observed height decay of the slug, as below this closure pressure is the third, and final, stage where drainage is at the usual low pressures controlled by the permeability of the rock just around the packed-off interval of the borehole. We develop the most general mathematical model for this phenomenon using conservation of mass of the water in the borehole-fracture system for each separate pressure stage. This model of fracture behaviour, with separate flow stages and pressure regimes, explains the physical basis for the increasingly used technique in high-pressure injection tests of measuring the so-called re-opening pressure at the sharp change in slope of the rate of change of water height vs. height. Analytic and numerical solutions yield the rate at which the water height drops during the slug test and comparison with high-pressure test data confirms the accuracy of the model. Analysis of data from Yucca Mountain, Nevada, gives values for the fracture opening pressures that indicate the presence of extensional stresses of significant magnitude and of zones of pre-existing open fracture networks in the vicinity of the proposed repository site.
Copyright 2015 John Bruce Davies, Ph.D.. All rights reserved.
ph: 805 886 9856
john