Mô Hình McCarty Trong Xử Lý Nước Thải: R = fs·Rs + fa·Ra − Rd, Bảng 5.5 Và Tính Toán Hiệu Suất Sinh Khối 53%

Mô Hình Toán Học Hóa Học Của Quá Trình Sinh Hóa Trong Xử Lý Nước Thải

Trong kỹ thuật xử lý nước thải, hầu hết các thiết kế đều dựa trên thực nghiệm và kinh nghiệm vận hành. Tuy nhiên, có một công cụ lý thuyết mạnh mẽ ít được biết đến ở Việt Nam cho phép tính toán hoàn chỉnh mọi phản ứng sinh hóa trong hệ xử lý — bao gồm lượng sinh khối tạo ra, lượng O₂ hoặc NO₃⁻ tiêu thụ, và lượng CO₂ thải ra — chỉ từ các thông số đầu vào cơ bản.

Đó là mô hình McCarty (Perry McCarty, Stanford University), được GS. Lê Văn Cát trình bày trong chương 5 của giáo trình Xử Lý Nước Thải Giàu Hợp Chất Nitơ Và Phốtpho (2007) qua phương trình:

R = fs · Rs + fa · Ra − Rd (phương trình 5-15)

Kết hợp với Bảng 5.5 chứa toàn bộ các nửa phản ứng thành phần, mô hình này cho phép kỹ sư xây dựng phương trình sinh hóa hoàn chỉnh cho bất kỳ quá trình nào trong hệ xử lý nước thải. Bài viết này giải thích chi tiết nguyên lý mô hình, cấu trúc Bảng 5.5 và minh họa ứng dụng tính hiệu suất sinh khối 53% từ glucose.

1. Nguyên Lý Mô Hình McCarty — R = fs · Rs + fa · Ra − Rd

1.1. Ý Nghĩa Từng Thành Phần

R — Phản ứng tổng thể: Phương trình hóa học hoàn chỉnh mô tả toàn bộ quá trình sinh hóa — đầu vào (cơ chất, chất nhận điện tử, nguồn N) và đầu ra (sinh khối, sản phẩm cuối, nước).

Rd — Phản ứng cho điện tử (Donor reaction): Nửa phản ứng oxy hóa của cơ chất (chất cho điện tử). Ví dụ: oxy hóa glucose, acetate, methanol, H₂, NH₄⁺, S²⁻… Quy ước: viết theo chiều cơ chất → sản phẩm oxy hóa + H⁺ + e⁻

Ra — Phản ứng nhận điện tử (Acceptor reaction): Nửa phản ứng khử của chất nhận điện tử. Ví dụ: khử O₂, NO₃⁻, SO₄²⁻, CO₂… Quy ước: viết theo chiều chất nhận điện tử → sản phẩm khử

Rs — Phản ứng tổng hợp tế bào (Synthesis reaction): Nửa phản ứng tổng hợp sinh khối từ nguồn carbon và nguồn N. Công thức tế bào: C₅H₇O₂N

fs — Phần điện tử dành cho tổng hợp tế bào: Tỷ lệ mol điện tử (từ Rd) được Vi sinh vật dùng để tổng hợp sinh khối (đồng hóa).

fa — Phần điện tử dành cho năng lượng: Tỷ lệ mol điện tử (từ Rd) được dùng để thu năng lượng (dị hóa) thông qua Ra.

Ràng buộc quan trọng: fs + fa = 1

Dấu trừ trước Rd: Do Rd viết theo chiều oxy hóa (cơ chất → sản phẩm + H⁺ + e⁻), khi tổ hợp vào phương trình tổng R phải đảo chiều Rd. Dấu trừ trước Rd trong công thức R = fs·Rs + fa·Ra − Rd thể hiện điều này.

1.2. Ý Nghĩa Sinh Học Của fs Và fa

fs và fa là phân bổ nguồn điện tử giữa hai nhiệm vụ căn bản của vi sinh vật:

• Tồn tại = cần năng lượng → dùng fa để dị hóa
• Sinh trưởng = cần vật liệu tế bào → dùng fs để đồng hóa

