A 2-D fully nonlinear numerical wave tank (NWT) is developed to investigate the Bragg resonance scattered by submerged multi-array breakwaters for random waves. This model is based on a boundary integral equation method with linear element scheme. The incident random waves are generated by JONSWAP spectrum at one end of the wave tank. Results indicate that the reflection coefficient of random waves at the first peak of resonance is smaller than that of the regular waves, but the band width of frequency of Bragg effect has become wider and this advantage may compensate the peak reduction. Moreover, it shows that the Bragg reflection phenomenon will increase as the increase of the following parameters such as the relative bar height, the ratio of width between top and bottom of the bar, the number of bars, and the toe angle of submerged breakwaters.