Tổng quan nghiên cứu

Trong công tác quản lý và xử lý các khu vực bị ô nhiễm môi trường đất và nước ngầm, việc xác định chính xác độ tan của các hóa chất hữu cơ kỵ nước giữ vai trò then chốt trong việc dự báo khả năng lan truyền và đánh giá rủi ro sinh thái. Trên thực tế hiện trường, các chất ô nhiễm hiếm khi tồn tại ở dạng đơn chất tinh khiết mà thường hòa trộn phức tạp thành các pha hữu cơ không phân cực, điển hình như hỗn hợp dầu khoáng, hợp chất thơm bay hơi (BTEX), hydrocarbon thơm đa vòng (PAHs) và phenol. Dữ liệu từ các trạm quan trắc cho thấy, mặc dù độ tan lý thuyết của benzene đơn chất trong nước có thể đạt xấp xỉ 1.780 mg/l, nhưng nồng độ thực tế ghi nhận tại các tầng chứa nước bị ô nhiễm dầu thường chỉ dao động quanh mức 50 mg/l, tức giảm hơn 97% so với giá trị tra cứu trong các sổ tay hóa học.

Nghiên cứu này được thực hiện nhằm giải quyết khoảng trống giữa lý thuyết đơn chất và hành vi thực tế của các chất ô nhiễm trong hỗn hợp. Mục tiêu cụ thể của đề tài là khảo sát thực nghiệm sự thay đổi độ tan bão hòa trong nước của 6 hợp chất hữu cơ đại diện gồm Ethylbenzene, Toluene, Phenol, Anthracene, Naphthalene và Tetradecane ở cả dạng đơn chất lẫn hệ hỗn hợp đa cấu tử tại các mức nhiệt độ 5°C, 10°C và 20°C. Đề tài được triển khai trong khuôn khổ hợp tác học thuật giữa Trường Đại học Khoa học Tự nhiên (ĐHQG Hà Nội) và Viện Quản lý Chất thải & Xử lý Điểm ô nhiễm thuộc Đại học Kỹ thuật Dresden (CHLB Đức). Kết quả nghiên cứu cung cấp cơ sở dữ liệu thực nghiệm chuẩn xác, hỗ trợ đắc lực cho việc hiệu chỉnh các mô hình mô phỏng lan truyền chất ô nhiễm như SiWaPro DSS, nâng cao độ tin cậy của các dự báo thủy văn môi trường tại Việt Nam và quốc tế.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu vận dụng nền tảng nhiệt động học cân bằng pha lỏng - lỏng và rắn - lỏng để giải thích hành vi hòa tan của các hợp chất hữu cơ kỵ nước. Các khung lý thuyết và mô hình chính bao gồm:

  • Định luật Raoult và Cân bằng pha nhiệt động: Mô tả mối quan hệ giữa phân số mol của từng cấu tử trong pha hữu cơ và áp suất hơi riêng phần cũng như độ tan cân bằng của chúng trong pha nước. Đối với các dung dịch lý tưởng, độ tan của từng chất trong hỗn hợp tỷ lệ thuận trực tiếp với nồng độ mol của chính chất đó trong pha dầu.
  • Mô hình đóng góp nhóm UNIFAC và UNIFAC cải tiến: Phương pháp ước tính hệ số hoạt độ nhiệt động dựa trên cấu trúc các nhóm chức hóa học của phân tử. Mô hình cho phép tính toán sự sai lệch so với trạng thái lý tưởng trong pha hữu cơ và pha nước bằng cách phân tích tương tác năng lượng giữa các nhóm chức ở các mức nhiệt độ từ 5°C đến 20°C.
  • Hiệu ứng muối kết tủa và tác động đồng dung môi: Ứng dụng phương trình Setschenow để đánh giá ảnh hưởng của các ion khoáng và mô hình log-tuyến tính của Morris để phân tích hiện tượng tăng độ tan khi có sự hiện diện của các dung môi hữu cơ phân cực hòa tan một phần.

Các khái niệm trọng tâm xuyên suốt nghiên cứu gồm: độ tan bão hòa pha nước, hệ số hoạt độ nhiệt động, hệ số phân bố octanol - nước ($K_{ow}$), và cân bằng pha lỏng - lỏng (LLE).

Phương pháp nghiên cứu

Nghiên cứu được thiết kế thực nghiệm nghiêm ngặt với cỡ mẫu gồm 6 chất hữu cơ tinh khiết cao (Ethylbenzene ≥ 99%, Toluene ≥ 99.5%, Phenol ≥ 99.5%, Anthracene ≥ 99%, Naphthalene ≥ 99% và Tetradecane chuẩn sắc ký). Phương pháp chọn mẫu áp dụng kỹ thuật lấy mẫu song song với 2 đến 3 bình phản ứng độc lập cho mỗi điều kiện khảo sát; mỗi điểm dữ liệu được phân tích lặp lại 3 lần để tính giá trị trung bình, độ lệch chuẩn và khoảng tin cậy thống kê Student ở mức ý nghĩa 95%.

