Phân Hủy Nội Sinh Và Tốc Độ Sinh Trưởng Thực vg' — Phương Trình vp = −kp·X Trong Mô Hình Hóa Hệ Thống Xử Lý Nước Thải

22/09/2026
3 views

Tại Sao Sinh Khối Trong Bể Không Tăng Mãi Dù Luôn Có Thức Ăn?

Một nghịch lý thoạt nhìn khó hiểu: trong bể aerotank đang vận hành ổn định với COD đầu vào đều đặn, lượng bùn hoạt tính MLVSS không tăng vô hạn mà dừng lại ở một mức cân bằng nào đó. Tại sao vậy? Nếu vi sinh vật liên tục ăn và sinh sản, tại sao không có hiện tượng "bùng nổ vi sinh"?

Câu trả lời nằm ở một quá trình song hành nhưng ngược chiều: phân hủy nội sinh (endogenous decay). Trong bất kỳ quần thể vi sinh vật nào, lúc nào cũng tồn tại đồng thời hai nhóm: nhóm đang sinh trưởng và nhóm đang già đi, chết, tự phân hủy. Khi tốc độ sinh trưởng bằng tốc độ phân hủy → quần thể ổn định. Khi phân hủy vượt sinh trưởng → quần thể suy giảm.

GS. Lê Văn Cát (2007) mô tả quá trình này qua phương trình đơn giản nhưng sâu sắc:

vp = −kp × X (Phương trình 6-9)

Và tốc độ sinh trưởng thực — tính cả phân hủy nội sinh — là:

vg' = vg + vp = μm × X × S/(Ks + S) − kp × X (Phương trình 6-10)

Đây là nền tảng để tính lượng bùn dư thực tế, nhu cầu oxy thực tế, và thiết kế SRT tối ưu trong mọi hệ thống xử lý nước thải sinh học.
 

Hệ Thống Xử Lý Nước Thải


Phần 1: Bản Chất Và Cơ Chế Phân Hủy Nội Sinh

1.1 Vi Sinh Vật Già Và Chết Là Tất Yếu

Tập đoàn vi sinh trong hệ xử lý nước thải gồm nhiều thế hệ với tuổi khác nhau. Vi sinh tuổi cao suy thoái vì:

  • Tích lũy sản phẩm độc nội bào: tế bào già tích lũy các sản phẩm trao đổi chất không thể loại ra ngoài
  • Hư hỏng DNA và protein: các đột biến tích lũy theo thời gian, enzyme hoạt tính giảm
  • Mất cân bằng áp suất thẩm thấu: màng tế bào già kém bền vững hơn
  • Cạnh tranh nội bộ: vi sinh trẻ hoạt tính cao hơn cạnh tranh cơ chất tốt hơn vi sinh già

Khi tế bào chết, phần lớn vật chất bên trong bị thủy phân bởi enzyme ngoại bào → trở thành cơ chất cho vi sinh còn sống (quá trình "tái chế" vật chất trong hệ). Một phần nhỏ (~15–20% VSS) chuyển thành chất trơ không phân hủy được, tích lũy thành bùn già.

1.2 Phân Biệt Phân Hủy Nội Sinh Với Các Quá Trình Khác

Quá Trình

Cơ Chế

Điều Khiển Được?

Phân hủy nội sinh

Già hóa và chết tự nhiên

Không (chỉ giảm qua nhiệt độ thấp)

Xả bùn dư

Loại bỏ chủ động bùn từ hệ thống

Có (điều chỉnh SRT)

Phân hủy sinh học

Vi sinh phân hủy chất hữu cơ nước thải

Không (phụ thuộc cơ chất)

Lysis do predation

Động vật nguyên sinh ăn vi khuẩn

Một phần (giảm HRT)

Phân hủy nội sinh (endogenous decay) là quá trình sinh hóa tự nhiên không phải do vận hành chủ động — điều này phân biệt nó với việc xả bùn dư có kiểm soát.


Phần 2: Phương Trình vp = −kp·X — Phân Tích Và Giá Trị Thực Tế

2.1 Ý Nghĩa Vật Lý

vp (g VSS/L·ngày): tốc độ mất đi sinh khối do phân hủy nội sinh — đại lượng âm.

kp (ngày⁻¹): hằng số phân hủy nội sinh — tỷ lệ % sinh khối bị phân hủy mỗi ngày.

X (g VSS/L): nồng độ sinh khối vi sinh.

Phương trình vp = −kp × X mô tả động học bậc một: tốc độ phân hủy tỉ lệ thuận với lượng sinh khối hiện có — tương tự động học phân hủy phóng xạ. Đây là hệ quả tự nhiên: càng nhiều tế bào → càng nhiều tế bào già chết đi mỗi ngày.

