Tổng quan nghiên cứu

Lưu vực sông Nhuệ giữ vai trò huyết mạch trong hệ thống thủy lợi và thoát nước của vùng châu thổ sông Hồng với tổng diện tích tự nhiên 107.530 ha, trải dài 74 km qua địa phận thành phố Hà Nội (chiếm 82% diện tích) và tỉnh Hà Nam (chiếm 18% diện tích). Áp lực từ tốc độ đô thị hóa nhanh cùng sự phát triển mạnh mẽ của các cụm công nghiệp, làng nghề đã biến dòng sông thành nơi tiếp nhận hơn 450.000 m³/ngày đêm nước thải chưa qua xử lý đầy đủ, làm suy giảm đến 71% chất lượng môi trường nước. Trong các nhóm chất độc hại, kim loại nặng nổi lên như một hiểm họa tiềm tàng do tính bền vững cao, khả năng tích lũy sinh học và xu hướng lắng đọng sâu trong nền bùn đáy.

Vấn đề cốt lõi đặt ra là sự phân bố ô nhiễm không đồng đều theo dòng chảy và nguy cơ giải phóng ngược kim loại nặng vào cột nước khi điều kiện hóa lý thay đổi. Đề tài tập trung xác định hàm lượng tổng số của ba nguyên tố đồng (Cu), chì (Pb), kẽm (Zn) và làm rõ các dạng liên kết hóa học, khả năng hấp phụ đẳng nhiệt của lớp bùn đáy tầng mặt (0–30 cm). Nghiên cứu được triển khai thực địa từ tháng 11/2011 đến tháng 11/2012 tại 8 điểm quan trắc chiến lược trải dài từ Cống Liên Mạc (thượng nguồn) đến Cống Phủ Lý (hạ lưu). Kết quả đo đạc cho thấy nồng độ Zn đạt cực đại 1226 mg/kg và Pb đạt 156 mg/kg, vượt giới hạn an toàn sinh thái của Canada (PEL 2002) lần lượt là 3,89 lần và 1,70 lần. Nghiên cứu cung cấp hệ thống dữ liệu định lượng chuẩn xác, hỗ trợ đắc lực cho các cơ quan quản lý trong việc đánh giá rủi ro sinh thái và thiết lập quy hoạch kiểm soát ô nhiễm liên tỉnh.

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 dựa trên nền tảng lý thuyết hấp phụ bề mặt và mô hình đẳng nhiệt hấp phụ Freundlich ($Q_s = K_F \cdot C_e^n$). Mô hình này mô tả tương tác dị thể giữa các ion kim loại hòa tan và bề mặt hạt keo trầm tích, trong đó hệ số $K_F$ biểu thị dung lượng hấp phụ và chỉ số $n$ phản ánh cường độ ái lực liên kết. Cơ chế lưu giữ và chuyển hóa chất ô nhiễm còn được soi chiếu qua lý thuyết "Quả bom hóa học nổ chậm" (Chemical Time Bomb – CTB do Stigliani đề xuất năm 1991), cảnh báo nguy cơ kim loại nặng bị giải phóng ồ ạt khi độ chua hoặc điện thế oxy hóa khử (Eh) của môi trường biến động.

