Tổng quan nghiên cứu

Hoạt động khai khoáng kim loại màu tại Việt Nam đang đặt ra những thách thức môi trường cấp bách đối với hệ sinh thái đất và nguồn nước mặt. Tại tỉnh Thái Nguyên, mỏ chì kẽm Làng Hích thuộc huyện Đồng Hỷ là trung tâm khai khoáng lâu đời với lịch sử hơn 100 năm kể từ năm 1905. Các khảo sát trước đây ghi nhận biểu hiện ô nhiễm kim loại nặng đáng báo động khi hàm lượng asen trong đất mỏ từng vượt quy chuẩn tới 78 lần, đi kèm hàm lượng kẽm vượt 48 lần và cadimi vượt 96,6 lần. Vấn đề cốt lõi là việc chỉ đo lường hàm lượng tổng số chưa thể phản ánh chính xác độc tính sinh học cũng như mức độ linh động thực tế của độc chất vào chuỗi thức ăn.

Xuất phát từ thực tiễn đó, nghiên cứu tập trung phân tích các dạng liên kết hóa học của asen, từ đó xác định quy luật phân bố không gian, cơ chế giải phóng và mức độ rủi ro sinh thái của nguyên tố này. Đối tượng khảo sát gồm hệ thống mẫu đất bãi thải, đất canh tác nông nghiệp và trầm tích suối thoát nước mỏ tại xã Tân Long, huyện Đồng Hỷ, được thu thập vào tháng 11 năm 2018 tại tầng mặt từ 0 đến 20 cm. Ý nghĩa của đề tài được lượng hóa rõ nét thông qua việc tích hợp chỉ số tích lũy địa chất Igeo và chỉ số mã đánh giá rủi ro sinh thái RAC. Kết quả tạo lập bộ dữ liệu thực nghiệm chuẩn xác với độ tin cậy đạt 95%, hỗ trợ đắc lực cho các cơ quan quản lý trong việc hoạch định chính sách kiểm soát ô nhiễm và phục hồi hệ sinh thái nông nghiệp địa phươ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 dựa trên lý thuyết phân tích dạng hóa học, khẳng định độc tính sinh học và khả năng di động của kim loại nặng trong đất phụ thuộc vào dạng liên kết hơn là hàm lượng tổng số. Mô hình chiết tuần tự Tessier cải tiến phân chia asen thành 5 dạng liên kết: dạng trao đổi (F1) liên kết lỏng lẻo dễ hòa tan; dạng liên kết cacbonat (F2) nhạy cảm với sự suy giảm pH; dạng liên kết Fe-Mn oxit (F3) dễ bị khử khi thế oxy hóa khử hạ thấp; dạng liên kết chất hữu cơ (F4) nhạy cảm với điều kiện oxy hóa; và dạng cặn dư (F5) nằm trong cấu trúc mạng tinh thể silicat bền vững.

Về độc tính học, asen vô cơ độc hơn dạng hữu cơ khoảng 100 lần, và dạng As(III) độc gấp 60 lần As(V). Sự suy giảm thế oxy hóa khử từ 500 mV xuống 200 mV làm tăng độ hòa tan của asen lên gấp 13 lần. Đề tài ứng dụng chỉ số tích lũy địa chất Igeo theo Muller với hàm lượng nền tự nhiên trong vỏ Trái Đất là 1,5 mg/kg, kết hợp chỉ số mã rủi ro sinh thái RAC tính bằng tổng phần trăm hai dạng linh động F1 và F2 trên tổng số để định lượng nguy cơ môi trường.

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

Cỡ mẫu nghiên cứu gồm 10 mẫu đại diện, ký hiệu từ S1 đến S10, thu thập vào tháng 11 năm 2018 tại xã Tân Long. Phương pháp chọn mẫu là lấy mẫu hỗn hợp ngẫu nhiên tầng mặt ở độ sâu từ 0 đến 20 cm theo nguyên tắc đường chéo, gộp từ 4 đến 5 điểm phụ trong bán kính 5 mét với khối lượng 500 gam mỗi điểm, sau đó phơi khô tự nhiên và rây qua cỡ lỗ 2 mm để bảo đảm tính đại diện của đối tượng nghiên cứu.

Quy trình phá mẫu xác định hàm lượng tổng sử dụng lò vi sóng Mars 6 ở nhiệt độ 170 độ C và công suất 1200 W trong hỗn hợp axit nitric và axit clohydric đặc. Quy trình chiết tuần tự Tessier sử dụng muối amoni axetat 1M, axit axetic chỉnh pH về mức 5, dung dịch hydroxylamin clohydrit 0,04M trong axit axetic 25% ở nhiệt độ 95 độ C, và amoni axetat 3,2M trong axit nitric 20%.

