Horizontal well multi-cluster fracturing technology is crucial for the economic development of fractured shale reservoirs. The abundance of natural fractures in shale reservoirs significantly influences the propagation path of hydraulic fractures and determines the formation of complex fracture networks. To investigate the impact of natural weak planes on the geometric parameters of fractures in shale reservoirs, we first conducted tests on the mechanical characteristics of core samples from outcropping shale in the Weiyuan area using the indoor three-point bending method and digital image correlation (DIC) technology, providing data validation for subsequent numerical models. Secondly, considering the interaction between hydraulic and natural weak planes in three-dimensional space, we established a three-dimensional numerical model for horizontal well fracturing to simulate the synchronous competition and expansion of fractures in multi-cluster fracturing. Based on this foundation, we analyzed the influence of formation parameters and engineering parameters on the formation patterns of complex fracture networks. The results indicate that the difference in in situ stress is a significant factor affecting the selection of fracture propagation paths. As the in situ stress difference increases, it becomes more challenging to open natural fractures, leading to a reduced probability of activation of natural weak interfaces. When the cohesive strength of natural fractures is smaller, they are more likely to open and capture hydraulic fractures, thereby increasing shear slip length and fracture network area. Each fracturing stage has an optimal perforation density combination, where a higher density of perforations leads to reduced perforation pressure drop and weaker ability to mitigate inter-cluster stress interference. To achieve a comprehensive and balanced development of multi-clusters, the inter-cluster stress interference can be alleviated by increasing the perforation pressure drop. For dense perforation clusters, higher injection rates and viscosity can be employed to ensure the uniform development of multiple perforation clusters. This study provides new insights into predicting the formation of complex fracture networks in shale reservoirs and offers valuable guidance for optimizing hydraulic fracturing designs.