Wellbores are largely constructed during coal mining, shale gas production, and geothermal exploration. Studying the shape and size of the disturbed zone in surrounding rock is of great significance for wellbore stability control. In this paper, a theoretical model for elastic-plastic-damage analysis around a deep circular wellbore under non-uniform compression is proposed. Based on the elastoplastic softening constitutive model and Mohr-Coulomb strength criterion, the analytical expressions of stresses in the elastic, plastic and damaged zones around a circle wellbore are derived. Further, the boundary line equations among the three zones are obtained according to the conditions of stress continuity. Then, the influence rules of non-uniform in-situ stress and mechanical parameters on the stress distribution and plastic zone size in surrounding rock mass are analyzed. The plastic and the damaged zones are both approximately elliptical in shape. When the lateral stress coefficient of the in-situ stress field takes the value 1, the model degenerates into the Yuan Wenbo’s Solution. If the brittleness coefficient of the surrounding rock is 0, the model degenerates into the Kastner’s Equation. Finally, the results are compared with those under two special cases (in the elastoplastic softening rock under a uniform stress field, in the ideal elastoplastic rock under a non-uniform stress field) and a common approximation method (the perturbation method).