In recent decades, increasing attention and efforts are focused on developing accurate models to describe general power law frequency-dependent acoustic attenuation.[1][10][13][14][15][16][17] Most of these recent frequency-dependent models are established via the analysis of the complex wave number and are then extended to transient wave propagation.[18] The multiple relaxation model considers the power law viscosity underlying different molecular relaxation processes.[16] Szabo[8] proposed a time convolution integral dissipative acoustic wave equation.