The valence photoemission spectra of alkali metals exhibit multiple plasmon satellite structure. The calculated spectral functions within the GW approximation show only one plasmon satellite at too large binding energy. In this Letter we use the cumulant expansion approach to obtain the spectral functions of Na and Al from ab initio calculations including the effects of band structure. The GW spectral functions are dramatically improved and the positions of the multiple plasmon satellites are in very good agreement with experiment while their intensities cannot be explained from intrinsic effects only.