Bên cạnh đó, khung nghiên cứu tích hợp lý thuyết phân đoạn hóa học (chemical speciation) theo quy trình chiết phân đoạn liên tiếp. Khung lý thuyết này định nghĩa năm dạng tồn tại chính của kim loại: dạng trao đổi (F1), dạng liên kết cacbonat (F2), dạng liên kết oxit sắt - mangan (F3), dạng liên kết hữu cơ (F4) và dạng liên kết trong mạng tinh thể khoáng sét (F5). Dung tích trao đổi cation (CEC) cùng thành phần cơ giới hạt mịn đóng vai trò là các biến số trung gian kiểm soát tốc độ di chuyển và độ linh động sinh học của các cation kim loại nặng hóa trị hai.

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

  • Cỡ mẫu và thiết kế thu mẫu: Nghiên cứu chọn cỡ mẫu $N = 8$ vị trí đại diện dọc theo toàn tuyến sông Nhuệ, phân bố từ cống Liên Mạc, Cầu Diễn, Cầu Trắng, Cầu Tả Thanh Oai, Cầu Chiếc, Đập Đồng Quan, Cầu Nhật Tựu đến Cống Phủ Lý. Phương pháp lấy mẫu hỗn hợp tầng mặt (độ sâu 0–30 cm) được thực hiện bằng gầu ngoạm Wildco chuyên dụng, tuân thủ nghiêm ngặt tiêu chuẩn quốc gia TCVN 6631-13:2000. Cỡ mẫu này bảo đảm bao quát đầy đủ các vùng hợp lưu tiếp nhận nước thải đô thị, nước thải công nghiệp và nông nghiệp.
  • Phương pháp phân tích hóa lý: Chỉ tiêu pH được đo trong dung dịch KCl; hàm lượng chất hữu cơ (CHC) xác định theo phương pháp Walkley-Black; dung tích trao đổi cation (CEC) định lượng bằng phương pháp Scheffer; thành phần cấp hạt cơ giới phân tích theo phương pháp ống hút kết hợp cột lắng. Hàm lượng tổng số của Cu, Pb, Zn được phá mẫu bằng hỗn hợp cường thủy ($HNO_3 : HCl = 1 : 3$) và định lượng trên máy quang phổ hấp thụ nguyên tử ngọn lửa (F-AAS). Quy trình chiết phân đoạn 5 bước của Salbu (1998) được lựa chọn để phân tách độ linh động của từng kim loại.
  • Xử lý thống kê: Mối tương quan giữa các thông số hóa lý (CHC, CEC, cấp hạt sét) và hàm lượng kim loại được phân tích thông qua hệ số tương quan Pearson ($r$) với độ tin cậy $p < 0,05$. Tiến trình nghiên cứu được triển khai xuyên suốt trong 12 tháng, kết hợp giữa khảo sát thực địa và phân tích trong phòng thí nghiệm chuyên chuẩn.

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

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

  1. Mức độ tích lũy kẽm (Zn) ở mức báo động đỏ trên toàn tuyến: Hàm lượng Zn tổng số dao động từ 277 mg/kg đến 1226 mg/kg, với giá trị trung bình toàn lưu vực đạt 648 mg/kg. Toàn bộ 7/8 điểm quan trắc đều vượt tiêu chuẩn chất lượng trầm tích Canada (315 mg/kg). Điểm cực đại ghi nhận tại Cầu Tả Thanh Oai (1226 mg/kg), vượt ngưỡng an toàn tới 389%, tiếp đến là Cầu Chiếc Hiền Giang (887 mg/kg, vượt 282%) và Cầu Trắng (763 mg/kg, vượt 242%).
  2. Ô nhiễm chì (Pb) cục bộ tại các khu vực đô thị và làng nghề: Nồng độ Pb dao động từ 65 mg/kg đến 156 mg/kg. Trong đó, 50% số điểm lấy mẫu vượt tiêu chuẩn Canada (91,3 mg/kg). Điểm ô nhiễm chì cao nhất là Cầu Diễn với 156 mg/kg (vượt tiêu chuẩn 170%), Cầu Tả Thanh Oai đạt 105 mg/kg (vượt 115%) và Cầu Chiếc đạt 100 mg/kg (vượt 110%).
  3. Hàm lượng đồng (Cu) duy trì dưới ngưỡng rủi ro: Nồng độ Cu ghi nhận trong khoảng 74,8 mg/kg đến 130,0 mg/kg, đạt đỉnh tại Cầu Tả Thanh Oai (130 mg/kg). So sánh với tiêu chuẩn Canada (197 mg/kg), toàn bộ 100% mẫu phân tích đều nằm trong giới hạn cho phép, chưa có dấu hiệu gây độc sinh thái cấp tính.
  4. Phân bố dạng liên kết hóa học: Dữ liệu chiết liên tiếp cho thấy phần lớn kim loại tồn tại ở dạng trơ trong mạng tinh thể khoáng (F5) và dạng liên kết oxit Fe-Mn (F3) hoặc hữu cơ (F4). Tại Cầu Diễn, dạng trơ F5 của Pb đạt 120,72 mg/kg và Zn đạt 431,28 mg/kg. Dạng trao đổi dễ linh động (F1) chiếm tỷ lệ nhỏ (khoảng 0,5% đến 2,5%), tuy nhiên dạng liên kết oxit Fe-Mn (F3) của Zn chiếm tới 216,66 mg/kg tại Cầu Diễn và 265,60 mg/kg tại Tả Thanh Oai, tiềm ẩn nguy cơ tái hòa tan cao.

Thảo luận kết quả

