Tổng quan nghiên cứu

Lưu vực sông Cầu là một trong những lưu vực sông trọng điểm tại miền Bắc Việt Nam, có diện tích lưu vực tự nhiên khoảng 6.030 km² trải rộng trên địa bàn 6 tỉnh, thành phố bao gồm Bắc Kạn, Thái Nguyên, Bắc Giang, Bắc Ninh, Vĩnh Phúc và Hà Nội. Đây là không gian sinh sống của hơn 8,5 triệu người với mật độ dân số bình quân 557 người/km², đồng thời là trung tâm công nghiệp, nông nghiệp và khai khoáng phát triển năng động bậc nhất vùng kinh tế trọng điểm Bắc Bộ. Nguồn nước mặt trên dòng chính sông Cầu dài 288,5 km giữ vai trò huyết mạch trong việc cung cấp nước cho sinh hoạt, công nghiệp, thủy sản và tưới tiêu nông nghiệp.

Tuy nhiên, tài nguyên nước mặt trên lưu vực đang phải đối mặt với áp lực kép về số lượng và chất lượng. Sự biến động sâu sắc của khí hậu gió mùa khiến phân bố dòng chảy không đều theo thời gian: khoảng 76% đến 80% tổng lượng nước tập trung vào 5 tháng mùa mưa (từ tháng V đến tháng IX), trong khi 8 tháng mùa kiệt chỉ đóng góp 20% đến 24% tổng lượng dòng chảy năm, dẫn đến tình trạng khan hiếm nước nghiêm trọng vào các tháng I, II và III. Song song với đó, sự phát triển của hơn 2.000 doanh nghiệp sản xuất công nghiệp, 27 khu công nghiệp tập trung cùng hàng trăm cơ sở khai thác khoáng sản và làng nghề đã xả thải lượng lớn chất ô nhiễm hữu cơ, chất rắn lơ lửng ra sông mà chưa qua xử lý triệt để, làm suy thoái chất lượng nước vùng trung và hạ lưu.

Luận văn thạc sĩ chuyên ngành Thủy văn học của tác giả Mai Thị Ngân Anh hướng tới mục tiêu đánh giá toàn diện hiện trạng tài nguyên nước mặt cả về trữ lượng và chất lượng; lượng hóa mâu thuẫn sử dụng nước thông qua mô phỏng cân bằng nước; phân tích lan truyền chất ô nhiễm dọc dòng chính và đề xuất giải pháp quản lý tổng hợp tài nguyên nước bền vững cho giai đoạn 2010 - 2020 và định hướng tương lai.

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 cân bằng nước lưu vực (Water Balance Theory) kết hợp nguyên lý bảo toàn khối lượng và chuyển động dòng chảy một chiều không ổn định trong mạng lưới sông ngòi tự nhiên. Về mặt phân tích chuỗi dữ liệu, đề tài ứng dụng lý thuyết xác suất thống kê thủy văn theo phân bố Pearson loại III và lý thuyết Trêbưsep để xác định tính đồng pha, chu kỳ dao động nhiều năm của mưa và dòng chảy.

Các khái niệm chuyên ngành then chốt được chuẩn hóa bao gồm: môđun dòng chảy năm (l/s/km²), hệ số biến sai dòng chảy Cv, lưu lượng kiệt thiết kế, nhu cầu cấp nước nông nghiệp ứng với tần suất bảo đảm P=75%, vùng xáo trộn pha trộn ô nhiễm (mixing zone), chỉ số oxy hòa tan (DO), nhu cầu oxy hóa sinh học (BOD5) và nhu cầu oxy hóa cacbon sinh hóa (CBOD).

Hệ thống mô hình toán thủy văn và chất lượng nước tích hợp gồm 4 công cụ chuyên sâu:

  • Mô hình thủy văn thông số tập trung NAM: Mô phỏng quá trình mưa - dòng chảy mặt, dòng chảy ngầm và bổ cập ẩm tầng đất.
  • Mô hình cân bằng nước MIKE BASIN: Mô phỏng và phân bổ tài nguyên nước đa mục tiêu giữa các ngành kinh tế theo cấu trúc nút và sơ đồ hệ thống.
  • Mô hình chất lượng nước QUAL2K: Mô phỏng biến đổi chất lượng nước một chiều (1D) dọc theo dòng chính với các thông số SS, DO, CBOD nhanh và CBOD chậm.
  • Mô hình CORMIX: Phân tích thủy lực pha trộn và trường lan truyền nồng độ ô nhiễm cục bộ tại các cửa xả công nghiệp ven bờ.

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

