Tổng quan nghiên cứu

Theo thống kê từ Quỹ Nhi đồng Liên Hiệp Quốc và các tổ chức môi trường quốc tế, mỗi ngày có khoảng 2 triệu tấn chất thải sinh hoạt xả trực tiếp ra hệ thống sông ngòi và đại dương. Tại Việt Nam, áp lực từ tốc độ công nghiệp hóa và đô thị hóa nhanh chóng đang tạo ra những thách thức nặng nề đối với chất lượng nước mặt. Tỉnh Bắc Ninh là một trong những địa phương có mật độ dân số cao hàng đầu cả nước với 1.404 người/km² và sở hữu tỷ trọng làng nghề lớn khi chiếm 18% tổng số làng nghề toàn quốc. Mặc dù các làng nghề tiểu thủ công nghiệp đóng góp 75% đến 80% giá trị sản xuất công nghiệp ngoài quốc doanh của tỉnh, hoạt động xả thải chưa qua xử lý đã làm suy thoái nghiêm trọng các nguồn nước tiếp nhận.

Vấn đề cốt lõi đặt ra là các phương pháp quan trắc truyền thống trước đây chủ yếu dựa vào việc phân tích từng chỉ tiêu riêng lẻ tại các điểm cố định rồi đối chiếu giới hạn cho phép. Cách tiếp cận này bộc lộ hạn chế lớn khi không phản ánh được bức tranh toàn diện về chất lượng nước theo không gian và diễn biến theo thời gian trên toàn thủy vực. Luận văn tập trung giải quyết bài toán này với mục tiêu ứng dụng Hệ thống Thông tin Địa lý (GIS) kết hợp Chỉ số Chất lượng Nước (WQI) để đánh giá và mô hình hóa phân bố chất lượng nước mặt hệ thống sông Cầu khu vực thành phố Bắc Ninh và các huyện phụ cận.

Phạm vi nghiên cứu bao quát mạng lưới sông Cầu và các kênh rạch liên quan qua 4 quý trong năm 2017. Công trình có ý nghĩa thực tiễn to lớn khi xây dựng cơ sở dữ liệu không gian tích hợp 32 điểm quan trắc với 9 thông số hóa lý và vi sinh. Kết quả lượng hóa trực quan giúp phát hiện kịp thời các vùng ô nhiễm nghiêm trọng có chỉ số WQI xuống mức báo động từ 5 đến 25, cung cấp công cụ đắc lực cho các nhà quản lý trong việc hoạch định chính sách bảo vệ tài nguyên nướ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 được xây dựng trên nền tảng lý thuyết phân tích không gian và địa thống kê trong GIS, kết hợp với các mô hình đánh giá chất lượng môi trường hiện đại. Trọng tâm lý thuyết bao gồm:

  • Lý thuyết nội suy không gian: Ứng dụng quy luật tương quan khoảng cách địa lý thông qua các thuật toán nội suy tất định như Nghịch đảo khoảng cách có trọng số (IDW), hàm Spline và phương pháp địa thống kê Kriging. Trong đó, Kriging khai thác cấu trúc biến thiên không gian thông qua hàm bán phương sai (Semivariogram) để tối thiểu hóa phương sai sai số ước tính.
  • Khung đánh giá Chỉ số Chất lượng Nước (WQI): Áp dụng phương pháp tính toán chuẩn hóa do Tổng cục Môi trường ban hành theo Quyết định số 879/QĐ-TCMT, kết hợp hệ thống quy chuẩn kỹ thuật quốc gia QCVN 08-MT:2015/BTNMT. Mô hình này chuyển đổi các giá trị nồng độ phức tạp thành một chỉ số tích hợp duy nhất chạy từ 0 đến 100.

Hệ thống khái niệm then chốt được vận hành trong nghiên cứu gồm:

  • Chất lượng nước mặt: Trạng thái thành phần hóa học, vật lý và sinh học của nước tồn tại trên bề mặt lục địa.
  • Hệ thống Thông tin Địa lý (GIS): Hệ thống máy tính thu thập, lưu trữ, xử lý, phân tích và trực quan hóa dữ liệu không gian và thuộc tính.
  • Chỉ số WQI: Thang đo không thứ nguyên thể hiện mức độ trong sạch hay ô nhiễm tổng thể của nguồn nước.
  • Sai số trung phương (RMSE): Thước đo định lượng dùng để đánh giá độ chính xác và kiểm định sự phù hợp của mô hình nội suy.

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

