Tổng quan nghiên cứu

Khai thác và chế biến khoáng sản kim loại màu đóng vai trò quan trọng trong phát triển công nghiệp, song hoạt động này cũng phát thải hàng triệu tấn đất đá thải và bùn quặng ra môi trường. Tại Việt Nam, mỏ chì kẽm Làng Hích thuộc xã Tân Long, huyện Đồng Hỷ, tỉnh Thái Nguyên có tổng diện tích quy hoạch 320,7 ha, đi vào hoạt động khai thác từ năm 1987. Trải qua hơn 30 năm vận hành, các khu vực bãi thải rộng 3,1 ha và hệ thống suối lân cận đang chịu áp lực ô nhiễm kim loại nặng nghiêm trọng, đặc biệt là Cadmi (Cd). Cadmi là nguyên tố độc tính cao, có chu kỳ bán hủy sinh học trong cơ thể người kéo dài từ 20 đến 30 năm, gây tổn thương thận, giòn xương và rối loạn chức năng enzyme. Nhiều khảo sát trước đây ghi nhận hàm lượng Cd tại mỏ vượt quy chuẩn môi trường từ 96,6 đến 185 lần.

Luận văn thạc sĩ Hóa học của học viên Nguyễn Thị Ánh Tuyết, chuyên ngành Hóa Phân tích (mã số 8.18), thực hiện dưới sự hướng dẫn của Tiến sĩ Vương Trường Xuân và Phó Giáo sư Yang Guangming tại Trường Đại học Khoa học – Đại học Thái Nguyên. Đề tài được tài trợ bởi Bộ Giáo dục và Đào tạo theo mã số B2020-TNA-15. Nghiên cứu hướng đến ba mục tiêu cốt lõi: tối ưu hóa quy trình chiết tuần tự để phân lập 5 dạng liên kết của Cd; đánh giá quy luật phân bố hàm lượng tổng số và các dạng hóa học của Cd trên 22 mẫu đất bãi thải, trầm tích suối và đất nông nghiệp thu thập vào tháng 11/2018; đồng thời xác định mức độ rủi ro sinh thái thông qua các chỉ số tích lũy địa chất và đánh giá rủi ro chuyên sâu. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn trong việc cung cấp cơ sở dữ liệu khoa học tin cậy để địa phương quản lý rủi ro và quy hoạch sử dụng đất an toàn.

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 được xây dựng trên nền tảng lý thuyết phân tích dạng hóa học (Speciation Analysis) và độc học môi trường. Độc tính, độ linh động cũng như khả năng tích lũy sinh học của Cadmi trong đất phụ thuộc chủ yếu vào dạng tồn tại hóa học thay vì tổng nồng độ. Theo quy trình Tessier cải tiến, Cadmi được phân chia thành 5 phân đoạn liên kết:

  1. Dạng trao đổi (F1): Liên kết tĩnh điện yếu trên bề mặt khoáng sét, dễ dàng giải phóng vào pha nước khi lực ion thay đổi.
  2. Dạng liên kết carbonat (F2): Rất nhạy cảm với sự suy giảm pH của môi trường đất.
  3. Dạng liên kết sắt-mangan oxit (F3): Không bền và bị khử giải phóng ion kim loại trong điều kiện yếm khí.
  4. Dạng liên kết chất hữu cơ (F4): Bị phá vỡ khi xảy ra quá trình oxy hóa sinh hóa.
  5. Dạng cặn dư (F5): Nằm sâu trong mạng lưới cấu trúc silicat nguyên sinh bền vững, không có khả năng tan vào nước.

Để định lượng rủi ro sinh thái nhân tạo, nghiên cứu áp dụng đồng thời ba mô hình đánh giá quốc tế:

  • Chỉ số tích lũy địa chất ($I_{geo} = \log_2[C_n / (1,5 \times B_n)]$), sử dụng giá trị nền vỏ Trái đất $B_n = 0,2$ mg/kg để phân loại 7 cấp độ ô nhiễm từ không ô nhiễm ($I_{geo} \le 0$) đến rất nghiêm trọng ($I_{geo} \ge 5$).
  • Nhân tố gây ô nhiễm cá nhân ($ICF = (F1 + F2 + F3 + F4) / F5$), xác định tỷ lệ giữa các dạng linh động so với dạng trơ cặn dư.
  • Mã đánh giá rủi ro ($RAC = [(F1 + F2) / \text{Tổng}] \times 100%$), chia thành 4 mức độ từ thấp (<10%), trung bình (10-30%), cao (30-50%) đến rất cao (>50%).

