Tổng quan nghiên cứu

Lưu vực sông Cầu giữ vị trí chiến lược đặc biệt quan trọng đối với sự phát triển kinh tế - xã hội của vùng kinh tế trọng điểm Bắc Bộ, bao phủ diện tích hơn 6.030 km² trải dài qua 5 tỉnh gồm Bắc Kạn, Thái Nguyên, Vĩnh Phúc, Bắc Giang, Bắc Ninh và thành phố Hà Nội. Với dân số toàn lưu vực đạt hơn 2,838 triệu người (năm 2007) và mật độ trung bình 487 người/km², trong đó cư dân nông thôn chiếm tới 80,01%, tài nguyên nước sông Cầu đang phải gánh vác áp lực cấp nước khổng lồ cho sinh hoạt, công nghiệp và diện tích sản xuất nông nghiệp lên tới 194.554 ha. Sự phân bố dòng chảy không đều theo không gian và thời gian tạo nên mâu thuẫn gay gắt giữa nguồn nước đến và nhu cầu sử dụng. Trong khi mùa mưa (từ tháng V đến tháng IX) tập trung 75% đến 80% tổng lượng mưa năm thì mùa kiệt (từ tháng X đến tháng IV năm sau) chỉ chiếm 20% đến 25% lượng dòng chảy, dẫn đến tình trạng cạn kiệt nguồn nước mặt nghiêm trọng tại nhiều tiểu khu hạ du.

Trước thực trạng khai thác thiếu đồng bộ và sự suy giảm chất lượng nước do tốc độ đô thị hóa, công nghiệp hóa nhanh chóng, nghiên cứu tập trung vào mục tiêu định lượng hóa cán cân cung - cầu nước trên toàn lưu vực. Bằng việc ứng dụng công nghệ mô hình hóa tiên tiến MIKE BASIN kết hợp mô đun mưa - dòng chảy NAM trên nền tảng Hệ thống Thông tin Địa lý (ArcView GIS), đề tài tiến hành phân tích hiện trạng cân bằng nước năm 2007 và dự báo bức tranh phân bổ nước theo quy hoạch phát triển đến năm 2015. Kết quả nghiên cứu cung cấp cơ sở khoa học chính xác, lượng hóa tỷ lệ thiếu hụt nước tại các nút tiêu thụ, qua đó hỗ trợ các nhà quản lý xây dựng chiến lược điều tiết công trình và phân bổ nguồn nước bền vững cho toàn 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 được xây dựng dựa trên nguyên lý cân bằng nước hệ thống, xuất phát từ định luật bảo toàn vật chất trong thủy văn học. Phương trình cân bằng nước tổng quát được thiết lập dưới dạng quan hệ giữa tổng lượng nước đến, lượng nước đi và sự biến thiên lượng trữ nước trong một đơn vị thời gian:

$$(X + Z_1 + Y_1 + W_1) - (Z_2 + Y_2 + W_2) = \Delta U$$

Trong đó, $X$ là lượng mưa rơi trên lưu vực; $Z_1, W_1, Y_1$ lần lượt là lượng nước ngưng tụ, dòng chảy ngầm và dòng chảy mặt chảy vào; $Z_2, Y_2, W_2$ là lượng bốc hơi, dòng chảy mặt và dòng ngầm chảy ra; $\Delta U$ là độ biến thiên lượng trữ nước. Khái niệm hệ thống nguồn nước được xem xét như một chỉnh thể phức hợp gồm ba thành phần cốt lõi: tài nguyên nước tự nhiên, hệ thống công trình khai thác điều tiết và các hộ tiêu thụ nước (nông nghiệp, công nghiệp, sinh hoạt, môi trường).

