Tổng quan nghiên cứu

Lưu vực sông Thu Bồn thuộc tỉnh Quảng Nam là một trong những khu vực chịu tác động thiên tai khốc liệt nhất tại miền Trung Việt Nam. Với tổng diện tích tự nhiên 10.406 km2, địa bàn tỉnh sở hữu nhiều tiềm năng phát triển kinh tế biển và văn hóa nhưng thường xuyên đối mặt với hiểm họa lũ lụt nghiêm trọng. Thống kê giai đoạn 2003 - 2007 cho thấy tổng thiệt hại kinh tế do thiên tai gây ra tại Quảng Nam chiếm trung bình khoảng 6,26% tổng GDP toàn tỉnh. Điển hình, cơn bão số 9 năm 2009 đã gây tổn thất tài sản ước tính hơn 500 tỷ đồng, cướp đi sinh mạng của nhiều người dân và tàn phá cơ sở hạ tầng ven sông.

Vấn đề nghiên cứu trọng tâm là sự hình thành các đợt lũ đặc biệt lớn và hiện tượng ngập úng kéo dài trên diện rộng tại vùng đồng bằng hạ lưu. Mục tiêu cụ thể của công trình nhằm xác định nguyên nhân cơ bản, phân tích các yếu tố tự nhiên lẫn nhân sinh tác động đến dòng chảy lũ, ứng dụng công nghệ mô hình hóa tiên tiến để xây dựng bản đồ phân vùng ngập lụt và dự báo diện ngập tương ứng với các kịch bản biến đổi khí hậu.

Phạm vi không gian nghiên cứu trải rộng trên toàn bộ lưu vực sông Thu Bồn thuộc địa phận tỉnh Quảng Nam, giới hạn từ 14°57'10" đến 16°03'50" vĩ độ Bắc và 107°12'50" đến 108°44'20" kinh độ Đông, tập trung sâu vào vùng đồng bằng hạ du như thị xã Hội An, huyện Điện Bàn, Đại Lộc và Duy Xuyên. Kết quả nghiên cứu có ý nghĩa khoa học và 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 quy hoạch đô thị, bảo vệ hơn 200 di tích lịch sử và hỗ trợ Chiến lược Quốc gia phòng chống và giảm nhẹ thiên tai với mục tiêu giảm thiểu trên 15% thiệt hại về người và tài sản trong khu 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 trên nền tảng lý thuyết tương tác hệ thống tự nhiên - nhân sinh và lý thuyết thủy động lực học dòng chảy trong mạng lưới sông ngòi tự nhiên. Ba khái niệm cốt lõi được chuẩn hóa bao gồm:

  1. Lũ lụt (Flood): Hiện tượng dòng chảy vượt quá khả năng thoát nước của lòng dẫn, dâng cao và nhấn chìm các vùng đất trũng thấp ven sông.
  2. Đồng bằng ngập lụt (Floodplain): Vùng đất bồi tích ven bờ chịu tác động ngập định kỳ theo các cấp tần suất lũ.
  3. Đỉnh lũ và tần suất lũ (Peak Flood & Return Period): Mực nước và lưu lượng cực đại quan trắc được trong trận lũ, xác định chu kỳ lặp lại thông qua phương trình tần suất P = 1/RI.

Mô hình nghiên cứu tích hợp phương trình Saint-Venant một chiều (1D) giải bài toán dòng chảy không ổn định kết hợp mô hình số hóa độ cao DEM (Digital Elevation Model) trên nền tảng hệ thống thông tin địa lý (GIS). Khung tiếp cận này cho phép liên kết chặt chẽ giữa quá trình tích tụ bồi tích lòng dẫn, đặc điểm biến dạng bờ sông với cơ chế thoát lũ vùng cửa sông ven biển.

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

