Tích Phân Phương Trình Động Học, Phản Ứng Bậc Không Và Bốn Mô Hình Kinh Nghiệm Trong Xử Lý Nước Thải Sinh Học

Từ Phương Trình Vi Phân Đến Công Cụ Thiết Kế Thực Tiễn

Các phương trình động học Monod và phân hủy nội sinh mô tả tốc độ tức thời (rate equations, vi phân). Trong thiết kế và vận hành thực tế, kỹ sư cần biết tổng lượng cơ chất loại bỏ sau một khoảng thời gian t — điều này đòi hỏi tích phân các phương trình vi phân.

GS. Lê Văn Cát (2007) trình bày quá trình tích phân phương trình tốc độ tiêu thụ cơ chất (6-14) và dẫn ra một kết quả quan trọng: trong điều kiện thực tế của nhiều bể xử lý sinh học, phương trình động học phức tạp dạng Monod có thể đơn giản hóa về phản ứng bậc không — dạng đơn giản nhất nhưng thực dụng nhất cho tính toán nhanh.

Bên cạnh phương trình Monod có cơ sở lý thuyết, nhiều phương trình kinh nghiệm khác (6-17 đến 6-20) được đề xuất từ quan sát thực tế. Mỗi phương trình phù hợp với một tình huống cụ thể và có giá trị thực dụng riêng, dù thiếu tính chặt chẽ lý thuyết như Monod.

Bài viết này phân tích sâu quá trình dẫn xuất từ vi phân đến tích phân, ý nghĩa của phản ứng bậc không, và hướng dẫn sử dụng bốn phương trình kinh nghiệm 6-17 đến 6-20 trong thiết kế và vận hành nhà máy XLNT tại Việt Nam.

Phần 1: Tích Phân Phương Trình Động Học Monod

1.1 Phương Trình Gốc Cần Tích Phân

Từ phương trình tốc độ tiêu thụ cơ chất 6-14:

vsu = dS/dt = −k × X × S / (Ks + S) (Phương trình 6-14)

Đây là phương trình vi phân thường (ODE) bậc nhất, có thể tách biến:

(Ks + S)/S × dS = −k × X × dt

Giả thiết quan trọng: X = hằng số trong suốt thời gian phản ứng. Giả thiết này đúng khi:

• Bể batch (closed system) trong thời gian ngắn (vài giờ đầu)
• X bổ sung (thêm bùn) hoặc phân hủy nội sinh bỏ qua so với tốc độ tiêu thụ cơ chất
• Hệ thống liên tục ở steady state (X không đổi theo thời gian)

1.2 Giải Tích Phân — Phương Trình 6-15

Tách biến và tích phân hai vế từ 0 đến t:

∫[(Ks + S)/S] dS = −k × X × ∫dt

Phân rã phân số: (Ks + S)/S = Ks/S + 1

∫[Ks/S + 1] dS = Ks × ln(S) + S + C

Áp dụng điều kiện ban đầu: t = 0 thì S = S₀; t = t thì S = S:

Ks × [ln(S₀) − ln(S)] + (S₀ − S) = k × X × t

Ks × ln(S₀/S) + S₀ − S = k × X × t (Phương trình 6-15)

Đây là nghiệm chính xác của phương trình Monod với X = hằng số. Phương trình này là phương trình ẩn (implicit equation) — không thể giải trực tiếp ra S theo t mà phải giải bằng phương pháp lặp (iterative) hoặc đồ thị.

1.3 Kiểm Tra Phương Trình 6-15 — Giới Hạn Thời Gian Dài

Khi t lớn (S → 0 hoàn toàn): Ks × ln(S₀/0) + S₀ − 0 = k × X × t → vế trái → ∞

Điều này cho thấy: cần t → ∞ để loại bỏ hoàn toàn cơ chất — phù hợp với thực tế (không bao giờ đạt S = 0).

