Từ Gián Đoạn Sang Liên Tục — Sự Khác Biệt Cơ Bản

Trong kỹ thuật xử lý nước thải, bể phản ứng gián đoạn (batch reactor) và bể phản ứng liên tục (continuous reactor) là hai phương thức vận hành căn bản. Bể gián đoạn (SBR) đã được phân tích với phương trình PT 7-1 đến PT 7-7, trong đó nồng độ chỉ thay đổi theo thời gian trong một thể tích kín không có dòng chảy vào ra trong quá trình phản ứng.

Khi chuyển sang bể phản ứng liên tục khuấy trộn đều (CSTR — Continuous Stirred Tank Reactor), bức tranh thay đổi hoàn toàn: nước thải liên tục chảy vào với lưu lượng Q (m³/h) mang theo nồng độ C₀, đồng thời nước đã xử lý liên tục chảy ra mang theo nồng độ C thấp hơn. Tính liên tục này đưa vào phương trình một số hạng bổ sung, làm cho mối quan hệ C = f(t) phức tạp hơn nhiều so với trường hợp gián đoạn — như GS. Lê Văn Cát đã nhấn mạnh trong công trình "Xử Lý Nước Thải Giàu Hợp Chất Nitơ và Phốtpho" (NXB Khoa học Tự nhiên và Công nghệ, 2007).

Bài viết này phân tích chi tiết quá trình thiết lập phương trình thăng bằng chất PT 7-8, tổng hợp thành dạng chuẩn PT 7-9 với hệ số α, và dẫn xuất nghiệm quá độ PT 7-10 mô tả sự biến thiên nồng độ C(t) trong giai đoạn hệ thống chưa đạt trạng thái ổn định.
 

Phần 1: Nguyên Lý Thăng Bằng Chất Trong CSTR — Cơ Sở Thiết Lập PT 7-8

1.1 Hệ Thống Vật Lý Của CSTR

Xét bể phản ứng liên tục khuấy trộn đều với các thông số:

• Thể tích bể: V (m³) — không đổi
• Lưu lượng dòng vào = dòng ra: Q (m³/h) — điều kiện cân bằng thể tích
• Nồng độ đầu vào: C₀ (g/m³ hay mg/L) — không đổi (giả thiết)
• Nồng độ trong bể = nồng độ đầu ra: C (g/m³) — thay đổi theo thời gian t
• Tốc độ phản ứng tiêu thụ cơ chất: vi (g/m³·h) — giá trị dương theo quy ước tiêu thụ

Điều kiện khuấy trộn đều hoàn toàn (perfect mixing) đảm bảo: nồng độ tại mọi điểm trong bể tại mỗi thời điểm đều bằng nhau và bằng nồng độ đầu ra. Đây là điều kiện lý tưởng hóa CSTR.

1.2 Phương Trình Thăng Bằng Chất Tổng Quát — PT 7-8

Áp dụng nguyên lý bảo toàn khối lượng cho cấu tử nguyên liệu trong bể trong khoảng thời gian dt:

Biến động lượng chất trong bình = Tăng do đầu vào − Giảm do đầu ra − Giảm do phản ứng tiêu thụ

Lượng chất trong bể tại thời điểm t: V·C (g) Biến động trong dt: d(V·C)/dt = V·dC/dt (vì V = const)

Vế phải:

• Vào: Q·C₀ (g/h)
• Ra: Q·C (g/h)
• Tiêu thụ bởi phản ứng: vi·V (g/h)

Phương trình thăng bằng chất:

V·(dC/dt) = Q·C₀ − Q·C − vi·V (PT 7-8)

Hay viết gọn:

dC/dt = (Q/V)·C₀ − (Q/V)·C − vi

1.3 Ý Nghĩa Từng Số Hạng Trong PT 7-8

Số hạng Q·C₀/V: Tốc độ nạp cơ chất vào bể trên đơn vị thể tích (g/m³·h). Đây là "nguồn" liên tục bổ sung cơ chất.

Số hạng −Q·C/V: Tốc độ thải cơ chất qua dòng ra trên đơn vị thể tích. Đây là "bể chứa" hút cơ chất ra khỏi bể, tốc độ tỷ lệ với nồng độ C hiện tại.

Số hạng −vi: Tốc độ tiêu thụ cơ chất bởi phản ứng hóa học hoặc sinh học. Trong xử lý nước thải, đây là tốc độ phân hủy hữu cơ hay nitrat hóa.

