The validity of the cubic law for laminar flow of fluids through open fractures consisting of parallel planar plates has been established by others over a wide range of conditions with apertures ranging down to a minimum of 0.2 µm. The law may be given in simplified form by Q/Δh = C(2b)3, where Q is the flow rate, Δh is the difference in hydraulic head, C is a constant that depends on the flow geometry and fluid properties, and 2b is the fracture aperture. The validity of this law for flow in a closed fracture where the surfaces are in contact and the aperture is being decreased under stress has been investigated at room temperature by using homogeneous samples of granite, basalt, and marble. Tension fractures were artificially induced, and the laboratory setup used radial as well as straight flow geometries. Apertures ranged from 250 down to 4µm, which was the minimum size that could be attained under a normal stress of 20 MPa. The cubic law was found to be valid whether the fracture surfaces were held open or were being closed under stress, and the results are not dependent on rock type. Permeability was uniquely defined by fracture aperture and was independent of the stress history used in these investigations. The effects of deviations from the ideal parallel plate concept only cause an apparent reduction in flow and may be incorporated into the cubic law by replacing C by C/ƒ. The factor ƒ varied from 1.04 to 1.65 in these investigations. The model of a fracture that is being closed under normal stress is visualized as being controlled by the strength of the asperities that are in contact. These contact areas are able to withstand significant stresses while maintaining space for fluids to continue to flow as the fracture aperture decreases. The controlling factor is the magnitude of the aperture, and since flow depends on (2b)3, a slight change in aperture evidently can easily dominate any other change in the geometry of the flow field. Thus one does not see any noticeable shift in the correlations of our experimental results in passing from a condition where the fracture surfaces were held open to one where the surfaces were being closed under stress.

/pdf/validity-of-cubic-law-for-fluid-flow-in-a-deformable-rock-2ynq4id3hq.pdf

VALIDITY OF CUBIC LAW FOR FLUID FLOW IN A DEFORMABLE ROCK FRACTURE - eScholarship

A review of recent developments concerning the structure, mechanics and fluid flow properties of fault zones

Heat carried by ground water serves as a tracer to identify surface water infiltration, flow through fractures, and flow patterns in ground water basins. Temperature measurements can be analyzed for recharge and discharge rates, the effects of surface warming, interchange with surface water, hydraulic conductivity of streambed sediments, and basin-scale permeability. Temperature data are also used in formal solutions of the inverse problem to estimate ground water flow and hydraulic conductivity. The fundamentals of using heat as a ground water tracer were published in the 1960s, but recent work has significantly expanded the application to a variety of hydrogeological settings. In recent work, temperature is used to delineate flows in the hyporheic zone, estimate submarine ground water discharge and depth to the salt-water interface, and in parameter estimation with coupled ground water and heat-flow models. While short reviews of selected work on heat as a ground water tracer can be found in a number of research papers, there is no critical synthesis of the larger body of work found in the hydrogeological literature. The purpose of this review paper is to fill that void and to show that ground water temperature data and associated analytical tools are currently underused and have not yet realized their full potential.

http://www.geology.wisc.edu/~andy/GW_Heat%20as%20a%20Tracer.pdf

Heat as a ground water tracer.

Ground-penetrating radar and its use in sedimentology: principles, problems and progress

