Приведены результаты математического моделирования кондуктивного тепломассообмена в приповерхностных геотермальных системах. Теплообменник системы описывается как вертикальный канал в непроницаемом породном массиве. Для описания температуры пород используется уравнение теплопроводности в цилиндрических координатах. На границе пород и теплообменника выполняются условия III рода.