Tổng quan nghiên cứu

Trong các thủy vực tự nhiên, khoảng 50% đến 90% tổng lượng kim loại nặng có xu hướng liên kết chặt chẽ với các hạt keo lơ lửng và tích tụ lâu dài trong lớp trầm tích đáy. Lưu vực sông Nhuệ - Đáy là một trong những hệ thống thủy văn huyết mạch tại miền Bắc Việt Nam, nhưng đồng thời cũng là điểm nóng về ô nhiễm môi trường nước mặt. Khu vực hạ lưu sông Đáy, đoạn chảy qua các tỉnh Hà Nam, Nam Định và Ninh Bình, đang chịu sức ép cực kỳ nặng nề từ quá trình công nghiệp hóa và đô thị hóa. Sự hiện diện của 176 làng nghề tại tỉnh Hà Nam, 128 làng nghề tại tỉnh Nam Định, cùng 9 khu công nghiệp với tổng diện tích gần 2.000 ha trên địa bàn tỉnh Nam Định đã và đang gia tăng nguy cơ phát thải kim loại nặng độc hại ra môi trường tiếp nhận.

Mục tiêu cốt lõi của nghiên cứu là đánh giá định lượng mức độ ô nhiễm của bốn kim loại nặng chính gồm Chì (Pb), Cadimi (Cd), Đồng (Cu) và Crom (Cr) trong trầm tích mặt, đồng thời xác định mức độ rủi ro sinh thái tiềm năng tại khu vực hạ lưu sông Đáy. Phạm vi nghiên cứu trải dài trên 100 km đường sông từ huyện Kim Bảng (Hà Nam) qua Ninh Bình đến khu vực cửa Đáy (Nam Định), với chương trình khảo sát thực địa được triển khai đồng bộ tại 22 vị trí quan trắc trọng điểm vào tháng 4 năm 2018. Kết quả nghiên cứu cung cấp bộ chỉ số môi trường có độ tin cậy cao, đóng vai trò là cơ sở khoa học then chốt hỗ trợ các cơ quan quản lý nhà nước kiểm soát chất lượng môi trường nước, giảm thiểu rủi ro lan truyền độc chất qua chuỗi thức ăn thủy sinh và bảo vệ an toàn sức khỏe cộng đồng dân cư trong lưu vực.

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 hệ thống lý thuyết địa hóa học môi trường và độc học sinh thái để đánh giá sự tích lũy của kim loại nặng trong môi trường trầm tích. Hai mô hình định lượng chủ đạo được tích hợp gồm:

Chỉ số tích lũy địa chất ($I_{geo}$) do Muller đề xuất năm 1961, xác định mức độ ô nhiễm thông qua việc so sánh nồng độ kim loại quan trắc được với hàm lượng nền tự nhiên trong vỏ Trái đất, có tính đến hệ số hiệu chỉnh thạch học 1,5 nhằm loại bỏ sự biến thiên địa chất tự nhiên.

Mô hình chỉ số rủi ro sinh thái tiềm năng ($RI$) của Lars Hakanson (1980), đánh giá tổng hợp tác động sinh thái dựa trên độc tính sinh học của từng nguyên tố. Hệ số độc tính được xác định nghiêm ngặt cho từng kim loại: Cadimi có hệ số độc tính bằng 30, Chì và Đồng bằng 5, Crom bằng 2.

Bên cạnh đó, nghiên cứu sử dụng hệ số rủi ro ($RQ$) dựa trên tỷ số giữa nồng độ quan trắc thực tế ($MEC$) và nồng độ ngụ ý không gây tác động ($PNEC$), cùng với mã đánh giá rủi ro ($RAC$) của Tessier dựa trên tỷ lệ phần trăm của các dạng liên kết trao đổi và cacbonat dễ giải phóng vào nước. Khái niệm trầm tích mặt (0 - 10 cm) được định nghĩa là tầng hoạt động sinh học cao, nơi diễn ra quá trình hấp phụ, trao đổi ion và giải phóng kim loại nặng giữa pha rắn và pha lỏng.

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

