This study employs the discrete element method for a series of coupled hydro-mechanical analyses aimed at investigating the hydraulic stimulation effect on a slate geothermal reservoir. Initially, a ubiquitous foliation model (UFM) is devised to characterize the mechanical properties of slate. The UFM consists of a foliation model for near-field discontinuity features and a ubiquitous model for far-field anisotropic rock behavior. The devised model undergoes validation through an established benchmark analysis considering isotropic conditions. The analysis is subsequently extended to investigate anisotropic shear response influenced by the foliation. In cases where the anisotropic angle equals 90°, new fractures initiate along the foliation, which leads to a concentrated and severe damage pattern. In contrast, models with the anisotropic angle of 45° and 135° exhibit minor failures due to foliation orientation hindering fracture propagation. Lastly, an anisotropic hydraulic-mechanical analysis is undertaken to assess the hydraulic stimulation influence and its consequential effects on a slate geothermal reservoir. The results indicate that higher injection rates amplify damage to the rock matrix and the propagation of fractures, which are concentrated between open fractures and adjacent foliations.