Abstract: Fracture propagation is affected by multi-metal-veins formed by geological diagenesis in shale during the hydraulic fracturing. However, the influence of multi-metal-veins on fractures propagation remains unclear. To solve the problem, based on the semi-circle bending (SCB) test and the extended finite element (XFEM) theory, the interaction between multi-metal-veins and fractures is investigated. The experimental results reveal that the fractures usually deflect at the upper or lower interfaces between metal veins and rocks (e.g. the specimen S-2), which is different from the propagation behavior of fractures in calcite veins. Meanwhile, the fracture toughness of the specimen S-1 is 24.40% higher than that of the specimen S-2, indicating that the increasing of total thickness of multiple metal veins increases the resistance to the fracture vertical propagation. The simulation results show that the increasing of the number, total thickness of veins, the modulus difference between veins and rock, the approach angle and the notch angle all increase the resistance of the fracture passing through metal veins. The maximum deviation distance (Dmax) of the fracture decreases with the number of veins, while thickness combination types of metal veins do not affect Dmax. The reduction of the notch angle leads to the more tortuous fracture propagation path. Finally, we propose a new comprehensive fracture network pattern. Fracture networks are divided into two categories, including orthogonal fracture networks and sub-orthogonal fracture networks, and then divided into six sub-categories further. The research results will provide reference for hydraulic fracturing of shale reservoirs containing multi-metal-veins.