Modeling fluid flow in a heterogeneous, fault-controlled hydrothermal system

Previous studies have shown that most hydrothermal systems discharging at the land surface are associated with faulting, and that the location, temperature and rate of discharge of these systems are controlled by the geometry and style of the controlling fault(s). Unfortunately, the transport of heat and fluid in fault-controlled hydrothermal systems is difficult to model realistically; although heterogeneity and anisotropy are assumed to place important controls on flow in faults, few data or observations are available to constrain the distribution of hydraulic properties within active faults. Here, analytical and numerical models are combined with geostatistical models of spatially varying hydraulic properties to model the flow of heat and fluid in the Borax Lake fault of south-east Oregon, USA. A geometric mean permeability within the fault of 7 x 10(-14) m(2) with 2x vertical/horizontal anisotropy in correlation length scale is shown to give the closest match to field observations. Furthermore, the simulations demonstrate that continuity of flow paths is an important factor in reproducing the observed behavior. In addition to providing some insight into possible spatial distributions of hydraulic properties at the Borax Lake site, the study highlights one potential avenue for integrating field observations with simulation results in order to gain greater understanding of fluid flow in faults and fault-controlled hydrothermal and petroleum reservoirs.