Để giải quyết bài toán cân bằng nước trên quy mô lớn, nghiên cứu lựa chọn bộ công cụ MIKE BASIN do Viện Thủy lực Đan Mạch (DHI) phát triển. MIKE BASIN vận hành trên cấu trúc mạng lưới dạng nhánh - nút (node-branch network) bao gồm 7 loại nút chức năng: nút biên nguồn (simple node), nút hứng nước (catchment node), nút chuyển nước (offtake node), nút tưới tiêu (irrigation node), nút cấp nước sinh hoạt - công nghiệp (water supply node), nút hồ chứa (reservoir node) và nút thủy điện (hydropower node). Đi kèm với đó là mô đun mưa - dòng chảy NAM, mô phỏng quá trình thủy văn liên tục qua hệ thống 5 bể chứa ẩm thẳng đứng: bể mặt ($U_{\max}$), bể tầng rễ cây ($L_{\max}$), bể nước ngầm tầng trên và tầng dưới, cùng hệ thống tham số dòng chảy mặt ($CQOF$), dòng sát mặt ($CQIF$) và dòng ngầm ($CKBF$).

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ừ mạng lưới quan trắc khí tượng thủy văn quốc gia trên lưu vực sông Cầu. Cỡ mẫu nghiên cứu bao gồm chuỗi số liệu thủy văn nhiều năm (1960–2002) tại 5 trạm thủy văn chủ lực: Thác Bưởi (22 năm số liệu lưu lượng liên tục), Thác Riềng (36 năm số liệu dòng chảy), Tân Cương (16 năm), Phú Cường (10 năm) và Giang Tiên (10 năm). Về dữ liệu khí tượng, chuỗi số liệu mưa ngày và bốc hơi 43 năm liên tục (1961–2002) được trích xuất từ 7 trạm khí tượng đại diện: Bắc Kạn, Định Hóa, Thái Nguyên, Hiệp Hòa, Bắc Giang, Bắc Ninh và Đại Từ. Phương pháp chọn mẫu chuỗi thời gian dài và bao phủ toàn diện từ thượng lưu, trung lưu đến hạ lưu đảm bảo nắm bắt đầy đủ các chu kỳ năm nhiều nước, năm trung bình và năm ít nước cực đoan ($P = 75%$ đến $90%$).

Phương pháp phân tích dựa trên sự tích hợp chặt chẽ giữa mô hình số trị và ArcView GIS. Bản đồ số hóa độ cao DEM và ranh giới tiểu lưu vực được sử dụng để phân chia toàn bộ sông Cầu thành các tiểu vùng tính toán chi tiết. Mô hình NAM được hiệu chỉnh và kiểm nghiệm tự động kết hợp thử sai dựa trên 4 hàm mục tiêu: cực tiểu hóa sai số tổng lượng dòng chảy, dạng đường quá trình, lưu lượng đỉnh lũ và lưu lượng kiệt. Độ hữu hiệu của mô hình được đánh giá thông qua chỉ số Nash-Sutcliffe ($R^2$) theo tiêu chuẩn của Tổ chức Khí tượng Thế giới (WMO), bảo đảm $R^2$ đạt từ 0,65 đến trên 0,85 (loại khá và tốt). Lý do lựa chọn giải pháp này là khả năng khôi phục chính xác chuỗi dòng chảy đến tại các lưu vực không có trạm đo và mô phỏng trực quan các kịch bản điều tiết hồ chứa phức tạp.

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

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

Thứ nhất, kết quả tính toán dòng chảy bằng mô hình NAM đã khôi phục thành công bức tranh thủy văn chi tiết tại các tiểu khu. Lưu lượng dòng chảy năm biến đổi mạnh theo không gian, với môđun dòng chảy năm bình quân toàn lưu vực đạt khoảng 21 đến 30 l/s/km². Sự phân phối dòng chảy theo mùa phản ánh rõ tính chất cực đoan: 4 tháng mùa lũ (tháng VI đến tháng IX) tạo ra lưu lượng đỉnh lũ lịch sử tại Thác Bưởi lên tới 3.490 m³/s, trong khi 8 tháng mùa kiệt chỉ đóng góp 20% đến 25% tổng lượng nước năm, khiến môđun dòng chảy kiệt nhỏ nhất ở nhiều nhánh sông tụt xuống dưới 1,5 l/s/km².

Thứ hai, cơ cấu nhu cầu sử dụng nước trên lưu vực bộc lộ sự áp đảo tuyệt đối của ngành nông nghiệp. Trong tổng diện tích đất canh tác 194.554 ha, diện tích trồng lúa và hoa màu chiếm 146.999 ha (76%), đòi hỏi khối lượng nước tưới khổng lồ tập trung từ tháng XII đến tháng IV năm sau. Các hộ tiêu thụ nước khác như công nghiệp (tập trung tại khu Gang thép Thái Nguyên, Sông Công, Nội Bài), sinh hoạt cho 2,838 triệu dân và nuôi trồng thủy sản (với 1.967 ha tại Hà Nội và 2.687 ha tại Thái Nguyên) tạo ra sức ép cạnh tranh nguồn nước gay gắt vào mùa khô.

