We propose a new methodology projected for the estimation of the adsorption energy distribution from the monolayer part of a single nitrogen adsorption isotherm determined at 77 K based on the lattice density functional theory (DFT) via the Aranovich-Donohue formalism. At first sight, the presented approach is computationally more difficult than a classical one. However, it is more flexible and comprehensible. Next, we developed a numerical program and used it for the estimation of the adsorption energy distribution from the experimental data on carbon black samples. The main nitrogen molecule-carbon black surface interaction energy can be estimated as approximate to7-8 kJ/mol, but the heterogeneity of the investigated materials differs significantly. Furthermore, we compare the results obtained from the lattice DFT via the Aranovich-Donohue formalism with the solution of the integral equation with the kernel represented by the classical monolayer localized Fowler-Guggenheim isotherm equation. The similarity between these two independent approaches is observed. The proposed methodology can be used for the investigation of the energetic heterogeneity of not only the carbonaceous materials but also the other "flat-surfaced" solids.