Tổng quan nghiên cứu

Thành phố Hà Nội nằm tại trung tâm vùng đồng bằng sông Hồng với tổng diện tích tự nhiên 334.470 ha trải rộng trên 29 quận, huyện, thị xã. Vị trí địa lý đặc thù khiến Thủ đô trở thành điểm hội tụ của các hệ thống sông lớn nhất miền Bắc trước khi đổ ra biển, đặc biệt là hợp lưu của ba con sông lớn gồm sông Đà, sông Thao và sông Lô tại khu vực Sơn Tây. Sông Hồng chảy qua Hà Nội với chiều dài khoảng 118 km, lưu lượng dòng chảy bình quân năm đạt 2.640 m³/s với tổng lượng nước trung bình khoảng 83,5 tỷ m³. Đồng thời, sông Đuống giữ vai trò là phân lưu chính với chiều dài 22 km chảy qua địa phận huyện Gia Lâm và Đông Anh, chuyển tải trung bình 27 tỷ m³ nước mỗi năm sang hệ thống sông Thái Bình.

Vấn đề nghiên cứu trọng tâm xuất phát từ thực trạng tốc độ đô thị hóa nhanh chóng khiến năng lực tiêu thoát nước của đô thị Hà Nội chỉ đáp ứng được khoảng 35% nhu cầu thực tế. Song song với đó, tình trạng lấn chiếm lòng dẫn và bãi sông diễn ra nghiêm trọng với 41.488 hộ dân sinh sống ngoài đê sông Hồng đoạn từ Thượng Cát đến Vạn Phúc và hơn 2.000 hộ ven sông Đuống. Cao trình các tuyến đê bối tự phát bị tôn cao vượt mức Báo động 3, làm co hẹp mặt cắt thoát lũ và gia tăng nguy cơ ngập úng nghiêm trọng cho toàn vùng hạ du khi xuất hiện các trận lũ lớn.

Mục tiêu cụ thể của luận văn là xây dựng cơ sở khoa học và thực tiễn vững chắc, ứng dụng các mô hình thủy văn và thủy lực hiện đại để đề xuất hệ thống giải pháp công trình và phi công trình nhằm tăng cường khả năng tiêu thoát nước của sông Hồng và sông Đuống. Nghiên cứu được hoàn thành vào tháng 9 năm 2012, hướng tới bảo vệ an toàn tuyệt đối cho trung tâm kinh tế - chính trị của cả nước với tần suất lũ thiết kế chu kỳ 250 năm (P = 0,4%) và 500 năm (P = 0,2%), đảm bảo mực nước lũ tại trạm Thủy văn Hà Nội không vượt quá cao trình 13,40 m.

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 thủy lực dòng chảy hở không ổn định biến đổi chậm và lý thuyết thủy văn lưu vực đồng bằng châu thổ. Chế độ thủy lực dòng chảy một chiều dọc theo mạng sông được mô tả bằng hệ phương trình vi phân đạo hàm riêng Saint-Venant gồm phương trình liên tục và phương trình động lượng, phản ánh sự biến đổi của lưu lượng và mực nước theo không gian và thời gian dưới tác động của gia tốc trọng trường, độ nhám lòng dẫn và lực ma sát đáy.

Đối với các vùng bãi tràn ngập lũ phức tạp hai bên bờ sông Hồng, lý thuyết thủy lực dòng chảy hai chiều phẳng được áp dụng nhằm mô tả chính xác sự tương tác dòng chảy ngang và dòng chủ lưu thông qua hệ phương trình vi phân đạo hàm riêng hai chiều có xét đến hệ số phân bố động lượng, lực Coriolis và ứng suất tiếp tuyến mặt cắt. Nghiên cứu định hình 4 khái niệm cốt lõi:

  • Chỉ giới thoát lũ: Ranh giới không gian bắt buộc phải giải tỏa và duy trì thông thoáng để đảm bảo chuyển tải an toàn lưu lượng lũ thiết kế.
  • Đê bối: Tuyến đê bao bảo vệ vùng bãi sông nông nghiệp với tiêu chuẩn khống chế cao trình ở mức Báo động 2 và cho phép tràn chủ động khi lũ vượt Báo động 3.
  • Tỷ lệ phân lưu mùa lũ: Phần trăm lưu lượng đỉnh lũ từ dòng chính sông Hồng được chuyển hướng sang các sông phân lưu như sông Đuống (28 - 30%) và sông Luộc (10 - 14%).
  • Dung tích cắt lũ hạ du: Khối tích phòng lũ chủ động tại các hồ chứa thủy điện thượng nguồn nhằm hạ thấp đỉnh lũ qua Hà Nội.

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