Nguồn dữ liệu của luận văn được thu thập thông qua chiến dịch lấy mẫu thực địa kéo dài 3 ngày từ ngày 14/4/2018 đến ngày 16/4/2018 tại 22 vị trí quan trắc đại diện dọc theo sông Đáy. Mẫu trầm tích mặt được thu thập ở độ sâu từ 0 đến 10 cm bằng thiết bị gầu ngoạm Peterson diện tích 600 $cm^2$ theo quy chuẩn TCVN 6663-13:2015 (ISO 5667-13:2011). Song song đó, 22 mẫu nước mặt được lấy ở độ sâu 0,1 m bằng thiết bị lấy mẫu đứng Wilco và axit hóa ngay bằng $HNO_3$ đặc để duy trì $pH < 2$. Phương pháp chọn mẫu theo mặt cắt dọc dòng chảy từ thượng nguồn Kim Bảng đến cửa biển nhằm phản ánh đầy đủ quy luật lan truyền chất ô nhiễm theo thủy lực sông.

Phương pháp phân tích hóa nghiệm được thực hiện tại phòng thí nghiệm đạt chuẩn Vilas 955. Mẫu trầm tích được xử lý theo phương pháp US-EPA 3050B bằng hỗn hợp $HNO_3$ và $H_2O_2$ 30% để phá mẫu hoàn toàn, sau đó phân tích định lượng trên thiết bị quang phổ hấp thụ nguyên tử ngọn lửa (F-AAS) Thermo Fisher Solaar M6. Lý do lựa chọn F-AAS cho mẫu trầm tích là nhờ độ chính xác cao và dải đo thích hợp với nồng độ kim loại ở mức $mg/kg$ (ppm). Đối với mẫu nước, do nồng độ kim loại hòa tan rất nhỏ (thường dưới 1 $\mu g/l$), nghiên cứu áp dụng kỹ thuật hấp thụ nguyên tử không ngọn lửa lò graphit (GF-AAS) theo tiêu chuẩn SMEWW 3111 B:2017 nhằm đạt độ nhạy phát hiện siêu vết ở mức 0,1 đến 10 $ng/ml$ (ppb). Dữ liệu sau đó được xử lý thống kê, phân tích tương quan Pearson và mô hình hồi quy tuyến tính trên phần mềm SPSS 20.0.

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 trên 22 mẫu trầm tích mặt và 22 mẫu nước mặt khu vực hạ lưu sông Đáy mang lại các phát hiện khoa học quan trọng:

Thứ nhất, hàm lượng của cả 4 kim loại nặng phân bố không đồng đều và có xu hướng gia tăng rõ rệt từ thượng nguồn về phía hạ lưu và cửa biển. Nồng độ tích lũy cao nhất tập trung tại các vị trí tiếp nhận nước thải công nghiệp và giao thoa thủy văn: trạm D2 (điểm hợp lưu với sông Nhuệ), trạm D8 (giao sông Vân, gần bãi tập kết cát TP. Ninh Bình) và trạm D10 (ngã ba sông Sắt, khu vực cảng than lớn).

Thứ hai, khi đối chiếu với Quy chuẩn kỹ thuật Quốc gia về chất lượng trầm tích (QCVN 43:2012/BTNMT), đa số các mẫu đều nằm trong ngưỡng giới hạn cho phép đối với trầm tích nước ngọt và nước mặn. Tuy nhiên, tính toán chỉ số tích lũy địa chất ($I_{geo}$) chỉ ra rằng kim loại Chì ($Pb$) và Cadimi ($Cd$) đã đạt mức ô nhiễm từ nhẹ đến trung bình ($I_{geo} > 0$) tại các trạm chịu tác động mạnh từ làng nghề đúc kim loại và bến bãi giao thông thủy (D5, D8, D10).

Thứ ba, hàm lượng kim loại nặng trong trầm tích cao gấp từ 1.000 đến 10.000 lần so với nồng độ hòa tan trong nước mặt. Sự chênh lệch này khẳng định trầm tích đóng vai trò như một bể tích lũy chất ô nhiễm lâu dài của toàn bộ lưu vực.

