It has been observed that the shale gas production modeled with conventional simulators/models is much lower than actually observed field data. Generally reservoir and/or stimulated reservoir volume (SRV) parameters are modified (without much physical support) to match production data. One of the important parameters controlling flow is the effective permeability of the intact shale. In this project we aim to model flow in shale nano pores by capturing the physics behind the actual process. For the flow dynamics, in addition to Darcy flow, the effects of slippage at the boundary of pores and Knudsen diffusion have been included. For the gas source, the compressed gas stored in pore spaces, gas adsorbed at pore walls and gas diffusing from the kerogen have been considered. To imitate the actual scenario, real gas has been considered to model the flow. Partial differential equations were derived capturing the physics and finite difference method was used to solve the coupled differential equations numerically. The contribution of Knudsen diffusion and gas slippage, gas desorption and gas diffusion from kerogen to total production was studied in detail. It was seen that including the additional physics causes significant differences in pressure gradients and increases cumulative production. We conclude that the above effects should be considered while modeling and making production forecasts for shale gas reservoirs.