Giới Thiệu: Tại Sao Cần Nhiều Bể CSTR Nối Tiếp?

Trong thiết kế hệ thống xử lý nước thải (XLNT), một trong những câu hỏi thực tiễn quan trọng nhất là: làm thế nào để tăng hiệu suất xử lý mà không tốn quá nhiều thể tích bể? Câu trả lời không phải lúc nào cũng là xây một bể CSTR (Continuous Stirred Tank Reactor) thật lớn — vì như đã thấy qua PT 7-14, hiệu suất E = k·τ/(1 + k·τ) tăng rất chậm theo τ khi Da = k·τ đã lớn.

Giải pháp hiệu quả hơn nhiều: ghép nối tiếp nhiều bể CSTR để tạo ra chuỗi phản ứng bậc thang. Đầu ra của bể này là đầu vào của bể tiếp theo, nồng độ giảm dần qua mỗi bể. GS. Lê Văn Cát trong "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) đã xây dựng hệ phương trình mô tả chuỗi n bể CSTR, bắt đầu từ phương trình thăng bằng chất PT 7-22 cho bể thứ 2 trong chuỗi.

Bài viết này tập trung vào quá trình thiết lập phương trình, định nghĩa các tham số đặc trưng τ và H, và dẫn xuất PT 7-25 — phương trình vi phân nền tảng để tính nồng độ tại bể số 2.
 

Phần 1: Cấu Trúc Hệ Thống n Bể CSTR Nối Tiếp

1.1 Mô Tả Hệ Thống

Xét hệ thống gồm n bể CSTR nối tiếp với:

• Tổng thể tích: V (m³) — không đổi
• Thể tích mỗi bể: V/n (m³) — tất cả bằng nhau
• Lưu lượng dòng chảy: Q (m³/h) — qua toàn bộ chuỗi
• Nồng độ đầu vào bể 1: C₀ (mg/L)
• Nồng độ đầu ra bể i: Ci (mg/L)

Dòng chảy qua mỗi bể:

• Bể 1: vào C₀, ra C1
• Bể 2: vào C1, ra C2
• Bể i: vào Ci-1, ra Ci
• Bể n: vào Cn-1, ra Cn (= nồng độ đầu ra toàn hệ)

1.2 Đặc Tính Thủy Lực Mỗi Bể

Thời gian lưu thủy lực của mỗi bể riêng lẻ:

τ_mỗi bể = (V/n)/Q = V/(n·Q) = τ/n

Trong đó τ = V/Q là thời gian lưu thủy lực toàn hệ.

Lưu ý quan trọng: khi ghép nối tiếp n bể, mỗi bể có τ_riêng = τ/n, nhưng tổng τ của toàn chuỗi vẫn là τ = V/Q — giống như một bể CSTR đơn có cùng tổng thể tích. Điều này có nghĩa là việc chia nhỏ thành n bể không làm thay đổi tổng thời gian lưu nước, chỉ thay đổi phân bố nồng độ.

1.3 Tại Sao Bể 2 Là Bắt Đầu Phức Tạp?

Bể 1 trong chuỗi hoạt động như một CSTR thông thường — đã được phân tích đầy đủ bởi PT 7-8 đến PT 7-14. Bể 2 là nơi bắt đầu sự phức tạp: đầu vào của nó là C1 — một hàm thay đổi theo thời gian (kết quả từ bể 1), không phải hằng số C₀ như đầu vào của bể 1.

Đây là lý do phương trình thăng bằng chất cho bể 2 (PT 7-22) khó hơn PT 7-8 cho bể 1.

Phần 2: Phương Trình Thăng Bằng Chất PT 7-22 Cho Bể Thứ 2

2.1 Thiết Lập PT 7-22

Áp dụng nguyên lý bảo toàn khối lượng cho bể 2 (thể tích V/n):

Tích lũy = Vào − Ra − Tiêu thụ bởi phản ứng

(V/n)·dC2/dt = Q·C1 − Q·C2 − k·C2·(V/n)

(V/n)·dC2/dt = Q·C1 − Q·C2 (PT 7-22, bỏ qua phản ứng ban đầu để đơn giản)

Lưu ý: PT 7-22 như được viết chỉ gồm số hạng đầu vào Q·C1 và đầu ra Q·C2, chưa kể phản ứng. Phản ứng được tích hợp qua việc C1 đã là kết quả sau phản ứng ở bể 1.

2.2 Từ PT 7-22 Đến PT 7-23

Chia hai vế của PT 7-22 cho V/n:

