Algorithmic approach to discrete fracture network flow modeling in consideration of realistic connections in large-scale fracture networks

Analyzing rock mass seepage using the discrete fracture network (DFN) flow model poses challenges when dealing with complex fracture networks. This paper presents a novel DFN flow model that incorporates the actual connections of large-scale fractures. Notably, this model efficiently manages over 20,000 fractures without necessitating adjustments to the DFN geometry. All geometric analyses, such as identifying connected fractures, dividing the two-dimensional domain into closed loops, triangulating arbitrary loops, and refining triangular elements, are fully automated. The analysis processes are comprehensively introduced, and core algorithms, along with their pseudo-codes, are outlined and explained to assist readers in their programming endeavors. The accuracy of geometric analyses is validated through topological graphs representing the connection relationships between fractures. In practical application, the proposed model is employed to assess the water-sealing effectiveness of an underground storage cavern project. The analysis results indicate that the existing design scheme can effectively prevent the stored oil from leaking in the presence of both dense and sparse fractures. Furthermore, following extensive modification and optimization, the scale and precision of model computation suggest that the proposed model and developed codes can meet the requirements of engineering applications.