Nguồn dữ liệu đầu vào của luận văn được thu thập đồng bộ từ hai nguồn chính:

  • Dữ liệu thuộc tính: Kết quả quan trắc môi trường nước mặt tại 32 vị trí đại diện trên toàn thủy vực sông Cầu và các vùng tiếp nhận nước thải làng nghề, thực hiện đều đặn qua 4 quý của năm 2017 bởi các cơ quan chuyên môn tại Bắc Ninh.
  • Dữ liệu không gian: Bản đồ hiện trạng sử dụng đất tỷ lệ 1:10.000 trong hệ quy chiếu quốc gia VN-2000, được chiết tách lớp mặt nước để làm khung ranh giới nội suy.

Cỡ mẫu và phương pháp chọn mẫu: Cỡ mẫu nghiên cứu gồm 32 điểm quan trắc được bố trí theo phương pháp chọn mẫu phân tầng kết hợp chọn mẫu có chủ đích. Các điểm đo phân bố phủ khắp dòng chính sông Cầu, các nhánh sông liên huyện và các kênh xả thải chính từ cụm công nghiệp. Toàn bộ 32 điểm được phân chia khoa học thành hai tập dữ liệu độc lập: 26 điểm được sử dụng làm mẫu huấn luyện để chạy thuật toán nội suy không gian và 6 điểm độc lập được giữ lại làm mẫu kiểm định chéo độ chính xác mô hình.

Phương pháp phân tích và lý do lựa chọn: Quy trình phân tích gồm 9 thông số đại diện: pH, DO, TSS, COD, BOD5, N-NH4+, P-PO4(3-), độ đục và tổng Coliform. Phương pháp tính toán WQI theo Quyết định 879/QĐ-TCMT được lựa chọn vì tính pháp lý cao và phù hợp với đặc thù nguồn nước tiếp nhận tại Việt Nam. Về xử lý không gian, phương pháp Kriging được lựa chọn ưu tiên so với IDW sau khi thực hiện phân tích đối sánh sai số trung phương. Kriging khắc phục được nhược điểm gán trọng số khoảng cách đơn thuần của IDW, giúp tạo ra bề mặt biến thiên liên tục, mịn và có độ tin cậy thống kê cao nhất dọc theo dòng chảy sông ngòi. Timeline nghiên cứu bao quát trọn vẹn 12 tháng của năm 2017, phản ánh rõ nét sự thay đổi chất lượng nước theo mùa mưa và mùa khô.

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

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

Nghiên cứu đã lượng hóa chi tiết hiện trạng ô nhiễm nguồn nước mặt khu vực thành phố Bắc Ninh và các huyện phụ cận với 4 phát hiện quan trọng:

  • Chỉ số chất lượng nước WQI tổng thể ở mức báo động đỏ: Qua cả 4 quý năm 2017, chỉ số WQI tại hầu hết các trạm quan trắc chỉ dao động trong khoảng từ 5 đến 25. Đây là mức chất lượng nước rất xấu, nguồn nước bị ô nhiễm nặng nề và đòi hỏi phải có biện pháp xử lý khẩn cấp trong tương lai.
  • Ô nhiễm hữu cơ và dinh dưỡng vượt quy chuẩn nhiều lần: Nồng độ BOD5 và COD tại các điểm tiếp nhận nước thải vượt ngưỡng quy chuẩn QCVN 08-MT:2015/BTNMT cột A1 từ 5 đến 20 lần, thậm chí có thời điểm lấy mẫu ghi nhận mức vượt gần 50 lần. Hàm lượng Amoni (N-NH4+) có diễn biến đặc biệt nghiêm trọng, điển hình như tại khu vực ao thôn Đức Lý (xã Tam Đa, huyện Yên Phong) ghi nhận nồng độ Amoni vượt quy chuẩn cho phép tới 69 lần.
  • Ô nhiễm vi sinh và chất rắn lơ lửng tập trung cao tại vùng làng nghề: Các điểm xả thải tại làng nghề sản xuất giấy Phong Khê, bún Khắc Niệm và tái chế kim loại Văn Môn có hàm lượng TSS và tổng Coliform vượt giới hạn cho phép từ 6,4 đến 15 lần, cho thấy tải lượng vi khuẩn gây bệnh và cặn hữu cơ không qua lắng lọc là rất lớn.
  • Phân hóa không gian và mùa vụ rõ rệt: Chất lượng nước suy giảm mạnh nhất vào quý I và quý IV (mùa khô). Lưu lượng dòng chảy tự nhiên giảm từ 30% đến 40% trong giai đoạn này khiến khả năng tự làm sạch của sông Cầu suy giảm nghiêm trọng, làm nồng độ chất ô nhiễm tăng vọt so với mùa mưa (quý II và quý III).

