Tổng quan nghiên cứu

Lưu vực sông Nhuệ - Đáy có tổng diện tích tự nhiên đạt 7.665 km2, trong đó phần diện tích thuộc địa bàn thành phố Hà Nội chiếm 2.361 km2 với quy mô dân số trên 5,7 triệu người. Đây là trung tâm chính trị, văn hóa và kinh tế đặc biệt quan trọng của cả nước, đóng góp khoảng 280 nghìn tỷ đồng vào tổng sản phẩm nội địa (GRDP) năm 2011. Mật độ dân số tại các quận nội thành ở mức rất cao, tiêu biểu như quận Đống Đa lên tới 38.936 người/km2, tạo sức ép to lớn lên hệ thống hạ tầng đô thị. Vấn đề nghiên cứu trọng tâm xuất phát từ thực trạng ngập úng nội đồng diễn biến ngày càng phức tạp khi xuất hiện mưa lớn cực đoan. Khác với lũ quét hay ngập do tràn bờ ở miền Trung, khu vực Hà Nội được bao bọc bởi hệ thống đê sông kiên cố, vừa ngăn nước lũ sông Hồng nhưng đồng thời lại gây khó khăn lớn cho việc tiêu thoát úng nội đồng. Trong đợt ngập lịch sử năm 2008, lượng mưa năm tại trạm Hà Đông lên tới 2.977,9 mm đã khiến nhiều khu vực chìm sâu trong nước, làm ngưng trệ các hoạt động dân sinh.

Mục tiêu cụ thể của luận văn là ứng dụng bộ mô hình thủy văn - thủy lực MIKE FLOOD kết hợp công nghệ GIS để mô phỏng diễn biến ngập lụt theo các tần suất thiết kế 1%, 2%, 5%, 10% và thành lập bản đồ tính dễ bị tổn thương cho lưu vực sông Nhuệ - Đáy trên địa bàn Hà Nội. Phạm vi nghiên cứu sử dụng chuỗi dữ liệu khí tượng thủy văn quan trắc từ năm 1961 đến năm 2011, tập trung phân tích các đợt mưa lũ điển hình năm 2007 và 2008. Kết quả nghiên cứu mang ý nghĩa thực tiễn sâu sắc, cung cấp cơ sở dữ liệu định lượng phục vụ công tác quản lý rủi ro thiên tai, bảo vệ các mục tiêu phát triển kinh tế và hỗ trợ triển khai hiệu quả Quy hoạch chung xây dựng Thủ đô Hà Nội đến năm 2030 theo Quyết định số 1878/QĐ-TTg của Thủ tướng Chính phủ.

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 vận dụng khung lý thuyết đánh giá tính dễ bị tổn thương do thiên tai từ Ủy ban Liên chính phủ về Biến đổi Khí hậu (IPCC), Chiến lược Giảm nhẹ Thiên tai Quốc tế (ISDR 2004) và Viện Giáo dục Nước Quốc tế (UNESCO-IHE). Khái niệm tổn thương được chuẩn hóa thông qua mô hình hàm số: Tổn thương bằng Sự lộ diện trừ đi Khả năng chống chịu, kế thừa từ các công trình nghiên cứu của Villagrán de León (2006) và Đặng Đình Khá (2011). Năm khái niệm học thuật cốt lõi được áp dụng gồm: Hiểm họa ngập lụt (Hazard), Sự lộ diện của đối tượng (Exposure), Tính nhạy cảm (Sensitivity), Khả năng chống chịu của cộng đồng (Coping Capacity) và Bản đồ tính dễ bị tổn thương do lũ (Flood Vulnerability Map).

Về mặt toán học thủy động lực, mô hình mưa - dòng chảy NAM (Nedbør-Afstrømnings-Model) với cấu trúc 19 thông số được sử dụng để mô phỏng dòng chảy mặt, dòng sát mặt và dòng ngầm từ các bể chứa lưu vực. Dòng chảy trong mạng lưới sông một chiều được mô tả qua hệ phương trình Saint-Venant giải bằng sơ đồ sai phân ẩn 6 điểm Abbott trong mô-đun MIKE 11. Sự lan truyền dòng tràn trên mặt đất hai chiều được mô tả bằng hệ phương trình nước nông giải bằng phương pháp sai phân hữu hạn xen hướng ẩn ADI trong mô-đun MIKE 21.

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

