Tổng quan nghiên cứu

Lưu vực sông Lô chảy qua địa phận tỉnh Tuyên Quang có chiều dài khoảng 145 km với diện tích lưu vực trực tiếp đạt 2.090 km2 trên tổng diện tích tự nhiên toàn tỉnh là 5.868 km2. Hệ thống sông ngòi này đóng vai trò sống còn đối với an sinh và phát triển kinh tế, đóng góp tới 24,1% tổng lượng nước hàng năm cho toàn bộ lưu vực sông Hồng. Tuy nhiên, đặc điểm địa hình với hơn 73% là đồi núi dốc cùng mạng lưới sông suối dạng rẻ quạt khiến vùng hạ du thường xuyên đối mặt với các đợt lũ lụt ác liệt. Điển hình như trận lũ lịch sử tháng 8 năm 1971 ghi nhận lưu lượng đỉnh lũ lên đến 12.000 m3/s tại trạm Tuyên Quang, mực nước dâng cao 31,78 m làm ngập úng toàn bộ khu vực trung tâm đô thị từ 7 đến 9 ngày.

Mục tiêu cốt lõi của nghiên cứu là đánh giá toàn diện chế độ thủy văn dòng chảy mùa lũ, phân tích thực trạng ngập úng và xác định các nguyên nhân gây lũ trên hệ thống sông Lô đoạn qua tỉnh Tuyên Quang. Từ đó, đề tài xây dựng các phương án diễn toán thủy lực và đề xuất hệ thống giải pháp công trình kết hợp phi công trình đồng bộ. Phạm vi nghiên cứu tập trung vào không gian lưu vực sông Lô và sông Gâm thuộc tỉnh Tuyên Quang, xử lý chuỗi số liệu quan trắc khí tượng thủy văn kéo dài liên tục gần 50 năm từ năm 1960 đến năm 2010. Ý nghĩa thực tiễn của công trình thể hiện qua việc bảo vệ an toàn cho 36,5 km đê bao hiện hữu, bảo vệ hơn 70.195 ha đất sản xuất nông nghiệp và tạo hành lang pháp lý vững chắc cho mục tiêu tăng trưởng kinh tế công nghiệp đạt mức 27,2%/năm của địa phương.

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 các lý thuyết chuyên ngành kỹ thuật tài nguyên nước và quản lý thiên tai hiện đại:

  • Lý thuyết thủy động lực học dòng chảy hở một chiều không ổn định: Vận dụng hệ phương trình vi phân Saint-Venant gồm phương trình liên tục và phương trình động lượng để mô tả chính xác sự biến đổi mực nước và lưu lượng theo không gian và thời gian dọc theo các tuyến sông phân nhánh.
  • Lý thuyết quản lý rủi ro thiên tai và ngập lụt tích hợp: Tiếp cận bài toán phòng chống lũ theo hướng kết hợp hài hòa giữa kiểm soát dòng chảy bằng giải pháp công trình và thích ứng mềm bằng giải pháp phi công trình.
  • Lý thuyết động lực học lòng dẫn và bồi xói hạ lưu: Đánh giá sự suy giảm bùn cát lơ lửng và quá trình xói lở bờ sông sau khi xây dựng các đập thủy điện quy mô lớn.
  • Các khái niệm chuyên môn nền tảng được định lượng hóa gồm: Hành lang thoát lũ, Cao trình đỉnh đê thiết kế, Tần suất lũ thiết kế chu kỳ lặp 500 năm tương ứng tần suất 0,2%, Độ dốc thủy lực đáy sông lớn hơn 0,5% và Hệ số nhám lòng dẫn Manning.

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

