Abstract. The fracture distribution in basalt flows is a direct result of thermal processes. Thus basalt flows present a unique opportunity to characterize a nearly perfect deterministic system with its fundamental physical parameters. Fracture distribution data collected on cliff exposures of basalt flows near the Idaho National Engineering and Environmental Laboratory (INEEL) are combined with calculations of cooling rate and temperature distribution from a finite element model to construct a predictive methodology for fracture spacing. The methodology is based on an empirical power law relationship between inverse cooling rate and fracture spacing. The methodology may be applied to unexposed basalt flows of approximately elliptical cross section whose thickness and width are constrained only by geophysical or borehole data if sufficient fracture data on nearby exposed flows are available. The methodology aids waste remediation efforts at sites involving contaminant transport through fractured basalt, such as the INEEL and the Hanford site in Washington, as well as involving transport and fluid flow through volcanic or intrusive rocks where thermal processes are responsible for fracturing.
IntroductionColumnar tensile fractures in basalt flows are driven by thermal elastic strain due to differential cooling of the basalt, as opposed to regional tectonic extension. Given the low matrix permeability of basalt, characterizing the spatial distribution of thermal contraction fractures is the key to understanding the overall permeability structure of multiple basalt flows. The goal of the geological portion of the integrated study, and the subject of this paper, is the development of a methodology to predict fracture distribution in basalt flows when only limited knowledge of the flow geometry is available. Obtaining fracture information on basalt flows not exposed at the surface is challenging, as geophysical methods have resolution falling below that required to image the fractures directly and core fracture measurements sample a very small volume. The We integrate field data on basalt flow geometry and fracture spacing with numerical constraints on the thermal history during fracturing. The result provides a methodology to estimate fracture spacing as a function of cooling rate, controlled by flow thickness and geometry. Although this work is part of a study for hazardous waste containment, aspects are directly applicable to water, oil, or gas reservoirs in thermally fractured volcanic or intrusive rocks.