©2017. American Geophysical Union. All Rights Reserved. Acoustic borehole televiewer (BHTV) logs provide measurements of fracture attributes (orientations, thickness, and spacing) at depth. Orientation, censoring, and truncation sampling biases similar to those described for one-dimensional outcrop scanlines, and other logging or drilling artifacts specific to BHTV logs, can affect the interpretation of fracture attributes from BHTV logs. K-means, fuzzy K-means, and agglomerative clustering methods provide transparent means of separating fracture groups on the basis of their orientation. Fracture spacing is calculated for each of these fracture sets. Maximum likelihood estimation using truncated distributions permits the fitting of several probability distributions to the fracture attribute data sets within truncation limits, which can then be extrapolated over the entire range where they naturally occur. Akaike Information Criterion (AIC) and Schwartz Bayesian Criterion (SBC) statistical information criteria rank the distributions by how well they fit the data. We demonstrate these attribute analysis methods with a data set derived from three BHTV logs acquired from the high-temperature Rotokawa geothermal field, New Zealand. Varying BHTV log quality reduces the number of input data points, but careful selection of the quality levels where fractures are deemed fully sampled increases the reliability of the analysis. Spacing data analysis comprising up to 300 data points and spanning three orders of magnitude can be approximated similarly well (similar AIC rankings) with several distributions. Several clustering configurations and probability distributions can often characterize the data at similar levels of statistical criteria. Thus, several scenarios should be considered when using BHTV log data to constrain numerical fracture models.

©2017. American Geophysical Union. All Rights Reserved. Acoustic borehole televiewer (BHTV) logs provide measurements of fracture attributes (orientations, thickness, and spacing) at depth. Orientation, censoring, and truncation sampling biases similar to those described for one-dimensional outcrop scanlines, and other logging or drilling artifacts specific to BHTV logs, can affect the interpretation of fracture attributes from BHTV logs. K-means, fuzzy K-means, and agglomerative clustering methods provide transparent means of separating fracture groups on the basis of their orientation. Fracture spacing is calculated for each of these fracture sets. Maximum likelihood estimation using truncated distributions permits the fitting of several probability distributions to the fracture attribute data sets within truncation limits, which can then be extrapolated over the entire range where they naturally occur. Akaike Information Criterion (AIC) and Schwartz Bayesian Criterion (SBC) statistical information criteria rank the distributions by how well they fit the data. We demonstrate these attribute analysis methods with a data set derived from three BHTV logs acquired from the high-temperature Rotokawa geothermal field, New Zealand. Varying BHTV log quality reduces the number of input data points, but careful selection of the quality levels where fractures are deemed fully sampled increases the reliability of the analysis. Spacing data analysis comprising up to 300 data points and spanning three orders of magnitude can be approximated similarly well (similar AIC rankings) with several distributions. Several clustering configurations and probability distributions can often characterize the data at similar levels of statistical criteria. Thus, several scenarios should be considered when using BHTV log data to constrain numerical fracture models.

©2017. American Geophysical Union. All Rights Reserved. Analysis of fracture orientation, spacing, and thickness from acoustic borehole televiewer (BHTV) logs and cores in the andesite-hosted Rotokawa geothermal reservoir (New Zealand) highlights potential controls on the geometry of the fracture system. Cluster analysis of fracture orientations indicates four fracture sets. Probability distributions of fracture spacing and thickness measured on BHTV logs are estimated for each fracture set, using maximum likelihood estimations applied to truncated size distributions to account for sampling bias. Fracture spacing is dominantly lognormal, though two subordinate fracture sets have a power law spacing. This difference in spacing distributions may reflect the influence of the andesitic sequence stratification (lognormal) and tectonic faults (power law). Fracture thicknesses of 9–30 mm observed in BHTV logs, and 1–3 mm in cores, are interpreted to follow a power law. Fractures in thin sections (∼5 μm thick) do not fit this power law distribution, which, together with their orientation, reflect a change of controls on fracture thickness from uniform (such as thermal) controls at thin section scale to anisotropic (tectonic) at core and BHTV scales of observation. However, the ∼5% volumetric percentage of fractures within the rock at all three scales suggests a self-similar behavior in 3-D. Power law thickness distributions potentially associated with power law fluid flow rates, and increased connectivity where fracture sets intersect, may cause the large permeability variations that occur at hundred meter scales in the reservoir. The described fracture geometries can be incorporated into fracture and flow models to explore the roles of fracture connectivity, stress, and mineral precipitation/dissolution on permeability in such andesite-hosted geothermal systems.

©2018. American Geophysical Union. All Rights Reserved. Ultrasonic image logs acquired in the DFDP-2B borehole yield the first continuous, subsurface description of the transition from schist to mylonite in the hangingwall of the Alpine Fault, New Zealand, to a depth of 818 m below surface. Three feature sets are delineated. One set, comprising foliation and foliation-parallel veins and fractures, has a constant orientation. The average dip direction of 145° is subparallel to the dip direction of the Alpine Fault, and the average dip magnitude of 60° is similar to nearby outcrop observations of foliation in the Alpine mylonites that occur immediately above the Alpine Fault. We suggest that this foliation orientation is similar to the Alpine Fault plane at ∼1 km depth in the Whataroa valley. The other two auxiliary feature sets are interpreted as joints based on their morphology and orientation. Subvertical joints with NW-SE (137°) strike occurring dominantly above ∼500 m are interpreted as being formed during the exhumation and unloading of the Alpine Fault's hangingwall. Gently dipping joints, predominantly observed below ∼500 m, are interpreted as inherited hydrofractures exhumed from their depth of formation. These three fracture sets, combined with subsidiary brecciated fault zones, define the fluid pathways and anisotropic permeability directions. In addition, high topographic relief, which perturbs the stress tensor, likely enhances the slip potential and thus permeability of subvertical fractures below the ridges, and of gently dipping fractures below the valleys. Thus, DFDP-2B borehole observations support the inference of a large zone of enhanced permeability in the hangingwall of the Alpine Fault.

©2018. American Geophysical Union. All Rights Reserved. Ultrasonic image logs acquired in the DFDP-2B borehole yield the first continuous, subsurface description of the transition from schist to mylonite in the hangingwall of the Alpine Fault, New Zealand, to a depth of 818 m below surface. Three feature sets are delineated. One set, comprising foliation and foliation-parallel veins and fractures, has a constant orientation. The average dip direction of 145° is subparallel to the dip direction of the Alpine Fault, and the average dip magnitude of 60° is similar to nearby outcrop observations of foliation in the Alpine mylonites that occur immediately above the Alpine Fault. We suggest that this foliation orientation is similar to the Alpine Fault plane at ∼1 km depth in the Whataroa valley. The other two auxiliary feature sets are interpreted as joints based on their morphology and orientation. Subvertical joints with NW-SE (137°) strike occurring dominantly above ∼500 m are interpreted as being formed during the exhumation and unloading of the Alpine Fault's hangingwall. Gently dipping joints, predominantly observed below ∼500 m, are interpreted as inherited hydrofractures exhumed from their depth of formation. These three fracture sets, combined with subsidiary brecciated fault zones, define the fluid pathways and anisotropic permeability directions. In addition, high topographic relief, which perturbs the stress tensor, likely enhances the slip potential and thus permeability of subvertical fractures below the ridges, and of gently dipping fractures below the valleys. Thus, DFDP-2B borehole observations support the inference of a large zone of enhanced permeability in the hangingwall of the Alpine Fault.