Tổng quan nghiên cứu

Lưu vực sông Cầu đóng vai trò huyết mạch trong mạng lưới thủy văn phía Bắc Việt Nam với chiều dài 288,5 km, diện tích lưu vực xấp xỉ 6.030 km² và lưu lượng nước tự nhiên hàng năm đạt khoảng 4,5 tỷ m³. Nguồn nước này đang chịu áp lực ô nhiễm nặng nề từ 16.000 m³/ngày đêm nước thải ngành luyện kim và hơn 35.000 m³/ngày nước thải sản xuất giấy tại khu vực công nghiệp tỉnh Thái Nguyên. Ô nhiễm kim loại nặng trong môi trường nước không chỉ giới hạn ở pha hòa tan mà còn tích lũy với nồng độ cao trong trầm tích đáy, tạo ra nguy cơ giải phóng ngược vào chuỗi thức ăn thủy sinh khi điều kiện lý hóa thay đổi.

Đề tài tập trung giải quyết bài toán định lượng hàm lượng tổng và phân tích 5 dạng liên kết hóa học của bốn kim loại nặng điển hình gồm Đồng (Cu), Chì (Pb), Kẽm (Zn) và Cadimi (Cd). Mục tiêu then chốt là đánh giá biến thiên hàm lượng kim loại theo chiều sâu bồi lắng tại 3 cột trầm tích đại diện, xác định mức độ ô nhiễm thông qua chỉ số tích lũy địa chất (Igeo) và đánh giá rủi ro sinh thái bằng mã RAC (Risk Assessment Code).

Nghiên cứu được triển khai thực địa vào tháng 7 năm 2012 trên đoạn sông Cầu chảy qua thành phố Thái Nguyên, kéo dài từ cầu Gia Bảy đến đập Ba Đa. Kết quả cung cấp bộ thông số hóa học định lượng chuẩn xác, phục vụ trực tiếp cho công tác kiểm soát ô nhiễm của Ủy ban Bảo vệ môi trường lưu vực sông Cầu đối với địa bàn tập trung hơn 2,8 triệu cư dân sinh sống.

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 lý thuyết phân tích dạng hóa học của Tessier (1979) kết hợp với các hiệu chỉnh của Viện Hóa học (2010). Phân tích dạng phân tách kim loại thành 5 phân đoạn có độ bền liên kết và hoạt tính sinh thái khác nhau: dạng trao đổi (F1), dạng liên kết với cacbonat (F2), dạng liên kết với Fe-Mn oxit (F3), dạng liên kết với chất hữu cơ (F4) và dạng cặn dư silicat bền vững (F5). Độ linh động và độc tính của kim loại giảm dần từ phân đoạn F1 đến F5.

Mức độ rủi ro sinh thái được xác định thông qua hai mô hình quốc tế chuẩn hóa:

  • Chỉ số tích lũy địa chất ($I_{geo}$) do Muller thiết lập, so sánh nồng độ thực tế với giá trị nền vỏ Trái Đất theo thang chuẩn của Hamilton (Cu: 55 mg/kg, Pb: 15 mg/kg, Zn: 70 mg/kg, Cd: 0,1 mg/kg).
  • Chỉ số đánh giá rủi ro sinh học RAC (Risk Assessment Code) tính bằng tổng tỷ lệ phần trăm của hai pha linh động ($F1 + F2$), phân chia rủi ro theo 5 mức từ không đáng kể (dưới 1%) đến rất cao (trên 50%).

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

Đối tượng khảo sát gồm 3 cột trầm tích nguyên dạng thu thập tại các vị trí trọng điểm: cầu treo Huống Thượng (SC01, độ sâu 80 cm), đập Ba Đa (SC02, độ sâu 90 cm) và cầu treo Oánh (SC03, độ sâu 110 cm). Mẫu được thu bằng thiết bị ống phóng PVC chuyên dụng gia tải trọng lực, bảo quản lạnh ở nhiệt độ -4°C và phân đoạn đều đặn 10 cm mỗi lớp, tạo thành 28 mẫu phân tích độc lập.