Vi sinh vật tối ưu hóa tỷ lệ fs/fa theo điều kiện môi trường. Khi nguồn dinh dưỡng dồi dào và cơ chất dễ phân hủy: fs cao (sinh trưởng nhiều). Khi thiếu dinh dưỡng hoặc cơ chất khó phân hủy: fa cao (tiêu tốn nhiều năng lượng để “xử lý” cơ chất, ít sinh trưởng).

Giá trị fs điển hình theo loại phản ứng:

• Carbohydrat + O₂ (dị dưỡng hiếu khí, dễ phân hủy): fs = 0,60–0,72
• Carbohydrat + NO₃⁻ (khử nitrat): fs = 0,53–0,65 (thấp hơn vì ít năng lượng hơn)
• NH₄⁺ + O₂ (tự dưỡng nitrat hóa): fs = 0,08–0,14 (rất thấp vì ít năng lượng)
• Acetate + O₂: fs = 0,55–0,65
• Methanol + O₂: fs = 0,45–0,55

2. Bảng 5.5 — Thư Viện Nửa Phản Ứng Sinh Hóa

2.1. Cấu Trúc Bảng 5.5

Bảng 5.5 trong giáo trình GS. Lê Văn Cát chứa đầy đủ các nửa phản ứng thành phần chuẩn hóa theo 1 mol điện tử. Đây là yếu tố then chốt: tất cả nửa phản ứng đều được viết dạng chứa “H⁺ + e⁻” với hệ số 1, cho phép cộng trừ trực tiếp khi tổ hợp.

Nhóm 1 — Phản ứng tổng hợp tế bào (Rs): 2 biến thể

Rs với NH₄⁺ làm nguồn N (số 1):

(1/20)CO₂ + (1/20)HCO₃⁻ + (1/20)NH₄⁺ + H⁺ + e⁻ → (1/20)C₅H₇O₂N + (9/20)H₂O

Rs với NO₃⁻ làm nguồn N (số 2):

(1/28)NO₃⁻ + (5/28)CO₂ + (29/28)H⁺ + e⁻ → (1/28)C₅H₇O₂N + (11/28)H₂O

Việc có 2 biến thể Rs phản ánh thực tế: nếu môi trường không có NH₄⁺, vi khuẩn có thể dùng NO₃⁻ làm nguồn N cho tổng hợp tế bào (tốn thêm năng lượng để khử NO₃⁻ → NH₄⁺ → tế bào).

Nhóm 2 — Phản ứng nhận điện tử (Ra): 4 loại chất nhận

STT	Chất nhận điện tử	ΔG° (kJ/mol e⁻)
3	O₂: 0,25O₂ + H⁺ + e⁻ → 0,5H₂O	−78,14
4	NO₃⁻: 0,2NO₃⁻ + 1,2H⁺ + e⁻ → 0,1N₂ + 0,6H₂O	−71,67
5	SO₄²⁻: 0,125SO₄²⁻ + 1,1875H⁺ + e⁻ → 0,0625H₂S + 0,0625HS⁻ + 0,5H₂O	+21,27
6	CO₂ (methanogenesis): 0,125CO₂ + H⁺ + e⁻ → 0,125CH₄ + 0,25H₂O	+24,11

Nhóm 3 — Phản ứng cho điện tử (Rd): Chất hữu cơ (dị dưỡng)

STT	Cơ chất	ΔG° (kJ/mol e⁻)
7	Nước thải sinh hoạt (C₁₀H₁₉O₃N đại diện)	+31,80
8	Protein/axit amin (C₁₆H₂₄O₅N đại diện)	+32,22
9	Carbohydrat (CH₂O đại diện)	+41,84
10	Dầu mỡ (C₈H₁₆O đại diện)	+27,61
11	Methanol (CH₃OH)	+37,51

Nhóm 4 — Phản ứng cho điện tử (Rd): Chất vô cơ (tự dưỡng)

STT	Cơ chất vô cơ	ΔG° (kJ/mol e⁻)
12	Fe²⁺ → Fe³⁺	−74,40
13	NO₂⁻ → NO₃⁻	−40,15
14	NH₄⁺ → NO₃⁻	−34,50
15	NH₄⁺ → NO₂⁻	−32,62
16–20	Các dạng oxy hóa S	+19–44

Nguồn: Bảng 5.5, GS. Lê Văn Cát, 2007