Nguồn dữ liệu địa hình phục vụ tính toán thủy lực được số hóa chi tiết từ 122 mặt cắt ngang sông đo đạc thực tế, gồm 62 mặt cắt sông Đáy, 26 mặt cắt sông Nhuệ, 20 mặt cắt sông Tích và 14 mặt cắt sông Hoàng Long. Dữ liệu khí tượng đầu vào bao gồm chuỗi số liệu mưa giờ và bốc hơi giai đoạn 1961 - 2010 tại trạm Láng, Ba Thá, Hưng Thi. Dữ liệu hiện trạng sử dụng đất và chỉ tiêu kinh tế - xã hội được trích xuất từ niên giám thống kê thành phố Hà Nội.

Nhằm lượng hóa khả năng chống chịu của người dân, nghiên cứu thực hiện điều tra xã hội học với cỡ mẫu 300 phiếu khảo sát. Phương pháp chọn mẫu phân tầng ngẫu nhiên (Stratified Random Sampling) được triển khai tại các khu vực dân cư trũng thấp ven sông thuộc huyện Chương Mỹ, Ứng Hòa, Quốc Oai, Ba Vì và quận Hà Đông để thu thập thông tin về nhận thức, công trình phòng hộ gia đình và khả năng phục hồi sau ngập úng.

Lý do lựa chọn giải pháp tích hợp MIKE FLOOD là nhằm khắc phục triệt để nhược điểm của mô hình 1D (không mô tả được trường vận tốc lan tỏa trên bề mặt đồng ruộng) và mô hình 2D độc lập (đòi hỏi chia lưới quá mịn gây quá tải thời gian tính toán). Bốn dạng liên kết không gian (liên kết tiêu chuẩn, liên kết bên, liên kết công trình và liên kết khô) cho phép mô phỏng chính xác quá trình tràn qua bờ đê bối. Mô hình được hiệu chỉnh với trận lũ tháng 8/1973, tháng 10/2007 và kiểm định với trận mưa lũ lịch sử tháng 11/2008.

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 thủy văn - thủy lực đạt độ chính xác rất cao trong quá trình mô phỏng. Chỉ số hiệu quả Nash-Sutcliffe trong giai đoạn hiệu chỉnh mô hình mưa - dòng chảy NAM đạt 71,1% tại trạm Ba Thá và 89,7% tại trạm Hưng Thi; trong giai đoạn kiểm định độc lập, chỉ số Nash đạt 91,0% tại trạm Ba Thá và 81,0% tại trạm Hưng Thi. Quá trình kiểm định trận ngập tháng 11/2008 cho thấy mực nước đỉnh lũ tính toán sai lệch dưới 0,15 m so với thực đo tại trạm Phủ Lý, Ba Thá và Hà Đông, mức độ trùng khớp diện tích ngập lụt đạt trên 88% so với vết lũ thực tế.

Thứ hai, bản đồ hiểm họa thể hiện rõ sự phân bố ngập lụt theo các cấp tần suất. Ứng với trận mưa thiết kế tần suất 1% (chu kỳ lặp lại 100 năm), khu vực ngập sâu trên 1,0 m tập trung lớn tại vùng trũng phía Tây Nam như huyện Ứng Hòa, Mỹ Đức, Chương Mỹ và dọc hành lang sông Tích. Thời gian ngập kéo dài từ 3 đến hơn 7 ngày do đáy sông Duy Tiên, sông La Khê bị bồi lắng cao hơn thiết kế từ 2,0 m đến 3,0 m, làm giảm khả năng tự tiêu thoát nước ra sông Đáy.