Thứ tư, đánh giá chỉ số rủi ro sinh thái tiềm năng ($RI$) cho thấy khoảng 80% các điểm trạm có mức độ rủi ro sinh thái thấp ($RI < 110$). Mặc dù nồng độ tổng số của Cadimi ($Cd$) trong trầm tích thấp hơn nhiều so với Đồng hay Crom, nhưng $Cd$ lại chiếm tới hơn 60% tổng giá trị chỉ số rủi ro sinh thái $RI$ tại các điểm quan trắc do sở hữu hệ số độc tính sinh học cao vượt trội ($T_r = 30$).

Thảo luận kết quả

Sự biến thiên hàm lượng kim loại trong trầm tích hạ lưu sông Đáy phản ánh sự kết hợp chặt chẽ giữa các hoạt động nhân sinh và quá trình thủy động lực học. Nguyên nhân cốt lõi gây ô nhiễm $Pb$ và $Cu$ tại các trạm thượng nguồn và trung lưu xuất phát từ nước thải chưa qua xử lý triệt để của các làng nghề cơ khí, mạ điện (chiếm 11,71% tổng làng nghề Nam Định) và hoạt động khai thác đá vôi, chế biến khoáng sản tại Thanh Liêm (Hà Nam). Hoạt động giao thông thủy mật độ cao và các bãi tập kết than tại Ninh Bình cũng đóng góp đáng kể lượng bụi kim loại lắng đọng.

Về mặt thủy lực, khi dòng chảy đổ về khu vực cửa Đáy (từ trạm D19 đến D22), lưu tốc dòng nước giảm mạnh kết hợp với hiện tượng keo tụ muối khoáng khi gặp môi trường nước mặn, làm gia tăng tốc độ lắng đọng của các hạt hữu cơ mang kim loại nặng xuống đáy sông.

Dữ liệu phân bố không gian có thể được minh họa trực quan thông qua biểu đồ cột thể hiện nồng độ kim loại tại 22 trạm quan trắc kết hợp bảng ma trận tương quan Pearson. Kết quả phân tích thống kê cho thấy tương quan thuận rất chặt chẽ giữa $Cu$ và $Pb$ với hệ số tương quan $r > 0,65$, chứng minh hai kim loại này có chung nguồn gốc phát thải từ hoạt động gia công kim loại và xả thải công nghiệp.

So sánh với các nghiên cứu quốc tế, mức độ ô nhiễm tại hạ lưu sông Đáy thấp hơn đáng kể so với sông Châu Giang (Trung Quốc) nơi hàm lượng $Cr$ lên tới 1.656 $mg/kg$, nhưng có xu hướng tích lũy sinh thái tương đồng với các nghiên cứu tại cửa sông Soài Rạp và hạ lưu sông Cu Đê ở Việt Nam, nơi $Cd$ và $Cu$ luôn là các nhân tố sinh thái cần được kiểm soát hàng đầu.

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

Dựa trên kết quả định lượng mức độ ô nhiễm và đánh giá rủi ro sinh thái, bốn nhóm giải pháp quản trị môi trường cấp bách được đề xuất như sau:

Thứ nhất, siết chặt kiểm soát và quy chuẩn hóa nguồn thải công nghiệp - làng nghề. Sở Tài nguyên và Môi trường các tỉnh Hà Nam, Nam Định, Ninh Bình cần phối hợp liên tỉnh thanh tra định kỳ 176 cơ sở làng nghề và 9 khu công nghiệp trong lưu vực; đặt mục tiêu giảm thiểu từ 25% đến 30% tổng tải lượng kim loại nặng xả thải trực tiếp vào dòng sông trong giai đoạn 2024–2026.

Thứ hai, nâng cấp hạ tầng xử lý nước thải tập trung. Yêu cầu 100% các khu, cụm công nghiệp nằm dọc hành lang sông Đáy (như KCN Khánh Phú, Gián Khẩu, Hòa Mạc) phải hoàn thiện hệ thống xử lý nước thải đạt QCVN 40:2011/BTNMT trước khi xả ra môi trường, hoàn thành nâng cấp công nghệ trong thời hạn 24 tháng dưới sự giám sát của Ban Quản lý các Khu công nghiệp.