Quy trình phân tích hóa lý chuyên sâu được lựa chọn dựa trên đặc tính bay hơi và độ tan của từng nhóm chất:

  • Nhóm VOCs (Ethylbenzene, Toluene): Phân tích bằng hệ thống Sắc ký khí tích hợp bộ lấy mẫu không gian hơi tự động (GC-FID/Headspace) trên cột mao quản Agilent DB-624 (30m x 0.32mm), điều kiện ủ nhiệt 70°C trong 180 phút, đạt giới hạn phát hiện 5 µg/l cho Ethylbenzene và 75 µg/l cho Toluene.
  • Nhóm PAHs (Anthracene, Naphthalene): Phân tích bằng Sắc ký lỏng hiệu năng cao (HPLC WATERS 2695) sử dụng đầu dò PDA tại bước sóng 220 nm và đầu dò huỳnh quang với bước sóng kích thích 235 - 280 nm, phát xạ 340 - 430 nm trên cột chuyên dụng MZ-PAH C18 (5 µm, 25 cm).
  • Nhóm Phenol và Alkane: Phenol được định lượng bằng quang phổ hấp thụ phân tử UV-Vis ở bước sóng cực đại 271 nm; Tetradecane ở nồng độ vết được chiết làm giàu 2 lần bằng n-Heptane qua cột Na2SO4/Florisil và định lượng trên hệ thống GC-FID Shimadzu 2010 với chất chuẩn nội Dodecane.

Lý do lựa chọn tổ hợp phương pháp này nhằm đảm bảo tính chọn lọc cao, loại trừ hoàn toàn hiện tượng nhiễu nền khi phân tích các hỗn hợp đa cấu tử có nồng độ chênh lệch từ mức vài µg/l đến hàng nghìn mg/l.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình phân tích thực nghiệm đã mang lại 4 phát hiện quan trọng về động học và cân bằng hòa tan:

  1. Thời gian đạt trạng thái cân bằng hòa tan: Các hợp chất bay hơi (VOCs) như Ethylbenzene và Toluene đạt trạng thái bão hòa ổn định trong pha nước sau 5 ngày, trong khi các hợp chất đa vòng (PAHs) như Anthracene và Naphthalene cần 7 ngày. Dữ liệu theo dõi nồng độ Anthracene sau 7 ngày đạt 16.32 µg/l, sau 14 ngày đạt 17.01 µg/l và sau 28 ngày duy trì ở mức 16.96 µg/l, chứng minh hệ đạt trạng thái cân bằng bền vững với độ lệch dưới 4%.
  2. Độ tan đơn chất ở các mức nhiệt độ: Tại 20°C, độ tan thực nghiệm của Ethylbenzene đạt 180.52 mg/l (tương ứng phân số mol $3.54 \times 10^{-5}$), hoàn toàn tương thích với giá trị chuẩn 181 mg/l do cơ sở dữ liệu quốc tế IUPAC-NIST khuyến nghị. Nồng độ Toluene bão hòa ở 20°C đạt 553.42 mg/l (phân số mol $10.81 \times 10^{-5}$), nằm trọn trong dải y văn quốc tế từ 518 đến 608 mg/l.
  3. Ảnh hưởng của chế độ khuấy trộn: Thử nghiệm đối chứng cho thấy mẫu Ethylbenzene có khuấy liên tục cho nồng độ biểu kiến trong nước cao hơn từ 8% đến 14% so với mẫu không khuấy (đạt 180.19 - 190.41 mg/l so với ~163.78 mg/l). Nguyên nhân được xác định là do lực khuấy cơ học tạo ra các vi giọt dung môi lơ lửng có kích thước từ 2 đến 300 µm, nhưng trạng thái tĩnh vẫn phản ánh chính xác độ tan bão hòa nhiệt động học.
  4. Sự sụt giảm độ tan trong hệ hỗn hợp: Khi phối trộn trong hỗn hợp đa cấu tử (như hệ Ethylbenzene - Toluene - Anthracene - Naphthalene), độ tan trong nước của từng đơn chất suy giảm mạnh theo tỷ lệ phân số mol của chúng trong pha hữu cơ. Độ tan của Anthracene giảm từ mức đơn chất 16.32 µg/l xuống dưới 5 µg/l khi có mặt các dung môi hữu cơ đồng hành.

Thảo luận kết quả

