Hiệu Suất Quá Độ E(t) Theo PT 7-12 Và Ảnh Hưởng Của Lưu Lượng Dòng Đến CSTR: Phân Tích Động Học Và Ứng Dụng Khởi Động Hệ Thống XLNT

Khoảng Cách Giữa Khởi Động Và Trạng Thái Ổn Định

Trong vòng đời vận hành của một trạm xử lý nước thải (XLNT) tại Việt Nam, có những khoảnh khắc đặc biệt quan trọng khi hệ thống chưa đạt trạng thái ổn định: lần đầu khởi động sau xây dựng, khởi động lại sau bảo trì định kỳ, phục hồi sau sự cố điện hay máy bơm, hoặc điều chỉnh lưu lượng khi tải trọng thay đổi. Trong tất cả các tình huống này, câu hỏi thực tiễn là: hiệu suất xử lý hiện tại là bao nhiêu, và bao giờ hệ thống đạt mức xử lý theo tiêu chuẩn?

Câu trả lời nằm trong phương trình hiệu suất quá độ PT 7-12 — mô tả E(t) từ thời điểm bắt đầu vận hành đến khi đạt trạng thái ổn định. Đây là phương trình đặc trưng cho chế độ quá độ (transient regime) của bể phản ứng liên tục khuấy trộn đều (CSTR), khác biệt hoàn toàn so với hiệu suất trạng thái ổn định E_∞ theo PT 7-14.

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) đã chỉ ra tường minh ảnh hưởng của tốc độ dòng đến hiệu suất: “Tốc độ dòng càng lớn thì τ càng nhỏ (V = const), τ⁻¹ càng cao, hiệu suất phản ứng giảm.” Phân tích sâu điều này và phương trình PT 7-12 là nội dung trọng tâm của bài viết.

Phần 1: Dẫn Xuất Phương Trình Hiệu Suất Quá Độ PT 7-12

1.1 Từ Nghiệm Quá Độ PT 7-10 Đến Hiệu Suất E(t)

Từ nghiệm quá độ PT 7-10 (với điều kiện ban đầu C_in = C₀, tức là bể ban đầu chứa nước chưa xử lý):

C(t) = Q·C₀/(V·α) · (1 − e^(−α·t)) + C₀·e^(−α·t)

Trong đó α = k + Q/V = k + 1/τ

Sử dụng Q/V = 1/τ:

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

Tính hiệu suất E(t) = 1 − C(t)/C₀:

E(t) = 1 − [1/(k·τ + 1) · (1 − e^(−α·t)) + e^(−α·t)]

= 1 − 1/(k·τ+1) + e^(−α·t)/(k·τ+1) − e^(−α·t)

= k·τ/(k·τ+1) − e^(−α·t) · [1 − 1/(k·τ+1)]

= k·τ/(k·τ+1) − e^(−α·t) · k·τ/(k·τ+1)

E(t) = [k·τ/(k·τ+1)] · (1 − e^(−(1/τ+k)·t)) (PT 7-12 — dạng triển khai)

Hay viết theo ký hiệu gốc:

E(t) = [1 − 1/(1+k·τ)] · (1 − e^(−(k+1/τ)·t)) (PT 7-12)

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

Tại t = 0: E(0) = E_∞ · (1 − e⁰) = E_∞ · 0 = 0 ✓ (ban đầu chưa có xử lý)

Khi t → ∞: E(∞) = E_∞ · (1 − 0) = E_∞ = k·τ/(1+k·τ) ✓ (đạt trạng thái ổn định PT 7-14)

Dạng tổng quát: E(t) = E_∞ · (1 − e^(−α·t)) — tích của giá trị tiệm cận E_∞ và hàm tăng dần (1 − e^(−α·t))

1.3 Phân Tích Cấu Trúc PT 7-12

PT 7-12 là tích của hai thành phần:

Thành phần 1: E_∞ = k·τ/(1+k·τ) — đây là trần hiệu suất (efficiency ceiling) mà hệ thống tiệm cận. Phụ thuộc hoàn toàn vào k và τ — không phụ thuộc thời gian.

Thành phần 2: (1 − e^(−α·t)) — đây là hàm tốc độ tiến đến trần. Phụ thuộc vào α = k + 1/τ — đặc trưng tốc độ đạt trạng thái ổn định.

Điều kỳ thú: cùng một α có thể tạo ra từ các cặp (k, τ) khác nhau, nhưng mỗi cặp cho E_∞ khác nhau. Ví dụ:

→ Hai hệ thống có cùng α (ví dụ α = 1,1 giờ⁻¹) có tốc độ quá độ như nhau nhưng E_∞ rất khác nhau. Điều này có nghĩa là tốc độ đạt trạng thái ổn định không quyết định chất lượng trạng thái ổn định.

Phần 2: Ảnh Hưởng Của Lưu Lượng Dòng Q Đến Hiệu Suất E

2.1 Cơ Chế Ảnh Hưởng Của Q

GS. Lê Văn Cát phân tích: “Tốc độ dòng càng lớn thì τ càng nhỏ (V = const), τ⁻¹ = Q/V càng cao, hiệu suất phản ứng giảm.”

Định lượng hóa: Với V cố định, Q ảnh hưởng đến cả E_∞ và tốc độ quá độ qua α:

Khi Q tăng n lần (τ giảm n lần):

• E_∞ = k·τ/(1+k·τ) = k·(τ/n) / (1+k·τ/n) — giảm rõ rệt khi Da = k·τ/n giảm
• α = k + n/τ — tăng, tức là đạt trạng thái ổn định nhanh hơn, nhưng trạng thái ổn định đó tệ hơn

Nghịch lý thú vị: Q tăng → hệ đạt trạng thái ổn định nhanh hơn (α lớn hơn) NHƯNG trạng thái ổn định đó cho E_∞ thấp hơn.

2.2 Ví Dụ Định Lượng — Bể Aerotank Khi Tải Trọng Thủy Lực Thay Đổi

Bể aerotank điển hình: V = 2.000 m³, k = 0,2 giờ⁻¹, C₀ = 200 mg BOD/L

Kịch bản A: Vận hành bình thường, Q = 250 m³/h

• τ = 2000/250 = 8 giờ
• α = 0,2 + 1/8 = 0,325 giờ⁻¹
• E_∞ = (0,2×8)/(1+0,2×8) = 1,6/2,6 = 61,5%
• t_98% = 4/0,325 ≈ 12,3 giờ

Kịch bản B: Mùa mưa, Q = 500 m³/h (tăng gấp đôi)

• τ = 2000/500 = 4 giờ
• α = 0,2 + 1/4 = 0,45 giờ⁻¹
• E_∞ = (0,2×4)/(1+0,2×4) = 0,8/1,8 = 44,4%
• t_98% = 4/0,45 ≈ 8,9 giờ

Kịch bản C: Hạ tải, Q = 125 m³/h (giảm một nửa)

• τ = 2000/125 = 16 giờ
• α = 0,2 + 1/16 = 0,2625 giờ⁻¹
• E_∞ = (0,2×16)/(1+0,2×16) = 3,2/4,2 = 76,2%
• t_98% = 4/0,2625 ≈ 15,2 giờ

→ Mùa mưa: Q tăng gấp đôi → E_∞ giảm từ 61,5% xuống 44,4% (giảm 17%). Mùa khô hạ tải: E_∞ tăng lên 76,2%. Sự dao động lớn này là nguyên nhân chính gây khó khăn cho công tác quản lý chất lượng đầu ra quanh năm.

2.3 Bất Đẳng Thức Q < Q_max — Điều Kiện Duy Trì E Tối Thiểu

Từ PT 7-14, để duy trì E ≥ E_min (ví dụ E_min = 85% theo QCVN):

k·τ/(1+k·τ) ≥ E_min → k·τ ≥ E_min/(1-E_min)

k·V/Q ≥ E_min/(1-E_min)

Q ≤ Q_max = k·V·(1-E_min)/E_min

Ví dụ: V = 2.000 m³, k = 0,2 giờ⁻¹, E_min = 85%: Q_max = 0,2 × 2.000 × 0,15/0,85 = 400 × 0,1765 = 70,6 m³/h

→ Bể chỉ có thể xử lý tối đa Q = 70,6 m³/h để đạt E = 85%. Đây là công suất thiết kế thực của bể đối với yêu cầu xử lý cụ thể.

Thực tế: kỹ sư vận hành cần biết Q_max cho từng chỉ tiêu (BOD, NH₄-N, TP) và điều chỉnh không cho Q vượt giới hạn tương ứng.

Phần 3: Động Học Quá Độ E(t) — Ứng Dụng Trong Thực Tiễn Vận Hành

3.1 Thời Gian Đạt Ngưỡng An Toàn E_safe < E_∞

Trong thực tế, kỹ sư muốn biết: sau bao lâu kể từ khi khởi động, hệ thống đạt mức E_safe đủ để bắt đầu xả nước thải ra môi trường mà không vi phạm QCVN?

Từ PT 7-12: E(t) = E_∞ · (1 − e^(−α·t)) ≥ E_safe