Thứ ba, bản đồ tính dễ bị tổn thương tổng hợp được thành lập dựa trên việc chồng xếp bản đồ hiểm họa (trọng số độ sâu ngập 50%, thời gian ngập 30%, vận tốc dòng chảy 20%), bản đồ lộ diện sử dụng đất (khu dân cư đô thị có trọng số nhạy cảm cao nhất 0,85 so với đất nông nghiệp 0,45) và bản đồ khả năng chống chịu. Kết quả chỉ ra rằng diện tích chịu tổn thương ở mức rất cao và cao chiếm khoảng 28,5% diện tích nghiên cứu, tập trung tại các vùng đô thị hóa nhanh và khu dân cư ngoài đê bối.

Thảo luận kết quả

Nguyên nhân dẫn đến mức độ tổn thương ngập lụt cao tại lưu vực sông Nhuệ - Đáy là sự kết hợp giữa yếu tố tự nhiên và phát triển đô thị. Địa hình khu vực nghiêng dần từ Tây Bắc xuống Đông Nam với độ cao phổ biến từ +5 m đến +20 m, tạo thành các thung lũng tụ nước. Cùng với đó, quá trình bê tông hóa làm thu hẹp không gian trữ nước của hơn 80 hồ đầm tự nhiên (như hồ Yên Sở 145,6 ha, hồ Linh Đàm 67,0 ha), làm giảm đáng kể khả năng tự điều tiết khi mưa lớn.

So với các nghiên cứu trước đây chỉ đánh giá rủi ro dựa trên độ sâu ngập đơn thuần, cách tiếp cận tích hợp khả năng chống chịu giúp lượng hóa chính xác mức độ tổn hại thực tế. Các cộng đồng ven sông có kinh nghiệm ứng phó tốt giúp giảm chỉ số tổn thương cục bộ, trong khi các khu đô thị mới thiếu kinh nghiệm đối mặt với rủi ro cao hơn. Kết quả dữ liệu được trực quan hóa hiệu quả thông qua hệ thống bản đồ GIS đa lớp, kết hợp biểu đồ cột so sánh diện tích ngập theo 4 cấp tần suất và bảng ma trận chéo thể hiện tương quan giữa độ sâu ngập và từng loại hình sử dụng đất.

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

Nâng cấp hạ tầng trạm bơm và hoàn thiện hệ thống tiêu úng: Tiến hành đầu tư nâng tổng công suất tiêu thoát nước toàn hệ thống thêm 35% đến 40%, ưu tiên mở rộng cụm công trình đầu mối trạm bơm Yên Sở (công suất 90 m3/s), trạm bơm Yên Nghĩa (công suất 120 m3/s) và cải tạo cống Liên Mạc. Mục tiêu cụ thể là kiểm soát thời gian ngập úng đô thị dưới 12 giờ sau các trận mưa có cường độ trên 200 mm/24h. Timeline thực hiện: Giai đoạn 2013 - 2020. Chủ thể thực hiện: Sở Xây dựng phối hợp Sở Nông nghiệp và Phát triển nông thôn Hà Nội.

Nạo vét toàn diện và khơi thông dòng chảy các trục sông chính: Triển khai nạo vét lòng dẫn trên toàn tuyến 74 km sông Nhuệ, 250 km sông Đáy cùng các sông nhánh gồm sông La Khê (6,8 km), sông Vân Đình (11,8 km) và sông Duy Tiên (21 km), hạ cao độ đáy sông về đúng cao trình thiết kế để tăng 25% lưu lượng thoát lũ tự nhiên. Kiên quyết giải tỏa các điểm lấn chiếm dòng chảy và rác thải hai bên bờ sông. Timeline thực hiện: Giai đoạn 2013 - 2018. Chủ thể thực hiện: Ban Quản lý Dự án Đầu tư Xây dựng Công trình Nông nghiệp và PTNT Hà Nội.

