The production decline model for conventional gas reservoirs is widely used for dynamic analysis of shale gas reservoirs at home and abroad, but due to some special features of the occurrence of the shale gas involving mainly adsorption state and free state, the regressive analysis methods of conventional gas reservoirs are no longer valid. Based on the traditional dual-porosity model, a mathematical model with consideration of the stress sensitivity and the sorption effects is established by introducing the non-dimensional stress sensitive coefficient and defining two non-dimensional adsorption parameters--dimensionless Langmuir volume and dimensionless Langmuir pressure. The numerical solution of the production is obtained together with the production decline Blasingame type curves through the numerical difference and Newton iterative methods. It is shown that the higher the stress sensitivity, the smaller the non-dimensional production, the production integral and the production decline integral derivative will be. The larger the adsorption of the shale gas matrix, the more delayed the point will be where the non-dimensional production decline integral derivative curve becomes concave down. At the same time, with the help of the shale gas well production data, the shale gas reservoir adsorption constants, other formation parameters, and the stress sensitive effect caused by the decline in the production can be obtained by typical curves.