1. Cần giữ lại đại lượng nào khi pha loãng
Một mẫu nước sau khi hòa thêm nước sạch có thể cho nồng độ thấp hơn dù lượng chất đang theo dòng chảy không giảm. Vì vậy, khi đánh giá một phương án gom dòng trong nhà máy, cần ghi đồng thời lưu lượng, nồng độ và tải lượng. Chỉ một con số mg/L không mô tả được tổng lượng chất đi qua ranh giới khảo sát.
Tài liệu NPDES của EPA trình bày cân bằng lưu lượng nhân nồng độ cho mô hình trộn hoàn toàn ở trạng thái ổn định; tài liệu cũng phân biệt trường hợp chất bảo toàn với trường hợp cần mô hình phản ứng hoặc trộn không hoàn toàn.[1] Bài này dùng nguyên lý cân bằng đó cho các số liệu tự xây dựng. Mọi nồng độ đích đều là giả định tính toán, không phải giới hạn xả thải và không phải đề xuất dùng pha loãng thay cho xử lý.
2. Đặt giả thiết và kiểm tra đơn vị
Xét một chất hòa tan bảo toàn trong khoảng thời gian khảo sát: không phản ứng, không bay hơi, không hấp phụ và không kết tủa. Hai dòng nước loãng có thể tích cộng được, cùng được đo trong một cơ sở thời gian. Gọi Q là lưu lượng m³/h, C là nồng độ mg/L và L là tải lượng kg/h.
L = Q C / 1.000
Q_tổng = Q_1 + Q_2
L_tổng = L_1 + L_2
C_trộn = (Q_1 C_1 + Q_2 C_2)/(Q_1 + Q_2)
Hệ số 1/1.000 xuất hiện vì một m³ bằng 1.000 L và một kg bằng một triệu mg. Khi Q dùng m³/ngày, kết quả cùng công thức sẽ là kg/ngày. Việc đổi nhãn giờ sang ngày mà không đổi số làm tải lượng sai 24 lần. Cũng cần phân biệt L, ký hiệu tải lượng trong bài, với chữ L biểu diễn đơn vị lít.
3. Hai dòng cùng chảy, một dòng vẫn có nồng độ nền
Giả định dòng thứ nhất Q_1 = 40 m³/h, C_1 = 600 mg/L. Dòng thứ hai Q_2 = 60 m³/h, C_2 = 20 mg/L. Đây là hai dòng đồng thời, ổn định, trộn hoàn toàn trước vị trí lấy mẫu.
L_1 = 40 × 600/1.000 = 24,0 kg/h
L_2 = 60 × 20/1.000 = 1,2 kg/h
L_tổng = 25,2 kg/h; Q_tổng = 100 m³/h
C_trộn = 252 mg/L
Nồng độ giảm từ 600 xuống 252 mg/L, nhưng tải lượng sau điểm nhập dòng lớn hơn tải lượng của riêng dòng thứ nhất. Phần tăng 1,2 kg/h đến từ nồng độ nền của dòng thứ hai. Nếu bỏ qua C_2, người tính sẽ dự báo 240 mg/L và thiếu đúng phần tải lượng này.
Nếu thực sự C_2 = 0, nồng độ 240 mg/L đi cùng tải lượng 24 kg/h. Khi đó không có một kg chất nào bị loại bỏ trong mô hình. Lấy 1 − 240/600 = 60% rồi gọi là “hiệu suất xử lý” sẽ đánh đồng sự tăng lượng nước với sự giảm lượng chất.
4. Giải ngược lưu lượng và đọc giới hạn
Chọn C_đích = 300 mg/L chỉ để khảo sát đại số. Với C_1 = 600 và C_2 = 20 mg/L, lưu lượng dòng thứ hai cần để đạt đúng nồng độ đó trong mô hình là:
Q_2 = Q_1 (C_1 − C_đích)/(C_đích − C_2)
Q_2 = 40 × 300/280 = 42,8571 m³/h
Kết quả hữu hạn, dương đòi hỏi nồng độ đích nằm giữa hai nồng độ đầu vào. Khi C_đích tiến gần 20 mg/L từ phía trên, Q_2 tăng không giới hạn. Đích nhỏ hơn 20 mg/L không thể đạt bằng cách trộn hai dòng này. Đây là giới hạn của phương trình, trước cả khi xét khả năng cấp nước, kích thước bể hoặc quy định áp dụng.
Tại Q_2 = 60 m³/h, đạo hàm của nồng độ theo lưu lượng thêm bằng (C_2 − C_trộn)/Q_tổng = −2,32 (mg/L)/(m³/h). Giá trị này chỉ là độ nhạy cục bộ. Không thể nhân nó với một thay đổi lưu lượng rất lớn rồi coi quan hệ tuyến tính trên toàn miền.
5. Bể có tồn chứa sẽ không đổi nồng độ tức thời
Giờ đưa hỗn hợp 252 mg/L vào một bể khuấy trộn lý tưởng thể tích không đổi V = 250 m³, lưu lượng vào và ra cùng 100 m³/h. Ban đầu bể chứa nước ở C_0 = 600 mg/L. Khối lượng trong bể là M = VC/1.000 kg. Bài tập cân bằng của TU Delft cũng đặt bài toán bể trộn có lưu lượng vào bằng ra và nồng độ đầu khác nồng độ cấp vào.[2] Với bộ số riêng ở đây, ta xây dựng nghiệm từ tích lũy bằng lượng vào trừ lượng ra.
Với chất bảo toàn trong ví dụ, hai số hạng phản ứng bằng không. Vì đầu ra có nồng độ bằng nồng độ đồng nhất trong bể, phương trình và nghiệm là:
V dC/dt = Q (C_vào − C)
τ = V/Q = 2,5 h
C(t) = 252 + 348 exp(−t/2,5)
Sau 2 h, C = 408,3665 mg/L; sau 4 h, C = 322,2600 mg/L. Thời gian để giảm đến 300 mg/L bằng 2,5 ln(348/48) = 4,9525 h. Bể không đạt trạng thái mới sau đúng một thời gian lưu; tại t = τ, phần chênh nồng độ còn lại bằng exp(−1) của chênh lệch ban đầu.
6. Kiểm toán khối lượng qua giai đoạn quá độ
Ban đầu M_0 = 150 kg. Trong 2 h, lượng chất đi vào bằng 25,2 × 2 = 50,4 kg. Cuối khoảng đó, M_2 = 250 × 408,3665/1.000 = 102,0916 kg. Từ cân bằng tích phân, tổng lượng đi ra là:
M_ra = M_0 + M_vào − M_2
M_ra = 150 + 50,4 − 102,0916 = 98,3084 kg
Lượng trong bể giảm 47,9084 kg, nhưng đây là lượng được vận chuyển ra ngoài nhiều hơn lượng đi vào. Mô hình không có cơ chế tiêu hủy. Cũng không được lấy nồng độ cuối 408,3665 mg/L nhân toàn bộ 200 m³ nước ra trong 2 h: đầu ra đã biến đổi liên tục từ 600 mg/L. Phải tích phân hoặc dùng cân bằng tồn chứa như trên.
Một kiểm tra khác là lượng vào 50,4 kg cộng lượng giảm tồn chứa 47,9084 kg phải bằng lượng ra 98,3084 kg. Cách đối chiếu này giúp phát hiện báo cáo bỏ quên nước và chất còn nằm trong bể vào đầu hoặc cuối kỳ.
7. Dữ liệu biến thiên theo giờ cần ghép đúng cặp
Xét một trường hợp riêng, không phải hai dòng đồng thời ở mục 3. Một đường ống vận hành 1 h ở 20 m³/h, 900 mg/L, rồi 1 h ở 80 m³/h, 100 mg/L. Tổng thể tích là 100 m³ và tổng lượng chất là 18 + 8 = 26 kg. Nồng độ trung bình theo thể tích bằng 260 mg/L.
Trung bình cộng của hai nồng độ là 500 mg/L. Nhân số đó với lưu lượng trung bình 50 m³/h và 2 h sẽ cho 50 kg, cao hơn 26 kg thực sự trong ví dụ. Sai lệch xuất hiện vì thời kỳ lưu lượng cao lại có nồng độ thấp. Cặp Q và C phải giữ cùng dấu thời gian; trung bình của tích không mặc nhiên bằng tích của hai trung bình.
8. Điểm lấy mẫu và mức độ trộn
Nghiên cứu của USGS về chất đánh dấu sông nêu rằng, khi mặt cắt chưa trộn hoàn toàn, nồng độ dùng với tổng lưu lượng phải là trung bình mặt cắt có trọng số theo lưu lượng.[3] Một mẫu lấy đúng vùng dòng đậm đặc hoặc vùng nước nền sẽ không tự đại diện cho toàn mặt cắt.
Trong bể thực, mô hình khuấy lý tưởng còn cần đối chiếu với đường đi ngắn, vùng ít trao đổi và thể tích hoạt động. Thể tích hình học 250 m³ không chứng minh toàn bộ thể tích tham gia trộn như giả định. Nếu có phản ứng hoặc chuyển chất sang pha rắn, cần bổ sung số hạng tương ứng và đo các dòng liên quan.
Một bảng tính có thể kiểm tra được nên giữ nguyên dữ liệu Q, C, khoảng thời gian, thể tích tồn đầu–cuối, vị trí mẫu và cơ sở đổi đơn vị. Kết quả cần báo cáo cả kg/h hoặc kg/kỳ cùng mg/L. Nhờ vậy người đọc biết nồng độ thay đổi do lượng nước, do vận chuyển tồn chứa hay do một cơ chế loại bỏ đã được xác định.
Nguồn tham khảo
[1] NPDES Permit Writers’ Manual, Chapter 6 — Simple mass-balance equation — U.S. Environmental Protection Agency, September 2010. Truy cập ngày 09/11/2026.
[2] 1.2 The mass balance: examples — Diluting within a stirred tank reactor; unsteady state mass balance — TU Delft OpenCourseWare, The Basics of Transport Phenomena. Truy cập ngày 09/11/2026.
[3] Prediction of Traveltime and Longitudinal Dispersion in Rivers and Streams, WRI 96-4013 — Harvey E. Jobson, U.S. Geological Survey, 1996. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Một mẫu nước sau khi hòa thêm nước sạch có thể cho nồng độ thấp hơn dù lượng chất đang theo dòng chảy không giảm. Vì vậy, khi đánh giá một phương án gom dòng trong nhà máy, cần ghi đồng thời lưu lượng, nồng độ và tải lượng. Chỉ một con số mg/L không mô tả được tổng lượng chất đi qua ranh giới khảo sát.
Tài liệu NPDES của EPA trình bày cân bằng lưu lượng nhân nồng độ cho mô hình trộn hoàn toàn ở trạng thái ổn định; tài liệu cũng phân biệt trường hợp chất bảo toàn với trường hợp cần mô hình phản ứng hoặc trộn không hoàn toàn.[1] Bài này dùng nguyên lý cân bằng đó cho các số liệu tự xây dựng. Mọi nồng độ đích đều là giả định tính toán, không phải giới hạn xả thải và không phải đề xuất dùng pha loãng thay cho xử lý.
2. Đặt giả thiết và kiểm tra đơn vị
Xét một chất hòa tan bảo toàn trong khoảng thời gian khảo sát: không phản ứng, không bay hơi, không hấp phụ và không kết tủa. Hai dòng nước loãng có thể tích cộng được, cùng được đo trong một cơ sở thời gian. Gọi Q là lưu lượng m³/h, C là nồng độ mg/L và L là tải lượng kg/h.
L = Q C / 1.000
Q_tổng = Q_1 + Q_2
L_tổng = L_1 + L_2
C_trộn = (Q_1 C_1 + Q_2 C_2)/(Q_1 + Q_2)
Hệ số 1/1.000 xuất hiện vì một m³ bằng 1.000 L và một kg bằng một triệu mg. Khi Q dùng m³/ngày, kết quả cùng công thức sẽ là kg/ngày. Việc đổi nhãn giờ sang ngày mà không đổi số làm tải lượng sai 24 lần. Cũng cần phân biệt L, ký hiệu tải lượng trong bài, với chữ L biểu diễn đơn vị lít.
3. Hai dòng cùng chảy, một dòng vẫn có nồng độ nền
Giả định dòng thứ nhất Q_1 = 40 m³/h, C_1 = 600 mg/L. Dòng thứ hai Q_2 = 60 m³/h, C_2 = 20 mg/L. Đây là hai dòng đồng thời, ổn định, trộn hoàn toàn trước vị trí lấy mẫu.
L_1 = 40 × 600/1.000 = 24,0 kg/h
L_2 = 60 × 20/1.000 = 1,2 kg/h
L_tổng = 25,2 kg/h; Q_tổng = 100 m³/h
C_trộn = 252 mg/L
Nồng độ giảm từ 600 xuống 252 mg/L, nhưng tải lượng sau điểm nhập dòng lớn hơn tải lượng của riêng dòng thứ nhất. Phần tăng 1,2 kg/h đến từ nồng độ nền của dòng thứ hai. Nếu bỏ qua C_2, người tính sẽ dự báo 240 mg/L và thiếu đúng phần tải lượng này.
Nếu thực sự C_2 = 0, nồng độ 240 mg/L đi cùng tải lượng 24 kg/h. Khi đó không có một kg chất nào bị loại bỏ trong mô hình. Lấy 1 − 240/600 = 60% rồi gọi là “hiệu suất xử lý” sẽ đánh đồng sự tăng lượng nước với sự giảm lượng chất.
4. Giải ngược lưu lượng và đọc giới hạn
Chọn C_đích = 300 mg/L chỉ để khảo sát đại số. Với C_1 = 600 và C_2 = 20 mg/L, lưu lượng dòng thứ hai cần để đạt đúng nồng độ đó trong mô hình là:
Q_2 = Q_1 (C_1 − C_đích)/(C_đích − C_2)
Q_2 = 40 × 300/280 = 42,8571 m³/h
Kết quả hữu hạn, dương đòi hỏi nồng độ đích nằm giữa hai nồng độ đầu vào. Khi C_đích tiến gần 20 mg/L từ phía trên, Q_2 tăng không giới hạn. Đích nhỏ hơn 20 mg/L không thể đạt bằng cách trộn hai dòng này. Đây là giới hạn của phương trình, trước cả khi xét khả năng cấp nước, kích thước bể hoặc quy định áp dụng.
Tại Q_2 = 60 m³/h, đạo hàm của nồng độ theo lưu lượng thêm bằng (C_2 − C_trộn)/Q_tổng = −2,32 (mg/L)/(m³/h). Giá trị này chỉ là độ nhạy cục bộ. Không thể nhân nó với một thay đổi lưu lượng rất lớn rồi coi quan hệ tuyến tính trên toàn miền.
5. Bể có tồn chứa sẽ không đổi nồng độ tức thời
Giờ đưa hỗn hợp 252 mg/L vào một bể khuấy trộn lý tưởng thể tích không đổi V = 250 m³, lưu lượng vào và ra cùng 100 m³/h. Ban đầu bể chứa nước ở C_0 = 600 mg/L. Khối lượng trong bể là M = VC/1.000 kg. Bài tập cân bằng của TU Delft cũng đặt bài toán bể trộn có lưu lượng vào bằng ra và nồng độ đầu khác nồng độ cấp vào.[2] Với bộ số riêng ở đây, ta xây dựng nghiệm từ tích lũy bằng lượng vào trừ lượng ra.
Với chất bảo toàn trong ví dụ, hai số hạng phản ứng bằng không. Vì đầu ra có nồng độ bằng nồng độ đồng nhất trong bể, phương trình và nghiệm là:
V dC/dt = Q (C_vào − C)
τ = V/Q = 2,5 h
C(t) = 252 + 348 exp(−t/2,5)
Sau 2 h, C = 408,3665 mg/L; sau 4 h, C = 322,2600 mg/L. Thời gian để giảm đến 300 mg/L bằng 2,5 ln(348/48) = 4,9525 h. Bể không đạt trạng thái mới sau đúng một thời gian lưu; tại t = τ, phần chênh nồng độ còn lại bằng exp(−1) của chênh lệch ban đầu.
6. Kiểm toán khối lượng qua giai đoạn quá độ
Ban đầu M_0 = 150 kg. Trong 2 h, lượng chất đi vào bằng 25,2 × 2 = 50,4 kg. Cuối khoảng đó, M_2 = 250 × 408,3665/1.000 = 102,0916 kg. Từ cân bằng tích phân, tổng lượng đi ra là:
M_ra = M_0 + M_vào − M_2
M_ra = 150 + 50,4 − 102,0916 = 98,3084 kg
Lượng trong bể giảm 47,9084 kg, nhưng đây là lượng được vận chuyển ra ngoài nhiều hơn lượng đi vào. Mô hình không có cơ chế tiêu hủy. Cũng không được lấy nồng độ cuối 408,3665 mg/L nhân toàn bộ 200 m³ nước ra trong 2 h: đầu ra đã biến đổi liên tục từ 600 mg/L. Phải tích phân hoặc dùng cân bằng tồn chứa như trên.
Một kiểm tra khác là lượng vào 50,4 kg cộng lượng giảm tồn chứa 47,9084 kg phải bằng lượng ra 98,3084 kg. Cách đối chiếu này giúp phát hiện báo cáo bỏ quên nước và chất còn nằm trong bể vào đầu hoặc cuối kỳ.
7. Dữ liệu biến thiên theo giờ cần ghép đúng cặp
Xét một trường hợp riêng, không phải hai dòng đồng thời ở mục 3. Một đường ống vận hành 1 h ở 20 m³/h, 900 mg/L, rồi 1 h ở 80 m³/h, 100 mg/L. Tổng thể tích là 100 m³ và tổng lượng chất là 18 + 8 = 26 kg. Nồng độ trung bình theo thể tích bằng 260 mg/L.
Trung bình cộng của hai nồng độ là 500 mg/L. Nhân số đó với lưu lượng trung bình 50 m³/h và 2 h sẽ cho 50 kg, cao hơn 26 kg thực sự trong ví dụ. Sai lệch xuất hiện vì thời kỳ lưu lượng cao lại có nồng độ thấp. Cặp Q và C phải giữ cùng dấu thời gian; trung bình của tích không mặc nhiên bằng tích của hai trung bình.
8. Điểm lấy mẫu và mức độ trộn
Nghiên cứu của USGS về chất đánh dấu sông nêu rằng, khi mặt cắt chưa trộn hoàn toàn, nồng độ dùng với tổng lưu lượng phải là trung bình mặt cắt có trọng số theo lưu lượng.[3] Một mẫu lấy đúng vùng dòng đậm đặc hoặc vùng nước nền sẽ không tự đại diện cho toàn mặt cắt.
Trong bể thực, mô hình khuấy lý tưởng còn cần đối chiếu với đường đi ngắn, vùng ít trao đổi và thể tích hoạt động. Thể tích hình học 250 m³ không chứng minh toàn bộ thể tích tham gia trộn như giả định. Nếu có phản ứng hoặc chuyển chất sang pha rắn, cần bổ sung số hạng tương ứng và đo các dòng liên quan.
Một bảng tính có thể kiểm tra được nên giữ nguyên dữ liệu Q, C, khoảng thời gian, thể tích tồn đầu–cuối, vị trí mẫu và cơ sở đổi đơn vị. Kết quả cần báo cáo cả kg/h hoặc kg/kỳ cùng mg/L. Nhờ vậy người đọc biết nồng độ thay đổi do lượng nước, do vận chuyển tồn chứa hay do một cơ chế loại bỏ đã được xác định.
Nguồn tham khảo
[1] NPDES Permit Writers’ Manual, Chapter 6 — Simple mass-balance equation — U.S. Environmental Protection Agency, September 2010. Truy cập ngày 09/11/2026.
[2] 1.2 The mass balance: examples — Diluting within a stirred tank reactor; unsteady state mass balance — TU Delft OpenCourseWare, The Basics of Transport Phenomena. Truy cập ngày 09/11/2026.
[3] Prediction of Traveltime and Longitudinal Dispersion in Rivers and Streams, WRI 96-4013 — Harvey E. Jobson, U.S. Geological Survey, 1996. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