So sánh với bể gián đoạn PT 7-1 (-dC/dt = vi): bể gián đoạn chỉ có số hạng phản ứng, không có số hạng đầu vào/ra. Đây là sự khác biệt cơ bản.

Phần 2: Đơn Giản Hóa Với Bậc Nhất Và Tham Số α — PT 7-9

2.1 Giả Thiết Phản Ứng Bậc Nhất

Giả thiết phản ứng tiêu thụ cơ chất tuân theo bậc nhất: vi = k·C (trong đó k > 0 là hằng số tốc độ phản ứng, giờ⁻¹).

Thay vi = k·C vào PT 7-8:

V·dC/dt = Q·C₀ − Q·C − k·C·V

Chia hai vế cho V:

dC/dt = (Q/V)·C₀ − (Q/V)·C − k·C

dC/dt = (Q/V)·C₀ − [(Q/V) + k]·C

Đặt α = k + Q/V, PT 7-8 trở thành:

dC/dt = −α·C + (Q/V)·C₀ (PT 7-9)

Đây là phương trình vi phân tuyến tính bậc nhất với hệ số hằng, dạng chuẩn: dC/dt + α·C = (Q/V)·C₀

2.2 Ý Nghĩa Vật Lý Của Tham Số α = k + Q/V

Tham số α có đơn vị h⁻¹ và là hằng số tốc độ tổng hợp của quá trình phản ứng liên tục:

α = k + Q/V = k + 1/τ

Trong đó τ = V/Q là thời gian lưu thủy lực (HRT), và 1/τ = Q/V là nghịch đảo HRT — tốc độ "pha loãng" hay "thải" cơ chất qua dòng ra.

Diễn giải: Trong bể liên tục, cơ chất C bị loại bỏ bởi hai cơ chế đồng thời:

1. Phản ứng hóa/sinh học với tốc độ k·C
2. Dòng ra mang cơ chất ra ngoài với tốc độ (Q/V)·C = C/τ

Tổng tốc độ loại bỏ = (k + 1/τ)·C = α·C

Đây là điểm khác biệt căn bản so với bể gián đoạn: trong bể gián đoạn, cơ chất chỉ bị loại bỏ bởi phản ứng (tốc độ k·C). Trong CSTR, còn có thêm "hao hụt" do dòng ra với tốc độ C/τ.

2.3 Tầm Quan Trọng Của Tỷ Lệ k/τ⁻¹

So sánh hai thành phần của α:

Trường hợp k >> 1/τ (phản ứng nhanh, τ lớn): α ≈ k, hệ thống hành xử gần như bể gián đoạn — thời gian lưu nước đủ dài để phản ứng hoàn thành.

Trường hợp k << 1/τ (phản ứng chậm, τ nhỏ): α ≈ 1/τ, hệ thống bị khống chế bởi thời gian lưu nước — cơ chất bị "đẩy" ra khỏi bể trước khi phản ứng kịp hoàn thành. Đây là tình trạng "washout" (rửa trôi) nguy hiểm trong thực tế vận hành.

Trường hợp k ≈ 1/τ: Cả hai cơ chế đóng góp tương đương — thiết kế cần cân bằng cẩn thận.

Ví dụ thực tế: Với k_nitrat hóa = 0,2 giờ⁻¹ (ở 20°C) và τ = 8 giờ (HRT điển hình bể aerotank):

• k = 0,2 giờ⁻¹
• 1/τ = 0,125 giờ⁻¹
• α = 0,2 + 0,125 = 0,325 giờ⁻¹

Hai thành phần đóng góp gần bằng nhau — không thể bỏ qua thành phần dòng chảy.

Phần 3: Dẫn Xuất Nghiệm Quá Độ PT 7-10

3.1 Phương Pháp Giải Phương Trình Vi Phân PT 7-9

Phương trình PT 7-9 là phương trình vi phân tuyến tính bậc nhất phi thuần nhất:

dC/dt + α·C = (Q/V)·C₀

Giải bằng phương pháp nghiệm tổng quát (variation of parameters hoặc integrating factor):

Bước 1 — Nghiệm thuần nhất (vế phải = 0): dC_h/dt + α·C_h = 0 → C_h = A·e^(-α·t)

Bước 2 — Nghiệm riêng (hằng số): Thử C_p = const: dC_p/dt = 0 → α·C_p = (Q/V)·C₀ → C_p = (Q/V)·C₀/α = Q·C₀/(V·α)

Bước 3 — Nghiệm tổng quát: C(t) = C_h + C_p = A·e^(-α·t) + Q·C₀/(V·α)

