Tổng quan nghiên cứu

Sông Nhuệ có tổng chiều dài 74 km, trong đó đoạn chảy qua địa phận thành phố Hà Nội dài 61,5 km với diện tích lưu vực tự nhiên khoảng 1.075 km2, giữ vai trò huyết mạch trong việc cấp nước tưới tiêu nông nghiệp và tiêu thoát nước cho hơn 6 triệu dân vùng hữu ngạn sông Hồng. Tuy nhiên, tốc độ đô thị hóa nhanh chóng cùng áp lực phát triển kinh tế đã khiến dòng sông phải tiếp nhận lưu lượng nước thải khổng lồ trung bình khoảng 5,4 m3/s từ 99.927 cơ sở sản xuất kinh doanh cùng các khu dân cư tập trung chưa qua xử lý đồng bộ. Dữ liệu quan trắc môi trường cho thấy chất lượng nước suy thoái nghiêm trọng: nồng độ oxy hòa tan (DO) xuống rất thấp, dao động trong khoảng từ 2,6 đến 3,1 mg/l, trong khi nồng độ COD và BOD5 vượt quy chuẩn kỹ thuật quốc gia QCVN 08:2008/BTNMT loại B1 lần lượt tới 37,1 lần và 54,1 lần. Chỉ số chất lượng nước WQI tại nhiều đoạn sông giảm sâu xuống mức 25 đến 50, đe dọa trực tiếp đến an ninh nguồn nước và hệ sinh thái thủy vực.

Mục tiêu then chốt của luận văn là ứng dụng bộ công cụ mô hình toán MIKE 11 để mô phỏng động lực dòng chảy, đánh giá khả năng lan truyền ô nhiễm và dự báo diễn biến các thông số chất lượng nước cốt lõi gồm DO, BOD5 và tổng Nitơ (T-N) trên 6 phân đoạn sông Nhuệ đến năm 2015 và 2020. Nghiên cứu thực hiện trong phạm vi không gian từ cống Liên Mạc đến cống Phủ Lý, sử dụng chuỗi số liệu thủy văn và nguồn thải năm 2008 làm nền tảng kiểm định, kết hợp định hướng Quy hoạch tổng thể phát triển kinh tế - xã hội Thủ đô Hà Nội. Kết quả đạt được cung cấp cơ sở khoa học định lượng với mức hiệu quả mô phỏng từ 88% đến 97%, hỗ trợ đắc lực cho các nhà quản lý trong việc kiểm soát tải lượng xả thải và bảo vệ môi trường nước mặt.

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 thủy động lực học dòng chảy một chiều không ổn định trong lòng dẫn hở thông qua hệ phương trình Saint-Venant, bao gồm phương trình liên tục và phương trình động lượng. Hệ phương trình này mô tả chính xác sự biến thiên của lưu lượng dòng chảy Q và mực nước h theo thời gian t và chiều dài lòng dẫn x dưới tác động của trọng trường, độ dốc đáy và lực ma sát thành lòng dẫn. Bên cạnh đó, quá trình lan truyền và biến đổi chất hòa tan trong nước tuân theo phương trình truyền tải - khuếch tán một chiều và định luật Fick, mô tả hai cơ chế chuyển vận chính là chuyển động đối lưu do vận tốc dòng chảy và khuếch tán rối do gradient chênh lệch nồng độ chất ô nhiễm.

Hệ thống mô hình toán tích hợp MIKE 11 do Viện Thủy lực Đan Mạch (DHI) phát triển được lựa chọn làm công cụ tính toán chủ đạo. Mô hình vận hành đồng bộ 3 mô đun chuyên sâu: mô đun thủy động lực (HD) giải hệ phương trình vi phân bằng sơ đồ sai phân ẩn 6 điểm Abbott-Ionescu; mô đun truyền tải - khuếch tán (AD) mô phỏng động thái dịch chuyển của các chất ô nhiễm; và mô đun sinh thái (Ecolab) tính toán tới 13 thông số chất lượng nước ở 6 cấp độ sinh thái. Khung lý thuyết này gắn liền với các khái niệm cốt lõi: Oxy hòa tan (DO), Nhu cầu oxy sinh hóa (BOD5), Tổng Nitơ (T-N), Chỉ số chất lượng nước (WQI) và Cơ chế tự làm sạch của thủy vực tiếp nhận.

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

Nguồn dữ liệu đầu vào của mô hình được thu thập toàn diện với cỡ mẫu gồm 43 mặt cắt ngang chi tiết dọc theo 74 km dòng chính sông Nhuệ, kết hợp 60 mặt cắt mở rộng trên hệ thống sông Đáy liên kết. Dữ liệu biên thủy lực được thiết lập từ chuỗi số liệu quan trắc mực nước và lưu lượng tại 12 trạm thủy văn cùng 7 trạm khí tượng trên lưu vực. Nguồn thải được thống kê chi tiết theo từng đoạn sông từ tài liệu hiện trạng của 8 quận, huyện thuộc thành phố Hà Nội.

