WesterntechVN – Từ Nghiệm Đặc Biệt Đến Công Thức Tổng Quát

Trong kỹ thuật phản ứng gián đoạn (batch reactor), các phương trình bậc không (PT 7-2), bậc nhất (PT 7-3) và bậc hai (PT 7-4) là những trường hợp đặc biệt của một công thức tổng quát hơn. Khi bậc phản ứng n không nguyên (n = 0,5; 0,7; 1,5…) hoặc khi muốn kiểm tra nhiều giả thuyết về bậc phản ứng cùng lúc, cần đến phương trình nghiệm tổng quát bậc n — PT 7-5.

Song song với việc mô tả sự biến đổi nồng độ C(t), kỹ sư thiết kế hệ thống xử lý nước thải thực tế không chỉ quan tâm đến giá trị nồng độ tuyệt đối mà còn cần đánh giá mức độ xử lý đạt được — đây chính là vai trò của hiệu suất phản ứng E theo PT 7-6 và PT 7-7.

Hai bộ công thức này — PT 7-5 về nghiệm bậc n và PT 7-6, 7-7 về hiệu suất E — cùng nhau tạo nên công cụ thiết kế hoàn chỉnh để xác định thời gian lưu nước (HRT) cần thiết trong bể gián đoạn nhằm đạt mục tiêu xử lý cụ thể theo tiêu chuẩn Việt Nam.

Phần 1: Phương Trình Nghiệm Bậc N Tổng Quát — PT 7-5

1.1 Dẫn Xuất Phương Trình

Với bậc phản ứng n tổng quát (n ≠ 1), phương trình tốc độ PT 7-1 trở thành:

-dC/dt = k · Cⁿ

Tách biến và tích phân với điều kiện t = 0 thì C = C₀:

∫(từ C₀ đến C) C^(-n) · dC = -k · ∫(từ 0 đến t) dt

C^(1-n)/(1-n) – C₀^(1-n)/(1-n) = -k·t

Nhân hai vế với (1-n):

C^(1-n) – C₀^(1-n) = -(1-n)·k·t = (n-1)·k·t

Suy ra:

1/C^(n-1) = 1/C₀^(n-1) + (n-1)·k·t (PT 7-5)

Hay viết dạng tường minh:

C(t) = [C₀^(1-n) – (n-1)·k·t]^[1/(1-n)] (với n ≠ 1)

1.2 Kiểm Tra Tính Nhất Quán Với Các Trường Hợp Đặc Biệt

PT 7-5 phải thu gọn về PT 7-2 và PT 7-4 khi n = 0 và n = 2:

n = 0: PT 7-5 cho: 1/C^(-1) = 1/C₀^(-1) + (-1)·k·t → C = C₀ – kt ✓ (PT 7-2)

n = 2: PT 7-5 cho: 1/C^1 = 1/C₀^1 + (1)·k·t → 1/C = 1/C₀ + kt ✓ (PT 7-4)

n = 1: PT 7-5 không áp dụng trực tiếp (chia cho 0), nhưng giới hạn khi n → 1 cho kết quả PT 7-3: C = C₀·e^(-kt) ✓

1.3 Đặc Điểm Của Phản Ứng Bậc Phân Số (0 < n < 1)

Bậc phản ứng phân số thường gặp trong:

Xử lý sinh học hỗn hợp cơ chất: Khi nước thải chứa nhiều loại hợp chất hữu cơ, mỗi loại có tốc độ phân hủy riêng, tổng hợp quan sát được cho hành vi bậc phân số. Ví dụ: nước thải công nghiệp thực phẩm tại Đồng Nai thường có bậc phản ứng quan sát được n ≈ 0,6-0,8 đối với COD tổng.

Phân hủy trong đất và trầm tích: Các chất ô nhiễm (dư lượng thuốc trừ sâu, hydrocarbons) trong đất tại các khu vực nông nghiệp Đồng bằng sông Cửu Long thường phân hủy với bậc phân số do quá trình khuếch tán trong lỗ rỗng đất chi phối.

Bậc Langmuir-Hinshelwood: Trong oxy hóa xúc tác quang (TiO₂/UV) xử lý thuốc nhuộm, tốc độ tuân theo cơ chế Langmuir-Hinshelwood, tương đương bậc 0,5-0,8 trong điều kiện nồng độ trung bình.

1.4 Xác Định Bậc n Từ Thực Nghiệm — Phương Pháp Log-Log

Từ PT 7-1: vi = k·Cⁿ

