Sublevel caving (SLC) is a mass mining method based upon the utilization of gravity flow of blasted ore and caved waste rock. Drilling and blasting are the initial and the major impact on primary fragmentation and later material flow characteristics. The heterogeneity of rock mass has a great impact on materials’ mechanical properties and plays a significant role in rock fragmentation by blasting. In this work, a numerical model for sublevel caving based on 3D Voronoi tessellations was developed to section rock mass into arbitrarily shaped regions for modelling the heterogeneous rock mass. A coupled finite element and discrete element method was employed to investigate both the fragmentation due to blasting and the gravity flow of ore in a heterogeneous orebody. In the coupled model, the caved waste rock was modelled with loose discrete particles while the orebody to be blasted was simulated with bonded particles. The detonation of explosives was described with a particle blast method and a finite element model was used to model the remaining orebody. The results indicated that the fragmentation at the upper part of the ring is coarse while it is fine at the lower part of the ring. The heterogeneity of rock mass induced some boulders after blasting.