Tốc độ suy giảm 1%: nếu kp = 0,01 ngày⁻¹, thì mỗi ngày 1% sinh khối bị phân hủy. Sau 100 ngày không có sinh trưởng: X còn lại = X₀ × e^(−0,01×100) = X₀ × 0,368 = 36,8%.

2.2 Giá Trị kp Theo Điều Kiện Và Loại Vi Sinh

Vi Sinh / Điều Kiện

kp (ngày⁻¹)

Ghi Chú

Vi khuẩn dị dưỡng, hiếu khí, 20°C

0,05–0,12

Giá trị thường dùng: 0,08 ngày⁻¹

Vi khuẩn dị dưỡng, hiếu khí, 30°C

0,08–0,18

Tăng theo nhiệt độ

Vi khuẩn dị dưỡng, thiếu khí

0,03–0,08

Thấp hơn hiếu khí 30–50%

Nitrosomonas (AOB)

0,04–0,08

Dùng 0,05 ngày⁻¹

Nitrobacter (NOB)

0,02–0,05

Thấp hơn AOB

Methanogens

0,01–0,03

Rất chậm — bền vững

Anammox

0,003–0,010

Cực chậm — đặc điểm quan trọng

Nhận xét quan trọng:

  • Hiếu khí (kp = 0,08) phân hủy nhanh hơn yếm khí (methanogens, kp = 0,02) khoảng 4 lần
  • Anammox (kp = 0,005) cực kỳ bền vững → lý do khi khởi động lại hệ Anammox sau nghỉ dài vẫn có thể phục hồi

2.3 Ảnh Hưởng Nhiệt Độ Lên kp

Kp cũng thay đổi theo nhiệt độ, nhưng với hệ số θ nhỏ hơn μm: kp(T) = kp(20°C) × θ_kp^(T-20) với θ_kp ≈ 1,04–1,06

Tại 30°C: kp(30°C) = kp(20°C) × 1,05^10 = kp(20°C) × 1,63

Nghĩa là vào mùa hè (nước thải 30°C), phân hủy nội sinh nhanh hơn 63% so với mùa đông (20°C). Cùng với việc μm cũng tăng, nhưng kp tăng ít hơn → vg' = vg − kp × X tăng ròng → bể sản xuất nhiều bùn hơn vào mùa hè.


Phần 3: Phương Trình Tốc Độ Sinh Trưởng Thực vg' — 6-10 Và 6-11

3.1 Phương Trình 6-10 — Dạng Khai Triển

vg' = μm × X × S/(Ks + S) − kp × X

Đây là tốc độ tăng sinh khối thực tế trong bể — là hiệu giữa hai quá trình đối nghịch:

  • μm × X × S/(Ks + S): tốc độ sinh trưởng thuần từ Monod (vg)
  • kp × X: tốc độ mất đi do phân hủy nội sinh (|vp|)

Phương trình 6-10 có thể được biến đổi: vg' = X × [μm × S/(Ks + S) − kp] = X × [μ − kp] = X × μ'

Trong đó μ' = μ − kp là hằng số phát triển riêng thực (sẽ phân tích trong bài tiếp theo).

3.2 Phương Trình 6-11 — Dạng Kết Nối Với Cơ Chất

Thay vg = −Y × vsu vào phương trình 6-10: vg' = −Y × vsu − kp × X (Phương trình 6-11)

Đây là dạng hữu ích nhất để tính toán vì:

  • vsu có thể tính từ tải BOD và điều kiện bể
  • Y và kp là hằng số (tra bảng hoặc đo thực nghiệm)
  • X là biến cần theo dõi vận hành (đo MLVSS)

Ví dụ tính toán (aerotank ổn định, không xả bùn, chỉ có phân hủy nội sinh):

