Thermally induced damage has an important influence on rock mechanics and engineering, especially for high-level radioactive waste repositories, geological carbon storage, underground coal gasification, and hydrothermal systems. Additionally, the wide application of geothermal heat requires knowledge of the geothermal conditions of reservoir rocks at elevated temperature. However, few methods to date have been reported for investigating the micro-mechanics of specimens at elevated temperatures. Therefore, this paper uses a cluster model in particle flow code in two dimensions (PFC^2D) to simulate the uniaxial compressive testing of Australian Strathbogie granite at various elevated temperatures. The peak strength and ultimate failure mode of the granite specimens at different elevated temperatures obtained by the numerical methods are consistent with those obtained by experimentation. Since the tensile force is always concentrated around the boundary of the crystal, cracks easily occur at the intergranular contacts, especially between the b-b and b-k boundaries where less intragranular contact is observed. The intergranular and intragranular cracking of the specimens is almost constant with increasing temperature at low temperature, and then it rapidly and linearly increases. However, the inflection point of intergranular micro-cracking is less than that of intragranular cracking. Intergranular cracking is more easily induced by a high temperature than intragranular cracking. At an elevated temperature, the cumulative micro-crack counts curve propagates in a stable way during the active period, and it has no unstable crack propagation stage. The micro-cracks and parallel bond forces in the specimens with elevated temperature evolution and axial strain have different characteristics than those at lower temperature. More branch fractures and isolated wider micro-cracks are generated with increasing temperature when the temperature is over 400 °C. Therefore, the total number of cracks is almost constant when the temperature is below 400 °C; next, it linearly increases when the temperature is over 400 °C. This trend is the same as that observed by experimentation.

Drilling vertical and horizontal wellbores in the shale reservoir may trigger the in-situ stress release around the wellbore walls and change the original stress equilibrium state, leading the wellbores to instability. This stress change in the wellbore corresponds to the stress paths of confining pressure unloading and axial stress loading under laboratory conditions. In this paper, according to the conventional triaxial compression test results, laboratory experiments and DEM simulations by PFC^2D were conducted to deeply study the strength, failure, strain energy evolution, and micro-crack damage mechanism of shale specimens under confining pressure unloading conditions. The shale specimens at different bedding inclinations were tested under different initial axial stress levels and confining pressure unloading rates, with fixed initial unloading confining pressure. This research revealed that confining pressure unloading induces greater plastic deformation, more micro-crack damage and strain energy dissipation, and a more complex failure pattern. The strain energy dissipation and dilatation under confining pressure unloading conditions are mainly induced by the generation and accumulation of tensile cracks. Moreover, the unloading rate has a significant effect on the mechanical properties, and the high unloading rate enhances the failure strength and induces more strain energy dissipation and micro tensile cracks. For the wellbore drilling in shale formations, when the buried depth and vertical stress are fixed, the lower the lateral stress is, the easier it is to form tensile failure around the wellbore wall in the drilling process, and the more induced fractures will be generated in the formation around the wellbore.