Calcite veins, as weak structural planes in shale reservoirs, jointly influence hydraulic fracturing behavior along with the significant in-situ stress differences in northern Guizhou. Based on CT scanning and numerical simulation, this study constructs a three-dimensional heterogeneous model to elucidate the coupling mechanism between the horizontal stress differential coefficient ( \(K_{h}\) ) and the inclination angle that calcite veins have on the propagation of fluid fractures. The results indicate that, due to the angular effect of the principal stress and the failure axis, the failure stress of shales shows a tendency to initially increase and then decrease with the increase in calcite vein inclination. Meanwhile, the weak-plane effect of calcite veins determines the initial fracture propagation direction, allowing hydraulic fractures to preferentially extend along the vein strike before branching and diffusing into the shale matrix. A significant negative correlation is observed between \(K_{h}\) and the fracture complexity index ( \(F_{CI}\) ) during hydraulic fracturing: as \(K_{h}\) increases, \(F_{CI}\) decreases, leading to a more simplified fracture morphology. Furthermore, the dynamic evolution of the acoustic emission b-value demonstrates the link between magnitude distribution and hydraulic fracture complexity, with bigger b-value peaks corresponding to lower fracture complexity. This study looks at the coupling effect of geological weak planes and the in-situ stress field on hydraulic fracturing, giving theoretical support for fracture design and complex fracture network optimization in northern Guizhou shale reservoirs.