Quy trình phân tích áp dụng phương pháp quang phổ hấp thụ nguyên tử (AAS) trên thiết bị Perkin Elmer AAS-3300:

  • Hàm lượng tổng được vô cơ hóa hoàn toàn bằng 30 ml hỗn hợp nước cường thủy ($HNO_3:HCl = 1:3$) trên bếp cách cát ở 80°C.
  • Quy trình chiết 5 phân đoạn sử dụng lần lượt các hệ dung môi chọn lọc: $CH_3COONH_4$ 1M (pH 7, lắc 1 giờ), $CH_3COONH_4$ 1M axit hóa bằng $CH_3COOH$ (pH 5, lắc 5 giờ), $NH_2OH\cdot HCl$ 0,04M trong $CH_3COOH$ 25% (lắc 5 giờ), $CH_3COONH_4$ 3,2M trong $HNO_3$ 20% (lắc 30 phút) và cặn dư phá bằng cường thủy.
  • Kỹ thuật ngọn lửa (F-AAS) với hỗn hợp khí $C_2H_2$/không khí (tỷ lệ 2/10 lít/phút) được dùng để đo Cu (324,8 nm), Pb (283,3 nm) và Zn (213,9 nm). Kỹ thuật lò graphit không ngọn lửa (GF-AAS) kết hợp đèn D2 bổ chính nền được chỉ định cho Cd tại 283,3 nm với nhiệt độ nguyên tử hóa 1650°C để đạt độ nhạy tới mức microgam/kg.
  • Ảnh hưởng của nền được triệt tiêu bằng cách xây dựng 5 bộ đường chuẩn riêng biệt cho 5 hệ dung dịch chiết, đạt hệ số tương quan $R^2 > 0,999$ và giới hạn phát hiện phương pháp đạt dưới 0,035 mg/L đối với F-AAS.

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

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

Kết quả phân tích 28 phân đoạn mẫu trầm tích cho thấy sự tích lũy kim loại nặng phân hóa rõ rệt giữa các nguyên tố:

  • Cadimi (Cd) là kim loại có mức độ ô nhiễm cao nhất, với hàm lượng tổng dao động từ khoảng 0,45 mg/kg đến trên 6,5 mg/kg, vượt từ 4,5 đến hơn 65 lần so với giá trị nền tự nhiên (0,1 mg/kg). Hàm lượng Kẽm (Zn) dao động trong khoảng 120 đến 460 mg/kg, vượt ngưỡng ô nhiễm nghiêm trọng theo tiêu chuẩn US EPA (200 mg/kg).
  • Đồng (Cu) và Chì (Pb) ghi nhận hàm lượng trung bình lần lượt là 25 - 98 mg/kg và 20 - 85 mg/kg, tập trung chủ yếu ở mức ô nhiễm nhẹ đến trung bình theo thang đo Canada (2002).
  • Về dạng tồn tại, Cd và Zn thể hiện hoạt tính sinh học rất cao khi tổng hàm lượng ở dạng trao đổi (F1) và cacbonat (F2) chiếm từ 32% đến 56% tổng lượng Cd, và từ 24% đến 41% tổng lượng Zn. Ngược lại, Cu và Pb chủ yếu cố định trong pha hữu cơ F4 (25% - 40%) và pha cặn dư F5 (45% - 62%).
  • Phân bố theo chiều sâu bồi lắng chỉ ra rằng các lớp mặt (0 - 30 cm) có nồng độ kim loại cao hơn từ 1,6 đến 2,7 lần so với các tầng trầm tích đáy sâu (70 - 110 cm), phản ánh sự gia tăng phát thải công nghiệp trong các giai đoạn gần đây.

Thảo luận kết quả

Hàm lượng Cd và Zn cao ở các phân đoạn linh động gắn liền với nguồn tiếp nhận nước thải trực tiếp từ suối Cam Giá (khu gang thép Lưu Xá) và suối Văn Dương (khu công nghiệp Sông Công), nơi các chỉ số BOD5 vượt hơn 2 lần và amoni vượt 16 lần so với QCVN. Ở điều kiện pH nước mặt trung tính dao động từ 7,1 đến 7,2, sự chuyển hóa qua lại giữa pha rắn và pha hòa tan diễn ra chậm, nhưng khi xảy ra hiện tượng axit hóa cục bộ hoặc xáo trộn cơ học, lượng lớn Cd ở dạng F1 và F2 sẽ hòa tan ngược vào nước, gây độc trực tiếp cho hệ thủy sinh.

