This work is focused on developing a computational model to predict the production rate and pressure evolution of natural gas from unconventional reservoirs, particularly shale gas deposits. The model is based on the principle of conservation of mechanical energy and was developed from the transient solution of Bernoulli’s equation. This solution was obtained by computing the pressure evolution in the well resulting from the combined action of extracting the free gas and of gasification from kerogen. The transient behavior of gas production by hydraulic fracturing was calculated by numerically integrating Bernoulli’s equation. The curves representing gas flow evolution were considered as a series of stepwise steady states under a constant gas flow rate, similar to the pressure–time curves. These time steps were connected by instantaneous drops in pressure or gas flow rates. On the other hand, the delayed release of the adsorbed and dissolved gas in the kerogen was accurately calculated by introducing a semi-empirical gas pressure source term into the gas well pressure equation. The effect of this source is to gradually increase the gas pressure in the reservoir, emulating the gas release mechanisms from the organic matter. Model validation was based on production data from the unconventional reservoirs Eagle Ford, U.S.A., and Burgos basin, México. The initial measured gas production rate was used to determine a global friction factor of the gas flowing out from soil cracks and ducts. Additionally, measured production rate data were used to determine the coefficients of the source term function. Pearson correlation coefficients of 0.97 and 0.96 were obtained for Eagle Ford and Burgos basins data, respectively. In contrast, the corresponding coefficients calculated from the traditional Arps’ model were 0.89 and 0.5, respectively. The present pressure source model (PSM) represents a new approach to characterize the process of gas production from unconventional reservoirs, proving to be accurate in forecasting both the gas flow rate and pressure evolution during gas production. The postulated pressure source term was shown to mimic the desorption and diffusion kinetics, which release free gas from the kerogen.