Phương pháp nghiên cứu kết hợp chặt chẽ giữa thu thập số liệu thực địa, phân tích thống kê toán học và mô phỏng số trị trên máy tính. Dữ liệu địa hình sử dụng bình đồ đo đạc tỷ lệ 1/5.000 năm 2000 kết hợp số liệu cập nhật mặt cắt ngang sông Hồng năm 2006. Chuỗi số liệu khí tượng thủy văn được thu thập liên tục qua nhiều thập kỷ tại các trạm quan trắc tiêu biểu như trạm Hà Nội (lượng mưa 1 ngày lớn nhất 211,5 mm), Gia Lâm (216,6 mm), Đông Anh (217,2 mm) cùng tài liệu quan trắc dòng chảy từ các trạm thủy văn Sơn Tây, Hòa Bình, Yên Bái và Phù Ninh.

Mô hình thủy lực một chiều MIKE 11 được lựa chọn để mô phỏng toàn bộ mạng lưới sông Hồng - sông Thái Bình gồm 12 tuyến sông chính với tổng chiều dài hàng trăm kilômét, kết hợp mô hình thủy văn NAM tính toán lưu lượng gia nhập từ 16 tiểu lưu vực khu giữa. Mô hình MIKE 21 HD với lưới ô vuông kích thước 50 x 50 m được sử dụng để tính toán chi tiết chế độ thủy lực hai chiều cho đoạn sông Hồng qua Hà Nội. Phương pháp số sử dụng thuật toán sai phân hữu hạn sơ đồ ẩn với phương pháp quét kép (Double Sweep) giúp đảm bảo độ ổn định và chính xác cao. Toàn bộ mô hình được hiệu chỉnh và kiểm nghiệm thực tế qua trận lũ lịch sử tháng 8/1996 với mức độ tương thích cao về đường quá trình mực nước và lưu lượng đỉnh lũ.

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

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

Phân tích thủy văn và thủy lực hệ thống sông đã chỉ ra 3 phát hiện khoa học quan trọng:

Thứ nhất, tổ hợp dòng chảy lũ sông Hồng chịu sự chi phối áp đảo từ sông Đà với tỷ lệ đóng góp lượng lũ 8 ngày lớn nhất về Sơn Tây đạt trung bình 49,2% (từng đạt đỉnh 68,8% vào năm 1964), tiếp theo là sông Lô đóng góp 28,2% và sông Thao đóng góp 19,0%. Về tỷ lệ phân phối dòng chảy châu thổ, sông Đuống gánh vác 28 - 30% lưu lượng mùa lũ từ sông Hồng (tương đương 6.865 m³/s tại Cửa Đuống trong trận lũ thiết kế), đóng vai trò then chốt trong việc giảm tải áp lực cho hạ du sông Hồng.

Thứ hai, quá trình bồi lắng lòng dẫn và lấn chiếm hành lang bãi sông diễn ra nghiêm trọng làm co hẹp cục bộ các mặt cắt thoát lũ trọng điểm, đặc biệt tại khu vực Cầu Thăng Long (bề rộng lòng sông thu hẹp còn 910 m kể cả trụ cầu), Cầu Chương Dương (lòng chính 500 m) và đoạn Long Biên - Bồ Đề. Sự tồn tại của 19 tuyến đê bối bờ tả dài 144,14 km và 19 tuyến đê bối bờ hữu dài 78,51 km với cao trình tiệm cận Báo động 3 đã cản trở nghiêm trọng không gian thoát lũ tự nhiên.