Nghiên cứu sử dụng chuỗi dữ liệu quan trắc lịch sử 52 năm (giai đoạn 1957 - 2009) thu thập từ 21 trạm đo mưa, 12 trạm khí tượng bề mặt và 16 trạm thủy văn trên toàn lưu vực sông Cầu. Trong đó, có 10 trạm quan trắc lưu lượng tiêu biểu như Thác Riềng, Thác Bưởi, Gia Bảy, Tân Cương, Giang Tiên, Cầu Mai và Phú Cường.

Phương pháp chọn mẫu dữ liệu được thực hiện theo nguyên tắc mạng lưới trạm đại biểu, bao phủ đồng đều 3 vùng địa hình đặc trưng: vùng núi cao thượng lưu (Bắc Kạn), vùng trung du đồi bát úp (Thái Nguyên) và vùng đồng bằng hạ lưu (Bắc Ninh, Bắc Giang, Vĩnh Phúc, Hà Nội). Nghiên cứu ứng dụng phương pháp đa giác Thiessen để tính toán lượng mưa bình quân lưu vực với trọng số trạm được xác định chi tiết (Vĩnh Yên 0,12; Phúc Yên 0,11; Bắc Kạn 0,10; Tam Đảo 0,10; Hiệp Hòa 0,09; Thái Nguyên 0,07; Đại Từ 0,07).

Lý do lựa chọn tổ hợp phương pháp này xuất phát từ thực trạng chuỗi đo đạc lưu lượng tại một số trạm quan trọng (như Thác Bưởi) bị gián đoạn sau năm 1996 do điều kiện địa hình phức tạp. Việc sử dụng mô hình NAM cho phép kéo dài chuỗi số liệu dòng chảy mặt một cách khoa học từ chuỗi mưa ngày liên tục. Mô hình MIKE BASIN, QUAL2K và CORMIX giúp giải quyết bài toán quản lý nước đa ngành, lượng hóa chính xác cả tải lượng ô nhiễm diện rộng lẫn nồng độ tích tụ tại các vùng nước tiếp nhận xung yếu.

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

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

Thứ nhất, lượng mưa trung bình nhiều năm trên lưu vực sông Cầu đạt 1.664 mm, dao động trong khoảng từ 1.186 mm đến 2.144 mm. Phân bố mưa theo không gian có sự phân hóa mạnh mẽ: khu vực núi cao Tam Đảo đạt giá trị lớn nhất 2.459 mm/năm, phần trung lưu đạt khoảng 1.900 mm/năm và hạ lưu chỉ đạt từ 1.400 đến 1.600 mm/năm. Lượng bốc hơi tiềm năng bình quân dao động từ 541 mm/năm (Tam Đảo) đến 1.011 mm/năm (Hà Nội, Vĩnh Yên đạt 998 mm/năm, Thái Nguyên đạt 975 mm/năm).

Thứ hai, đặc trưng dòng chảy năm có môđun bình quân toàn lưu vực đạt 21,4 l/s/km², với sự chênh lệch lớn giữa mùa lũ và mùa kiệt. Mùa lũ kéo dài từ tháng VI đến tháng IX, với đỉnh lũ lịch sử ghi nhận tại trạm Thác Bưởi lên tới 3.490 m³/s. Ngược lại, trong mùa cạn, môđun dòng chảy kiệt nhỏ nhất trên toàn lưu vực rơi xuống mức dưới 1,0 l/s/km². Điển hình trên sông Công, lưu lượng nhỏ nhất tại trạm Núi Hồng chỉ đạt 0,16 m³/s và sông Cầu Mai chỉ đạt 0,02 m³/s. Việc hồ Núi Cốc (dung tích 175 triệu m³) vận hành từ năm 1978 làm hạ lưu sông Công suy giảm dòng chảy tự nhiên nghiêm trọng trong mùa khô.

Thứ ba, kết quả tính toán cân bằng nước bằng mô hình MIKE BASIN cho thấy tổng nhu cầu sử dụng nước nông nghiệp với tần suất thiết kế P=75% chiếm trên 70% tổng lượng nước tiêu thụ. Tại các nút cấp nước thuộc khu vực hạ lưu đập Thác Huống và vùng hạ lưu Núi Cốc, vào các tháng I, II và III xảy ra tình trạng thiếu hụt nước cục bộ từ 18% đến 32% so với nhu cầu thực tế theo các kịch bản phát triển kinh tế 2010 - 2020.