Phương pháp phân tích công cụ được lựa chọn là phổ khối plasma cảm ứng cao tần ICP-MS trên thiết bị Nexion 2000. Kỹ thuật ICP-MS được lựa chọn nhờ độ nhạy vượt trội, khả năng phân tích lượng siêu vết ở đồng vị As-75, với giới hạn phát hiện LOD đạt 0,001 ppb và giới hạn định lượng LOQ đạt 0,003 ppb. Độ tin cậy được kiểm soát qua mẫu chuẩn trầm tích quốc tế MESS-4 cho hiệu suất thu hồi đạt 92,3%, cùng độ thu hồi tổng năm phân đoạn đạt từ 93,11% đến 121,83%.

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

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

Nghiên cứu ghi nhận bốn phát hiện khoa học nổi bật về hiện trạng địa hóa tại khu vực mỏ:

Thứ nhất, hàm lượng asen tổng số dao động trong biên độ rất rộng từ 6,1 mg/kg đến 78,3 mg/kg. Mẫu trầm tích suối S8 ghi nhận nồng độ cao nhất đạt 78,3 mg/kg, vượt 4,6 lần ngưỡng quy chuẩn kỹ thuật quốc gia về chất lượng trầm tích nước ngọt QCVN 43-MT:2017/BTNMT là 17 mg/kg, đồng thời vượt 3,1 lần giới hạn đất công nghiệp. Mẫu bãi thải S1 đạt 42,7 mg/kg, vượt 1,7 lần giới hạn đất công nghiệp theo QCVN 03-MT:2015/BTNMT là 25 mg/kg.

Thứ hai, các mẫu đất ruộng ven suối (S9 và S10) có hàm lượng asen lần lượt là 11,4 mg/kg và 11,8 mg/kg, cao hơn từ 48% đến 93% so với các thửa ruộng gần bãi thải khô (S2 đến S6 dao động từ 6,1 mg/kg đến 7,7 mg/kg), dù tất cả mẫu đất nông nghiệp vẫn dưới ngưỡng quy chuẩn 15 mg/kg.

Thứ ba, dạng cặn dư F5 chiếm ưu thế ở các mẫu từ S1 đến S8 với tỷ lệ từ 40,3% đến 66,0%, nhưng suy giảm mạnh xuống mức 19,3% ở mẫu S10 và 20,8% ở mẫu S9.

Thứ tư, dạng liên kết cacbonat F2 chiếm tỷ trọng cao bất thường, dao động từ 23,8% đến 41,4%, đạt đỉnh tại mẫu đất ruộng S10 với 41,4% và bãi thải S1 với 37,0%. Dạng liên kết Fe-Mn oxit F3 cũng tăng mạnh từ mức 7,0% ở bãi thải lên tới 39,6% tại mẫu ruộng ven suối S9.

Thảo luận kết quả

Sự tích tụ asen cực đại tại trầm tích suối S8 (78,3 mg/kg), vượt xa bãi thải trên cạn S1 (42,7 mg/kg), xuất phát từ việc nước thải tuyển khoáng xả thẳng ra tự nhiên sau khi chỉ qua lắng lọc sơ bộ, khiến hạt mịn giàu asen lắng đọng lâu năm tại vùng trũng thủy văn. Đất ruộng S9 và S10 tích tụ asen cao hơn nhóm S2 đến S6 là hệ quả từ tập quán lấy nước suối tưới tiêu kết hợp bụi quặng phát tán từ xe vận chuyển. Tỷ lệ dạng liên kết cacbonat F2 chiếm tới 41,4% phản ánh sát đặc trưng địa chất mỏ Làng Hích cấu tạo từ nham tầng đá vôi cacbonat, phù hợp với nhận định của một nghiên cứu gần đây về quặng chì kẽm cộng sinh cacbonat tại miền Bắc Việt Nam.

So với các vùng mỏ lân cận, nồng độ asen tại Làng Hích (6,1 đến 78,3 mg/kg) thấp hơn mỏ chì kẽm Chợ Đồn tại tỉnh Bắc Kạn (489 đến 5630 mg/kg) và mỏ Núi Pháo (13,55 đến 122,1 mg/kg). Dù vậy, rủi ro tiềm ẩn tại Làng Hích lại rất đáng ngại vì chỉ số RAC ở đa số các mẫu đều vượt mức 30%, thuộc phân cấp rủi ro sinh thái cao. Dạng F2 rất kém bền, sẽ lập tức giải phóng asen hòa tan khi môi trường đất bị chua hóa.