Thứ ba, mô phỏng cân bằng nước hiện trạng năm 2007 bằng MIKE BASIN chỉ ra sự chênh lệch lớn giữa kịch bản có và không có công trình thủy lợi. Khi không có công trình điều tiết, hầu hết các tiểu khu trung và hạ lưu (Thái Nguyên, Bắc Giang, Bắc Ninh) đều rơi vào trạng thái thiếu nước nghiêm trọng trong mùa cạn với lượng thiếu hụt hàng chục triệu m³. Sự hiện diện của các hồ chứa hiện hữu như hồ Núi Cốc (dung tích 175 triệu m³), hồ Đại Lải (25 triệu m³) và hồ Xạ Hương (12,7 triệu m³) đã cắt giảm đáng kể tình trạng thiếu hụt, nâng tỷ lệ đáp ứng nước cho vùng đồng bằng lên trên 80%.

Thứ tư, phân tích phương án quy hoạch năm 2015 cho thấy việc gia tăng nhu cầu nước đô thị và công nghiệp sẽ tiếp tục gây áp lực lên hệ thống. Tuy nhiên, việc đưa hồ chứa Văn Lăng vào vận hành cùng với việc nâng cấp các trạm bơm và hệ thống chuyển nước đã giúp tái lập cân bằng dòng chảy, giảm mức độ căng thẳng nước tại các tiểu vùng trọng điểm xuống dưới mức báo động 10%.

Thảo luận kết quả

Nguyên nhân cốt lõi dẫn đến tình trạng thiếu nước cục bộ trên lưu vực sông Cầu xuất phát từ sự bất đối xứng giữa chế độ mưa tự nhiên và mùa vụ sản xuất. Địa hình dốc trên 10‰ ở thượng lưu làm dòng chảy lũ rút nhanh, trong khi hạ lưu có độ dốc đáy chỉ 0,1‰ và chịu tác động của thủy triều (biên độ triều tại Đáp Cầu từ 0,3 đến 0,8 m), gây khó khăn cho việc tự chảy lấy nước ngọt. Thêm vào đó, diện tích đất trống đồi núi trọc chiếm tới 39% diện tích đất lâm nghiệp trong quá khứ đã làm suy giảm khả năng giữ nước tự nhiên của thảm phủ thực vật, gia tăng hệ số biến sai dòng chảy ($C_v$ dao động từ 0,25 đến 0,35).

So sánh với các nghiên cứu cân bằng nước sử dụng mô hình WEAP hay IQQM tại lưu vực sông Mê Kông và sông Đáy, mô hình MIKE BASIN trên lưu vực sông Cầu thể hiện ưu thế vượt trội trong việc gắn kết dữ liệu không gian GIS với các nút tính toán chi tiết. Kết quả tính toán cân bằng nước được trình bày trực quan thông qua:

  • Sơ đồ hóa mạng lưới sông và các nút điều tiết phân bố theo lưu vực nhánh;
  • Biểu đồ phân bố cơ cấu sử dụng nước giữa các ngành kinh tế;
  • Bản đồ phân vùng thiếu nước GIS phân cấp theo màu sắc rõ rệt giữa hai kịch bản có và không có công trình điều tiết.

Các bảng tổng hợp cán cân nước theo từng tháng tại các tiểu khu đã chỉ ra chính xác thời điểm căng thẳng nhất rơi vào các tháng I, II và III, là căn cứ không thể thay thế cho công tác lập lịch xả nước hồ chứa.

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