Các kết quả thực nghiệm chỉ ra rằng sự thay đổi độ tan trong hỗn hợp chịu sự chi phối chặt chẽ bởi tương tác liên phân tử bên trong pha hữu cơ hơn là tương tác trong pha nước. Khi các chất ô nhiễm kỵ nước cùng tồn tại, pha hữu cơ đóng vai trò như một dung môi lý tưởng hoặc bán lý tưởng, làm giảm hoạt độ nhiệt động của từng cấu tử truyền khối sang pha nước.

Dữ liệu nghiên cứu có thể được hệ thống hóa qua biểu đồ tương quan nhiệt độ - độ tan, minh chứng rằng độ tan của hầu hết các hydrocacbon thơm tăng đều theo nhiệt độ từ 5°C lên 20°C. Sự sai khác giữa số liệu thực tế ngoài hiện trường bãi chôn lấp (nơi nồng độ benzene chỉ đạt 50 mg/l) so với giá trị lý thuyết đơn chất đã được sáng tỏ: nồng độ chất ô nhiễm giải phóng vào nước ngầm bị khống chế bởi thành phần pha dầu ô nhiễm theo đúng định luật Raoult và hệ số hoạt độ UNIFAC. Việc so sánh giữa giá trị tính toán lý thuyết và giá trị đo đạc sắc ký cho thấy sai số chỉ dao động trong khoảng 5% đến 12%, khẳng định tính khả thi của việc ứng dụng mô hình nhiệt động vào mô phỏng môi trường.

Đề xuất và khuyến nghị

Dựa trên các kết luận thực nghiệm, nghiên cứu đưa ra 4 khuyến nghị hành động cụ thể cho công tác quản lý và xử lý ô nhiễm:

  1. Tích hợp mô hình nhiệt động UNIFAC vào phần mềm mô phỏng: Đề xuất các viện nghiên cứu và cơ quan quản lý tài nguyên tích hợp thuật toán tính toán hệ số hoạt độ hỗn hợp vào các công cụ mô phỏng lan truyền nước ngầm (như SiWaPro DSS) nhằm kéo giảm sai số dự báo phát tán chất ô nhiễm xuống dưới 10% trong vòng 12 tháng tới.
  2. Chuẩn hóa quy trình quan trắc và lấy mẫu hiện trường: Cục Kiểm soát ô nhiễm môi trường và các trung tâm quan trắc cần ban hành hướng dẫn kỹ thuật lấy mẫu nước ngầm tại các khu vực rò rỉ xăng dầu, quy định rõ việc sử dụng thiết bị vi tiêm qua màng ngăn và duy trì thời gian lắng tĩnh tối thiểu 5 ngày để loại bỏ hoàn toàn sai số do vi giọt dầu phân tán.
  3. Thiết lập cơ sở dữ liệu quốc gia về độ tan hỗn hợp chất ô nhiễm: Khuyến nghị Bộ Tài nguyên và Môi trường chủ trì xây dựng bộ thông số cân bằng pha cho các tổ hợp ô nhiễm đặc thù (BTEX, PAHs, Phenol, dầu khoáng) trong dải nhiệt độ nhiệt đới từ 10°C đến 30°C trong giai đoạn 2026 - 2028.
  4. Ứng dụng công nghệ xử lý đất ô nhiễm dựa trên cân bằng pha: Các đơn vị xử lý môi trường cần áp dụng giải pháp rửa giải tăng cường bằng chất hoạt động bề mặt (Surfactant Flushing) có tính toán đến phân số mol pha dầu, nhằm tối ưu hóa hiệu suất bóc tách hydrocacbon nặng lên trên 85% trong thời gian xử lý từ 6 đến 18 tháng.

Đối tượng nên tham khảo luận văn

Nội dung và phương pháp luận của luận văn mang lại giá trị thực tiễn cao cho 4 nhóm đối tượng:

  1. Nhà nghiên cứu và Giảng viên ngành Kỹ thuật Môi trường / Hóa học: Nguồn tài liệu tham khảo chuyên sâu về nhiệt động học môi trường, kỹ thuật sắc ký chuyên biệt cho VOCs, PAHs và phương pháp tính toán hệ số hoạt độ UNIFAC.
  2. Kỹ sư tư vấn Đánh giá Tác động Môi trường (ĐTM) và Mô hình hóa: Cung cấp thông số đầu vào chuẩn xác để thiết lập mô hình lan truyền chất ô nhiễm trong tầng chứa nước ngầm, phục vụ công tác thẩm định dự án công nghiệp.
  3. Cơ quan quản lý nhà nước về Môi trường và Địa chất: Cơ sở khoa học để xây dựng quy chuẩn kỹ thuật quốc gia về ngưỡng đánh giá rủi ro tại các bãi chôn lấp chất thải và khu công nghiệp lọc hóa dầu.
  4. Kỹ thuật viên phòng thí nghiệm phân tích môi trường: Quy trình chi tiết về chuẩn bị mẫu, xử lý màng ngăn, tối ưu hóa thông số GC-FID Headspace và HPLC huỳnh quang cho các hợp chất hữu cơ khó phân tích.