Dữ liệu nghiên cứu phản ánh quy luật phân bố không gian rõ rệt: mức độ ô nhiễm tăng vọt từ thượng nguồn (Cống Liên Mạc có Cu 75 mg/kg, Pb 92 mg/kg, Zn 323 mg/kg) đến khu vực trung lưu tiếp nhận nước thải sông Tô Lịch (Tả Thanh Oai có Zn lên tới 1226 mg/kg), sau đó giảm dần về phía hạ lưu (Cống Phủ Lý còn Zn 277 mg/kg, Pb 65 mg/kg) nhờ quá trình tự làm sạch và pha loãng thủy lực. Sự tích tụ vượt trội của Zn và Pb tại Cầu Diễn và Tả Thanh Oai bắt nguồn trực tiếp từ nguồn xả của cụm công nghiệp Phú Minh, xí nghiệp dệt len Hà Đông, cơ khí Tân Hòa và làng nghề Lụa Vạn Phúc (nơi xả khoảng 1.000 m³/ngày nước thải chứa nhiều hóa chất nhuộm và kim loại).

Mức độ tích lũy kim loại nặng tại các vị trí trọng điểm (mg/kg):
+---------------------+-------+--------+--------+
| Vị trí lấy mẫu      | Đồng  |  Chì   |  Kẽm   |
|                     | (Cu)  |  (Pb)  |  (Zn)  |
+---------------------+-------+--------+--------+
| T1: Cống Liên Mạc   | 75,0  |  92,0  |  323,0 |
| T2: Cầu Diễn        | 107,0 | 156,0  |  695,0 |
| T3: Cầu Trắng       | 101,0 |  98,0  |  763,0 |
| T4: Tả Thanh Oai    | 130,0 | 105,0  | 1226,0 |
| T5: Cầu Chiếc       | 102,0 | 100,0  |  887,0 |
| Tiêu chuẩn Canada   | 197,0 |  91,3  |  315,0 |
+---------------------+-------+--------+--------+

Phân tích tương quan Pearson chứng minh vai trò then chốt của các đặc tính hóa lý bùn đáy:

  • Hàm lượng chất hữu cơ (CHC dao động từ 2,96% đến 16,63%) có tương quan thuận cực kỳ chặt chẽ với Zn ($r = 0,95$) và Cu ($r = 0,94$).
  • Cấp hạt sét mịn ($< 0,002\text{ mm}$, chiếm từ 20% đến 35%) tương quan chặt chẽ với Zn ($r = 0,81$) và Cu ($r = 0,80$).
  • Dung tích trao đổi cation (CEC từ 15,25 đến 27,50 cmol/kg) tương quan cao với Zn ($r = 0,81$) và Cu ($r = 0,69$).