Nghiên cứu kế thừa và xử lý chuỗi số liệu khí tượng thủy văn đo đạc liên tục hơn 32 năm (1976 - 2007) tại các trạm quan trắc tiêu biểu như Nông Sơn, Thạnh Mỹ, Ái Nghĩa, Câu Lâu và Giao Thủy. Cỡ mẫu nghiên cứu bao gồm toàn bộ số liệu đo đạc của 78 phụ lưu cấp I đến cấp IV với chiều dài dòng nhánh lớn hơn 10 km trên tổng diện tích hứng nước 10.350 km2. Phương pháp chọn mẫu tập trung vào các trận lũ lịch sử điển hình năm 1964, 1998, 1999 và 2007 để làm dữ liệu chuẩn kiểm định.

Phương pháp khảo sát thực địa được triển khai dọc theo 205 km chiều dài dòng chính và các phân lưu Cửa Đại, Cửa Hàn, Cửa Lở nhằm đo đạc trắc dọc, trắc ngang và điều tra cao trình vết lũ thực tế. Phương pháp viễn thám và GIS sử dụng ảnh vệ tinh SPOT độ phân giải cao kết hợp bản đồ chuyên đề tỷ lệ 1/25.000 để chiết xuất các lớp thông tin về độ dốc, địa mạo và hiện trạng sử dụng đất.

Lý do lựa chọn bộ công cụ mô hình thủy lực MIKE 11 (DHI) kết hợp MIKE 11 GIS là nhờ thuật toán sai phân hữu hạn 6 điểm Abbott - Ionescu ổn định, có khả năng mô phỏng chính xác mạng lưới phân lưu phức tạp, kiểm soát chế độ thủy triều ven bờ và liên kết trực tiếp với dữ liệu không gian DEM để xuất bản đồ ngập lụt đa kịch bản với độ tin cậy vượt trội so với các mô hình truyền thống trước đây.

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

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

Quá trình phân tích dữ liệu thực nghiệm và kết quả mô phỏng số đã mang lại các phát hiện cốt lõi sau:

Thứ nhất, lượng mưa cực đoan đóng vai trò kích hoạt chính. Tổng lượng mưa trung bình năm toàn lưu vực đạt 2.978 mm, tương ứng khối lượng khoảng 30,7 tỷ m3 nước. Mùa lũ kéo dài từ tháng 10 đến tháng 12 tập trung tới 65% - 80% tổng lượng nước cả năm, riêng hai tháng 10 và 11 đã chiếm 40% - 50%. Các tổ hợp thời tiết nguy hiểm gồm bão kết hợp không khí lạnh gây ra khoảng 73% các trận lũ lớn lịch sử. Tại tâm mưa Trà My, lượng mưa tháng 10/1981 từng đạt kỷ lục 1.894 mm, lượng mưa ngày ứng với tần suất 1% chạm ngưỡng 900 mm.

Thứ hai, đặc điểm sơn văn tạo hiệu ứng truyền lũ cực nhanh. Địa hình đồi núi chiếm hơn 80% diện tích lưu vực với độ dốc sườn trung bình 25,5%, nhiều nơi sườn dốc trên 35 độ. Hệ số uốn khúc của sông chính Thu Bồn đạt 1,86 và sông Bung đạt 2,02. Cấu trúc mạng sông hình nan quạt gom dòng chảy từ 3 nhánh lớn (sông Bung, dòng chính Thu Bồn và sông Tĩnh Yên) khiến lũ tập trung nhanh, cường suất mực nước dâng lên tới 1 - 2 mét/giờ, tạo đỉnh lũ đơn hoặc lũ đa đỉnh rất nhọn.

Thứ ba, sự tắc nghẽn tiêu thoát nước tại hạ du làm trầm trọng hóa mức độ ngập lụt. Đồng bằng hạ lưu rất nhỏ hẹp, chỉ chiếm khoảng 20% diện tích lưu vực. Lưu lượng đỉnh lũ lớn nhất đo được tại trạm Nông Sơn lên đến 10.815 m3/s trong trận lũ năm 2007 và đạt 18.250 m3/s trong trận lũ lịch sử năm 1964. Hệ thống cồn cát ven biển cao 6 - 10 m cùng các tuyến đường sắt Bắc - Nam, Quốc lộ 1A chạy vuông góc hướng thoát lũ đã tạo thành những con đê chắn nhân tạo. Thêm vào đó, doi cát bồi lấp phía ngoài Cửa Đại hình thành từ năm 1983 làm giảm mạnh tiết diện thoát nước ra biển.

