An imbibition process of water into a matrix is required to investigate the influences of large-volume fracturing fluids on gas production of unconventional formations. Slip flow has been recognized by recent studies as a major mechanism of fluid transport in nanotubes. For nanopores in shale, a slip boundary is nonnegligible in the imbibition process. In this study, we established an analytic equation of spontaneous imbibition considering slip effects in capillaries. A spontaneous imbibition model that couples the analytic equation considering the slip effect was constructed based on fractal theory. We then used a model for various conditions, such as slip boundary, pore structure, and fractal dimension of pore tortuosity, to capture the imbibition characteristics considering the slip effect. A dynamic contact angle was integrated into the modeling. Results of our study verify that the slip boundary influences water imbibition significantly. The imbibition speed is significantly improved when slip length reaches the equivalent diameter of a tube. Therefore, disregarding the slip effect will underestimate the imbibition speed in shale samples.