Các giá trị định lượng được đối chiếu trực tiếp với Quy chuẩn kỹ thuật quốc gia QCVN 03:2008/BTNMT (quy định ngưỡng Cd trong đất nông nghiệp là 2,0 mg/kg và đất công nghiệp là 10,0 mg/kg) cùng các tiêu chuẩn quốc tế từ Hà Lan và Ba Lan.

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

Nghiên cứu khảo sát trên tập dữ liệu gồm 22 mẫu đại diện thu thập tại hiện trường tầng đất mặt ở độ sâu từ 0 đến 20 cm vào tháng 11/2018. Mạng lưới lấy mẫu bao gồm 6 mẫu bãi thải mỏ (BT1 đến BT6), 5 mẫu trầm tích suối tiếp nhận nước mỏ (TTS1 đến TTS5), 9 mẫu đất canh tác lúa và ngô xung quanh bãi thải, cùng 2 mẫu đất ruộng ven suối (R_S1 và R_S2). Tại mỗi điểm, 4 đến 5 mẫu phụ được lấy ngẫu nhiên trong bán kính 5 m, gộp thành 500 g mẫu hỗn hợp, sấy khô tự nhiên, nghiền mịn và rây qua rây 2 mm theo đúng quy chuẩn phân tích thổ nhưỡng.

Quy trình phân tích hàm lượng tổng số được thực hiện bằng hệ thống phá mẫu vi sóng MARS 6 với hỗn hợp cường thủy tỷ lệ 3:1 (15 mL HCl đặc và 5 mL $HNO_3$ đặc) trên 1,0000 g mẫu khô qua 4 giai đoạn gia nhiệt đạt đỉnh 180°C trong 10 phút. Quy trình chiết 5 phân đoạn Tessier cải tiến sử dụng các hóa chất phân tích tinh khiết cao: $CH_3COONH_4$ 1M (chiết F1), $CH_3COONH_4$ 1M axit hóa pH 5 bằng $CH_3COOH$ (chiết F2), $NH_2OH\cdot HCl$ 0,04M trong $CH_3COOH$ 25% ở 95°C (chiết F3), $CH_3COONH_4$ 3,2M trong $HNO_3$ 20% (chiết F4) và phá cặn F5 bằng hỗn hợp axit mạnh.

Hàm lượng Cd được xác định trên hệ thống khối phổ plasma cảm ứng Agilent 7900 ICP-MS đo đồng vị Cd-111 tại công suất cao tần 1600 W, lưu lượng khí mang 0,7 L/phút và dòng khí heli 4,3 mL/phút. Phương pháp đạt độ tuyến tính tuyệt đối với đường chuẩn từ 1 đến 100 ppb có phương trình $Y = 988,64X + 5,7233$ ($R^2 = 1,0000$). Giới hạn phát hiện (LOD) đạt 0,00244 ppb và giới hạn định lượng (LOQ) đạt 0,0072 ppb. Độ đúng của phương pháp được kiểm chứng qua mẫu chuẩn trầm tích quốc tế MESS-4, ghi nhận độ thu hồi đạt 92,11% (hàm lượng đo được 0,26 ± 0,04 ppm so với giá trị chứng chỉ 0,28 ppm), khẳng định độ tin cậy vượt trội của toàn bộ dữ liệu thực nghiệm.

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 22 mẫu đất và trầm tích tại mỏ chì kẽm Làng Hích chỉ ra mức độ ô nhiễm diện rộng với các số liệu cụ thể:

  • Đất bãi thải quặng chịu ô nhiễm nặng: Nồng độ Cd tổng số dao động từ 10,028 mg/kg đến 61,660 mg/kg, vượt từ 1,0 lần đến hơn 6,1 lần so với ngưỡng quy định cho đất công nghiệp tại QCVN 03:2008/BTNMT (10 mg/kg). Mẫu BT1 (khu vực trũng thấp) và mẫu BT6 (ngay miệng cống thoát thải) ghi nhận nồng độ tích lũy cao nhất toàn bãi thải.
  • Trầm tích suối là túi chứa kim loại độc hại: Toàn bộ 5 mẫu trầm tích suối tiếp nhận nước thải mỏ (TTS1 đến TTS5) có hàm lượng Cd tổng số cực cao, trải dài từ 51,536 mg/kg đến 84,695 mg/kg, vượt quy chuẩn đất công nghiệp từ 5,1 đến 8,4 lần và cao hơn cả nhiều vị trí trên bãi thải.
  • Lan truyền mạnh sang đất canh tác nông nghiệp: Nồng độ Cd trong đất trồng lúa và ngô gần bãi thải nằm trong khoảng 1,107 mg/kg đến 3,111 mg/kg, trong đó có 6 trên 8 mẫu (tỷ lệ 75%) vượt quá giới hạn cho phép của đất nông nghiệp (1,5 - 2,0 mg/kg). Đặc biệt, hai mẫu ruộng ven suối R_S1 và R_S2 có hàm lượng Cd tăng vọt lên mức 6,582 mg/kg và 13,445 mg/kg, vượt từ 3,3 đến 6,7 lần ngưỡng an toàn.
  • Ưu thế của các dạng linh động sinh học: Tỷ lệ thu hồi tổng 5 phân đoạn so với hàm lượng phá mẫu tổng đạt từ 87,30% đến 118,68%. Tỷ trọng hai dạng linh động cao nhất là dạng trao đổi (F1) và dạng carbonat (F2) chiếm từ 20% đến trên 45% tổng hàm lượng Cd tại nhiều vị trí đất ruộng và trầm tích, tạo ra chỉ số rủi ro sinh thái RAC vượt mức 30% đến 50%.

