The morphological characteristics and distribution range of the plastic zone around tunnels, influenced by the true three-dimensional high geostress environment in deep engineering, constitute fundamental reference data for both tunnel design and the control of surrounding rock stability. Based on the GZZ three-dimensional rock mass strength criteria and the complex variable method, a novel semi-analytical method for solving the plastic zone of a circular tunnel subjected to non-hydrostatic stress is presented. The conformal mapping is employed to intricately transform the complex boundaries of plastic zones characterized by diverse shapes onto the unit circle boundary within the image plane. This framework proves invaluable for the precise analytical process of the complex plastic zone morphology and the plastic influence range within the surrounding rock under complex stress conditions. In contrast to traditional analytical methods, the efficacy and versatility of our proposed approach are substantiated. Notably, it demonstrates unique advantages, particularly in the context of non-hydrostatic stress fields and varied rock mass parameters, with respect to predicting both surrounding rock stability and plastic zone morphology. Through the analysis of numerous cases, the study reveals the systematic patterns of influence that in-situ stress and surrounding rock characteristic parameters have on the shape and maximum depth of the plastic zone. The results indicate that (1) there are notable disparities in the morphology and extent of the plastic zone between the GZZ and GHB criteria; (2) the transition of the plastic zone in surrounding rock from U-shaped to V-shaped is significantly influenced by the geostress-to-strength ratio Kσ; (3) there is a notable linear correlation observed between the maximum plastic zone thickness rmax and both Kσ and the stress concentration factor KSCF, importantly, this correlation is independent of far-field stress and the brittleness characteristics of the surrounding rock. These findings provide a novel discriminative method for the rapid analysis of surrounding rock stability in deep circular tunnels.