Nghiên cứu tích hợp hệ phương pháp đa ngành dựa trên cơ sở dữ liệu quan trắc thực địa quy mô lớn:

  • Nguồn dữ liệu đầu vào: Cỡ mẫu quan trắc khí tượng thủy văn gồm chuỗi số liệu liên tục 45 năm từ năm 1961 đến năm 2005 tại 3 trạm khí tượng Tuyên Quang, Chiêm Hóa, Hàm Yên và 5 trạm đo mưa Na Hang, Sơn Dương, Hàm Yên, Chiêm Hóa, Tuyên Quang. Dữ liệu hình thái sông được trích xuất từ hàng trăm mặt cắt thực đo cùng thông số vận hành của hồ thủy điện Tuyên Quang có dung tích 2.240 triệu m3.
  • Phương pháp chọn mẫu: Phương pháp chọn mẫu toàn diện theo chuỗi thời gian liên tục 50 năm nhằm bảo đảm tính đại diện thống kê cực trị của chuỗi dòng chảy lũ và loại bỏ các sai số đo đạc cục bộ.
  • Phương pháp phân tích và mô hình hóa: Ứng dụng mô hình thủy lực thủy động lực học MIKE 11 do Viện Thủy lực Đan Mạch phát triển kết hợp với phần mềm hệ thống thông tin địa lý MapInfo để số hóa mạng sông và biên tập bản đồ ngập lụt.
  • Lý do lựa chọn mô hình: MIKE 11 có khả năng tính toán chính xác mạng lưới sông tự nhiên phức tạp gồm hệ thống sông Lô, sông Gâm và sông Hồng với nhiều công trình kiểm soát như hồ chứa Hòa Bình, Thác Bà và Tuyên Quang. Mô hình được hiệu chỉnh dựa trên trận lũ thực tế tháng 8 năm 1996 và kiểm định độc lập qua trận lũ lớn tháng 9 năm 1985 nhằm bảo đảm độ tin cậy tuyệt đối trước khi diễn toán các kịch bản tương lai.

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

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

Quá trình mô phỏng thủy lực và phân tích thực địa đã làm sáng tỏ 4 phát hiện khoa học quan trọng:

  • Phát hiện về nguồn sinh lũ chính: Sông Gâm là nhánh sông hợp lưu lớn nhất, đóng góp xấp xỉ 35% tổng lượng nước của toàn bộ lưu vực sông Lô tại ngã ba Khe Lau với lưu lượng trung bình năm đạt 442,76 m3/s, tương đương tổng lượng nước hàng năm 14.016 triệu m3. Sự tập trung dòng chảy nhanh từ sông Gâm tạo nên các đợt lũ kép có đỉnh nhọn, tốc độ dòng chảy cực lớn từ 3 m/s đến 5 m/s và chênh lệch biên độ mực nước mùa lũ lên tới 14,6 m tại trạm Tuyên Quang.
  • Phát hiện về khả năng phòng hộ của đê điều: Hệ thống 36,5 km đê bao hiện tại mới chỉ đáp ứng mức ngăn lũ tương ứng cao độ mực nước cos +30,0 m tại thành phố Tuyên Quang. Kết quả mô phỏng trận lũ thiết kế tần suất 0,2% cho thấy nếu không có sự phối hợp cắt lũ của hồ Tuyên Quang và hồ Thác Bà, mực nước lũ sẽ vượt qua cao trình đỉnh đê từ 1,2 m đến 1,8 m trên diện rộng.
  • Phát hiện về nguy cơ xói lở hạ lưu sau đập: Công trình thủy điện Tuyên Quang với chiều cao đập 125 m đã giữ lại phần lớn lượng phù sa bùn cát, làm thay đổi mạnh mẽ chế độ thủy lực dòng chảy hạ du. Tình trạng xói lở bờ sông diễn ra đặc biệt nghiêm trọng tại các xã Thái Long, An Khang và Sầm Dương. Tại thôn Thái Thịnh, hiện tượng sạt lở đã tạo nên vết nứt đê kéo dài 130 m với bề rộng nứt từ 0,1 m đến 0,2 m, đe dọa trực tiếp đến an toàn các cống tiêu Sài Lĩnh và Bè Cạn.
  • Phát hiện về tác động của biến đổi khí hậu: Kịch bản biến đổi khí hậu dự báo lượng mưa cực đoan gia tăng, làm lưu lượng đỉnh lũ thiết kế 0,2% tăng thêm từ 8,5% đến 12% so với kịch bản nền, đòi hỏi mở rộng hành lang thoát lũ tại các điểm co thắt qua huyện Hàm Yên và thành phố Tuyên Quang. Toàn bộ các kết quả này được trực quan hóa qua biểu đồ đường quá trình mực nước thực đo so với tính toán mô phỏng và bảng ma trận lưu lượng đỉnh lũ giữa các kịch bản.

