The liquid slippage behavior due to molecular interactions at fluid-solid (F-S) interfaces is of great importance to the transport of shale oil and clay water. A mesoscopic single-phase lattice Boltzmann method (LBM), based on a continuous and exponentially decaying F-S interaction force and a midgrid bounce-back boundary condition, is proposed to be responsible for the apparent liquid slippage. The F-S interaction force is established at the particle level and thus can be readily extended to porous media. When it is repulsive (attractive), the phenomena of positive (negative) slip lengths and fluid slip (damping) are successfully recovered. This model is validated by the velocity profiles on hydrophobic and hydrophilic surfaces in a benchmark microchannel flow experiment. The slip length is found to be independent of shear rate (its constituents including body force, pore diameter, and kinematic viscosity), but dependent on pore geometry (smaller in porous media than in capillary tubes). Both slip length and permeability ratio follow a power law relationship with interaction parameters (strength and decay length) in capillary tubes and porous media. The permeability ratio estimated analytically with the slip length considered agrees well with that calculated from the LBM simulations, except for the fluid slip in porous media with a significant overestimation. The estimated permeability ratio indicates that it increases (decreases) nonlinearly as the pore diameter decreases, suggesting the great importance of the F-S interaction particularly for thin capillary tubes and microporous media.