2.2. Quy Ước Dấu Trong Bảng 5.5

Tất cả nửa phản ứng trong Bảng 5.5 được viết theo chiều:

Chất oxy hóa + H⁺ + e⁻ → Chất khử (chiều khử)

Khi dùng làm Ra (phản ứng nhận điện tử trong R = fs·Rs + fa·Ra − Rd): dùng trực tiếp chiều này.

Khi dùng làm Rd (phản ứng cho điện tử): đảo chiều — đây là lý do có dấu “−” trước Rd trong công thức.

3. Ví Dụ Tính Toán: Hiệu Suất Sinh Khối 53% Từ Oxy Hóa Glucose

GS. Lê Văn Cát trình bày ví dụ chi tiết tính hiệu suất hình thành sinh khối khi oxy hóa glucose (carbohydrat) với oxy.

3.1. Điều Kiện Bài Toán

• Cơ chất (Rd): Carbohydrat CH₂O (glucose) — phản ứng số 9 trong Bảng 5.5
• Chất nhận điện tử (Ra): O₂ — phản ứng số 3
• Nguồn N (Rs): NH₄⁺ — phản ứng số 1
• fs = 0,70 (70% điện tử dùng cho tổng hợp tế bào)
• fa = 1 − fs = 0,30 (30% điện tử dùng cho năng lượng)

3.2. Các Nửa Phản Ứng Thành Phần

Rs (phản ứng 1 — tổng hợp tế bào với NH₄⁺):

(1/20)CO₂ + (1/20)HCO₃⁻ + (1/20)NH₄⁺ + H⁺ + e⁻ → (1/20)C₅H₇O₂N + (9/20)H₂O

Ra (phản ứng 3 — khử O₂):

(1/4)O₂ + H⁺ + e⁻ → (1/2)H₂O

Rd (phản ứng 9 — oxy hóa carbohydrat — viết theo chiều ngược trong Bảng 5.5):

(1/4)CO₂ + H⁺ + e⁻ → (1/4)CH₂O + (1/4)H₂O (chiều trong bảng)

Đảo chiều để dùng làm Rd:

(1/4)CH₂O + (1/4)H₂O → (1/4)CO₂ + H⁺ + e⁻

3.3. Tổ Hợp Theo R = fs·Rs + fa·Ra − Rd

Bước 1: Tính fs·Rs = 0,70·Rs

0,70 × [(1/20)CO₂ + (1/20)HCO₃⁻ + (1/20)NH₄⁺ + H⁺ + e⁻ → (1/20)C₅H₇O₂N + (9/20)H₂O]

= 0,035CO₂ + 0,035HCO₃⁻ + 0,035NH₄⁺ + 0,70H⁺ + 0,70e⁻ → 0,035C₅H₇O₂N + 0,315H₂O

Bước 2: Tính fa·Ra = 0,30·Ra

0,30 × [(1/4)O₂ + H⁺ + e⁻ → (1/2)H₂O]

= 0,075O₂ + 0,30H⁺ + 0,30e⁻ → 0,150H₂O

Bước 3: Tính −Rd (đảo chiều phản ứng 9)

−Rd có nghĩa: đảo chiều phương trình Rd ban đầu:

Rd gốc (chiều cho e⁻):

(1/4)CH₂O + (1/4)H₂O → (1/4)CO₂ + H⁺ + e⁻

Khi viết −Rd trong tổ hợp: đưa Rd về vế trái của phản ứng tổng, tức là cộng:

• (1/4)CH₂O + (1/4)H₂O = (1/4)CO₂ + H⁺ + e⁻ → đưa vế trái thành đầu vào của R

Hay viết: 0,25CH₂O + 0,25H₂O → 0,25CO₂ + H⁺ + e⁻ (chiều Rd thực tế)

3.4. Phương Trình Tổng R (Phản Ứng 5-17 Trong Giáo Trình)

Tổ hợp ba thành phần, cân bằng H⁺ và e⁻ (triệt tiêu nhau):

Đầu vào:

• Từ Rd: 0,25CH₂O + 0,25H₂O
• Từ fs·Rs: 0,035HCO₃⁻ + 0,035NH₄⁺
• Từ fa·Ra: 0,075O₂

Đầu ra:

