We present a kinetic theory neuronal network model that captures how first and second order statistics of population activity depend on certain statistical patterns of connectivity. To capture pairwise correlations among neurons in each population, we model the activity of each population by a joint population density ρ(x1, x2, t), where (x1, x2) represents the state of any pair of neurons in the population. We derive the evolution equations of ρ(x1, x2, t) using the Chapman-Kolmolgorov formalism based on the assumption that the input to any pair of neurons is a multivariate Poisson point process. We then derive the network equations by coupling population equations. The coupling scheme is based on (1) the first and second order statistics of connectivity, (2) the approximation of delayed correlation as instantaneous correlation in each population activity, and (3) the inference of higher order statistics of population activity from lower order statistics via a conditional independence method. We test the validity of our network equations by comparing the numerical solutions to the Monte Carlo simulation under various network configurations. Finally, we discuss the success and failure of the implementation of our kinetic theory network model.