Để trực quan hóa các dữ liệu này, biểu đồ cột xếp chồng 100% là phương án tối ưu giúp thể hiện rõ nét sự dịch chuyển tỷ trọng của 5 dạng liên kết qua từng vị trí. Đồng thời, biểu đồ mạng nhện đa trục sẽ đối sánh hiệu quả mức độ rủi ro giữa hai chỉ số Igeo và RAC, trong khi một bảng ma trận tương quan địa hóa sẽ minh chứng mối liên hệ thuận giữa độ chua của đất và khả năng giải phóng ion asen tự do.

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

Dựa trên kết quả định lượng, bốn giải pháp hành động cụ thể được đề xuất nhằm kiểm soát triệt để nguy cơ ô nhiễm:

Thứ nhất, nâng cấp công nghệ xử lý nước thải tuyển mỏ bằng hệ thống keo tụ tạo bông hóa học kết hợp oxy hóa xúc tác, nhằm hạ nồng độ asen trong nước thải trước khi xả ra môi trường xuống dưới ngưỡng 0,05 mg/L theo quy chuẩn nước thải công nghiệp loại A. Thời hạn hoàn thành trong vòng 12 tháng, do Ban Giám đốc Công ty Cổ phần Kim loại màu Thái Nguyên chủ trì thi công và vận hành.

Thứ hai, thiết lập hành lang cách ly nguồn nước tưới nông nghiệp tại xã Tân Long, đồng thời xây dựng 2 hồ lắng sinh học đa ngăn với tổng dung tích 5000 m3 để xử lý sơ bộ nguồn nước trước khi đưa vào kênh mương nội đồng. Thời gian triển khai trong 18 tháng, dưới sự điều phối của Phòng Tài nguyên và Môi trường huyện Đồng Hỷ cùng Ủy ban nhân dân xã Tân Long.

Thứ ba, áp dụng biện pháp ổn định hóa học đất canh tác bằng cách bón vôi nông nghiệp kết hợp khoáng bentonite định kỳ với liều lượng từ 2 đến 3 tấn trên một hecta mỗi năm, nhằm duy trì độ pH đất luôn ổn định trên mức 6,5, qua đó cố định vững chắc trên 85% lượng asen ở dạng liên kết cacbonat F2, ngăn chặn nguy cơ xâm nhập vào nông sản. Kế hoạch được thực hiện từ năm 2022 đến 2025, do Trạm Khuyến nông huyện Đồng Hỷ hướng dẫn trực tiếp cho các hộ nông dân.

Thứ tư, triển khai mô hình phục hồi sinh thái bằng thực vật thông qua việc trồng cây dương xỉ Pteris vittata dọc 3 km hành lang bờ suối và chân bãi thải, nhằm hấp thu và làm giảm từ 30% đến 45% lượng asen di động trong tầng đất mặt sau 24 tháng theo dõi. Dự án do Sở Khoa học và Công nghệ tỉnh Thái Nguyên cấp kinh phí và phối hợp cùng các viện nghiên cứu chuyên ngành.

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

Công trình nghiên cứu mang giá trị học thuật và ứng dụng cao đối với bốn nhóm đối tượng chuyên môn:

Nhóm thứ nhất là các nhà khoa học và học viên cao học chuyên ngành Hóa phân tích, Khoa học Môi trường và Địa hóa học: Luận văn cung cấp quy trình chiết tuần tự 5 bước cải tiến Tessier chuẩn xác, bộ thông số vận hành thiết bị ICP-MS tối ưu với giới hạn phát hiện 0,001 ppb, cùng phương pháp đánh giá độ tin cậy bằng mẫu chuẩn MESS-4 đạt hiệu suất 92,3%.

Nhóm thứ hai là các kỹ sư môi trường và nhà quản trị mỏ khoáng sản: Đây là tài liệu kỹ thuật giá trị để đánh giá tác động môi trường thực chất, nhận diện cơ chế rủi ro của asen dạng cacbonat chiếm tới 41,4%, từ đó thiết kế quy trình xử lý nước thải và bãi lưu giữ quặng đuôi an toàn hơn.

Nhóm thứ ba là cơ quan quản lý nhà nước như Sở Tài nguyên và Môi trường tỉnh Thái Nguyên: Kết quả cung cấp bộ chỉ số Igeo và RAC định lượng trên 10 vị trí quan trắc, làm căn cứ khoa học vững chắc để quy hoạch vùng canh tác an toàn và xây dựng kế hoạch bảo vệ lưu vực sông Cầu.

Nhóm thứ tư là cán bộ khuyến nông và cộng đồng người dân tại xã Tân Long: Giúp nhận biết rõ nguy cơ khi dùng nước suối nhiễm asen để canh tác lúa, từ đó chủ động bón vôi cân bằng pH đất và lựa chọn giống cây trồng ít mẫn cảm nhằm bảo vệ sức khỏe cộng đồng.

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