• Từ fs·Rs: 0,035C₅H₇O₂N + 0,315H₂O
• Từ Rd → 0,25CO₂ (đưa về phải)
• Từ Rd → phần CO₂ còn lại sau khử trừ: 0,25CO₂ − 0,035CO₂ = 0,215CO₂

Phương trình tổng (GS. Lê Văn Cát viết):

0,25CH₂O + 0,035HCO₃⁻ + 0,035NH₄⁺ + 0,075O₂ → 0,035C₅H₇O₂N + 0,11CO₂ + 0,21H₂O (phản ứng 5-17)

3.5. Tính Hiệu Suất Sinh Khối 53%

Lượng carbohydrat (CH₂O) tiêu thụ:

n(CH₂O) = 0,25 mol × M(CH₂O) = 0,25 × 30 = 7,5 g

Lượng sinh khối (C₅H₇O₂N) tạo ra:

n(C₅H₇O₂N) = 0,035 mol × M(C₅H₇O₂N) = 0,035 × 113 = 3,96 g

Hiệu suất hình thành tế bào:

H = 3,96/7,5 = 0,528 ≈ 53%

GS. Lê Văn Cát kết luận: “hiệu suất hình thành tế bào là 53% (3,96/7,5).”

Liên hệ với hệ số năng suất Y:

Y = 3,96 g VSS / 7,5 g CH₂O

Nếu COD của CH₂O = 1,07 g COD/g CH₂O (tính theo C₆H₁₂O₆ + 6O₂ → 6CO₂ + 6H₂O):

Y = 3,96 g VSS / (7,5 × 1,07 g COD) = 3,96/8,03 = 0,49 g VSS/g COD

Kết quả 0,49 g VSS/g COD hoàn toàn phù hợp với giá trị Y = 0,4–0,6 thực nghiệm thường dùng trong thiết kế aerotank.

4. Ứng Dụng Mô Hình McCarty Cho Các Quá Trình Khác

4.1. Nitrat Hóa — Tự Dưỡng Chemoautotroph Với fs Thấp

Tính hiệu suất sinh khối trong nitrat hóa NH₄⁺ → NO₂⁻ bởi Nitrosomonas, với:

• Rd: số 15 trong Bảng 5.5 — oxy hóa NH₄⁺ → NO₂⁻
• Ra: số 3 — khử O₂
• Rs: số 1 — tổng hợp tế bào với NH₄⁺
• fs = 0,10 (giá trị điển hình cho tự dưỡng hóa năng)

fs·Rs (0,10·Rs):

0,005CO₂ + 0,005HCO₃⁻ + 0,005NH₄⁺ + 0,10H⁺ + 0,10e⁻ → 0,005C₅H₇O₂N + 0,045H₂O

fa·Ra (0,90·Ra):

0,225O₂ + 0,90H⁺ + 0,90e⁻ → 0,45H₂O

−Rd (đảo phản ứng 15): Phản ứng 15 trong bảng (chiều khử): (1/6)NO₂⁻ + (4/3)H⁺ + e⁻ → (1/6)NH₄⁺ + (1/3)H₂O

Đảo chiều: (1/6)NH₄⁺ + (1/3)H₂O → (1/6)NO₂⁻ + (4/3)H⁺ + e⁻

Tổ hợp R = 0,10·Rs + 0,90·Ra − Rd (với 1 mol e⁻):

(1/6)NH₄⁺ + 0,225O₂ → (1/6)NO₂⁻ + 0,005C₅H₇O₂N + 0,005HCO₃⁻ + sản phẩm

Tính Y:

• Sinh khối: 0,005 mol × 113 g/mol = 0,565 g sinh khối
• NH₄⁺-N oxy hóa: (1/6) mol × 14 g/mol = 2,33 g N
• Y = 0,565/2,33 = 0,24 g VSS/g NH₄⁺-N (rất thấp vì fs thấp)

So với giá trị thực nghiệm Y(Nitrosomonas) = 0,04–0,13 g VSS/g N — kết quả từ mô hình hơi cao do mô hình McCarty chưa tính phân hủy nội sinh, nhưng cho thấy xu hướng đúng.

4.2. Khử Nitrat Dị Dưỡng Với Methanol

