A novel CDEM method was used for numerical simulation of hydraulic fracturing.
3D initiation and propagation of cracks were investigated using heterogeneous model.
Material heterogeneity has considerable effect on crack initiation.
Crack propagation is controlled primarily by geostress difference.