Lấy logarithm hai vế:

log(vi) = log(k) + n·log(C)

Vẽ đồ thị log(vi) theo log(C) — nếu tuyến tính:

• Độ dốc = n (bậc phản ứng)
• Giao điểm trục tung = log(k)

Quy trình thực hành:

1. Đo nồng độ Ci tại các thời điểm ti
2. Tính vi ≈ -(Ci+1 – Ci-1)/(ti+1 – ti-1) (đạo hàm số trung tâm)
3. Vẽ log(vi) vs log(Ci) và hồi quy tuyến tính
4. Đọc n từ độ dốc

Phần 2: Hiệu Suất Phản Ứng E — PT 7-6 Và Ý Nghĩa Công Nghệ

2.1 Định Nghĩa Và Ý Nghĩa Của E

Hiệu suất phản ứng E (conversion efficiency, hay degree of conversion) được định nghĩa theo PT 7-6:

E = (C₀ – C) / C₀ = 1 – C/C₀ (PT 7-6)

Trong đó:

• C₀: nồng độ cơ chất đầu vào (mg/L hoặc g/m³)
• C: nồng độ cơ chất tại thời điểm t hoặc nồng độ đầu ra Ce
• E: hiệu suất phản ứng (vô thứ nguyên, giá trị từ 0 đến 1, hay biểu diễn %)

Ý nghĩa vật lý: E = 0,9 nghĩa là 90% lượng cơ chất ban đầu đã được chuyển hóa — chỉ còn lại 10% C₀ trong pha lỏng.

Lưu ý quan trọng: Trong xử lý nước thải thực tế tại Việt Nam, E có thể được định nghĩa theo nhiều chỉ tiêu khác nhau tùy mục đích:

• E_COD = (COD_vào – COD_ra) / COD_vào → hiệu quả xử lý carbon hữu cơ
• E_NH4 = (NH₄-N_vào – NH₄-N_ra) / NH₄-N_vào → hiệu quả nitrat hóa
• E_TP = (TP_vào – TP_ra) / TP_vào → hiệu quả khử photpho

Mỗi giá trị E có ý nghĩa riêng và phản ánh một quá trình vi sinh vật khác nhau trong hệ thống.

2.2 Tiêu Chuẩn E Theo Quy Chuẩn Việt Nam

Để đạt các tiêu chuẩn xả thải, cần E tối thiểu:

QCVN 14:2008/BTNMT (Nước thải sinh hoạt — Cột A):

• BOD₅ ≤ 30 mg/L: Với C₀ = 200 mg/L → E_BOD ≥ 85%
• COD ≤ 75 mg/L: Với C₀ = 400 mg/L → E_COD ≥ 81%
• NH₄-N ≤ 5 mg/L: Với C₀ = 40 mg/L → E_N ≥ 87,5%
• Tổng P ≤ 6 mg/L: Với C₀ = 8 mg/L → E_P ≥ 25%

QCVN 40:2011/BTNMT (Nước thải công nghiệp — Cột A):

• BOD₅ ≤ 30 mg/L: Với C₀ = 1.000 mg/L (chế biến thực phẩm) → E_BOD ≥ 97%
• COD ≤ 75 mg/L: Với C₀ = 3.000 mg/L → E_COD ≥ 97,5%
• Tổng N ≤ 20 mg/L: Với C₀ = 60 mg/L → E_N ≥ 67%
• Tổng P ≤ 4 mg/L: Với C₀ = 15 mg/L → E_P ≥ 73%

Đây là các mục tiêu E cụ thể mà kỹ sư thiết kế phải đạt được.

Phần 3: Phương Trình Hiệu Suất Phản Ứng Bậc Nhất — PT 7-7

3.1 Dẫn Xuất PT 7-7 Từ PT 7-3 Và PT 7-6

Kết hợp PT 7-3 (C = C₀·e^(-kt)) và PT 7-6 (E = 1 – C/C₀):

E = 1 – C/C₀ = 1 – C₀·e^(-kt)/C₀

E = 1 – e^(-k·t) (PT 7-7)

3.2 Phân Tích Đặc Tính Của E(t) Theo PT 7-7

Tính chất toán học:

• t = 0: E = 1 – e⁰ = 0 (hiệu suất ban đầu = 0, chưa có xử lý)
• t → ∞: E → 1 – 0 = 1 (hiệu suất tối đa = 100%, tiệm cận, không bao giờ đạt đúng 1)
• Hàm đơn điệu tăng, lõm xuống (dE/dt > 0, d²E/dt² < 0)
• Tốc độ tăng E chậm dần theo thời gian