Ngược lại, chì (Pb) thể hiện mức tương quan trung bình đến yếu với CHC ($r = 0,45$), cấp hạt sét ($r = 0,35$) và CEC ($r = 0,22$). Hiện tượng này được giải thích là do Pb có khả năng tạo kết tủa khó tan dưới dạng hợp chất sulfua ($PbS$) hoặc cacbonat ($PbCO_3$) trong môi trường yếm khí khử mạnh, ít phụ thuộc thuần túy vào các vị trí hấp phụ hữu cơ hay khoáng sét. Đường đẳng nhiệt Freundlich xác nhận trầm tích tại Tả Thanh Oai và Cầu Diễn có dung lượng hấp phụ kim loại lớn nhất nhờ tỷ lệ keo sét và mùn hữu cơ cao.

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

  1. Thực hiện nạo vét cơ học bùn đáy có kiểm soát: Sở Tài nguyên và Môi trường thành phố Hà Nội cần phối hợp cùng các đơn vị chuyên môn tiến hành nạo vét khẩn cấp lớp trầm tích ô nhiễm tầng 0–30 cm tại 3 điểm nóng trọng điểm gồm Cầu Diễn, Cầu Trắng và Cầu Tả Thanh Oai. Mục tiêu nạo vét hướng tới loại bỏ 45% đến 60% khối lượng Zn và Pb đang tích tụ trong lòng dẫn, hoàn thành giai đoạn 1 trong khung thời gian 2024–2026. Quá trình hút bùn phải áp dụng kỹ thuật quây chắn kín nhằm ngăn ngừa phát tán ngược chất lơ lửng chứa độc tố về phía hạ lưu.
  2. Hiện đại hóa trạm xử lý nước thải công nghiệp và làng nghề: Ủy ban nhân dân thành phố Hà Nội cần ban hành chế tài bắt buộc các cơ sở sản xuất tại Cụm công nghiệp Phú Minh, làng lụa Vạn Phúc và xã Minh Khai đầu tư hệ thống tách dòng và xử lý kim loại nặng cục bộ. Chỉ tiêu xả thải phải bảo đảm giảm tối thiểu 80% nồng độ ion Cu, Pb, Zn trước khi hòa vào mạng lưới thoát nước chung, áp dụng quy chuẩn nghiêm ngặt từ quý 2/2025.
  3. Ứng dụng giải pháp xử lý sinh học (Phytoremediation): Triển khai trồng các dải thực vật thủy sinh ngập nước như cỏ Vetiver (Chrysopogon zizanioides) và lục bình tại các đoạn sông mở rộng, vùng đất ngập nước tự nhiên ven sông thuộc huyện Thanh Trì và Thường Tín. Mục tiêu sinh học đặt ra là hấp thu và cố định từ 20% đến 35% lượng ion kim loại di động trong nước và lớp bùn lắng, duy trì môi trường pH ổn định từ 6,8 đến 7,4 nhằm triệt tiêu rủi ro kích hoạt "quả bom hóa học", thực hiện liên tục trong giai đoạn 2025–2028.
  4. Xây dựng mạng lưới quan trắc chất lượng trầm tích liên tỉnh: Tổng cục Môi trường chủ trì liên kết Sở Tài nguyên và Môi trường Hà Nội và Hà Nam thiết lập 8 trạm giám sát định kỳ 3 tháng/lần đối với cả pha nước và pha trầm tích. Ứng dụng hệ thống cảm biến đo đạc tự động các thông số Eh, pH, DO để dự báo sớm nguy cơ giải phóng kim loại nặng khi mùa kiệt xuất hiện dòng chảy thấp.

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

  • Cơ quan quản lý nhà nước và nhà hoạch định chính sách: Cán bộ công tác tại Cục Quản lý Tài nguyên Nước, Chi cục Bảo vệ Môi trường Hà Nội và Hà Nam có thể sử dụng bộ số liệu thực địa chi tiết để xây dựng bản đồ phân vùng ô nhiễm, phục vụ trực tiếp cho Đề án bảo vệ môi trường lưu vực sông Nhuệ - sông Đáy giai đoạn mới.
  • Kỹ sư và chuyên gia công nghệ môi trường: Đội ngũ chuyên môn phụ trách thiết kế hệ thống xử lý nước thải công nghiệp, đơn vị tư vấn xử lý bùn thải đô thị có thể khai thác các thông số đẳng nhiệt hấp phụ Freundlich ($K_F, n$) để tối ưu hóa quá trình tính toán vật liệu hấp phụ bùn khoáng.
  • Giảng viên và nhà nghiên cứu học thuật: Các nhà khoa học thuộc lĩnh vực Địa hóa môi trường, Khoa học Đất và Độc học sinh thái sử dụng quy trình chiết 5 bước phân đoạn Salbu cùng ma trận tương quan Pearson trong công trình làm tài liệu tham khảo giảng dạy và phát triển đề tài nghiên cứu chuyên sâu.
  • Học viên cao học và sinh viên chuyên ngành: Người học thuộc các ngành Khoa học Môi trường, Kỹ thuật Môi trường và Quản lý Môi trường tìm thấy trong luận văn một phương pháp luận chuẩn mực về quy hoạch mạng lưới quan trắc, kỹ thuật phân tích AAS và xử lý dữ liệu thống kê môi trường.

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

Tại sao hàm lượng kẽm (Zn) trong trầm tích sông Nhuệ lại cao vượt trội so với đồng (Cu) và chì (Pb)?

Hàm lượng Zn dao động từ 323 mg/kg đến 1226 mg/kg, vượt chuẩn Canada tới 3,89 lần tại Tả Thanh Oai. Nguyên nhân do lưu vực tiếp nhận lượng lớn nước thải xi mạ, đúc kim loại, sản xuất cao su từ các cụm công nghiệp và làng nghề dọc sông. Đồng thời, Zn có ái lực rất cao với chất hữu cơ ($r = 0,95$) và cấp hạt sét ($r = 0,81$), dẫn tới xu hướng tích tụ áp đảo trong lớp bùn đáy.

Dạng liên kết hóa học nào của kim loại nặng trong trầm tích tiềm ẩn nguy cơ gây độc lớn nhất?

