/* Chat */

Phương Trình vsu = −(μm·X·S)/(Y·(Ks+S)) — Phân Tích Bốn Yếu Tố Quyết Định Tốc Độ Loại Bỏ Cơ Chất Trong Bể Sinh Học

16/09/2026 9 lượt xem quantri

WesterntechVN – Tốc Độ Xử Lý Phụ Thuộc Vào Điều Gì?

Một câu hỏi thực tiễn thường gặp trong vận hành nhà máy xử lý nước thải: Tại sao cùng một bể sinh học aerotank, lúc thì xử lý tốt, lúc thì COD đầu ra vượt chuẩn? Câu trả lời không chỉ đến từ một yếu tố đơn lẻ mà từ sự tương tác phức tạp của ít nhất bốn biến số khác nhau.

Phương trình 6-6 trong giáo trình GS. Lê Văn Cát (2007) biểu diễn tất cả bốn yếu tố đó trong một phương trình duy nhất:

vsu = −(μm × X × S) / (Y × (Ks + S)) (Phương trình 6-6)

Đây là phương trình “bản đồ” cho kỹ sư vận hành: mỗi biến số tương ứng với một đòn bẩy kiểm soát — tăng X, tăng S, tăng μm (qua nhiệt độ/pH), hoặc giảm Ks (qua chọn giống vi sinh thích nghi). Bài viết này phân tích chi tiết từng yếu tố và cách tối ưu hóa trong thực tiễn Việt Nam.

Bể Sinh Học

Phần 1: Nguồn Gốc Phương Trình 6-6 — Kết Hợp Monod Và Quan Hệ Sinh Khối-Cơ Chất

1.1 Dẫn Xuất Từ Các Phương Trình Cơ Bản

Phương trình 6-6 là kết quả tự nhiên từ hai phương trình đã biết:

Từ Monod (6-3): vg = μm × X × S / (Ks + S)

Từ quan hệ sinh khối-cơ chất (6-5): vg = −Y × vsu → vsu = −vg/Y

Thay thế: vsu = −(1/Y) × μm × X × S / (Ks + S) = −μm × X × S / (Y × (Ks + S))

Phương trình 6-6 cho thấy tốc độ loại bỏ cơ chất phụ thuộc vào:

  • X: mật độ sinh khối vi sinh
  • S: nồng độ cơ chất (nồng độ ô nhiễm)
  • μm: tốc độ sinh trưởng tối đa (hoạt tính vi sinh)
  • Y: hiệu suất sinh khối (bao nhiêu sinh khối tạo ra từ mỗi đơn vị cơ chất)
  • Ks: hằng số bán bão hòa (ái lực của vi sinh với cơ chất)

1.2 Tái Viết Dưới Dạng Tốc Độ Tiêu Thụ Riêng

Phương trình 6-6 có thể biểu diễn dưới dạng khác, chia cả hai vế cho X:

q ≡ |vsu|/X = (μm/Y) × S/(Ks + S) = k × S/(Ks + S)

Trong đó q là tốc độ tiêu thụ cơ chất riêng (specific substrate utilization rate, g COD/g VSS·ngày). Dạng này tách bạch ảnh hưởng của X ra khỏi ảnh hưởng của S — hữu ích để so sánh hoạt tính vi sinh giữa các hệ thống khác nhau.

Phần 2: Yếu Tố X — Nồng Độ Sinh Khối Và Cách Tối Ưu

2.1 X Tỉ Lệ Thuận Với vsu — Cơ Sở Của Công Nghệ Bùn Hoạt Tính

Từ phương trình 6-6: |vsu| ∝ X

Tức là: tăng gấp đôi MLVSS → tốc độ xử lý COD tăng gấp đôi (khi S và các điều kiện khác không đổi).

Đây chính là nguyên lý cơ bản của công nghệ bùn hoạt tính: duy trì mật độ vi sinh cao trong bể → tốc độ xử lý cao → thể tích bể nhỏ hơn.

So sánh:

  • Ao hồ ổn định (stabilization pond): X ≈ 50–200 mg VSS/L; HRT = 20–60 ngày
  • Aerotank thông thường: X ≈ 2.000–4.000 mg VSS/L; HRT = 4–8 giờ
  • Aerotank tải cao (high-rate): X ≈ 4.000–8.000 mg VSS/L; HRT = 1–3 giờ
  • MBR (Membrane Bioreactor): X ≈ 8.000–15.000 mg VSS/L; HRT = 0,5–2 giờ

