Содержание
1. Физические основы мембранного разделения
Мембранное разделение — процесс разделения жидких смесей под действием перепада давления с использованием полупроницаемых мембран, которые селективно пропускают растворитель (воду) и задерживают растворённые вещества (соли, органику, микроорганизмы).
1.1. Явление осмоса и обратного осмоса
Осмос — самопроизвольный переход растворителя через полупроницаемую мембрану из области с меньшей концентрацией растворённого вещества в область с большей концентрацией. Движущей силой является разность химических потенциалов растворителя по обе стороны мембраны.
Давление, при котором осмотический поток прекращается, называется осмотическим давлением (π). Для разбавленных растворов оно описывается уравнением Вант-Гоффа:
π = i × C × R × T
где:
- i — изотонический коэффициент (число частиц при диссоциации: для NaCl i = 2, для CaCl₂ i = 3);
- C — молярная концентрация, моль/л;
- R — универсальная газовая постоянная, 0.0831 л·бар/(моль·К);
- T — абсолютная температура, К.
Обратный осмос — процесс, при котором к концентрированному раствору прикладывается давление P > π, что заставляет растворитель двигаться в обратном направлении (от концентрированного раствора к разбавленному). Это основа работы систем обратного осмоса (RO).
1.2. Уравнения массопереноса
Поток пермеата (воды) через мембрану описывается уравнением:
Jw = A × (ΔP - Δπ)
где:
- Jw — удельный поток воды (flux), л/(м²·ч) или LMH;
- A — коэффициент проницаемости мембраны для воды, л/(м²·ч·бар);
- ΔP — приложенный перепад давления, бар;
- Δπ — разность осмотических давлений по обе стороны мембраны, бар.
Поток солей через мембрану:
Js = B × ΔC
где:
- Js — удельный поток соли, г/(м²·ч);
- B — коэффициент проницаемости мембраны для соли, л/(м²·ч);
- ΔC — разность концентраций соли по обе стороны мембраны, г/л.
Селективность (солеразделение) мембраны:
R = (1 - Cp / Cf) × 100%
где:
- R — селективность, %;
- Cp — концентрация соли в пермеате, мг/л;
- Cf — концентрация соли в питании, мг/л.
2. Классификация мембранных процессов
| Параметр | Микрофильтрация (MF) | Ультрафильтрация (UF) | Нанофильтрация (NF) | Обратный осмос (RO) |
|---|---|---|---|---|
| Размер пор, нм | 50 – 500 | 2 – 50 | 0.5 – 2 | < 0.5 |
| MWCO (молекулярная масса отсечки), Да | 100 000 – 1 000 000 | 1 000 – 100 000 | 150 – 1 000 | < 100 |
| Рабочее давление, бар | 0.5 – 2 | 1 – 5 | 5 – 15 | 10 – 80 |
| Удаляемые вещества | Взвеси, бактерии, простейшие | Вирусы, коллоиды, белки, полимеры | Двухвалентные ионы (Ca²⁺, Mg²⁺, SO₄²⁻), органика MW > 200 | Все ионы, включая одновалентные (Na⁺, Cl⁻), органика MW > 100 |
| Селективность по NaCl, % | 0 | 0 | 20 – 80 | 95 – 99.8 |
| Типичный flux, LMH | 100 – 500 | 50 – 200 | 15 – 30 | 15 – 25 |
| Применение в водоподготовке | Предподготовка RO, удаление мутности | Предподготовка RO, удаление SDI | Умягчение, удаление органики, частичное обессоливание | Глубокое обессоливание, получение питьевой и ультрачистой воды |
Нанофильтрация селективна преимущественно к двухвалентным ионам (Ca²⁺, Mg²⁺, SO₄²⁻) и органике с MW > 200 Да. Одновалентные ионы (Na⁺, Cl⁻) проходят через NF-мембрану частично (селективность 20–80%). Это позволяет работать при более низком давлении (5–15 бар против 10–80 бар для RO) и получать пермеат с умеренной минерализацией, что полезно для питьевой воды.
3. Ключевые параметры расчёта мембранных систем
3.1. Удельная производительность (flux)
Flux — объём пермеата, проходящего через единицу площади мембраны в единицу времени:
J = Qp / Smembrane
где:
- J — flux, л/(м²·ч) или LMH (литров на м² в час);
- Qp — производительность по пермеату, л/ч;
- Smembrane — активная площадь мембраны, м².
Типичные диапазоны flux для различных мембран (при 25 °C):
| Тип мембраны | Типичный flux, LMH | Максимальный flux, LMH | Примечание |
|---|---|---|---|
| RO для солоноватых вод (BW) | 18 – 25 | 30 | При превышении — ускоренный fouling |
| RO для морской воды (SW) | 14 – 20 | 25 | Ниже из-за высокого осмотического давления |
| NF | 15 – 30 | 40 | Зависит от селективности |
| UF | 50 – 150 | 250 | При dead-end режиме |
Превышение максимально допустимого flux приводит к ускоренному уплотнению осадка на поверхности мембраны (concentration polarization) и необратимому fouling. Производители мембран строго нормируют максимальный flux для каждой модели — его нельзя превышать ни при каких условиях.
3.2. Степень извлечения (Recovery)
Recovery — доля питательной воды, преобразованной в пермеат:
Y = Qp / Qf × 100%
где:
- Y — recovery, %;
- Qp — расход пермеата, м³/ч;
- Qf — расход питания (feed), м³/ч.
Расход концентрата:
Qc = Qf - Qp = Qf × (1 - Y)
3.3. Коэффициент концентрации (Concentration Factor, CF)
CF показывает, во сколько раз концентрация солей в концентрате выше, чем в питании:
CF = 1 / (1 - Y)
| Recovery, % | CF | Примечание |
|---|---|---|
| 50 | 2.0 | Типично для SWRO |
| 75 | 4.0 | Типично для BWRO |
| 85 | 6.7 | Высокое значение, риск солеотложения |
| 90 | 10.0 | Только для чистой воды или с антискалантом |
3.4. Концентрация солей в концентрате
Средняя концентрация солей в концентрате рассчитывается с учётом селективности мембраны:
Cc = Cf × CF × (1 - R/100) / (1 - Y × R/100)
Для упрощённого расчёта при высокой селективности (R > 95%):
Cc ≈ Cf × CF
Для точного расчёта осмотического давления применяется средняя концентрация между входом и выходом корпуса:
Cavg = (Cf + Cc) / 2
4. Расчёт осмотического давления
4.1. Упрощённый расчёт через TDS
Для пресных и солоноватых вод при TDS < 10 г/л осмотическое давление приближённо рассчитывается:
π ≈ 0.0075 × TDS × (1 + 0.02 × (T - 25))
где:
- π — осмотическое давление, бар;
- TDS — общая минерализация, мг/л;
- T — температура воды, °C.
Более точная формула (при T = 25 °C):
π ≈ 0.0075 × TDS
Пример: TDS = 2000 мг/л → π ≈ 0.0075 × 2000 = 15 бар
4.2. Точный расчёт через ионный состав
Для точного расчёта применяется уравнение Вант-Гоффа с учётом ионного состава:
π = R × T × Σ(ij × Cj)
где:
- R = 0.0831 л·бар/(моль·К);
- T — температура, К;
- ij — изотонический коэффициент j-го иона;
- Cj — молярная концентрация j-го иона, моль/л.
| Ион | Молярная масса, г/моль | Изотонический коэффициент i |
|---|---|---|
| Na⁺ | 22.99 | 1 |
| K⁺ | 39.10 | 1 |
| Ca²⁺ | 40.08 | 1 (эффективно) |
| Mg²⁺ | 24.31 | 1 (эффективно) |
| Cl⁻ | 35.45 | 1 |
| SO₄²⁻ | 96.06 | 1 (эффективно) |
| HCO₃⁻ | 61.02 | 1 |
| NaCl (полностью) | 58.44 | 2 |
| CaCl₂ (полностью) | 110.98 | 3 |
Осмотическое давление прямо пропорционально абсолютной температуре. При снижении температуры с 25 °C до 5 °C осмотическое давление снижается на 7%, но вязкость воды возрастает на 52%, что требует увеличения рабочего давления для поддержания того же flux.
5. Расчёт требуемого рабочего давления
5.1. Уравнение рабочего давления
Минимальное рабочее давление, необходимое для получения пермеата с заданным flux:
P = πavg + Jw / A + ΔPloss
где:
- P — рабочее давление, бар;
- πavg — среднее осмотическое давление на стороне концентрата, бар;
- Jw — заданный flux, л/(м²·ч);
- A — коэффициент проницаемости мембраны, л/(м²·ч·бар);
- ΔPloss — потери давления в корпусе и трубопроводах, бар (обычно 0.3–1.0).
С учётом давления пермеата (обычно 0.5–2.0 бар для подачи потребителю):
Ppump = P + Ppermeate
5.2. Практическая формула для BWRO
Для систем обратного осмоса солоноватых вод (BWRO) применяется упрощённая формула:
Ppump ≈ πc,avg + (J / A) + 1.5
где 1.5 бар — суммарный запас на потери и давление пермеата.
5.3. Практическая формула для SWRO
Для систем обратного осмоса морской воды (SWRO):
Ppump ≈ πc,avg + (J / A) + 2.0
Типичные рабочие давления:
| Тип воды | TDS, мг/л | πfeed, бар | πc,avg, бар | Рабочее давление, бар |
|---|---|---|---|---|
| Пресная (водопровод) | 200 – 500 | 0.15 – 0.4 | 0.6 – 1.6 | 5 – 10 |
| Солоноватая (скважина) | 1 000 – 5 000 | 0.75 – 3.8 | 3 – 15 | 10 – 20 |
| Морская вода | 30 000 – 45 000 | 22 – 34 | 44 – 68 | 55 – 80 |
6. Конфигурация ступеней и корпусов
6.1. Архитектура мембранной системы
Мембранная система состоит из:
- Мембранных элементов — стандартные спирально-навитые элементы диаметром 4" или 8" и длиной 40" (1016 мм) или 36" (914 мм).
- Корпусов (pressure vessels) — цилиндрические сосуды высокого давления, в которые устанавливается от 1 до 8 элементов последовательно.
- Ступеней (stages) — группы корпусов, концентрат которых последовательно подаётся на следующую группу для дополнительного извлечения.
- Потоков (passes) — группы корпусов, пермеат которых последовательно подаётся на следующую группу для дополнительного солеразделения.
6.2. Правило конфигурации ступеней (15/10/5)
Для предотвращения ускоренного fouling на последних элементах корпуса (где концентрация солей максимальна) применяется правило ограничения flux:
| Тип воды | Максимальный flux на последнем элементе, LMH | Правило |
|---|---|---|
| Пресная вода (TDS < 1000 мг/л) | 18 – 22 | Не строго |
| Солоноватая вода (1000 < TDS < 10 000 мг/л) | 15 | Правило 15 LMH |
| Морская вода (TDS > 10 000 мг/л) | 10 | Правило 10 LMH |
| Высокий риск fouling | 5 | Правило 5 LMH |
6.3. Типовые конфигурации
Конфигурация записывается в формате «A:B:C», где A — количество корпусов в первой ступени, B — во второй, C — в третьей.
| Тип воды | Recovery, % | Типовая конфигурация | Обоснование |
|---|---|---|---|
| Пресная вода | 50 – 75 | 2:1 или 3:2 | Низкий риск fouling, достаточно одной ступени |
| Солоноватая вода | 60 – 80 | 3:2:1 или 4:2:1 | Ограничение flux на последних элементах |
| Морская вода | 35 – 50 | 5:3 или 6:4:2 | Высокое осмотическое давление, низкий recovery |
| Высокий риск fouling | 50 – 65 | 4:2:1 или 5:3:2 | Правило 5 LMH, увеличенная площадь мембран |
6.4. Расчёт количества элементов и корпусов
Шаг 1. Требуемая площадь мембран:
Stotal = Qp / J
где Qp — требуемая производительность по пермеату, л/ч; J — заданный flux, л/(м²·ч).
Шаг 2. Количество элементов:
Nelements = Stotal / Selement
где Selement — активная площадь одного элемента:
- 8" × 40" элемент: S ≈ 37 м² (для BW) или 400 ft²
- 4" × 40" элемент: S ≈ 8.5 м²
Шаг 3. Количество элементов в корпусе:
Обычно 6 или 7 элементов в корпусе длиной 40". При длине 36" — 6 элементов.
Шаг 4. Количество корпусов:
Nvessels = Nelements / Nper_vessel
Шаг 5. Распределение корпусов по ступеням:
Согласно выбранной конфигурации (A:B:C) и правилу flux.
7. Расчёт удельного энергопотребления (SEC)
7.1. Определение SEC
SEC (Specific Energy Consumption) — количество электроэнергии, затрачиваемое на производство 1 м³ пермеата:
SEC = Ppump / (36 × ηpump × ηmotor × Y)
где:
- SEC — удельное энергопотребление, кВт·ч/м³;
- Ppump — рабочее давление насоса, бар;
- ηpump — КПД насоса (0.65–0.85);
- ηmotor — КПД двигателя (0.90–0.95);
- Y — recovery (в долях единицы).
7.2. Теоретический минимум SEC
Минимальная энергия, необходимая для разделения, определяется термодинамикой:
SECmin = πavg / (36 × Y)
Для типичной BWRO-системы (πavg = 10 бар, Y = 0.75):
SECmin = 10 / (36 × 0.75) = 0.37 кВт·ч/м³
Реальное SEC с учётом КПД насоса и двигателя:
SECreal = 0.37 / (0.75 × 0.92) = 0.54 кВт·ч/м³
7.3. Типичные значения SEC
| Тип системы | TDS, мг/л | SEC, кВт·ч/м³ |
|---|---|---|
| BWRO (солоноватая вода) | 1 000 – 5 000 | 0.5 – 2.0 |
| SWRO (морская вода) | 30 000 – 45 000 | 2.5 – 4.5 |
| SWRO с энергорекуперацией | 30 000 – 45 000 | 2.0 – 3.0 |
| NF | 500 – 2 000 | 0.3 – 1.0 |
| UF | < 1 000 | 0.1 – 0.3 |
7.4. Энергорекуперация для SWRO
Для морских систем применяются устройства рекуперации энергии из концентрата:
- Турбодетандеры (Hydraulic Turbines): возвращают 50–60% энергии концентрата. Снижение SEC на 15–20%.
- Pressure Exchangers (PX): изобарные рекуператоры, возвращают 90–96% энергии. Снижение SEC на 30–40%. Стандарт для современных SWRO.
8. Пример расчёта системы обратного осмоса
Исходные данные:
- Производительность по пермеату: Qp = 5 м³/ч = 5000 л/ч
- Состав питания: TDS = 3000 мг/л, Fe²⁺ = 0.1 мг/л, жёсткость = 8 мг-экв/л
- Температура воды: 15 °C
- Требуемое качество пермеата: TDS < 100 мг/л
- Целевой recovery: 75%
Шаг 1. Проверка лимитирующих факторов:
- LSI в концентрате (при Y = 0.75, CF = 4): требуется расчёт по методике LSI. Принимаем LSIconcentrate ≈ +1.5 → требуется антискалант.
- Fe²⁺ = 0.1 мг/л < 0.05 мг/л (лимит для RO без предварительной очистки) — требуется предварительное обезжелезивание до < 0.05 мг/л.
- SDI₁₅ должен быть < 5 (требуется UF-предподготовка).
Шаг 2. Расчёт осмотического давления:
πfeed = 0.0075 × 3000 = 22.5 бар
Cc = 3000 × 4 = 12 000 мг/л
πc = 0.0075 × 12 000 = 90 бар
πavg = (22.5 + 90) / 2 = 56.25 бар
Шаг 3. Выбор мембраны:
Для солоноватой воды с TDS = 3000 мг/л применяются BW-мембраны. Принимаем стандартную 8" × 40" мембрану с:
- A = 5.5 л/(м²·ч·бар) (при 25 °C)
- Селективность R = 99.5%
- Selement = 37 м²
- Максимальный flux = 25 LMH
Шаг 4. Температурная коррекция:
TCF при 15 °C ≈ 0.80 (по данным производителя).
A15°C = 5.5 × 0.80 = 4.4 л/(м²·ч·бар)
Шаг 5. Задание flux:
Принимаем J = 20 LMH (с запасом от лимита 25 LMH).
Шаг 6. Требуемая площадь мембран:
Stotal = 5000 / 20 = 250 м²
Шаг 7. Количество элементов:
Nelements = 250 / 37 = 6.76 → принимаем 7 элементов
Шаг 8. Конфигурация:
7 элементов при 6 элементах в корпусе → 2 корпуса в одной ступени. Конфигурация: 2:1 (2 корпуса в первой ступени, 1 корпус во второй).
Шаг 9. Расчёт рабочего давления:
P = πavg + J / A15°C + 1.5 = 56.25 + 20 / 4.4 + 1.5 = 56.25 + 4.55 + 1.5 = 62.3 бар
Это очень высокое давление для BWRO. Причина — высокая минерализация (3000 мг/л) и низкая температура (15 °C). Принимаем рабочее давление 65 бар.
Шаг 10. Подбор насоса высокого давления:
Qfeed = Qp / Y = 5 / 0.75 = 6.67 м³/ч
Ppump = 65 бар = 650 м вод. ст.
Мощность насоса: N = (Q × P) / (3600 × η) = (6.67 × 650) / (3600 × 0.70) = 1.72 кВт
Принимаем насос мощностью 2.2 кВт с запасом.
Шаг 11. Расчёт SEC:
SEC = 65 / (36 × 0.70 × 0.92 × 0.75) = 65 / 17.39 = 3.74 кВт·ч/м³
Это высокое значение SEC обусловлено высокой минерализацией и низкой температурой. Для снижения SEC можно:
- Увеличить температуру воды (подогрев).
- Снизить recovery (увеличится расход концентрата, но снизится давление).
- Применить NF-предподготовку для снижения TDS питания.
Шаг 12. Проверка качества пермеата:
Cp = Cf × (1 - R/100) = 3000 × (1 - 0.995) = 15 мг/л
15 мг/л < 100 мг/л — качество пермеата соответствует требованию.
9. Ограничения и лимитирующие факторы
9.1. Лимитирующие параметры для RO
| Параметр | Лимит для BWRO | Лимит для SWRO | Обоснование |
|---|---|---|---|
| SDI₁₅ | < 5 | < 3 | Коллоидное загрязнение мембран |
| Мутность | < 1 NTU | < 0.5 NTU | Механическое загрязнение |
| Fe²⁺ + Fe³⁺ | < 0.05 мг/л | < 0.01 мг/л | Окисление и осаждение на мембране |
| Mn²⁺ | < 0.05 мг/л | < 0.01 мг/л | Осаждение MnO₂ |
| Остаточный Cl₂ | < 0.1 мг/л | 0 мг/л | Окисление полиамидного слоя |
| LSI в концентрате | < +1.8 (с антискалантом) | < +0.5 | Карбонатное солеотложение |
| SiO₂ в концентрате | < 100–120 мг/л | < 80 мг/л | Полимеризация силикатов |
| Индекс насыщения CaSO₄ | < 1.0 (с антискалантом) | < 0.5 | Гипсовое солеотложение |
| Температура | 5 – 45 °C | 5 – 45 °C | Ограничение материала мембраны |
| pH (при работе) | 2 – 11 | 2 – 11 | Гидролиз полиамида |
| pH (при CIP) | 1 – 13 (кратковременно) | 1 – 13 (кратковременно) | Только при промывке |
9.2. Concentration Polarization
Концентрационная поляризация — повышение концентрации солей у поверхности мембраны по сравнению с объёмом потока. Коэффициент поляризации:
β = Csurface / Cbulk
Типичные значения β = 1.1–1.3. При β > 1.5 резко возрастает риск солеотложения и fouling.
Снижение β достигается:
- Увеличением скорости тангенциального потока (cross-flow velocity).
- Применением турбулизаторов (spacers) в мембранном элементе.
- Снижением flux.
10. Типовые ошибки проектирования
- Расчёт только по исходной воде без учёта концентрата. LSI, осмотическое давление и концентрация солей рассчитываются по параметрам концентрата (с учётом CF), а не питания. Это основная ошибка, приводящая к солеотложению.
- Превышение максимального flux. Работа при flux выше лимита производителя (25 LMH для BW, 20 LMH для SW) приводит к ускоренному fouling и сокращению срока службы мембран в 2–3 раза.
- Игнорирование температурной коррекции. Подбор мембран и насосов по 25 °C без проверки при 5 °C. Зимой производительность падает на 30–40%, требуется увеличение давления или снижение flux.
- Отсутствие предварительной очистки. Подача воды с Fe²⁺ > 0.05 мг/л или SDI₁₅ > 5 на RO без предподготовки приводит к необратимому fouling за 3–6 месяцев.
- Неправильная конфигурация ступеней. Применение конфигурации 1:1 для солоноватой воды при recovery 75% приводит к превышению flux на последних элементах и ускоренному fouling.
- Занижение давления насоса. Расчёт давления по πfeed, а не по πc,avg. Реальное давление в 2–3 раза выше, что приводит к недостаточной производительности.
- Отсутствие дозирования антискаланта. При LSIconcentrate > 0 или индексе насыщения CaSO₄ > 0.5 обязательно дозирование антискаланта. Без него солеотложение гарантировано за 3–6 месяцев.
- Отсутствие дехлорирования. Подача воды с остаточным Cl₂ > 0.1 мг/л на полиамидные мембраны приводит к необратимому окислению активного слоя. Требуется дозирование NaHSO₃ или угольный фильтр.
- Игнорирование кремния. При SiO₂ > 50 мг/л в питании и recovery 75% концентрация SiO₂ в концентрате превышает 200 мг/л, что приводит к силикатному fouling, не удаляемому стандартными CIP.
- Отсутствие резервирования. Один насос высокого давления без резерва. При выходе насоса из строя — остановка всей системы. Для ответственных объектов требуется N+1 резервирование.
- Неправильный подбор CIP-системы. Отсутствие CIP-системы или её неправильная конфигурация приводит к невозможности восстановления мембран и их преждевременной замене (через 1–2 года вместо 5–7 лет).
- Завышение recovery без расчёта. Попытка получить recovery 85% для воды с высокой минерализацией без проверки LSI, индекса гипса и кремния приводит к катастрофическому fouling.
11. Нормативная и справочная база
- ГОСТ Р 57060-2016 «Мембраны обратного осмоса и нанофильтрации. Методы контроля».
- ГОСТ Р 56237-2014 «Вода питьевая. Отбор проб на станциях водоподготовки и в трубопроводных распределительных системах».
- СП 31.13330.2021 «Водоснабжение. Наружные сети и сооружения».
- ASTM D5472-18 «Standard Test Method for Water Permeability of Reverse Osmosis Membranes».
- ASTM D4189-07 «Standard Test Method for Silt Density Index (SDI) of Water».
- DuPont FilmTec™ Reverse Osmosis Membranes Technical Manual. — разделы Formulas and Calculations, System Design Guidelines.
- Hydranautics Technical Service Bulletins. — разделы Design Guidelines, Cleaning Procedures.
- Toray Membrane Technical Data. — спецификации мембран и методики расчёта.
- Crittenden, J.C., et al. MWH's Water Treatment: Principles and Design. 3rd ed. Wiley, 2012. — разделы 14 (Membrane Processes), 15 (Reverse Osmosis).
- Metcalf & Eddy, Inc. Wastewater Engineering: Treatment and Resource Recovery. 5th ed. McGraw-Hill, 2014. — раздел по мембранным процессам.
Смотри также
- Расчёт индекса насыщения Ланжелье (LSI) для мембранных систем водоподготовки
- Индекс Ларсона-Сколда (LR) и оценка коррозионной агрессивности воды
- Индекс стабильности Ризнара (RSI)
- Расчёт производительности систем водоподготовки
- Гидравлический расчёт систем водоподготовки
- Расчёт реагентного хозяйства систем водоподготовки
- Калькулятор дозирования антискаланта