Thảo luận kết quả

Sự tích tụ nồng độ Cadmi cao đột biến tại bãi thải và lòng suối bắt nguồn trực tiếp từ hoạt động khai thác, tuyển rửa quặng sunfua (grenokit CdS thường cộng sinh cùng sphalerit ZnS). Quá trình phong hóa oxy hóa ngoài trời kết hợp nước mưa đã hòa tan và rửa trôi các hạt khoáng mịn xuống các vùng trũng thấp như vị trí lấy mẫu BT1 và BT6. Lượng nước thải xả trực tiếp không qua xử lý hoàn chỉnh ra dòng suối địa phương qua hàng chục năm đã biến hệ thống suối thành bồn trĩu lắng đọng trầm tích ô nhiễm nghiêm trọng.

Hiện tượng đất ruộng lúa ven suối R_S1 và R_S2 bị ô nhiễm nặng hơn gấp nhiều lần so với đất ruộng sát bãi thải được minh chứng bởi hai con đường phát tán:

  1. Thói quen canh tác lâu năm của nông dân địa phương khi bơm dẫn trực tiếp nguồn nước mặt từ dòng suối ô nhiễm vào đồng ruộng để tưới tiêu.
  2. Tuyến đường vận chuyển quặng chạy sát các thửa ruộng này, khiến bụi quặng chứa Cd rơi vãi liên tục trong thời gian dài.

Khi trực quan hóa dữ liệu qua biểu đồ cột xếp chồng (stacked bar chart) phân bố phần trăm 5 dạng hóa học, các mẫu đất ruộng thể hiện hàm lượng dạng cặn dư F5 rất thấp, trong khi dạng liên kết hữu cơ F4 và carbonat F2 chiếm tỷ lệ chủ đạo. Chỉ số tích lũy địa chất $I_{geo}$ tại bãi thải và suối đều đạt mức cấp 5 và cấp 6, thuộc nhóm ô nhiễm nghiêm trọng. Đồng thời, hệ số ô nhiễm cá nhân $ICF > 6$ và mã rủi ro $RAC > 30%$ khẳng định Cadmi tại khu vực Làng Hích có độ hoạt động hóa học cao, dễ dàng chuyển dịch vào chuỗi thức ăn thông qua hệ rễ cây lúa và ngô.

Kết quả này hoàn toàn nhất quán với các nghiên cứu của Chu Thị Thu Hà (hàm lượng Cd đạt 11,34 - 61,04 mg/kg) và Bùi Thị Kim Anh (hàm lượng Cd dao động 67,4 - 80,9 mg/kg tại bãi thải Làng Hích). So sánh quốc tế cho thấy nồng độ Cd linh động tại đây tương đương với các điểm nóng khai khoáng tại Vũ Hán, Trung Quốc (dạng F1 chiếm 19,8% - 25,2%) và Ba Lan, cho thấy tính cấp bách trong việc can thiệp xử lý môi trường.

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