Câu hỏi thường gặp

Tại sao nồng độ benzene ghi nhận tại hiện trường bãi ô nhiễm lại thấp hơn nhiều so với độ tan đơn chất lý thuyết?
Tại các điểm ô nhiễm thực tế, benzene không tồn tại đơn độc mà hòa tan trong pha dầu khoáng cùng hàng trăm hợp chất hữu cơ khác. Phân số mol của benzene trong pha dầu bị pha loãng, làm giảm hoạt độ nhiệt động truyền khối sang pha nước xuống mức khoảng 50 mg/l so với độ tan bão hòa đơn chất lý thuyết là khoảng 1.780 mg/l.

Thời gian cần thiết để các chất ô nhiễm hữu cơ đạt cân bằng hòa tan trong pha nước là bao lâu?
Nghiên cứu chứng minh các hợp chất hữu cơ dễ bay hơi (VOCs) như Ethylbenzene đạt trạng thái cân bằng hòa tan sau 5 ngày, trong khi các hydrocarbon thơm đa vòng (PAHs) như Anthracene cần 7 ngày để nồng độ trong nước đạt mức bão hòa ổn định 16.32 µg/l và không thay đổi suốt 28 ngày sau đó.

Kỹ thuật lấy mẫu nước có ảnh hưởng như thế nào đến độ chính xác của phép đo độ tan?
Việc khuấy trộn mạnh hoặc lấy mẫu không cẩn thận có thể cuốn theo các vi giọt dung môi lơ lửng (kích thước 2 - 300 µm), làm nồng độ đo được tăng ảo từ 8% đến 14%. Việc sử dụng vi tiêm Hamilton qua màng cao su silicone trong điều kiện mẫu tĩnh giúp loại bỏ hoàn toàn hạt cặn và giọt dầu.

Phương pháp sắc ký nào phù hợp nhất để xác định các hợp chất PAHs ở nồng độ vết trong nước?
Phương pháp HPLC sử dụng kết hợp đầu dò quang tuyến mảng (PDA) tại bước sóng 220 nm và đầu dò huỳnh quang với bước sóng phát xạ 340 - 430 nm là giải pháp tối ưu, cho phép định lượng chính xác Anthracene và Naphthalene ở mức giới hạn phát hiện dưới 1 µg/l.

Mô hình UNIFAC cải tiến hỗ trợ gì trong việc đánh giá sự lan truyền của chất ô nhiễm?
Mô hình UNIFAC phân tích các nhóm chức phân tử để tính toán chính xác hệ số hoạt độ của từng cấu tử trong hỗn hợp phức tạp. Nhờ đó, các chuyên gia có thể dự báo độ tan thực tế mà không cần tiến hành thực nghiệm tốn kém cho mọi tổ hợp hóa chất phát sinh tại hiện trường.

Kết luận

  • Luận văn đã xác định thành công độ tan bão hòa trong nước của 6 hợp chất hữu cơ nguy hại (Ethylbenzene, Toluene, Phenol, Anthracene, Naphthalene, Tetradecane) ở dải nhiệt độ 5°C, 10°C và 20°C.
  • Dữ liệu thực nghiệm của Ethylbenzene (180.52 mg/l) và Toluene (553.42 mg/l) ở 20°C có độ tương thích cao với cơ sở dữ liệu quốc tế IUPAC-NIST, khẳng định độ tin cậy của quy trình phân tích.
  • Nghiên cứu làm sáng tỏ cơ chế sụt giảm độ tan của các cấu tử trong hệ hỗn hợp đa thành phần, giải thích thỏa đáng sự chênh lệch nồng độ chất ô nhiễm giữa thực tế hiện trường và lý thuyết đơn chất.
  • Đóng góp bộ thông số thực nghiệm quan trọng hỗ trợ chuẩn hóa công cụ tính toán thủy văn và mô phỏng lan truyền chất ô nhiễm trong đất và nước ngầm tại Việt Nam.
  • Trong giai đoạn tiếp theo (2026 - 2028), các nghiên cứu cần tiếp tục mở rộng khảo sát trên các nền mẫu nước ngầm thực tế có độ mặn và hàm lượng khoáng cao; các đơn vị quản lý môi trường được khuyến khích áp dụng ngay các phát hiện này vào quy trình đánh giá rủi ro hiện trường.