Thảo luận kết quả

Kết quả mô phỏng qua mô hình MIKE 11 GIS cho thấy sự ăn khớp chặt chẽ giữa mực nước tính toán và mực nước thực đo tại trạm Câu Lâu với sai số biên dưới 5%, chứng minh độ tin cậy cao của phương pháp mô hình hóa. Dữ liệu ngập lụt được trực quan hóa sinh động thông qua biểu đồ đường cong quá trình thủy văn (hydrograph) lưu lượng - mực nước theo thời gian và bảng phân cấp diện tích ngập theo 4 dải độ sâu: dưới 0,5 m; 0,5 - 1,0 m; 1,0 - 2,0 m và trên 2,0 m.

So sánh với các nghiên cứu trước đây như mô hình DHM hay bản đồ ngập lụt xây dựng trên nền địa hình tỷ lệ 1/500.000 của những năm 1996 - 2000, nghiên cứu này đã khắc phục hoàn toàn tính chắp vá dữ liệu. Việc tích hợp mô hình độ cao số DEM chi tiết tỷ lệ 1/25.000 giúp bóc tách rõ ranh giới ngập của từng xã, thôn thuộc huyện Điện Bàn, Duy Xuyên và thành phố Hội An.

Nguyên nhân ngập sâu không chỉ do lượng nước từ thượng nguồn đổ về đột ngột mà còn do hiện tượng nâng cao đáy sông tại các đoạn uốn khúc cổ và tác động chặn dòng của vòm nâng hiện đại Hòa Tiến, khiến phân lưu Vu Gia chuyển hướng dòng chảy mạnh về phía Cửa Đại thay vì đổ về Cửa Hàn như trước năm 1993. Khi áp dụng các kịch bản biến đổi khí hậu với lượng mưa tăng lần lượt 0,7%, 1,5% và 2,5%, diện tích ngập lụt trên mức báo động III mở rộng thêm từ 8,5% đến 18,2%, đặt ra thách thức sống còn cho hạ tầng ven biển.

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

Dựa trên kết quả mô phỏng và phân tích nguyên nhân tai biến, 4 nhóm giải pháp chiến lược được đề xuất triển khai cụ thể:

  1. Nạo vét và chỉnh trị luồng lạch cửa sông: Tiến hành nạo vét khẩn cấp khoảng 15 - 20 km luồng dẫn tại Cửa Đại, đoạn sông Câu Lâu và hệ thống sông Cổ Cò; mục tiêu gia tăng lưu lượng thoát lũ thêm 25% - 30% và hạ thấp mực nước ngập cục bộ tại Hội An từ 0,3 - 0,5 m; timeline thực hiện trong giai đoạn 2013 - 2016; chủ thể thực hiện là Sở Nông nghiệp và Phát triển nông thôn phối hợp Ban Quản lý Dự án Thủy lợi tỉnh Quảng Nam.

  2. Cải tạo nâng cấp khẩu độ thoát lũ công trình giao thông: Xây dựng bổ sung cầu cạn và mở rộng khẩu độ thoát nước tại các điểm xung yếu dọc tuyến Quốc lộ 1A và đường sắt Thống Nhất; target metric nâng khả năng tiêu thoát lũ đồng bằng thêm 20%; timeline hoàn thành trước năm 2018; chủ thể thực hiện do Sở Giao thông Vận tải Quảng Nam phối hợp với Tổng công ty Đường sắt Việt Nam và Bộ Giao thông Vận tải triển khai.

  3. Xây dựng hệ thống cảnh báo sớm thời gian thực: Thiết lập hệ thống hỗ trợ ra quyết định ứng dụng mô hình MIKE FLOOD WATCH tích hợp trạm đo mưa tự động; mục tiêu kéo dài thời gian dự báo trước khi lũ về lên 12 - 24 giờ với độ chính xác đạt trên 85%; timeline vận hành giai đoạn 2013 - 2015; chủ thể thực hiện là Ban Chỉ huy Phòng chống thiên tai và Tìm kiếm cứu nạn tỉnh Quảng Nam.

  4. Tăng cường quản lý rừng đầu nguồn và quy hoạch hành lang thoát lũ: Tái lập quy hoạch phân vùng sử dụng đất ven sông, nghiêm cấm xây dựng công trình cản lũ trái phép, đồng thời nâng tỷ lệ che phủ rừng phòng hộ đầu nguồn sông Tranh và sông Bung lên trên 65%; timeline thực hiện liên tục từ 2013 đến 2020; chủ thể thực hiện do Sở Tài nguyên và Môi trường cùng Chi cục Kiểm lâm tỉnh chịu trách nhiệm giám sát.

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