New observations of the internal structure of the San Gabriel fault (SGF) are combined with previous characterizations of the Punchbowl fault (PF) to evaluate possible explanations for the low frictional strength and seismic characteristics of the San Andreas fault (SAF). The SGF and PF are ancient, large-displacement faults of the SAF system exhumed to depths of 2 to 5 km. These fault zones are internally zoned; the majority of slip was confined to the cores of principal faults, which typically consist of a narrow layer (less than tens of centimeters) of ultracataclasite within a zone of foliated cataclasite several meters thick. Each fault core is bounded by a zone of damaged host rock of the order of 100 m thick. Orientations of subsidiary faults and other fabric elements imply that (1) the maximum principal stress was oriented at large angles to principal fault planes, (2) strain was partitioned between simple shear in the fault cores and nearly fault-normal contraction in the damaged zones and surrounding host rock, and (3) the principal faults were weak. Microstructures and particle size distributions in the damaged zone of the SGF imply deformation was almost entirely cataclastic and can be modeled as constrained comminution. In contrast, cataclastic and fluid-assisted processes were significant in the cores of the faults as shown by pervasive syntectonic alteration of the host rock minerals to zeolites and clays and by folded, sheared, and attenuated cross-cutting veins of laumontite, albite, quartz, and calcite. Total volume of veins and neocrystallized material reaches 50% in the fault core, and vein structure implies episodic fracture and sealing with time-varying and anisotropic permeability in the fault zone. The structure of the ultracataclasite layer reflects extreme slip localization and probably repeated reworking by particulate flow at low effective stresses. The extreme slip localization reflects a mature internal fault structure resulting from a positive feedback between comminution and transformation weakening. The structural, mechanical, and hydrologic characteristics of the Punchbowl and San Gabriel faults support the model for a weak San Andreas based on inhomogeneous stress and elevated pore fluid pressures contained within the core of a seismogenic fault. Elevated fluid pressures could be repeatedly generated in the core of the fault by a combination of processes including coseismic dilatancy and creation of fracture permeability, fault-valve behavior to recharge the fault with fluid, post-seismic self-sealing of fracture networks to reduce permeability and trap fluids, and time-dependent compaction of the core to generate high pore pressure. The localized slip and fluid-saturated conditions are wholly compatible with additional dynamic weakening by thermal pressurization of fluids during large seismic slip events, which can help explain both the low average strength of the San Andreas and seismogenic characteristics such as large stress relief. In addition, such a dynamic weakening mechanism is expected only in mature fault zones and thus could help explain the apparent difference in strength of large-displacement faults from smaller-displacement, subsidiary seismogenic faults.

Internal structure and weakening mechanisms of the San Andreas Fault

Fault zone architecture and related permeability structures form primary controls on fluid flow in upper-crustal, brittle fault zones. We develop qualitative and quantitative schemes for evaluating fault-related permeability structures by using results of field investigations, laboratory permeability measurements, and numerical models offlow within andnearfaultzones.Thequalitativeschemecomparesthepercentageofthetotalfaultzone width composed of fault core materials (e.g., anastomosing slip surfaces, clay-rich gouge, cataclasite,andfaultbreccias)tothepercentageofsubsidiarydamagezonestructures(e.g., kinematically related fracture sets, small faults, and veins). A more quantitative scheme is developed to define a set of indices that characterize fault zone architecture and spatial variability.Thefaultcoreanddamagezonearedistinctstructuralandhydrogeologicunits that reflect the material properties and deformation conditions within a fault zone. Whether a fault zone will act as a conduit, barrier, or combined conduit-barrier system is controlled by the relative percentage of fault core and damage zone structures and the inherent variability in grain scale and fracture permeability. This paper outlines a frameworkforunderstanding,comparing,andcorrelatingthefluidflowpropertiesoffaultzones in various geologic settings.

/pdf/fault-zone-architecture-and-permeability-structure-ygkia82fr8.pdf

Fault zone architecture and permeability structure

https://works.bepress.com/james_evans/139/download/

Permeability of fault-related rocks, and implications for hydraulic structure of fault zones

