A series of numerical model tests were performed to investigate the behaviour of the anisotropic rock surrounding circular excavations under high confining pressures. The aim was to provide information on the formation of fractures and failure around deep level rock tunnels under controlled conditions. Solid cubes containing a circular hole were confined to a vertical pressure with same as the confinement in the horizontal directions. In this modeling, the inhomogeneous rock is generated by using Weibull parameters which are related to the microstructural properties determined by crack size distribution and grain size. The fracture angle is assumed to be 45o. The observed failure zone around the excavation was simulated using both the maximum tensile strain criterion and Mohr-Coulomb criterion respectively (as the damage threshold). And RFPA (Realistic Failure Process Analysis) code was used as the calculating tool in this modelling, three opening modes are simulated and compared. Computational model predictions that include crack propagation and failure modes of rock show a good agreement with those of the observation in site. It is pointed out that the damage evolution of EDZ strongly depends on the inhomogeneous, the excavation mode, anisotropic property, and the various loading conditions. Concerning the existence of a weak plane, the amount of displacement at the side wall of the tunnel was quite large, since the shear deformation occurred in EDZ. The model is implemented in RFPA code and is able to represent the change in fracture patterns between the solid and jointed parts. This provides confidence for the application of the numerical model to the design of rock tunnels at great depth.