Understanding the flow characteristics of shale gas especially in nanopores is extremely important for the exploitation. Here, we perform molecular dynamics (MD) simulations to investigate the hydrodynamics of methane in nanometre-sized slit pores. Using equilibrium molecular dynamics (EMD), the static properties including density distribution and self-diffusion coefficient of the confined methane are firstly analyzed. For a 6 nm slit pore, it is found that methane molecules in the adsorbed layer diffuse more slowly than those in the bulk. Using nonequilibrium molecular dynamics (NEMD), the pressure-driven flow behavior of methane in nanopores is investigated. The results show that velocity profiles manifest an obvious dependence on the pore width and they translate from parabolic flow to plug flow when the width is decreased. In relatively large pores (6 – 10 nm), the parabolic flow can be described by the Navier-Stokes (NS) equation with appropriate boundary conditions because of its slip flow characteristic. Based on this equation, corresponding parameters such as viscosity and slip length are determined. Whereas, in small pores (∼ 2 nm), the velocity profile in the center exhibits a uniform tendency (plug flow) and that near the wall displays a linear increase due to the enhanced mechanism of surface diffusion. Furthermore, the profile is analyzed and fitted by a piecewise function. Under this condition, surface diffusion is found to be the root of this anomalous flow characteristic, which can be negligible in large pores. The essential tendency of our simulation results may be significant for revealing flow mechanisms at nanoscale and estimating the production accurately.