Thảo luận kết quả

Nguyên nhân ngập úng tại Tuyên Quang mang tính chất tổng hợp giữa yếu tố tự nhiên và hoạt động nhân sinh. Mùa mưa kéo dài từ tháng 5 đến tháng 9 chiếm tới 75% đến 80% tổng lượng mưa cả năm, trong đó tháng mưa cao nhất chiếm 20%. Khi mực nước sông Lô dâng cao, các ngòi tiêu tự nhiên như Ngòi Là, Ngòi Cơi, Ngòi Chả bị nghẽn dòng chảy ngược vào nội đồng. Tuyến đê Cấp Tiến dài 7,5 km dù đã bảo vệ được 400 ha hoa màu nhưng các cống tiêu trọng yếu như Ngòi Cát, Ngòi Liễn, Ngòi Khổng chưa được đầu tư đồng bộ, dẫn đến tình trạng ngập úng kéo dài.

So sánh với các nghiên cứu quốc tế trên sông Vàng tại Trung Quốc hay hệ thống đập kiểm soát lũ tại lưu vực sông Itajaí-Açu ở Brazil, việc phòng chống lũ không thể chỉ phụ thuộc hoàn toàn vào giải pháp tôn cao đê mà bắt buộc phải kết hợp vận hành dung tích phòng lũ của các hồ chứa thượng nguồn với việc bảo vệ không gian thoát lũ tự nhiên. Kết quả nghiên cứu khẳng định tính đúng đắn của việc xác định chỉ giới thoát lũ, tạo nền tảng khoa học vững chắc để địa phương vừa kiểm soát rủi ro thiên tai vừa hiện thực hóa mục tiêu chuyển dịch cơ cấu nông nghiệp, nâng tỷ trọng ngành chăn nuôi lên 45% vào năm 2020.

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

Nhằm giải quyết triệt để nguy cơ ngập lụt và sạt lở bờ sông, nghiên cứu đề xuất 4 nhóm giải pháp chiến lược:

  • Nâng cấp và kiên cố hóa hệ thống đê điều: Tiến hành gia cố toàn bộ 36,5 km đê bao sông Lô đạt cao trình chống lũ an toàn cos +32,0 m; hoàn thành dứt điểm 7,5 km tuyến đê Cấp Tiến và xây mới đồng bộ các cống tiêu Ngòi Cát, Ngòi Liễn, Ngòi Khổng trước năm 2015. Chủ thể thực hiện là Sở Nông nghiệp và Phát triển nông thôn tỉnh Tuyên Quang phối hợp với Ủy ban nhân dân các huyện Sơn Dương, Yên Sơn.
  • Xây dựng công trình kè chỉnh trị chống xói lở bờ sông: Xử lý khẩn cấp đoạn nứt đê dài 130 m tại xã Sầm Dương; triển khai xây dựng 21,8 km kè bờ sông suối xung yếu gồm 13 km kè bảo vệ khu hành chính Lâm Bình, 3 km kè Bình Xa, 3 km kè Minh Hương và 2,8 km kè Yên Phú trong giai đoạn 2012-2016. Chủ thể thực hiện là Chi cục Thủy lợi và Ban Quản lý dự án Đầu tư xây dựng công trình nông nghiệp tỉnh Tuyên Quang.
  • Tối ưu hóa quy trình điều tiết liên hồ chứa: Vận hành khoa học dung tích phòng lũ 2.240 triệu m3 của hồ Tuyên Quang phối hợp nhịp nhàng với hồ Thác Bà và hồ Hòa Bình, cắt giảm từ 25% đến 30% lưu lượng đỉnh lũ chảy về hạ du khi xảy ra các trận lũ tần suất 0,2%. Chủ thể thực hiện là Công ty Thủy điện Tuyên Quang dưới sự điều phối của Ban Chỉ đạo Quốc gia về Phòng chống thiên tai.
  • Hoàn thiện công tác quản lý phi công trình và cắm mốc hành lang thoát lũ: Ban hành quy định quản lý chỉ giới thoát lũ sông Lô qua thành phố Tuyên Quang và huyện Hàm Yên, hoàn thành cắm mốc 100% các đoạn sông trọng điểm trước năm 2014; thiết lập mạng lưới quan trắc cảnh báo sớm tự động đạt độ bao phủ tin cậy trên 95% số hộ dân ven sông. Chủ thể thực hiện là Sở Tài nguyên và Môi trường cùng Ủy ban nhân dân thành phố Tuyên Quang.

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

