Methane needle-plane discharge has practical application prospect and scientific research significance since methane conversion heavy oil hydrogenation is formed by coupling methane needle-plane discharge with heavy oil hydrogenation, which can achieve high-efficient heavy oil hydrogenation and increase the yields of high value-added light olefins. In this paper, a two-dimensional fluid model is built up for numerically simulating the methane needle-plane discharge plasma at atmospheric pressure. Spatial and axial distributions of electric intensity, electron temperature and particle densities are obtained. Reaction yields are summarized and crucial pathways to produce various kinds of charged and neutral particles are found out. Simulation results indicate that axial evolutions of CH3+ and CH4+ densities, electric intensity and electron temperature are similar and closely related. The CH5+ and C2H5+ densities first increase and then decrease along the axial direction. The CH3 and H densities have nearly identical spatial and axial distributions. Particle density distributions of CH2, C2H4 and C2H5 are obviously different in the area near the cathode but comparatively resemblant in the positive column region. The CH3+ and CH4+ are produced by electron impact ionizations between electrons and CH4. The CH5+ and C2H5+ are respectively generated by molecular impact dissociations between CH3+ and CH4 and between CH4+ and CH4. Electron impact decomposition between electrons and CH4 is a dominated reaction to produce CH3, CH2, CH and H. The reactions between CH2 and CH4 and between electrons and C2H4 are critical pathways to produce C2H4 and C2H2, respectively. In addition, the yields of electron impact decomposition reactions between electrons and CH4 and reactions between CH2 and CH4 account for 52.15% and 47.85% of total yields of H2 respectively.