1 − e^(−α·t) ≥ E_safe/E_∞

e^(−α·t) ≤ 1 − E_safe/E_∞

t_safe = −(1/α) · ln(1 − E_safe/E_∞)

Ví dụ thực tế:

Bể aerotank mới khởi động (V = 1.000 m³, Q = 100 m³/h, k = 0,2 giờ⁻¹, C₀ = 200 mg BOD/L):

• τ = 10 giờ; α = 0,2 + 0,1 = 0,3 giờ⁻¹
• E_∞ = 2/(1+2) = 66,7%
• Yêu cầu E_safe = 85%? → Không thể! E_∞ = 66,7% < 85%. Cần điều chỉnh thiết kế trước.
• Yêu cầu E_safe = 60%: t_safe = −(1/0,3) · ln(1 − 60/66,7) = −3,33 · ln(0,10) = 7,7 giờ

→ Hệ thống cần khoảng 8 giờ vận hành trước khi đạt 60% hiệu suất — tức là đến 20:00 nếu khởi động lúc 12:00. Đây là thông tin thiết thực để lên kế hoạch khởi động vận hành.

3.2 Quá Trình Điều Chỉnh Sau Thay Đổi Lưu Lượng

Một tình huống quan trọng: hệ thống đang vận hành ở trạng thái ổn định với Q₁, sau đó đột ngột chuyển sang Q₂ (tăng hoặc giảm tải). Hệ thống sẽ điều chỉnh từ E_∞₁ về E_∞₂.

Phương trình mô tả quá trình chuyển tiếp (từ điều kiện ban đầu C = C_∞₁ khi đang ở trạng thái ổn định cũ):

E(t) = E_∞₂ − (E_∞₂ − E_∞₁) · e^(−α₂·t)

Trong đó α₂ = k + Q₂/V và E_∞₁, E_∞₂ là hiệu suất trạng thái ổn định tương ứng với Q₁ và Q₂.

Ứng dụng: Điều tiết Q khi cần tăng E tạm thời

Kỹ sư vận hành có thể giảm Q tạm thời (giảm tải) khi cần tăng E nhanh (ví dụ trước khi thanh tra môi trường):

• Giảm Q từ 250 xuống 150 m³/h → τ tăng từ 8 lên 13,3 giờ
• E_∞_mới = 0,2×13,3/(1+0,2×13,3) = 2,66/3,66 = 72,7% (tăng từ 61,5%)
• Hệ thống đạt 90% E_∞_mới sau: t = −(1/α₂)·ln(0,1) = −(1/0,275)·ln(0,1) = 8,4 giờ

Tuy nhiên: giảm Q tạm thời đồng nghĩa tích tụ nước thải chưa xử lý ở đầu vào — chỉ dùng khi có bể điều hòa đầy đủ.

3.3 Phân Tích Hiệu Suất Tức Thời Trong Vận Hành Ca

Nhiều trạm XLNT tại Việt Nam vận hành theo 3 ca (6-14h, 14-22h, 22-6h). Tải trọng thay đổi giữa các ca. Dùng PT 7-12, kỹ sư ca trực có thể ước tính E(t) thực tế tại bất kỳ thời điểm nào trong ca mình:

Công thức đơn giản hóa (tuyến tính hóa gần đúng):

Khi α·t không quá lớn (α·t < 0,5): 1 − e^(−α·t) ≈ α·t

E(t) ≈ E_∞ · α · t (chỉ đúng khi hệ mới khởi động hoặc sau thay đổi lớn)

Khi α·t > 3 (hệ đã gần ổn định): E(t) ≈ E_∞ · (1 − e^(−α·t)) ≈ E_∞ (sai số < 5%)

Trong thực tế: sau 3-4 lần hằng số thời gian 1/α, hệ thống xem như đã ổn định. Kỹ thuật viên vận hành chỉ cần nhớ: t_ổn định ≈ 4/α = 4/(k + 1/τ).

Phần 4: So Sánh Chế Độ Quá Độ CSTR Với Bể Gián Đoạn SBR

4.1 Bảng So Sánh Toàn Diện

4.2 Khi Nào Ưu Tiên SBR, Khi Nào Ưu Tiên CSTR?

Ưu tiên SBR (bể gián đoạn):

• Lưu lượng nhỏ (Q < 2.000 m³/ngày) — phổ biến ở khu dân cư, bệnh viện, khách sạn tại Việt Nam
• Yêu cầu E cao (> 90%) — SBR đạt được bằng thời gian phản ứng dài mà không cần thể tích bể quá lớn
• Yêu cầu xử lý N và P trong cùng một bể — SBR đa pha linh hoạt hơn CSTR
• Tải trọng dao động mạnh — chu kỳ SBR dễ điều chỉnh

