В связи с интенсификацией строительства в районах распространения лессовых просадочных грунтов вопросы прогнозирования развития процессов подтопления являются актуальной задачей, т.к. могут привести к аварийному замачиванию, неравномерному подъему горизонта грунтовых вод, изменению напряженно-деформированного состояния грунтов и, соответственно, к потере пригодности эксплуатации здания или сооружения. Лессовый грунт обладает ярко выраженной фильтрационной анизотропией. Просадка, фильтрация воды происходят в условиях неполного водонасыщения. Появление новых компьютерных технологий позволяет совершенствовать методы математического моделирования и разрабатывать математические модели численными методами, достоверно отражающими внутрипочвенные процессы. Настоящая статья посвящена совершенствованию математической модели задачи влагопереноса для неоднородных фильтрационно-анизотропных лессовых грунтов с учетом их структурных особенностей.