This study investigates the dynamic response of fractional poroviscoelastic layered media subjected to a moving harmonic load. The material viscosity includes two parts: the flow-dependent viscosity between the skeleton and pore water and the skeleton viscosity. Firstly, the governing equations of poroelastic media are derived based on Biot's theory and then transformed into ordinary differential equations by using the double Fourier integral transform. By incorporating fractional calculus theory and the dynamic correspondence principle, the applicability of these equations is extended to fractional poroviscoelastic media. The extended precise integration method (PIM) is employed to obtain solutions. Subsequently, the fractional Zener model and the presented method are validated. Finally, parametric studies are conducted to analyze the effects of fractional order, moving velocity, load frequency, and stratification on the dynamic behaviors of the media.