<p>Simulation of geological fractured media remains a significant research focus in computational geomechanics, with complexities arising from crack interactions and multi-physical coupling. The consideration of three-dimensional (3D) multiple cracks further complicates this problem. To address this issue, we propose a novel 3D peridynamic framework for simulating thermal-induced cracking in fractured media. This framework is built upon a thermo-mechanical peridynamic (TM-PD) model, grounded in thermoelasticity and peridynamic theory, and incorporates the influences of multiple fractures and non-uniform temperature distributions. To facilitate efficient simulations of quasi-static cracking, an adaptive dynamic relaxation method with a modified damage criterion is developed. Then, numerical accuracy and convergence are assessed through benchmark tests, demonstrating the efficiency of the TM-PD model. Characteristics of cracking for various fracture patterns under non-uniform temperature distributions are studied, with a quantitative analysis of variations in kinetic energy, damage accumulation, and deformation at crack tips. Finally, a complex fractured medium containing multiple random/parallel fractures is constructed for further exploration. Simulation results reveal that cracking and damage evolutions under non-uniform temperature conditions significantly differ from those under uniform temperature conditions. These disparities are attributed to the complex interactions, coalescence, and penetration among fractures.</p>