Nội dung và kết quả tính toán của luận văn đem lại giá trị ứng dụng thiết thực cho 4 nhóm đối tượng trọng tâm:

  • Cơ quan quản lý nhà nước và cơ quan hoạch định quy hoạch: Sở Nông nghiệp và Phát triển nông thôn, Sở Kế hoạch và Đầu tư, Sở Xây dựng sử dụng cơ sở dữ liệu thủy lực để thẩm định quy hoạch phát triển đô thị, xây dựng kết cấu hạ tầng giao thông và phân vùng sử dụng đất an toàn trước rủi ro ngập úng.
  • Ban Chỉ huy Phòng chống thiên tai và Tìm kiếm cứu nạn các cấp: Ứng dụng bản đồ ngập lụt và các kịch bản lũ cực hạn tần suất 0,2% để xây dựng phương án ứng phó khẩn cấp, phân vùng sơ tán dân cư an toàn tại các vùng trũng thấp ven sông Lô và sông Gâm.
  • Kỹ sư thiết kế thủy lợi, thủy điện và đơn vị tư vấn xây dựng: Khai thác trực tiếp các thông số kỹ thuật về lưu lượng đỉnh lũ, mực nước lớn nhất và vận tốc dòng chảy xói lở từ 3 m/s đến 5 m/s để tính toán kết cấu ổn định thân đê, mố kè và cống dưới đê.
  • Các nhà khoa học, giảng viên, học viên cao học và sinh viên chuyên ngành Kỹ thuật tài nguyên nước, Thủy văn, Biến đổi khí hậu: Sử dụng luận văn như một tài liệu tham khảo mẫu mực về quy trình thiết lập, hiệu chỉnh và kiểm định mô hình thủy lực MIKE 11 cho các lưu vực sông miền núi phức tạp.

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

  • Nguyên nhân cơ bản nào gây ngập lụt nghiêm trọng tại tỉnh Tuyên Quang? Tình trạng ngập lụt xuất phát từ lượng mưa lớn tập trung 75% đến 80% trong mùa lũ, kết hợp với lượng nước sông Gâm chiếm xấp xỉ 35% lưu lượng sông Lô hợp lưu tại Khe Lau. Đồng thời, hệ thống cống dưới đê chưa đồng bộ và lòng dẫn bị co hẹp khiến nước lũ tiêu thoát chậm, gây ngập úng từ 7 đến 9 ngày như trận lũ lịch sử năm 1971.

  • Mô hình thủy lực MIKE 11 đóng vai trò kỹ thuật gì trong nghiên cứu? Mô hình MIKE 11 thực hiện diễn toán một chiều dòng chảy không ổn định trên toàn bộ hệ thống sông Lô, sông Gâm và sông Hồng. Sau khi được kiểm định chính xác qua các trận lũ lịch sử năm 1985 và năm 1996, mô hình cung cấp các thông số mực nước và lưu lượng lũ ứng với tần suất 0,2% làm căn cứ thiết kế công trình đê kè.

  • Hồ thủy điện Tuyên Quang tác động như thế nào đến hạ du sông Lô? Với dung tích 2.240 triệu m3 và chiều cao đập 125 m, hồ Tuyên Quang có vai trò quyết định trong việc cắt giảm đỉnh lũ cho hạ du. Tuy nhiên, việc giữ lại phù sa bùn cát trong lòng hồ đã làm thay đổi động lực dòng chảy, gây xói lở mạnh bờ sông tại xã Thái Long, Sầm Dương với vết nứt đê dài 130 m.

  • Tần suất lũ thiết kế 0,2% trong luận văn có ý nghĩa gì đối với an toàn đê điều? Lũ tần suất 0,2% là đợt lũ kiểm tra tương ứng với chu kỳ lặp lại 500 năm một lần theo tiêu chuẩn phòng lũ quốc gia. Việc tính toán mực nước lũ ứng với kịch bản này giúp các kỹ sư xác định chính xác cao trình đỉnh đê cần nâng cấp lên mức cos +32,0 m nhằm ngăn chặn nguy cơ vỡ đê thảm khốc.

  • Những giải pháp phi công trình trọng tâm nào được đề xuất thực hiện? Giải pháp phi công trình tập trung vào việc xác định và cắm mốc hành lang thoát lũ sông Lô qua thành phố Tuyên Quang và huyện Hàm Yên. Bên cạnh đó, nghiên cứu kiến nghị số hóa bản đồ rủi ro ngập lụt, kiểm soát việc xây dựng xâm lấn bãi bồi và tăng cường độ che phủ rừng phòng hộ ổn định trên mức 60%.

