Multi-Factor Analysis and Graded Remediation Strategy for Goaf Stability in Underground Metal Mines: Fluid–Solid Coupling Simulation and Genetic Algorithm-Based Optimization Approach

To ensure the green, safe, and efficient extraction of mineral resources and promote sustainability, the stability of mined-out areas has become a critical factor affecting safe production and ecological restoration in underground metal mines. The instability of underground goafs poses a significant threat to mine safety, especially when irregular excavation patterns interact with high ground stress, exacerbating instability risks. Most existing studies lack a systematic and multidisciplinary integrated framework for comprehensive evaluation and management. This paper proposes a trinity research system of “assessment–optimization–governance”, integrating theoretical analysis, three-dimensional fluid–solid coupling numerical simulation, and a filling sequence optimization method based on genetic algorithms. An analysis of data measured from 243 pillars and 49 goafs indicates that approximately 20–30% of the pillars have a factor of safety (FoS) below 1.0, signaling immediate instability risks; additionally, 58% do not meet the threshold for long-term stability (FoS ≥ 1.5). Statistical and spatial analyses highlight that pillar width-to-height ratio (W/H) and cross-sectional area significantly influence stability; when W/H exceeds 1.5, FoS typically surpasses 2.0. Numerical simulations reveal pore water pressures of 1.4–1.8 MPa in deeper goafs, substantially reducing effective stress and accelerating plastic zone expansion. Stability classification categorizes the 49 goafs into 7 “poor”, 37 “moderate”, and 5 “good” zones. A genetic algorithm-optimized filling sequence prioritizes high-risk area remediation, reducing maximum principal stress by 60.96% and pore pressure by 28.6%. Cemented waste rock filling applied in high-risk areas, complemented by general waste rock filling in moderate-risk areas, significantly enhances overall stability. This integrated method provides a scientific foundation for stability assessment and dynamic remediation planning under complex hydrogeological conditions, offering a risk-informed and scenario-specific application of existing tools that improves engineering applicability.