Nhiều nhà máy xử lý nước thải tại Việt Nam bổ sung methanol làm nguồn carbon cho bể anoxic. Tính R:

• Rd: số 11 — oxy hóa methanol CH₃OH
• Ra: số 4 — khử NO₃⁻
• Rs: số 1 — tổng hợp tế bào với NH₄⁺
• fs = 0,50

Giá trị ΔG° tổng:

ΔG° = 0,50 × ΔG°(Rs) + 0,50 × ΔG°(Ra) − ΔG°(Rd) = 0,50 × ~20 + 0,50 × (−71,67) − 37,51 = 10 − 35,84 − 37,51 = −63,35 kJ/mol e⁻

ΔG° < 0 → phản ứng tự xảy ra. Năng lượng thu được −63,35 kJ/mol e⁻ là cơ sở sinh trưởng của vi khuẩn khử nitrat khi dùng methanol.

5. Từ Mô Hình McCarty Đến Thiết Kế Hệ Xử Lý Thực Tế

5.1. Xác Định Nhu Cầu Oxy Dựa Trên fs

Từ phương trình R đã tính, có thể xác định hệ số tiêu thụ O₂:

Trong ví dụ glucose: 0,075 mol O₂ được tiêu thụ để oxy hóa 0,25 mol CH₂O = 7,5 g CH₂O.

a = lượng O₂ / lượng cơ chất = (0,075 × 32) / 7,5 = 2,4/7,5 = 0,32 g O₂/g CH₂O

Hay theo COD: 0,32/(1,07 g COD/g CH₂O) = 0,30 g O₂/g COD

Tổng nhu cầu O₂ gồm cả phân hủy nội sinh (bₐ):

O₂ tổng = a × ΔCOD + 1,42 × bₐ × X × V

5.2. Tác Động Của fs Đến Lượng Bùn Dư

Kết quả từ mô hình McCarty xác nhận nguyên lý thiết kế quan trọng:

• fs cao (= Y cao) → nhiều sinh khối tạo ra → nhiều bùn dư → chi phí xử lý bùn cao
• fs thấp (= Y thấp) → ít bùn → nhưng vi sinh vật sinh trưởng chậm hơn, cần SRT dài hơn

Điều này dẫn đến quyết định thiết kế chiến lược:

Công nghệ	fs điển hình	Y (g VSS/g COD)	Lượng bùn	Đặc điểm
Aerotank SRT ngắn (3–5 ngày)	0,60–0,70	0,40–0,55	Nhiều	Dễ vận hành, nhiều bùn
Aerotank SRT dài (15–30 ngày)	0,45–0,55	0,25–0,35	Ít hơn	Khó vận hành hơn, ít bùn
MBR (SRT > 20 ngày)	0,40–0,50	0,20–0,30	Rất ít	Chi phí màng cao
UASB (Methanogens)	0,05–0,08	0,02–0,05	Rất ít	Chậm, ít bùn, sinh biogas

5.3. Ứng Dụng Tại Nhà Máy Xử Lý Nước Thải Công Nghiệp Việt Nam

Ví dụ thực tế: Nhà máy xử lý nước thải khu công nghiệp Bình Dương

Nước thải hỗn hợp từ nhiều ngành công nghiệp:

• COD vào = 800 mg/L (đại diện: nước thải sinh hoạt, C₁₀H₁₉O₃N)
• Q = 10.000 m³/ngày
• Yêu cầu: COD đầu ra < 75 mg/L (QCVN 40:2011 cột B)

Áp dụng mô hình McCarty với Rd = nước thải sinh hoạt (số 7, ΔG°=+31,80), Ra = O₂ (số 3), Rs = NH₄⁺ (số 1), fs = 0,65:

Từ đó tính:

• Y = ~0,50 g VSS/g COD
• a (O₂/COD) = ~0,52 g O₂/g COD
• Bùn dư = 0,50 × (800−75) mg/L × 10.000 m³/ngày × 10⁻³ = 3.625 kg VSS/ngày
• O₂ cần = 0,52 × 725 mg/L × 10.000 m³/ngày × 10⁻³ = 3.770 kg O₂/ngày

Với thổi khí hiệu suất 1,5 kg O₂/kWh:

Điện thổi khí = 3.770/1,5 = 2.513 kWh/ngày = 104 kW công suất trung bình

