The interaction between multiple fractures is important in the analysis of rock fracture propagation, fracture network evolution and stability and integrity of rocks under hydro-mechanical (HM) coupling conditions. At present, modeling the mechanical behavior of multiple fractures is still challenging. Under the condition of multiple fractures, the opening, closing, sliding, propagation and penetration of fractures become more complicated. In order to simulate the HM coupling behavior of multi-fracture system, the paper presents a novel numerical scheme, including mesh reconstruction and topology generation algorithm, to efficiently and accurately represent fractures and their propagation process, and a potential function-based algorithm to address contact problem. The fracture contact algorithm does not need to set contact pairs and thus is suitable for complex contact situations from small to large deformations induced by HM loading. The topology of fracture interfaces is constructed by the dynamic adding algorithm, which makes the mesh reconstruction more rapid in the modeling of fracturing process, especially in the case of multiple fractures intersections. The numerical scheme is implemented in CASRock, a self-developed numerical code, to simulate the propagation process of rock fractures and the interaction of multiple fractures under the condition of HM coupling. To verify the suitability of the code, a series of tests were performed. The code was then applied to simulate hydraulic fracture propagation and fracture interactions caused by fluid injection. The ability of this method to study fracture propagation, multi-fracture interaction and fracture network evolution under hydro-mechanical coupling conditions is demonstrated.

Download PDF:

Keywords: CASRock, Multi-fracture interaction, Hydro-mechanical (HM) coupling, Efficient grid reconstruction and topology generation, Distributed contact computation