Abstract: The saturation seepage in the surrounding rock of the deep vertical shaft under high water pressure condition is generally in the non-Darcy seepage state, taking this state of deep vertical shaft grouting and waterproofing combined with steel lining support scheme as the research object, based on the Izbash non-Darcy seepage model, thick-walled cylinder theory and continuous medium elasticity theory, with the help of computer symbolic computing software, the displacement field and stress field of the well perimeter under the action of grouting and steel lining support are analytically solved. The analytical solution and numerical solution are compared by using the finite element software ABAQUS, and the results are basically consistent. Based on the obtained analytical solution, we further studied the effect of the change of non-Darcy seepage coefficient and the change of slurry-steel liner joint support parameters on the support effect. The design optimization method of grouting-steel liner joint support is proposed and calculated by controlling the amount of water inrush. The results showed that the larger the non-Darcy seepage coefficient, the greater the reduction of circumferential stress at the steel well wall with the increase of steel liner thickness and slurry radius. It has a positive effect to enhance the grouting - steel lining joint support effect by the greater non-Darcy seepage coefficient.