Bảo vệ nghiêm ngặt không gian trữ nước mặt và kiểm soát quy hoạch: Ban hành quy định bảo tồn nguyên trạng diện tích mặt nước của hơn 80 hồ đầm trên địa bàn, đặc biệt là các hồ điều hòa lớn như hồ Suối Hai (1.000 ha), hồ Đồng Sương (203,1 ha) và hồ Xuân Khanh (104 ha). Yêu cầu các dự án khu đô thị mới bắt buộc dành tối thiểu 20% quỹ đất cho diện tích thấm nước tự nhiên và hồ cảnh quan. Timeline thực hiện: Hàng năm trong giai đoạn 2013 - 2025. Chủ thể thực hiện: Sở Tài nguyên và Môi trường cùng Sở Quy hoạch - Kiến trúc Hà Nội.

Tăng cường hệ thống cảnh báo sớm và nâng cao năng lực cộng đồng: Lắp đặt hệ thống giám sát mực nước tự động tại 15 vị trí khí tượng thủy văn then chốt, tích hợp bản đồ tổn thương vào hệ thống hỗ trợ ra quyết định phòng chống thiên tai nhằm phát bản tin cảnh báo trước ngập lụt ít nhất 6 giờ. Tổ chức tập huấn kỹ năng ứng phó cho 100% cán bộ và người dân tại các xã thuộc vùng trũng thấp có nguy cơ tổn thương cao. Timeline thực hiện: Giai đoạn 2013 - 2017. Chủ thể thực hiện: Ban Chỉ huy Phòng chống Thiên tai và Tìm kiếm Cứu nạn TP Hà Nội.

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

Cơ quan quản lý nhà nước và nhà hoạch định chính sách: Ban Chỉ huy Phòng chống Thiên tai, Sở Nông nghiệp & PTNT, Sở Xây dựng và UBND các cấp sử dụng bản đồ tổn thương để xác định các trọng điểm ngập lụt, xây dựng kịch bản di dời dân cư và phân bổ vốn đầu tư công trình phòng hộ trong quy hoạch tổng thể phát triển Thủ đô.

Chuyên gia tư vấn quy hoạch đô thị và kỹ sư công trình thủy lợi: Các viện quy hoạch xây dựng và đơn vị thiết kế hạ tầng khai thác hệ thống 122 mặt cắt sông cùng chuỗi số liệu thủy lực tần suất 1% đến 10% để tính toán cao trình san nền, khẩu độ cầu cống và bố trí hành lang thoát lũ an toàn cho các khu đô thị mới.

Giảng viên, nghiên cứu sinh 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, Quản lý tài nguyên nước, Biến đổi khí hậu và Hệ thống thông tin địa lý (GIS) tham khảo quy trình liên kết mô hình NAM, MIKE 11, MIKE 21 và phương pháp luận thành lập bản đồ tổn thương làm tài liệu giảng dạy và nghiên cứu chuyên sâu.

Các tổ chức xã hội và cộng đồng dân cư địa phương: Các tổ chức phi chính phủ và người dân tại vùng trũng thuộc lưu vực 2.361 km2 khai thác thông tin rủi ro để chủ động xây dựng kế hoạch phòng tránh thiên tai tại chỗ, cải tạo nhà ở thích ứng và nâng cao năng lực tự phục hồi kinh tế sau ngập lụt.

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

Tại sao lưu vực sông Nhuệ - Đáy trên địa bàn Hà Nội thường xuyên bị ngập úng kéo dài khi mưa lớn?

Khu vực nghiên cứu có địa hình trũng thấp dạng lòng chảo với cao độ trung bình từ +5 m đến +20 m, chịu lượng mưa năm lớn đạt trung bình 1.608,19 mm tập trung chủ yếu vào mùa hè. Hệ thống đê kiên cố ngăn lũ sông Hồng nhưng cản trở tiêu thoát nội đồng, trong khi các trục sông chính như sông Nhuệ (74 km) và sông Đáy (250 km) bị bồi lắng nặng làm suy giảm khả năng thoát nước tự nhiên.

Ưu điểm vượt trội của mô hình MIKE FLOOD trong nghiên cứu ngập lụt là gì?

