The intricate morphology of porous media can considerably affect the progression of penetration and the diffusion of grouting slurry. In this study, a Newtonian fluid’s refined seepage motion equation was proposed to incorporate the impact of tortuosity on the grouting process into gravel soil, and the column-hemispherical penetration grouting mechanism was analyzed using the method of theoretical analysis. Utilizing secondary-development programming techniques, a numerical simulation program was developed with COMSOL Multiphysics to model the penetration grouting mechanism of a Newtonian fluid in a column-hemispherical pattern, considering a medium’s tortuosity. The penetration grouting process of Newtonian cement into gravel soils was then simulated. Finally, the theoretical analysis, experimental values, and numerically simulated values were compared.The findings suggested that incorporating the tortuosity of porous media is more efficacious in depicting the penetration and diffusion behavior for Newtonian fluid grouting in porous media, as compared to omitting the tortuosity. The findings of this study contribute to a better understanding of grouting engineering in porous media strata, guiding practical design and construction.