Dựa trên các phát hiện thực nghiệm về mức độ ô nhiễm và dạng tồn tại của Cadmi, các biện pháp can thiệp cần được triển khai đồng bộ:

  1. Nâng cấp hệ thống xử lý nước thải mỏ đạt chuẩn xả thải: Ban lãnh đạo Công ty TNHH MTV Kim loại màu Thái Nguyên cần đầu tư xây dựng hệ thống lắng kết tủa hóa lý kết hợp hồ sinh học lọc bùn đáy trước quý IV/2021. Mục tiêu kiểm soát 100% dòng nước thải phát sinh từ khu vực chế biến đạt nồng độ Cd dưới 0,05 mg/L theo QCVN 40:2011/BTNMT trước khi xả ra lưu vực suối.
  2. Cô lập và phủ xanh sinh học bãi thải quặng: Sở Tài nguyên và Môi trường tỉnh Thái Nguyên chủ trì phối hợp với doanh nghiệp xây dựng bờ kè bê tông kiên cố ngăn rửa trôi trên diện tích 3,1 ha bãi thải trong giai đoạn 2021–2023. Đồng thời, triển khai trồng các loài thực vật siêu tích lũy kim loại như cỏ Vetiver và cây sậy nhằm cố định sinh học (phytostabilization), giảm tối thiểu 60% lượng bụi và ion kim loại phát tán ra khu dân cư.
  3. Chuyển đổi cơ cấu cây trồng và nguồn nước tưới: Ủy ban nhân dân huyện Đồng Hỷ ban hành văn bản cấm tuyệt đối việc sử dụng nguồn nước suối quanh mỏ để tưới tiêu hoa màu; đồng thời hỗ trợ người dân chuyển đổi 100% diện tích đất ruộng ven suối (khoảng 15 ha) sang trồng cây lấy gỗ, cây cảnh hoặc cây công nghiệp không dùng làm thực phẩm trong vòng 12 tháng.
  4. Thiết lập trạm quan trắc định kỳ dạng hóa học kim loại: Chi cục Bảo vệ Môi trường tỉnh Thái Nguyên định kỳ lấy mẫu và phân tích dạng hóa học Cd 2 lần/năm tại 20 điểm trọng yếu quanh khu vực mỏ, công bố rộng rãi chỉ số rủi ro sinh thái RAC cho người dân trước ngày 30/6 hàng năm.

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

Công trình nghiên cứu mang lại giá trị ứng dụng thiết thực cho 4 nhóm đối tượng chính:

  • Học viên, nghiên cứu sinh và giảng viên chuyên ngành Hóa Phân tích, Hóa Môi trường: Tài liệu cung cấp quy trình chiết tuần tự 5 phân đoạn Tessier cải tiến chuẩn hóa và kỹ thuật định lượng vết Cd bằng ICP-MS với giới hạn phát hiện đạt mức dưới ppb ($LOD = 0,00244$ ppb).
  • Cán bộ quản lý tại Sở Tài nguyên và Môi trường, Chi cục Bảo vệ Môi trường: Ứng dụng phương pháp luận kết hợp ba chỉ số $I_{geo}$, $ICF$ và $RAC$ để đánh giá rủi ro sinh thái thực tế tại các điểm khai khoáng, thay thế cho cách tiếp cận hàm lượng tổng đơn thuần trong cấp phép xả thải.
  • Doanh nghiệp khai thác và chế biến khoáng sản kim loại màu: Cung cấp số liệu đầu vào chuẩn xác để lập báo cáo đánh giá tác động môi trường (ĐTM), thiết kế quy hoạch bãi thải và phương án phục hồi môi trường sau đóng cửa mỏ.
  • Cán bộ khuyến nông và chính quyền cơ sở huyện Đồng Hỷ: Cơ sở dữ liệu khoa học tin cậy giúp địa phương hoạch định ranh giới vùng đất an toàn cho sản xuất nông nghiệp, bảo vệ sức khỏe cộng đồng nông thôn.

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

Tại sao phân tích dạng hóa học của Cadmi lại quan trọng hơn việc chỉ đo tổng hàm lượng?
Hàm lượng tổng chỉ phản ánh mức độ hiện diện thô của kim loại. Phân tích dạng hóa học giúp tách biệt các phân đoạn linh động sinh học (F1, F2 chiếm từ 20% đến 45% tại Làng Hích) có khả năng hòa tan gây độc trực tiếp cho cây trồng, khác biệt hoàn toàn với dạng cặn dư bền vững F5 không gây hại tức thời.