Thứ ba, sự tham gia của hệ thống liên hồ chứa thượng nguồn gồm Hòa Bình (dung tích cắt lũ 4,9 tỷ m³), Thác Bà (0,45 tỷ m³) và Tuyên Quang (1,0 tỷ m³) theo Quyết định số 80/2007/QĐ-TTg cho phép cắt giảm đáng kể đỉnh lũ tự nhiên 42.000 m³/s tại Sơn Tây. Kết quả tính toán thủy lực xác lập lưu lượng lũ thiết kế qua trạm Thủy văn Hà Nội là 20.578 m³/s tương ứng mực nước 13,40 m, và tại trạm Thượng Cát là 27.538 m³/s tương ứng mực nước 14,65 m.

Thảo luận kết quả

Hiện tượng ngập úng gia tăng tại Hà Nội không chỉ do lượng mưa cực đoan tập trung 85% vào mùa mưa mà còn do tác động cộng hưởng của hiện tượng sụt lún kiến tạo trẻ Holoxen với tốc độ hạ thấp địa hình lên tới 8,5 mm/năm. Thêm vào đó, diện tích mặt nước điều hòa đô thị bị san lấp hơn 30% kể từ năm 1995 đã làm suy giảm nghiêm trọng khả năng tự tiêu thoát nước cục bộ.

So sánh với các phương án quy hoạch truyền thống vốn chủ yếu dựa vào việc nâng cao thân đê chính, giải pháp tích hợp điều tiết liên hồ chứa cắt lũ kết hợp phân lũ qua cống Vân Cốc và đập Đáy (với 6 cửa van cung kích thước 33,75 m mỗi cửa) mang lại tính khả thi và an toàn vượt trội. Dữ liệu thủy lực được hệ thống hóa chi tiết qua bảng quan hệ mực nước - lưu lượng thiết kế tại các mặt cắt:

Vị trí mặt cắt Tuyến sông Lưu lượng lũ thiết kế (m³/s) Mực nước lũ thiết kế (m)
Thượng Cát Sông Hồng 27.538 14,65
Cầu Thăng Long Sông Hồng 27.471 14,22
Trạm TV Hà Nội Sông Hồng 20.578 13,40
Cửa Đuống Sông Đuống 6.865 13,34
Cầu Phù Đổng Sông Đuống 6.853 11,76

Bản đồ phân bố vận tốc dòng chảy 2D từ mô hình MIKE 21 thể hiện rõ hiện tượng xoáy cuộn và nghẽn dòng tại các doi cát, bãi bồi như bãi Minh Châu, Võng La, Tự Nhiên, khẳng định tính cấp thiết phải nạo vét và chỉnh trị hình học lòng dẫn.

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

