Geologic data are examined to evaluate whether certain geologic characteristics occur in higher abundance or higher magnitude along production geothermal wells relative to non-productive wells. We perform 3D geologic mapping, 3D stress modeling, and fault-slip modeling to estimate fourteen different geologic factors that are hypothesized to control or correlate with well productivity. The geologic factors are; heat, fault-damage zone thickness, distance from active faults, fault intersection/termination density, fault curvature, slip tendency of faults, dilation tendency of faults, dilation resulting from modeled fault slip, normal stress reduction resulting from modeled fault slip, Coulomb shear stress increase resulting from modeled fault slip, the summed thickness of ‘favorable’ lithologies within a borehole, the summed thickness of fault damage zones in favorable lithologies within a borehole, the distance along the borehole to the nearest geologic contact, and the thickness of individual stratigraphic units. These geologic factors are quantified along fifty wells at Brady geothermal system, including twelve production wells and thirty-one non-productive wells. Results indicate that geologic factors such as stress changes associated with faulting, nearness to and thickness of fault zones, distance from geologic contacts, and heat occur in higher magnitude or higher abundance along production wells relative to non-productive wells. These geologic factors may play an important role in controlling the locations and distribution of fluid circulation in geothermal fields.