So sánh với các công bố khoa học tại lưu vực sông Nhuệ - Đáy (Vũ Đức Lợi, 2010), mức độ tích lũy Cd tại sông Cầu đoạn qua Thái Nguyên tương đương vùng hạ lưu sông Đáy, nhưng tỷ lệ dạng trao đổi F1 cao hơn khoảng 8%. Biến thiên dữ liệu phân tích được hệ thống hóa qua 28 biểu đồ cột nồng độ theo độ sâu và 5 bảng ma trận tương quan Pearson, cho thấy mối liên kết thuận chặt chẽ ($r > 0,82$) giữa hàm lượng Zn và Cd, chứng minh hai nguyên tố này cùng bắt nguồn từ hoạt động luyện kim màu và mạ điện công nghiệp.

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

  1. Nâng cấp và kiểm soát trạm xử lý nước thải công nghiệp: Ban Quản lý các Khu công nghiệp và Sở Tài nguyên và Môi trường tỉnh Thái Nguyên cần yêu cầu 100% cơ sở luyện kim, cán thép tại khu vực Lưu Xá và cụm công nghiệp Cao Ngạn lắp đặt hệ thống khử kim loại nặng đạt chuẩn cột A QCVN 40:2011/BTNMT, hướng tới mục tiêu giảm 85% tải lượng phát thải Cu, Pb, Zn, Cd ra sông Cầu trước năm 2028.
  2. Triển khai nạo vét và xử lý trầm tích tầng mặt có kiểm soát: Ủy ban Nhân dân tỉnh Thái Nguyên chỉ đạo thực hiện nạo vét chọn lọc lớp trầm tích ô nhiễm nghiêm trọng ở độ sâu 0 - 30 cm tại 3 điểm nóng thủy vực gồm hạ lưu cầu Gia Bảy, khu vực đập Ba Đa và cầu Huống Thượng trong giai đoạn 2027 - 2029, xử lý đóng rắn an toàn bùn thải nhằm giảm 60% nguy cơ giải phóng Cd dạng linh động.
  3. Thiết lập mạng lưới quan trắc phân tích dạng định kỳ: Sở Tài nguyên và Môi trường tỉnh Thái Nguyên phối hợp với Viện Hóa học triển khai chương trình giám sát chất lượng trầm tích định kỳ 6 tháng/lần, bổ sung chỉ số RAC và $I_{geo}$ vào hệ thống chỉ tiêu quan trắc lưu vực trước năm 2027.
  4. Xây dựng quy chuẩn kỹ thuật địa phương về ngưỡng rủi ro trầm tích: Bộ Tài nguyên và Môi trường cần sớm ban hành bộ hướng dẫn đánh giá chất lượng trầm tích nước ngọt nội địa dựa trên tỷ lệ phần trăm dạng chiết linh động ($F1 + F2 < 15%$), thay vì chỉ kiểm soát hàm lượng tổng số đơn thuần, áp dụng thí điểm cho toàn lưu vực sông Cầu vào năm 2029.

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

  • Học viên cao học và nghiên cứu sinh chuyên ngành Hóa phân tích, Khoa học Môi trường: Tiếp cận phương pháp luận chiết liên tục 5 bước cải tiến Tessier và kỹ thuật kiểm soát hiệu ứng nền phức tạp trên máy quang phổ hấp thụ nguyên tử F-AAS, GF-AAS.
  • Cán bộ quản lý môi trường tại các Sở Tài nguyên và Môi trường lưu vực sông Cầu: Sử dụng số liệu thực chứng và bộ chỉ số $I_{geo}$, RAC để xây dựng báo cáo hiện trạng môi trường, thẩm định đánh giá tác động môi trường (ĐTM) cho các dự án công nghiệp.
  • Kỹ sư công nghệ môi trường tại các nhà máy luyện kim, mạ điện và chế biến khoáng sản: Tham khảo đặc tính phân đoạn hóa học của bùn thải để cải tiến công nghệ kết tủa hóa lý, trung hòa pH và thu hồi kim loại tái chế.
  • Các tổ chức nghiên cứu bảo tồn nguồn nước và phát triển bền vững: Khai thác dữ liệu lịch sử bồi lắng trầm tích 28 phân đoạn để mô hình hóa xu hướng lan truyền chất ô nhiễm và xây dựng dự án phục hồi sinh thái sông ngòi.

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

