Fracture Assessment of Quasi-brittle Rock Simulated by Modified Discrete Element Method

Abstract

New developments of an in-house hybrid code, named Modified Discrete Element Method (MDEM) are presented in the paper. The new developments are on the treatment of pre-existing and propagating fractures in quasi-brittle materials. These developments are the embedment of Linear Elastic Fracture Mechanics (LEFM) and elastic-softening crack band model -based methodologies in the MDEM and their application in lab and reservoir scale. Using the first methodology, MDEM can calculate stress intensity factors, 𝐾�I and 𝐾�II using the internal contact forces of particles. 𝐾�I and 𝐾�II are calculated independent of boundary conditions and geometrical configuration with acceptable accuracy level. The methodology has been also used in reservoir scale to study the rupture likelihood of faults and fractures due to fluid injection. This methodology enables the code to model mode I and mode II failures and propagation direction based on the fracturing model proposed by Rao et al. (Int J Rock Mech Min Sci 40(3): 355–375, 2003). Using the second methodology, the MDEM can model nonlinear behavior of quasi-brittle materials including or excluding preexisting cracks based on fracture energy. A model was verified against an experiment of a three point bend test with a notch. The numerically obtained force-crack mouth opening curve was reasonably comparable to the experimental test. The analysis was repeated for three other mesh sizes and the results are less mesh size dependent. Finally, it was shown that MDEM has the potential in studying fracture mechanics of quasi-brittle materials both in lab and large-scale investigations.