The hydraulic fracturing process is a prominent example of fracture network evolution under stress. However, the interactions between hydraulic fractures and natural fracture networks, along with the connectivity evolution of the resultant fracture networks, require more research. This research incorporates discrete fracture networks to characterize subsurface structures and employs the Discrete Element ‐ Lattice Boltzmann Method to simulate the hydraulic fracturing process. The dynamic evolution of subsurface structures, including the initiation of hydraulic fractures and their interaction with natural fractures, is systematically investigated. Results indicate that natural fractures significantly impact fracture initiation, propagation, and connectivity evolution, which in turn affects fluid production. Fracture strength is key for the interaction, and hydraulic fractures tend to propagate along weak natural fractures with low resistance. Fracture strength variability determines the final fracture networks, with low‐strength fractures breaking due to the altered in‐situ stress and forming local clusters. High injection rates and fluid viscosity result in a large pressure buildup and exaggerate the influential region. A multi‐cluster system is thus formed during the hydraulic fracturing process, and its connectivity can be well quantified with a novel connectivity metric. In low‐permeable reservoirs, fracture clusters connected to production wells can contribute instantly, while local clusters may contribute to production from a long‐term perspective. Injection rate, fluid viscosity, fracture orientation, and clustering effects have consistent positive correlations with the total connectivity and production. Heterogeneity has a weak positive correlation with fluid production, while a moderate negative correlation with total connectivity.