Tổng quan nghiên cứu
Lũ lụt tại khu vực miền Trung Việt Nam luôn là thách thức nghiêm trọng đối với sự phát triển bền vững do đặc thù sông ngòi ngắn, độ dốc lớn và tốc độ truyền lũ nhanh. Trong giai đoạn từ năm 1977 đến năm 2000, tại trạm Gia Vòng trên sông Bến Hải đã ghi nhận khoảng 110 trận lũ với tần suất trung bình 5,33 trận mỗi năm, gây thiệt hại nặng nề về người và tài sản cho 162.092 cư dân sinh sống trong khu vực. Trước bối cảnh biến đổi khí hậu diễn biến phức tạp, việc định lượng hóa mức độ rủi ro thiên tai trở thành yêu cầu cấp thiết nhằm thay thế các giải pháp ứng phó thụ động bằng chiến lược quản lý rủi ro chủ động.
Đề tài tập trung nghiên cứu, thử nghiệm và so sánh các công thức tính toán chỉ số đánh giá tính dễ bị tổn thương do lũ lụt (Flood Vulnerability Index - FVI) trên toàn bộ phạm vi không gian lưu vực sông Bến Hải thuộc tỉnh Quảng Trị, bao gồm 44 xã, phường, thị trấn trải rộng trên địa bàn 3 huyện Vĩnh Linh, Gio Linh và Cam Lộ. Mục tiêu cốt lõi của công trình là thiết lập bộ tiêu chí đánh giá chuẩn hóa, áp dụng thuật toán trọng số thống kê để tính toán chính xác mức độ tổn thương cho từng đơn vị hành chính cơ sở. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn khi cung cấp bộ công cụ định lượng hỗ trợ các nhà quản lý xây dựng phương án phòng chống lụt bão, quy hoạch không gian kinh tế - xã hội hợp lý và tối ưu hóa nguồn lực cứu trợ, góp phần giảm thiểu hơn 30% nguy cơ thiệt hại tiềm tàng cho cộng đồng đị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 vận dụng khung lý thuyết đánh giá tính dễ bị tổn thương của Ủy ban Liên chính phủ về Biến đổi Khí hậu (IPCC) và Tổ chức Giáo dục, Khoa học và Văn hóa Liên Hợp Quốc - Viện Giáo dục Tài nguyên Nước (UNESCO-IHE). Theo các lý thuyết này, tính dễ bị tổn thương do lũ lụt được xác lập qua hàm số của ba thành phần cơ bản: Độ phơi nhiễm (Exposure - E), Tính nhạy cảm (Sensitivity - S) và Khả năng phục hồi thích ứng (Resilience/Adaptive Capacity - R).
Khái niệm độ phơi nhiễm phản ánh mức độ tiếp xúc trực tiếp của con người, tài sản và cơ sở hạ tầng trước hiểm họa thiên tai; tính nhạy cảm thể hiện các đặc trưng kinh tế - xã hội làm gia tăng mức độ tổn thất; trong khi khả năng phục hồi biểu thị năng lực ứng phó, khắc phục và thích nghi của cộng đồng. Trên nền tảng đó, nghiên cứu tích hợp mô hình đánh giá của Connor và Hiroki (năm 2005) với phương trình FVI = C + H + S - M (bao gồm các yếu tố khí hậu, thủy văn địa hình, kinh tế xã hội và giải pháp ứng phó) cùng mô hình đa chiều cải tiến của Balica (giai đoạn 2007-2012) gồm 4 phân hệ chính: xã hội, kinh tế, môi trường và vật lý theo cấu trúc tương tác FVI = (E x S) / R.
Phương pháp nghiên cứu
Nghiên cứu tiếp cận cỡ mẫu toàn diện gồm 44 đơn vị hành chính cấp xã thuộc 3 huyện Vĩnh Linh, Gio Linh và Cam Lộ trên lưu vực sông Bến Hải (N = 44). Phương pháp chọn mẫu tổng thể được áp dụng nhằm đảm bảo tính bao quát, không bỏ sót bất kỳ tiểu vùng sinh thái hay vùng trũng ngập cục bộ nào dọc theo dòng chảy chính. Nguồn dữ liệu thu thập bao gồm chuỗi số liệu quan trắc khí tượng thủy văn hơn 40 năm tại trạm Gia Vòng, số liệu niên giám thống kê năm 2016 của Cục Thống kê tỉnh Quảng Trị, bản đồ địa hình và hiện trạng sử dụng đất.
Để xử lý các biến số có đơn vị đo lường khác nhau, phương pháp chuẩn hóa dữ liệu theo thang điểm từ 0 đến 1 của Balica cùng chuẩn hóa Connor và Hiroki được áp dụng. Nghiên cứu lựa chọn thuật toán phân tích trọng số Iyengar - Sudarshan kết hợp phương pháp phân tích thứ bậc AHP. Lý do lựa chọn tổ hợp phương pháp này là nhằm loại bỏ hoàn toàn yếu tố định kiến chủ quan của chuyên gia, bảo đảm trọng số của từng tiêu chí phản ánh chính xác mức độ biến thiên khách quan của dữ liệu thống kê thực tế. Toàn bộ kết quả tính toán sau đó được phân cấp mức độ tổn thương dựa trên hàm phân bố xác suất Beta và số hóa không gian trên nền tảng phần mềm MapInfo.
Kết quả nghiên cứu và thảo luận
Những phát hiện chính
Thứ nhất, nghiên cứu chỉ ra hiện trạng lớp phủ thực vật trên lưu vực sông Bến Hải bị suy thoái nghiêm trọng, ảnh hưởng trực tiếp đến khả năng điều tiết dòng chảy mặt. Diện tích trảng cây bụi chiếm tới 68,38% (tương đương 194 km²), trong khi diện tích rừng tự nhiên giàu chỉ chiếm vỏn vẹn 2,75% (7,802 km²), làm giảm khả năng giữ nước đầu nguồn và gia tăng hệ số dòng chảy lũ vượt mức 0,7.
Thứ hai, sự phân hóa mật độ dân số và mức độ phơi nhiễm giữa các tiểu vùng diễn ra vô cùng gay gắt. Mật độ dân số trung bình toàn lưu vực đạt 148 người/km², nhưng tại đô thị trung tâm như thị trấn Hồ Xá con số này lên tới 1.674 người/km², cao gấp hơn 23 lần so với xã miền núi Vĩnh Trường (72 người/km²), khiến chỉ số tổn thương xã hội tại vùng hạ lưu luôn ở mức báo động.
Thứ ba, quy luật lũ lụt có tính tập trung cao độ theo thời gian. Các trận lũ lớn xuất hiện dày đặc trong tháng 10 (chiếm 44,8% tổng số trận) và tháng 11 (chiếm 31,8% tổng số trận). Biên độ lũ dao động từ 2m đến 15m với cường suất lũ lên cực nhanh, đạt từ 100 đến 3.000 cm/giờ, tạo áp lực khổng lồ lên hệ thống đê điều và cơ sở hạ tầng giao thông.
Thứ tư, việc thử nghiệm mô hình Balica cho thấy độ nhạy vượt trội so với mô hình Connor và Hiroki, phân loại rõ rệt 44 xã thành 5 cấp độ tổn thương từ rất thấp đến rất cao, trong đó hơn 35% số xã vùng đồng bằng ven biển rơi vào nhóm có mức độ tổn thương nghiêm trọng nhất.
Thảo luận kết quả
Nguyên nhân chính dẫn đến sự chênh lệch chỉ số tổn thương giữa các khu vực xuất phát từ sự kết hợp bất lợi giữa địa hình đồi bát úp phía Tây dốc mạnh về phía Đông và vùng đồng bằng hạ lưu trũng thấp chỉ có cao trình từ 1 đến 2,5m. Khi xuất hiện các tổ hợp thời tiết cực đoan như bão kết hợp không khí lạnh, lượng mưa tăng vọt khiến thời gian truyền lũ chỉ mất từ 10 đến 24 giờ, trong khi ảnh hưởng của thủy triều khiến thời gian ngập kéo dài từ 2 đến 4 ngày.
Về mặt trực quan hóa, toàn bộ dữ liệu phân cấp tổn thương FVI được mô hình hóa rõ nét qua hệ thống bản đồ chuyên đề GIS, sử dụng thang màu tương phản từ xanh lục (tổn thương rất thấp) đến đỏ thẫm (tổn thương rất cao). Kết hợp với các biểu đồ cột thể hiện tương quan trọng số giữa 4 thành phần kinh tế, xã hội, môi trường và vật lý, các nhà quản lý có thể nhanh chóng nhận diện các túi ngập xung yếu.
So với các nghiên cứu trước đây tại lưu vực sông Vu Gia - Thu Bồn hay sông Thạch Hãn, việc ứng dụng công thức Balica kết hợp trọng số Iyengar - Sudarshan trong luận văn đã cải thiện độ chính xác phân vùng rủi ro thêm khoảng 15%, chứng minh tính ưu việt của việc tiếp cận đa tiêu chí định lượng trong quản lý thiên tai.
Đề xuất và khuyến nghị
Thứ nhất, nâng cấp và kiên cố hóa hệ thống công trình thoát lũ tại vùng hạ lưu. Ban Chỉ huy Phòng chống thiên tai và Tìm kiếm cứu nạn tỉnh Quảng Trị cần chủ trì dự án nạo vét lòng sông Bến Hải, gia cố 100% các tuyến đê bao xung yếu tại huyện Gio Linh và Cam Lộ, đặt mục tiêu giảm 25% thời gian ngập úng cục bộ trong giai đoạn 2020 - 2025.
Thứ hai, phục hồi và mở rộng diện tích rừng phòng hộ đầu nguồn. Sở Nông nghiệp và Phát triển nông thôn phối hợp với Ủy ban nhân dân các huyện miền núi triển khai chương trình trồng mới và làm giàu 50 km² rừng tự nhiên, phấn đấu nâng tỷ lệ che phủ rừng chất lượng cao từ 2,75% lên mức 12% trước năm 2030 nhằm tăng cường khả năng giữ nước và cắt lũ tự nhiên.
Thứ ba, hiện đại hóa mạng lưới quan trắc và cảnh báo sớm. Chi cục Thủy văn khu vực cần lắp đặt bổ sung 5 trạm đo mưa và mực nước tự động dọc lưu vực sông Bến Hải, ứng dụng công nghệ truyền tin thời gian thực để nâng cao độ chính xác dự báo và kéo dài thời gian cảnh báo trước lũ lên tối thiểu 6 giờ cho các khu dân cư trũng thấp.
Thứ tư, chuẩn hóa quy hoạch không gian sống và nâng cao năng lực tự thích ứng của cộng đồng. Chính quyền 44 xã thuộc lưu vực cần tổ chức tập huấn kỹ năng ứng phó thiên tai cho 100% hộ dân thuộc diện nguy cơ cao, đồng thời kiểm soát nghiêm ngặt quy hoạch xây dựng, nghiêm cấm san lấp các vùng thoát lũ tự nhiên có cao trình dưới 2m.
Đối tượng nên tham khảo luận văn
Thứ nhất, các nhà quản lý và cơ quan hoạch định chính sách phòng chống thiên tai tại tỉnh Quảng Trị. Nhóm đối tượng này có thể sử dụng trực tiếp bản đồ phân vùng FVI để phân bổ nguồn ngân sách cứu trợ, xây dựng kịch bản sơ tán dân cư sát với thực tế địa bàn và ban hành các quy định quản lý rủi ro thiên tai cấp cơ sở.
Thứ hai, các nhà khoa học, giảng viên và nghiên cứu sinh chuyên ngành Thủy văn học, Biến đổi khí hậu và Quản lý tài nguyên nước. Luận văn cung cấp khung phương pháp luận hoàn chỉnh về việc tích hợp thuật toán trọng số khách quan Iyengar - Sudarshan với mô hình FVI của Balica, đóng vai trò như một tài liệu tham khảo học thuật chuẩn mực cho các đề tài nghiên cứu lưu vực sông miền Trung.
Thứ ba, các kỹ sư quy hoạch đô thị và hạ tầng nông thôn. Dữ liệu chi tiết về cao trình ngập, biên độ lũ từ 2m đến 15m và thời gian duy trì ngập lụt là căn cứ kỹ thuật quan trọng để thiết kế nền móng công trình, hệ thống tiêu thoát nước và mạng lưới giao thông bền vững trước thiên tai.
Thứ tư, các tổ chức phi chính phủ và các định chế phát triển quốc tế hoạt động trong lĩnh vực ứng phó biến đổi khí hậu. Bản phân tích định lượng 44 xã là tài liệu thẩm định giá trị giúp xây dựng các dự án tài trợ sinh kế chống chịu bão lũ và nâng cao năng lực cộng đồng có địa chỉ rõ ràng.
Câu hỏi thường gặp
Lưu vực sông Bến Hải có những đặc trưng tự nhiên nào làm gia tăng rủi ro tổn thương do lũ lụt? Lưu vực có địa hình dốc đứng từ dãy Trường Sơn đổ ra biển, thảm phủ rừng giàu chỉ chiếm 2,75% trong khi đất cây bụi chiếm tới 68,38%. Lượng mưa tập trung cao độ trong các tháng 9, 10, 11 kết hợp với nền đất feralit tầng dày và thủy triều hạ lưu làm lũ dâng nhanh và rút rất chậm.
Mô hình Balica mang lại ưu thế gì nổi bật so với công thức Connor và Hiroki? Mô hình Balica tiếp cận toàn diện 4 thành phần kinh tế, xã hội, môi trường, vật lý qua 3 giai đoạn trước, trong và sau lũ. Mô hình này giúp lượng hóa chi tiết mức độ tổn thương của 44 xã trên lưu vực thay vì chỉ đưa ra bức tranh tổng quát như công thức tuyến tính 11 chỉ số của Connor và Hiroki.
Tại sao nghiên cứu lại áp dụng thuật toán tính trọng số Iyengar và Sudarshan? Thuật toán Iyengar và Sudarshan tính toán trọng số dựa hoàn toàn trên phương sai và độ biến thiên thống kê thực tế của dữ liệu thu thập. Phương pháp này triệt tiêu sự thiên lệch cảm tính của chuyên gia, đảm bảo kết quả phân cấp FVI cho từng xã đạt độ tin cậy và khách quan khoa học cao nhất.
Yếu tố mật độ dân số tác động như thế nào đến chỉ số tổn thương tại địa phương? Mật độ dân số cao làm gia tăng trực tiếp độ phơi nhiễm trước thiên tai. Điển hình tại thị trấn Hồ Xá với mật độ 1.674 người/km², rủi ro tổn thất về người và tài sản khi mực nước lũ vượt báo động 3 lớn hơn rất nhiều so với các khu vực miền núi thưa dân như xã Vĩnh Trường (72 người/km²).
Kết quả nghiên cứu của luận văn có thể chuyển giao và ứng dụng thực tế ra sao? Bản đồ phân vùng tổn thương FVI và bộ cơ sở dữ liệu số hóa có thể tích hợp trực tiếp vào hệ thống quản lý rủi ro thiên tai cấp tỉnh. Chính quyền địa phương căn cứ vào đó để lập bản đồ sơ tán chi tiết cho 162.092 người dân và xây dựng danh mục công trình phòng chống lũ ưu tiên.
Kết luận
- Luận văn đã hệ thống hóa và làm sáng tỏ khung lý thuyết đánh giá tính dễ bị tổn thương do lũ lụt dựa trên các mô hình tiên tiến của IPCC, UNESCO-IHE, Connor và Hiroki cùng mô hình đa thành phần của Balica.
- Thiết lập thành công bộ cơ sở dữ liệu đa chiều bao gồm các chỉ tiêu khí tượng thủy văn, địa hình, thảm phủ và kinh tế - xã hội cho toàn bộ 44 xã thuộc lưu vực sông Bến Hải.
- Ứng dụng xuất sắc thuật toán trọng số khách quan Iyengar - Sudarshan, chứng minh tính vượt trội của công thức Balica trong việc phản ánh chân thực độ nhạy cảm và khả năng chống chịu lũ lụt tại miền Trung.
- Xây dựng bản đồ phân cấp tính dễ bị tổn thương FVI trên nền tảng MapInfo, chỉ ra chính xác các trọng điểm nguy cơ cao tại khu vực đồng bằng trũng thấp ven sông và vùng đô thị hóa cao.
- Đề xuất hệ thống 4 nhóm giải pháp công trình và phi công trình đồng bộ, tạo luận cứ khoa học vững chắc phục vụ công tác quy hoạch phát triển bền vững lưu vực sông Bến Hải.
Trong giai đoạn 2020 - 2025, hướng nghiên cứu tiếp theo cần tích hợp thêm các kịch bản biến đổi khí hậu cập nhật và mô hình hóa thủy lực 2D để nâng cao hơn nữa độ phân giải không gian ngập lụt. Hãy liên hệ và khai thác toàn văn tài liệu luận văn tại Thư viện Đại học Quốc gia Hà Nội để nắm bắt trọn vẹn phương pháp luận định lượng rủi ro thiên tai chuyên sâu này.