Fracture Expansion and Closure in Overburden: Mechanisms Controlling Dynamic Water Inflow to Underground Reservoirs in Shendong Coalfield

The construction of underground reservoirs in coal goafs is an innovative technology to alleviate the coal–water conflict in arid mining areas of northwest China. However, its widespread application is constrained by the challenge of accurately predicting water inflow, which fluctuates significantly due to the dynamic “expansion–closure” behavior of mining-induced fractures. This study focuses on the Shendong mining area, where repeated multi-seam mining occurs, and employs a coupled Finite Discrete Element Method (FDEM) and Computational Fluid Dynamics (CFD) numerical model, combined with in situ tests such as drilling fluid loss and groundwater level monitoring, to quantify the evolution of overburden fractures and their impact on reservoir water inflow during mining, 8 months post-mining, and after 7 years. The results demonstrate that the height of the water-conducting fracture zone decreased from 152 m during mining to 130 m after 7 years, while fracture openings in the key aquifer and aquitard were reduced by over 50%. This closure process caused a dramatic decline in water inflow from 78.3 m3/h to 2.6 m3/h—a reduction of 96.7%. The CFD-FDEM simulations showed a deviation of only 10.6% from field measurements, confirming fracture closure as the dominant mechanism driving inflow attenuation. This study reveals how fracture closure shifts water flow patterns from vertical to lateral recharge, providing a theoretical basis for optimizing the design and sustainable operation of underground reservoirs.