Ưu tiên CSTR liên tục:

• Lưu lượng lớn và liên tục (Q > 5.000 m³/ngày) — khu đô thị, KCN lớn
• Yêu cầu vận hành đơn giản, ít van và van tự động
• Nước thải đã tương đối đồng đều về thành phần và lưu lượng
• Hệ thống A²/O hay chuỗi CSTR nối tiếp cho E tổng hợp cao

4.3 Ứng Dụng Tại Trạm XLNT Đô Thị Việt Nam

Điển hình là trạm XLNT Yên Sở (Hà Nội), Bình Hưng (TP.HCM) và Hòa Xuân (Đà Nẵng) — đều vận hành theo nguyên lý CSTR liên tục với dòng chảy Q hàng nghìn m³/giờ. Mùa mưa, khi Q tăng đột biến, hiệu suất E_∞ giảm theo PT 7-14 — đây là thách thức vận hành thực tế quan trọng.

Phần 5: Phân Tích Nhạy Cảm — Tham Số Nào Ảnh Hưởng Nhiều Nhất Đến E(t)?

5.1 Phân Tích Đạo Hàm Riêng Của E(t) Theo Các Biến

Từ PT 7-12: E(t) = [k·τ/(1+k·τ)] · [1 − e^(−(k+1/τ)·t)]

Đạo hàm theo k: ∂E/∂k = [τ/(1+k·τ)²] · [1−e^(−α·t)] + [k·τ/(1+k·τ)] · t · e^(−α·t)

Cả hai số hạng dương → E tăng khi k tăng, tác động qua cả E_∞ và α.

Đạo hàm theo τ (hay V): ∂E/∂τ = [k/(1+k·τ)²] · [1−e^(−α·t)] − [k·τ/(1+k·τ)] · (1/τ²) · t · e^(−α·t)

Số hạng thứ nhất dương (tăng E_∞), số hạng thứ hai âm (giảm α, làm chậm đạt ổn định). Tác động net: tăng τ luôn tăng E(t) ở trạng thái ổn định, nhưng tốc độ tăng E(t) trong quá độ chậm hơn.

Hàm ý thực tế:

• Để tăng E nhanh trong ngắn hạn (ví dụ khởi động khẩn cấp): ưu tiên tăng k (cấy bùn giống, điều chỉnh nhiệt độ, bổ sung dinh dưỡng vi lượng)
• Để tăng E trong dài hạn (thiết kế bền vững): ưu tiên tăng τ (thể tích bể) — tác động ổn định và không phụ thuộc điều kiện vận hành biến động

5.2 Ảnh Hưởng Của Nhiệt Độ Đến E(t) Và α

Nhiệt độ ảnh hưởng đến k: k_T = k₂₀ × (1,04)^(T-20)

Nhiệt độ không ảnh hưởng trực tiếp đến τ = V/Q (thông số thủy lực, không phụ thuộc nhiệt độ).

Do đó: α_T = k_T + 1/τ — nhiệt độ cao tăng k_T, tăng α_T, làm hệ đạt trạng thái ổn định nhanh hơn VÀ E_∞ cao hơn.

Ví dụ: Cùng bể, cùng Q, chỉ khác nhiệt độ:

• T = 15°C: k = 0,164 giờ⁻¹, τ = 8h → α = 0,289 giờ⁻¹, E_∞ = 56,8%, t_98% = 13,8h
• T = 25°C: k = 0,219 giờ⁻¹, τ = 8h → α = 0,344 giờ⁻¹, E_∞ = 63,7%, t_98% = 11,6h
• T = 30°C: k = 0,266 giờ⁻¹, τ = 8h → α = 0,391 giờ⁻¹, E_∞ = 68,1%, t_98% = 10,2h

Từ mùa đông (15°C) sang mùa hè (30°C): E_∞ tăng 11,3%, thời gian ổn định rút ngắn 3,6 giờ. Tại TP.HCM nhiệt độ quanh năm 27-32°C là điều kiện thuận lợi so với Hà Nội có mùa đông lạnh.

Phần 6: Tối Ưu Hóa Vận Hành Dựa Trên PT 7-12

6.1 Xác Định Q Tối Ưu Trong Quá Trình Khởi Động

Mục tiêu: đạt E ≥ E_yêu cầu sớm nhất có thể trong giai đoạn khởi động.

Điều kiện cần: E_∞ ≥ E_yêu cầu → Q ≤ Q_max = k·V·(1-E_yêu cầu)/E_yêu cầu

