Blasingame production decline curve analysis (PDCA) is crucial for extracting natural gas from fractured tight sandstone reservoirs (FTSRs). The natural gas flow in FTSRs is complex due to high-velocity non-Darcy (hvnD) flow, threshold pressure gradient (TPG), and stress sensitivity of reservoirs. These factors pose a challenge in solving nonlinear equations. In this paper, we establish a mathematical model for the flow of natural gas in FTSRs. The fracture is characterized by a hvnD flow coupled using an embedded discrete fracture model (EDFM). The gas flow in the matrix is coupled with the TPG. Stress sensitivity has been considered in both fractures and matrix. The finite volume method (FVM) was employed to discretize and solve the nonlinear equations in the model. The pressure-normalized rate algorithm and gas flow model were used to present a Blasingame production decline curve for fractured horizontal wells in FTSRs. The Blasingame production decline curve can be used to identify natural gas flow regimes, matrix petrophysical properties, and fracture parameters. The effect of the hvnD flow in fractures, the lvnD flow in the matrix (TPG), stress sensitivity, fracture parameters, and production schedule on the Blasingame production decline curve and pressure profile are discussed in detail. This study provides scientific and theoretical support for predicting gas wells performance and the gas production decline law in FTSRs.