Tăng X từ 3.000 lên 12.000 mg VSS/L (tương đương chuyển từ aerotank sang MBR) → tốc độ xử lý tăng 4 lần → thể tích bể giảm 4 lần. Đây là lý do MBR tốn ít đất hơn aerotank truyền thống.

2.2 Giới Hạn Tăng X — Vấn Đề Lắng Và Chi Phí Năng Lượng

Tăng X không phải vô hạn. Giới hạn thực tiễn:

Công Nghệ MLVSS Tối Đa Giới Hạn Chính
Aerotank (lắng trọng lực) 4.000–5.000 mg/L Khả năng lắng của bể lắng thứ cấp
SBR 4.000–6.000 mg/L Hiệu suất lắng trong pha Settle
MBBR (màng sinh học) 10.000–15.000 mg/L Bề mặt giá thể (carrier)
MBR 12.000–20.000 mg/L Tắc màng (membrane fouling)
Aerotank MBBR kết hợp 6.000–10.000 mg/L Cân bằng giữa lắng và màng

Thực tiễn Việt Nam: Tại các khu công nghiệp miền Nam, nhiều nhà máy đang nâng cấp từ aerotank (MLSS 3.000 mg/L) lên MBBR (MLSS 6.000 mg/L) để tăng công suất mà không cần xây thêm bể, tiết kiệm 40–60% chi phí đầu tư so với xây bể mới.

2.3 Kiểm Soát X Trong Vận Hành — Cân Bằng Sinh Khối

Trong bể liên tục ổn định, X được kiểm soát qua tốc độ xả bùn dư:

X = Y_net × SRT × (Q/V) × ΔCOD / (1 + b_a × SRT)

Để duy trì X mục tiêu, kỹ sư điều chỉnh lưu lượng bơm bùn dư Qw hàng ngày, dựa trên đo MLVSS hàng ngày.

Phần 3: Yếu Tố S — Nồng Độ Cơ Chất Và Đặc Tính Phi Tuyến

3.1 Quan Hệ Phi Tuyến Giữa S Và |vsu|

Khác với X (tuyến tính), ảnh hưởng của S lên |vsu| là phi tuyến (nonlinear) theo hàm Monod:

|vsu|/X = k × S/(Ks + S)

Phân tích theo vùng:

  • S >> Ks (S >> 60 mg BOD/L): S/(Ks + S) ≈ 1 → |vsu|/X ≈ k → không phụ thuộc S nữa
  • S ≈ Ks: S/(Ks + S) ≈ 0,5 → |vsu|/X ≈ k/2
  • S << Ks (S << 60 mg BOD/L): S/(Ks + S) ≈ S/Ks → |vsu|/X ≈ k × S/Ks (tỉ lệ tuyến tính với S)

Hệ quả quan trọng cho thiết kế bể:

  • Đầu bể (S cao): tốc độ tiêu thụ cơ chất rất nhanh, gần mức tối đa k
  • Cuối bể (S thấp): tốc độ giảm mạnh, cần thể tích lớn để đạt tiêu chuẩn đầu ra nghiêm ngặt

Đây là lý do tại sao bể plug-flow (PFR) thường hiệu quả hơn bể khuấy hoàn toàn (CSTR) trong xử lý nước thải:

  • CSTR: S trong bể = S_đầu_ra (thấp) → tốc độ thấp toàn bể → cần thể tích lớn hơn
  • PFR: S giảm dần từ S_đầu_vào đến S_đầu_ra → tốc độ cao ở đầu, thấp ở cuối → thể tích nhỏ hơn cho cùng hiệu suất

3.2 Tác Động Của Tải Trọng Sốc (Shock Loading)

Khi COD đầu vào đột ngột tăng (S tăng):

  • S vẫn trong vùng S >> Ks → |vsu| tăng tỉ lệ thuận với X (không phụ thuộc S nhiều)
  • Vấn đề không phải thiếu tốc độ mà thiếu X (MLVSS không đủ)!

Tải trọng sốc hiệu quả được xác định bởi F/M = Q×S/(V×X) = tải trọng thực / khả năng tiêu thụ. Khi F/M > 0,5 kg BOD/kg MLVSS·ngày → nguy cơ bùn phân tán (dispersed growth), COD đầu ra tăng.

Giải pháp thực tế: bể điều lưu (equalization tank) trước aerotank để san bằng tải trọng, đặc biệt tại các nhà máy công nghiệp hoạt động theo ca (shift operation).