Phương pháp chọn mẫu phân tầng theo không gian địa lý được áp dụng để chia sông Nhuệ thành 6 phân đoạn đặc thù dựa trên điều kiện thủy lực, vị trí tiếp nhận nước thải và công năng sử dụng nước. Nhu cầu sử dụng nước và lưu lượng nước thải phát sinh được ước tính dựa trên tiêu chuẩn TCXDVN 33:2006 với định mức cấp nước sinh hoạt từ 120 đến 180 lít/người/ngày đêm và tỷ lệ phát thải bằng 80% lượng nước cấp. Phương pháp mô hình hóa số trị MIKE 11 được lựa chọn nhờ khả năng mô phỏng vượt trội đối với mạng lưới sông ngòi phân nhánh phức tạp, đảm bảo độ ổn định tính toán và độ chính xác cao. Toàn bộ chuỗi số liệu năm 2008 được sử dụng để hiệu chỉnh, kiểm định mô hình trước khi tiến hành chạy dự báo chất lượng nước cho các mốc thời gian năm 2015 và năm 2020.

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

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

Thứ nhất, mô hình MIKE 11 thể hiện mức độ tương thích vượt trội trong việc tái hiện chế độ thủy lực và chất lượng nước sông Nhuệ. Kết quả kiểm định mô hình thủy lực với dữ liệu thực đo cho thấy sai số tuyệt đối về mực nước chỉ dao động từ 0,15 m đến 0,35 m. Mức hiệu quả mô phỏng chất lượng nước đạt độ tin cậy từ 88% đến 97% đối với thông số oxy hòa tan DO và nhu cầu oxy sinh hóa BOD5 trong đợt quan trắc tháng 3/2008.

Thứ hai, nghiên cứu phát hiện sự phân hóa ô nhiễm nghiêm trọng giữa các phân đoạn sông. Đoạn 2 (từ Cầu Hà Đông đến Cầu Tó dài 4,34 km) và Đoạn 3 (từ Cầu Tó đến Cầu Chiếc dài 17,14 km) là hai khu vực ô nhiễm nặng nề nhất. Nồng độ BOD5 tại đây luôn vượt tiêu chuẩn quy định từ 4 đến 6 lần do tiếp nhận trực tiếp toàn bộ nước thải sinh hoạt chưa qua xử lý từ sông Tô Lịch và các đô thị trung tâm.

Thứ ba, kịch bản dự báo đến năm 2015 và 2020 chỉ ra xu hướng suy thoái nước mặt gia tăng mạnh mẽ. Do tốc độ tăng dân số cơ học đạt 1% đến 2%/năm (dân số dự kiến đạt 9,1 triệu người vào năm 2030), lưu lượng nước thải sinh hoạt tăng từ 15% đến 25%. Nồng độ DO tại trung lưu sông Nhuệ tiếp tục giảm sâu dưới ngưỡng 2,0 mg/l, trong khi hàm lượng T-N và BOD5 dự báo tăng từ 20% đến 35% so với năm 2008 nếu không có biện pháp xử lý bổ sung.

Thảo luận kết quả

Nguyên nhân cốt lõi dẫn đến tình trạng ô nhiễm nghiêm trọng trên sông Nhuệ là do sự mất cân đối gay gắt giữa lưu lượng xả thải khổng lồ khoảng 5,4 m3/s và lưu lượng dòng chảy tự nhiên trong mùa kiệt. Khả năng tự làm sạch của dòng sông gần như bị tê liệt hoàn toàn khi hàm lượng chất hữu cơ vượt quá ngưỡng tự phân hủy sinh học. Các kết quả mô phỏng của luận văn hoàn toàn tương đồng với nghiên cứu của Nguyễn Thanh Sơn năm 2011 (nồng độ DO đạt 2,89 mg/l, BOD5 vượt chuẩn loại A từ 4 đến 6 lần) và nghiên cứu của Nguyễn Thị Nga năm 2012 (DO dao động từ 2,6 đến 3,1 mg/l). Độ chính xác mô phỏng của mô hình (88% đến 97%) cũng tương đương với các nghiên cứu ứng dụng MIKE 11 trên lưu vực sông Cầu (96,4% đến 99,1%) và sông Sài Gòn - Đồng Nai (90% đến 94%) của nhóm tác giả Trần Hồng Thái.