Thứ ba, tiến hành nạo vét và xử lý trầm tích tại các điểm nóng bồi lắng. Chi cục Thủy lợi và Công ty Quản lý Đô thị các địa phương cần triển khai kế hoạch nạo vét lớp trầm tích mặt từ 0 đến 15 cm tại khu vực ngã ba sông Vân (D8), cảng than Ninh Bình (D10) và vùng bồi lắng cửa Đáy trong các tháng mùa cạn (tháng 11 đến tháng 4 năm sau). Hoạt động này nhằm loại bỏ triệt để nguồn giải phóng kim loại ngược vào nước, giảm thiểu 40% nguy cơ ô nhiễm thứ cấp.

Thứ tư, thiết lập mạng lưới quan trắc tự động liên tỉnh. Bộ Tài nguyên và Môi trường cần chủ trì lắp đặt từ 5 đến 8 trạm quan trắc tự động liên tục đo đạc chất lượng nước mặt và trầm tích tại các vị trí giáp ranh địa giới hành chính trong vòng 12 đến 18 tháng, phục vụ cảnh báo sớm các sự cố ô nhiễm kim loại nặng xuyên biên giới.

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

Nội dung và kết quả nghiên cứu của luận văn mang giá trị ứng dụng thực tiễn cao cho 4 nhóm đối tượng chuyên môn:

Nhóm cán bộ quản lý môi trường và nhà hoạch định chính sách: Luận văn cung cấp cơ sở dữ liệu định lượng về tải lượng ô nhiễm 4 kim loại nặng ($Pb, Cd, Cu, Cr$) tại 22 vị trí trọng yếu, hỗ trợ đắc lực cho việc rà soát và thực thi Đề án bảo vệ môi trường lưu vực sông Nhuệ - sông Đáy do Thủ tướng Chính phủ phê duyệt.

Nhóm chuyên gia, nghiên cứu sinh và học viên ngành Khoa học Môi trường: Tài liệu là mẫu hình nghiên cứu chuẩn mực về việc tích hợp phương pháp phân tích thực nghiệm US-EPA 3050B, kỹ thuật quang phổ AAS cùng các mô hình đánh giá rủi ro sinh thái hiện đại như $I_{geo}$, $RI$, $RQ$ và $RAC$.

Nhóm ban quản lý khu công nghiệp và kỹ sư công nghệ môi trường: Cung cấp thông số nền chất lượng trầm tích để phục vụ công tác lập Báo cáo Đánh giá tác động môi trường (ĐTM) và thiết kế hệ thống xử lý nước thải công nghiệp phù hợp với đặc thù tiếp nhận của nguồn nước sông Đáy.

Nhóm đơn vị nuôi trồng và quản lý thủy sản: Giúp các cơ sở nuôi thủy sản tại khu vực hạ lưu và cửa sông (từ trạm D11 đến D18) nắm bắt được nguy cơ phơi nhiễm $Cd$ và $Pb$, từ đó xây dựng kế hoạch khoanh vùng bảo vệ an toàn cho các vùng nuôi trồng lồng bè.

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

Tại sao nghiên cứu môi trường nước lại tập trung đánh giá trầm tích đáy thay vì chỉ phân tích nước mặt?

Trầm tích đáy giữ vai trò là bể chứa lâu dài, nơi tích lũy từ 50% đến 90% tổng lượng kim loại nặng trong thủy vực thông qua các quá trình lắng đọng cơ học và hấp phụ keo. Việc phân tích tầng trầm tích mặt 0 - 10 cm phản ánh chính xác lịch sử tích lũy ô nhiễm liên tục theo thời gian, trong khi nồng độ kim loại hòa tan trong nước mặt biến động rất mạnh theo mùa và lưu lượng xả thải tức thời.

Kim loại nặng nào gây ra rủi ro sinh thái tiềm năng lớn nhất tại khu vực hạ lưu sông Đáy?

Cadimi ($Cd$) là nguyên tố đóng góp lớn nhất vào mức độ rủi ro sinh thái, chiếm trên 60% tổng chỉ số rủi ro sinh thái tiềm năng $RI$ tại các điểm quan trắc. Nguyên nhân là do $Cd$ có hệ số độc tính sinh thái rất cao ($T_r = 30$), có khả năng can thiệp vào các phản ứng enzyme tế bào và tích lũy sinh học mạnh trong gan, thận của thủy sinh vật.

