In this paper, a multi-relaxation-time lattice Boltzmann (LB) model for nanoscale liquid flow is established to investigate the liquid flow characteristics in nanoporous media. The slip length and effective viscosity obtained by molecular dynamics (MD) simulations are adopted to account for the nanoscale effect. First, the LB model for water flow in nanopores is built and water flow characteristics in nanoporous media are investigated. The fluid-solid interaction force shows significant influence on water flow in nanoporous media. The water flux through nanoporous media increases with the decrease of fluid-solid interaction force. In hydrophobic nanoporous media, the nanoscale effect can decrease the gap between water flow resistance in large pores and that in small pores, making the velocity distribution more uniform. In addition, the end effect caused by the bending of streamlines can induce significant additional flow resistance. Neglecting the end effect can greatly overestimate water flow ability. The pore structure also has significant influence on water flow in nanoporous media. With the increase of specific interfacial length, the nanoscale effect increases. Finally, the LB model for oil (octane) flow in quartz nanopores is also established by incorporating the MD simulation results. Oil flow simulation in quartz nanoporous media shows that the conclusions obtained for water flow are also applicable for oil flow.