Kiểm tra chiều: khi S₀ = 200 mg/L, S = 20 mg/L, Ks = 60 mg/L, k × X × t = k × X × t: Vế trái = 60 × ln(200/20) + 200 − 20 = 60 × ln(10) + 180 = 60 × 2,303 + 180 = 138,2 + 180 = 318,2 mg/L k × X × t = 318,2 → với k = 11,1 gCOD/gVSS·ngày; X = 3 gVSS/L → t = 318,2/(11,1 × 3) = 9,55 giờ (HRT ≈ 10 giờ).

Phần 2: Điều Kiện Ks << S₀ Và Phản Ứng Bậc Không

2.1 Bỏ Qua Số Hạng Ks × ln(S₀/S)

Trong phần lớn các trường hợp thực tế, nồng độ COD đầu vào S₀ >> Ks:

• Ks (BOD) ≈ 60 mg BOD/L
• S₀ điển hình = 200–1.000 mg COD/L
• Tỷ lệ S₀/Ks = 3–17 lần

Khi S₀ >> Ks → Ks × ln(S₀/S) << S₀ − S (trong giai đoạn đầu khi S vẫn còn lớn)

Bỏ qua số hạng Ks × ln(S₀/S) trong phương trình 6-15:

S = S₀ − k × X × t (Phương trình 6-16)

2.2 Ý Nghĩa Của Phản Ứng Bậc Không

Phương trình S = S₀ − k × X × t là phản ứng bậc không (zero-order reaction) — tốc độ giảm cơ chất không phụ thuộc vào nồng độ S mà chỉ phụ thuộc vào k × X:

|dS/dt| = k × X = hằng số

Đồ thị S theo t là đường thẳng dốc xuống với độ dốc = −k × X.

Điều kiện vật lý để phản ứng bậc không đúng: Đây đúng khi S >> Ks — tức là khi nồng độ cơ chất trong bể còn cao hơn nhiều so với Ks. Tại giai đoạn cuối khi S tiếp cận Ks, bậc không không còn đúng nữa và cần dùng phương trình 6-15 đầy đủ.

2.3 Ứng Dụng Thực Tiễn Của Phương Trình Bậc Không

Phương trình 6-16 được ứng dụng phổ biến vì tính đơn giản:

Ứng dụng 1 — Tính thời gian phản ứng: Biết S₀, mục tiêu S_e, k và X: t = (S₀ − S_e) / (k × X)

Ví dụ: S₀ = 300 mg COD/L; S_e = 30 mg COD/L; k = 11,1 gCOD/gVSS·ngày × (1 ngày/24 giờ) = 0,463 gCOD/gVSS·giờ; X = 3.000 mg VSS/L = 3,0 gVSS/L: t = (300 − 30) / (0,463 × 3,0) = 270 / 1,389 = 194 phút ≈ 3,25 giờ

→ HRT cần thiết ≈ 3,25 giờ (với giả thiết S > Ks trong suốt quá trình — cần kiểm tra: S_e = 30 mg < Ks = 60 mg → không thỏa mãn S >> Ks nữa → cần dùng phương trình 6-15 để chính xác hơn).

Ứng dụng 2 — Vẽ đường cong S theo t trong batch test: Trong batch test (thí nghiệm mẻ), vẽ COD theo thời gian trong giai đoạn đầu (khi S >> Ks) → đường thẳng → độ dốc = −k × X → tính k từ độ dốc và X đã biết.

Ứng dụng 3 — So sánh hoạt tính bùn từ các nguồn khác nhau: Cùng X và cùng thời gian → So sánh (S₀ − S_t): mẫu nào giảm nhiều hơn → hoạt tính (k) cao hơn.

2.4 Khi Nào Bậc Không Không Còn Đúng?

Phương trình 6-16 không còn đúng khi S tiếp cận Ks. Cụ thể, sai số < 10% khi S > 5×Ks:

• Với Ks = 60 mg BOD/L: áp dụng bậc không khi S > 300 mg BOD/L
• Với KN = 0,74 mg NH₄-N/L: áp dụng bậc không khi S > 3,7 mg NH₄-N/L (thực ra hầu hết các bể nitrification có NH₄ trong bể > 5 mg/L ở tải trung bình nên bậc không vẫn áp dụng được)

Phần 3: Bốn Phương Trình Kinh Nghiệm — Khi Nào Dùng Cái Nào?

