Nonlinear seepage behavior within rock fractures represents a critical and actively researched challenge in underground engineering, energy exploitation, and environmental sciences. Through the integration of nonlinear seepage theory with coupled numerical simulations of fracture flow and matrix flow, this study systematically investigates the synergistic mechanisms governing the influence of filling particles, tortuous fractures, and porous matrices on fluid transport within fracture–porous matrix seepage systems. Key findings reveal that: (1) Horizontal fractures continuously receive fluid influx from the surrounding porous matrix, where the flow field maintains remarkable symmetry, with a critical matrix height-to-fracture aperture ratio regulating streamline divergence and convergence at the fracture outlet; (2) The flow field within horizontal fractures undergoes substantial transformation when the Reynolds number exceeds a critical threshold, while maintaining stable flow patterns and -ΔP-Q relationships below this value, demonstrating a distinct inertial-controlled flow regime transition; (3) Tortuous fracture geometries induce localized vortex formation and significant velocity fluctuations, particularly in the front and rear dip-angle zones, substantially enhancing fluid exchange efficiency compared to horizontal configurations; (4) The volumetric flow rate exhibits a non-monotonic relationship with inclination angle, peaking at approximately 36°, while a synergistic effect between fracture inclination and infill particle diameter systematically modulates pressure-drop-flow-rate relationships, with a critical d/h = 0.5 threshold distinguishing fundamentally different flow behaviors. These findings provide quantitative criteria for predicting nonlinear seepage in practical engineering scenarios involving complex fracture networks and filling materials, offering significant implications for risk assessment and drainage design in deep underground projects.
A Thu, study studied this question.