Abstract: As an alternative means of experiments at sea, reverberation signal simulation supports the studies on sonar design and anti-reverberation techniques. In this paper, scattering acoustic field of seabed in shallow water is studied by means of normal-mode method, which is treated as a transfer function formed by ocean transmission network. On this basis, scattering region is divided into many cells according to azimuth angles, and the phase differences between seabed cells are taken into account. Then, a simulation method of reverberation sequence received by horizontal array is put forward. Finally, the method is illustrated in the case of shallow water with negative gradient and the results show that the spectrum characteristics, statistical characteristics and spatial correlation characteristics of simulation signals are in accordance with theoretical predictions, which are useful for the research on anti-reverberation techniques.