Thảo luận kết quả

Nguyên nhân chính dẫn đến tình trạng ô nhiễm nghiêm trọng bắt nguồn từ việc phát triển kinh tế thiếu kiểm soát môi trường. Với mật độ dân cư đông đúc 1.404 người/km², lưu lượng nước thải sinh hoạt chưa qua xử lý xả thẳng vào các kênh mương là rất lớn. Đồng thời, các làng nghề đóng góp tới hơn 30% giá trị sản xuất công nghiệp toàn tỉnh nhưng phần lớn chưa xây dựng trạm xử lý nước thải tập trung. Nước thải từ quá trình ngâm tẩy giấy, cô đúc kim loại và chế biến nông sản chứa nhiều hóa chất độc hại được xả trực tiếp ra dòng chảy.

So sánh với các nghiên cứu tương tự trên lưu vực sông Sài Gòn - Đồng Nai và lưu vực sông Ogun - Osun tại Nigeria, có thể thấy tính quy luật tương đồng: sự kết hợp giữa nước thải sinh hoạt đô thị và chất thải công nghiệp luôn khiến các chỉ tiêu BOD5, Amoni và Coliform tăng đột biến tại các vùng hạ lưu và điểm giao cắt kênh xả.

Dữ liệu nghiên cứu được trực quan hóa sinh động thông qua hệ thống bản đồ chuyên đề GIS đa tầng và biểu đồ thống kê. Bản đồ phân vùng chất lượng nước phân tách rõ 5 dải màu theo tiêu chuẩn WQI, trong đó gam màu đỏ và da cam bao trùm phần lớn diện tích sông Cầu đoạn qua thành phố Bắc Ninh. Đồ thị sai số trung phương (RMSE) so sánh giữa Kriging và IDW qua 4 quý chứng minh rõ ràng đường cong sai số của Kriging luôn thấp hơn IDW, khẳng định Kriging là mô hình toán tối ưu để mô phỏng sự lan truyền ô nhiễm trên mặt nước.

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