Hàm lượng Cadmi tại mỏ chì kẽm Làng Hích vượt quy chuẩn cho phép ở mức độ nào?
Nồng độ Cd tổng số trong đất bãi thải dao động 10,03 - 61,66 mg/kg (vượt quy chuẩn đất công nghiệp 1,0 - 6,1 lần); trầm tích suối đạt 51,54 - 84,70 mg/kg (vượt 5,1 - 8,4 lần); đất ruộng ven suối đạt 6,58 - 13,45 mg/kg (vượt quy chuẩn đất nông nghiệp 3,3 - 6,7 lần).

Phương pháp chiết tuần tự Tessier cải tiến trong luận văn gồm những bước nào?
Quy trình gồm 5 giai đoạn liên tiếp trên 1 g mẫu: (F1) chiết dạng trao đổi bằng $CH_3COONH_4$ 1M; (F2) chiết dạng carbonat bằng $CH_3COONH_4$ pH 5; (F3) chiết dạng Fe-Mn oxit bằng $NH_2OH\cdot HCl$ 0,04M ở 95°C; (F4) chiết dạng hữu cơ bằng $CH_3COONH_4$ trong $HNO_3$ 20%; (F5) phá cặn bằng cường thủy.

Nguyên nhân nào khiến đất ruộng lúa ven suối bị ô nhiễm nặng hơn ruộng sát bãi thải?
Đất ruộng ven suối chịu tác động kép từ việc nông dân dẫn trực tiếp nguồn nước suối ô nhiễm (trầm tích chứa tới 84,70 mg/kg Cd) để tưới tiêu trong nhiều năm, kết hợp với bụi quặng rơi vãi từ các đoàn xe vận tải khoáng sản chạy dọc tuyến đường dân sinh lân cận.

Thiết bị ICP-MS mang lại ưu điểm gì trong phép đo vết kim loại Cadmi?
Hệ thống Agilent 7900 ICP-MS phân tích đồng vị Cd-111 có độ nhạy vượt trội, đạt giới hạn phát hiện cực thấp $LOD = 0,00244$ ppb, khoảng tuyến tính hoàn hảo ($R^2 = 1,0000$) và độ thu hồi trên mẫu chuẩn MESS-4 đạt 92,11%, hạn chế tối đa sai số nền so với phương pháp F-AAS truyền thống.

Kết luận

  • Luận văn xây dựng thành công quy trình phân tích dạng hóa học của Cadmi bằng phương pháp chiết Tessier cải tiến kết hợp định lượng ICP-MS với độ tin cậy cao, độ thu hồi đạt 92,11% trên mẫu chuẩn quốc tế MESS-4.
  • Đất bãi thải mỏ Làng Hích ghi nhận ô nhiễm nặng với hàm lượng Cd tổng số từ 10,028 mg/kg đến 61,660 mg/kg, trong khi trầm tích suối đạt đỉnh ô nhiễm lên tới 84,695 mg/kg.
  • Đã phát hiện 75% số mẫu đất ruộng canh tác gần bãi thải và 100% mẫu ruộng ven suối vượt ngưỡng quy chuẩn kỹ thuật quốc gia QCVN 03:2008/BTNMT cho đất nông nghiệp.
  • Các chỉ số $I_{geo}$, $ICF$ và $RAC$ đồng loạt cảnh báo mức độ rủi ro sinh thái cao đến rất nghiêm trọng do tỷ lệ dạng trao đổi và carbonat linh động chiếm tỷ trọng lớn trong đất.
  • Nghiên cứu đề xuất lộ trình xử lý nước thải mỏ hoàn thành trước quý IV/2021 và chuyển đổi 15 ha cơ cấu cây trồng nông nghiệp ven suối trong vòng 12 tháng.

Công trình đóng góp bộ dữ liệu thực nghiệm toàn diện đầu tiên về dạng hóa học của Cadmi tại mỏ kẽm chì Làng Hích, mở đường cho các giải pháp cố định hóa sinh học trong tương lai. Quý độc giả và các nhà khoa học có thể liên hệ Trung tâm Học liệu và Công nghệ thông tin – Đại học Thái Nguyên để tiếp cận toàn văn công trình nghiên cứu giá trị này.