В работе используются статистические уравнения и формулы для неоднородных многокомпонентных систем, полученные в рамках двухуровневого молекулярно-статистического подхода, который базируется на одновременном применении метода коррелятивных функций Боролюбова – Борна – Грина – Ивона (ББГКИ) и метода условных распределений Ротта, а также метода термодинамических функционалов плотности. В частности, в развиваемом статистическом подходе однокомпонентная система с тепловыми вакансиями рассматривается как гипотетическая двухкомпонентная, которая состоит из реальных частиц сорта а и фиктивных частиц сорта b. Фиктивные частицы используются для учета вкладов от тепловых вакансий в кристаллическом состоянии вещества. Выражения для унарной и бинарной функций распределения молекул неоднородной системы, определяющие ее свободную энергию, в данной работе используются для изучения микроструктуры однородной бинарной системы, которая описывается в рамках классической решеточной модели с парным взаимодействием между атомами сортов а и b. Проведено сопоставление выражения для чисел заполнения в двухуровневом молекулярно-статистическом подходе с ранее полученными вероятностными функциями, которые описывают ближний порядок в решеточных однородных бинарных сплавах с учетом взаимодействия только ближайших соседей.