Bước 4 — Áp dụng điều kiện đầu: tại t = 0, C = C_in (nồng độ ban đầu trong bể khi bắt đầu vận hành liên tục): C_in = A + Q·C₀/(V·α) → A = C_in - Q·C₀/(V·α)

Nghiệm tổng quát:

C(t) = Q·C₀/(V·α) · (1 - e^(-α·t)) + C_in · e^(-α·t) (PT 7-10)

Trong trường hợp đặc biệt quan trọng: khi bắt đầu vận hành từ bể sạch (C_in = C₀, tức là ban đầu bể chứa nước thải cùng nồng độ đầu vào), ta có C_in = C₀:

C(t) = (Q·C₀)/(V·α) · (1 - e^(-α·t)) + C₀ · e^(-α·t) (PT 7-10, dạng thực hành)

3.2 Phân Tích Cấu Trúc Nghiệm PT 7-10

Nghiệm PT 7-10 gồm hai số hạng với ý nghĩa vật lý rõ ràng:

Số hạng thứ nhất: Q·C₀/(V·α) · (1 - e^(-α·t))

• Tại t = 0: số hạng này = 0
• Khi t → ∞: số hạng này → Q·C₀/(V·α) = C_∞ (nồng độ trạng thái ổn định)
• Đây là phần "tăng dần" tiến đến trạng thái ổn định

Số hạng thứ hai: C₀ · e^(-α·t)

• Tại t = 0: số hạng này = C₀
• Khi t → ∞: số hạng này → 0
• Đây là phần "suy giảm" của điều kiện ban đầu

Kiểm tra tính đúng đắn:

• t = 0: C(0) = 0 + C₀ = C₀ ✓ (đúng điều kiện đầu)
• t → ∞: C(∞) = Q·C₀/(V·α) + 0 = C_∞ ✓ (tiệm cận về trạng thái ổn định)

3.3 So Sánh PT 7-10 (CSTR) Với PT 7-3 (Batch)

GS. Lê Văn Cát nhấn mạnh: so với PT 7-3 (batch: C = C₀·e^(-kt)), PT 7-10 "phức tạp hơn nhiều, gây ra bởi tính liên tục của quá trình phản ứng." Sự phức tạp thể hiện ở:

Đặc điểm	Batch PT 7-3	CSTR PT 7-10
Số số hạng	1 (thuần suy giảm)	2 (suy giảm + tăng dần)
Hằng số trong hàm mũ	k (phản ứng đơn thuần)	α = k + 1/τ (tổng hợp)
C khi t→∞	0 (tiêu thụ hết)	Q·C₀/(V·α) > 0 (trạng thái ổn định ≠ 0)
Ảnh hưởng C₀	Chỉ qua điều kiện đầu	Qua cả điều kiện đầu và trạng thái ổn định
Ảnh hưởng Q	Không có (Q = 0 trong batch)	Qua α và Q/V·C₀

Trong batch, nồng độ tiệm cận về 0 khi t → ∞ (lý thuyết). Trong CSTR, nồng độ tiệm cận về một giá trị hữu hạn C_∞ > 0 — đây là trạng thái ổn định, không thể loại hết cơ chất trong bể liên tục.

Phần 4: Dạng Viết Lại PT 7-11 Với Thời Gian Lưu Thủy Lực τ

4.1 Thay Thế Q/V = 1/τ Vào PT 7-10

Sử dụng định nghĩa τ = V/Q (thời gian lưu thủy lực, HRT), ta có Q/V = 1/τ và:

α = k + Q/V = k + 1/τ

Thay vào PT 7-10:

C(t) = (1/τ)·C₀ / (k + 1/τ) · (1 - e^(-α·t)) + C₀·e^(-α·t)

Đây chính là PT 7-11:

C(t) = C₀ · [1/(k·τ + 1)] · (1 - e^(-(1/τ+k)·t)) + C₀ · e^(-(1/τ+k)·t) (PT 7-11)

4.2 Ý Nghĩa Của Tích k·τ

Tích k·τ = k·(V/Q) là thông số vô thứ nguyên quan trọng (Damköhler number Da):

Da = k·τ

Da đặc trưng cho tỷ lệ giữa tốc độ phản ứng và tốc độ dòng chảy:

• Da >> 1: phản ứng nhanh so với dòng chảy → hiệu suất cao
• Da << 1: phản ứng chậm so với dòng chảy → hiệu suất thấp
• Da = 1: điểm cân bằng, cần τ = 1/k để đạt E = 50% ở trạng thái ổn định

