In this paper, considering on the real threshold pressure gradient (TPG), the mathematical model for low velocity non-Darcy unstable gas flow was established for the tight gas reservoir. Based on the low velocity non-Darcy unstable gas flow model, the gas flow model in heterogeneity formation was developed, and the analytical solution was obtained on the condition of constant production for inner boundary. Then we calculated the model numerically by Matlab programming. Its analyzed that the pressure distribution and the influencing factors such as production rate, producing time, reservoir thickness and the threshold pressure gradient. We compared the heterogeneity formation with the homogeneous formation by the heterogeneity low velocity non-Darcy unstable gas flow model, which shows great difference and illustrate the significance of research on heterogeneity.