Nhằm kiểm soát và từng bước phục hồi chất lượng nước mặt trên địa bàn tỉnh Bắc Ninh, nghiên cứu đề xuất 4 nhóm giải pháp trọng tâm:

  • Xây dựng và hoàn thiện trạm xử lý nước thải tập trung: Sở Tài nguyên và Môi trường tỉnh Bắc Ninh cùng Ban Quản lý các Khu công nghiệp cần đôn đốc và đưa vào vận hành 100% hệ thống xử lý nước thải đạt tiêu chuẩn cột A theo QCVN 40:2011/BTNMT tại các làng nghề trọng điểm như Phong Khê, Văn Môn và Khắc Niệm trước quý IV năm 2025, hướng tới mục tiêu nâng chỉ số WQI trung bình toàn vùng từ dưới 25 lên trên 50 (mức vàng).
  • Thiết lập mạng lưới trạm quan trắc tự động kết hợp WebGIS: Ủy ban nhân dân tỉnh Bắc Ninh chỉ đạo Trung tâm Quan trắc Môi trường triển khai lắp đặt thêm 15 trạm quan trắc tự động liên tục trên sông Cầu và các kênh tiêu chính trong giai đoạn 2024 - 2026. Dữ liệu cần được truyền trực tiếp về hệ thống WebGIS thời gian thực 24/7 với sai số dưới 5% nhằm phát hiện kịp thời các hành vi xả thải lén lút.
  • Tăng cường kiểm tra, cưỡng chế và di dời cơ sở gây ô nhiễm nặng: Thanh tra Sở Tài nguyên và Môi trường phối hợp cùng Công an Môi trường tổ chức thanh tra định kỳ 100% doanh nghiệp và cơ sở sản xuất; xử phạt nghiêm các vi phạm xả thải vượt tiêu chuẩn từ 2 lần trở lên, đồng thời hoàn thành kế hoạch di dời 80% cơ sở tái chế nhôm và giấy gây ô nhiễm nặng ra khỏi khu dân cư trước năm 2027.
  • Ứng dụng giải pháp sinh học và kiểm soát nguồn thải nông nghiệp: Sở Nông nghiệp và Phát triển nông thôn chủ trì triển khai mô hình hồ sinh học, bãi lọc ngập nước nhân tạo và tổ chức nạo vét định kỳ trên 45 km kênh mương thoát nước chính; phấn đấu giảm thiểu 35% lượng nitơ và photpho phát tán từ chất thải chăn nuôi và phân bón hóa học trong giai đoạn 2024 - 2028.

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

Nội dung và phương pháp luận của luận văn mang lại giá trị thiết thực cho 4 nhóm đối tượng sau:

  • Cơ quan quản lý nhà nước về tài nguyên và môi trường: Cán bộ tại Sở Tài nguyên và Môi trường, Chi cục Bảo vệ Môi trường và Ủy ban nhân dân các cấp có thể sử dụng quy trình tích hợp GIS - WQI để lập quy hoạch phân vùng chức năng nguồn nước, xây dựng bản đồ cảnh báo rủi ro môi trường và quản lý cấp phép xả thải cho hơn 30 khu và cụm công nghiệp trên địa bàn.
  • Các nhà khoa học, giảng viên và học viên cao học: Giảng viên và học viên thuộc các chuyên ngành Kỹ thuật Trắc địa - Bản đồ, Khoa học Môi trường, Quản lý Đất đai và Địa tin học có thể khai thác mô hình nội suy địa thống kê Kriging trên 32 điểm đo làm tài liệu tham khảo chuẩn mực để phát triển các mô hình lan truyền ô nhiễm và nghiên cứu chuyên sâu.
  • Chủ đầu tư hạ tầng khu công nghiệp và doanh nghiệp xử lý nước thải: Các kỹ sư môi trường và ban quản lý dự án có thể dựa vào số liệu ô nhiễm hữu cơ và kim loại nặng trong luận văn để tính toán chính xác công suất thiết kế trạm xử lý nước thải tập trung, tối ưu hóa vị trí cửa xả nhằm giảm tải áp lực cho nguồn nước tiếp nhận.
  • Các tổ chức phi chính phủ và cộng đồng dân cư: Các tổ chức xã hội bảo vệ môi trường và người dân lưu vực sông Cầu có nguồn dữ liệu khoa học tin cậy để nâng cao nhận thức cộng đồng, chủ động các giải pháp bảo vệ sức khỏe và tham gia giám sát việc chấp hành pháp luật môi trường của các cơ sở sản xuất.

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

  1. Tại sao nghiên cứu lại lựa chọn thuật toán Kriging thay vì IDW để nội suy chất lượng nước? Kriging là phương pháp địa thống kê tối ưu dựa trên hàm bán phương sai, cho phép lượng hóa cấu trúc tương quan không gian và giảm thiểu phương sai sai số ước tính. Kết quả kiểm tra chéo trên 6 điểm kiểm định độc lập cho thấy sai số trung phương của Kriging nhỏ hơn rõ rệt so với IDW, phản ánh liên tục và chính xác hơn hiện trạng lan truyền ô nhiễm trên 32 trạm quan trắc.

  2. Thang điểm WQI từ 5 đến 25 ghi nhận trong nghiên cứu phản ánh mức độ nguy hại ra sao? Theo Quyết định số 879/QĐ-TCMT của Tổng cục Môi trường, thang điểm WQI từ 0 đến 25 được gán màu đỏ, đại diện cho chất lượng nước ô nhiễm rất nặng. Nguồn nước ở mức này hoàn toàn không thể sử dụng cho sinh hoạt, nuôi trồng thủy sản hay tưới tiêu do nồng độ các chất hữu cơ BOD5, COD và Amoni vượt ngưỡng an toàn từ 5 đến 69 lần.

  3. Những nguồn thải chính nào làm suy thoái chất lượng nước mặt tại Bắc Ninh? Nguồn thải chính bắt nguồn từ nước thải sản xuất chưa qua xử lý của hơn 60 làng nghề truyền thống và các cụm công nghiệp, điển hình là làng giấy Phong Khê và tái chế nhôm Văn Môn. Bên cạnh đó, mật độ dân cư cao 1.404 người/km² xả lượng lớn nước thải sinh hoạt chứa Coliform và TSS trực tiếp ra sông ngòi cũng là tác nhân gây ô nhiễm nặng.

  4. Dữ liệu bản đồ không gian trong nghiên cứu được thu thập và xử lý từ nguồn nào? Dữ liệu không gian lớp thủy hệ được chiết tách trực tiếp từ bản đồ hiện trạng sử dụng đất tỷ lệ 1:10.000 trên hệ tọa độ quốc gia VN-2000 do cơ quan chuyên trách cung cấp. Nền địa lý này kết hợp với tọa độ thực địa của 32 điểm đo giúp xây dựng mô hình không gian chuẩn xác trên các phần mềm ArcGIS và Microstation V8i.

  5. Công thức tính WQI theo Quyết định 879/QĐ-TCMT tích hợp những nhóm thông số nào? Mô hình WQI của Tổng cục Môi trường tích hợp 9 thông số chia làm 4 nhóm: nhóm thông số hữu cơ và dinh dưỡng (DO, BOD5, COD, NH4+, PO4(3-)), nhóm thông số vật lý (TSS, độ đục), nhóm vi sinh (tổng Coliform) và nhóm chỉ số pH. Các chỉ số phụ được chuẩn hóa và tính toán gộp thành một thang điểm tổng hợp từ 0 đến 100.

