Hydrodynamic forces on small diameter subsea pipelines and cables placed near seabed are important for their on-bottom stability design. In offshore environments, these pipelines are usually subjected to extreme wave conditions. The present study investigates hydrodynamic forces acting on a pipeline near a flat seabed subjected to a wave-induced boundary layer flow. The Keulegan-Carpenter numbers (KC) of the wave-induced boundary layer flow are 20, 140 and 200, defined based on the pipeline diameter (D), the maximum velocity of the undisturbed near-bed orbital velocity (Uw) and the period of the incoming oscillatory flow (Tw). Reynolds number is 1 × 104 based on Uw and D. A seabed roughness ratio ks/D (ks is the Nikuradse equivalent sand roughness) of up to 0.1 and different gap ratios of G/D = 0.05~0.5 between the pipeline and the seabed are considered. Numerical simulations have been carried out based on two-dimensional (2D) Unsteady Reynolds-Averaged Navier-Stokes (URANS) equations combined with the k-ω Shear Stress Transport (SST) turbulence model. A preliminary one-dimensional (1D) simulation is carried out to obtain a fully developed wave-induced boundary layer velocity profile, which is used as inlet flow for the 2D simulations. The numerical model is validated against the experimental data reported by Sumer et al. (1991) at KC = 10. Influences of KC, ks/D and G/D on the hydrodynamic forces and the surrounding flows are discussed in detail.