Confined phase behavior plays a critical role in predicting production from shale reservoirs. In this work, a pseudo-potential lattice Boltzmann method is applied to directly model the phase equilibrium of fluids in nanopores. First, vapor-liquid equilibrium is simulated by capturing the sudden jump on simulated adsorption isotherms in a capillary tube. In addition, effect of pore size distribution on phase equilibrium is evaluated by using a bundle of capillary tubes of various sizes. Simulated coexistence curves indicate that an effective pore size can be used to account for the effects of pore size distribution on confined phase behavior. With simulated coexistence curves from pore-scale simulation, a modified equation of state is built and applied to model the thermodynamic phase diagram of shale oil. Shifted critical properties and suppressed bubble points are observed when effects of confinement is considered. The compositional simulation shows that both predicted oil and gas production will be higher if the modified equation of state is implemented. Results are compared with those using methods of capillary pressure and critical shift.