1. Tại sao cần xác định 5 dạng liên kết kim loại thay vì chỉ phân tích hàm lượng tổng?
Hàm lượng tổng chỉ phản ánh mức độ tập trung kim loại mà không thể hiện độc tính thực tế. Các dạng trao đổi (F1) và cacbonat (F2) chiếm tới hơn 50% ở Cadimi rất dễ tan vào nước khi môi trường thay đổi, trong khi dạng cặn dư (F5) của Chì và Đồng lên tới trên 55% lại hoàn toàn trơ và an toàn đối với sinh vật.

2. Kỹ thuật GF-AAS có ưu thế gì so với F-AAS khi phân tích Cadimi?
Cadimi trong mẫu dịch chiết có nồng độ siêu vết ở mức vài microgam/lít, nằm dưới giới hạn phát hiện của kỹ thuật ngọn lửa F-AAS. Kỹ thuật lò graphit GF-AAS với chương trình gia nhiệt tro hóa 320°C và nguyên tử hóa 1650°C nâng độ nhạy lên gấp 100 đến 1000 lần, đảm bảo phép đo đạt độ chính xác cao.

3. Chỉ số rủi ro RAC phản ánh nguy cơ thực tế như thế nào tại sông Cầu?
Chỉ số RAC đo lường tổng tỷ lệ phần trăm hai dạng linh động ($F1 + F2$). Kết quả nghiên cứu ghi nhận RAC của Cd vượt ngưỡng 30% - 50%, xếp vào mức rủi ro cao đến rất cao. Điều này chứng minh Cd tại sông Cầu có khả năng phát tán mạnh vào chuỗi thức ăn thủy sản.

4. Yếu tố nào gây ra sự khác biệt về nồng độ kim loại giữa các tầng trầm tích?
Sự khác biệt phản ánh cường độ hoạt động công nghiệp qua từng thời kỳ. Các lớp mặt từ 0 đến 30 cm có nồng độ Cd và Zn cao gấp 2 đến 3 lần tầng đáy 80 - 110 cm do sự phát triển mạnh của các cơ sở luyện kim, cán thép và mạ điện tại Thái Nguyên từ những năm 2000 trở lại đây.

5. Vì sao cần khảo sát đường chuẩn trên 5 hệ nền dịch chiết riêng biệt?
Các dung môi chiết như $CH_3COONH_4$, $NH_2OH\cdot HCl$ hay cường thủy chứa nồng độ muối và axit rất cao, gây hiện tượng hấp thụ phân tử và tán xạ ánh sáng. Việc lập 5 đường chuẩn tương ứng từng hệ nền giúp loại bỏ triệt để sai số nền, đảm bảo hệ số tương quan đạt $R^2 > 0,999$.

Kết luận

  • Đề tài đã xây dựng và chuẩn hóa thành công quy trình phân tích định lượng hàm lượng tổng và 5 dạng phân đoạn của Cu, Pb, Zn, Cd trong 28 mẫu trầm tích sông Cầu bằng kỹ thuật F-AAS và GF-AAS.
  • Kết quả khẳng định Cd và Zn là hai tác nhân ô nhiễm chính với hàm lượng vượt từ 4 đến 65 lần mức nền địa chất, trong đó Cd có chỉ số RAC trên 30%, tiềm ẩn nguy cơ rủi ro sinh thái rất cao.
  • Đồng và Chì chủ yếu phân bố ở các dạng liên kết hữu cơ (F4) và cặn dư silicat (F5) bền vững với tỷ lệ trên 60%, ít có khả năng giải phóng ngược vào nguồn nước tự nhiên.
  • Biến thiên nồng độ theo chiều sâu phản ánh xu hướng tích lũy ô nhiễm gia tăng mạnh ở các lớp trầm tích tầng mặt (0 - 30 cm) do tác động của dòng thải công nghiệp luyện kim địa phương.
  • Trong giai đoạn 2027 - 2030, các cơ quan quản lý cần ưu tiên triển khai nạo vét cục bộ các điểm nóng ô nhiễm và áp dụng quy chuẩn phân tích dạng vào mạng lưới quan trắc môi trường lưu vực sông Cầu. Để khai thác chi tiết số liệu thực nghiệm và quy trình phân tích hóa học, quý độc giả có thể liên hệ tham khảo toàn văn luận văn tại thư viện chuyên ngành.