Công trình nghiên cứu mang giá trị ứng dụng thực tiễn cao, đặc biệt hữu ích cho 4 nhóm đối tượng sau:

  1. Cơ quan quản lý nhà nước và Ban Chỉ huy Phòng chống thiên tai các cấp: Cung cấp bản đồ nguy cơ ngập lụt chi tiết để hoạch định phương án sơ tán dân cư kịp thời, chuẩn bị phương tiện cứu hộ tại các vùng cô lập và phân bổ ngân sách phòng chống lụt bão hàng năm đạt hiệu quả tối ưu.

  2. Kỹ sư và nhà quy hoạch đô thị, giao thông, thủy lợi: Sử dụng các thông số thủy lực đỉnh lũ và độ sâu ngập làm căn cứ tính toán khẩu độ cầu cống, xác định cốt nền xây dựng an toàn cho các khu đô thị mới tại Điện Bàn, Hội An và khu kinh tế mở Chu Lai.

  3. Giảng viên, nghiên cứu sinh và học viên chuyên ngành Môi trường, Thủy văn, GIS: Làm tài liệu tham khảo mẫu mực về quy trình tích hợp mô hình thủy động lực học MIKE 11 với GIS trong phân tích tai biến tự nhiên, cung cấp nguồn số liệu thực đo phong phú phục vụ các đề tài nghiên cứu mở rộng.

  4. Các tổ chức bảo tồn di sản văn hóa và phát triển du lịch: Ứng dụng bản đồ cảnh báo mức độ ngập để xây dựng phương án kè chắn, chống xói lở móng di tích cổ tại phố cổ Hội An và khu di tích Mỹ Sơn, đảm bảo an toàn cho du khách và bảo vệ hơn 1.000 công trình kiến trúc gỗ có giá trị lịch sử.

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

Đâu là nguyên nhân cốt lõi gây ra các trận lũ lịch sử tại hạ lưu sông Thu Bồn?
Nguyên nhân cốt lõi là sự kết hợp đồng thời giữa các hình thế thời tiết cực đoan (bão, áp thấp nhiệt đới tương tác với không khí lạnh chiếm 73% các đợt lũ lớn) và đặc điểm địa hình dốc trên 25,5% ở thượng lưu. Địa hình gom nước dạng nan quạt làm lưu lượng đỉnh lũ tại Nông Sơn có thể vượt 10.815 m3/s, dồn nhanh xuống vùng hạ lưu nhỏ hẹp.

Tại sao mô hình MIKE 11 kết hợp GIS lại được ưu tiên ứng dụng trong nghiên cứu này?
Mô hình MIKE 11 sở hữu thuật toán giải hệ phương trình Saint-Venant 1D hoàn chỉnh, mô phỏng chuẩn xác mạng sông phân lưu phức tạp và tương tác triều. Khi tích hợp với GIS trên bản đồ tỷ lệ 1/25.000, mô hình giúp tự động hóa việc xuất bản đồ phân vùng ngập lụt trực quan, đạt sai số kiểm định dưới 5% so với mực nước thực đo tại trạm Câu Lâu.