Trong cấu trúc báo cáo, diễn biến ô nhiễm được trực quan hóa rõ nét qua các biểu đồ phân bố nồng độ chất theo chiều dài dòng sông và bảng thống kê sai số tại các trạm then chốt như Cống Liên Mạc, Cầu Hà Đông, Cầu Tó, Cầu Chiếc, Đồng Quan và Cống Thần, giúp phân định rõ ranh giới các điểm nóng ô nhiễm phục vụ công tác điều hành thực địa.

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

Thứ nhất, đầu tư xây dựng và hoàn thiện hệ thống thu gom, xử lý nước thải sinh hoạt đô thị tập trung. Ủy ban nhân dân thành phố Hà Nội phối hợp với Sở Xây dựng khẩn trương đẩy nhanh tiến độ các nhà máy xử lý nước thải quy mô lớn dọc lưu vực, hướng tới mục tiêu thu gom và xử lý đạt quy chuẩn từ 70% đến 80% tổng lượng nước thải đô thị (tương đương khoảng 4,0 m3/s đến 4,5 m3/s), kiểm soát nồng độ BOD5 đầu ra dưới 30 mg/l trong giai đoạn 2015 đến 2020.

Thứ hai, điều tiết thủy lực chủ động và bổ cập nguồn nước sạch thau rửa dòng chảy. Công ty TNHH MTV Đầu tư phát triển thủy lợi Sông Nhuệ cần phối hợp vận hành tối ưu cống Liên Mạc để tiếp dẫn nguồn nước từ sông Hồng với lưu lượng duy trì từ 25 m3/s đến 30 m3/s trong suốt mùa kiệt. Giải pháp này giúp tăng cường khả năng pha loãng, nâng nồng độ DO toàn tuyến lên trên 4,0 mg/l, đáp ứng tiêu chuẩn QCVN 08:2008/BTNMT cột B1 phục vụ tưới tiêu nông nghiệp.

Thứ ba, thanh tra và bắt buộc kiểm soát triệt để các nguồn thải công nghiệp, làng nghề. Sở Tài nguyên và Môi trường thành phố Hà Nội cùng Cảnh sát Môi trường triển khai kiểm tra liên ngành, yêu cầu 100% các khu, cụm công nghiệp và bệnh viện tại 8 quận, huyện ven sông phải vận hành hệ thống xử lý nước thải đạt chuẩn trước khi xả thải, phấn đấu giảm 90% lượng tải lượng ô nhiễm hữu cơ và hóa chất độc hại từ 99.927 cơ sở sản xuất kinh doanh trước năm 2020.

Thứ tư, thiết lập hệ thống quan trắc tự động và ứng dụng mô hình MIKE 11 trong cảnh báo sớm. Chi cục Bảo vệ Môi trường Hà Nội cần tích hợp mô hình số MIKE 11 với hệ thống bản đồ số GIS và mạng lưới trạm quan trắc tự động liên tục, rút ngắn thời gian cảnh báo sự cố ô nhiễm nước xuống dưới 2 giờ và nâng cao độ chính xác dự báo lan truyền chất ô nhiễm đạt trên 90% từ năm 2020.

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

Nhóm cơ quan quản lý nhà nước về tài nguyên và môi trường (Bộ Tài nguyên và Môi trường, Sở Tài nguyên và Môi trường Hà Nội và Hà Nam): Sử dụng cơ sở dữ liệu định lượng và kết quả phân vùng tải lượng ô nhiễm để hoạch định chính sách cấp phép xả thải, đánh giá sức chịu tải của sông Nhuệ và ban hành các quy định bảo vệ môi trường lưu vực sông liên tỉnh.

Nhóm đơn vị khai thác công trình thủy lợi và nông nghiệp (Công ty Khai thác công trình thủy lợi Sông Nhuệ): Ứng dụng các kịch bản thủy lực từ 43 mặt cắt thực đo để tối ưu hóa quy trình đóng mở cống Liên Mạc, cống Lương Cổ, chủ động điều tiết nguồn nước tưới an toàn cho hơn 50.000 ha đất canh tác nông nghiệp, hạn chế tối đa rủi ro đưa nước ô nhiễm vào đồng ruộng.

Nhóm nghiên cứu, giảng viên và học viên cao học chuyên ngành Môi trường, Thủy văn: Tham khảo khung phương pháp luận kết hợp 3 mô đun HD, AD và Ecolab của MIKE 11 với độ tin cậy từ 88% đến 97%. Đây là tài liệu học thuật mẫu mực cho các đề tài nghiên cứu mô hình hóa chất lượng nước mặt và tính toán lan truyền chất ô nhiễm trên mạng lưới sông ngòi phức tạp.

Nhóm doanh nghiệp, ban quản lý khu công nghiệp và đơn vị tư vấn môi trường: Nắm rõ hiện trạng và giới hạn tiếp nhận nồng độ BOD5, T-N trên 6 phân đoạn sông, từ đó thiết kế quy chuẩn công nghệ trạm xử lý nước thải cục bộ đạt tiêu chuẩn xả thải loại A hoặc B theo quy định pháp luật.

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