Fault zones in the upper crust are typically composed of complex fracture networks and discrete zones of comminuted and geochemically altered fault rocks. Determining the patterns and rates of fluid flow in these distinct structural discontinuities is a three-dimensional problem. A series of numerical simulations of fluid flow in a set of three-dimensional discrete fracture network models aids in identifying the primary controlling parameters of fault-related fluid flow, and their interactions, throughout episodic deformation. Four idealized, but geologically realistic, fault zone architectural models are based on fracture data collected along exposures of the Stillwater Fault Zone in Dixie Valley, Nevada and geometric data from a series of normal fault zones in east Greenland. The models are also constrained by an Andersonian model for mechanically compatible fracture networks associated with normal faulting. Fluid flow in individual fault zone components, such as a fault core and damage zone, and full outcrop scale model domains are simulated using a finite element routine. Permeability contrasts between components and permeability anisotropy within components are identified as the major controlling factors in fault-related fluid flow. Additionally, the structural and hydraulic variations in these components are also major controls of flow at the scale of the full model domains. The four models can also be viewed as a set of snapshots in the mechanical evolution of a single fault zone. Changes in the hydraulic parameters within the models mimic the evolution of the permeability structure of each model through a single deformation cycle. The model results demonstrate that small changes in the architecture and hydraulic parameters of individual fault zone components can have very large impacts, up to five orders of magnitude, on the permeability structure of the full model domains. Closure of fracture apertures in each fault zone magnifies the magnitude and orientation of permeability anisotropy in ways that are closely linked to the implicitly modeled deformation. Changes in fault zone architecture can cause major changes in permeability structure that, in turn, significantly impact the magnitude and patterns of fluid flux and solute transport both within and near the fault zone. Inferences derived from the model results are discussed in the context of the mechanical strength of an evolving fault zone, fault zone sealing mechanisms which control the conduit-barrier systematics of a fault zone as a flow system, and how these processes are related to fluid flow in natural fault zones.

Fault Zone Architecture and Fluid Flow: Insights from Field Data and Numerical Modeling

We determined the structure and permeability variations of a 4 km-long normal fault by integrating surface mapping with data from five boreholes drilled through the fault (borehole to tens of meters scale). The Big Hole fault outcrops in the Jurassic Navajo Sandstone, central Utah. A total of 363.2 m of oriented drill core was recovered at two sites where fault displacement is 8 and 3-5 m. The main fault core is a narrow zone of intensely comminuted grains that is a maximum of 30 cm thick and is composed of low-porosity amalgamated deformation bands that have slip surfaces on one or both sides. Probe permeameter measurements showed a permeability decline from greater than 2000 to less than 0.1 md as the fault is approached. Whole-core analyses showed that fault core permeability is less than 1 md and individual deformation band permeability is about 1 md. Using these data, we calculated the bulk permeability of the fault zone. Calculated transverse permeability over length scales of 5-10 m is 30-40 md, approximately 1-4% the value of the host rock. An inverse power mean calculation (representing a fault array with complex geometry) yielded total fault-zone permeabilities of 7-57 md. The bulk fault-zone permeability is most sensitive to variations in fault core thickness, which exhibits the greatest variability of the fault components.

/pdf/structural-heterogeneity-and-permeability-in-faulted-eolian-51l6rx5tdx.pdf

Structural heterogeneity and permeability in faulted eolian sandstone: Implications for subsurface modeling of faults

Groundwater flow in high-relief mountainous terrain is modeled to examine how geology, surface topography, climate, and regional heat flow influence the pattern and magnitude of an advective thermal disturbance. A numerical procedure is used to estimate the position of the water table within constraints provided by the available infiltration rate and the permeability of the mountain massif. Results show that where the water table is located at depth within the mountain massif, the rate of groundwater recharge, rather than permeability, is the appropriate factor to characterize the potential for an advective disturbance. Thermal disturbances should be assessed within the framework of the two-dimensional character of the flow system because interpretations based on one-dimensional models are prone to significant error. Modeling of site-specific systems can exploit the existence of a permeability “window,” where temperatures in discharge areas and in springs discharging from fracture zones reach peak values. Temperatures of thermal springs reflect the complex interaction between flow within the mountain massif and flow through permeable fracture zones. Consequently, simple calculations relating the geothermal gradient, depth of maximum groundwater circulation, and discharge temperature improperly represent the physics of the process.

/pdf/the-influence-of-groundwater-flow-on-thermal-regimes-in-19h129ugct.pdf

Craig B. Forster

Papers

Fault zone architecture and permeability structure

Permeability of fault-related rocks, and implications for hydraulic structure of fault zones

Fault Zone Architecture and Fluid Flow: Insights from Field Data and Numerical Modeling

Structural heterogeneity and permeability in faulted eolian sandstone: Implications for subsurface modeling of faults

The influence of groundwater flow on thermal regimes in mountainous terrain: A model study