dC2/dt = (n·Q/V)·C1 − (n·Q/V)·C2

dC2/dt = (n·Q/V)·C1 − (n·Q/V)·C2 (PT 7-23)

Đặt n·Q/V = n/τ (vì τ = V/Q), PT 7-23 có thể viết:

dC2/dt = (n/τ)·C1 − (n/τ)·C2

Đây là phương trình vi phân tuyến tính bậc nhất phi thuần nhất với biến C2 và hàm kích thích C1(t) — chưa biết dạng tường minh.

2.3 So Sánh PT 7-23 Với PT 7-9 (CSTR Đơn)

Nhắc lại PT 7-9 cho CSTR đơn:

dC/dt = −α·C + (Q/V)·C₀ = −(k + 1/τ)·C + C₀/τ

PT 7-23 cho bể 2 trong chuỗi:

dC2/dt = −(n/τ)·C2 + (n/τ)·C1(t)

Điểm tương đồng: Cả hai đều là ODE tuyến tính bậc nhất với số hạng kích thích (forcing term).

Điểm khác biệt:

• PT 7-9: hệ số là α = k + 1/τ, kích thích là C₀/τ = hằng số
• PT 7-23: hệ số là n/τ, kích thích là (n/τ)·C1(t) = hàm thời gian biến đổi

Chính sự biến đổi của C1(t) làm cho PT 7-23 phức tạp hơn PT 7-9.

Phần 3: Tham Số τ Và H — Hai Thước Đo Thời Gian Không Thể Thiếu

3.1 Thời Gian Lưu Thủy Lực Toàn Hệ τ = V/Q

τ = V/Q là thời gian lưu thủy lực (HRT) của toàn bộ chuỗi n bể — thời gian một phần tử nước thải trung bình lưu lại trong toàn hệ thống nếu không có phản ứng.

τ có đơn vị giờ và là thông số thiết kế cơ bản:

• τ nhỏ: lưu lượng lớn hoặc bể nhỏ → thời gian tiếp xúc ngắn → hiệu suất thấp
• τ lớn: lưu lượng nhỏ hoặc bể lớn → thời gian tiếp xúc dài → hiệu suất cao
• τ = thể tích/lưu lượng = V/Q — hoàn toàn do thiết kế quyết định

3.2 Thời Gian Vô Thứ Nguyên H = t/τ

GS. Lê Văn Cát định nghĩa tỉ lệ:

H = t/τ

H là thời gian vô thứ nguyên (dimensionless time):

• H = 0: lúc bắt đầu (t = 0)
• H = 1: khi t = τ (một thời gian lưu đã trôi qua)
• H = 2: khi t = 2τ (hai thời gian lưu đã trôi qua)
• H < 1: giai đoạn quá độ sớm
• H >> 1: giai đoạn tiệm cận trạng thái ổn định

Ý nghĩa vật lý của H: H biểu thị tỉ lệ giữa thời gian phản ứng thực tế t và thời gian lưu thủy lực τ của toàn hệ. Nó cũng tương ứng với tỉ lệ thời gian lưu tại một bể đơn so với thời gian lưu toàn hệ — vì thời gian lưu mỗi bể là τ/n, sau khi nước đi qua i bể: t_tích lũy = i × (τ/n), tức là H_sau bể i = i/n.

Tại sao dùng H thay vì t? Việc vô thứ nguyên hóa thời gian bằng τ cho phép so sánh các hệ thống có τ khác nhau trên cùng một trục — kết quả có tính tổng quát, không phụ thuộc vào giá trị cụ thể của V, Q.

3.3 Mối Liên Hệ n·H Và Thể Tích Đã Xử Lý

Tích n·H = n·t/τ = t/(τ/n) = t/τ_bể có ý nghĩa: số lần thời gian lưu của một bể đơn đã trôi qua.

Khi t = τ (H = 1): n·H = n — chính xác n lần thời gian lưu một bể đã qua.

Phần 4: Biểu Thức C1(t) — PT 7-24

4.1 C1 Từ Phản Ứng Bậc Nhất Ở Bể 1

Bể 1 là CSTR với thể tích V/n, lưu lượng Q, đầu vào C₀ hằng số, phản ứng bậc nhất k. Áp dụng kết quả PT 7-10 cho CSTR đơn (với C_in = C₀):

C1(t) = Q·C₀/(n·Q·(1/τ + k)/n·V·(1/τ + k)) ...

Nhưng GS. Lê Văn Cát chọn cách tiếp cận đơn giản hơn: sử dụng trực tiếp kết quả phản ứng bậc nhất (PT 7-3 — tương tự batch với thời gian = t, hằng số = k):

