Depleted reservoirs are widely used for underground gas storage due to their advantages of low construction cost and easy development. Under the in uence of complex geological conditions and frequent operations, the integrity of the wells in depleted reservoirs is prone to failure, which would potentially lead to gas leakage. In this study, by using a nite element based computational uid dynamics model, we have developed evaluation criteria for assessing the severity of the occurred wellbore integrity failure and the risk of the un-occurred wellbore integrity failures respectively to identify hazardous zones potentially prone to wellbore integrity failure. The study results indicate that the gas storage wellbore integrity failure is prone to occur inside the wellbore structure in the direction of the minimum ground stress near the lower boundary of the formation interlayer. The wellbore integrity failure hazardous zones are mainly concentrated at the formation interlayer boundaries. The practical guidelines and solutions derived from current research results can be used for environmental protection work to avoid pollution problems caused by natural gas leakage.