Nhằm giải quyết triệt để bài toán thiếu hụt nguồn nước và nâng cao hiệu quả quản lý tài nguyên nước lưu vực sông Cầu, 4 nhóm giải pháp trọng tâm cần được thực thi đồng bộ:

  • Đẩy nhanh tiến độ xây dựng và tối ưu hóa vận hành hồ chứa đa mục tiêu: Hoàn thiện xây dựng và đưa hồ chứa Văn Lăng vào khai thác trước năm 2015, kết hợp nâng cấp, nạo vét hệ thống kênh dẫn của các hồ Núi Cốc, Đại Lải nhằm nâng tổng dung tích điều tiết hữu ích toàn lưu vực lên thêm 15% đến 20%. Đơn vị thực hiện: Bộ Nông nghiệp và Phát triển nông thôn phối hợp cùng Ủy ban nhân dân tỉnh Thái Nguyên và Vĩnh Phúc.
  • Chuyển đổi cơ cấu cây trồng và ứng dụng công nghệ tưới tiết kiệm nước: Rà soát, chuyển đổi khoảng 15% diện tích trồng lúa kém hiệu quả tại các vùng gò đồi sang cây rau màu ít tiêu thụ nước; áp dụng công nghệ tưới phun mưa, tưới nhỏ giọt nhằm cắt giảm 20% định mức dùng nước nông nghiệp đến năm 2018. Đơn vị thực hiện: Sở Nông nghiệp và Phát triển nông thôn các tỉnh trong lưu vực.
  • Tăng cường bảo vệ rừng đầu nguồn và phục hồi thảm phủ sinh thái: Triển khai các dự án trồng rừng phòng hộ tại các huyện Chợ Đồn, Bạch Thông (Bắc Kạn) và Định Hóa, Võ Nhai (Thái Nguyên), nâng tỷ lệ che phủ rừng tự nhiên lên trên 60% vào năm 2020 nhằm tăng khả năng tích trữ nước ngầm mùa kiệt. Đơn vị thực hiện: Chi cục Kiểm lâm và Sở Tài nguyên và Môi trường các tỉnh Bắc Kạn, Thái Nguyên.
  • Thiết lập hệ thống quan trắc tự động và quy trình vận hành liên hồ chứa theo thời gian thực: Lắp đặt mạng lưới đo mưa và mực nước tự động tại 13 nhánh sông chính, ban hành quy chế điều tiết liên hồ trước mùa khô năm 2016 nhằm duy trì dòng chảy môi trường tối thiểu sau công trình đạt trên 1,5 m³/s. Đơn vị thực hiện: Ủy ban Quản lý Lưu vực sông Cầu phối hợp Trung tâm Khí tượng Thủy văn Quốc gia.

Đố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 lại giá trị thực tiễn và học thuật cao cho 4 nhóm đối tượng chính:

  • Cơ quan quản lý nhà nước và nhà hoạch định chính sách: Cán bộ Cục Quản lý Tài nguyên nước, Viện Quy hoạch Thủy lợ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 tại 5 tỉnh lưu vực sông Cầu và Hà Nội sử dụng số liệu làm căn cứ cấp phép khai thác nguồn nước, phê duyệt quy hoạch thủy lợi và xây dựng quy chế bảo vệ môi trường nước lưu vực.
  • Kỹ sư quy hoạch thủy lợi và chuyên gia mô hình hóa: Các kỹ sư công trình thủy, chuyên gia tư vấn thiết kế tìm thấy trong luận văn quy trình chuẩn hóa về thiết lập sơ đồ mạng lưới MIKE BASIN, phương pháp nội suy số liệu mưa - dòng chảy bằng NAM và kỹ thuật tích hợp ArcView GIS.
  • Giảng viên, nghiên cứu sinh và sinh viên ngành Thủy văn - Môi trường: Giảng viên và người học thuộc các chuyên ngành Thủy văn học, Kỹ thuật Tài nguyên nước, Khoa học Môi trường tại Đại học Khoa học Tự nhiên, Đại học Thủy lợi tham khảo làm tài liệu giảng dạy, học tập chuyên sâu về phương pháp tính toán cân bằng nước lưu vực sông phức tạp.
  • Đơn vị quản lý, khai thác công trình thủy nông và cấp nước đô thị: Các Công ty Khai thác Công trình Thủy lợi (như Thủy nông Sông Cầu, Núi Cốc) và các nhà máy nước sạch sử dụng dữ liệu dự báo thiếu hụt nước để xây dựng kế hoạch tích nước, vận hành máy bơm và tối ưu hóa lịch phân phối nước theo từng thời đoạn mùa khô.

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

Mô hình MIKE BASIN có ưu điểm gì nổi bật so với các mô hình cân bằng nước khác khi áp dụng cho sông Cầu?
MIKE BASIN tích hợp trực tiếp trên nền tảng ArcView GIS, cho phép biểu diễn trực quan mạng lưới sông dạng nhánh - nút và xử lý linh hoạt 7 loại nút tính toán. Đối với lưu vực sông Cầu có mạng lưới 68 sông suối và địa hình phức tạp, mô hình giúp mô phỏng chính xác quan hệ giữa lượng nước đến tự nhiên và nhu cầu dùng nước đa mục tiêu tại từng tiểu vùng.