Dù dạng trao đổi (F1) chỉ chiếm 0,5% đến 2,5%, nhưng dạng liên kết oxit Fe-Mn (F3) và dạng hữu cơ (F4) mới là mối đe dọa tiềm tàng lớn nhất. Khi sông Nhuệ chuyển sang điều kiện yếm khí hoặc pH hạ thấp, lượng Zn liên kết Fe-Mn lên tới 265,60 mg/kg tại Cầu Tả Thanh Oai sẽ bị phân hủy và tái hòa tan cực nhanh vào nguồn nước mặt.

Tại sao chỉ phân tích chất lượng nước mặt là chưa đủ để đánh giá ô nhiễm sông Nhuệ?

Nước mặt chịu ảnh hưởng lớn bởi biến động lưu lượng dòng chảy và quá trình pha loãng theo mùa, làm nồng độ kim loại hòa tan có thể thấp giả tạo. Trong khi đó, hơn 85% kim loại nặng xâm nhập vào sông có xu hướng lắng đọng và tích tụ bền vững trong trầm tích đáy với nồng độ gấp hàng nghìn lần trong nước, phản ánh chính xác lịch sử xả thải dài hạn.

Yếu tố lý hóa nào quyết định mạnh nhất đến khả năng lưu giữ kim loại nặng tại sông Nhuệ?

Chất hữu cơ (CHC) là yếu tố chi phối mạnh mẽ nhất với hệ số tương quan $r = 0,95$ đối với Zn và $r = 0,94$ đối với Cu. Trầm tích tại Cầu Tả Thanh Oai có hàm lượng CHC cao nhất (16,63%), tạo ra mạng lưới phức chất chelat bền vững với ion kim loại, kết hợp dung tích trao đổi cation cao (27,5 cmol/kg) để giữ chặt các chất ô nhiễm.

Phương pháp phân đoạn Salbu có ưu điểm gì so với phương pháp xác định hàm lượng tổng số?

Phương pháp xác định tổng số bằng cường thủy chỉ cung cấp tổng nồng độ kim loại mà không làm rõ mức độ độc tính. Quy trình chiết liên tiếp 5 bước của Salbu bóc tách chính xác 5 dạng hóa học từ linh động đến bền vững, giúp đánh giá định lượng khả năng sinh khả dụng và nguy cơ dịch chuyển thực tế của ion kim loại vào chuỗi thức ăn thủy sinh.

Kết luận

  • Luận văn đã định lượng toàn diện hàm lượng Cu, Pb, Zn tại 8 điểm đại diện dọc 74 km dòng sông Nhuệ, ghi nhận nồng độ Zn cực đại đạt 1226 mg/kg (vượt tiêu chuẩn Canada 3,89 lần) và Pb đạt 156 mg/kg (vượt tiêu chuẩn 1,7 lần).
  • Xác lập cơ chế địa hóa kiểm soát ô nhiễm thông qua ma trận tương quan Pearson: chất hữu cơ (2,96–16,63%) và cấp hạt sét (20–35%) đóng vai trò là tác nhân lưu giữ chủ đạo với hệ số tương quan $r > 0,80$ đối với Cu và Zn.
  • Làm rõ phân bố hóa học 5 pha theo Salbu: dạng trơ trong mạng khoáng F5 và dạng liên kết oxit Fe-Mn F3 chiếm tỷ trọng lớn nhất, đóng vai trò nguồn dự trữ phát tán kim loại tiềm ẩn.
  • Khẳng định tính cấp thiết của việc kiểm soát ô nhiễm tại khu vực trung lưu đô thị (đoạn từ Cầu Diễn đến Cầu Tả Thanh Oai) – nơi chịu tải lượng nước thải sinh hoạt và công nghiệp cao nhất toàn lưu vực.
  • Đề xuất lộ trình xử lý môi trường đồng bộ từ nạo vét cơ học 45–60% bùn đáy ô nhiễm giai đoạn 2024–2026 đến triển khai quan trắc tự động và xử lý sinh học liên tỉnh giai đoạn 2025–2028.

Công trình là tài liệu học thuật có giá trị thực tiễn cao cho công tác quản trị lưu vực sông. Các nhà quản lý, kỹ sư và nghiên cứu sinh quan tâm có thể khai thác sâu các bộ thông số thực nghiệm để phát triển các mô hình kiểm soát ô nhiễm trầm tích bền vững cho các dòng sông đô thị tại Việt Nam.