Điểm uốn và tối ưu kinh tế: Tốc độ thay đổi hiệu suất: dE/dt = k·e^(-kt). Tốc độ này giảm liên tục — điều này có nghĩa là mỗi đơn vị thời gian thêm vào mang lại mức tăng E ngày càng nhỏ hơn. Đây là cơ sở kinh tế để lựa chọn HRT tối ưu (không nên kéo dài HRT vô hạn mà chỉ cần đủ để đạt E_yêu cầu).

3.3 Đảo Ngược PT 7-7: Tính Thời Gian Từ Mục Tiêu E

Từ PT 7-7: E = 1 – e^(-kt)

Đảo ngược để tính thời gian phản ứng cần thiết khi biết E mục tiêu và k:

e^(-kt) = 1 – E

-kt = ln(1 – E)

t_yêu cầu = -ln(1 – E) / k = ln[1/(1-E)] / k

Đây là công thức thiết kế cơ bản cho bể gián đoạn với phản ứng bậc nhất.

Bảng tra nhanh: Với k = 0,2 giờ⁻¹ (ví dụ điển hình nước thải sinh hoạt, 25°C):

Nhận xét: Để tăng E từ 90% lên 99% (loại thêm 9%), cần thêm 11,52 giờ — gần bằng thời gian đã tốn cho 90% đầu tiên. Điều này phản ánh đặc tính “lợi suất giảm dần” đặc trưng của phản ứng bậc nhất.

3.4 Hiệu Suất E Cho Phản Ứng Bậc Khác Nhất

Để so sánh toàn diện, biểu diễn E theo t cho các bậc phản ứng khác nhau:

Bậc 0 (từ PT 7-2): E = (C₀ – C)/C₀ = kt/C₀ (tuyến tính, tăng đều đến E = 1 tại t* = C₀/k)

Bậc 1 (PT 7-7): E = 1 – e^(-kt) (tiệm cận, tăng nhanh ban đầu rồi chậm dần)

Bậc 2 (từ PT 7-4): E = k·C₀·t / (1 + k·C₀·t) (dạng hyperbola, tăng nhanh ban đầu rồi chậm rất rõ)

Bậc n tổng quát (từ PT 7-5): E = 1 – [1 + (n-1)·k·C₀^(n-1)·t]^(-1/(n-1))

So sánh tốc độ đạt E = 80% (ví dụ, k·C₀ = 1 đơn vị/giờ):

• Bậc 0: t = 0,8 giờ (nhanh nhất)
• Bậc 1: t = 1,609 giờ
• Bậc 2: t = 4,0 giờ (chậm nhất, đặc biệt ở nồng độ thấp)

Phần 4: Tích Hợp PT 7-5, 7-6, 7-7 Vào Quy Trình Thiết Kế Bể Gián Đoạn

4.1 Quy Trình Thiết Kế 5 Bước

Bước 1: Xác định yêu cầu đầu vào và đầu ra

• Đặc trưng nước thải: C₀ (nồng độ đầu vào), Q (lưu lượng ngày)
• Tiêu chuẩn đầu ra: Ce (nồng độ tối đa cho phép theo QCVN)
• Tính E mục tiêu: E = 1 – Ce/C₀

Bước 2: Xác định bậc phản ứng và hằng số k từ thực nghiệm

• Thực hiện thí nghiệm batch với nước thải thực tế
• Vẽ đồ thị tuyến tính hóa (C vs t / ln(C) vs t / 1/C vs t / log(v) vs log(C))
• Xác định n và k từ hồi quy tuyến tính

Bước 3: Tính thời gian phản ứng t_rx từ E mục tiêu

• Bậc 1: t_rx = -ln(1-E)/k = ln[C₀/Ce]/k
• Bậc n: t_rx = [C₀^(1-n) – Ce^(1-n)] / [(n-1)·k]
• Thêm hệ số an toàn SF = 1,2-1,5 → t_thiết kế = SF · t_rx

Bước 4: Tính chu kỳ vận hành và thể tích bể

• Thời gian một chu kỳ: T_ck = t_nạp + t_rx + t_lắng + t_tháo
  • t_nạp: 0,5-1 giờ
  • t_lắng: 1,5-3 giờ (phụ thuộc tính lắng của bùn)
  • t_tháo: 0,5-1 giờ
• Số chu kỳ/ngày: n_ck = 24/T_ck
• Thể tích bể: V = Q / n_ck