MIKE FLOOD tích hợp đồng bộ mô hình 1D MIKE 11 (tính toán 122 mặt cắt lòng dẫn) và mô hình 2D MIKE 21 (mô phỏng dòng tràn trên diện tích 2.361 km2). Nhờ 4 loại liên kết thủy lực linh hoạt, mô hình tái hiện chính xác trường vận tốc và độ sâu ngập trên bề mặt mà vẫn tối ưu hóa tài nguyên tính toán của máy tính.

Tính dễ bị tổn thương do lũ lụt được lượng hóa bằng phương pháp nào trong luận văn?

Tổn thương được xác định thông qua hàm số: Tổn thương bằng Sự lộ diện trừ đi Khả năng chống chịu. Trong đó, sự lộ diện được thiết lập từ bản đồ hiểm họa ngập lụt (tích hợp độ sâu, thời gian ngập và vận tốc dòng chảy) kết hợp với bản đồ sử dụng đất, còn khả năng chống chịu được định lượng qua 300 phiếu điều tra xã hội học thực tế tại cộng đồng.

Độ tin cậy của mô hình tính toán được chứng minh qua những số liệu nào?

Mô hình mưa - dòng chảy NAM đạt chỉ số hiệu quả Nash-Sutcliffe từ 81,0% đến 91,0% trong bước kiểm nghiệm. Khi mô phỏng kiểm tra trận lũ lịch sử tháng 11/2008, đường quá trình mực nước tính toán bám sát số liệu thực đo tại các trạm Phủ Lý, Ba Thá, Hà Đông với sai số đỉnh lũ dưới 0,15 m và độ trùng khớp diện tích ngập trên 88%.

Kết quả nghiên cứu hỗ trợ như thế nào cho công tác quy hoạch phát triển Thủ đô Hà Nội?

Luận văn cung cấp bộ bản đồ phân vùng rủi ro ngập lụt chi tiết ứng với các tần suất thiết kế 1%, 2%, 5%, 10%, phục vụ trực tiếp cho việc triển khai Quyết định số 1878/QĐ-TTg của Thủ tướng Chính phủ về Quy hoạch chung Thủ đô Hà Nội đến năm 2030, hỗ trợ định hướng không gian xây dựng an toàn và bảo vệ hơn 80 hồ đầm điều hòa.

Kết luận

Nội dung và giá trị khoa học của luận văn được tổng kết qua 5 điểm cốt lõi:

  • Xây dựng thành công hệ cơ sở dữ liệu thủy văn, khí tượng và địa hình với 122 mặt cắt ngang trên 4 tuyến sông chính thuộc lưu vực sông Nhuệ - Đáy.
  • Ứng dụng xuất sắc bộ công cụ mô hình hóa MIKE FLOOD (NAM, MIKE 11, MIKE 21) với chỉ số kiểm định Nash đạt từ 81,0% đến 91,0%, khẳng định độ chuẩn xác cao.
  • Thành lập bộ bản đồ hiểm họa ngập lụt ứng với 4 kịch bản tần suất thiết kế (1%, 2%, 5%, 10%) và bản đồ phân vùng tính dễ bị tổn thương chi tiết.
  • Lượng hóa rõ nét mối quan hệ giữa yếu tố hiểm họa ngập, tính nhạy cảm sử dụng đất và năng lực chống chịu thực tế của hơn 5,7 triệu cư dân đô thị.
  • Đề xuất đồng bộ 4 nhóm giải pháp công trình và phi công trình khả thi nhằm bảo vệ an toàn cho khu vực kinh tế tạo ra 280 nghìn tỷ đồng GRDP.

Luận văn là công trình khoa học có giá trị lý luận và ứng dụng thực tiễn cao trong lĩnh vực thủy văn và quản lý rủi ro thiên tai đô thị. Trong giai đoạn 2013 - 2025, các đơn vị quản lý cần tiếp tục cập nhật dữ liệu không gian theo các kịch bản biến đổi khí hậu mới. Các nhà quản lý, nhà quy hoạch và cộng đồng nghiên cứu hãy chủ động tham khảo và ứng dụng các phát hiện từ luận văn này nhằm nâng cao hiệu quả phòng chống ngập úng và hướng tới phát triển đô thị bền vững!