Kết luận

  • Luận văn đã phân tích toàn diện đặc điểm khí tượng thủy văn lưu vực sông Lô dài 145 km qua tỉnh Tuyên Quang, làm rõ cơ chế tập trung lũ từ sông Gâm chiếm 35% lưu lượng toàn dòng.
  • Thiết lập và kiểm định thành công mô hình thủy động lực học MIKE 11 trên chuỗi số liệu thực đo gần 50 năm, diễn toán chính xác các kịch bản lũ thiết kế tần suất 0,2% có xét đến biến đổi khí hậu.
  • Đánh giá sâu sắc hiện trạng 36,5 km đê bao, chỉ rõ các điểm xung yếu sạt lở bờ sông tại Sầm Dương dài 130 m và tình trạng hạ thấp lòng dẫn sau khi vận hành hồ Tuyên Quang dung tích 2.240 triệu m3.
  • Đề xuất hệ thống giải pháp công trình và phi công trình đồng bộ, ưu tiên hoàn thiện 7,5 km đê Cấp Tiến, xây dựng 21,8 km kè bảo vệ và quản lý nghiêm ngặt hành lang thoát lũ.
  • Kết quả nghiên cứu là tài liệu khoa học nền tảng phục vụ đắc lực cho công tác quy hoạch phát triển kinh tế - xã hội, bảo vệ vững chắc thành quả sản xuất nông nghiệp và ổn định đời sống nhân dân vùng hạ du tỉnh Tuyên Quang.

Để nâng cao năng lực phòng chống thiên tai và ứng phó hiệu quả với biến đổi khí hậu trong giai đoạn tới, các cơ quan quản lý và đơn vị chuyên môn cần sớm triển khai áp dụng các kết quả tính toán thủy lực của đề tài vào thực tiễn quản lý quy hoạch xây dựng và vận hành công trình đê điều tại địa phương.