Hạ tầng giao thông như Quốc lộ 1A và đường sắt ảnh hưởng thế nào đến ngập lụt?
Các tuyến đường sắt và Quốc lộ 1A chạy dọc đồng bằng vuông góc với hướng thoát lũ tự nhiên, kết hợp hệ thống cồn cát cao 6 - 10 m đóng vai trò như đê chắn. Khẩu độ cống cầu hiện hữu không đủ tiêu thoát lượng nước lũ khổng lồ, làm tăng thời gian ngập úng nội đồng thêm 1 - 2 ngày và nâng cao mực nước ngập từ 0,3 - 0,7 m.

Kịch bản biến đổi khí hậu ảnh hưởng ra sao đến diện ngập lụt vùng nghiên cứu?
Theo kết quả mô phỏng kịch bản biến đổi khí hậu trong luận văn, khi lượng mưa tăng từ 0,7% đến 2,5%, diện tích ngập sâu trên 1 m tại hạ lưu sông Thu Bồn mở rộng thêm khoảng 8,5% đến 18,2%. Vùng chịu tác động mạnh nhất tập trung tại các khu vực trũng thấp ven sông thuộc huyện Duy Xuyên, Điện Bàn và thành phố Hội An.

Đóng góp mới nổi bật nhất của luận văn so với các nghiên cứu trước đây là gì?
Luận văn đã khắc phục hạn chế của các nghiên cứu giai đoạn 1996 - 2000 vốn dùng bản đồ tỷ lệ nhỏ 1/500.000. Tác giả đã liên kết dữ liệu quan trắc 32 năm với mô hình số độ cao DEM chi tiết, định lượng hóa thành công diện ngập cho các kịch bản lũ năm 1999, 2007 và dự báo biến đổi khí hậu phục vụ quy hoạch thực tiễn.

Kết luận

  • Luận văn đã phân tích toàn diện mối quan hệ tương tác giữa các yếu tố tự nhiên và nhân sinh tác động đến nguy cơ tai biến lũ lụt tại hạ lưu sông Thu Bồn, Quảng Nam.
  • Khẳng định tổ hợp bão, áp thấp nhiệt đới và không khí lạnh tạo ra lượng mưa cực đại trên 900 mm/ngày là nguyên nhân trực tiếp kích hoạt lũ lụt nghiêm trọng.
  • Ứng dụng thành công mô hình thủy lực MIKE 11 kết hợp GIS trên chuỗi số liệu quan trắc hơn 32 năm để xây dựng bản đồ phân vùng ngập lụt chuẩn xác với sai số dưới 5%.
  • Dự báo chi tiết diện tích ngập lụt ứng với các kịch bản biến đổi khí hậu, chỉ ra nguy cơ mở rộng vùng ngập thêm 8,5% - 18,2% khi lượng mưa gia tăng.
  • Đề xuất hệ thống 4 nhóm giải pháp công trình và phi công trình đồng bộ, xác định rõ lộ trình thực hiện đến năm 2020 và đơn vị chịu trách nhiệm cụ thể.

Đóng góp chính của luận văn là thiết lập cơ sở khoa học định lượng vững chắc cho công tác quản lý rủi ro thiên tai và quy hoạch không gian kinh tế - xã hội bền vững tại tỉnh Quảng Nam. Bước tiếp theo cần triển khai ngay trong giai đoạn tới là số hóa hệ thống bản đồ ngập lụt đưa lên nền tảng WebGIS để chính quyền và hơn 1,4 triệu người dân địa phương chủ động tra cứu, ứng phó thiên tai. Các cơ quan quản lý và đơn vị nghiên cứu hãy tham khảo, áp dụng ngay mô hình này vào thực tiễn quản lý lưu vực sông nhằm bảo vệ an toàn tính mạng và di sản văn hóa trước biến đổi khí hậu.