Crack micro-geometries and tribological properties pose an important impact on the elastic characteristics of fractured rocks. Numerical simulation as a promising way for this issue still faces some challenges. With the rapid development of computers and computational techniques, discrete-based numerical approaches with desirable properties have been increasingly developed, but few attempts to consider the particle surface roughness in a lattice type model. For this purpose, an integrated numerical scheme accounting rough contact deformation is developed by coupling modified LSM and DFN modeling for predicting the effective mechanical properties of a realistic outcrop. Smooth joint logic is introduced to consider contact and slip behaviors at fracture surfaces and a modified contact relation to estimating the normal force-displacement from rough contact deformation. Improved constitutive laws are developed and employed for rock matrix and rough fracture surface and implemented in the modified LSM. Complex fracture networks presented by DFNs are automatically extracted based on the gradient Hough transform algorithm. This developed framework is validated by classic equivalent medium theories. It shows the model could be used to emulate naturally-fractured media and to quantitatively investigate the effects of fracture attributes and micro-scale surface roughness on the compression mechanism.
Download