Kết luận

  • Luận văn đã ứng dụng thành công công nghệ GIS và thuật toán địa thống kê Kriging để mô hình hóa toàn diện bức tranh phân bố chất lượng nước mặt sông Cầu qua 4 quý năm 2017.
  • Xác định hiện trạng chất lượng nước mặt khu vực thành phố Bắc Ninh và các huyện lân cận đang ở mức ô nhiễm nghiêm trọng với chỉ số WQI dao động rất thấp từ 5 đến 25.
  • Chỉ rõ các tâm điểm ô nhiễm hữu cơ, dinh dưỡng và vi sinh tại các làng nghề sản xuất giấy, bún và đúc kim loại với nồng độ chất ô nhiễm vượt quy chuẩn từ 5 đến 69 lần.
  • Chứng minh phương pháp nội suy Kriging có độ chính xác và độ tin cậy vượt trội so với IDW trong việc mô phỏng sự phân bố không gian của các thông số môi trường nước.
  • Cung cấp luận cứ khoa học và giải pháp kỹ thuật có tính ứng dụng cao phục vụ công tác quy hoạch, quản lý và kiểm soát ô nhiễm lưu vực sông.

Trong lộ trình giai đoạn 2024 - 2030, việc đầu tư hạ tầng xử lý nước thải và hiện đại hóa mạng lưới trạm quan trắc tự động là yêu cầu cấp thiết. Các nhà quản lý, nhà nghiên cứu và doanh nghiệp hãy chủ động ứng dụng dữ liệu và phương pháp luận từ luận văn để triển khai ngay những hành động thiết thực nhằm bảo vệ và phục hồi tài nguyên nước mặt quý giá của địa phương.