Numerical integration method for seabed response to water waves