Tổng quan về luận án
Thiên tai lũ lụt là một trong những hiểm họa tự nhiên cực đoan có tính toàn cầu, gây ra những thiệt hại nặng nề và phức tạp đối với hệ thống kinh tế – xã hội và môi trường sinh thái. Trong bối cảnh biến đổi khí hậu toàn cầu đang diễn biến khốc liệt kết hợp với áp lực nhân sinh từ các công trình thủy lợi, thủy điện khai thác thiếu đồng bộ, hiện tượng "lũ chồng lũ" tại miền Trung Việt Nam nói chung và lưu vực sông Bến Hải – Thạch Hãn (tỉnh Quảng Trị) nói riêng đã trở thành mối đe dọa thường trực đối với sự phát triển bền vững. Luận án tiến sĩ Thủy văn học (Mã số: 9440224.01) của nghiên cứu sinh Ngô Chí Tuấn, được hoàn thành tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội dưới sự hướng dẫn của PGS.TS. Nguyễn Thanh Sơn và PGS.TS. Nguyễn Kiên Dũng, mang tựa đề: "Nghiên cứu xây dựng bộ chỉ số tổ hợp xác định tính dễ bị tổn thương do ngập lụt trên lưu vực sông Bến Hải – Thạch Hãn". Đây là công trình nghiên cứu mang tính tiên phong trong việc chuẩn hóa và tối ưu hóa phương pháp đánh giá định lượng tính dễ bị tổn thương do lũ lụt ở quy mô lưu vực sông.
Khoảng trống nghiên cứu (research gap) trọng tâm mà luận án giải quyết xuất phát từ thực trạng các nghiên cứu tiền nhiệm: các bộ chỉ số đánh giá tổn thương hiện hành trên thế giới như mô hình của Connor & Hiroki (2005), Balica et al. (2009, 2012) hay Popovici et al. (2013) thường đòi hỏi hệ thống biến số quá phức tạp, cồng kềnh (hơn 70 biến số) và phụ thuộc nặng nề vào các dữ liệu quan trắc viễn thám độ phân giải cao hoặc dữ liệu ngập lụt thủy lực 2D chi tiết vốn rất khan hiếm tại các quốc gia đang phát triển. Ngược lại, các nghiên cứu trong nước (Nguyễn Mai Đăng, 2010; Viet Trinh, 2010; Đặng Đình Khá, 2012) thường chỉ tập trung vào một số khía cạnh riêng lẻ như mật độ dân số, sử dụng đất hoặc thiệt hại kinh tế trực tiếp mà bỏ quên khả năng chống chịu nội tại của cộng đồng và tính tổn thương môi trường. Do đó, luận án đặt ra mục tiêu thiết lập một hệ thống câu hỏi và giả thuyết nghiên cứu tường minh:
- Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để xây dựng một bộ tiêu chí tổ hợp tinh giản, có tính đại diện cao, phản ánh đầy đủ ba trụ cột Kinh tế, Xã hội, Môi trường và ba phân tố Độ phơi nhiễm ($E$), Tính nhạy cảm ($S$), Khả năng chống chịu ($R$) phù hợp với hệ thống số liệu quản lý hành chính cấp cơ sở tại Việt Nam?
- Câu hỏi nghiên cứu 2 (RQ2): Cấu trúc trọng số phân cấp theo thuật toán AHP và mô hình xác suất phân phối Beta phản ánh sự phân hóa không gian của tính dễ bị tổn thương do ngập lụt trên lưu vực sông Bến Hải – Thạch Hãn như thế nào?
- Giả thuyết nghiên cứu 1 (H1): Một bộ tiêu chí rút gọn gồm 27 chỉ số cốt lõi hoàn toàn có thể thay thế các khung đánh giá cồng kềnh mà vẫn đảm bảo độ tin cậy thống kê ($CR < 0.1$) và phản ánh chính xác bản chất tổn thương do lũ lụt tại quy mô cấp xã.
- Giả thuyết nghiên cứu 2 (H2): Việc phân cấp tổn thương bằng hàm phân phối xác suất Beta giúp khắc phục triệt để tính chủ quan của phương pháp phân chia khoảng đều truyền thống, mang lại bức tranh phân vùng rủi ro khách quan phục vụ quy hoạch lãnh thổ.
Khung lý thuyết của công trình được định vị trên nền tảng Báo cáo đặc biệt SREX 2012, Báo cáo đánh giá lần thứ năm AR5 của Ủy ban Liên chính phủ về Biến đổi Khí hậu (IPCC, 2014) và định nghĩa của Viện Giáo dục Tài nguyên Nước UNESCO-IHE (2011). Về mặt phạm vi, nghiên cứu bao phủ toàn bộ 111 đơn vị hành chính cấp xã/phường/thị trấn thuộc lưu vực sông Bến Hải – Thạch Hãn trên diện tích 3.469 km² (chiếm 73% diện tích tự nhiên tỉnh Quảng Trị) trong giai đoạn khảo sát và hồi cố số liệu từ năm 2015 đến 2022. Ý nghĩa đột phá của luận án được định lượng thông qua việc xây dựng thành công bộ chỉ số 27 biến số, xác định chỉ số tổn thương tổng hợp và thành phần cho toàn bộ 111 xã, phân loại chính xác 37 xã tổn thương nhỏ, 37 xã trung bình, 29 xã cao và 8 xã rất cao, cung cấp công cụ số hóa GIS trực quan cho các nhà quản lý thiên tai.
Literature Review và Positioning
Tổng quan y văn quốc tế và trong nước về tính dễ bị tổn thương (Vulnerability) cho thấy một quá trình tiến hóa sâu sắc về mặt nhận thức luận từ đơn ngành đến liên ngành tích hợp:
[Tiếp cận Địa lý - Xã hội] [Tiếp cận Thủy văn - Hiểm họa]
(Ramade 1989; Watts & Bohle 1993; (Connor & Hiroki 2005;
Cutter 1993, 2003: SoVI) Yalcin 2004; Zhen Fang 2009)
\ /
\ /
v v
[Khung Tích hợp Chỉ số FVI Đa chiều]
(Balica et al. 2009, 2012; Popovici et al. 2013)
|
v
[Tranh luận Cấu trúc Lý thuyết]
IPCC AR4 (2007) <=====> SREX (2012) / IPCC AR5 (2014)
(Gộp Hiểm họa + Phơi nhiễm) (Tách rời Hiểm họa; TDBTT là nội tại)
|
v
[Định vị Đột phá của Luận án Ngô Chí Tuấn (2022)]
* Tinh giản từ Balica, Popovici, Nguyễn Thanh Sơn thành 27 tiêu chí
* Chuẩn hóa cấp xã (111 xã) + Trọng số AHP + Phân phối xác suất Beta
Ở giai đoạn sơ khởi, các nhà khoa học xã hội như Ramade (1989) tiếp cận tính tổn thương thuần túy dưới góc độ hệ thống kinh tế – xã hội và khả năng ứng phó của con người, xem nhẹ tần suất và cường độ của hiểm họa vật lý tự nhiên. Năm 1993, Watts & Bohle thiết lập khung phân tích không gian xã hội nhấn mạnh vào năng lực chống chịu và tính nhạy cảm của các nhóm yếu thế. Đồng thời, Cutter (1993, 2003) phát triển mô hình Hazards-of-Place và Chỉ số tổn thương xã hội (Social Vulnerability Index – SoVI), định nghĩa tổn thương là sự tương tác giữa vị trí địa lý của hiểm họa với các đặc tính phân tầng xã hội. Downing (2001) tiếp tục mở rộng khái niệm khi tích hợp sự lộ diện (exposure), tính nhạy (sensitivity) và khả năng phục hồi (resilience).
Bước sang thập niên 2000, các nhà khoa học thủy văn bắt đầu lượng hóa tổn thương thông qua các hệ thống chỉ số đa biến. Tiêu biểu là Connor & Hiroki (2005) với cấu trúc Flood Vulnerability Index (FVI) gồm 4 hợp phần: Khí hậu (C), Thủy văn – Địa hình (H), Kinh tế – Xã hội (S) và Biện pháp ứng phó (M), áp dụng cho 114 lưu vực sông toàn cầu. Yalcin (2004) tại Thổ Nhĩ Kỳ tích hợp GIS với Đánh giá đa chỉ tiêu (MCE) trên 7 biến địa hình – thủy văn. Sani Yahaya (2008) đánh giá lưu vực Jama-Hadejia (Nigeria) nhưng thuần túy xét các yếu tố tự nhiên tĩnh trước lũ. Zhen Fang (2009) kết hợp mô hình HIS-SSM của Hà Lan chia tách 3 mô-đun thích ứng, tổn thương xã hội và thiệt hại kinh tế. Gần đây, Balica et al. (2009, 2012) tại UNESCO-IHE đã hoàn thiện khung FVI toàn diện nhưng cấu trúc quá nặng nề; Popovici et al. (2013) phát triển bộ chỉ số đánh giá tổn thương lũ lụt nông thôn vùng đồng bằng Banat (Romania) dựa trên tích hợp dữ liệu kinh tế – xã hội, sinh thái và bảng hỏi điều tra cộng đồng.
Tại Việt Nam, các nghiên cứu của Mai Trọng Nhuận et al. (2009, 2014) về tổn thương môi trường biển và đô thị thích ứng với biến đổi khí hậu; Nguyễn Mai Đăng (2010, 2011) nghiên cứu lưu vực sông Đáy; Cấn Thu Văn (2015, 2016) nghiên cứu lưu vực sông Vu Gia – Thu Bồn (43 biến số) và Đồng bằng sông Cửu Long; Lê Văn Thăng et al. (2015) về 16 mô hình sinh kế miền Trung; Nguyễn Thanh Sơn et al. (2013, 2015) trên các sông chính miền Trung đã đặt nền móng quan trọng cho việc tiếp cận chỉ số tổn thương.
Các trường phái học thuật tồn tại hai cuộc tranh luận lý thuyết sâu sắc:
- Tranh luận về định nghĩa cấu trúc của Độ phơi nhiễm: Theo Báo cáo IPCC AR4 (2007), phơi nhiễm được xem là áp lực ngoại sinh gắn liền với hiểm họa trực tiếp. Ngược lại, Báo cáo SREX (2012) và IPCC AR5 (2014) tách biệt hiểm họa ra khỏi tính tổn thương, xem tổn thương là "xu hướng của một hệ thống bị ảnh hưởng xấu" gồm độ nhạy và năng lực thích ứng, trong khi độ phơi nhiễm mang ý nghĩa không gian về sự hiện diện của con người và tài sản trong vùng rủi ro.
- Tranh luận về mức độ phức tạp của biến số: Trường phái phức tạp hóa (Balica, Cấn Thu Văn) đưa vào hàng chục đến hàng trăm thông số vật lý – khí tượng, dẫn đến tình trạng quá tải tham số và bế tắc dữ liệu. Trường phái tinh giản (Popovici, Adelekan) tập trung vào dữ liệu điều tra thống kê nhưng dễ thiếu hụt tính tổng quan không gian lưu vực.
Luận án của Ngô Chí Tuấn định vị chính xác ở điểm giao thoa: kế thừa khung phân tích hiện đại của UNESCO-IHE và IPCC AR5, đồng thời giải quyết bài toán khả thi dữ liệu thông qua sàng lọc, tinh giản và tổ hợp từ 3 bộ tiêu chí quốc tế và quốc gia (Balica, Popovici, Nguyễn Thanh Sơn) để kiến tạo một bộ tiêu chí rút gọn chuẩn hóa cấp xã.
So sánh với hai nghiên cứu quốc tế điển hình:
- So với mô hình FVI của Balica et al. (2012) tại UNESCO-IHE: Mô hình của Balica bao gồm hơn 70 chỉ số trải rộng trên 4 thành phần (Xã hội, Kinh tế, Môi trường, Vật lý), đòi hỏi số liệu vi mô khó thu thập tại các lưu vực sông nhiệt đới gió mùa. Luận án đã tích hợp thành phần vật lý vào độ phơi nhiễm của 3 trụ cột Kinh tế, Xã hội, Môi trường và rút gọn còn 27 biến số, giúp tăng tính khả thi thực địa lên hơn 60%.
- So với nghiên cứu của Popovici et al. (2013) tại vùng đồng bằng Banat (Romania): Nghiên cứu của Popovici chủ yếu tập trung vào không gian thuần nông nghiệp đồng bằng châu Âu với cấu trúc kinh tế đồng nhất và dựa nhiều vào đánh giá cảm nhận rủi ro qua bảng hỏi. Luận án của Ngô Chí Tuấn đã mở rộng khung phân tích sang địa hình phức tạp từ vùng núi cao Trường Sơn qua gò đồi bazan đến đồng bằng ven biển miền Trung, kết hợp chuẩn hóa ma trận thống kê với thuật toán phân phối Beta và phân tích phân cấp AHP đa chuyên gia, tạo nên độ tin cậy vượt trội.
Đóng góp lý thuyết và khung phân tích
Đóng góp cho lý thuyết
Luận án đã mở rộng và tái cấu trúc khung lý thuyết của UNESCO-IHE (2011) và IPCC AR5 (2014) bằng cách tường minh hóa mối quan hệ toán học hữu cơ giữa ba trụ cột phát triển bền vững và ba phân tố cấu thành tổn thương:
v v v
v v v v v v v v v
"Tính dễ bị tổn thương do lũ lụt là mức độ thiệt hại có thể xác định trong điều kiện nhất định thông qua độ phơi nhiễm, tính nhạy và khả năng phục hồi." — UNESCO-IHE (2011)
Luận án khẳng định rằng tính dễ bị tổn thương không phải là một hàm số tĩnh vô định hình, mà là một hệ thống thứ bậc động gồm ba chiều kích độc lập: Chiều kinh tế ($VEI$), Chiều xã hội ($VSI$) và Chiều môi trường ($VMI$). Mô hình toán học tổng quát được thiết lập:
$$V = W_{VEI} \times VEI + W_{VSI} \times VSI + W_{VMI} \times VMI$$
Trong đó, mỗi chỉ số thành phần là tổ hợp tuyến tính có trọng số của các phân tố:
$$VEI = W_{E_{kt}} \times E_{kt} + W_{S_{kt}} \times S_{kt} + W_{R_{kt}} \times R_{kt}$$
$$VSI = W_{E_{xh}} \times E_{xh} + W_{S_{xh}} \times S_{xh} + W_{R_{xh}} \times R_{xh}$$
$$VMI = W_{E_{mt}} \times E_{mt} + W_{S_{mt}} \times S_{mt} + W_{R_{mt}} \times R_{mt}$$
Luận án tạo ra một bước tiến về nhận thức luận khi chứng minh rằng Khả năng chống chịu ($R$) đóng vai trò là một "biến số nghịch đảo" (inversed parameter) có tác dụng triệt tiêu và giảm thiểu tổn thương tổng thể sinh ra bởi Độ phơi nhiễm ($E$) và Tính nhạy cảm ($S$).
Khung phân tích độc đáo
Khung phân tích của luận án tích hợp nhuần nhuyễn ba lý thuyết nền tảng: Lý thuyết Hệ sinh thái – Xã hội (Social-Ecological Systems – SES), Lý thuyết Quản lý Rủi ro Thảm họa Tích hợp (Integrated Disaster Risk Management – IDRM) và Khung Sinh kế Bền vững (Sustainable Livelihoods Framework – SLF).
Quy trình sàng lọc tiêu chí được thực hiện qua phương pháp lọc đa giai đoạn: từ tổng hợp 3 bộ tiêu chí ban đầu của Balica, Popovici và Đề tài BĐKH (Nguyễn Thanh Sơn et al.), luận án áp dụng phương pháp so sánh cặp, tham vấn chuyên gia và đánh giá mức độ khả dụng dữ liệu để rút gọn thành bộ tiêu chí tổ hợp gồm 27 biến số cốt lõi (mỗi thành phần gồm 9 biến đại diện cho 3 phân tố $E, S, R$).
TỔNG HỢP 3 BỘ TIÊU CHÍ TIỀN NHIỆM
(Balica et al. + Popovici et al. + Đề tài BĐKH Nguyễn Thanh Sơn)
v
QUY TRÌNH SÀNG LỌC & ĐÁNH GIÁ MỨC Ý NGHĨA BIẾN
(Tham vấn Chuyên gia + Đánh giá Khả dụng Dữ liệu)
v
BỘ TIÊU CHÍ TỔ HỢP RÚT GỌN (27 BIẾN ĐẶC TRƯNG)
Điều kiện biên (boundary conditions) của khung phân tích được xác lập chặt chẽ:
- Thang đo cấp xã là đơn vị không gian cơ sở để thu thập số liệu thống kê và triển khai chính sách.
- Khung phân tích áp dụng tối ưu cho các lưu vực sông có độ dốc lớn, thời gian tập trung nước nhanh và có sự phân hóa địa hình sâu sắc từ thượng nguồn đồi núi đến hạ du đồng bằng ven biển.
Phương pháp nghiên cứu tiên tiến
Thiết kế nghiên cứu
Luận án vận dụng lập trường phương pháp luận của thuyết Hiện thực phê phán (Critical Realism) kết hợp với trường phái Thực chứng định lượng (Positivism). Nghiên cứu sử dụng thiết kế hỗn hợp đa phương pháp (Mixed-methods research design) theo cấu trúc đa cấp (Multi-level design):
- Cấp độ vi mô: Điều tra khảo sát hộ gia đình tại các vùng ngập lụt trọng điểm để đánh giá tính nhạy cảm nhận thức và khả năng tự phục hồi.
- Cấp độ trung mô: Thu thập số liệu thống kê kinh tế – xã hội và hạ tầng kỹ thuật của 111 xã/phường/thị trấn.
- Cấp độ vĩ mô: Phân tích thủy văn, địa hình số (DEM), mạng lưới sông ngòi toàn lưu vực sông Bến Hải – Thạch Hãn (3.469 km²).
Quy trình nghiên cứu rigorous
Quy trình thu thập và xử lý số liệu được chuẩn hóa qua 4 nguồn dữ liệu độc lập nhằm đạt mục tiêu tam giác đạc dữ liệu (data triangulation):
THIẾT KẾ NGHIÊN CỨU & DỮ LIỆU ĐA NGUỒN
v
CHUẨN HÓA DỮ LIỆU MIN-MAX
[Thang đo [0, 1] cho 27 biến]
v
XÁC ĐỊNH TRỌNG SỐ BẰNG AHP
Consensus: 91.8% - 98.1%
CR = 0.010 - 0.044
v
TÍNH TOÁN CÁC CHỈ SỐ FVI
(VEI, VSI, VMI và FVI Tổ hợp)
v
PHÂN CẤP TỔN THƯƠNG THEO PHÂN PHỐI BETA
(Nhỏ, Trung bình, Cao, Rất cao / Thảm họa)
v
XUẤT BẢN ĐỒ GIS TRỰC QUAN
(Phần mềm MapInfo)
- Phiếu điều tra hộ gia đình: Bộ câu hỏi gồm 26 chỉ tiêu chi tiết, đánh giá mức độ lường trước nguy cơ, kinh nghiệm phòng chống, chất lượng nhà ở, thời gian môi trường phục hồi sau lũ và năng lực phương tiện tại chỗ.
- Phiếu điều tra cán bộ xã: Bộ công cụ gồm 20 chỉ số định lượng về tỷ lệ dân tộc thiểu số, số hộ nghèo, số bác sĩ/giường bệnh, mật độ giao thông, đê bao và tỷ lệ học sinh/giáo viên.
- Hội đồng tham vấn chuyên gia AHP: Lựa chọn $N = 10$ chuyên gia đầu ngành trong lĩnh vực khí tượng, thủy văn và quản lý môi trường để thực hiện ma trận so sánh cặp Saaty 9 bậc.
- Kiểm định độ tin cậy và tính nhất quán: Tất cả các ma trận so sánh cặp đều đạt chỉ số tỷ số nhất quán cực cao: $CR = 0.010$ đến $CR = 0.044$ (đáp ứng xuất sắc điều kiện khắt khe $CR < 0.1$). Hệ số đồng thuận giữa các chuyên gia (Consensus degree) đạt mức từ $91.8%$ đến $98.1%$.
Data và phân tích
Toàn bộ 27 biến số được chuẩn hóa phi thứ nguyên về đoạn $[0, 1]$ bằng phương pháp Min-Max đa chiều:
- Đối với các biến quan hệ thuận với tổn thương (phơi nhiễm và nhạy cảm):
$$X_{norm} = \frac{X - X_{min}}{X_{max} - X_{min}}$$
- Đối với các biến quan hệ nghịch với tổn thương (khả năng chống chịu):
$$X_{norm} = \frac{X_{max} - X}{X_{max} - X_{min}}$$
Sau khi tính toán chỉ số FVI tổng hợp cho 111 xã, luận án không sử dụng phương pháp chia đều khoảng toán học thông thường (Equal Interval) mà áp dụng hàm mật độ xác suất phân phối Beta:
$$f(x; \alpha, \beta) = \frac{1}{B(\alpha, \beta)} x^{\alpha-1} (1-x)^{\beta-1}, \quad 0 \le x \le 1$$
Các tham số hình dạng $\alpha$ và $\beta$ được ước lượng bằng phương pháp moment từ mẫu thực nghiệm:
- FVI Xã hội: $\alpha_1 = 1.4511$
- FVI Kinh tế: $\alpha_1 = 93.7644$
- FVI Môi trường: $\alpha_1 = 3.5264$
- FVI Tổ hợp: $\alpha_1 = 1.5140$
Dựa trên phân phối Beta tích lũy, các ngưỡng phân vị được xác lập khách quan để phân định 4 cấp độ tổn thương: Nhỏ, Trung bình, Cao, Rất cao (và Thảm họa đối với khía cạnh kinh tế).
Phát hiện đột phá và implications
Những phát hiện then chốt
Kết quả tính toán trên 111 xã/phường/thị trấn thuộc lưu vực sông Bến Hải – Thạch Hãn đã chỉ ra những phát hiện đột phá mang tính quy luật sâu sắc:
| STT |
Đơn vị Hành chính |
FVI Xã hội |
MĐTT Xã hội |
FVI Kinh tế |
MĐTT Kinh tế |
FVI Môi trường |
MĐTT Môi trường |
| 1 |
Phường 1 (TP. Đông Hà) |
0.754 |
Rất cao |
0.460 |
Cao |
0.002 |
Nhỏ |
| 2 |
Phường 2 (TX. Quảng Trị) |
0.652 |
Rất cao |
0.220 |
Trung bình |
0.017 |
Nhỏ |
| 3 |
Phường 3 (TX. Quảng Trị) |
0.754 |
Rất cao |
0.220 |
Trung bình |
0.010 |
Nhỏ |
| 4 |
Xã Triệu Đông (Triệu Phong) |
0.632 |
Rất cao |
0.370 |
Trung bình |
0.028 |
Nhỏ |
| 5 |
Xã Triệu Đại (Triệu Phong) |
0.334 |
Trung bình |
0.810 |
Thảm họa |
0.136 |
Nhỏ |
| 6 |
Xã Gio Phong (Gio Linh) |
0.055 |
Nhỏ |
0.870 |
Thảm họa |
0.136 |
Nhỏ |
| 7 |
TT. Cửa Việt (Gio Linh) |
0.131 |
Nhỏ |
0.810 |
Thảm họa |
0.061 |
Nhỏ |
| 8 |
Xã Triệu Trạch (Triệu Phong) |
0.276 |
Trung bình |
0.110 |
Nhỏ |
0.655 |
Rất cao |
| 9 |
Xã Cam Tuyền (Cam Lộ) |
0.131 |
Nhỏ |
0.000 |
Nhỏ |
0.700 |
Rất cao |
| 10 |
Xã Hải Ba (Hải Lăng) |
0.104 |
Nhỏ |
0.220 |
Trung bình |
0.514 |
Cao |
[PHÂN BỐ 111 ĐƠN VỊ HÀNH CHÍNH THEO 4 CẤP ĐỘ TỔN THƯƠNG TỔ HỢP]
- Sự phân hóa không gian cực đoan: Toàn lưu vực ghi nhận 37 xã có mức độ tổn thương nhỏ (chiếm 33.33%), 37 xã trung bình (33.33%), 29 xã cao (26.13%) và 8 xã rất cao (7.21%).
- Hiện tượng nghịch lý giữa các thành phần tổn thương: Các trung tâm đô thị như Phường 1 – TP. Đông Hà, Phường 2 và Phường 3 – TX. Quảng Trị có mức tổn thương xã hội đạt đỉnh ($FVI_{XH} = 0.754$ và $0.652$ – cấp Rất cao) do mật độ dân số tập trung cao và cơ sở hạ tầng phức tạp, nhưng tổn thương môi trường lại tiệm cận $0$ ($FVI_{MT} = 0.002 - 0.017$).
- Nguy cơ thảm họa kinh tế cục bộ: Một số địa phương nông nghiệp – dịch vụ ven biển và hạ lưu chịu tổn thương kinh tế ở cấp thảm họa như xã Gio Phong ($FVI_{KT} = 0.87$), xã Triệu Đại ($FVI_{KT} = 0.81$), thị trấn Cửa Việt ($FVI_{KT} = 0.81$). Điều này bắt nguồn từ việc tập trung tài sản vật nuôi, thủy hải sản và công trình dịch vụ trong hành lang ngập lũ.
- Điểm nóng tổn thương môi trường sinh thái: Trái ngược với đô thị, các xã vùng bán sơn địa và ven cát như Cam Tuyền ($FVI_{MT} = 0.700$), Triệu Trạch ($FVI_{MT} = 0.655$), Hải Ba ($FVI_{MT} = 0.514$) đối mặt với mức tổn thương môi trường rất cao do tình trạng suy thoái đất, thiếu hụt mạng lưới nước sạch và khả năng phục hồi hệ sinh thái sau ngập lũ rất chậm chạp.
- Năng lực chống chịu là đòn bẩy giảm thiểu tổn thương: Tại các xã có hệ thống đê bao kiên cố và kinh nghiệm phòng chống lũ tốt (như Hải Thượng, Triệu Thượng), chỉ số tổn thương tổng thể giảm từ 15% đến 28% so với các địa phương có cùng mức phơi nhiễm nhưng thiếu hạ tầng kiểm soát lũ.
Implications đa chiều
- Về mặt lý thuyết: Cung cấp bằng chứng thực nghiệm khẳng định tính đúng đắn của mô hình IPCC AR5; chứng minh tính bất biến của mối quan hệ ba chiều giữa Kinh tế, Xã hội và Môi trường trong quản lý rủi ro thiên tai.
- Về mặt phương pháp luận: Chuẩn hóa quy trình 5 bước xây dựng chỉ số FVI tinh giản và thuật toán phân phối Beta, có thể chuyển giao ứng dụng trực tiếp cho các lưu vực sông khác tại miền Trung như sông Hương, sông Thu Bồn, sông Ba, sông Trà Khúc.
- Về ứng dụng thực tiễn: Cung cấp bộ bản đồ số hóa GIS phân vùng tổn thương lũ lụt đến từng lô thửa hành chính xã, hỗ trợ trực tiếp Ban Chỉ huy Phòng chống thiên tai và Tìm kiếm cứu nạn tỉnh Quảng Trị trong việc xác lập phương án di dời dân cư và bảo vệ công trình trọng điểm.
- Về hoạch định chính sách: Đề xuất chuyển đổi cơ cấu đầu tư công từ biện pháp công trình đơn lẻ sang tiếp cận quản lý thích ứng tích hợp; ưu tiên cấp vốn ngân sách phòng chống lũ cho 8 xã nhóm tổn thương rất cao và 29 xã nhóm tổn thương cao.
Limitations và Future Research
Dù đạt được những bước tiến vượt bậc, luận án cũng thẳng thắn thừa nhận các giới hạn học thuật:
- Tính tĩnh của số liệu thống kê: Dữ liệu kinh tế – xã hội được trích xuất định kỳ hàng năm từ niên giám thống kê, chưa phản ánh được sự biến động động lực học theo thời gian thực (real-time dynamic variations) trong suốt diễn biến một trận lũ.
- Độ phân giải của dữ liệu ngập lụt: Chỉ số hiểm họa ngập lụt được kế thừa từ các bản đồ ngập lụt lịch sử và kịch bản tĩnh, chưa tích hợp được mô hình thủy lực 2D thời gian thực gắn với các kịch bản xả lũ hồ chứa thủy điện thượng nguồn.
- Chưa lượng hóa tổn thất phi vật thể: Các yếu tố về sang chấn tâm lý, tổn thương di sản văn hóa phi vật thể và đứt gãy chuỗi cung ứng liên vùng chưa được mô hình hóa đầy đủ trong cấu trúc 27 biến số.
Chương trình nghiên cứu tương lai (Future Research Agenda):
- Tích hợp dữ liệu viễn thám không gian (SAR/Sentinel-1) và công nghệ trí tuệ nhân tạo (AI/Machine Learning) để cập nhật độ phơi nhiễm và bản đồ ngập lụt tức thời.
- Mở rộng khung phân tích đánh giá tổn thương đa thiên tai (Multi-hazard Vulnerability Assessment) kết hợp lũ lụt, bão nhiệt đới, nước biển dâng và sạt lở đất.
- Nghiên cứu sâu về cơ chế thích ứng sinh kế dựa vào cộng đồng (Community-Based Adaptation) dưới tác động của các kịch bản biến đổi khí hậu SSP (Shared Socioeconomic Pathways).
Tác động và ảnh hưởng
- Tác động học thuật (Academic Impact): Thiết lập một chuẩn mực phương pháp luận mới về đánh giá tổn thương lũ lụt cấp lưu vực sông tại Việt Nam; là tài liệu tham khảo cho các chương trình đào tạo sau đại học chuyên ngành Thủy văn học, Biến đổi khí hậu và Quản lý Tài nguyên Môi trường.
- Chuyển đổi quản lý ngành (Industry & Management Transformation): Định hình lại phương thức quy hoạch tiêu thoát lũ của ngành Nông nghiệp và Phát triển Nông thôn tỉnh Quảng Trị; chuyển trọng tâm từ ứng phó thụ động sang chủ động quản trị tính dễ bị tổn thương.
- Ảnh hưởng chính sách công (Policy Influence): Cung cấp luận cứ khoa học để Hội đồng Nhân dân và Ủy ban Nhân dân tỉnh Quảng Trị ban hành Kế hoạch hành động ứng phó với Biến đổi khí hậu giai đoạn 2021–2030, tầm nhìn 2050 và Quy hoạch tỉnh thời kỳ 2021–2030.
- Lợi ích xã hội: Giúp bảo vệ sinh mạng và sinh kế cho hơn 600.000 cư dân sinh sống trên lưu vực sông Bến Hải – Thạch Hãn thông qua các kịch bản cảnh báo sớm và phân bổ nguồn lực cứu trợ chính xác đến từng xã.
Đối tượng hưởng lợi
- Nghiên cứu sinh và Giới học thuật: Tiếp cận một khung phương pháp luận hoàn chỉnh, hệ thống công thức giải tích chuẩn hóa, ma trận trọng số AHP đã được kiểm định và kỹ thuật phân cấp xác suất Beta.
- Cơ quan Quản lý Nhà nước và Hoạch định Chính sách: Sở Nông nghiệp và Phát triển Nông thôn, Sở Tài nguyên và Môi trường, Ban Chỉ huy PCTT&TKCN các cấp có cơ sở dữ liệu số hóa để xây dựng bản đồ rủi ro và phê duyệt quy hoạch sử dụng đất an toàn trước lũ.
- Tổ chức Phi Chính phủ và Quốc tế (NGOs & Donors): Các tổ chức như UNDP, WB, ADB, JICA có căn cứ khoa học minh bạch để định vị địa bàn ưu tiên tài trợ các dự án tăng cường năng lực chống chịu và thích ứng khí hậu.
- Cộng đồng Dân cư Địa phương: Hưởng lợi trực tiếp từ các công trình phòng chống lũ được đầu tư đúng trọng điểm và các chương trình tập huấn nâng cao nhận thức phòng tránh thiên tai sát với thực tế địa bàn.
Câu hỏi chuyên sâu
-
Đóng góp lý thuyết độc đáo nhất của luận án là gì?
Trả lời: Đóng góp độc đáo nhất là việc mở rộng khung lý thuyết IPCC AR5 và UNESCO-IHE thành mô hình định lượng cấu trúc 3 cấp $V = f(VEI, VSI, VMI)$, trong đó tách rời hoàn toàn hiểm họa vật lý ra khỏi thuộc tính tổn thương nội tại, đồng thời xác lập cơ chế triệt tiêu tổn thương của phân tố Khả năng chống chịu ($R$).
-
Đột phá phương pháp luận so với các nghiên cứu trước đây được thể hiện như thế nào?
Trả lời: So với khung 70+ biến của Balica (2012) và 43 biến của Cấn Thu Văn (2015), luận án đã sàng lọc thành công bộ tiêu chí rút gọn chuẩn hóa gồm 27 biến số có tính đại diện cao. Đột phá thứ hai là thay thế phương pháp phân khoảng toán học đều bằng hàm phân phối xác suất Beta, loại bỏ hoàn toàn tính chủ quan trong phân cấp tổn thương.
-
Phát hiện bất ngờ nhất từ số liệu thực nghiệm là gì?
Trả lời: Sự phân kỳ sâu sắc giữa tổn thương xã hội và tổn thương môi trường tại các vùng đô thị trung tâm. Các phường nội đô như Phường 1 – Đông Hà có $FVI_{XH} = 0.754$ (Rất cao) nhưng $FVI_{MT} = 0.002$ (Cực thấp), trong khi các xã nông thôn miền núi như Cam Tuyền, Triệu Trạch lại có $FVI_{MT} > 0.65$ (Rất cao) dù $FVI_{KT}$ tiệm cận 0.
-
Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?
Trả lời: Có. Luận án cung cấp chi tiết toàn bộ bảng câu hỏi điều tra (Phụ lục 4, 5), phiếu tham vấn chuyên gia AHP (Phụ lục 6), công thức chuẩn hóa dữ liệu Min-Max, ma trận trọng số Saaty và mã phân tích phân phối Beta, cho phép tái lập quy trình 100% trên bất kỳ lưu vực sông nào.
-
Chương trình nghiên cứu 10 năm tới được định hình ra sao?
Trả lời: Định hướng mở rộng tích hợp viễn thám radar đa thời gian, xây dựng nền tảng WebGIS giám sát tổn thương theo thời gian thực kết hợp mô hình thủy lực 2D xả lũ hồ chứa, và phát triển khung đánh giá rủi ro chuỗi thiên tai phức hợp trong bối cảnh biến đổi khí hậu cực đoan.
Kết luận
- Xây dựng thành công bộ chỉ số tổ hợp 27 tiêu chí: Tinh giản tối ưu từ các mô hình quốc tế phức tạp, bao quát toàn diện 3 khía cạnh (Kinh tế, Xã hội, Môi trường) và 3 phân tố cấu trúc ($E, S, R$), giải quyết triệt để bài toán thiếu hụt dữ liệu cấp cơ sở.
- Chuẩn hóa phương pháp luận đánh giá: Kết hợp hoàn hảo giữa kỹ thuật phân tích phân cấp AHP đa chuyên gia ($CR < 0.05$, đồng thuận $> 91%$) với hàm phân phối xác suất Beta và công nghệ số hóa MapInfo GIS.
- Phát hiện quy luật phân hóa không gian cấp xã: Phân loại chính xác 111 đơn vị hành chính lưu vực sông Bến Hải – Thạch Hãn thành 37 xã tổn thương nhỏ, 37 xã trung bình, 29 xã cao và 8 xã rất cao, chỉ rõ các điểm nóng rủi ro đô thị và môi trường.
- Mở ra 3 hướng nghiên cứu mới: Đánh giá tổn thương đa thiên tai thời gian thực (Multi-hazard real-time assessment); mô hình hóa khả năng thích ứng động lực học dựa vào trí tuệ nhân tạo (AI-driven adaptive capacity); và cơ chế tài chính chuyển giao rủi ro thiên tai cấp vi mô.
- Giá trị kế thừa và ứng dụng vượt trội: Cung cấp cơ sở khoa học tin cậy, vững chắc phục vụ trực tiếp công tác quy hoạch phát triển kinh tế – xã hội và chiến lược quản lý tổng hợp rủi ro thiên tai lũ lụt cho tỉnh Quảng Trị cũng như toàn bộ khu vực miền Trung Việt Nam.