3.1 Phương Trình 6-17 — Bậc Không Đơn Giản

vsu = −k (Phương trình 6-17)

Đây là dạng đơn giản nhất: tốc độ suy giảm cơ chất không phụ thuộc vào X, S, hay bất kỳ điều kiện nào khác. Tích phân trực tiếp:

S = S₀ − k × t

Khi nào dùng: Trường hợp quan sát thực tế cho thấy S giảm tuyến tính theo t với tốc độ gần như hằng số, bất kể thay đổi X (mật độ bùn) hay S. Điều này xảy ra khi:

• Quá trình bị giới hạn bởi yếu tố vật lý (ví dụ tốc độ khuếch tán, hòa tan khí) chứ không phải yếu tố sinh hóa
• Hệ thống đã đạt điều kiện bão hòa (X và S đều >> Ks) → tốc độ phụ thuộc vào yếu tố khác
• Mô hình hóa đơn giản để ước tính nhanh tải trọng xử lý

Ví dụ ứng dụng: Quá trình loại bỏ mùi (odor removal) trong bể biofilter khi nồng độ H₂S thấp nhưng ổn định — tốc độ loại bỏ ổn định không phụ thuộc lớn vào nồng độ H₂S.

3.2 Phương Trình 6-18 — Bậc Một Theo S

vsu = −k × S (Phương trình 6-18)

Tốc độ tiêu thụ tỉ lệ với nồng độ cơ chất, không phụ thuộc X. Nghiệm tích phân:

S = S₀ × e^(−k × t)

Đây là phản ứng bậc một cổ điển — đường cong suy giảm dạng exponential.

Khi nào dùng: Phương trình 6-18 phù hợp khi:

• Nồng độ cơ chất thấp (S << Ks) → Monod đơn giản về bậc một: vsu ≈ −(k/Ks) × S = −k’ × S
• Mô hình hóa quá trình loại bỏ chất ô nhiễm vi lượng (micropollutants) như thuốc trừ sâu, hormone, kháng sinh tồn dư — thường ở nồng độ µg/L, << Ks
• Dùng phổ biến trong kiểm soát dịch tễ: mô hình tử vong vi sinh vật gây bệnh (pathogen inactivation)

Ví dụ tính toán: Loại bỏ kháng sinh ciprofloxacin (S₀ = 100 µg/L) trong bùn hoạt tính với k’ = 0,5 L/gVSS·ngày; X = 3.000 mgVSS/L: k_eff = k’ × X = 0,5 × 3 = 1,5 ngày⁻¹ S(t) = 100 × e^(−1,5t) → sau HRT = 1 ngày: S = 100 × e^(−1,5) = 100 × 0,223 = 22 µg/L (loại bỏ 78%).

3.3 Phương Trình 6-19 — Bậc Hai (Tỉ Lệ X và S)

vsu = −k × X × S (Phương trình 6-19)

Đây là dạng đơn giản hóa của Monod khi S << Ks: Monod → vsu ≈ −(k/Ks) × X × S → đổi hằng số thành k × X × S.

Phương trình 6-19 khác phương trình Monod đầy đủ ở chỗ không có hệ số (Ks + S) ở mẫu. Đây là xấp xỉ bậc một theo S và bậc một theo X.

Nghiệm tích phân (với X = const): S = S₀ × e^(−k × X × t)

Khi nào dùng:

• S << Ks (nồng độ cơ chất thấp hơn Ks nhiều)
• Đặc biệt phù hợp cho nitrification giai đoạn cuối (S = NH₄ gần tiêu chuẩn đầu ra, KN = 0,74 mg/L)
• Mô hình hóa hệ thống bể nối tiếp (tanks-in-series): mỗi bể ứng với S khác nhau, nhưng khi S << Ks thì tất cả đều áp dụng bậc một theo S

Ví dụ: Bể nitrification SBR, pha react cuối (NH₄ = 2 mg/L, KN = 0,74 → S/KN = 2,7 → không nhỏ hơn Ks nhiều, nhưng gần đủ để dùng 6-19 như xấp xỉ).

3.4 Phương Trình 6-20 — Bậc Hai Có Chuẩn Hóa S₀