Ưu điểm nổi bật của mô hình MIKE 11 trong việc đánh giá chất lượng nước sông Nhuệ là gì? Mô hình MIKE 11 tích hợp đồng bộ 3 mô đun HD, AD và Ecolab, sử dụng sơ đồ sai phân ẩn 6 điểm Abbott-Ionescu giúp mô phỏng chính xác chế độ thủy lực trên mạng sông phân nhánh phức tạp dài 74 km với 43 mặt cắt. Độ chính xác kiểm định đạt từ 88% đến 97%, vượt trội so với các mô hình một chiều truyền thống.

Phân đoạn nào trên sông Nhuệ bị ô nhiễm nghiêm trọng nhất và nguyên nhân do đâu? Đoạn sông từ Cầu Hà Đông đến Cầu Chiếc dài hơn 21 km (thuộc đoạn 2 và đoạn 3) ô nhiễm nặng nhất. Nguyên nhân do tiếp nhận trực tiếp nước thải sinh hoạt chưa qua xử lý từ sông Tô Lịch và khu vực nội thành với lưu lượng khoảng 5,4 m3/s, khiến DO giảm sâu dưới 2,0 mg/l và BOD5 vượt chuẩn hơn 50 lần.

Nước sông Nhuệ hiện tại có đủ điều kiện phục vụ tưới tiêu nông nghiệp hay không? Nước sông Nhuệ không đáp ứng tiêu chuẩn QCVN 08:2008/BTNMT cột B1 dùng cho tưới tiêu. Các thông số COD và BOD5 vượt quy chuẩn từ 37,1 đến 54,1 lần, hàm lượng hợp chất Nitơ (Amoni, Nitrat) rất cao, tiềm ẩn nguy cơ tích tụ độc chất trong đất và ảnh hưởng nghiêm trọng đến cây trồng.

Dữ liệu đầu vào phục vụ dự báo chất lượng nước đến năm 2015 và 2020 dựa trên cơ sở nào? Dự báo dựa trên Quyết định số 1878/QĐ-TTg phê duyệt Quy hoạch chung xây dựng Thủ đô Hà Nội đến năm 2030, kết hợp định mức tiêu chuẩn cấp nước TCXDVN 33:2006 (120 đến 180 lít/người/ngày) và số liệu quy hoạch phát triển kinh tế - xã hội của 8 quận, huyện ven sông để tính toán tải lượng xả thải.

Biện pháp kỹ thuật nào giúp cải thiện nhanh chóng chất lượng nước sông Nhuệ trong mùa khô? Biện pháp cấp bách là điều tiết cống Liên Mạc để lấy nước từ sông Hồng bổ cập vào sông Nhuệ với lưu lượng từ 25 m3/s đến 30 m3/s trong mùa kiệt nhằm tăng cường pha loãng chất ô nhiễm, kết hợp nạo vét bùn lắng và vận hành các trạm xử lý nước thải cục bộ.

Kết luận

  • Luận văn đã ứng dụng thành công hệ thống mô hình MIKE 11 (tích hợp các mô đun HD, AD, Ecolab) để mô phỏng động lực dòng chảy và diễn biến chất lượng nước trên 74 km sông Nhuệ với độ tin cậy đạt từ 88% đến 97%.
  • Nghiên cứu đã phân đoạn chi tiết 6 khúc sông và định lượng chính xác áp lực từ các nguồn thải với lưu lượng khoảng 5,4 m3/s, xác định rõ đoạn Cầu Hà Đông - Cầu Chiếc là điểm nóng ô nhiễm hữu cơ vượt quy chuẩn trên 50 lần.
  • Thiết lập thành công kịch bản dự báo các thông số DO, BOD5 và T-N đến năm 2015 và năm 2020 dựa trên dữ liệu quy hoạch phát triển kinh tế - xã hội và tỷ lệ phát sinh nước thải sinh hoạt chiếm 80% nước cấp.
  • Đóng góp hệ thống giải pháp đồng bộ kết hợp giữa công trình (xây dựng nhà máy xử lý nước thải, bổ cập 25 đến 30 m3/s nước sông Hồng) và phi công trình (quan trắc tự động, kiểm soát nguồn thải công nghiệp) nhằm phục hồi chất lượng nước trước năm 2025.
  • Đề xuất các cơ quan chức năng khẩn trương số hóa công tác quản lý nguồn nước trong giai đoạn 2016 đến 2020 và khuyến nghị các nhà nghiên cứu tiếp tục mở rộng mô hình thủy động lực 2 chiều tại các điểm hợp lưu sông Đáy - Nhuệ.