C1 = C₀·e^(-n·Q·t/V) = C₀·e^(-n·t/τ) = C₀·e^(-n·H) (PT 7-24)

4.2 Diễn Giải PT 7-24

PT 7-24 nói rằng: đầu ra của bể 1 (C1) biến thiên theo hàm mũ giảm dần với hệ số n/τ.

Tại t = 0: C1 = C₀·e⁰ = C₀ (đầu vào) Khi t → ∞: C1 → 0 (tiệm cận, nếu phản ứng hoàn toàn) Khi t = τ/n (H = 1/n): C1 = C₀·e⁻¹ ≈ 0,368·C₀

Tốc độ suy giảm C1 tỷ lệ với n: Với n bể, hệ số mũ là n/τ — tốc độ suy giảm nhanh gấp n lần so với CSTR đơn (hệ số 1/τ). Điều này phản ánh thực tế: mỗi bể có thể tích V/n nhỏ hơn → HRT ngắn hơn → dòng ra thay đổi nhanh hơn.

4.3 Kiểm Tra Tính Nhất Quán

Với n = 1 (1 bể): C1 = C₀·e^(-t/τ) — phản ứng bậc nhất trong bể đơn ✓

Với n → ∞ (PFR lý tưởng): C1 = C₀·e^(-∞·t/τ) → 0 ngay lập tức (trừ t = 0) — nước vào bể 1 "vô cùng nhỏ" bị xử lý ngay gần như tức thì ✓

Phần 5: Phương Trình Vi Phân Đầy Đủ Cho Bể 2 — PT 7-25

5.1 Thay C1 Vào PT 7-23

Thay PT 7-24 (C1 = C₀·e^(-n·H) = C₀·e^(-n·t/τ)) vào PT 7-23:

dC2/dt = (n/τ)·C₀·e^(-n·t/τ) − (n/τ)·C2 (PT 7-25)

Hay viết dạng chuẩn ODE:

dC2/dt + (n/τ)·C2 = (n/τ)·C₀·e^(-n·t/τ)

Đây là phương trình vi phân tuyến tính bậc nhất phi thuần nhất với hàm kích thích (n/τ)·C₀·e^(-n·t/τ) — hàm mũ giảm dần.

5.2 Phân Tích Cấu Trúc PT 7-25

GS. Lê Văn Cát nhận xét: "Phương trình 7-25 có dạng chính tắc (chuẩn) như phương trình 7-9." Điểm giống nhau:

Đặc điểm	PT 7-9 (CSTR đơn)	PT 7-25 (Bể 2 trong chuỗi)
Loại PT	ODE tuyến tính phi thuần nhất	ODE tuyến tính phi thuần nhất
Hệ số C	α = k + 1/τ	n/τ
Hàm kích thích	(Q/V)·C₀ = C₀/τ = hằng số	(n/τ)·C₀·e^(-n·t/τ) = biến đổi
Nghiệm	Tổng hàm mũ	Tổng hàm mũ (phức tạp hơn)

Điểm khác biệt then chốt: hàm kích thích của PT 7-25 là hàm mũ giảm, không phải hằng số. Điều này tạo ra cộng hưởng đặc biệt: hệ số của C2 trong PT 7-25 cũng là n/τ — bằng với hệ số mũ của hàm kích thích. Đây là trường hợp đặc biệt trong lý thuyết ODE gọi là "resonance" (cộng hưởng tham số), dẫn đến nghiệm có dạng đặc biệt.

5.3 Phương Pháp Giải PT 7-25

PT 7-25 thuộc loại: dy/dt + a·y = b·e^(-a·t) (hệ số mũ hàm kích thích = hệ số của y → cộng hưởng)

Phương pháp giải: Biến đổi biến u = C2·e^(n·t/τ):

dC2/dt = (du/dt)·e^(-n·t/τ) − (n/τ)·u·e^(-n·t/τ)

Thay vào PT 7-25:

(du/dt)·e^(-n·t/τ) − (n/τ)·u·e^(-n·t/τ) + (n/τ)·u·e^(-n·t/τ) = (n/τ)·C₀·e^(-n·t/τ)

du/dt = (n/τ)·C₀

Tích phân: u = (n/τ)·C₀·t + A

Quay lại: C2 = u·e^(-n·t/τ) = [(n/τ)·C₀·t + A]·e^(-n·t/τ)

Điều kiện đầu t = 0, C2 = C₀ (bể ban đầu chứa đầy nước thải):