Kỹ thuật quang phổ hấp thụ nguyên tử nào được sử dụng để phân tích mẫu trầm tích và mẫu nước?

Nghiên cứu ứng dụng kỹ thuật ngọn lửa (F-AAS) trên thiết bị Thermo Fisher Solaar M6 để phân tích mẫu trầm tích sau khi phá mẫu bằng axit vì hàm lượng kim loại ở ngưỡng $mg/kg$ (ppm). Đối với mẫu nước mặt, kỹ thuật không ngọn lửa lò graphit (GF-AAS) được lựa chọn để phát hiện lượng siêu vết ở cấp độ $ng/ml$ (ppb) theo tiêu chuẩn SMEWW 3111 B:2017.

Các nguồn thải nhân tạo nào đóng góp chính vào ô nhiễm kim loại nặng tại khu vực nghiên cứu?

Nguồn thải chủ yếu đến từ nước thải chưa qua xử lý của 176 làng nghề tại Hà Nam và 128 làng nghề tại Nam Định, đặc biệt là các làng nghề đúc kim loại, cơ khí và mạ điện. Ngoài ra, nước chảy tràn từ các bãi tập kết than, cát tại TP. Ninh Bình và nước thải từ 9 khu công nghiệp ven sông là nguồn gia tăng đáng kể nồng độ $Pb$ và $Cu$.

Tại sao hàm lượng kim loại nặng lại có xu hướng tăng cao tại khu vực cửa Đáy?

Tại vùng cửa sông tiếp giáp biển, tốc độ dòng chảy giảm đột ngột kết hợp với quá trình hòa trộn giữa nước ngọt và nước mặn gây ra hiện tượng keo tụ hóa lý. Quá trình này đẩy nhanh tốc độ kết tủa và lắng đọng của các hạt phù sa lơ lửng chứa kim loại nặng xuống đáy sông, làm tăng lượng tích lũy kim loại tại các trạm từ D19 đến D22.

Kết luận

  • Luận văn đã định lượng thành công hàm lượng 4 kim loại nặng ($Pb, Cd, Cu, Cr$) tại 22 vị trí quan trắc hạ lưu sông Đáy, xác định nồng độ trong trầm tích cao gấp 1.000 đến 10.000 lần so với nước mặt.
  • Chỉ số tích lũy địa chất ($I_{geo}$) xác nhận khu vực nghiên cứu chịu tác động nhiễm bẩn $Pb$ và $Cd$ từ nhẹ đến trung bình tại các điểm giao cắt thủy văn và vùng làng nghề.
  • Đánh giá rủi ro sinh thái chỉ ra Cadimi ($Cd$) là nhân tố đóng góp trên 60% tổng rủi ro sinh thái tiềm năng ($RI$), đòi hỏi phải được ưu tiên kiểm soát nghiêm ngặt.
  • Chuẩn hóa thành công quy trình quan trắc kết hợp phương pháp US-EPA 3050B với kỹ thuật phân tích F-AAS và GF-AAS đạt chuẩn kiểm định AOAC.
  • Đề xuất hệ thống 4 nhóm giải pháp quản lý đồng bộ từ xử lý nước thải tại nguồn, nạo vét 0 - 15 cm bùn đáy đến thiết lập từ 5 đến 8 trạm quan trắc tự động liên tỉnh.

Đóng góp lớn nhất của nghiên cứu là cung cấp bức tranh toàn cảnh định lượng và hệ thống bản đồ hóa về hiện trạng tích lũy kim loại nặng trong trầm tích hạ lưu sông Đáy. Kế hoạch tiếp theo cần duy trì tần suất quan trắc định kỳ 6 tháng một lần trong giai đoạn 2024–2028 để theo dõi biến động môi trường. Độc giả quan tâm và các nhà nghiên cứu hãy tham khảo toàn văn luận văn thạc sĩ khoa học môi trường để ứng dụng chi tiết phương pháp luận và số liệu vào các đề án quản trị lưu vực sông thực tiễn.