Wellbore breakout is one of the critical issues in drilling due to the fact that the related problems result in additional costs and impact the drilling scheme severely. However, the majority of such wellbore breakout analyses were based on continuum mechanics. In addition to failure in intact rocks, wellbore breakouts can also be initiated along natural discontinuities, e.g. weak planes and fractures. Furthermore, the conventional models in wellbore breakouts with uniform distribution fractures could not reflect the real drilling situation. This paper presents a fully coupled hydro-mechanical model of the SB-X well in the Tarim Basin, China for evaluating wellbore breakouts in heavily fractured rocks under anisotropic stress states using the distinct element method (DEM) and the discrete fracture network (DFN). The developed model was validated against caliper log measurement, and its stability study was carried out by stress and displacement analyses. A parametric study was performed to investigate the effects of the characteristics of fracture distribution (orientation and length) on borehole stability by sensitivity studies. Simulation results demonstrate that the increase of the standard deviation of orientation when the fracture direction aligns parallel or perpendicular to the principal stress direction aggravates borehole instability. Moreover, an elevation in the average fracture length causes the borehole failure to change from the direction of the minimum in-situ horizontal principal stress (i.e. the direction of wellbore breakouts) towards alternative directions, ultimately leading to the whole wellbore failure. These findings provide theoretical insights for predicting wellbore breakouts in heavily fractured rocks.