WesterntechVN – Trong xử lý nước thải sinh học hiếu khí, ba hằng số động học quan trọng nhất của vi khuẩn dị dưỡng — hằng số phân hủy nội sinh kp, hằng số tốc độ tiêu thụ cơ chất riêng k, và hằng số bão hòa Ks — đều phụ thuộc nhiệt độ theo cách riêng của chúng.
Bài viết này phân tích chi tiết phương trình PT 6-23: kp,T = kp,20 × (1,04)^(T-20) và PT 6-24: kT = k20 × (1,09)^(T-20) cùng với Ks,T = Ks,20 × (1,09)^(T-20), so sánh các hệ số hiệu chỉnh θ = 1,04 và θ = 1,09, và trình bày bảng tính toán cụ thể cho hai điều kiện khí hậu đặc trưng của Việt Nam: Hà Nội mùa đông (15°C) và TP.HCM quanh năm (30°C).
1. Tại sao mỗi hằng số động học có hệ số hiệu chỉnh nhiệt độ khác nhau?
Trong mô hình động học Monod mở rộng có phân hủy nội sinh, quá trình sinh học của vi khuẩn dị dưỡng được mô tả bởi bốn hằng số chính:
- μm: tốc độ phát triển riêng cực đại (d⁻¹ hoặc h⁻¹)
- Ks: hằng số bão hòa cơ chất (mg COD/L hoặc mg BOD/L)
- Y: hệ số sinh khối (g VSS/g COD), thường được coi là không phụ thuộc nhiệt độ đáng kể
- kp: hằng số phân hủy nội sinh (d⁻¹)
Trong thực hành thiết kế, thay vì dùng μm riêng biệt, kỹ sư thường dùng hằng số k = μm/Y (hằng số tốc độ tiêu thụ cơ chất riêng cực đại, g COD/(g VSS·d)), do đó hệ thống có thể viết lại theo tốc độ tiêu thụ cơ chất:
vsu = -k × X × S / (Ks + S)
Câu hỏi đặt ra: tại sao kp lại có hệ số nhiệt độ θ = 1,04 (giống như hệ số của bùn hoạt tính tổng thể), còn k và Ks lại có θ = 1,09 (cao hơn đáng kể)?
Câu trả lời nằm ở bản chất sinh hóa của từng quá trình:
Phân hủy nội sinh (endogenous decay) là quá trình vi khuẩn tự oxy hóa chất dự trữ nội bào và thành phần tế bào khi thiếu cơ chất từ bên ngoài. Đây là quá trình tương đối chậm, ít phụ thuộc vào hệ enzyme phức tạp và điều kiện môi trường, do đó nhạy cảm với nhiệt độ ở mức độ thấp hơn (θ = 1,04).
Tiêu thụ cơ chất và ái lực cơ chất (Ks) là các quá trình liên quan trực tiếp đến hệ enzyme oxy hóa khử ngoại bào (extracellular enzymes), cơ chế vận chuyển qua màng, và chuỗi phản ứng sinh hóa cần xúc tác enzyme ở tốc độ cao. Các quá trình này nhạy cảm hơn với nhiệt độ (θ = 1,09).
Điều thú vị là Ks,T cũng có cùng hệ số θ = 1,09 với kT. Ks là thông số phản ánh ái lực của enzyme với cơ chất: Ks nhỏ = ái lực cao, vi khuẩn “bắt” được cơ chất ngay cả khi nồng độ cơ chất thấp. Khi nhiệt độ tăng, ái lực enzyme thường tăng (Ks giảm), dẫn đến Ks,T giảm khi T > 20°C và tăng khi T < 20°C.
2. Phương trình PT 6-23: kp,T = kp,20 × (1,04)^(T-20)
2.1 Phân tích phương trình
kp,T = kp,20 × (1,04)^(T-20) — (PT 6-23)
Đây là phương trình hiệu chỉnh nhiệt độ cho hằng số phân hủy nội sinh kp. Hệ số θ = 1,04 nghĩa là mỗi 1°C tăng, kp tăng 4%, và mỗi 1°C giảm, kp giảm ~3,85% (phép tính nghịch đảo).
Giá trị kp,20 tham khảo cho vi khuẩn dị dưỡng:
Trong tài liệu kỹ thuật quốc tế, kp,20 (còn ký hiệu là kd hoặc b) nằm trong khoảng:
- Điều kiện hiếu khí: kp,20 = 0,05-0,10 d⁻¹ (giá trị điển hình 0,06-0,08 d⁻¹)
- Điều kiện thiếu khí (anoxic): kp,20 = 0,03-0,04 d⁻¹ (thường lấy ~60-65% giá trị hiếu khí)
- Điều kiện kị khí: kp,20 = 0,01-0,03 d⁻¹
Giá trị này phản ánh thực tế: phân hủy nội sinh phụ thuộc vào nồng độ oxy hòa tan và điều kiện oxi hóa khử. Trong bể thiếu oxy, vi khuẩn tiêu thụ chất dự trữ nội bào chậm hơn.
2.2 Tính toán kp,T theo nhiệt độ
Lấy kp,20 = 0,07 d⁻¹ (giá trị điển hình cho bùn hoạt tính hiếu khí):
Bảng kp,T theo nhiệt độ (kp,20 = 0,07 d⁻¹):
| Nhiệt độ T (°C) | (1,04)^(T-20) | kp,T (d⁻¹) | % thay đổi so với 20°C |
|---|---|---|---|
| 35°C | 1,04^15 = 1,800 | 0,126 | +80,0% |
| 30°C | 1,04^10 = 1,480 | 0,104 | +48,0% |
| 25°C | 1,04^5 = 1,217 | 0,085 | +21,7% |
| 20°C | 1,04^0 = 1,000 | 0,070 | 0% |
| 15°C | 1,04^(-5) = 0,822 | 0,058 | -17,8% |
| 12°C | 1,04^(-8) = 0,731 | 0,051 | -26,9% |
| 10°C | 1,04^(-10) = 0,676 | 0,047 | -32,4% |
Ý nghĩa thực tế: Khi nhiệt độ xuống 15°C (Hà Nội mùa đông), kp giảm từ 0,07 xuống còn 0,058 d⁻¹. Điều này có vẻ tốt — vi khuẩn ít bị tự phân hủy hơn. Tuy nhiên, cần đặt trong bối cảnh tổng thể: đồng thời tốc độ phát triển (thể hiện qua k) cũng giảm mạnh hơn, dẫn đến nồng độ sinh khối khó duy trì ở mức thiết kế.
3. Phương trình PT 6-24: kT = k20 × (1,09)^(T-20)
3.1 Hệ số θ = 1,09 — Nhạy cảm nhiệt độ cao
kT = k20 × (1,09)^(T-20) — (PT 6-24)
So sánh θ = 1,09 với θ = 1,04:
- Mỗi 10°C thay đổi: với θ = 1,04, hệ số là 1,04^10 = 1,48; với θ = 1,09, hệ số là 1,09^10 = 2,37
- Tức là k nhạy cảm với nhiệt độ gần gấp đôi so với kp
Điều này có ý nghĩa thiết kế quan trọng: khi nhiệt độ giảm từ 20°C xuống 10°C, kp chỉ giảm 32% nhưng k giảm đến 58%. Tốc độ phát triển vi khuẩn giảm nhanh hơn nhiều so với tốc độ tự phân hủy, dẫn đến khó duy trì nồng độ bùn (MLSS) trong bể.
Giá trị k20 tham khảo:
k20 = μm,20 / Y
Với các giá trị điển hình:
- μm,20 = 5-20 d⁻¹ (vi khuẩn dị dưỡng, hiếu khí)
- Y = 0,4-0,6 g VSS/g COD
- k20 = μm/Y = thường 8-30 g COD/(g VSS·d)
Trong thực hành thiết kế, nhiều tài liệu dùng k20 = 5-8 d⁻¹ cho điều kiện wastewater sinh hoạt.
3.2 Bảng tính kT theo nhiệt độ
Lấy k20 = 6,0 d⁻¹ (giá trị tham khảo cho nước thải sinh hoạt, điều kiện hiếu khí):
| Nhiệt độ T (°C) | (1,09)^(T-20) | kT (d⁻¹) | % thay đổi |
|---|---|---|---|
| 35°C | 1,09^15 = 3,642 | 21,9 | +264% |
| 30°C | 1,09^10 = 2,367 | 14,2 | +137% |
| 25°C | 1,09^5 = 1,539 | 9,2 | +53,9% |
| 20°C | 1,09^0 = 1,000 | 6,0 | 0% |
| 15°C | 1,09^(-5) = 0,650 | 3,9 | -35,0% |
| 12°C | 1,09^(-8) = 0,502 | 3,0 | -49,8% |
| 10°C | 1,09^(-10) = 0,422 | 2,5 | -57,8% |
Nhận xét quan trọng: Tại T = 30°C (điều kiện miền Nam Việt Nam), kT = 14,2 d⁻¹ — cao hơn 2,37 lần so với 20°C. Điều này có nghĩa là bể xử lý bùn hoạt tính hoạt động ở miền Nam có thể được thiết kế với thể tích nhỏ hơn đáng kể so với thiết kế tiêu chuẩn ở 20°C.
4. Hằng số bão hòa Ks,T = Ks,20 × (1,09)^(T-20)
4.1 Ks và ý nghĩa sinh học
Ks là hằng số Michaelis-Menten, đại diện cho nồng độ cơ chất S mà ở đó tốc độ tiêu thụ đạt một nửa tốc độ cực đại. Ks nhỏ biểu thị ái lực cao của enzyme với cơ chất — vi khuẩn hoạt động hiệu quả ngay cả khi nồng độ cơ chất thấp.
Ks,T = Ks,20 × (1,09)^(T-20)
Khi T > 20°C: (1,09)^(T-20) > 1 → Ks,T > Ks,20 → ái lực enzyme giảm đi
Điều này nghe có vẻ mâu thuẫn — nhiệt độ cao lại làm giảm ái lực enzyme? Thực ra đây là hiện tượng đã được xác nhận thực nghiệm: ở nhiệt độ cao, enzyme hoạt động với tốc độ cực đại nhanh hơn (kT tăng mạnh) nhưng đồng thời cấu hình enzyme cũng linh hoạt hơn, dẫn đến mức độ “bắt” cơ chất không còn chặt chẽ như ở nhiệt độ thấp (Ks tăng).
Khi T < 20°C: Ks,T < Ks,20 → ái lực enzyme tăng. Ở nhiệt độ thấp, enzyme “bắt” chặt cơ chất hơn, nhưng vì kT thấp nên dù bắt được cơ chất thì cũng xử lý chậm.
4.2 Tác động của Ks,T đến nồng độ cơ chất ra (Se)
Trong bể bùn hoạt tính hoàn toàn hỗn hợp (CSTR) ở trạng thái ổn định, nồng độ cơ chất ra Se được xác định bởi:
Se = Ks × (kp + 1/SRT) / (kT – kp – 1/SRT) × (1 + kp × SRT) / SRT
(Phương trình được rút ra từ cân bằng sinh khối và cơ chất trong bể CSTR)
Nhìn vào phương trình này, khi nhiệt độ thay đổi:
- kT giảm (T giảm) → tử số thay đổi ít, mẫu số giảm → Se tăng (xử lý kém hơn)
- Ks,T giảm (T giảm) → có xu hướng giảm Se (bù một phần cho sự giảm kT)
- kp,T giảm (T giảm) → mẫu số tăng một chút (bù thêm)
Sự tương tác phức tạp này cho thấy khi T giảm, ảnh hưởng lớn nhất vẫn là từ kT (θ = 1,09) chứ không phải từ Ks,T hoặc kp,T. Do đó, thiết kế cần tập trung vào việc bù đắp sự suy giảm của kT khi nhiệt độ thấp.
4.3 Bảng tổng hợp Ks,T
Lấy Ks,20 = 60 mg BOD/L (giá trị điển hình cho nước thải sinh hoạt):
| T (°C) | (1,09)^(T-20) | Ks,T (mg BOD/L) | Nhận xét |
|---|---|---|---|
| 30°C | 2,367 | 142 | Ái lực rất thấp, cần S cao để đạt tốc độ tối đa |
| 25°C | 1,539 | 92 | Ái lực trung bình |
| 20°C | 1,000 | 60 | Chuẩn thiết kế |
| 15°C | 0,650 | 39 | Ái lực tăng, nhưng kT thấp hơn nhiều |
| 10°C | 0,422 | 25 | Ái lực rất cao, nhưng tốc độ cực đại rất thấp |
5. Phân tích tổng hợp cho hai vùng khí hậu Việt Nam
5.1 Kịch bản Hà Nội mùa đông — T = 15°C
Đây là điều kiện bất lợi nhất mà hệ thống xử lý nước thải đô thị phía Bắc phải đối mặt. Lấy các giá trị chuẩn ở 20°C:
- kp,20 = 0,07 d⁻¹
- k20 = 6,0 d⁻¹
- Ks,20 = 60 mg COD/L
Tại T = 15°C:
kp,15 = 0,07 × (1,04)^(15-20) = 0,07 × 0,822 = 0,058 d⁻¹
k15 = 6,0 × (1,09)^(15-20) = 6,0 × 0,650 = 3,90 d⁻¹
Ks,15 = 60 × (1,09)^(15-20) = 60 × 0,650 = 39 mg COD/L
Hậu quả đối với thiết kế SRT:
SRT tối thiểu để tránh bùn bị rửa trôi (washout):
SRT_min = 1 / (k15/Ks,15 × Se_max – kp,15)
Khi nhiệt độ giảm xuống 15°C, k15/Ks,15 = 3,90/39 = 0,100 (d·mg COD/L)⁻¹. So sánh với ở 20°C: k20/Ks,20 = 6,0/60 = 0,100 (d·mg COD/L)⁻¹.
Tỷ số k/Ks không thay đổi khi θ_k = θ_Ks! Đây là một kết quả quan trọng: vì cả k và Ks đều có cùng hệ số θ = 1,09, tỷ số k/Ks = k20/Ks,20 = const, không phụ thuộc nhiệt độ.
Ý nghĩa: nồng độ cơ chất ra Se ở trạng thái ổn định không thay đổi theo nhiệt độ nếu SRT được giữ nguyên. Tuy nhiên, tốc độ đạt trạng thái ổn định chậm hơn khi T giảm, và SRT cần thiết để duy trì hiệu quả cao hơn (vì kp không đổi tương quan như k).
Điều chỉnh SRT cần thiết tại Hà Nội mùa đông:
Nếu thiết kế ở 20°C với SRT = 10 ngày, thì tại 15°C, để đạt cùng mức độ xử lý, cần:
Sử dụng phương trình cân bằng sinh khối trong CSTR: 1/SRT = Y × k × S/(Ks + S) – kp
Để duy trì Se không đổi, cần SRT_T thỏa mãn: 1/SRT_T = Y × kT × Se/(Ks,T + Se) – kp,T
Vì k/Ks không đổi, nhưng kp,15 = 0,058 < kp,20 = 0,070, nên 1/SRT_15 = Y × k15 × Se/(Ks,15 + Se) – 0,058
Với k15 thấp hơn k20 (3,90 vs 6,0), cần Se cao hơn một chút hoặc SRT dài hơn để cân bằng. Kết quả thực tế: SRT_15 tăng khoảng 30-50% so với SRT_20.
5.2 Kịch bản TP. Hồ Chí Minh — T = 30°C
Đây là điều kiện thuận lợi nhất cho xử lý sinh học hiếu khí. Tại T = 30°C:
kp,30 = 0,07 × (1,04)^(30-20) = 0,07 × 1,480 = 0,104 d⁻¹
k30 = 6,0 × (1,09)^(30-20) = 6,0 × 2,367 = 14,2 d⁻¹
Ks,30 = 60 × (1,09)^(30-20) = 60 × 2,367 = 142 mg COD/L
Phân tích: kT tăng gấp 2,37 lần, nhưng kp cũng tăng gấp 1,48 lần. Bùn phân hủy nội sinh nhanh hơn đáng kể ở 30°C — cần tăng tỷ lệ bùn tuần hoàn hoặc điều chỉnh SRT để bù.
Tỷ số k/Ks tại 30°C: k30/Ks,30 = 14,2/142 = 0,100 (d·mg COD/L)⁻¹ — đúng bằng tỷ số tại 20°C, xác nhận tính bất biến nhiệt độ của k/Ks.
Hàm ý cho SRT tại TP.HCM:
Kp,30 = 0,104 d⁻¹ cao gấp 1,48 lần kp,20. Do đó, dù hoạt tính vi khuẩn mạnh hơn nhiều, SRT không thể giảm tùy tiện vì cần đủ thời gian để bù cho phân hủy nội sinh cao hơn.
Ước tính SRT_30 thực tế cho nước thải sinh hoạt (xử lý BOD): khoảng 5-7 ngày (ngắn hơn so với 10 ngày ở 20°C), nhưng phải kiểm tra điều kiện nitrat hóa nếu cần loại bỏ nitơ.
6. So sánh hai hệ số θ = 1,04 và θ = 1,09 trong bối cảnh thực tế
6.1 Tác động tích lũy theo thời gian
Khi nhiệt độ thay đổi từ từ theo mùa, hiệu ứng của hai hệ số tích lũy và tương tác:
Tỷ số kT/kp,T theo nhiệt độ:
- T = 30°C: kT/kp,T = (k20 × 2,367) / (kp,20 × 1,480) = 2,367/1,480 × k20/kp,20 = 1,599 × k20/kp,20
- T = 20°C: kT/kp,T = k20/kp,20 (chuẩn)
- T = 15°C: kT/kp,T = (k20 × 0,650) / (kp,20 × 0,822) = 0,791 × k20/kp,20
- T = 10°C: kT/kp,T = (k20 × 0,422) / (kp,20 × 0,676) = 0,624 × k20/kp,20
Tỷ số kT/kp,T giảm đáng kể khi nhiệt độ giảm. Khi tỷ số này giảm, hiệu quả tổng hợp của quá trình xử lý giảm và SRT cần thiết tăng lên. Đây là lý do tại sao hệ bùn hoạt tính ở miền Bắc vào mùa đông cần SRT dài hơn nhiều so với mùa hè.
6.2 Thiết kế thực hành: bảng tra cứu kỹ sư
Bảng dưới đây tổng hợp hệ số điều chỉnh cho ba hằng số tại các nhiệt độ điển hình của Việt Nam:
| Nhiệt độ | Hà Nội đông | Hà Nội xuân/thu | Chuẩn TK | HCM mùa mưa | HCM mùa khô |
|---|---|---|---|---|---|
| T (°C) | 15 | 22 | 20 | 28 | 32 |
| θkp = (1,04)^(T-20) | 0,822 | 1,082 | 1,000 | 1,369 | 1,605 |
| θk = (1,09)^(T-20) | 0,650 | 1,188 | 1,000 | 2,172 | 3,140 |
| θKs = (1,09)^(T-20) | 0,650 | 1,188 | 1,000 | 2,172 | 3,140 |
Bảng này cho phép kỹ sư nhanh chóng tra cứu hệ số điều chỉnh mà không cần tính toán lại từ đầu. Ví dụ: tại Hà Nội tháng 1 (T = 15°C), k chỉ bằng 65% giá trị thiết kế chuẩn — cần bù bằng SRT dài hơn hoặc tăng thể tích bể.
7. Ứng dụng trong tính toán thiết kế bùn hoạt tính
7.1 Phương trình tính lượng bùn sản sinh theo nhiệt độ
Lượng bùn sản sinh thực tế trong bể bùn hoạt tính:
Px = Y’ × Qv × (S0 – Se) / 1000 (kg VSS/ngày)
Trong đó Y’ = Y – kp,T × Y × SRT / (1 + kp,T × SRT) là hệ số sinh khối thực tế.
Khi T giảm:
- kp,T giảm → Y’ tăng nhẹ → bùn sản sinh tăng một chút
- Tuy nhiên, vì k/Ks không đổi, Se không đổi nhiều
Kết quả: lượng bùn sản sinh thay đổi không nhiều theo nhiệt độ trong điều kiện Việt Nam (15-32°C), chủ yếu thay đổi qua kp,T theo PT 6-23.
7.2 Tính lượng oxy cần thiết theo nhiệt độ
Nhu cầu oxy trong bể bùn hoạt tính:
O2 = Qv × (S0 – Se) × (1 – 1,42 × Y’) / 1000 (kg O2/ngày)
Khi T tăng → Y’ giảm nhẹ (vì kp,T tăng, phân hủy nội sinh nhiều hơn) → (1 – 1,42Y’) tăng → nhu cầu oxy tăng.
Đây là lý do tại sao hệ thống cung cấp khí (aeration system) cần được thiết kế cho điều kiện nhiệt độ cao nhất (mùa hè miền Nam hoặc miền Bắc), không phải điều kiện trung bình năm.
8. Lưu ý thực hành và kinh nghiệm vận hành
8.1 Giám sát nhiệt độ nước thải thực tế
Nhiệt độ nước thải sinh hoạt thường cao hơn nhiệt độ không khí 3-7°C do nhiệt từ hoạt động tắm rửa, giặt giũ. Ở Hà Nội vào tháng 1, nhiệt độ không khí có thể 12-13°C nhưng nhiệt độ nước thải thường 16-18°C. Do đó, kỹ sư không nên lấy nhiệt độ không khí làm cơ sở thiết kế mà cần đo trực tiếp nhiệt độ nước thải.
8.2 Điều chỉnh vận hành theo mùa
Kinh nghiệm vận hành các trạm xử lý nước thải ở miền Bắc Việt Nam cho thấy:
Mùa đông (tháng 11-3):
- Tăng SRT từ 8-10 ngày (hè) lên 12-15 ngày
- Duy trì MLSS cao hơn (3.000-4.000 mg/L thay vì 2.500-3.000 mg/L)
- Giảm tỷ lệ rút bùn dư 20-30%
- Theo dõi chặt chẽ BOD đầu ra — dấu hiệu cảnh báo sớm nhất
Mùa hè (tháng 5-9):
- Có thể giảm SRT về 8-10 ngày
- Tăng tỷ lệ rút bùn để kiểm soát MLSS
- Chú ý nguy cơ sưng bùn (sludge bulking) do filamentous bacteria phát triển mạnh ở nhiệt độ cao
8.3 Lỗi thiết kế phổ biến liên quan đến nhiệt độ
Lỗi 1: Dùng số liệu thiết kế ở 20°C cho công trình miền Bắc Thiếu hệ số an toàn nhiệt độ dẫn đến bể quá nhỏ, không đạt chuẩn xả thải vào mùa đông.
Lỗi 2: Không phân biệt θ = 1,04 và θ = 1,09 Áp dụng cùng một hệ số cho tất cả các hằng số dẫn đến tính toán sai lệch, đặc biệt trong hệ thống có kết hợp nitrat hóa.
Lỗi 3: Không tính hiệu ứng Ks,T tăng ở nhiệt độ cao Ở T = 32°C, Ks,T = 142 mg COD/L — rất cao. Nếu nồng độ COD đầu vào sau hệ sơ cấp thấp (< 200 mg/L), phần lớn nồng độ này nằm trong vùng tuyến tính của Monod (S < Ks), dẫn đến hiệu quả xử lý thực tế thấp hơn lý thuyết.
9. Kết luận
PT 6-23 (θ = 1,04 cho kp) và PT 6-24 (θ = 1,09 cho k và Ks) phản ánh hai mức độ nhạy cảm nhiệt độ khác nhau của các quá trình sinh hóa cơ bản trong bùn hoạt tính. Đặc biệt quan trọng là phát hiện rằng k và Ks có cùng hệ số θ = 1,09, khiến tỷ số k/Ks bất biến theo nhiệt độ — từ đó nồng độ cơ chất đầu ra Se không thay đổi đáng kể khi SRT được duy trì, mặc dù tốc độ đạt cân bằng thay đổi theo nhiệt độ.
Đối với thiết kế công trình xử lý nước thải tại Việt Nam, cần áp dụng đồng thời cả ba phương trình hiệu chỉnh nhiệt độ (kp,T, kT, Ks,T) và xét đến điều kiện nhiệt độ cực đoan của từng vùng — đặc biệt là Hà Nội mùa đông (15°C) và TP.HCM mùa khô (32°C). Bài toán nhiệt độ không đơn thuần là tra cứu bảng mà đòi hỏi hiểu rõ cơ chế tương tác giữa các hằng số động học để ra quyết định thiết kế đúng đắn.