Phần 4: Yếu Tố μm — Hoạt Tính Vi Sinh Và Các Điều Kiện Ảnh Hưởng

4.1 Nhiệt Độ — Ảnh Hưởng Theo Arrhenius

μm thay đổi mạnh theo nhiệt độ theo phương trình Arrhenius đơn giản hóa: μm(T) = μm(20°C) × θ^(T-20)

với θ = 1,04–1,08 cho vi khuẩn dị dưỡng và θ = 1,06–1,12 cho Nitrosomonas.

Vì vsu ∝ μm, nhiệt độ tác động trực tiếp đến tốc độ xử lý:

Tại nhiệt độ 30°C (điển hình miền Nam Việt Nam): μm(30°C)/μm(20°C) = 1,07^10 ≈ 1,97 — xử lý nhanh gần gấp đôi so với điều kiện chuẩn 20°C!

Tại nhiệt độ 15°C (điển hình miền Bắc mùa đông): μm(15°C)/μm(20°C) = 1,07^(-5) ≈ 0,71 — chậm hơn 29% so với 20°C.

Hàm ý thiết kế: Nhà máy tại Hà Nội cần thể tích bể lớn hơn hoặc SRT dài hơn nhà máy cùng công suất tại TP. Hồ Chí Minh, vì nhiệt độ thấp hơn làm μm giảm → cần bù đắp bằng thể tích hoặc thời gian.

4.2 pH — Ảnh Hưởng Lên Enzyme

pH ảnh hưởng đến hoạt tính enzyme vi sinh thông qua việc thay đổi trạng thái ion hóa của nhóm hoạt động:

Nhóm Vi Sinh pH Tối Ưu pH Giảm 50% Hoạt Tính
Vi khuẩn dị dưỡng 6,5–8,0 < 5,5 hoặc > 9,5
Nitrosomonas 7,5–8,5 < 6,0 hoặc > 9,0
Nitrobacter 7,0–8,0 < 5,5 hoặc > 9,5
Methanogens 6,8–7,4 < 6,0 hoặc > 8,0
Anammox 6,8–8,0 < 6,0 hoặc > 9,0

Thực tiễn tại Việt Nam: Nước thải công nghiệp chế biến thực phẩm (đường mía, cà phê, tinh bột) thường có pH thấp (4–5) → cần trung hòa trước khi vào bể sinh học. Không điều chỉnh pH → μm giảm nghiêm trọng → |vsu| giảm → COD đầu ra vượt chuẩn.

4.3 Chất Ức Chế — Tác Nhân Làm Giảm μm

Nhiều chất trong nước thải công nghiệp ức chế vi sinh vật, hiệu lực qua việc giảm μm:

  • Kim loại nặng (Zn, Cu, Cr, Cd): ức chế enzyme tại nồng độ mg/L; Zn ức chế 50% tại 0,5 mg/L
  • Ammonia tự do (NH₃): ức chế Nitrosomonas tại > 150 mg NH₃-N/L; ức chế methanogens tại > 300 mg/L
  • H₂S: độc với vi khuẩn hiếu khí tại > 200 mg/L dưới dạng H₂S không ion hóa
  • Kháng sinh: ức chế vi khuẩn tại µg/L, đặc biệt nguy hiểm trong nước thải bệnh viện và dược phẩm
  • Dung môi hữu cơ (toluene, benzene): ức chế màng tế bào tại g/L

Khi ngộ độc xảy ra: μm giảm → |vsu| giảm → COD tăng đột ngột dù MLVSS vẫn bình thường. Đây là dấu hiệu phân biệt ngộ độc với washout (MLVSS giảm khi washout).

Phần 5: Yếu Tố Y Và Ks — Đặc Trưng Sinh Học Không Điều Chỉnh Được Trực Tiếp

5.1 Y — Hiệu Suất Năng Lượng Của Vi Sinh

Từ phương trình 6-6: |vsu| ∝ μm/Y = k (hệ số tiêu thụ cơ chất riêng)

Khi Y nhỏ (methanogens, Y = 0,025), k = μm/Y lớn hơn tương đối so với μm — vi sinh tiêu thụ nhiều cơ chất hơn trên mỗi đơn vị sinh khối tạo thành. Nhưng vì X cũng rất nhỏ trong hệ yếm khí, tổng |vsu| = k × X × f(S) vẫn thấp hơn hiếu khí.

