In order to study the influence of wave acceleration asymmetry on seabed response, a numerical model of wave-induced seabed response was developed based on the Biot theory of poroelasticity, and the model was validated against the measured data from flume experiments. A method of calculating wave pressure boundary conditions on the seabed surface, which considers wave acceleration asymmetry, is introduced. Effects of wave acceleration asymmetry on pore water pressures, effective stresses, and soil liquefaction of the seabed are discussed based on numerical simulation. Results show that wave acceleration asymmetry considerably reduces the peak value of the positive pore water pressure in the surface soil layer of the seabed, modifies the spatial distribution pattern of the pore water pressure, and eventually leads to decreases in both the depth and width of seabed liquefaction.