Модель ASM3: як математика описує біологічну очистку стічних вод
Оглядова стаття про модель активного мулу №3 (IWA): концепція внутрішньоклітинного накопичення, 13 змінних і 12 процесів, матриця Гуєра, кінетика, двоступенева нітрифікація, симультанна нітри-денітрифікація й калібрування — простою інженерною мовою.
Автор — Кобилко Ігор Васильович, інженер-технолог · Опубліковано · Оновлено
Навіщо взагалі моделювати очисні споруди
Біологічна очистка стічних вод — це керована екосистема: у резервуарі з активним мулом мільярди мікроорганізмів «поїдають» органіку й перетворюють азот. Проблема в тому, що зазирнути всередину цього процесу напряму не можна — видно лише вхід і вихід. Тому ще з кінця 1980-х інженери описують те, що відбувається в аеротенку, математичними моделями: вони дозволяють «прокрутити» роботу споруди на комп’ютері наперед — до того, як залито перший кубометр бетону чи змінено режим аерації на діючому об’єкті.
Модель відповідає на дуже практичні запитання: якого об’єму потрібен реактор, скільки подавати повітря, який тримати вік мулу, наскільки чистою вийде вода за піковим навантаженням і де можна зекономити енергію. ASM3 (Activated Sludge Model No. 3) — одна з базових моделей для задач, де головні гравці — вуглець і азот.
Це оглядова інженерна стаття про саму модель. Усі числові параметри наведено як типові значення за замовчуванням IWA — у реальних розрахунках їх калібрують під конкретний стік і споруду.
Сімейство моделей активного мулу IWA
ASM3 — не окрема розробка, а частина сімейства уніфікованих моделей, які Міжнародна асоціація води (IWA, раніше IAWQ) створювала як спільну «мову» для інженерів усього світу.
- ASM1 (Henze та ін., 1987) — перша уніфікована модель: окиснення органіки, нітрифікація й денітрифікація на основі концепції «загибелі–регенерації».
- ASM2 / ASM2d (1995/1999) — додано біологічне й хімічне видалення фосфору (фосфор-акумулюючі організми, PAO).
- ASM3 (Gujer, Henze, Mino, van Loosdrecht, 1999/2000) — переглянута базова модель вуглецю й азоту, побудована на концепції внутрішньоклітинного накопичення субстрату.
ASM3 не замінює ASM2d (там ключовий фосфор), а є сучаснішою альтернативою ASM1 для задач вуглецю й азоту. Канонічний опис усіх чотирьох моделей — у звіті IWA Scientific and Technical Report No. 9 (Henze та ін., 2000).
Чому знадобилась ASM3: слабкі місця ASM1
ASM1 добре послужила, але за роки практики виявились структурні обмеження, які й спонукали переписати модель:
- Переплетені процеси. Через механізм «загибелі–регенерації» продукти загибелі бактерій знову йдуть у їжу — окиснення й розпад взаємопов’язані, тож роздільно відкалібрувати їхні швидкості важко.
- Складно розділити дихання. У виміряній швидкості споживання кисню (OUR) непросто відокремити внесок росту від внеску ендогенного дихання.
- Ігнорування накопичення. Експериментально доведено: за наявності легкого субстрату мікроорганізми спершу запасають його у вигляді внутрішніх полімерів, а вже потім використовують для росту. ASM1 цього не описує.
- Погана придатність до циклічних режимів — SBR, чергування «надлишок–нестача» субстрату (feast–famine).
ASM3 переформульовує метаболізм так, щоб зробити процеси послідовними, а не циклічними, і явно ввести внутрішнє накопичення.
Головна ідея ASM3: клітина спершу запасає, потім росте
Центральна ідея звучить просто: увесь легкозасвоюваний субстрат S_S спочатку накопичується у вигляді внутрішньоклітинного продукту X_STO, і лише потім цей запас витрачається на ріст біомаси. Замість замкненого циклу ASM1 виходить лінійний ланцюг:
Ендогенне дихання приходить на зміну «загибелі–регенерації»: біомаса й накопичений запас окиснюються з витратою кисню (аеробно) або нітрату/нітриту (аноксидно), даючи інертний твердий залишок. Завдяки цьому швидкості накопичення, росту й дихання калібруються незалежно, легше інтерпретувати респірометрію, а модель точніше відтворює динаміку в реакторах із чергуванням надлишку й нестачі субстрату. Автотрофи-нітрифікатори в базовій ASM3 накопичення не мають — вони ростуть просто на амонійному азоті.
13 змінних стану: з чого «зроблена» модель
Модель стежить за 13 величинами. Розчинені компоненти позначають літерою S, тверді (завислі) — X. Одиниці переважно у ХСК (гCOD/м³), для азоту — гN/м³, для лужності — моль HCO₃⁻/м³.
| Позначення | Що це | Тип |
|---|---|---|
| S_O2 | Розчинений кисень (як від’ємний ХСК) | розчинений |
| S_I | Інертна розчинена органіка | розчинений |
| S_S | Легкозасвоюваний субстрат | розчинений |
| S_NH4 | Амонійний азот (NH₄⁺ + NH₃) | розчинений |
| S_N2 | Газоподібний азот (продукт денітрифікації) | розчинений |
| S_NOX | Сума нітрату й нітриту (один пул) | розчинений |
| S_ALK | Лужність (бікарбонатна) — контроль pH і заряду | розчинений |
| X_I | Інертна завись органіки | твердий |
| X_S | Повільно-біорозкладний субстрат (гідролізується) | твердий |
| X_H | Гетеротрофна біомаса | твердий |
| X_STO | Внутрішньоклітинний накопичений запас | твердий |
| X_A | Автотрофна (нітрифікуюча) біомаса | твердий |
| X_TS | Загальна завись твердих речовин (перерахунок у мг/л) | твердий |
12 процесів: що відбувається і що керує швидкістю
Усе, що відбувається в мулі, ASM3 зводить до 12 процесів. Їхню стехіометрію й кінетику подають у вигляді матриці Петерсена (Гуєра): рядки — процеси, стовпці — компоненти, на перетині — стехіометричні коефіцієнти, а праворуч — вираз швидкості. Такий табличний формалізм гарантує збереження маси за ХСК, азотом і зарядом: нічого не «виникає» і не «зникає». Швидкість кожного процесу — це максимальна константа, помножена на набір перемикальних (Моно) множників виду S/(K+S), які «вмикають» процес лише за наявності потрібних реагентів. Нижче — 8 ключових процесів; швидкості записано спрощено (позначення множників — у ключі під статтею).
| Процес | Що відбувається | Від чого залежить швидкість (Моно-множники) |
|---|---|---|
| Гідроліз | X_S → S_S: повільний субстрат стає доступним | k_H · [X_S/X_H ⁄ (K_X + X_S/X_H)] · X_H |
| Аеробне накопичення X_STO | S_S → X_STO, споживання O₂ | k_STO · [O₂/(K_O2+O₂)] · [S_S/(K_S+S_S)] · X_H |
| Аноксидне накопичення X_STO | S_S → X_STO, споживання NOX (денітрифікація) | k_STO · η_NOX · [K_O2/(K_O2+O₂)] · [NOX/(K_NOX+NOX)] · [S_S/(K_S+S_S)] · X_H |
| Аеробний ріст X_H | X_STO → X_H, споживання O₂ | μ_H · [O₂/(K_O2+O₂)] · [NH4/(K_NH4+NH4)] · [X_STO/X_H ⁄ (K_STO+…)] · X_H |
| Аноксидний ріст X_H | X_STO → X_H, денітрифікація | μ_H · η_NOX · [K_O2/(K_O2+O₂)] · [NOX/(K_NOX+NOX)] · [NH4…] · [X_STO…] · X_H |
| Аеробне ендогенне дихання X_H | окиснення біомаси киснем → X_I | b_H,O2 · [O₂/(K_O2+O₂)] · X_H |
| Ріст автотрофів = нітрифікація | S_NH4 → S_NOX, ріст X_A, споживання O₂ і лужності | μ_A · [O₂/(K_A,O2+O₂)] · [NH4/(K_A,NH4+NH4)] · [ALK…] · X_A |
| Ендогенне дихання X_A | окиснення автотрофів киснем → X_I | b_A,O2 · [O₂/(K_A,O2+O₂)] · X_A |
| …і ще: аноксидні ендогенні дихання X_H, X_STO, X_A | симетричні аеробним — окиснення нітратом замість кисню | b·η_NOX · [K_O2/(K_O2+O₂)] · [NOX/(K_NOX+NOX)] · X |
Стехіометрична матриця Петерсена як теплокарта
Класично матрицю Гуєра друкують числами (стехіометричні коефіцієнти), і читати її важко. Тут ми показуємо її наочно: рядки — процеси, стовпці — компоненти, а колір клітинки одразу каже, що процес із цим компонентом робить — виробляє (зелена, «+») чи споживає (червона, «−»). Порожня клітинка — компонент у процесі не бере участі. На телефоні таблиця гортається вбік.
| Процес ↓ / Компонент → | SO2 | SS | SNH4 | SNOX | SN2 | SALK | XS | XSTO | XH | XA | XI |
|---|---|---|---|---|---|---|---|---|---|---|---|
| 1. Гідроліз | + | + | + | − | |||||||
| 2. Аеробне накопичення X_STO | − | − | + | ||||||||
| 3. Аноксидне накопичення X_STO | − | − | + | + | + | ||||||
| 4. Аеробний ріст X_H | − | − | − | − | + | ||||||
| 5. Аноксидний ріст X_H | − | − | + | + | − | + | |||||
| 6. Ендогенне дихання X_H | − | + | + | − | + | ||||||
| 7. Нітрифікація (ріст X_A) | − | − | + | − | + | ||||||
| 8. Ендогенне дихання X_A | − | + | + | − | + |
Читається як «хто чим живиться»: нітрифікація (рядок 7) жадібно споживає кисень, амоній і лужність, зате виробляє нітрат і нову автотрофну біомасу; ендогенне дихання (рядки 6, 8) навпаки повертає амоній у воду й нарощує інертний залишок X_I. Саме ця «бухгалтерія» плюс збереження маси й робить ASM3 придатною для інженерного розрахунку.
Ключ позначень
У всіх таблицях і схемах діють єдині позначення:
Кінетика: чим задають швидкості
Швидкості будуються на функціях Моно та перемикальних (switching) функціях виду S/(K+S): вони «вмикають» процес за наявності реагента й «вимикають» за його відсутності. Наприклад, аеробний ріст гетеротрофів множиться на наявність кисню, запасу X_STO, амонію (як джерела азоту для біомаси) та лужності. Для аноксидних процесів кисневий множник замінюють на інгібувальний K_O2/(K_O2+S_O2), додають наявність нітрату й емпіричний коефіцієнт редукції η_NOX (< 1) — анаеробний метаболізм повільніший.
| Параметр | Типове значення (20 °C) | Що описує | Калібрується? |
|---|---|---|---|
| k_H | ≈ 3 гCOD·(гCOD)⁻¹·д⁻¹ | константа гідролізу X_S → S_S | так, чутливий |
| k_STO | ≈ 5 гCOD·(гCOD)⁻¹·д⁻¹ | швидкість накопичення S_S у X_STO | так |
| μ_H | ≈ 2 д⁻¹ | макс. швидкість росту гетеротрофів | так, чутливий |
| μ_A | ≈ 1 д⁻¹ | макс. швидкість росту автотрофів (нітрифікація) | так, критичний |
| Y_STO,O2 | ≈ 0,85 гCOD/гCOD | вихід запасу X_STO із субстрату (аеробно) | рідше |
| Y_H,O2 | ≈ 0,63 гCOD/гCOD | вихід гетеротрофної біомаси (аеробно) | рідше |
| Y_A | ≈ 0,24 гCOD/гN | вихід автотрофної біомаси на 1 г окисненого N | рідше |
| K_O2 | ≈ 0,2 гO₂/м³ | напівнасичення за киснем (аеробні процеси) | інколи |
| K_A,NH4 | ≈ 1 гN/м³ | напівнасичення нітрифікації за амонієм | так |
| b_H,O2 | ≈ 0,2 д⁻¹ | аеробне ендогенне дихання X_H | так |
Наголосимо ще раз: ці числа — орієнтир. Для конкретних стоків (особливо промислових — молочних, харчових) їх калібрують експериментально. Нітрифікація особливо чутлива до температури — це критично для холодних періодів.
Нітрифікація: один етап чи два
У базовій ASM3 нітрифікація спрощена до одного етапу: амоній одразу перетворюється на об’єднаний пул S_NOX однією групою автотрофів X_A. Це прийнятно, коли проміжний нітрит не накопичується. Фізично ж нітрифікація йде у два кроки:
- NH₄⁺ → NO₂⁻ — амоній-окиснювальні бактерії (AOB, напр. Nitrosomonas).
- NO₂⁻ → NO₃⁻ — нітрит-окиснювальні бактерії (NOB, напр. Nitrobacter/Nitrospira).
Коли важливе саме накопичення нітриту (часткова нітрифікація, nitritation/anammox, низький DO, короткий вік мулу, підвищена температура), базову модель розширюють до двоступеневої (ASM3-2N): вводять окремі групи AOB і NOB зі своїми параметрами. Це дозволяє відтворити «нітритний шунт» і його енергетичні переваги — менші витрати кисню й вуглецю. Такі розширення докладно опрацьовані у Koch та ін. (2000).
Симультанна нітри-денітрифікація й економія повітря
Окремий клас задач — симультанна нітри-денітрифікація (СНД), коли нітрифікація й денітрифікація йдуть одночасно в одному реакторі за низького розчиненого кисню (як правило, < 1,5 мг/л). Механізм — градієнт кисню усередині флокули мулу або біоплівки:
- Зовнішня (аеробна) оболонка флокули — нітрифікація, яка йде навіть за помірного об’ємного DO.
- Внутрішнє (аноксидне) ядро — денітрифікація гетеротрофами за рахунок накопиченого запасу X_STO.
Практична вигода СНД відчутна: не потрібні окрема аноксидна зона й внутрішня рециркуляція, а витрати на аерацію падають. Класичний експериментальний опис — Pochana & Keller (1999). У моделях цей ефект відтворюють або через просторовий підмодель дифузії всередині флокули (підхід Wanner & Gujer, середовище AQUASIM), або емпіричною корекцією швидкостей денітрифікації як функції DO. «Точкова» ASM3 без такого розширення градієнта сама не відтворює.
Розширення ASM3
За два десятиліття навколо ASM3 виросла ціла родина розширень під конкретні технології:
- ASM3 + BioP (EAWAG bio-P module, Rieger та ін., 2001) — додає фосфор-акумулюючі організми й біологічне видалення фосфору, зберігаючи storage-концепцію.
- Двоступенева нітрифікація (ASM3-2N) — окремі AOB і NOB для режимів із нітритом.
- Поєднання з моделями біоплівки — ASM3 як кінетичне ядро в 1-D моделі біоплівки (Wanner–Gujer); застосовується для MBBR, IFAS, біофільтрів, реакторів із плаваючими носіями, де важливі дифузійні градієнти.
- Модифікації для промислових стоків з високим ХСК/N (харчові, молочні) — корекція виходів, гідролізу, врахування легко- й важкозасвоюваних фракцій.
- Температурні й pH-розширення — врахування вільного аміаку й азотистої кислоти як інгібіторів нітрифікації.
Приклад сучасної української модифікації — модель ASM3P-2SND (В. С. Рожков, НУВГП): розширена ASM3 на 19 процесів і 16 компонентів, де азот обліковують окремо за нітритом (NO₂) і нітратом (NO₃), а швидкості коригують на дифузійні обмеження всередині флокули для коректного опису симультанної нітри-денітрифікації. Це показує, у якому напрямі базову ASM3 доопрацьовують під конкретні режими; наведені тут числа загальної ASM3 до цієї модифікації не належать.
Головний виклик практики: калібрування
ASM3 має 20+ кінетичних і стехіометричних параметрів, які неможливо надійно визначити з одного набору вимірювань. Це проблема ідентифікованості: різні набори параметрів можуть однаково добре описувати ті самі дані. Тому калібрування завжди супроводжують аналізом чутливості. Стандартний робочий процес:
- Характеризація стоку — фракціонування ХСК (S_S, S_I, X_S, X_I) та респірометрія (вимір OUR).
- Аналіз чутливості — спершу калібрують параметри, до яких вихід найчутливіший; решту лишають за замовчуванням.
- Поетапне калібрування (протоколи BIOMATH, STOWA, WERF): спочатку стаціонарні показники (вік мулу, продукція мулу), потім динамічні (добові профілі NH₄, NO₃, OUR).
- Бенчмарки BSM1/BSM2 — стандартні платформи для порівняння стратегій керування.
Система рівнянь ASM3 «жорстка» (швидкі й повільні процеси співіснують), тож її розв’язують неявними чисельними солверами (BDF) у середовищах GPS-X, WEST, SIMBA, BioWin, AQUASIM або відкритих реалізаціях на Python/MATLAB. SBR моделюють як послідовність фаз: наповнення → реакція → відстій → декантація.
Де ASM-моделювання в практиці МВК
МВК проєктує й обслуговує локальні очисні споруди для виробництв — насамперед харчових і молокопереробних, де стік «важкий»: висока органіка, жири, стрибки pH і навантаження. Саме для таких об’єктів моделі сімейства ASM корисні як інженерний інструмент розрахунку, а не абстрактна теорія.
Що дає модель на етапі проєктування й оптимізації
- Підбір об’єму й режиму реактора. Для аеротенка чи SBR модель допомагає оцінити потрібний об’єм, вік мулу й тривалість фаз під фактичне навантаження стоку, а не за усередненою нормою.
- Економія аерації. Аерація — найбільша стаття енерговитрат очисних. Розрахунок режиму низького DO зі СНД показує, де можна зменшити подачу повітря без втрати якості очищення за азотом.
- Прогноз за пікових навантажень. Молочне виробництво дає різкі коливання ХСК — модель допомагає перевірити, чи витримає біологія залпові скиди й закисання.
Важлива чесна засторога: модель — не «чарівна кнопка». Її результат вартий рівно стільки, скільки якісним був аналіз стоку й калібрування. Тому на практиці ASM-розрахунок працює в парі з лабораторним контролем і телеметрією: діючу лінію зручно вести під моніторингом зі смартфона, звіряючи фактичні DO, витрату й якість із тим, що передбачала модель.
Що ще почитати за темою
Біологічна очистка на практиці
- Очищення стічних вод молокозаводуФлотація, аеротенк-усереднювач, двоступенева біологія й ППС-доочистка — реальна технологічна лінія під «важкий» молочний стік.
- Очисні споруди під ключПроєктування, монтаж і сервіс локальних очисних для виробництв.
- Пінополістирольні фільтриПлаваюча засипка школи НУВГП — доочистка біологічно очищеної води.
Рішення МВК
- Промислова водопідготовкаЗагальний огляд технологій і бар’єрний принцип побудови схеми.
- Розумна ХВО: моніторинг зі смартфонаОнлайн-контроль параметрів лінії — DO, витрата, стан вузлів і аварійні сигнали.
- Локальні очисні спорудиНапрям МВК зі стічних вод — від жировловлювачів до біологічного блоку.
Інженерія й інструменти
- Як ми проєктуємо станцію3D-компонування в реальному приміщенні ще до монтажу.
- Онлайн-калькулятори водопідготовкиЖорсткість, SDI-тест та інші інженерні розрахунки.
- Як працює градирняЩе один розбір «об’єкт як механізм» — оборотне охолодження й реагентні програми.
- Gujer W., Henze M., Mino T., van Loosdrecht M. Activated Sludge Model No. 3 // Water Science and Technology. — 1999. — Vol. 39(1). — P. 183–193. DOI: 10.2166/wst.1999.0039.
- Henze M., Gujer W., Mino T., van Loosdrecht M. Activated Sludge Models ASM1, ASM2, ASM2d and ASM3. IWA Scientific and Technical Report No. 9. — London: IWA Publishing, 2000 (канонічний опис моделі, матриця Гуєра, типові параметри за замовчуванням). DOI: 10.2166/9781780402369.
- Koch G., Kühni M., Gujer W., Siegrist H. Calibration and validation of ASM3 for biological nitrogen removal // Water Research. — 2000. — Vol. 34(14). — P. 3580–3590 (калібрування, двоступенева нітрифікація).
- Rieger L., Koch G., Kühni M., Gujer W., Siegrist H. The EAWAG Bio-P module for ASM3 // Water Research. — 2001. — Vol. 35(16). — P. 3887–3903 (розширення для фосфору).
- Pochana K., Keller J. Study of factors affecting simultaneous nitrification and denitrification (SND) // Water Science and Technology. — 1999. — Vol. 39(6). — P. 61–68.
- Wanner O., Gujer W. A multispecies biofilm model // Biotechnology and Bioengineering. — 1986. — Vol. 28. — P. 314–328 (модель біоплівки, просторові градієнти).
- Рожков В. С. Дисертація (НУВГП): модель ASM3P-2SND — розширення ASM3 на 19 процесів і 16 компонентів з окремим обліком нітриту/нітрату та дифузійними обмеженнями (приклад української модифікації; числа загальної ASM3 звідси не запозичено).
- Значення параметрів наведено як типові (default) за IWA STR No. 9; у прикладних розрахунках підлягають калібруванню під конкретний стік і споруду.