vsu = −k × X × S / S₀ (Phương trình 6-20)

Phương trình này có thêm S₀ ở mẫu — điều chỉnh tốc độ theo mức độ “tương đối” của cơ chất so với giá trị ban đầu.

Đây là phương trình kinh nghiệm thuần túy, không có cơ sở lý thuyết rõ ràng từ Monod. Tuy nhiên nó phản ánh một quan sát thực tế: hệ vi sinh thích nghi với tải COD cao hơn (S₀ lớn) thường loại bỏ cơ chất nhanh hơn — bùn đã “học” cách xử lý nước thải nồng độ đó.

Nghiệm tích phân: S = S₀ × e^(−k × X × t / S₀) = S₀ × e^(−k’ × t) với k’ = k × X / S₀

Khi nào dùng:

• Nước thải có S₀ biến động lớn (ví dụ nước thải làng nghề theo mùa vụ)
• Mô hình hóa hệ thống batch có tải trọng đầu mẻ khác nhau
• Nghiên cứu sự thích nghi bùn với nồng độ nước thải cụ thể

Phần 4: So Sánh Bốn Mô Hình — Lựa Chọn Phù Hợp Với Điều Kiện Thực Tế

4.1 Bảng Tóm Tắt Bốn Phương Trình Kinh Nghiệm

Phương Trình	Dạng	Nghiệm (S theo t)	Điều Kiện Áp Dụng	Ứng Dụng Điển Hình
6-17: vsu = −k	Bậc 0	S = S₀ − k × t	Không phụ thuộc X, S	Bão hòa, giới hạn vật lý
6-18: vsu = −k × S	Bậc 1 theo S	S = S₀ × e^(−kt)	S << Ks, X không đổi	Micropollutants, inactivation
6-19: vsu = −k × X × S	Bậc 2	S = S₀ × e^(−kXt)	S << Ks, X quan trọng	Nitrification cuối bể
6-20: vsu = −k × X × S/S₀	Biến đổi	S = S₀ × e^(−kXt/S₀)	Thích nghi tải, S₀ biến động	Nước thải mùa vụ

4.2 Phương Pháp Chọn Mô Hình Từ Dữ Liệu Thực Nghiệm

Khi có dữ liệu S theo t từ batch test, chọn mô hình phù hợp bằng phân tích hồi quy:

Bước 1: Vẽ đồ thị S theo t → đường thẳng → bậc không (6-17 hoặc đặc biệt hóa 6-16)

Bước 2: Nếu không thẳng, vẽ ln(S) theo t → đường thẳng → bậc một (6-18 hoặc 6-19 nếu kiểm tra vai trò X)

Bước 3: Thay đổi X trong batch test (X₁, X₂, X₃) → nếu k_eff ∝ X → phương trình 6-19; nếu k_eff không phụ thuộc X → phương trình 6-18

Bước 4: Thay đổi S₀ → nếu k_eff ∝ 1/S₀ → phương trình 6-20

4.3 Hạn Chế Của Tất Cả Mô Hình — Nhận Định Của GS. Lê Văn Cát

GS. Lê Văn Cát (2007) nhận định thẳng thắn: “Tính chặt chẽ của các phương trình động học và lời giải của các hệ xử lý vi sinh không cao vì trong hệ xảy ra đồng thời nhiều quá trình chuyển hóa được thực hiện bởi nhiều chủng loại vi sinh với tốc độ rất khác nhau.”

Đây là thực tế quan trọng cần thẳng thắn thừa nhận:

• Bùn hoạt tính gồm hàng nghìn chủng vi sinh với k và Ks hoàn toàn khác nhau
• “k” đo được là hằng số đặc trưng trung bình của cả quần thể — không đặc trưng riêng cho bất kỳ chủng nào
• Các thông số k, Ks đo tại phòng thí nghiệm với nước thải nhân tạo có thể sai lệch đáng kể so với điều kiện thực (nước thải thực có hàng trăm thành phần phức tạp hơn)

Vì vậy, nguyên tắc sử dụng mô hình động học trong thiết kế thực tế tại Việt Nam:

• Dùng mô hình để hướng thiết kế, không phải dự đoán chính xác tuyệt đối
• Áp dụng hệ số an toàn đủ lớn (SF = 3–10 tùy loại vi sinh và mức độ tin cậy của tham số)
• Calibrate mô hình với dữ liệu thực từ nhà máy tương tự (cùng loại nước thải, cùng khí hậu)
• Kết hợp mô hình với kinh nghiệm vận hành thực tế — mô hình là công cụ hỗ trợ, không thay thế được phán đoán kỹ thuật

Phần 5: Ứng Dụng Tổng Hợp — Từ Mô Hình Đến Thiết Kế Bể Phản Ứng Dạng Mẻ

5.1 Thiết Kế Bể SBR Dựa Trên Phương Trình 6-15 Và 6-16

Bể SBR (Sequencing Batch Reactor) vận hành theo mẻ — thích hợp để áp dụng phương trình tích phân:

Yêu cầu: Xử lý 1.000 m³ nước thải/ngày, COD = 500 mg/L, yêu cầu đầu ra COD < 50 mg/L.

Dữ liệu: k = 11,1 gCOD/gVSS·ngày = 0,463 gCOD/gVSS·giờ; X = 3.500 mg VSS/L = 3,5 gVSS/L; Ks = 60 mg/L

Bước 1: Kiểm tra điều kiện bậc không (S₀/Ks = 500/60 = 8,3 → có thể dùng 6-16 trong giai đoạn đầu khi S > 150 mg/L)

Bước 2: Tính t từ S = 500 → 150 mg/L (giai đoạn bậc không): t₁ = (500 − 150) / (0,463 × 3,5) = 350/1,620 = 216 phút

Bước 3: Tính t từ S = 150 → 50 mg/L bằng phương trình 6-15 (vì S ≈ Ks): Ks × ln(150/50) + 150 − 50 = k × X × t₂ 60 × 1,099 + 100 = 1,620 × t₂ (giờ → đổi sang phút: 1,620 × 60 = 97,2 mg/L·phút) 65,9 + 100 = 97,2 × t₂ (phút) → t₂ = 165,9/97,2 = 171 phút

Tổng thời gian phản ứng: t_react = 216 + 171 = 387 phút ≈ 6,5 giờ

Cộng thêm thời gian Fill (30 phút) + Settle (60 phút) + Decant (30 phút) = tổng chu kỳ = 507 phút ≈ 8,5 giờ

Số chu kỳ/ngày = 24 × 60 / 507 = 2,8 chu kỳ/ngày

Thể tích bể = 1.000 / 2,8 = 357 m³ (cần 2 bể song song, mỗi bể 357 m³ vận hành xen kẽ).

Kết Luận

Hành trình từ phương trình vi phân Monod đến phương trình tích phân thực dụng phản ánh một quá trình quan trọng trong kỹ thuật: chuyển đổi mô hình lý thuyết phức tạp thành công cụ tính toán đơn giản, đủ chính xác cho mục đích thiết kế.

Ba điểm chính:

Một, phương trình 6-15 là nghiệm chính xác (với X = hằng số). Khi S₀ >> Ks, đơn giản hóa về phương trình 6-16 (bậc không) — đủ chính xác cho giai đoạn đầu khi COD còn cao. Giai đoạn cuối (S gần Ks) cần dùng 6-15 đầy đủ.

Hai, bốn phương trình kinh nghiệm 6-17 đến 6-20 phản ánh sự đa dạng trong hành vi thực tế. Mỗi dạng phù hợp với một điều kiện cụ thể: bậc không (bão hòa), bậc một theo S (S << Ks), bậc hai (S << Ks, X quan trọng), và phương trình chuẩn hóa S₀ (tải biến động).

Ba, tất cả mô hình đều có hạn chế vì bùn hoạt tính là hỗn hợp phức tạp của hàng nghìn chủng vi sinh. Nguyên tắc thực tiễn: dùng mô hình để hướng thiết kế, áp dụng hệ số an toàn đủ lớn, và luôn kiểm tra với dữ liệu thực tế từ công trình tương tự.