A fast iterative approach of generalized forward-backward method with the spectrum acceleration algorithm (GFBM/SAA) is developed for solving electromagnetic field integral equation and is applied to numerical simulation of radar surveillance of the ship target and oceanic clutters in two-dimensional (2-D) model. Randomly rough surface is realized by the Monte Carlo method using the 1-D Pierson-Moskowitz spectrum. Numerical results of bistatic and backscattering from the ship target and oceanic clutters demonstrate the functional dependencies upon the situation of radar surveillance such as the incidence and observation angles, radar altitude, ship RCS, and other oceanic conditions.