Thứ tư, kết quả mô phỏng QUAL2K và CORMIX xác định nước mặt sông Cầu đoạn chảy qua thành phố Thái Nguyên và hạ lưu đập Thác Huống bị ô nhiễm hữu cơ và chất rắn lơ lửng rõ rệt. Nồng độ DO tại một số đoạn giảm sâu xuống dưới 3,5 mg/l, trong khi nồng độ BOD5 và SS tại các điểm tiếp nhận nước thải từ Khu gang thép Lưu Xá, Nhà máy cán thép Gia Sàng, mỏ sắt Trại Cau và Nhà máy giấy Hoàng Văn Thụ vượt ngưỡng quy chuẩn cho phép từ 1,5 đến 3,8 lần.

Thảo luận kết quả

Nguyên nhân cốt lõi dẫn đến mâu thuẫn tài nguyên nước trên lưu vực sông Cầu là sự suy giảm khả năng tự điều tiết của lưu vực. Độ che phủ rừng tại hai tỉnh thượng nguồn Bắc Kạn và Thái Nguyên đã giảm xuống mức 32% (so với mức 34,6% của giai đoạn trước), khiến diện tích đất trống đồi núi trọc tăng lên 35%, làm gia tăng tốc độ lũ mùa mưa và làm cạn kiệt nguồn sinh thủy mùa khô.

Dữ liệu nghiên cứu được biểu diễn trực quan thông qua đường lũy tích sai chuẩn dòng chảy và mưa, thể hiện rõ chu kỳ dao động 11 năm và các bội số của 11 năm với các pha nhiều nước, ít nước kéo dài từ 4 đến 24 năm. Biểu đồ diễn tiến nồng độ CBOD nhanh, CBOD chậm và DO dọc theo 288,5 km chiều dài sông Cầu trong mô hình QUAL2K đã làm sáng tỏ cơ chế tự làm sạch bị quá tải tại các điểm hợp lưu tiếp nhận nước thải công nghiệp.

So với các công trình nghiên cứu trên lưu vực sông Nhuệ - Đáy, tốc độ lan truyền ô nhiễm trên sông Cầu có tính phân vùng rõ rệt hơn: vùng thượng nguồn có chất lượng nước tương đối tốt, nhưng vùng trung lưu chịu tác động nặng nề bởi chất thải luyện kim và khai khoáng, trong khi hạ lưu lại chịu tải lượng ô nhiễm từ làng nghề và đô thị với hệ số tích tụ ô nhiễm nhanh hơn 1,3 lần vào mùa kiệt.

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

Nhằm giải quyết triệt để tình trạng thiếu hụt nguồn nước và kiểm soát suy thoái môi trường trên lưu vực sông Cầu, luận văn đưa ra 4 nhóm giải pháp trọng tâm:

  1. Hiện đại hóa mạng lưới quan trắc khí tượng thủy văn và chất lượng nước: Đầu tư xây dựng bổ sung 15 trạm quan trắc tự động liên tục đo lưu lượng và các chỉ tiêu DO, pH, SS, COD tại các vị trí xung yếu (trạm Thác Bưởi, Gia Bảy, ngã ba sông Công và sông Cà Lồ) hoàn thành trước năm 2026. Cơ quan chủ trì: Ủy ban nhân dân 6 tỉnh lưu vực phối hợp cùng Bộ Tài nguyên và Môi trường. Mục tiêu: Giảm sai số dự báo dòng chảy mùa kiệt xuống dưới 8% và cung cấp dữ liệu cảnh báo ô nhiễm tức thời.
  2. Thiết lập quy trình vận hành điều tiết liên hồ chứa và hệ thống thủy nông: Rà soát, ban hành quy trình xả dòng chảy môi trường tối thiểu từ hồ Núi Cốc (duy trì tối thiểu 2,5 m³/s) và đập dâng Thác Huống (duy trì tối thiểu 5,0 m³/s) trong giai đoạn các tháng I đến tháng IV hàng năm, thực hiện từ giai đoạn 2026 - 2028. Cơ quan chủ trì: Bộ Nông nghiệp và Phát triển Nông thôn phối hợp với Công ty TNHH MTV Khai thác Thủy lợi Thái Nguyên. Mục tiêu: Giảm tỷ lệ thiếu hụt nước tưới vụ Đông Xuân từ 32% xuống dưới 5%.
  3. Kiểm soát nghiêm ngặt nguồn thải công nghiệp và áp dụng công nghệ xử lý tại nguồn: Bắt buộc 100% các cơ sở luyện kim, khai thác khoáng sản (Khu gang thép Lưu Xá, mỏ Trại Cau, Gia Sàng) và các cụm công nghiệp hoàn thành hệ thống xử lý nước thải đạt QCVN 40:2011/BTNMT cột A trước năm 2027. Cơ quan chủ trì: Cục Kiểm soát Ô nhiễm Môi trường và Sở Tài nguyên và Môi trường các tỉnh. Mục tiêu: Cắt giảm 40% tải lượng BOD5 và 50% tổng lượng chất rắn lơ lửng xả vào dòng chính sông Cầu.
  4. Tăng cường phục hồi rừng phòng hộ đầu nguồn và hành lang bảo vệ nguồn nước: Thực hiện đề án trồng mới và khoanh nuôi tái sinh 12.000 ha rừng phòng hộ tại thượng nguồn Bắc Kạn, Thái Nguyên và vùng đệm Tam Đảo trong giai đoạn 2026 - 2030, nâng tỷ lệ che phủ rừng toàn lưu vực từ 32% lên trên 45%. Cơ quan chủ trì: Chi cục Kiểm lâm các tỉnh. Mục tiêu: Tăng khả năng trữ nước ngầm và nâng môđun dòng chảy kiệt thêm 15% đến 20%.

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