Ví dụ: Bể hiếu khí xử lý BOD với k = 0,2 giờ⁻¹ và τ = 8 giờ → Da = 0,2 × 8 = 1,6

Trạng thái ổn định: C_∞ = C₀/(1 + Da) = C₀/2,6 ≈ 0,385·C₀

→ E_∞ = 1 - 1/2,6 ≈ 61,5% — cần tăng τ hoặc k để đạt tiêu chuẩn xả thải cao hơn.

Phần 5: Đặc Điểm Quá Độ Của CSTR — Ứng Dụng Trong Khởi Động Hệ Thống

5.1 Tốc Độ Đạt Trạng Thái Ổn Định

Từ PT 7-10, nồng độ đạt "đủ gần" trạng thái ổn định (98%) khi:

e^(-α·t) ≈ 0,02 → α·t ≈ 4 → t_98% ≈ 4/α = 4/(k + 1/τ)

Đây là thời gian khởi động thực tế của CSTR.

Ví dụ: với α = 0,325 giờ⁻¹ (như ví dụ ở trên): t_98% = 4/0,325 ≈ 12,3 giờ

Hệ thống cần khoảng 12-13 giờ vận hành liên tục để đạt trạng thái ổn định từ điều kiện ban đầu C = C₀. Đây là thông tin quan trọng cho kỹ thuật viên vận hành trạm XLNT khi khởi động hệ thống sau bảo trì.

5.2 Ảnh Hưởng Của Tải Trọng Đột Biến

Khi nồng độ đầu vào thay đổi đột ngột từ C₀ lên C₀' (shock load), PT 7-10 mô tả quá trình điều chỉnh của bể:

C(t) = C_∞' · (1 - e^(-α·t)) + C_∞_cũ · e^(-α·t)

Trong đó C_∞' = Q·C₀'/(V·α) là trạng thái ổn định mới.

Khả năng "giảm chấn" của CSTR: Nhờ thể tích V lớn, bể CSTR đệm tải trọng đột biến hiệu quả. Nồng độ trong bể không tăng ngay lên C₀' mà tăng từ từ theo hàm mũ với hằng số thời gian 1/α. Đây là ưu điểm nổi bật của CSTR so với PFR (dòng đẩy) trong xử lý nước thải công nghiệp có tải trọng dao động mạnh.

Thực tế tại Việt Nam: Các KCN chế biến thực phẩm ở Đồng Nai, Bình Dương thường có tải trọng COD dao động mạnh theo ca sản xuất (ví dụ: buổi sáng 500 mg/L, buổi trưa 2.000 mg/L khi xả nước vệ sinh). Bể CSTR với τ = 12-24 giờ sẽ "trung bình hóa" tải trọng, bảo vệ quá trình vi sinh khỏi sốc tải trọng.

5.3 Ứng Dụng PT 7-10 Trong Tối Ưu Hóa Quá Trình Khởi Động

Trong giai đoạn khởi động (start-up) hệ thống XLNT mới tại Việt Nam, kỹ sư thường muốn rút ngắn thời gian đạt trạng thái ổn định mà vẫn đảm bảo chất lượng đầu ra. PT 7-10 gợi ý các chiến lược:

Chiến lược 1: Tăng α bằng cách giảm τ ban đầu Khởi động với Q cao hơn thiết kế (τ nhỏ hơn) → α lớn hơn → đạt trạng thái ổn định nhanh hơn. Sau đó điều chỉnh Q về giá trị thiết kế. Tuy nhiên, cần đảm bảo vi sinh vật đủ thời gian thích nghi.

Chiến lược 2: Cấy bùn giống (seed sludge) Điền bể với bùn hoạt tính từ hệ thống đang vận hành → C_in thấp hơn C₀ ngay từ đầu → hệ thống đạt E mục tiêu sớm hơn. Trong PT 7-10, C_in nhỏ làm giảm số hạng thứ hai, rút ngắn thời gian quá độ.

Chiến lược 3: Chạy tuần hoàn nội bộ Trong giai đoạn đầu, giảm Q đến mức tối thiểu → τ lớn → α nhỏ hơn nhưng C_∞ = Q·C₀/(V·α) thấp hơn. Kỹ sư chọn Q phù hợp để đạt trạng thái ổn định với C_∞ đạt tiêu chuẩn ngay khi hệ thống ổn định lần đầu.

Phần 6: Ứng Dụng Thực Tế PT 7-8 Đến PT 7-10 Tại Các Công Trình Việt Nam

6.1 Bể Aerotank Thông Thường — CSTR Điển Hình