Vì sao không thể điều chỉnh Y? Y là đặc trưng sinh hóa cơ bản của vi sinh — phản ánh hiệu quả chuyển hóa năng lượng từ cơ chất. Kỹ sư không thể “ép” vi sinh tạo ít sinh khối hơn trên cùng lượng cơ chất (trừ khi thay đổi kiểu hô hấp — hiếu khí sang yếm khí). Chỉ có thể điều chỉnh Y_net qua SRT (phân hủy nội sinh).

5.2 Ks — Ái Lực Với Cơ Chất Và Chất Lượng Nước Đầu Ra

Ks ảnh hưởng đến |vsu| theo hai cách:

Cách 1 — Ảnh hưởng đến tốc độ tại nồng độ S thấp: Khi S < Ks: |vsu|/X ≈ k × S/Ks — vi sinh có Ks thấp hơn xử lý hiệu quả hơn ở nồng độ cơ chất thấp.

Cách 2 — Ảnh hưởng đến nồng độ đầu ra S_e: S_e tỉ lệ thuận với Ks (phân tích trong bài Monod). Vi sinh thích nghi (acclimatized) với nước thải cụ thể thường phát triển Ks thấp hơn → S_e thấp hơn → nước đầu ra sạch hơn.

Ứng dụng: Thay vì dùng bùn mồi từ nhà máy xử lý nước thải đô thị để xử lý nước thải công nghiệp đặc thù (ví dụ nước thải cà phê giàu polyphenol), nên dùng bùn đã thích nghi với nước thải cùng loại — Ks thấp hơn → hiệu suất xử lý cao hơn ngay từ đầu.

Phần 6: Ứng Dụng Tính Thể Tích Bể Từ Phương Trình 6-6

6.1 Phương Trình Thiết Kế Bể CSTR (Aerotank Khuấy Hoàn Toàn)

Trong bể CSTR ở steady state, cân bằng cơ chất cho toàn bể:

Q × (S₀ − Se) = |vsu| × V = [μm × X × Se / (Y × (Ks + Se))] × V

Kết hợp với cân bằng sinh khối (X × V/SRT = Y_net × Q × (S₀ − Se)):

Se = Ks × (1 + b_a × SRT) / (SRT × (μm − b_a) − 1)

Thể tích bể: V = Q × (S₀ − Se) × Y × (Ks + Se) / (μm × X × Se)

Ví dụ tính toán (nước thải đô thị, Q = 10.000 m³/ngày):

  • S₀ = 200 mg BOD/L; Se = 20 mg BOD/L
  • μm = 5 ngày⁻¹ (ở 28°C); Ks = 70 mg BOD/L; Y = 0,45; b_a = 0,08; SRT = 10 ngày
  • X = Y × SRT × (S₀−Se) × Q / (V × (1 + b_a × SRT)) ← cần giải hệ phương trình

Cách đơn giản hơn: HRT = V/Q = Y × SRT × (S₀−Se) / ((1 + b_a × SRT) × X_mong_muốn)

Với X = 3.000 mg VSS/L; Y_net = 0,23: HRT = 0,23 × (200−20) / (3.000/1000) = 0,23 × 180 / 3 = 13,8 giờ ≈ 14 giờ V = HRT × Q = 14/24 × 10.000 = 5.833 m³

6.2 So Sánh Thể Tích CSTR Và PFR

Với cùng hiệu suất xử lý (S₀ → Se):

  • CSTR: tốc độ xử lý = k × X × Se/(Ks + Se) (thấp — dùng tốc độ ở Se)
  • PFR: tốc độ xử lý tích phân từ S₀ đến Se — trung bình cao hơn vì S cao hơn trong phần lớn bể

Tỷ số thể tích: V_CSTR/V_PFR = [ln(S₀/Se) + Ks×(S₀−Se)/(Se×S₀)] × Se/S₀ × (Ks + S₀)/(Ks + Se)

Với S₀/Se = 10, Ks >> Se: V_CSTR/V_PFR ≈ ln(10) × Se/Ks = 2,3 × Se/Ks

Khi Se << Ks (hiệu suất cao), V_CSTR >> V_PFR — bể PFR cần ít thể tích hơn đáng kể để đạt cùng tiêu chuẩn đầu ra nghiêm ngặt. Đây là cơ sở ưu tiên dùng bể sinh học dài có nhiều ngăn (plug-flow approximation) hoặc SBR thay vì aerotank khuấy đơn thuần.

16/09/2026 9 lượt xem quantri

Có thể bạn quan tâm

Zalo
/* Chat plugin */