Giả sử: X = 3.000 mg VSS/L; vsu = −200 mg BOD/L·ngày; Y = 0,45; kp = 0,08 ngày⁻¹

  • Sinh trưởng thuần (vg): Y × |vsu| = 0,45 × 200 = 90 mg VSS/L·ngày
  • Phân hủy nội sinh (|vp|): kp × X = 0,08 × 3.000 = 240 mg VSS/L·ngày
  • Tốc độ thực (vg'): 90 − 240 = −150 mg VSS/L·ngày (âm! sinh khối đang giảm)

Kết quả âm này có nghĩa: bể đang tiêu thụ BOD nhưng vẫn giảm sinh khối vì phân hủy nội sinh vượt sinh trưởng. Điều này xảy ra khi tải BOD quá thấp so với lượng bùn hiện có — bùn đang "ăn chính mình" (endogenous respiration phase).

3.3 Điều Kiện Cân Bằng Sinh Khối (Steady State)

Ở điều kiện ổn định trong hệ thống liên tục (bể có xả bùn dư), vg' thực (sau khi tính thêm xả bùn) = 0:

vg_thực_vận_hành = vg' − (X × V / SRT) / V = 0

vg' = X/SRT

Thay phương trình 6-10: X × μ' = X/SRT → μ' = 1/SRT

Đây là phương trình thiết kế cơ bản: SRT được xác định bởi hằng số phát triển riêng thực μ', không phải μ. Khi thiết kế hệ nitrification cần SRT đủ để μ'_Nitrosomonas ≥ 1/SRT_thiết_kế.


Phần 4: Hệ Quả Của Phân Hủy Nội Sinh Đối Với Vận Hành

4.1 Nhu Cầu Oxy Từ Hô Hấp Nội Sinh

Phân hủy nội sinh không chỉ giảm MLVSS mà còn tiêu thụ oxy. Từ phản ứng: C₅H₇O₂N + 5O₂ → 5CO₂ + NH₄⁺ + 2H₂O

Cứ 1 g VSS bị phân hủy nội sinh cần 1,42 g O₂.

Tốc độ tiêu thụ oxy từ hô hấp nội sinh: OUR_nội_sinh = 1,42 × kp × X (g O₂/L·ngày)

Ví dụ (bể aerotank 5.000 m³, X = 3.500 mg VSS/L, kp = 0,08 ngày⁻¹): OUR_nội_sinh = 1,42 × 0,08 × 3.500 = 396 mg O₂/L·ngày Tổng O₂ nội sinh = 396 × 5.000 = 1.980.000 g O₂/ngày = 1.980 kg O₂/ngày

Đây là lượng oxy đáng kể — trong nhiều trường hợp chiếm 20–40% tổng nhu cầu oxy của bể. Nếu không tính hô hấp nội sinh → thiết kế máy thổi khí sẽ thiếu công suất!

Tổng nhu cầu oxy (AOR): AOR = O₂_oxy_hóa_COD + O₂_nội_sinh = Q × ΔCOD × (1 − p × Y_net) + 1,42 × kp × Xa × V

4.2 Bùn Trơ Tích Lũy Từ Phân Hủy Nội Sinh

Không phải toàn bộ sinh khối bị phân hủy nội sinh được tái sử dụng. Một phần f_e (thường 0,15–0,20 g VSS trơ/g VSS chết) chuyển thành chất trơ sinh học X_e:

Tốc độ tạo X_e = f_e × kp × X

Với f_e = 0,15; kp = 0,08; X = 3.500 mg/L: Tốc độ tạo bùn trơ = 0,15 × 0,08 × 3.500 = 42 mg X_e/L·ngày

Sau thời gian dài (SRT = 20 ngày), tổng X_e trong bể: X_e / X_total ≈ f_e × kp × SRT / (1 + f_e × kp × SRT) = 0,15 × 0,08 × 20 / (1 + 0,24) = 0,24/1,24 = 19,4% tổng VSS

Gần 20% VSS trong bể là chất trơ không đóng góp vào xử lý COD — "vốn chết" trong hệ thống. SRT càng dài, tỷ lệ bùn trơ càng cao.

4.3 Tác Động Của kp Đến Lượng Bùn Dư

Phương trình sinh khối dư (Px) có tính đến phân hủy nội sinh:

Px = [Y × Q × ΔCOD − kp × X_a × V] / (p × 1000)

Với Y_net = Y/(1 + kp × SRT) đã tích hợp sẵn kp.

So sánh có và không có phân hủy nội sinh (Q = 10.000 m³/ngày, ΔCOD = 300 mg/L, Y = 0,45, kp = 0,08, SRT = 15 ngày):

  • Không có kp (lý thuyết thuần túy): Px = 0,45 × 10.000 × 0,30 = 1.350 kg VSS/ngày
  • Có kp (thực tế): Y_net = 0,45/(1 + 0,08 × 15) = 0,45/2,2 = 0,205 Px = 0,205 × 10.000 × 0,30 = 615 kg VSS/ngày — chỉ bằng 45,6% lý thuyết!

Kết luận: phân hủy nội sinh làm giảm bùn dư thực tế gần một nửa so với lý thuyết không có kp — đây là lý do tại sao SRT dài (extended aeration) tạo rất ít bùn dư.


Phần 5: So Sánh kp Giữa Các Công Nghệ Và Điều Kiện Khí Hậu Việt Nam

5.1 Extended Aeration — Giải Pháp Giảm Bùn Tối Đa

Extended aeration (thổi khí kéo dài) dùng SRT = 20–30 ngày để tối đa hóa phân hủy nội sinh → tối thiểu hóa bùn dư:

Y_net (SRT = 25 ngày) = 0,45/(1 + 0,08 × 25) = 0,45/3,0 = 0,15 g VSS/g COD

So với SRT = 5 ngày: Y_net = 0,45/1,4 = 0,32 g VSS/g COD

Extended aeration tạo ra ít hơn 53% lượng bùn so với aerotank SRT ngắn. Đánh đổi: tiêu thụ nhiều O₂ hơn (hô hấp nội sinh tăng) và thể tích bể lớn hơn.

Phù hợp với Việt Nam: Extended aeration đặc biệt thích hợp cho các khu dân cư nhỏ và đô thị vừa (công suất < 5.000 m³/ngày) do:

  • Không cần bể phân hủy bùn riêng (bùn đã được phân hủy ngay trong bể)
  • Bùn ra khỏi hệ đã ổn định, có thể làm phân bón trực tiếp
  • Vận hành đơn giản, ít sự cố

Nhà máy XLNT Bình Phước, Tây Ninh, nhiều thị xã miền Tây Nam Bộ đang áp dụng extended aeration thành công.

5.2 Ảnh Hưởng Nhiệt Độ Cao Tại Miền Nam Việt Nam

Tại nhiệt độ nước thải 30°C:

  • kp(30°C) = 0,08 × 1,05^10 = 0,130 ngày⁻¹ (tăng 63% so với 20°C)
  • μm(30°C) = 5,0 × 1,07^10 = 9,9 ngày⁻¹ (tăng 98%)

Vì μm tăng nhiều hơn kp, tỷ lệ kp/μm giảm từ 0,08/5,0 = 0,016 xuống 0,130/9,9 = 0,013 → vi sinh hiệu quả hơn (ít "lãng phí" cho phân hủy hơn). Nhiệt độ cao có lợi cho xử lý sinh học theo nghĩa này.

Tuy nhiên, kp cao hơn vẫn có nghĩa là bùn dư ít hơn ở nhiệt độ cao → nhà máy tại miền Nam sẽ ít bùn dư hơn nhà máy tại miền Bắc cùng quy mô và SRT.


Phần 6: Đo Lường kp Tại Hiện Trường

6.1 Phương Pháp Hô Hấp Nội Sinh (Endogenous Respiration Rate)

kp đo trực tiếp bằng cách đo tốc độ tiêu thụ oxy khi không có cơ chất bên ngoài:

  1. Lấy bùn hoạt tính, rửa sạch loại bỏ COD hòa tan
  2. Đặt vào bình kín có đo oxy (respirometer)
  3. Không thêm cơ chất → tốc độ tiêu thụ O₂ = OUR_endogenous
  4. kp = OUR_endogenous / (1,42 × X)

Ví dụ thực tế: OUR_endogenous = 30 mg O₂/L·giờ = 720 mg O₂/L·ngày; X = 4.500 mg VSS/L: kp = 720 / (1,42 × 4.500) = 720/6.390 = 0,113 ngày⁻¹ (cao, bùn khỏe mạnh ở 30°C)

6.2 Phương Pháp Cân Bằng Khối (Mass Balance)

Đơn giản hơn nhưng kém chính xác hơn — tính từ dữ liệu vận hành:

kp = (Y × Q × ΔCOD/V − Px/V) / X

Trong đó Px là lưu lượng bùn dư thực tế (đo được từ máy bơm bùn dư và phân tích VSS bùn).

Phương pháp này tích hợp tốt vào hệ thống SCADA của các nhà máy lớn — theo dõi kp theo thời gian thực để phát hiện bùn lão hóa (kp giảm bất thường → vi sinh đang bị ức chế hoặc "ngủ đông").


Kết Luận

Phân hủy nội sinh với phương trình vp = −kp × X là quá trình không thể tránh khỏi trong bất kỳ hệ thống vi sinh nào. Tích hợp nó vào mô hình qua vg' = vg + vp cho phép dự đoán chính xác hơn về lượng bùn dư, nhu cầu oxy, và điều kiện ổn định hệ thống.

Bốn điểm thực hành cốt lõi:

Một, kp = 0,05–0,12 ngày⁻¹ cho vi khuẩn dị dưỡng hiếu khí — luôn dùng kp khi tính Y_net, bùn dư, và nhu cầu oxy.

Hai, hô hấp nội sinh tiêu thụ 1,42 g O₂/g VSS — lên đến 40% tổng nhu cầu oxy của bể aerotank, không thể bỏ qua trong thiết kế máy thổi khí.

Ba, SRT dài (extended aeration, SRT = 20–30 ngày) tối đa hóa phân hủy nội sinh → giảm bùn dư 50% so với SRT ngắn — phù hợp với các nhà máy Việt Nam muốn đơn giản hóa xử lý bùn.

Bốn, đo kp tại hiện trường bằng OUR_endogenous là cách kiểm tra sức khỏe bùn hoạt tính tin cậy nhất — kp bất thường báo hiệu ngộ độc, ức chế, hay bùn lão hóa.

 

Bình luận facebook