Phân tích dạng hóa học của asen có ưu thế vượt trội gì so với việc đo tổng hàm lượng?
Tổng hàm lượng chỉ phản ánh mức tích tụ mà không thể hiện độc tính thực tế. Trong nghiên cứu này, dạng cặn dư F5 chiếm đến 66,0% hoàn toàn bất động và ít độc, trong khi dạng cacbonat F2 chiếm tới 41,4% lại rất dễ hòa tan khi đất bị chua hóa. Phân tích dạng giúp nhận diện chính xác mức độ độc hại thực chất đối với hệ sinh thái.

Độ nhạy và giới hạn định lượng của thiết bị ICP-MS trong nghiên cứu đạt mức cụ thể như thế nào?
Thiết bị ICP-MS Nexion 2000 được thiết lập đường chuẩn tuyến tính với hệ số xác định R2 bằng 1 trên dải nồng độ từ 1 đến 100 ppb. Phương pháp đạt giới hạn phát hiện LOD là 0,001 ppb và giới hạn định lượng LOQ là 0,003 ppb, cho phép định lượng chính xác hàm lượng vết và siêu vết của nguyên tố asen.

Vì sao mẫu trầm tích lòng suối S8 lại có mức ô nhiễm asen cao hơn nhiều so với đất bãi thải S1?
Vị trí suối S8 nằm ở địa hình trũng thấp, tiếp nhận trực tiếp dòng nước thải mỏ chưa qua xử lý hóa học triệt để. Quá trình bồi tụ cơ học kéo dài nhiều năm khiến bùn thải lắng đọng, nâng hàm lượng asen tại S8 lên tới 78,3 mg/kg, vượt xa mức 42,7 mg/kg ghi nhận tại bãi thải khô trên bờ S1.

Đất ruộng lúa xung quanh khu mỏ Làng Hích đã bị ô nhiễm asen theo quy chuẩn Việt Nam hay chưa?
Theo quy chuẩn kỹ thuật quốc gia QCVN 03-MT:2015/BTNMT về hàm lượng asen trong đất nông nghiệp là 15 mg/kg, tất cả 7 mẫu ruộng lúa khảo sát đều an toàn với hàm lượng từ 6,1 đến 11,8 mg/kg. Tuy nhiên, chỉ số rủi ro sinh thái RAC của các ruộng lúa lại ở mức trung bình đến cao, từ 24,1% đến 41,8%.

Chỉ số mã đánh giá rủi ro sinh thái RAC cảnh báo điều gì về môi trường đất tại mỏ Làng Hích?
Chỉ số RAC phản ánh tổng tỷ lệ phần trăm hai dạng linh động F1 và F2. Tại Làng Hích, có tới 6 trên 10 mẫu khảo sát vượt ngưỡng RAC 30%, thuộc phân cấp rủi ro sinh thái cao. Kết quả này cảnh báo asen rất dễ bị giải phóng hàng loạt vào nguồn nước nếu độ kiềm của đất và nước mặt suy giảm.

Kết luận

Luận văn đã giải quyết trọn vẹn các mục tiêu khoa học đề ra thông qua các kết quả cốt lõi:

  • Xây dựng thành công quy trình chiết tuần tự 5 dạng liên kết asen Tessier cải tiến bằng kỹ thuật ICP-MS với độ thu hồi đạt 92,3% trên mẫu chuẩn MESS-4.
  • Định lượng hàm lượng asen tổng số trong 10 mẫu dao động từ 6,1 đến 78,3 mg/kg, xác định các điểm nóng ô nhiễm tại đất bãi thải S1 và trầm tích suối S7, S8 vượt chuẩn quốc gia từ 1,7 đến 4,6 lần.
  • Khám phá đặc tính phân bố dạng hóa học với dạng liên kết cacbonat F2 chiếm tỷ trọng cao bất thường từ 23,8% đến 41,4%, là nguồn phát tán asen tiềm tàng khi đất bị axit hóa.
  • Đánh giá toàn diện rủi ro sinh thái qua chỉ số Igeo từ 1,2 đến 5,0 và chỉ số RAC vượt mức 30% ở đa số mẫu, khẳng định mức độ rủi ro môi trường ở cấp độ cao.
  • Đề xuất lộ trình xử lý môi trường giai đoạn 2022 đến 2025 kết hợp xử lý nước thải mỏ dưới 0,05 mg/L, bón vôi duy trì pH trên 6,5 và trồng cây dương xỉ hấp phụ kim loại.

Đóng góp nổi bật của luận văn là cung cấp bức tranh địa hóa môi trường chuẩn xác, khẳng định tính tất yếu của việc phân tích dạng trong quản lý rủi ro độc chất mỏ. Độc giả và các đơn vị nghiên cứu hãy cùng kết nối, chia sẻ và áp dụng ngay các giải pháp khoa học này vào thực tiễn bảo vệ tài nguyên môi trường đất bền vững!