Кожухотрубный охладитель: расчёт по LMTD и NTU–ε — пошаговый пример
Тепловой расчёт кожухотрубного охладителя — ключевой этап проектирования. Ошибка в 10% по площади поверхности означает либо перегрев масла или газа, либо переплату за лишний металл. В этой статье разберём два метода расчёта: классический LMTD и современный NTU–ε, покажем пошаговые примеры для маслоохладителя 500 кВт и охладителя агрессивного газа, а также дадим формулы, таблицы и онлайн-калькулятор.
1. Исходные данные для расчёта
Перед выбором метода и началом расчёта необходимо собрать полный набор исходных данных. Неполные данные — главная причина ошибочного подбора охладителей на практике.
Q = G × Cp × (T_in - T_out) Q — тепловая нагрузка (Вт или кВт) — сколько тепла нужно отвести
LMTD = (ΔT₁ - ΔT₂) / ln(ΔT₁ / ΔT₂) LMTD — среднелогарифмический температурный напор (К) — движущая сила теплообмена
G_mass = ρ × G_vol G — массовый расход охлаждаемой среды (кг/с); ρ — плотность (кг/м³)
Минимальный набор исходных данных для расчёта кожухотрубного теплообменника:
- Охлаждаемая среда: состав, расход G (м³/ч или кг/ч), температура входа T₁_in и выхода T₁_out
- Охлаждающая среда: тип (вода/воздух), температура входа T₂_in и допустимая T₂_out
- Рабочее давление обеих сред
- Допустимый перепад давления (ΔP) по каждому контуру
- Теплофизические свойства сред: Cp, ρ, μ, λ при рабочих температурах
- Требования к материалам (агрессивность, нормативная база)
По этим данным находят тепловую нагрузку Q, а затем — требуемую площадь поверхности теплообмена. Подробнее об устройстве и типах охладителей читайте в статье «Кожухотрубный охладитель: принцип работы».
2. Метод LMTD: среднелогарифмический температурный напор
Метод LMTD (Logarithmic Mean Temperature Difference) — базовый инженерный инструмент теплового расчёта. Он применим, когда известны все четыре температуры: вход/выход горячей и холодной сред.
LMTD = (ΔT₁ − ΔT₂) / ln(ΔT₁ / ΔT₂) ΔT₁ = T_hot_in − T_cold_out (разность температур на «горячем» конце)
ΔT₂ = T_hot_out − T_cold_in (разность температур на «холодном» конце)
Формула справедлива для чистого противотока. При прямотоке ΔT₁ = T_hot_in − T_cold_in и ΔT₂ = T_hot_out − T_cold_out
Пример расчёта LMTD для маслоохладителя
Масло охлаждается с 80°C до 45°C, охлаждающая вода нагревается с 20°C до 35°C. Схема — противоток.
- ΔT₁ = 80 − 35 = 45 К (горячий конец: вход масла / выход воды)
- ΔT₂ = 45 − 20 = 25 К (холодный конец: выход масла / вход воды)
- LMTD = (45 − 25) / ln(45/25) = 20 / ln(1,8) = 20 / 0,588 = 34,0 К
Для противоточной схемы LMTD всегда выше, чем для прямотока при тех же температурах. Это означает меньшую площадь и более компактный аппарат. Проектируйте охладители в противотоке — это стандарт для кожухотрубных теплообменников.
Если ΔT₁ ≈ ΔT₂ (разность не более 5%), в знаменателе возникает неопределённость. В этом частном случае LMTD ≈ ΔT₁ = ΔT₂ (арифметическое среднее).
3. Поправочный коэффициент F на многоходовое течение
Реальные кожухотрубные аппараты часто имеют 2, 4 или 6 ходов по трубному пространству. Часть потока при этом движется в прямотоке, снижая эффективность. Для учёта этого в расчётную формулу вводят поправочный коэффициент F:
LMTD_eff = F × LMTD_counterflow F — безразмерный коэффициент, 0 < F ≤ 1,0
F рассчитывают по параметрам R и S:
R = (T₁_in − T₁_out) / (T₂_out − T₂_in) S = (T₂_out − T₂_in) / (T₁_in − T₂_in) Диаграммы F(R, S) приведены в TEMA и ГОСТ 27590
Если расчётное значение F < 0,75, конструкцию следует пересмотреть: увеличить число секций, выбрать чисто противоточную схему или применить метод NTU-ε. Работа при низком F означает неэффективное использование поверхности и высокую чувствительность к отклонениям параметров.
Типовые значения F для кожухотрубных охладителей:
- 1 ход по трубкам (противоток) — F = 1,00
- 2 хода по трубкам — F = 0,88–0,97 (зависит от R и S)
- 4 хода по трубкам — F = 0,82–0,95
- Перекрёстный ток (1 ход) — F = 0,75–0,90
4. Коэффициент теплопередачи k: от чего зависит, типовые значения
Коэффициент теплопередачи k объединяет все тепловые сопротивления — плёнки теплоотдачи с обеих сторон, стенку трубки и загрязнения:
1/k = 1/α₁ + Rf₁ + δ/λ_стен + Rf₂ + 1/α₂ α₁, α₂ — коэффициенты теплоотдачи (Вт/м²·К) с трубной и межтрубной стороны
δ/λ_стен — тепловое сопротивление стенки трубки (δ — толщина, λ — теплопроводность)
Rf₁, Rf₂ — коэффициенты загрязнения (м²·К/Вт) по TEMA
На коэффициент теплоотдачи α влияют скорость потока, вязкость, теплопроводность среды и геометрия (диаметр трубок, шаг перегородок). Для повышения α применяют: увеличение скорости, оребрение, турбулизаторы.
| Охлаждаемая среда | Охладитель | k, Вт/м²·К | Примечание |
|---|---|---|---|
| Вода (чистая) | Вода | 800–1 500 | Стандартные скорости 1–2 м/с |
| Масло (вязкое, >50 сСт) | Вода | 150–280 | Лимитирует вязкостная плёнка масла |
| Масло (маловязкое, <20 сСт) | Вода | 280–450 | Многоходовая компоновка |
| Воздух / газ при атм. давлении | Вода | 30–60 | Межтрубный газ, поперечное обтекание |
| Природный газ (5–20 бар) | Вода | 60–120 | Рост Re с давлением |
| Конденсат / горячая вода | Вода | 600–1 200 | Зависит от скорости и ΔT |
| Аммиак жидкий | Вода | 400–700 | Высокая λ аммиака |
| Кислота (H₂SO₄, HCl разб.) | Вода | 250–500 | Материал — нержавейка 316L / титан |
| Углеводороды C₅–C₈ | Вода | 300–550 | Зависит от вязкости фракции |
Для предварительного расчёта используйте нижнюю границу диапазона k — это заложит запас. Точное значение k получают из расчёта по методу Bell-Delaware (межтрубное пространство) и уравнения Sieder-Tate или Gnielinski (трубное пространство). Программа HTRI делает это автоматически.
5. Учёт загрязнений: коэффициент Rf по TEMA
Загрязнение поверхности — главная причина снижения производительности охладителя в процессе эксплуатации. Стандарт TEMA (Tubular Exchanger Manufacturers Association) устанавливает нормативные значения коэффициентов загрязнения Rf для проектирования.
k_clean = 1 / (1/α₁ + δ/λ + 1/α₂) k_fouled = 1 / (1/k_clean + Rf₁ + Rf₂) Площадь рассчитывают по k_fouled — с учётом будущих загрязнений
Запас поверхности: oversurface = (k_clean / k_fouled − 1) × 100%
Без учёта запаса на загрязнение аппарат не выйдет на расчётный режим охлаждения уже через 6–12 месяцев эксплуатации — по мере накопления биообрастания, накипи и отложений масла. Это критично для компрессорных маслоохладителей и охладителей оборотной воды.
| Среда | Rf, м²·К/Вт | Тип загрязнения | Примечание |
|---|---|---|---|
| Речная вода | 0,000 20 | Биообрастание, накипь | Минимальная скорость 0,9 м/с |
| Морская вода (<52°C) | 0,000 10 | Биообрастание | Трубки CuNi или титан |
| Оборотная охлаждающая вода | 0,000 20 | Накипь, биообрастание | Контроль ХВП обязателен |
| Питьевая водопроводная вода | 0,000 17 | Карбонатная накипь | Умягчение или ингибиторы |
| Дистиллят / деминерализат | 0,000 09 | Минимальное | pH контроль |
| Пар насыщенный | 0,000 09 | Масляные плёнки | Сепаратор пара на входе |
| Масло (минеральное) | 0,000 18 | Коксование, шламы | Трубная сторона |
| Нефтепродукты лёгкие | 0,000 20 | Полимеризация | T < 120°C |
| Нефтепродукты тяжёлые | 0,000 50 | Коксование, асфальтены | T > 150°C — риск закоксовки |
| Природный / попутный газ | 0,000 20 | Конденсат, механика | Сепаратор на входе |
| Аммиак жидкий | 0,000 10 | Масляные плёнки | Маслоотделитель обязателен |
Для охладителей с агрессивными средами (хлориды, кислоты) загрязнение дополняется коррозионным износом. Выбор материала трубок подробно рассмотрен в статье «Материалы и коррозия охладителей».
Сигнал к промывке — рост перепада давления на 25–40% от начального при неизменных расходах. Периодичность химической чистки: при нормальной водоподготовке — 1 раз в 2–3 года, при жёсткой воде — 1 раз в год. Подробнее о качестве воды — в статье «Качество воды для конденсаторов».
6. Требуемая площадь поверхности: F = Q / (k × LMTD)
После определения LMTD, k (с учётом загрязнений) и поправочного коэффициента F рассчитывают требуемую площадь поверхности теплообмена:
F_req = Q / (k_fouled × F_corr × LMTD) F_req — требуемая площадь поверхности (м²)
Q — тепловая нагрузка (Вт)
k_fouled — коэффициент теплопередачи с учётом загрязнений (Вт/м²·К)
F_corr — поправочный коэффициент на схему течения (0,75–1,00)
LMTD — среднелогарифмический температурный напор (К)
К расчётной площади добавляют конструктивный запас 15–25% для стандартных условий и 20–30% для агрессивных сред:
- F_design = F_req × 1,20 (стандарт для маслоохладителей)
- F_design = F_req × 1,25–1,30 (газы, морская вода, агрессивные среды)
Число трубок рассчитывают исходя из стандартных диаметров (25×2, 38×2 мм) и длины секции (2, 3, 4, 6 м). Затем проверяют скорость в трубках (рекомендуется 0,8–2,5 м/с) и перепад давления. Подбор по гидравлике описан в статье «Гидравлика в кожухотрубных аппаратах».
7. Метод NTU–ε: когда применять
Метод NTU–ε (Number of Transfer Units — Effectiveness) решает «обратную задачу»: когда известны площадь, k и расходы, но неизвестны выходные температуры. Также он точнее при сложных схемах течения.
NTU = k × F / C_min где C_min = min(G₁×Cp₁, G₂×Cp₂) — меньший поток теплоёмкости (Вт/К)
C* = C_min / C_max Для противотока эффективность ε:
ε = [1 − exp(−NTU×(1−C*))] / [1 − C*×exp(−NTU×(1−C*))] При C* = 1 (равные потоки теплоёмкости): ε = NTU / (NTU + 1)
Q = ε × C_min × (T_hot_in − T_cold_in) При NTU > 2 метод LMTD с поправочным коэффициентом F даёт погрешность более 15% из-за нелинейности температурных профилей. Используйте NTU–ε при: малых температурных напорах (LMTD < 10 К), многоходовых схемах, неизвестных выходных температурах, перекрёстном токе.
Метод NTU–ε особенно полезен при подборе кожухотрубных конденсаторов с переменной нагрузкой и в случаях, когда расход охладителя варьируется в зависимости от сезона.
Алгоритм NTU для проектировочного расчёта (нахождения площади):
- Задать целевую ε = Q_required / (C_min × ΔT_max)
- Из формулы ε(NTU, C*) найти NTU
- F = NTU × C_min / k — требуемая площадь
8. Пошаговый пример расчёта маслоохладителя 500 кВт
Рассмотрим полный тепловой расчёт маслоохладителя компрессора методом LMTD. Все числа — реальные параметры типичной установки.
Охладитель масла компрессора 500 кВт
Охладитель природного газа 8 бар, 320 кВт
Проверка по методу NTU для маслоохладителя
Убедимся, что метод LMTD применим (NTU < 2):
- C_масло = (53 800/3 600) × 2 000 = 29 900 Вт/К
- C_вода = G_воды × 4 186. При ΔT воды = 16 К и Q = 500 кВт: G_воды = 500 000/(4 186×16) = 7,47 кг/с → C_вода = 31 270 Вт/К
- C_min = 29 900 Вт/К (масло); C* = 29 900/31 270 = 0,956
- NTU = k×F/C_min = 280×69/29 900 = 0,646
- NTU = 0,646 < 2 → метод LMTD применим, погрешность < 3%
9. Пример расчёта охладителя газа (агрессивная среда, запас 1,2)
Рассмотрим охладитель сероводородсодержащего газа — типичный случай для нефтехимии, где требуются повышенные запасы и специальные материалы.
Для газов с H₂S свыше 10 мг/м³ трубки из углеродистой стали подвержены водородному растрескиванию (HIC/SSC). Применяйте трубки из нержавейки 316L или специальных сталей с ограниченной твёрдостью по NACE MR0175. Расчёт без учёта этого фактора недействителен.
Исходные данные: кислый газ (H₂S 0,5%), давление 12 бар, расход 8 500 нм³/ч, охлаждение с 160°C до 55°C, охлаждающая оборотная вода 28°C/45°C.
- Q = G×ρ×Cp×ΔT = (8 500/3 600) × 8,2 × 2 350 × (160−55) = 565 кВт
- LMTD: ΔT₁ = 160−45 = 115 К; ΔT₂ = 55−28 = 27 К → LMTD = (115−27)/ln(115/27) = 88/1,447 = 60,8 К
- k_clean для газ(12 бар)–вода: α_газ = 110, α_вода = 5 500, стенка ст.316L δ=2 мм (λ=16): 1/k = 1/110 + 0,002/16 + 1/5 500 = 0,00909+0,000125+0,000182 = 0,00940 → k_clean = 106 Вт/м²·К
- k_fouled: Rf_газ = 0,0002; Rf_вода_оборот = 0,0002: 1/k_f = 0,00940+0,0004 = 0,00980 → k_fouled = 102 Вт/м²·К
- F_req = 565 000 / (102 × 1,0 × 60,8) = 91,0 м²
- Запас 25% (агрессивная среда, NACE): F_design = 91,0 × 1,25 = 113,8 м² → 115 м²
- Исполнение: ХПГ 1000-115 (плавающая головка, трубки 25×2 из AISI 316L, кожух 09Г2С)
Большой перепад температур между газом (160°C) и кожухом (∼40°C) создаёт температурное расширение, которое может разрушить трубные решётки. Для газовых охладителей с ΔT_конструктивным > 50 К применяют тип ХПГ (плавающая головка) или ХКГ (компенсатор на кожухе). Подробнее — в статье «Горизонтальный, вертикальный охладитель с компенсатором и плавающей головкой».
10. Программы расчёта: HTRI, Aspen EDR, COMPRESS
Для инженерного проектирования кожухотрубных охладителей применяют специализированные программы. Каждая имеет свою нишу применения.
HTRI Xchanger Suite
Наиболее точная программа теплового расчёта. Реализует метод Bell-Delaware для межтрубного пространства и учитывает перетоки, байпасные потоки и нерасчётные режимы. Применяется в нефтегазовой и нефтехимической промышленности. Лицензия — коммерческая, используется на крупных инжиниринговых предприятиях.
Aspen EDR (Exchanger Design & Rating)
Интегрируется с Aspen Plus и HYSYS. Позволяет оптимизировать конструкцию в рамках технологической схемы. Поддерживает расчёт по TEMA, ASME и PED. Особенно удобен при расчёте охладителей в составе технологических установок.
COMPRESS и PV Elite
Программы расчёта на прочность по ASME VIII Div.1 и ГОСТ 34233. Не выполняют тепловой расчёт, но обязательны для проверки прочности кожуха, трубных решёток и фланцев при высоких давлениях. Применяют совместно с HTRI или Aspen EDR.
Excel + ручной расчёт
Для предварительного подбора охладителя при известных k и простой схеме течения достаточно формул LMTD в Excel. Используйте онлайн-калькулятор ниже для быстрой оценки LMTD и площади поверхности.
HTRI или Aspen EDR необходимы при: давлении свыше 25 бар, двухфазных средах, конденсирующих парах, вязкости > 100 сСт, газах высокого давления, системах с несколькими охладителями в обвязке. Для стандартных маслоохладителей 50–200 кВт достаточно расчёта по LMTD с табличными значениями k.
11. Как передать ТЗ инженеру S22: чек-лист исходных данных
Чем полнее исходные данные, тем точнее и быстрее инженер S22 выполнит расчёт и предложит оптимальный охладитель. Воспользуйтесь чек-листом:
- Тепловая нагрузка Q (кВт) — или данные для её расчёта (G, Cp, ΔT)
- Охлаждаемая среда: название, состав (% компонентов), расход G (м³/ч или кг/ч)
- Температуры горячей среды: вход T₁_in и выход T₁_out (°C)
- Охлаждающая среда: тип (вода оборотная / питьевая / морская / воздух)
- Температуры холодной среды: доступная T₂_in и допустимый нагрев T₂_out (°C)
- Рабочее давление горячей и холодной стороны (бар изб.)
- Допустимый перепад давления ΔP (бар) для каждого контура
- Агрессивность: pH, хлориды (мг/л), H₂S (мг/м³), кислоты — при наличии
- Анализ охлаждающей воды (жёсткость, хлориды, взвесь)
- Нормативная база: ГОСТ / ASME / PED / НПАОП
- Ограничения площадки: габариты, горизонтальная/вертикальная установка
- Срок службы и периодичность обслуживания
Передать ТЗ можно через форму ниже — инженер свяжется в течение 2 часов. Типовые сроки расчёта: стандартный охладитель — 1 рабочий день, нестандартный (агрессивные среды, высокое давление) — 2–3 рабочих дня.
Также рекомендуем ознакомиться с маркировками охладителей ХНГ/ХНВ/ХКГ/ХКВ/ХПГ/ХПВ — это поможет грамотно сформулировать техническое задание.