In view of the multiple compound viscoelastic wave flaws, the discrete singular convolution scheme has been adopted to simulate this model. The detailed discrete formula and the treatment method of boundary condition has been proposed, with the validity of the method being verified by numerical examples. It is demonstrated that the discrete singular convolution method can efficiently deal with singularity and high order space derivatives in the multiple compound viscoelastic waves, with even a more controllable accuracy.