WesterntechVN – Từ thiết bị thí nghiệm liên tục có hồi lưu sinh khối, kỹ sư xây dựng phương trình thăng bằng chất để mô tả biến động sinh khối (PT 6-27).
Áp dụng điều kiện trạng thái ổn định (dX/dt = 0), từng bước rút gọn phương trình để thu được biểu thức 1/θc = Y·vsu/X – kp (PT 6-29). Tiếp theo, từ phương trình phi tuyến Monod tốc độ tiêu thụ cơ chất (PT 6-30, 6-31), tuyến tính hóa theo phương pháp Lineweaver-Burk (PT 6-33): τ·X/(So-S) = Ks/k × (1/S) + 1/k, cho phép xác định k và Ks từ đồ thị đường thẳng mà không cần phương pháp số phức tạp. Bài viết phân tích từng bước toán học, ý nghĩa của độ dốc và đoạn chắn trục, và cách áp dụng thực tế cho nước thải tại Việt Nam.
1. Xây dựng phương trình thăng bằng sinh khối — PT 6-27
1.1 Nguyên lý thăng bằng chất
Bất kỳ hệ phản ứng nào cũng tuân theo nguyên lý bảo toàn khối lượng. Cho bình phản ứng thể tích V, khuấy trộn lý tưởng, có hồi lưu bùn, phương trình thăng bằng sinh khối X viết theo nguyên lý tích lũy = vào – ra + sinh ra:
V · dX/dt = Q·Xo – [Qw·Xr + Qe·Xe] + V·vg’ — (PT 6-27)
Phân tích từng số hạng:
V · dX/dt: tốc độ thay đổi lượng sinh khối trong bình (mg VSS/ngày). Dương nghĩa là X đang tăng, âm nghĩa X đang giảm. Ở trạng thái ổn định, số hạng này bằng 0.
Q·Xo: lượng sinh khối đưa vào bình mỗi ngày từ dòng nạp. Trong hầu hết các ứng dụng, nước thải đầu vào không chứa vi sinh hoạt tính (Xo = 0), nên số hạng này triệt tiêu. Lưu ý: nếu thực nghiệm dùng nước thải thực tế có chứa vi sinh (ví dụ nước thải hòa hợp từ cống chung), cần đo và tính Xo.
[Qw·Xr + Qe·Xe]: lượng sinh khối rời khỏi hệ mỗi ngày theo hai đường: bùn thải (Qw·Xr) và nước ra (Qe·Xe). Đây là lượng “thải ra” kiểm soát tuổi bùn.
V·vg’: lượng sinh khối sinh ra mỗi ngày do quá trình sinh trưởng trong bình. vg’ là tốc độ sinh trưởng thực (net growth rate) tính theo PT 6-10, 6-11:
vg’ = μm·X·S/(Ks+S) – kp·X = -Y·vsu – kp·X (PT 6-10/6-11)
Dấu âm trước vsu vì vsu là tốc độ tiêu thụ cơ chất (âm, vì S giảm), nên -Y·vsu dương (X tăng), trừ kp·X là phân hủy nội sinh (âm đối với sinh khối).
1.2 Ý nghĩa vật lý của từng thành phần
Để hiểu sâu PT 6-27, hãy xem xét “vòng đời” của vi khuẩn trong hệ thí nghiệm:
Vi khuẩn được sinh ra trong bình phản ứng nhờ hấp thu và chuyển hóa cơ chất. Chúng tồn tại trong bình một thời gian, rồi rời đi qua hai con đường: theo nước ra (Qe·Xe) hoặc bị rút ra với bùn thải (Qw·Xr). Ngoài ra, một phần vi khuẩn “chết” tại chỗ do phân hủy nội sinh (-kp·X trong vg’).
PT 6-27 ghi lại đầy đủ tất cả dòng vào, ra và sinh/chết của vi khuẩn trong bình.
2. Áp dụng điều kiện trạng thái ổn định — Thu gọn về PT 6-29
2.1 Ba điều kiện đơn giản hóa
Ở trạng thái ổn định (steady state), áp dụng ba điều kiện:
Điều kiện 1: dX/dt = 0 — Mật độ vi sinh trong bình không thay đổi theo thời gian. Đây là yêu cầu cơ bản để hệ được gọi là “ổn định.”
Điều kiện 2: Xo = 0 — Nước thải đầu vào không chứa vi sinh hoạt tính (vi sinh trong bình đã thích nghi, hoạt tính cao hơn nhiều so với vi sinh từ nước thải tươi không qua acclimation).
Điều kiện 3: Qe = Q — Lưu lượng ra bằng lưu lượng vào (thể tích bình không đổi).
2.2 Rút gọn phương trình
Áp dụng ba điều kiện vào PT 6-27:
0 = 0 – [Qw·Xr + Qe·Xe] + V·vg’
V·vg’ = Qw·Xr + Qe·Xe
Chia cả hai vế cho V·X:
vg’/X = (Qw·Xr + Qe·Xe) / (V·X)
Vế phải chính là 1/θc (từ PT 6-26):
vg’/X = 1/θc
Thay vg’ bằng biểu thức PT 6-11:
(-Y·vsu – kp·X)/X = 1/θc
Sắp xếp lại:
1/θc = Y·(-vsu/X) – kp = Y·(|vsu|/X) – kp (vì vsu < 0, viết |vsu| cho tiện):
1/θc = Y · vsu/X – kp — (PT 6-29)
Lưu ý: Trong PT 6-29, vsu được hiểu là độ lớn tốc độ tiêu thụ cơ chất (giá trị dương), tức |vsu|.
2.3 Ý nghĩa của PT 6-29
Phương trình này là nền tảng của toàn bộ quy trình xác định hằng số động học từ thực nghiệm:
- Vế trái 1/θc: biết từ thực nghiệm (đo Qw, Xr, Qe, Xe, V, X)
- Vế phải Y·vsu/X – kp: chứa hai ẩn Y và kp cần xác định, plus vsu/X biết từ thực nghiệm
Quan hệ 1/θc = Y·(vsu/X) – kp là phương trình đường thẳng với vsu/X là biến số. Vẽ đồ thị 1/θc vs. vsu/X, đường thẳng có độ dốc = Y và tung độ gốc = -kp. Phương pháp xác định Y và kp này sẽ được trình bày chi tiết trong bài tiếp theo.
3. Tốc độ tiêu thụ cơ chất vsu từ thực nghiệm — PT 6-30
3.1 Cân bằng chất cho cơ chất
Tương tự PT 6-27 cho sinh khối, xây dựng phương trình thăng bằng chất cho cơ chất S:
V · dS/dt = Q·So – Qe·S + V·vsu
Ở trạng thái ổn định (dS/dt = 0), với Qe = Q:
0 = Q·So – Q·S + V·vsu
V·vsu = -Q·(So – S)
vsu = -Q·(So – S)/V = -(So – S)/τ — (PT 6-30)
Dấu âm xác nhận vsu < 0 (cơ chất giảm, S < So).
3.2 Tính toán vsu thực tế
Từ PT 6-30, kỹ sư chỉ cần đo ba thông số dễ đo: Q (lưu lượng), So (COD đầu vào), S (COD đầu ra), và biết V (thể tích bình):
vsu = -(So – S) × Q / V = -(So – S) / τ
Ví dụ thực tế:
Bình V = 10 L, Q = 5 L/ngày (τ = 2 ngày), So = 300 mg COD/L, S = 50 mg COD/L:
vsu = -(300 – 50)/2 = -125 mg COD/(L·ngày)
Hay |vsu| = 125 mg COD/(L·ngày) — tốc độ tiêu thụ COD là 125 mg/L mỗi ngày.
4. Kết hợp Monod với thực nghiệm — PT 6-31 và 6-32
4.1 Từ mô hình Monod đến phương trình thực nghiệm
Từ mô hình Monod (PT 6-8): vsu = -k·X·S/(Ks+S)
Kết hợp với PT 6-30: -k·X·S/(Ks+S) = -(So-S)/τ
k·X·S/(Ks+S) = (So-S)/τ — (PT 6-31)
Chia cả hai vế cho X:
k·S/(Ks+S) = (So-S)/(τ·X) — (PT 6-32)
Đây là phương trình liên kết mô hình Monod lý thuyết (vế trái, chứa k và Ks chưa biết) với số liệu thực nghiệm đo được (vế phải, chứa So, S, τ, X).
4.2 Vấn đề phi tuyến của PT 6-32
PT 6-32 có hai ẩn k và Ks nằm trong hàm phi tuyến k·S/(Ks+S). Không thể giải trực tiếp. Cần phương pháp biến đổi để tuyến tính hóa — đây là bước then chốt của phương pháp Lineweaver-Burk.
5. Tuyến tính hóa Lineweaver-Burk — PT 6-33
5.1 Biến đổi toán học
Lấy nghịch đảo cả hai vế của PT 6-32:
(Ks+S)/(k·S) = τ·X/(So-S)
Ks/(k·S) + 1/k = τ·X/(So-S)
Thay 1/S = u (biến mới):
(Ks/k)·u + 1/k = τ·X/(So-S) — (PT 6-33)
Viết lại dưới dạng hàm tuyến tính:
τ·X/(So-S) = (Ks/k)·(1/S) + 1/k
Đây là phương trình đường thẳng y = a·x + b với:
- y = τ·X/(So-S) (tính từ thực nghiệm)
- x = 1/S (tính từ thực nghiệm)
- a = độ dốc = Ks/k (cần xác định)
- b = đoạn chắn trục tung (intercept) = 1/k (cần xác định)
5.2 Sơ đồ xác định k và Ks từ đồ thị
Bước 1: Thu thập số liệu thực nghiệm tại ít nhất 5 mức θc khác nhau. Với mỗi mức θc, ghi lại S, τ, X, So.
Bước 2: Tính hai cột số liệu cho đồ thị:
- Trục hoành (x) = 1/S (L/mg)
- Trục tung (y) = τ·X/(So-S) (ngày·mg VSS/L / mg COD/L = ngày·mg VSS/mg COD)
Bước 3: Vẽ đồ thị x-y và vẽ đường hồi quy tuyến tính (có thể bằng Excel hoặc tay).
Bước 4: Đọc thông số đường thẳng:
- Intercept b = 1/k → k = 1/b
- Slope a = Ks/k → Ks = a/b = a·k
5.3 Ví dụ tính toán minh họa
Giả sử số liệu thực nghiệm từ 5 mức θc như sau (V = 10 L, Q = 5 L/ngày, So = 300 mg COD/L):
| θc (ngày) | S (mg COD/L) | X (mg VSS/L) | τ (ngày) |
|---|---|---|---|
| 2 | 85 | 900 | 2 |
| 4 | 55 | 1.500 | 2 |
| 6 | 40 | 2.000 | 2 |
| 8 | 32 | 2.350 | 2 |
| 10 | 28 | 2.600 | 2 |
Lưu ý: τ = V/Q = 10/5 = 2 ngày, không đổi vì V và Q cố định. θc thay đổi bằng cách điều chỉnh Qw.
Tính toán x = 1/S và y = τ·X/(So-S):
| θc | S | 1/S = x (L/mg) | τ·X/(So-S) = y (ngày·mg VSS/mg COD) |
|---|---|---|---|
| 2 | 85 | 0,01176 | 2×900/(300-85) = 1800/215 = 8,37 |
| 4 | 55 | 0,01818 | 2×1500/(300-55) = 3000/245 = 12,24 |
| 6 | 40 | 0,02500 | 2×2000/(300-40) = 4000/260 = 15,38 |
| 8 | 32 | 0,03125 | 2×2350/(300-32) = 4700/268 = 17,54 |
| 10 | 28 | 0,03571 | 2×2600/(300-28) = 5200/272 = 19,12 |
Hồi quy tuyến tính y = a·x + b từ 5 điểm trên (có thể tính bằng Excel):
Giả sử kết quả hồi quy cho: a = 350, b = 4,3 (các giá trị ví dụ)
Từ đó:
- k = 1/b = 1/4,3 = 0,233 mg COD/(mg VSS·ngày) = 0,233 d⁻¹
- Ks = a·k = 350 × 0,233 = 81,6 mg COD/L
5.4 Kiểm tra tính hợp lệ của kết quả
Sau khi tính k và Ks, cần kiểm tra:
- R² của đường hồi quy ≥ 0,95 (độ tương quan cao). Nếu R² < 0,90, số liệu có vấn đề (sai số đo lớn hoặc chưa đạt trạng thái ổn định).
- Ks > 0 và k > 0. Nếu intercept b < 0 hoặc slope a < 0, kết quả vô nghĩa — kiểm tra lại số liệu.
- So sánh với giá trị tài liệu. k thường 0,1-0,4 mg COD/(mg VSS·ngày) cho nước thải sinh hoạt. Ks thường 30-120 mg COD/L. Nếu kết quả thực nghiệm lệch quá xa, kiểm tra lại quy trình thí nghiệm.
6. Sai số hệ thống trong phương pháp Lineweaver-Burk
6.1 Vấn đề phân bố sai số
Phương pháp Lineweaver-Burk (nghịch đảo kép — double reciprocal) là kỹ thuật tuyến tính hóa cổ điển, phổ biến từ thập niên 1930. Tuy nhiên, có nhược điểm thống kê quan trọng:
Khi lấy nghịch đảo 1/S và τ·X/(So-S), các điểm có S nhỏ (nồng độ cơ chất đầu ra thấp) sẽ có 1/S lớn và bị khuếch đại sai số. Đặc biệt:
- Điểm θc lớn (S thấp, 1/S lớn): sai số đo S nhỏ thôi sẽ gây biến động lớn trên trục x
- Điểm θc nhỏ (S cao, 1/S nhỏ): ít nhạy cảm với sai số đo
Hệ quả: đường hồi quy tuyến tính thường bị kéo về phía các điểm θc cao (S thấp), có thể làm sai lệch ước tính k và Ks.
6.2 Biện pháp giảm thiểu sai số Lineweaver-Burk
Phân bổ S đồng đều trên trục 1/S. Thay vì chọn θc đều nhau (2, 4, 6, 8, 10 ngày), chọn θc sao cho 1/S đều nhau. Ví dụ: nếu dự kiến S trong khoảng 20-100 mg/L, thì 1/S từ 0,01 đến 0,05 L/mg, chọn 5 điểm: 1/S = 0,01; 0,02; 0,03; 0,04; 0,05, từ đó ngược lại tính θc cần thiết.
Dùng trọng số hồi quy (weighted regression). Gán trọng số cao hơn cho điểm có S cao (sai số tương đối nhỏ hơn), giảm ảnh hưởng của điểm S thấp.
Tăng số điểm dữ liệu. Dùng 7-10 điểm thay vì 5, giảm ảnh hưởng của điểm bất thường.
6.3 So sánh Lineweaver-Burk với các phương pháp tuyến tính hóa khác
Ngoài Lineweaver-Burk (1/v vs 1/S), còn có các phương pháp tuyến tính hóa khác cho mô hình Monod:
Eadie-Hofstee plot: v/S vs v → độ dốc = -k, tung độ gốc = μm. Phân bố sai số tốt hơn Lineweaver-Burk.
Hanes-Woolf plot: S/v vs S → độ dốc = 1/k, tung độ gốc = Ks/k. Phân bố sai số tốt nhất trong các phương pháp tuyến tính hóa.
Trong bối cảnh PT 6-33, đồ thị τ·X/(So-S) vs 1/S tương đương Lineweaver-Burk. Kỹ sư có thể thử thêm biểu đồ S·τ·X/(So-S)^… để so sánh và kiểm tra độ nhất quán.
7. Ứng dụng tính toán cho nước thải tại Việt Nam
7.1 Nước thải khu công nghiệp chế biến thực phẩm
Giả sử kết quả thực nghiệm tại một nhà máy chế biến thực phẩm tại Bình Dương:
| θc (ngày) | S (mg COD/L) | X (mg VSS/L) |
|---|---|---|
| 3 | 160 | 1.200 |
| 5 | 95 | 1.900 |
| 8 | 62 | 2.800 |
| 12 | 45 | 3.600 |
| 15 | 38 | 4.100 |
τ = 4 ngày, So = 900 mg COD/L.
Tính y = τ·X/(So-S) và x = 1/S:
| θc | y = 4×X/(900-S) | x = 1/S |
|---|---|---|
| 3 | 4×1200/740 = 6,49 | 0,00625 |
| 5 | 4×1900/805 = 9,44 | 0,01053 |
| 8 | 4×2800/838 = 13,37 | 0,01613 |
| 12 | 4×3600/855 = 16,84 | 0,02222 |
| 15 | 4×4100/862 = 19,03 | 0,02632 |
Hồi quy tuyến tính (Excel): y = 496 × x + 3,31, R² = 0,998
Từ đó:
- k = 1/3,31 = 0,302 mg COD/(mg VSS·ngày)
- Ks = 496/k = 496 × 0,302 = 149,8 mg COD/L ≈ 150 mg COD/L
Nhận xét: Ks = 150 mg COD/L — cao hơn nhiều so với nước thải sinh hoạt (Ks ≈ 60 mg COD/L), phản ánh sự phức tạp của cơ chất trong nước thải thực phẩm (protein phân tử lớn, lipid, chuỗi carbohydrate). Nếu kỹ sư dùng Ks = 60 mg COD/L từ tài liệu cho thiết kế này, sẽ đánh giá thấp hơn 2,5 lần khả năng xử lý thực tế và thiết kế bể không đủ thể tích.
7.2 Kiểm tra nhất quán với phương trình Monod
Sau khi có k và Ks, kiểm tra lại: tại θc = 5 ngày, S thực nghiệm = 95 mg COD/L, X = 1.900 mg/L.
Từ PT 6-32: (So-S)/(τ·X) = k·S/(Ks+S)
Vế phải: k·S/(Ks+S) = 0,302 × 95 / (150 + 95) = 28,69/245 = 0,117
Vế trái: (900-95)/(4×1900) = 805/7600 = 0,106
Sai số: (0,117 – 0,106)/0,106 × 100% = 10,4% — chấp nhận được (< 15%)
8. Quy trình thực hành xác định k và Ks từ A đến Z
8.1 Tóm tắt quy trình
- Lắp hệ thí nghiệm liên tục có hồi lưu bùn (V, Q, hệ bơm)
- Cho hệ khởi động với bùn seed từ trạm XLNT hiện hữu, chạy 2-3 tuần để ổn định
- Đặt θc = θc,1 (giá trị nhỏ nhất, ví dụ 2 ngày) bằng cách điều chỉnh Qw
- Vận hành đến trạng thái ổn định (3-5 × θc,1 ngày)
- Đo Q, So, S, X, Xe, Xr, Qw — lấy trung bình 3-5 ngày liên tiếp
- Tăng θc lên giá trị tiếp theo, lặp lại bước 4-5
- Thu thập đủ 5+ điểm dữ liệu tại các mức θc khác nhau
- Tính x = 1/S, y = τ·X/(So-S) cho mỗi điểm
- Hồi quy tuyến tính y = a·x + b
- Tính k = 1/b, Ks = a·k
- Kiểm tra R² ≥ 0,95 và so sánh với giá trị tham khảo
8.2 Lưu ý đặc thù cho thực hành tại Việt Nam
Nước thải sinh hoạt: Thành phần thay đổi mạnh theo giờ trong ngày (sáng, trưa, tối). Cần lấy mẫu tổng hợp 24 giờ (composite sample) thay vì mẫu đơn, hoặc dùng bình chứa đều hóa (equalization tank) trước khi bơm vào thiết bị thí nghiệm.
Mùa mưa/mùa khô: So mùa khô thường cao hơn So mùa mưa do pha loãng từ nước mưa chảy tràn. Nên thực nghiệm cả hai mùa nếu đây là nước thải sinh hoạt đô thị, hoặc cố định So bằng cách pha thêm nước thải cô đặc (concentrated waste).
Tiêu chuẩn phân tích COD: Đảm bảo dùng phương pháp COD dichromate (phương pháp chuẩn SMEWW 5220C) để kết quả nhất quán và so sánh được với tài liệu quốc tế.
9. Kết luận
Chuỗi suy luận toán học từ PT 6-27 (thăng bằng sinh khối) → PT 6-29 (1/θc = Y·vsu/X – kp) → PT 6-30 (vsu = -(So-S)/τ) → PT 6-32 (kết hợp Monod với thực nghiệm) → PT 6-33 (tuyến tính hóa Lineweaver-Burk) là một chuỗi logic chặt chẽ, cho phép xác định hai trong bốn hằng số động học quan trọng nhất là k và Ks từ số liệu thực nghiệm đơn giản.
Điểm mấu chốt: PT 6-33 biến phương trình phi tuyến Monod thành phương trình đường thẳng y = a·x + b, từ đó chỉ cần công cụ hồi quy tuyến tính cơ bản (Excel) là đủ để tính k và Ks. Tuy nhiên, cần hiểu rõ giới hạn của phương pháp Lineweaver-Burk (sai số khuếch đại tại S thấp) và áp dụng biện pháp kiểm tra chất lượng phù hợp để đảm bảo kết quả đáng tin cậy.