Bước 5: Kiểm tra và hiệu chỉnh theo nhiệt độ mùa đông

• k_đông = k₂₀ · (1,04)^(T_đông – 20) (miền Bắc: T_đông ≈ 15°C → k_đông ≈ 0,82·k₂₀)
• t_rx,đông = -ln(1-E)/k_đông → dài hơn
• Kiểm tra: T_ck,đông vẫn cho phép đủ số chu kỳ/ngày với thể tích đã chọn

4.2 Ví Dụ Thiết Kế Hoàn Chỉnh: Bể SBR Xử Lý Nước Thải Sinh Hoạt

Thông số đầu vào:

• Q = 500 m³/ngày
• C₀_BOD = 250 mg/L, C₀_NH₄-N = 45 mg/L
• Yêu cầu: Ce_BOD ≤ 30 mg/L, Ce_NH₄-N ≤ 5 mg/L (QCVN 14:2008/BTNMT, Cột A)

Tính E mục tiêu:

• E_BOD = 1 – 30/250 = 0,88 (88%)
• E_N = 1 – 5/45 = 0,889 (88,9%)

Xác định k từ thực nghiệm (T = 25°C):

• k_BOD = 0,20 giờ⁻¹ (từ thí nghiệm batch — phản ứng bậc 1)
• k_N = 0,15 giờ⁻¹ (nitrat hóa — bậc 1 theo NH₄-N)

Tính thời gian phản ứng:

• t_rx,BOD = -ln(1-0,88)/0,20 = -ln(0,12)/0,20 = 2,12/0,20 = 10,6 giờ
• t_rx,N = -ln(1-0,889)/0,15 = -ln(0,111)/0,15 = 2,20/0,15 = 14,7 giờ
• Khống chế bởi nitrat hóa: t_rx = 14,7 giờ
• Với SF = 1,2: t_thiết kế = 1,2 × 14,7 = 17,6 giờ

Chu kỳ vận hành và thể tích:

• T_ck = 0,5 + 17,6 + 2 + 0,5 = 20,6 giờ → chỉ 1 chu kỳ/ngày (quá dài)
• Giải pháp: thiết kế 2 bể SBR vận hành lệch pha nhau 10 giờ
• Thể tích mỗi bể: V = Q = 500 m³ (xử lý toàn bộ lưu lượng ngày)

Kiểm tra mùa đông Hà Nội (T = 15°C):

• k_N,đông = 0,15 × (1,04)^(15-25) = 0,15 × (1,04)^(-10) = 0,15 × 0,676 = 0,101 giờ⁻¹
• t_rx,N,đông = 2,20/0,101 = 21,8 giờ → T_ck = 0,5 + 21,8 + 2 + 0,5 = 24,8 giờ
• Kết luận: Vào mùa đông, hệ thống cần chạy gần như liên tục 24 giờ/chu kỳ, gần giới hạn. Cần xem xét giữ ấm nước thải hoặc thiết kế dự phòng.

Phần 5: Hiệu Suất E Trong Các Quá Trình Xử Lý Sinh Học Đặc Thù

5.1 E Trong Nitrat Hóa — Yếu Tố Nhiệt Độ Quyết Định

Nitrat hóa là quá trình có hằng số k nhạy cảm nhất với nhiệt độ trong số các quá trình vi sinh thông thường. Tốc độ tăng trưởng cực đại của Nitrosomonas:

μm,N = 0,47 · e^(0,098·(T-15)) (ngày⁻¹)

Dẫn đến k_N biến thiên mạnh:

Với E_N = 88,9% yêu cầu: thời gian phản ứng t = -ln(0,111)/k_N tỷ lệ nghịch với k_N. Tại 10°C, thời gian phản ứng cần thiết tăng gần 1,5 lần so với 20°C — đây là lý do hệ thống SBR xử lý nitơ tại Hà Nội cần thiết kế HRT “dư” đáng kể so với tính toán ở 20°C.

5.2 E Trong Khử Nitrat — Tầm Quan Trọng Của Nguồn Carbon

Khử nitrat (denitrification) tuân theo bậc phản ứng gần nhất (≈ bậc 0 đến bậc 1) tùy theo nguồn carbon:

Sử dụng carbon nội sinh (nước thải nghèo carbon): Tốc độ chậm, hành vi bậc không → E tuyến tính theo thời gian. Điển hình gặp ở nước thải đô thị Hà Nội và TP.HCM đã qua bể hiếu khí.