Mô đun NAM xử lý bài toán thiếu dữ liệu dòng chảy thực đo như thế nào?
Mô đun NAM sử dụng mô hình toán 5 bể chứa ẩm liên tục để chuyển hóa dữ liệu mưa và bốc hơi giai đoạn 1961–2002 thành dòng chảy lưu vực. Bằng việc hiệu chỉnh thông số tại các trạm có số liệu như Thác Bưởi, Tân Cương và kiểm nghiệm đạt chỉ số $R^2$ từ 0,65 đến trên 0,85, mô hình đã khôi phục tin cậy chuỗi dòng chảy đến tại các nhánh sông không có trạm quan trắc.

Ngành nào gây áp lực lớn nhất lên cán cân nước lưu vực sông Cầu trong mùa kiệt?
Ngành nông nghiệp tạo áp lực lớn nhất, chiếm trên 70% tổng nhu cầu nước toàn lưu vực phục vụ cho 146.999 ha đất trồng lúa và hoa màu. Cao điểm lấy nước tưới tập trung vào các tháng I, II và III, trùng hoàn toàn với thời kỳ dòng chảy sông Cầu suy giảm thấp nhất (chỉ chiếm dưới 25% lượng dòng chảy cả năm), gây căng thẳng nguồn nước gay gắt.

Hồ chứa Văn Lăng đóng vai trò gì trong kịch bản quy hoạch phát triển năm 2015?
Theo kết quả mô phỏng quy hoạch năm 2015, hồ chứa Văn Lăng đóng vai trò công trình điều tiết chiến lược, bổ sung lưu lượng trữ để cấp nước cho nông nghiệp và các khu công nghiệp mới phát triển. Sự tham gia của hồ Văn Lăng giúp giảm thiểu cơ bản diện tích thiếu nước tại các tiểu vùng hạ lưu, giữ tỷ lệ thiếu nước toàn hệ thống ở mức an toàn.

Làm thế nào để kết hợp giữa tính toán cân bằng nước và bảo vệ chất lượng nước sông Cầu?
Tính toán cân bằng nước là tiền đề xác định lưu lượng dòng chảy tối thiểu cần duy trì trong sông (dòng chảy môi trường). Kết hợp dòng chảy kiệt được mô phỏng từ MIKE BASIN với các quy chuẩn kiểm soát xả thải công nghiệp tại Thái Nguyên, Bắc Giang sẽ tăng cường khả năng tự làm sạch của sông Cầu, ngăn chặn nguy cơ ô nhiễm nguồn nước mặt.

Kết luận

  • Luận văn đã xây dựng thành công cơ sở dữ liệu số hóa hoàn chỉnh và thiết lập mô hình cân bằng nước hệ thống MIKE BASIN tích hợp mô đun NAM trên nền GIS cho 6.030 km² lưu vực sông Cầu.
  • Đã khôi phục và lượng hóa chi tiết chuỗi dòng chảy đến từ 43 năm dữ liệu khí tượng thủy văn, xác định rõ lưu lượng kiệt và đỉnh lũ lịch sử tại các nút trọng điểm.
  • Đánh giá toàn diện cán cân nước hiện trạng năm 2007, khẳng định vai trò sống còn của các công trình hồ chứa hiện hữu trong việc cắt giảm nguy cơ thiếu hụt nước mùa khô cho 2,838 triệu dân và 194.554 ha đất canh tác.
  • Mô phỏng và chứng minh tính khả thi của kịch bản quy hoạch năm 2015, trong đó công trình hồ chứa Văn Lăng là đòn bẩy then chốt giúp hạ thấp tỷ lệ thiếu nước toàn vùng xuống dưới 10%.
  • Đề xuất lộ trình tích hợp các giải pháp phi công trình (trồng rừng đầu nguồn, tưới tiết kiệm) và giải pháp công trình (vận hành liên hồ thời gian thực) hướng tới mục tiêu an ninh nguồn nước giai đoạn 2015–2020.

Để tiếp tục hoàn thiện và phát triển hướng nghiên cứu này, các cơ quan quản lý và nhà khoa học cần đẩy mạnh việc cập nhật kịch bản biến đổi khí hậu mới nhất vào mô hình MIKE BASIN, đồng thời mở rộng mô phỏng lan truyền chất lượng nước nhằm bảo vệ toàn diện tài nguyên nước lưu vực sông Cầu trong tương lai.