Tối ưu tốc độ đạt E_yêu cầu: E(t*) = E_yêu cầu

Từ PT 7-12: E_yêu cầu = E_∞ · (1 − e^(−α·t*))

t* = −(1/α)·ln(1 − E_yêu cầu/E_∞)

Để tối thiểu t*, cần tối thiểu −(1/α)·ln(1 − r) trong đó r = E_yêu cầu/E_∞ < 1.

Khi α tăng (Q tăng hoặc V giảm): 1/α giảm (tốt), nhưng E_∞ giảm → r tăng → −ln(1-r) tăng (xấu). Có điểm tối ưu.

Chiến lược thực tiễn:

Vận hành với Q thấp hơn thiết kế trong 24-48 giờ đầu khởi động (Q_khởi động = 50-60% Q_thiết kế), sau đó tăng dần Q về giá trị thiết kế theo từng bước 10-20%, kiểm tra E trước mỗi lần tăng.

6.2 Quy Trình Khởi Động Cụ Thể Cho Trạm XLNT Tại Việt Nam

Dựa trên phân tích PT 7-12, quy trình khởi động đề xuất:

Giai đoạn 1 (0-24 giờ): Chạy tuần hoàn nội

• Không nạp nước thải thô vào
• Cấy bùn từ trạm XLNT lân cận, mật độ X₀ ≥ 2.000 mg VSS/L
• Sục khí liên tục để acclimate vi sinh
• Mục tiêu: k đạt > 50% giá trị thiết kế

Giai đoạn 2 (24-72 giờ): Nạp dần nước thải

• Tăng Q từ 20% → 40% → 60% Q_thiết kế, mỗi bước 12 giờ
• Theo dõi: pH, DO, SV30, COD đầu ra sau mỗi bước
• Chưa xả ra nguồn tiếp nhận (tuần hoàn về đầu vào hoặc chứa trong bể điều hòa)

Giai đoạn 3 (72-120 giờ): Chạy thử ở 80% Q

• Q = 80% Q_thiết kế → τ_thực = 1,25·τ_thiết kế → E_∞thực > E∞_thiết kế
• Kiểm tra chất lượng đầu ra theo QCVN, đặc biệt BOD₅, NH₄-N
• Khi E ≥ E_yêu cầu trong 3 ngày liên tiếp: chuyển sang vận hành chính thức

Giai đoạn 4: Vận hành chính thức

• Q = 100% Q_thiết kế
• Theo dõi E_∞ theo PT 7-14: kiểm tra hàng ngày bằng đo nhanh COD (COD test kit)
• Duy trì E ≥ E_yêu cầu bằng điều chỉnh SRT (bơm xả bùn) và Q khi cần

6.3 Xử Lý Sự Cố Mất Điện — Phục Hồi Nhanh

Sau sự cố mất điện (ngừng sục khí, ngừng bơm) trong thời gian t_sự cố:

Bể không sục khí → vi khuẩn hiếu khí chết dần, k giảm về k_thực_sau = f(t_sự cố) Khi khởi động lại, hệ thống ở trạng thái ban đầu C = C_trước_sự_cố, k = k_thực_sau

Áp dụng PT 7-12 với C_in = C_trước_sự_cố và k = k_thực_sau để ước tính t_phục hồi.

Nếu t_sự cố < 4 giờ: k giảm < 20%, hệ phục hồi trong 6-8 giờ Nếu t_sự cố > 12 giờ: cần cấy bổ sung bùn giống, thời gian phục hồi 24-48 giờ

Kết Luận

Phương trình hiệu suất quá độ PT 7-12 — E(t) = E_∞ · (1 − e^(−α·t)) — là cầu nối giữa lý thuyết CSTR và thực tiễn vận hành. Nó mô tả đầy đủ diễn biến hiệu suất xử lý từ khi hệ thống khởi động (E = 0) đến khi đạt trạng thái ổn định (E = E_∞), với hằng số thời gian 1/α = 1/(k+1/τ).

Phân tích ảnh hưởng của lưu lượng Q xác nhận quan sát của GS. Lê Văn Cát: Q tăng làm giảm τ, giảm E_∞ theo PT 7-14, dẫn đến hiệu suất xử lý thấp hơn dù hệ đạt trạng thái ổn định nhanh hơn. Đây là nguyên nhân gốc của các vấn đề chất lượng đầu ra vào mùa mưa tại các trạm XLNT đô thị Việt Nam, và là cơ sở kỹ thuật để thiết kế bể điều hòa và xây dựng quy trình vận hành linh hoạt theo mùa.