Experimental Study on Damage–Seepage Coupling of Small Faults Under Mining-Induced Stress Paths Based on Fractal Grading Method

To reveal the damage–seepage coupling mechanism of delayed floor water inrush induced by small fault activation under mining-induced stress, a cubic cement mortar specimen containing a persistent small fault was prepared based on similarity theory. Systematic triaxial loading–seepage tests were conducted under different fault fracture zone particle gradations, fracture zone widths, and fault angles, with simultaneous monitoring of stress–strain behavior, acoustic emission (AE) characteristics, and seepage flow evolution. The results show that: ① The peak strength decreases with increasing fracture zone width, but increases with increasing Talbot gradation coefficient (a fractal grading method) and fault angle. The failure mode transitions from shear-dominated to tension–shear composite failure. The spatial localization of AE events corresponds well with macroscopic fracture surfaces, and the AE source amplitude is positively correlated with compressive strength. ② The seepage flow exhibits a nonlinear evolution pattern of “compaction stabilization—stepwise rise—plateau stabilization” during loading. In the early loading stage, compaction of the fracture zone causes a slight decrease in flow. Approaching peak strength, the initiation and propagation of through-going fractures create interconnected seepage channels, leading to a stepwise jump in flow. In the post-peak stage, accompanied by fine particle erosion and framework reconfiguration, the flow tends to stabilize. A larger fracture zone width, smaller gradation coefficient, and smaller fault angle result in a more significant post-peak seepage surge, with the maximum flow rate reaching 3.6 times that of the specimen with a 2 mm wide fracture zone. ③ Grey relational analysis indicates that the fault angle is the most sensitive factor affecting the risk of delayed water inrush (correlation degree 0.788), followed by particle gradation and fracture zone width. The study demonstrates that under monotonic loading conditions, the damage evolution and seepage response of small faults are jointly controlled by their geometric parameters and internal structure, with the fractal grading method effectively quantifying the role of particle gradation. The findings provide a theoretical basis for risk assessment of delayed water inrush from small faults in working faces above confined aquifers.