A new method based on generalized reflection and transmission (R/T) coefficients method is proposed to calculate the single seismic phase (SSP) of cylindrically mlultilayered media including liquid intelayer. The use of nor-malization facotrs and normalzed Lame coefficients makes the algorithm stable numerically. Using the modified R/T matrices, we derive the iterative expressions of generalized R/T matrices, and by using the iterative relation we determine the SSP of each interface and the full waveforms. To show the superiority of this new approach for investigating of reflection and transmission properties of cylindrically multivalyered media, we simulate the full waveforms and SSPs of cased hole model with annulus 1 (casing-cement interface) channelling (or, cross-flow). The generalized reflection coefficient spectra and SSPs of interfaces obtained show the propagation mechanism of each component of full waveform clearly.