Sử dụng methanol bổ sung: Phản ứng bậc nhất theo methanol → E = 1 – e^(-k_MeOH·t). k_MeOH = 0,3-0,7 giờ⁻¹ ở 20°C. Chi phí methanol thường là yếu tố kinh tế quyết định.

Sử dụng carbon nước thải thô (bypass): Tốc độ nhanh hơn nhiều (bậc nhất, k lớn hơn). Đây là giải pháp được ưu tiên trong thiết kế các trạm XLNT đô thị tại Việt Nam để tiết kiệm chi phí hóa chất.

5.3 E Trong Tích Lũy Photpho Sinh Học (EBPR)

Quá trình tích lũy photpho bởi vi khuẩn PAO (Phosphorus Accumulating Organisms) trong giai đoạn hiếu khí của bể gián đoạn tuân theo động học phức tạp hơn. Tuy nhiên, ở giai đoạn đầu khi lượng PHB (polyhidroxybutyrate) nội bào còn đủ, tốc độ hấp thụ P gần bậc nhất theo nồng độ P trong dung dịch:

v_P = k_P · C_P → E_P = 1 – e^(-k_P · t)

k_P = 0,05-0,15 giờ⁻¹ ở 20°C trong điều kiện PAO phát triển tốt.

Để đạt E_P = 73% (loại P từ 15 mg/L xuống 4 mg/L theo QCVN 40:2011): t = -ln(1-0,73)/k_P = -ln(0,27)/k_P = 1,31/k_P

Với k_P = 0,1 giờ⁻¹: t ≈ 13 giờ — đây là thời gian pha hiếu khí cần thiết trong bể SBR xử lý photpho sinh học.

Phần 6: Giới Hạn Thực Tế Của E Và Nồng Độ Đầu Ra Tối Thiểu

6.1 E_max < 100% — Nguyên Nhân Kỹ Thuật

Trong thực tế, E không thể đạt 100% vì nhiều lý do:

Nồng độ cân bằng (equilibrium concentration) S_min: Trong phản ứng nghịch (reversible reactions), có một nồng độ tối thiểu S_min dưới đó phản ứng không tiến xa hơn được. Trong bùn hoạt tính, S_min phụ thuộc vào SRT và hệ số phân rã bùn kp.

Tế bào vi sinh vật tự phân hủy: Quá trình phân rã tế bào (endogenous respiration) tạo ra sản phẩm có COD không thể phân hủy thêm được (non-biodegradable endogenous products).

Hợp chất trơ sinh học: Nước thải luôn chứa một phần COD không phân hủy sinh học (refractory COD). Tỷ lệ này thường 10-30% tổng COD đầu vào trong nước thải công nghiệp.

6.2 Từ E Tính Toán Đến E Đo Đạc Thực Tế

Kỹ sư vận hành cần hiểu rằng E tính từ PT 7-7 là E lý thuyết dựa trên mô hình bậc nhất đơn giản. E thực đo được tại công trình có thể thấp hơn vì:

• Hành vi thủy lực không lý tưởng (dead zones, short-circuiting)
• Biến động tải trọng trong ngày (lưu lượng và nồng độ thay đổi)
• Sự thay đổi hoạt tính vi sinh vật do thay đổi nhiệt độ, pH, chất độc bất thường
• Không đạt SRT thiết kế do sự cố bể lắng hoặc mất bùn

Do đó, trong thực tế thiết kế, thường nhân thêm hệ số an toàn SF = 1,2-1,5 vào thời gian phản ứng tính từ PT 7-7.

Kết Luận

Phương trình nghiệm bậc n tổng quát (PT 7-5) mở rộng khả năng mô hình hóa của kỹ sư XLNT ra ngoài phạm vi của ba bậc nguyên (0, 1, 2), bao gồm các trường hợp bậc phân số thường gặp trong thực tế nước thải hỗn hợp.

Hiệu suất phản ứng E (PT 7-6) là ngôn ngữ chung kết nối hai thế giới: thế giới toán học của phương trình nồng độ C(t) và thế giới thực tiễn của tiêu chuẩn xả thải QCVN. Phương trình PT 7-7 (E = 1 – e^(-kt) cho bậc nhất) cho phép đảo ngược bài toán: từ mục tiêu E yêu cầu, tính trực tiếp thời gian phản ứng cần thiết — nền tảng không thể thiếu của mọi quy trình thiết kế bể gián đoạn và hệ thống SBR tại Việt Nam.