Dựa trên kết quả mô phỏng thủy lực, luận văn đề xuất 4 nhóm giải pháp đồng bộ:

  • Thiết lập và cắm mốc chỉ giới thoát lũ thực địa: Tiến hành số hóa và cắm mốc chỉ giới thoát lũ trên thực địa với mật độ trung bình 100 m/mốc dọc hai bờ sông Hồng và sông Đuống. Ủy ban nhân dân thành phố Hà Nội chủ trì phối hợp với Sở Nông nghiệp và Phát triển nông thôn hoàn thành việc cắm mốc trong giai đoạn 2012 - 2015 nhằm ngăn chặn triệt để tình trạng tái lấn chiếm bãi sông.
  • Nạo vét lòng dẫn và hạ thấp đê bối cản lũ: Thực hiện nạo vét mở rộng mặt cắt thoát lũ tại các bãi sông trọng điểm gồm bãi Đại Mạch - Võng La, Hải Bối, Tàm Xá - Xuân Canh bờ tả và bãi Thượng Cát, Tứ Liên, Phúc Xá bờ hữu. Điều chỉnh cao trình các tuyến đê bối về đúng mức quy định tại Nghị định số 62/1999/NĐ-CP (không vượt quá mức Báo động 2), đảm bảo lòng sông thoát an toàn lưu lượng 27.538 m³/s trong giai đoạn 2015 - 2020.
  • Kiên cố hóa công trình đê điều và kè chỉnh trị: Nâng cấp đồng bộ 37,709 km đê cấp đặc biệt hữu Hồng và nâng cấp đê tả Đuống lên cấp I. Gia cố, cứng hóa mặt cắt đê kết hợp xử lý dứt điểm các điểm thẩm lậu, mạch sủi tại 84 tuyến kè và hệ thống cống dưới đê hiện có; duy trì 54,30 km dải tre chắn sóng bảo vệ mái đê trước mùa mưa bão.
  • Di dời dân cư vùng lòng dẫn và bảo vệ rừng phòng hộ: Tổ chức di dời 11.954 hộ dân nằm trong chỉ giới thoát lũ nguy hiểm (gồm 10.548 hộ ven sông Hồng và 1.406 hộ ven sông Đuống). Đồng thời, Bộ Nông nghiệp và Phát triển nông thôn phối hợp các tỉnh thượng nguồn đẩy mạnh trồng rừng phòng hộ, nâng tỷ lệ che phủ rừng từ mức 7,58% hiện tại lên trên 15% vào năm 2030 để giảm dòng chảy mặt lũ quét về hạ du.

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

Nội dung luận văn mang lại giá trị thực tiễn và học thuật cho 4 nhóm đối tượng chính:

  • Cơ quan quản lý nhà nước và Ban chỉ huy phòng chống thiên tai: Cung cấp cơ sở khoa học để xây dựng quy hoạch phân lũ, quản lý không gian bãi sông và lập phương án ứng phó khẩn cấp khi mực nước sông Hồng tại Hà Nội chạm ngưỡng 13,40 m.
  • Các đơn vị tư vấn thiết kế công trình thủy lợi - giao thông thủy: Sử dụng trực tiếp các thông số thủy lực, mực nước lũ thiết kế (14,65 m tại Thượng Cát; 13,34 m tại Cửa Đuống) để tính toán khẩu độ thoát nước cho cầu vượt sông, kè chống sạt lở và cống qua đê.
  • Chuyên gia quy hoạch đô thị và quản lý tài nguyên đất: Ứng dụng bản đồ chỉ giới thoát lũ để hoạch định ranh giới xây dựng các khu đô thị ven sông, đảm bảo không gian trữ nước tự nhiên và bảo vệ 334.470 ha diện tích Thủ đô khỏi nguy cơ ngập lụt.
  • Giảng viên, nghiên cứu sinh và sinh viên ngành Kỹ thuật tài nguyên nước: Nắm bắt phương pháp luận ghép nối mô hình thủy lực 1D MIKE 11 và 2D MIKE 21 HD trong giải quyết bài toán thủy lực mạng sông phức tạp có bãi tràn rộng lớn.

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

Tại sao đoạn sông Hồng qua Hà Nội lại là trọng điểm thoát lũ xung yếu nhất toàn lưu vực?

Đoạn sông qua Hà Nội là nơi hợp lưu trực tiếp của ba nhánh sông lớn gồm sông Đà, sông Thao và sông Lô. Toàn bộ lượng nước khổng lồ từ lưu vực miền núi phía Bắc dồn về Sơn Tây với lưu lượng lũ tự nhiên có thể vượt 42.000 m³/s, tạo áp lực thủy lực cực lớn lên hệ thống đê hữu Hồng bảo vệ vùng lõi nội thành có mật độ dân số cao.

Tỷ lệ phân lưu dòng chảy giữa sông Hồng và sông Đuống có ý nghĩa thế nào đối với công tác chống lũ?

Sông Đuống đóng vai trò là van xả lũ tự nhiên quan trọng nhất, tiếp nhận 28 - 30% tổng lượng dòng chảy mùa lũ từ sông Hồng (chuyển khoảng 27 tỷ m³ nước mỗi năm sang sông Thái Bình). Việc duy trì tỷ lệ phân lưu hợp lý giúp hạ thấp đáng kể mực nước lũ tại trạm Hà Nội nhưng đòi hỏi phải kiên cố hóa hệ thống đê sông Đuống dài 22 km.

