Production of gas out of low permeability shale packages is very recent in the Western Canadian Sedimentary Basin (WCSB). The process of gas release and production from shale gas sediments is not well understood. Because of adsorptive capacity of certain shale constituents, including organic carbon content,: coalbed methane models are sometimes being applied to model and simulate tight shale gas production behaviour. Alternatively, conventional Darcy flow models are sometimes applied to tight shale gas. However, neither of these approaches takes into account the differences in transport mechanisms in shale due to additional nanopore networks. Hence, the application of existing models for shale results in erroneous evaluation and predictions. Our analysis shows that a combination of a nanopore network connected to a micrometre pore network controls the gas flow in shale. Mathematical modelling of gas flow in nanopores is difficult since the standard assumption of no-slip boundary conditions in the Navier-Stokes equation breaks down at the nanometre scale, while the computational times of applicable molecular-dynamics (MD) codes become exorbitant. We found that the gas flow in nanopores of the shale can be modelled with a diffusive transport regime with a constant diffusion coefficient and negligible viscous effects. The obtained diffusion coefficient is consistent with the Knudsen diffusivity which supports the slip: boundary condition at the nanopore surfaces. This model can be used for shale gas evaluation and production optimization.