Nuclear magnetic resonance (NMR) spectroscopy is regarded as one of the most powerful and versatile analytical approaches to assure the quality of heparin preparations. In particular, it was recently demonstrated that by using ¹H NMR coupled with chemometrics heparin and low molecular weight heparin (LMWH) samples derived from three major animal species (porcine, ovine and bovine) can be differentiated [Y.B. Monakhova et al. J. Pharm. Anal. 149 (2018) 114–119]. In this study, significant improvement of existing chemometric models was achieved by switching to 2D NMR experiments (heteronuclear multiple-quantum correlation (HMQC) and diffusion-ordered spectroscopy (DOSY)). Two representative data sets (sixty-nine heparin and twenty-two LMWH) belonged to different batches and distributed by different commercial companies were investigated. A trend for animal species differentiation was observed in the principal component analysis (PCA) score plot built based on the DOSY data. A superior model was constructed using HMQC experiments, where individual heparin (LMWH) clusters as well as their blends were clearly differentiated. The predictive power of different classification methods as well as unsupervised techniques (independent components analysis, ICA) clearly proved applicability of the model for routine heparin and LMWH analysis. The switch from 1D to 2D NMR techniques provides a wealth of additional information, which is beneficial for multivariate modeling of NMR spectroscopic data for heparin preparations.