Neutral hydrogen (HI) intensity mapping with single-dish experiments is a powerful approach for probing cosmology in the post-reionization epoch. However, the presence of bright foregrounds over four orders of magnitude stronger than the HI signal makes its extraction highly challenging. While all methods perform well when assuming a Gaussian beam degraded to the worst resolution, most of them degrade significantly in the presence of a more realistic beam model. In this work, we investigate the performance of SDecGMCA. This method extends DecGMCA to spherical data, combining sparse component separation with beam deconvolution. Our goal is to evaluate this method in comparison with established foreground removal techniques, assessing its ability to recover the cosmological HI signal from single-dish intensity mapping observations under varying beam conditions. We use simulated HI signal and foregrounds, covering the frequency ranges relevant to MeerKAT and SKA-Mid. The foreground removal techniques tested fall into two main categories: model-fitting methods (polynomial and parametric) and blind source separation methods (PCA, ICA, GMCA, and SDecGMCA). Their effectiveness is evaluated based on the recovery of the HI angular and frequency power spectra under progressively more realistic beam conditions. While all methods perform adequately under a uniform degraded beam, SDecGMCA remains robust when frequency-dependent beam distortions are introduced. In the oscillating beam case, SDecGMCA suppresses the spurious spectral peak at and achieves accuracy at intermediate angular scales ( ), outperforming other methods. Beam inversion, however, remains intrinsically unstable beyond , setting a practical limit on the method.