How can we improve estimates of bulk fault zone hydraulic properties?

The fluid flow properties of faults are highly variable and spatially heterogeneous. We use numerical simulation of flow through field maps of detailed fault zone architecture to demonstrate that flow across the fault zone is controlled by connected high-permeability pathways, which are highly tortuous in mapped fault outcrops. Such small-scale, geometrically complex, fault zone architectural features can never be resolved for subsurface faults. Consequently, the key to prediction of subsurface bulk fault zone hydraulic properties is a statistical characterisation of the likelihood and frequency of such connected pathways. We demonstrate for a single architectural feature, the fault core, that thickness variation along strike can be well described by a spatially correlated random field with a spherical covariance structure. These data are from a single site in a specific lithology. To enable such statistics to be used to make predictions at other sites, a large number of similar datasets must be pooled. This will enable us to relate such spatial statistics to gross properties such as host rock lithology and fault throw, which are measurable for subsurface faults.