This study developed a production-forecasting model to replace the numerical simulation and the decline curve analysis using reservoir and hydraulic fracture data in Montney shale gas reservoir, Canada. A shale-gas production curve can be generated if some of the decline parameters such as a peak rate, a decline rate, and a decline exponent are properly estimated based on reservoir and hydraulic fracturing parameters. The production-forecasting model was developed to estimate five decline parameters of a modified hyperbolic decline by using significant reservoir and hydraulic fracture parameters which are derived through the simulation experiments designed by design of experiments and statistical analysis: (1) initial peak rate ( P hyp ), (2) hyperbolic decline rate ( D hyp ), (3) hyperbolic decline exponent ( b hyp ), (4) transition time ( T transition ), and (5) exponential decline rate ( D exp ). Total eight reservoir and hydraulic fracture parameters were selected as significant parameters on five decline parameters from the results of multivariate analysis of variance among 11 reservoir and hydraulic fracture parameters. The models based on the significant parameters had high predicted R 2 values on the cumulative production. The validation results on the 1-, 5-, 10-, and 30-year cumulative production data obtained by the simulation showed a good agreement: R 2 > 0.89 . The developed production-forecasting model can be also applied for the history matching. The mean absolute percentage error on history matching was 5.28% and 6.23% for the forecasting model and numerical simulator, respectively. Therefore, the results from this study can be applied to substitute numerical simulations for the shale reservoirs which have similar properties with the Montney shale gas reservoir.