Tài liệu luận văn mang lại giá trị học thuật và thực tiễn sâu sắc cho 4 nhóm đối tượng chính:

  1. Cơ quan quản lý nhà nước và hoạch định chính sách: Lãnh đạo và chuyên viên tại Sở Tài nguyên và Môi trường, Sở Nông nghiệp và Phát triển Nông thôn của 6 tỉnh thuộc lưu vực có thể sử dụng kết quả mô phỏng MIKE BASIN và QUAL2K làm căn cứ khoa học để phê duyệt quy hoạch tài nguyên nước, cấp phép xả thải và phân bổ hạn ngạch khai thác nước mặt.
  2. Đơn vị tư vấn thiết kế công trình thủy lợi, cấp thoát nước: Các kỹ sư thủy văn và công trình có thể khai thác trực tiếp bảng phân phối mưa thiết kế tần suất P=75%, chuỗi số liệu dòng chảy lũ cực đại 3.490 m³/s và đặc trưng dòng chảy kiệt để tính toán kích thước công trình lấy nước, cống phân lũ và trạm bơm cấp nước sinh hoạt.
  3. Giảng viên, nhà nghiên cứu và học viên cao học: Giảng viên và học viên chuyên ngành Thủy văn học, Kỹ thuật Tài nguyên Nước, Quản lý Môi trường có thể tham khảo quy trình tích hợp chuỗi 4 mô hình toán (NAM, MIKE BASIN, QUAL2K, CORMIX) như một tài liệu giảng dạy mẫu mực về mô hình hóa lưu vực sông.
  4. Ban quản lý khu công nghiệp và doanh nghiệp sản xuất: Đội ngũ quản lý môi trường tại các nhà máy luyện kim, khai khoáng và chế biến có thể vận dụng mô hình CORMIX để tính toán vùng xáo trộn pha trộn, tối ưu hóa vị trí và chế độ xả thải ven bờ nhằm đáp ứng quy định pháp luật môi trường.

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

Tại sao cần phải kết hợp đồng thời mô hình NAM và mô hình MIKE BASIN trong nghiên cứu này? Do mạng lưới trạm đo lưu lượng trên lưu vực sông Cầu bị thiếu hụt sau năm 1996 tại các vị trí xung yếu như Thác Bưởi, mô hình NAM đóng vai trò chuyển đổi chuỗi số liệu mưa 52 năm thành chuỗi dòng chảy mặt liên tục. Sau đó, mô hình MIKE BASIN tiếp nhận dữ liệu này để tính toán phân bổ dòng chảy cho các vùng dùng nước trên diện tích 6.030 km², giúp phản ánh chính xác bài toán cân bằng nước lưu vực.