6. Kiểm Chứng Năng Lượng Tự Do Gibbs Của Phản Ứng 5-17

Từ ví dụ glucose, GS. Lê Văn Cát tính:

ΔG°tổng = ΔG°(Ra) + ΔG°(Rd ngược) = −78,14 + 41,84 = −36,3 kJ/mol e⁻

(Lưu ý: phần Rs không đóng góp năng lượng vào ΔG° tổng vì Rs là phản ứng tổng hợp sinh khối — vi sinh vật không “thu” năng lượng từ đây mà tiêu tốn năng lượng để tổng hợp tế bào)

Giá trị −36,3 kJ/mol e⁻ xác nhận:

• Phản ứng tự xảy ra (ΔG° < 0)
• Cân bằng hóa học dịch mạnh về phía phải (tuyệt đối lớn)
• Tương đương năng lượng: −36,3 kJ/mol e⁻ × 24 mol e⁻/mol glucose = −871 kJ/mol glucose

So sánh với giá trị thực nghiệm: nhiệt đốt cháy glucose (có tính đến hiệu ứng nhiệt) ≈ −2.810 kJ/mol, nhưng trong điều kiện sinh học với hiệu suất chuyển hóa thực tế, năng lượng thu được ước tính −800 đến −900 kJ/mol glucose — phù hợp hoàn toàn.

7. Giới Hạn Và Lưu Ý Khi Áp Dụng Mô Hình McCarty

7.1. Mô Hình Cho Kết Quả Định Tính/Bán Định Lượng

GS. Lê Văn Cát nhấn mạnh: sử dụng bảng 5.5 và biểu thức 5-15 “có thể đánh giá một phản ứng sinh hóa trong xử lý nước thải (sơ bộ, mang tính chất định tính).”

Mô hình McCarty có các giới hạn:

• Không tính phân hủy nội sinh (bₐ)
• fs là tham số trung bình, thực tế thay đổi theo điều kiện SRT, pH, nhiệt độ
• Không phân biệt vi khuẩn sống/chết trong sinh khối
• Không tính đến sự cạnh tranh giữa các nhóm vi sinh vật

7.2. Khi Nào Dùng Mô Hình McCarty?

• Nên dùng: Ước tính sơ bộ Y, nhu cầu O₂, lượng bùn dư, kiểm tra cân bằng vật chất; xây dựng phương trình stoichiometry cho phản ứng sinh hóa mới
• Không dùng thay thế cho: Mô phỏng động học (ASM1, ASM2d, ADM1) khi thiết kế chi tiết hệ thống phức tạp

Trong thực tiễn tại Việt Nam, mô hình McCarty phù hợp nhất cho:

• Bước tính toán sơ bộ khi lập báo cáo ĐTM (đánh giá tác động môi trường)
• Giảng dạy tại các trường Đại học Bách Khoa, Đại học Xây dựng, Đại học Tài nguyên Môi trường
• Kiểm tra tính hợp lý của kết quả mô phỏng bằng phần mềm GPS-X, BioWin

8. Kết Luận

Mô hình McCarty R = fs·Rs + fa·Ra − Rd là công cụ lý thuyết đơn giản nhưng mạnh mẽ, cho phép kỹ sư môi trường xây dựng phương trình stoichiometry hoàn chỉnh cho bất kỳ quá trình sinh hóa nào trong xử lý nước thải:

• Với fs = 0,70 và glucose + O₂ + NH₄⁺, phương trình tổng hợp cho hiệu suất sinh khối 53% (Y = 0,49 g VSS/g COD) — phù hợp với thực nghiệm
• Bảng 5.5 cung cấp đầy đủ ΔG° cho 20 nửa phản ứng thành phần — tra cứu và tổ hợp linh hoạt cho mọi kịch bản xử lý
• Nguyên lý fs + fa = 1 và sự phân bổ điện tử giữa đồng hóa/dị hóa là cơ sở để hiểu và dự đoán lượng bùn, nhu cầu O₂, và tiềm năng năng lượng trong bất kỳ hệ xử lý nào

Đây là kiến thức sinh hóa môi trường nền tảng không thể thiếu cho kỹ sư xử lý nước thải chuyên nghiệp tại Việt Nam.