C₀ = A → A = C₀

C2(t) = [(n/τ)·C₀·t + C₀]·e^(-n·t/τ) = C₀·[1 + n·t/τ]·e^(-n·t/τ)

Viết theo H = t/τ:

C2(H) = C₀·(1 + n·H)·e^(-n·H)

(Đây là dạng đầy đủ bao gồm cả điều kiện đầu. Khi bỏ qua điều kiện đầu và chỉ xét phần đặc biệt liên quan đến kích thích từ bể 1, ta được kết quả PT 7-26.)

Phần 6: Ứng Dụng Thực Tiễn Hệ Chuỗi CSTR Tại Việt Nam

6.1 Hệ A²/O — Ví Dụ Thực Tế 3 Bể Nối Tiếp

Hệ thống A²/O (Anaerobic-Anoxic-Oxic) là ứng dụng phổ biến nhất của nguyên lý chuỗi CSTR nối tiếp tại các KCN và đô thị Việt Nam:

• Bể 1 (Yếm khí): Đầu vào C₀ (COD, TP), đầu ra C1 (giải phóng P, phân hủy một phần COD)
• Bể 2 (Thiếu khí): Đầu vào C1, đầu ra C2 (khử nitrat NO₃⁻ → N₂)
• Bể 3 (Hiếu khí): Đầu vào C2, đầu ra C3 (oxy hóa COD còn lại, nitrat hóa, hấp thụ P)

Mỗi bể là CSTR với điều kiện khác nhau (kị khí/thiếu khí/hiếu khí). Đầu ra bể trước là đầu vào bể sau — đúng cấu trúc PT 7-22 và 7-23.

Nhưng trong thực tế A²/O, các bể không có cùng thể tích: thường V_yếm khí : V_thiếu khí : V_hiếu khí = 1 : 1 : 3. Đây là điểm sai lệch so với giả thiết "mỗi bể cùng thể tích" trong PT 7-22 — cần điều chỉnh công thức cho từng bể có HRT khác nhau.

6.2 Xác Định τ Và H Trong Vận Hành Thực Tế

Tại trạm XLNT khu đô thị (ví dụ: Q = 5.000 m³/ngày = 208 m³/h, V_tổng = 2.500 m³):

τ = V/Q = 2.500/208 = 12,0 giờ

Khi vận hành được t = 6 giờ sau khởi động: H = t/τ = 6/12 = 0,5

Với n = 3 bể nối tiếp (mỗi bể V/3 ≈ 833 m³): n·H = 3 × 0,5 = 1,5

Các giá trị này sẽ được dùng trực tiếp để tính C1, C2, C3 theo PT 7-24, 7-26, 7-27.

6.3 Tầm Quan Trọng Của PT 7-22 Trong Phân Tích Quá Độ

Phân tích quá độ (transient analysis) — không phải trạng thái ổn định — đặc biệt quan trọng khi:

• Khởi động hệ thống mới: Nồng độ trong mỗi bể thay đổi theo thời gian từ C₀ về trạng thái ổn định khác nhau
• Shock load (sốc tải trọng): C₀ đột ngột tăng → C1, C2 đáp ứng theo chuỗi thời gian trễ
• Bảo trì một bể trong chuỗi: Khi bể i dừng hoạt động, toàn bộ chuỗi sau đó bị ảnh hưởng

PT 7-22 và 7-23 chính là công cụ toán học để phân tích các tình huống này — không thể chỉ dùng phương trình trạng thái ổn định PT 7-13.

Kết Luận

Phương trình thăng bằng chất PT 7-22 cho bể thứ 2 trong chuỗi n CSTR nối tiếp thiết lập nền tảng toán học cho toàn bộ phân tích hệ nối tiếp. Hai tham số then chốt — thời gian lưu thủy lực toàn hệ τ = V/Q và thời gian vô thứ nguyên H = t/τ — cho phép diễn đạt kết quả ở dạng tổng quát, không phụ thuộc vào quy mô cụ thể.

Phương trình vi phân PT 7-25 cho bể 2 có cấu trúc tương tự PT 7-9 (CSTR đơn) nhưng với hàm kích thích là hàm mũ giảm (thay vì hằng số) — phản ánh thực tế rằng đầu vào bể 2 không cố định mà thay đổi liên tục theo diễn biến ở bể 1. Đây là khởi đầu của chuỗi phương trình dẫn đến nghiệm tổng quát PT 7-27 cho tất cả n bể, được phân tích trong bài viết tiếp theo của loạt bài này.