1. Một phép chia cần hai định nghĩa rõ ràng
Thời gian lưu thường được ghi bằng thể tích chia lưu lượng thể tích. Nhưng thể tích có thể là toàn vỏ thiết bị, phần chứa chất lỏng, khoảng rỗng giữa hạt xúc tác hoặc thể tích thực tham gia dòng chảy. Lưu lượng lại có thể được báo ở nhiệt độ và áp suất thực, hoặc ở điều kiện quy chiếu. Dùng hai số không cùng cơ sở sẽ tạo ra một “thời gian” không đại diện cho chuyển động của pha đang xét.
Tài liệu IAEA phân biệt thời gian trung bình suy từ thể tích và lưu lượng với giá trị thực nghiệm thu từ phân bố thời gian lưu, đồng thời xem xét thể tích chết và vùng trao đổi chậm.[1] Trong bài này, các ví dụ đều là giả định để kiểm tra cơ sở tính, không phải kết luận về một thiết bị đang vận hành.
2. Từ thể tích vỏ tới lượng pha thực có mặt
Xét một bồn có thể tích hình học 10 m³. Mức chất lỏng chiếm vùng hình học 7 m³; vật rắn và chi tiết nằm ngập chiếm 1 m³. Nếu bỏ qua lỗ rỗng chứa lỏng bên trong vật rắn, lượng pha lỏng có mặt là 6 m³. Với dòng lỏng không nén được 2 m³/h, thời gian danh nghĩa theo pha lỏng bằng:
tau_l = V_l/Q_l = 6/2 = 3 h
Dùng toàn vỏ sẽ cho 5 h, cao hơn 66,666667%. Dùng vùng dưới mức mà quên thể tích chiếm chỗ sẽ cho 3,5 h. Cả hai lỗi xuất phát từ cách định nghĩa thể tích, không phải phép làm tròn.
Giả định thêm 0,6 m³ lỏng bị cô lập hoàn toàn, không trao đổi với phần còn lại trong thí nghiệm. Thể tích tham gia là 5,4 m³ và tỷ số tương ứng 2,7 h. Nếu vùng đó chỉ trao đổi rất chậm, cách thay đơn giản 6 thành 5,4 m³ có thể không mô tả đầy đủ phần đuôi của đáp ứng. Cần phân biệt “không vào được” với “vào được nhưng trở ra chậm”.
3. Khí phải dùng lưu lượng tại đúng điều kiện
Xét một thiết bị khí độc lập có khoảng khí 6 m³. Lưu lượng quy chiếu là 1.000 m³/h tại 273,15 K, 101.325 Pa và Z_n = 1. Tại thiết bị, giả định 450 K, áp suất tuyệt đối 2.000.000 Pa, hệ số nén Z_a = 0,9. Cùng một dòng mol cho:
Q_a = Q_n (Z_a/Z_n)(T_a/T_n)(p_n/p_a)
Q_a = 75,117381 m³/h
tau_a = 6/Q_a = 4,792499 min
Nếu chia trực tiếp 6 m³ cho 1.000 m³/h, kết quả chỉ 0,36 phút. Chênh lệch không phải khí thực đang lưu “lâu hơn thiết kế”; nó do đã ghép thể tích thực của thiết bị với lưu lượng ở điều kiện khác.
Mô hình quy đổi trên giả định không có thay đổi mol do phản ứng hoặc ngưng tụ giữa điểm quy chiếu và điểm xét. Nếu lưu lượng thực thay đổi dọc thiết bị, tỷ số V/Q tại cửa vào là một chỉ số danh nghĩa; thời gian đi theo một đường dòng còn phụ thuộc vận tốc cục bộ.
4. Space velocity cần ghi cả tử số và mẫu số
IUPAC mô tả space velocity trong xúc tác qua tốc độ cấp chất phản ứng tính trên khối lượng, thể tích hoặc diện tích xúc tác.[2] Do đó một con số mang đơn vị nghịch đảo giờ chưa đủ để biết thể tích nào đã được dùng; cần đọc định nghĩa cụ thể trong tài liệu thiết kế.
Ví dụ tự định nghĩa LHSV = Q_liquid,ref/V_bed, trong đó Q_liquid,ref = 2 m³/h là thể tích cấp ở điều kiện quy chiếu đã ghi và V_bed = 10 m³ là thể tích lớp hạt. LHSV bằng 0,2 h⁻¹; nghịch đảo bằng 5 h. Nếu lượng lỏng thật giữ trong lớp chỉ 1,2 m³ và dòng lỏng thực cũng bằng 2 m³/h, thời gian theo lượng giữ lỏng là 0,6 h, tức 36 phút.
Hai số 5 h và 0,6 h không mâu thuẫn: chúng có hai mẫu số thể tích khác nhau. Cũng không được lấy nghịch đảo WHSV có cơ sở khối lượng rồi gọi là thời gian dịch chuyển của chất lỏng nếu chưa có phép nối giữa lượng giữ, lưu lượng và khối lượng xúc tác.
5. Phân bố thời gian lưu chứa nhiều thông tin hơn giá trị trung bình
Trong phép thử xung tracer lý tưởng, chất đánh dấu đi cùng pha cần khảo sát, không phản ứng, không hấp phụ đáng kể và được đo đủ tại cửa ra. Với lưu lượng không đổi và tín hiệu nền đã trừ, đặt:
E(t) = C(t)/integral[C(t)dt]
t_mean = integral[tE(t)dt]
variance = integral[(t − t_mean)²E(t)dt]
Các tích phân bao toàn bộ đáp ứng từ lúc đưa xung vào đến khi thu đủ tín hiệu. Tài liệu thí nghiệm của Columbia trình bày việc lấy mômen tín hiệu tracer để suy thời gian trung bình và thể tích lượng giữ.[3] Hệ số khuếch đại tuyến tính chung của C(t) triệt tiêu trong tỷ số mômen, nhưng điều đó không chữa được tín hiệu bão hòa hoặc nền bị lệch.
Ví dụ rời rạc tự xây dựng có bốn nhóm tracer ra tại 1; 2; 3; 5 giờ, với tỷ phần khối lượng lần lượt 0,10; 0,20; 0,50; 0,20. Các tỷ phần cộng thành 1. Trung bình bằng 0,10 × 1 + 0,20 × 2 + 0,50 × 3 + 0,20 × 5 = 3 h. Phương sai bằng 1,4 h².
Có 30% tracer đã ra không muộn hơn giờ thứ hai, dù thời gian trung bình là 3 h. Một mô hình mọi phần tử đều ra đúng giờ thứ ba cũng có trung bình 3 h nhưng phương sai bằng 0. Vì thế không thể dùng một giá trị trung bình để khẳng định mọi phần tử đã ở trong thiết bị ít nhất từng ấy thời gian.
6. Bỏ đuôi tín hiệu có thể tạo ra thể tích chết giả
Giả sử phép đo của ví dụ trên dừng trước nhóm 5 h. Chỉ 80% tracer được thu. Nếu người xử lý tự chuẩn hóa phần đã thu thành 100%, thời gian trung bình biểu kiến trở thành:
t_apparent = (0,10 × 1 + 0,20 × 2 + 0,50 × 3)/0,80
t_apparent = 2,5 h
Với Q = 2 m³/h, thể tích suy ra chỉ 5 m³ thay vì 6 m³. Nếu diễn giải ngay thành phần chết, sẽ báo thiếu 16,666667% thể tích dù nguyên nhân trong bộ số này chỉ là bỏ 20% khối lượng tracer ở cuối.
IAEA cũng lưu ý giá trị trung bình không phù hợp có thể liên quan đến sai số lưu lượng, thể tích, tương tác tracer hoặc vùng trao đổi chậm.[1] Một tỷ số nhỏ hơn dự kiến chưa tự xác định cơ chế. Cần kiểm phần thu hồi và khoảng thời gian đo trước khi gán ý nghĩa vật lý cho tham số mô hình.
7. Đặt mốc thời gian và biên thử nghiệm đúng
Nếu tracer được đưa vào ở đầu một đoạn ống dẫn trước thiết bị, nhưng cảm biến đặt sau một đoạn ống khác, đáp ứng bao cả ống và thiết bị. Lượng giữ trong đường nối không tự biến mất khi xử lý số liệu. Với một độ trễ vận chuyển cố định 0,2 h được xác định riêng, trung bình đo 3,2 h có thể tách thành 3 h của phần cần xét và 0,2 h đường nối trong mô hình nối tiếp ấy.
Khi xung vào có độ rộng đáng kể, đường ra là kết quả chồng của đáp ứng thiết bị với dạng xung vào. Chỉ trừ một thời gian đỉnh không luôn đủ để suy toàn bộ phân bố. Tương tự, nếu lưu lượng thay đổi theo thời gian, nồng độ tại cửa ra phải được gắn với lượng tracer vận chuyển, thay vì chuẩn hóa nồng độ như khi Q không đổi.
8. Hồ sơ cần giữ lại những gì
Một báo cáo nên ghi thể tích hình học, thể tích từng pha, phần chiếm chỗ, cơ sở lưu lượng, nhiệt độ, áp suất tuyệt đối và phương pháp xác định hệ số nén. Với thử tracer, lưu thêm vị trí vào ra, dạng xung, tín hiệu nền, thời gian thu, phần thu hồi và quy tắc tích phân.
Thời gian lưu là dữ liệu đầu vào cho đánh giá tiếp xúc, phản ứng hoặc tách pha. Nó không tự thay thế động học hay mô hình truyền khối. Chọn đúng cơ sở giúp phép tính có ý nghĩa; phân bố thời gian lưu giúp nhìn thấy những khác biệt mà một phép chia thể tích cho lưu lượng chưa thể hiện.
Nguồn tham khảo
[1] Radiotracer Residence Time Distribution Method for Industrial and Environmental Applications, Training Course Series31 — International Atomic Energy Agency. Truy cập ngày 09/11/2026.
[2] Space velocity, S05779 — IUPAC Compendium of Chemical Terminology,5th edition,2025. Truy cập ngày 09/11/2026.
[3] ChE E3820y: Packed Bed Dynamics — Residence Time Distribution Theory — Columbia University, Chemical Engineering. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Thời gian lưu thường được ghi bằng thể tích chia lưu lượng thể tích. Nhưng thể tích có thể là toàn vỏ thiết bị, phần chứa chất lỏng, khoảng rỗng giữa hạt xúc tác hoặc thể tích thực tham gia dòng chảy. Lưu lượng lại có thể được báo ở nhiệt độ và áp suất thực, hoặc ở điều kiện quy chiếu. Dùng hai số không cùng cơ sở sẽ tạo ra một “thời gian” không đại diện cho chuyển động của pha đang xét.
Tài liệu IAEA phân biệt thời gian trung bình suy từ thể tích và lưu lượng với giá trị thực nghiệm thu từ phân bố thời gian lưu, đồng thời xem xét thể tích chết và vùng trao đổi chậm.[1] Trong bài này, các ví dụ đều là giả định để kiểm tra cơ sở tính, không phải kết luận về một thiết bị đang vận hành.
2. Từ thể tích vỏ tới lượng pha thực có mặt
Xét một bồn có thể tích hình học 10 m³. Mức chất lỏng chiếm vùng hình học 7 m³; vật rắn và chi tiết nằm ngập chiếm 1 m³. Nếu bỏ qua lỗ rỗng chứa lỏng bên trong vật rắn, lượng pha lỏng có mặt là 6 m³. Với dòng lỏng không nén được 2 m³/h, thời gian danh nghĩa theo pha lỏng bằng:
tau_l = V_l/Q_l = 6/2 = 3 h
Dùng toàn vỏ sẽ cho 5 h, cao hơn 66,666667%. Dùng vùng dưới mức mà quên thể tích chiếm chỗ sẽ cho 3,5 h. Cả hai lỗi xuất phát từ cách định nghĩa thể tích, không phải phép làm tròn.
Giả định thêm 0,6 m³ lỏng bị cô lập hoàn toàn, không trao đổi với phần còn lại trong thí nghiệm. Thể tích tham gia là 5,4 m³ và tỷ số tương ứng 2,7 h. Nếu vùng đó chỉ trao đổi rất chậm, cách thay đơn giản 6 thành 5,4 m³ có thể không mô tả đầy đủ phần đuôi của đáp ứng. Cần phân biệt “không vào được” với “vào được nhưng trở ra chậm”.
3. Khí phải dùng lưu lượng tại đúng điều kiện
Xét một thiết bị khí độc lập có khoảng khí 6 m³. Lưu lượng quy chiếu là 1.000 m³/h tại 273,15 K, 101.325 Pa và Z_n = 1. Tại thiết bị, giả định 450 K, áp suất tuyệt đối 2.000.000 Pa, hệ số nén Z_a = 0,9. Cùng một dòng mol cho:
Q_a = Q_n (Z_a/Z_n)(T_a/T_n)(p_n/p_a)
Q_a = 75,117381 m³/h
tau_a = 6/Q_a = 4,792499 min
Nếu chia trực tiếp 6 m³ cho 1.000 m³/h, kết quả chỉ 0,36 phút. Chênh lệch không phải khí thực đang lưu “lâu hơn thiết kế”; nó do đã ghép thể tích thực của thiết bị với lưu lượng ở điều kiện khác.
Mô hình quy đổi trên giả định không có thay đổi mol do phản ứng hoặc ngưng tụ giữa điểm quy chiếu và điểm xét. Nếu lưu lượng thực thay đổi dọc thiết bị, tỷ số V/Q tại cửa vào là một chỉ số danh nghĩa; thời gian đi theo một đường dòng còn phụ thuộc vận tốc cục bộ.
4. Space velocity cần ghi cả tử số và mẫu số
IUPAC mô tả space velocity trong xúc tác qua tốc độ cấp chất phản ứng tính trên khối lượng, thể tích hoặc diện tích xúc tác.[2] Do đó một con số mang đơn vị nghịch đảo giờ chưa đủ để biết thể tích nào đã được dùng; cần đọc định nghĩa cụ thể trong tài liệu thiết kế.
Ví dụ tự định nghĩa LHSV = Q_liquid,ref/V_bed, trong đó Q_liquid,ref = 2 m³/h là thể tích cấp ở điều kiện quy chiếu đã ghi và V_bed = 10 m³ là thể tích lớp hạt. LHSV bằng 0,2 h⁻¹; nghịch đảo bằng 5 h. Nếu lượng lỏng thật giữ trong lớp chỉ 1,2 m³ và dòng lỏng thực cũng bằng 2 m³/h, thời gian theo lượng giữ lỏng là 0,6 h, tức 36 phút.
Hai số 5 h và 0,6 h không mâu thuẫn: chúng có hai mẫu số thể tích khác nhau. Cũng không được lấy nghịch đảo WHSV có cơ sở khối lượng rồi gọi là thời gian dịch chuyển của chất lỏng nếu chưa có phép nối giữa lượng giữ, lưu lượng và khối lượng xúc tác.
5. Phân bố thời gian lưu chứa nhiều thông tin hơn giá trị trung bình
Trong phép thử xung tracer lý tưởng, chất đánh dấu đi cùng pha cần khảo sát, không phản ứng, không hấp phụ đáng kể và được đo đủ tại cửa ra. Với lưu lượng không đổi và tín hiệu nền đã trừ, đặt:
E(t) = C(t)/integral[C(t)dt]
t_mean = integral[tE(t)dt]
variance = integral[(t − t_mean)²E(t)dt]
Các tích phân bao toàn bộ đáp ứng từ lúc đưa xung vào đến khi thu đủ tín hiệu. Tài liệu thí nghiệm của Columbia trình bày việc lấy mômen tín hiệu tracer để suy thời gian trung bình và thể tích lượng giữ.[3] Hệ số khuếch đại tuyến tính chung của C(t) triệt tiêu trong tỷ số mômen, nhưng điều đó không chữa được tín hiệu bão hòa hoặc nền bị lệch.
Ví dụ rời rạc tự xây dựng có bốn nhóm tracer ra tại 1; 2; 3; 5 giờ, với tỷ phần khối lượng lần lượt 0,10; 0,20; 0,50; 0,20. Các tỷ phần cộng thành 1. Trung bình bằng 0,10 × 1 + 0,20 × 2 + 0,50 × 3 + 0,20 × 5 = 3 h. Phương sai bằng 1,4 h².
Có 30% tracer đã ra không muộn hơn giờ thứ hai, dù thời gian trung bình là 3 h. Một mô hình mọi phần tử đều ra đúng giờ thứ ba cũng có trung bình 3 h nhưng phương sai bằng 0. Vì thế không thể dùng một giá trị trung bình để khẳng định mọi phần tử đã ở trong thiết bị ít nhất từng ấy thời gian.
6. Bỏ đuôi tín hiệu có thể tạo ra thể tích chết giả
Giả sử phép đo của ví dụ trên dừng trước nhóm 5 h. Chỉ 80% tracer được thu. Nếu người xử lý tự chuẩn hóa phần đã thu thành 100%, thời gian trung bình biểu kiến trở thành:
t_apparent = (0,10 × 1 + 0,20 × 2 + 0,50 × 3)/0,80
t_apparent = 2,5 h
Với Q = 2 m³/h, thể tích suy ra chỉ 5 m³ thay vì 6 m³. Nếu diễn giải ngay thành phần chết, sẽ báo thiếu 16,666667% thể tích dù nguyên nhân trong bộ số này chỉ là bỏ 20% khối lượng tracer ở cuối.
IAEA cũng lưu ý giá trị trung bình không phù hợp có thể liên quan đến sai số lưu lượng, thể tích, tương tác tracer hoặc vùng trao đổi chậm.[1] Một tỷ số nhỏ hơn dự kiến chưa tự xác định cơ chế. Cần kiểm phần thu hồi và khoảng thời gian đo trước khi gán ý nghĩa vật lý cho tham số mô hình.
7. Đặt mốc thời gian và biên thử nghiệm đúng
Nếu tracer được đưa vào ở đầu một đoạn ống dẫn trước thiết bị, nhưng cảm biến đặt sau một đoạn ống khác, đáp ứng bao cả ống và thiết bị. Lượng giữ trong đường nối không tự biến mất khi xử lý số liệu. Với một độ trễ vận chuyển cố định 0,2 h được xác định riêng, trung bình đo 3,2 h có thể tách thành 3 h của phần cần xét và 0,2 h đường nối trong mô hình nối tiếp ấy.
Khi xung vào có độ rộng đáng kể, đường ra là kết quả chồng của đáp ứng thiết bị với dạng xung vào. Chỉ trừ một thời gian đỉnh không luôn đủ để suy toàn bộ phân bố. Tương tự, nếu lưu lượng thay đổi theo thời gian, nồng độ tại cửa ra phải được gắn với lượng tracer vận chuyển, thay vì chuẩn hóa nồng độ như khi Q không đổi.
8. Hồ sơ cần giữ lại những gì
Một báo cáo nên ghi thể tích hình học, thể tích từng pha, phần chiếm chỗ, cơ sở lưu lượng, nhiệt độ, áp suất tuyệt đối và phương pháp xác định hệ số nén. Với thử tracer, lưu thêm vị trí vào ra, dạng xung, tín hiệu nền, thời gian thu, phần thu hồi và quy tắc tích phân.
Thời gian lưu là dữ liệu đầu vào cho đánh giá tiếp xúc, phản ứng hoặc tách pha. Nó không tự thay thế động học hay mô hình truyền khối. Chọn đúng cơ sở giúp phép tính có ý nghĩa; phân bố thời gian lưu giúp nhìn thấy những khác biệt mà một phép chia thể tích cho lưu lượng chưa thể hiện.
Nguồn tham khảo
[1] Radiotracer Residence Time Distribution Method for Industrial and Environmental Applications, Training Course Series31 — International Atomic Energy Agency. Truy cập ngày 09/11/2026.
[2] Space velocity, S05779 — IUPAC Compendium of Chemical Terminology,5th edition,2025. Truy cập ngày 09/11/2026.
[3] ChE E3820y: Packed Bed Dynamics — Residence Time Distribution Theory — Columbia University, Chemical Engineering. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
