Ocean Bottom Seismometers (OBS) placed on the seafloor surface are utilized for measuring the ocean bottom seismic waves. The vibration of OBS excited by underwater noise on its surface may interfere with its measured results of seismic waves. In this particular study, an OBS was placed on the seabed, while ray acoustic theory was used to deduce the sound field distribution around the OBS. Then using this information, the analytical expression for the OBS vibration velocity was obtained in order to find various factors affecting its amplitude. The finite element computing software COMSOL Multiphysics® (COMSOL) was used to obtain the vibration response model of the OBS which was exposed to underwater noise. The vibration velocity for the OBS calculated by COMSOL agreed with the theoretical result. Moreover, the vibration velocity of OBS with different densities, shapes, and characters were investigated as well. An OBS with hemispherical shape, consistent average density as that of the seafloor, and a physical structure of double tank has displayed minimum amplitude of vibration velocity. The proposed COMSOL model predicted the impact of underwater noise while detecting the ocean bottom seismic waves with the OBS. In addition, it provides significant help for the design and optimization of an appropriate OBS.