Hệ thống hồ chứa thượng nguồn giúp cắt giảm đỉnh lũ cho Thủ đô Hà Nội ra sao?

Các hồ chứa Hòa Bình (cắt 4,9 tỷ m³), Tuyên Quang (1,0 tỷ m³) và Thác Bà (0,45 tỷ m³) thực hiện tích nước trước đỉnh lũ theo quy trình liên hồ. Nhờ đó, lưu lượng lũ lớn nhất qua trạm Thủy văn Hà Nội được khống chế ở mức 20.578 m³/s, giữ mực nước lũ thiết kế không vượt quá 13,40 m đối với trận lũ chu kỳ 250 năm.

Vì sao bắt buộc phải di dời gần 12.000 hộ dân sinh sống tại các bãi sông?

Khảo sát thực tế ghi nhận 41.488 hộ dân lấn chiếm bãi sông Hồng gây cản trở dòng chảy nghiêm trọng. Phương án tối ưu yêu cầu giải tỏa 11.954 hộ dân nằm trong chỉ giới thoát lũ nguy hiểm (10.548 hộ sông Hồng, 1.406 hộ sông Đuống) nhằm khôi phục mặt cắt ướt tự nhiên, đảm bảo an toàn tính mạng cho người dân khi các hồ thượng nguồn xả lũ khẩn cấp.

Ưu điểm nổi trội của việc kết hợp mô hình MIKE 11 và MIKE 21 HD trong nghiên cứu này là gì?

Mô hình 1D MIKE 11 cho phép tính toán nhanh mạng lưới thủy lực rộng lớn gồm 12 sông chính, trong khi mô hình 2D MIKE 21 HD với lưới ô vuông 50 x 50 m mô phỏng chính xác trường dòng chảy phức tạp, vận tốc tràn bãi và tương tác lòng dẫn - đê bối tại vùng trọng điểm Hà Nội.

Kết luận

  • Luận văn đã xác lập đầy đủ cơ sở khoa học và thực tiễn về chế độ thủy lực dòng chảy lũ trên hệ thống sông Hồng và sông Đuống đoạn qua thành phố Hà Nội.
  • Ứng dụng thành công liên kết mô hình thủy văn NAM, mô hình thủy lực 1D MIKE 11 và 2D MIKE 21 HD trên lưới phân giải 50 x 50 m để mô phỏng chính xác các kịch bản lũ cực đoan.
  • Xác định bộ thông số mực nước và lưu lượng lũ thiết kế chuẩn xác tại các vị trí xung yếu, tiêu biểu là 13,40 m (20.578 m³/s) tại trạm Hà Nội và 14,65 m (27.538 m³/s) tại Thượng Cát.
  • Đề xuất giải pháp toàn diện kết hợp cắm mốc chỉ giới thoát lũ 100 m/mốc, nạo vét bãi tràn, kiên cố hóa 37,709 km đê cấp đặc biệt và di dời 11.954 hộ dân lấn chiếm lòng sông.
  • Luận văn đóng góp khung phương pháp luận hiện đại phục vụ công tác quy hoạch đê điều, quản lý bãi sông và phòng chống thiên tai cho lưu vực sông Hồng - sông Thái Bình.

Lộ trình triển khai khuyến nghị chia làm ba giai đoạn: hoàn thành cắm mốc chỉ giới thoát lũ vào năm 2015; hoàn tất nạo vét lòng dẫn và nâng cấp hệ thống đê kè vào năm 2020; đạt mục tiêu phủ xanh rừng phòng hộ đầu nguồn trên 15% vào năm 2030. Các nhà quản lý đô thị, cơ quan thủy lợi và kỹ sư công trình cần nhanh chóng ứng dụng kết quả nghiên cứu này vào thực tiễn để bảo vệ vững chắc an toàn đê điều và phát triển bền vững Thủ đô Hà Nội.