Bể aerotank trong hệ thống bùn hoạt tính (activated sludge) là ứng dụng CSTR phổ biến nhất:

Thông số điển hình (XLNT sinh hoạt, TP.HCM):

• V = 2.000 m³
• Q = 6.000 m³/ngày = 250 m³/giờ
• τ = V/Q = 2.000/250 = 8 giờ
• C₀_BOD = 200 mg/L (đầu vào)
• k_BOD = 0,18 giờ⁻¹ (20°C, nước thải sinh hoạt)
• α = 0,18 + 1/8 = 0,18 + 0,125 = 0,305 giờ⁻¹

Thời gian đạt trạng thái ổn định: t_98% = 4/0,305 ≈ 13,1 giờ

Sau khởi động hoặc sau sự cố ngắn, bể cần khoảng 13 giờ để trở lại trạng thái vận hành bình thường.

Nồng độ BOD trạng thái ổn định: C_∞ = C₀/(1 + k·τ) = 200/(1 + 0,18×8) = 200/2,44 ≈ 82 mg/L

→ E_∞ = 1 - 82/200 ≈ 59% — không đạt tiêu chuẩn QCVN 14 (BOD ≤ 30 mg/L). Cần tăng τ hoặc tăng k (bằng cách duy trì SRT cao hơn, nhiệt độ phù hợp).

6.2 Bể Kỵ Khí UASB — CSTR Trong Điều Kiện Yếm Khí

Bể UASB (Upflow Anaerobic Sludge Blanket) tại các nhà máy bia, chế biến tinh bột cassava ở Bình Dương, Tây Ninh về cơ bản là CSTR với dòng đi lên từ dưới:

Thông số điển hình (nước thải sản xuất bia):

• COD đầu vào C₀ = 3.000 mg/L
• Nhiệt độ 35°C (điều chỉnh)
• k_COD_yếm khí = 0,05 giờ⁻¹ ở 35°C
• τ thiết kế = 24 giờ (HRT = 1 ngày)
• α = 0,05 + 1/24 = 0,05 + 0,042 = 0,092 giờ⁻¹

Nồng độ COD trạng thái ổn định: C_∞ = 3.000/(1 + 0,05×24) = 3.000/2,2 ≈ 1.364 mg/L

E_∞ = 1 - 1.364/3.000 ≈ 54,5%

→ Sau xử lý UASB, nước thải vẫn cần bể hiếu khí tiếp theo để đạt tiêu chuẩn. Đây là lý do thiết kế hệ thống 2 bậc (kỵ khí + hiếu khí) phổ biến trong xử lý nước thải công nghiệp thực phẩm.

6.3 Điều Chỉnh PT 7-10 Cho Nước Thải Có Tải Trọng Dao Động

Trong thực tế, C₀ không phải là hằng số — nó biến thiên theo giờ trong ngày. PT 7-8 vẫn đúng nhưng cần giải số với C₀ = C₀(t) biến đổi:

V·dC/dt = Q·C₀(t) − Q·C − k·C·V

Giải số bằng phương pháp Euler hoặc Runge-Kutta cho phép mô phỏng diễn biến nồng độ C(t) theo tải trọng thực tế. Đây là công cụ phân tích quan trọng trong thiết kế hệ thống cân bằng (equalization tank) trước bể CSTR — giúp giảm dao động C₀ và ổn định quá trình xử lý.

Kết Luận

Phương trình thăng bằng chất PT 7-8 là khởi điểm lý thuyết cho toàn bộ phân tích CSTR, phản ánh ba cơ chế đồng thời: nạp cơ chất từ dòng vào, thải cơ chất qua dòng ra, và tiêu thụ bởi phản ứng. Sự tổng hợp này tạo ra tham số α = k + 1/τ — hằng số tốc độ hiệu dụng đặc trưng cho hệ thống liên tục.

Nghiệm quá độ PT 7-10 mô tả đầy đủ diễn biến nồng độ C(t) từ lúc khởi động đến khi đạt trạng thái ổn định, với hằng số thời gian 1/α. Sự phức tạp của PT 7-10 so với PT 7-3 (batch) phản ánh bản chất tích hợp của CSTR — vừa là bể phản ứng vừa là thiết bị phân tách dòng chảy. Nắm vững đặc tính quá độ này là nền tảng để kỹ sư vận hành dự đoán hành vi hệ thống khi khởi động, khi gặp sốc tải trọng, và khi điều chỉnh lưu lượng tại các trạm XLNT thực tế tại Việt Nam.