Nguyên nhân chính dẫn đến hiện tượng thiếu nước mùa khô trên lưu vực sông Cầu là gì? Nguyên nhân chủ yếu do sự phân bố mưa không đều khi 8 tháng mùa khô chỉ chiếm 24% tổng lượng dòng chảy năm. Vào các tháng I đến tháng III, môđun dòng chảy kiệt giảm xuống dưới 1,0 l/s/km², trùng với thời kỳ tưới cao điểm của vụ Đông Xuân, dẫn đến tình trạng thâm hụt nước cục bộ từ 18% đến 32% tại khu vực hạ lưu đập Thác Huống và hồ Núi Cốc.

Mô hình QUAL2K và CORMIX đánh giá mức độ ô nhiễm nước sông Cầu như thế nào? Mô hình QUAL2K chỉ ra rằng nước sông Cầu bị ô nhiễm hữu cơ nặng khi chảy qua khu vực đô thị và công nghiệp Thái Nguyên, làm hàm lượng DO giảm xuống dưới 3,5 mg/l. Mô hình CORMIX chứng minh vùng hòa trộn ô nhiễm kéo dài từ 500 m đến 1.200 m quanh các cống xả của Khu gang thép Lưu Xá và mỏ Trại Cau, với nồng độ SS và BOD5 vượt tiêu chuẩn từ 1,5 đến 3,8 lần.

Hồ Núi Cốc có vai trò và tác động ra sao đối với dòng chảy hạ lưu sông Công và sông Cầu? Hồ chứa Núi Cốc với dung tích 175 triệu m³ giữ vai trò cung cấp nước tưới và sinh hoạt thiết yếu cho tỉnh Thái Nguyên. Tuy nhiên, việc tích nước từ năm 1978 làm cho đoạn hạ lưu sông Công sau đập bị suy giảm dòng chảy tự nhiên trầm trọng trong mùa khô, lưu lượng kiệt tại trạm Núi Hồng có thời điểm chỉ còn 0,16 m³/s, làm suy yếu khả năng tự làm sạch của lòng dẫn.

Phương pháp đa giác Thiessen được triển khai như thế nào để xác định mưa bình quân lưu vực? Nghiên cứu đã thiết lập mạng lưới đa giác Thiessen từ 14 trạm đo mưa đại biểu với các trọng số cụ thể như Vĩnh Yên 0,12; Phúc Yên 0,11; Bắc Kạn 0,10; Tam Đảo 0,10; Hiệp Hòa 0,09; Thái Nguyên 0,07. Phương pháp này giúp khắc phục sai số địa hình phức tạp, phản ánh chính xác lượng mưa bình quân lưu vực 1.664 mm/năm phục vụ tính toán thủy văn.

Kết luận

  • Luận văn đã xử lý đồng bộ chuỗi dữ liệu khí tượng thủy văn 52 năm (1957 - 2009) trên toàn bộ lưu vực sông Cầu diện tích 6.030 km², xác định lượng mưa trung bình năm đạt 1.664 mm và môđun dòng chảy năm đạt 21,4 l/s/km².
  • Thiết lập thành công chuỗi công cụ mô hình thủy văn NAM nhằm khôi phục và kéo dài dòng chảy mặt tại 10 trạm lưu lượng, tạo tiền đề dữ liệu đầu vào chuẩn xác cho mô phỏng cân bằng nước.
  • Ứng dụng mô hình MIKE BASIN lượng hóa nguy cơ thiếu hụt nguồn nước mặt lên tới 32% vào mùa kiệt giai đoạn 2010 - 2020 tại các nút dùng nước trọng điểm hạ lưu Thác Huống và Núi Cốc.
  • Phân tích chi tiết trường nồng độ các chất ô nhiễm hữu cơ (DO, BOD5, CBOD) và chất rắn lơ lửng dọc 288,5 km dòng chính bằng mô hình QUAL2K và lan truyền cục bộ CORMIX.
  • Xây dựng hệ thống 4 nhóm giải pháp công trình và phi công trình gắn với mốc thời gian 2026 - 2030, hỗ trợ trực tiếp cho công tác quy hoạch phát triển bền vững lưu vực sông Cầu.

Các hướng nghiên cứu tiếp theo cần mở rộng tích hợp các kịch bản biến đổi khí hậu toàn cầu và bài toán kiểm soát ô nhiễm nước thải xuyên biên giới tỉnh giai đoạn 2026 - 2035. Các nhà quản lý, chuyên gia quy hoạch và kỹ sư công trình quan tâm có thể tra cứu toàn văn luận văn thạc sĩ của tác giả Mai Thị Ngân Anh để khai thác toàn bộ hệ thống bảng biểu thống kê và bộ thông số mô hình chi tiết.