Multi-fractured horizontal wells (MFHW) is one of the most effective technologies to develop tight gas reservoirs. The gas seepage from tight formations in MFHW can be divided into three stages: early stage with high productivity, transitional stage with declined productivity, and final stage with stable productivity. Considering the characteristics and mechanisms of porous flows in different regions and at different stages, we derive three coupled equations, namely the equations of porous flow from matrix to fracture, from fracture to near wellbore region, and from new wellbore region to wellbore then an unstable productivity prediction model for a MFHW in a tight gas reservoir is well established. Then, the reliability of this new model, which considers the multi-fracture interference, is verified using a commercial simulator (CMG). Finally, using this transient productivity prediction model, the sensitivity of horizontal well's productivity to several relevant factors is analyzed. The results illustrate that threshold pressure gradient has the most significant influence on well productivity, followed by stress sensitivity, turbulence flow, and slippage flow. To summarize, the proposed model has demonstrated a potential practical usage to predict the productivity of multi-stage fractured horizontal wells and to analyze the effects of certain factors on gas production in tight gas reservoirs.Energies 2020, 13, 1066 2 of 20 of tight gas reservoirs has made considerable breakthrough and become the primary growth point for oil and gas production [1][2][3][4][5][6]. Because of its strong heterogeneity, complicated microscopic pore structure, and poor connectivity of the effective sand bodies, the porous flow mechanism of tight gas reservoirs is obviously distinct from those of conventional reservoirs. Moreover, because of poor reservoir physical properties, small discharging radius, lack of or low natural productivity, high development difficulty, and some other potential disadvantages, a MFHW with hydraulic stimulations is extensively used for developing tight gas reservoirs [7][8][9][10][11][12][13][14][15][16].To date, many scholars are dedicated to the transient productivity of fractured wells in unconventional hydrocarbon reservoirs, considering the different characteristics of nonlinear flows compared to multi-scale flows. The methods to calculate the transient productivity of fractured horizontal wells primarily include analytical and numerical solutions, among which the analytical solution includes the complex potential theory, conformal transformation, and the equivalent flow resistance method [14,17,18]. Research on the transient productivity of tight gas reservoirs is increasingly wide. Through an iterative algorithm, Zheng et al. [19] quantified the correlation between saturation and pressure of the infinitesimal coal. They used the Corey relative permeability model to describe relative gas/water permeability as a function of the pressure. By applying the inter-porosity flow equation based on a pseudo-pr...