Tổng quan luận án

Luận án tiến sĩ kỹ thuật với đề tài "Nghiên cứu phát triển mô hình trí tuệ nhân tạo trong phân vùng nguy cơ lũ quét ở Việt Nam" do nghiên cứu sinh Ngô Thị Phương Thảo thực hiện tại Bộ môn Đo ảnh và Viễn thám, Khoa Trắc địa - Bản đồ và Quản lý đất đai, Trường Đại học Mỏ - Địa chất, dưới sự hướng dẫn khoa học của TS. Nguyễn Quang Khánh (Trường Đại học Mỏ - Địa chất) và GS.TS. Bùi Tiến Diệu (University of South-Eastern Norway, Na Uy). Luận án thuộc chuyên ngành Kỹ thuật Trắc địa - Bản đồ, mã số 9520503, bảo vệ năm 2024 tại Hà Nội.

Tính cấp thiết và khoảng trống nghiên cứu

Lũ lụt và lũ quét là những hiểm họa thiên tai gây thiệt hại nghiêm trọng về tính mạng con người và tài sản trên quy mô toàn cầu cũng như tại Việt Nam. Theo các báo cáo quốc tế, Việt Nam nằm trong nhóm 10 quốc gia chịu tác động nặng nề nhất do biến đổi khí hậu; trong đó lũ lụt đã ảnh hưởng đến khoảng 71% dân số và 59% diện tích đất đai. Giai đoạn 1996–2015, tổng thiệt hại do bão và lũ ước tính chiếm khoảng 0,62% GDP quốc gia (tương đương 2,12 tỷ USD). Đặc biệt tại các tỉnh miền núi phía Bắc, hiện tượng lũ quét xảy ra thường xuyên với cường độ dữ dội, diễn biến nhanh (thường dưới 6 giờ), bắt nguồn từ các trận mưa lớn cực đoan kết hợp với địa hình phân cắt mạnh, độ dốc lớn và các hoạt động dân sinh như phá rừng, sử dụng đất thiếu quy hoạch.

Khoảng trống nghiên cứu được luận án chỉ ra gồm hai điểm chính:

  1. Hạn chế của các mô hình truyền thống: Các mô hình thủy văn - thủy lực mô phỏng quá trình mưa - dòng chảy đòi hỏi hệ thống chuỗi dữ liệu quan trắc dài hạn và mạng lưới trạm quan trắc dày đặc. Tại các vùng núi hiểm trở ở Việt Nam, mật độ trạm đo rất thưa thớt, khiến các mô hình truyền thống khó đạt độ chính xác cao trong cảnh báo không gian.
  2. Khó khăn trong thành lập bản đồ hiện trạng lũ quét: Việc điều tra vết lũ ngoài thực địa bằng GPS thủ công rất tốn kém và mất thời gian, trong khi ảnh viễn thám quang học thường bị hạn chế bởi mây mù dày đặc trong mùa mưa bão. Ảnh radar trước đây có chu kỳ lặp dài và độ phân giải thấp. Do đó, việc ứng dụng ảnh Radar SAR Sentinel-1 kết hợp với các thuật toán học máy và tối ưu hóa metaheuristic để giải quyết bài toán phân loại không gian ("on-off") là hướng tiếp cận cần thiết nhằm lấp đầy khoảng trống này.

Mục tiêu nghiên cứu

Mục tiêu tổng quát của luận án là xây dựng các mô hình trí tuệ nhân tạo có độ chính xác cao phục vụ phân vùng nguy cơ lũ quét, trên cơ sở ứng dụng ảnh viễn thám Radar Sentinel-1, kỹ thuật hệ thống thông tin địa lý (GIS) và địa thống kê; thực nghiệm áp dụng cho khu vực vùng núi tỉnh Lào Cai.

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu:
    1. Dữ liệu địa không gian đa thời gian từ ảnh viễn thám Radar Sentinel-1A phục vụ chiết xuất hiện trạng lũ quét.
    2. Cơ sở dữ liệu GIS gồm các yếu tố: địa hình, địa mạo, loại đất, địa chất, khí hậu và thủy văn.
    3. Thuật toán chiết xuất điểm lũ quét dựa trên kỹ thuật phát hiện thay đổi (Change Detection).
    4. Các mô hình trí tuệ nhân tạo, thuật toán máy học và kỹ thuật tối ưu hóa metaheuristic (FA-LM-ANN, PSO-ELM, FURIA-GA-Ensemble).
  • Phạm vi không gian: Tỉnh Lào Cai, trong đó khu vực nghiên cứu thực nghiệm tập trung tại hai huyện Bắc Hà và Bảo Yên với tổng diện tích 1.510,4 km² (tọa độ từ 22°5′ đến 22°40′ vĩ độ Bắc và từ 104°10′ đến 105°37′ kinh độ Đông; độ cao biến thiên từ 38,9 m đến 1.878,69 m; 85,4% diện tích có độ dốc từ 10° đến 40°).
  • Phạm vi khoa học: Các thuật toán xử lý ảnh viễn thám SAR đa thời gian, phân tích dữ liệu không gian trên GIS, kỹ thuật phân loại nhị phân ("on-off") và các mô hình học máy tích hợp tối ưu hóa toàn cục.

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu lũ quét trên thế giới

Tổng quan y văn quốc tế được luận án hệ thống thành ba nhóm mô hình chính:

  • Nhóm 1 - Các mô hình phân tích thống kê: Sử dụng chuỗi số liệu mực nước và lưu lượng từ các trạm quan trắc để xây dựng phương trình hồi quy và ngoại suy (như mô hình Copula, mô hình Gumbel hỗn hợp, mô hình tích hợp Gumbel–Hougaard Copula, phân bố tần suất thống kê). Nhược điểm lớn là đòi hỏi chuỗi số liệu quan trắc liên tục, dài hạn.
  • Nhóm 2 - Các mô hình mô phỏng mưa - dòng chảy: Mô phỏng các quá trình thủy văn - thủy lực dòng chảy mặt dựa trên hệ phương trình vi phân và các tham số vật lý (tiêu biểu như HEC-HMS/RAS, HYDROTEL, Wetspa, SWAT, MIKE, KINEROS2/AGWA, HiResFlood-UCI). Nhóm mô hình này gặp trở ngại lớn khi áp dụng tại các lưu vực thiếu trạm đo.
  • Nhóm 3 - Các mô hình dựa trên giả thuyết thống kê "on-off" (Binary Pattern Recognition): Tiếp cận bài toán phân vùng nguy cơ lũ quét dưới dạng bài toán phân loại nhị phân hai lớp (lũ quét và không lũ quét), dựa trên nguyên lý: khu vực có điều kiện địa lý - môi trường tương đồng với nơi từng xảy ra lũ quét trong quá khứ sẽ có nguy cơ cao xảy ra lũ quét. Các phương pháp đã phát triển gồm: phân tích thứ bậc AHP, hồi quy Logistic, trọng số Bayesian, hệ mờ nơ-ron (ANFIS), cây quyết định (Decision Trees), máy vectơ hỗ trợ (SVM), mạng nơ-ron nhân tạo nhiều lớp (MLPNN), rừng ngẫu nhiên (Random Forest).

Điểm luận các nghiên cứu lũ quét tại Việt Nam

Luận án điểm lại các công trình tiêu biểu trong nước theo trật tự thời gian và phương pháp:

  • Cao Đặng Dư và cộng sự: Điều tra thống kê lũ quét giai đoạn 1951–2001 trên toàn quốc, xác định lượng mưa là yếu tố kích hoạt quan trọng nhất và đề xuất phân vùng nguy cơ thành 4 cấp (rất cao, cao, trung bình, thấp).
  • Đào Minh Đức và cộng sự: Nghiên cứu nguy cơ lũ quét lưu vực suối Nghĩa Đô (Bảo Yên, Lào Cai) bằng phương pháp thống kê, chỉ ra ngưỡng mưa 24 giờ vượt quá 29 mm là điều kiện đủ kích hoạt lũ quét.
  • Đào Văn Thánh; Nguyễn Ngọc Việt: Sử dụng kiến thức chuyên gia kết hợp điều tra hiện trạng để gán trọng số cho các bản đồ địa hình, địa chất, khí hậu; phương pháp mang tính định tính và phụ thuộc chủ quan vào chuyên gia.
  • Dương Thị Lợi và Đặng Phương Lan; Nguyễn Ngọc Thạch: Ứng dụng quy trình phân tích thứ bậc (AHP) kết hợp GIS và bổ sung nhân tố dòng chảy tích lũy để thành lập bản đồ phân vùng nguy cơ tại Vĩnh Phúc, Bắc Kạn.
  • Hồ Tiến Chung và cộng sự: Dự báo lũ quét và sạt lở đất dọc Quốc lộ 32 (Yên Bái - Lai Châu) giai đoạn 2008–2010.
  • Lã Thanh Hà và cộng sự: Tái hiện trận lũ quét lịch sử tại suối Phà (Văn Chấn, Yên Bái) năm 2005 bằng mô hình thủy văn - thủy lực HEC-HMS và HEC-RAS tích hợp trong ArcView.
  • Ngô Thị Phượng và cộng sự: Nghiên cứu hiện trạng và dự báo nguy cơ trượt lở, lũ bùn đá tại Hoàng Su Phì, Xín Mần, Yên Minh (Hà Giang).
  • Nguyễn Khắc Hải: Thử nghiệm ứng dụng mạng nơ-ron nhân tạo mô phỏng lưu lượng lũ quét cho một số lưu vực sông nhỏ (sông Dinh, Nậm La, sông Vệ).
  • Viện Khoa học Khí tượng Thủy văn và Biến đổi khí hậu: Thực hiện đề án điều tra, phân vùng và cảnh báo lũ quét miền núi Việt Nam giai đoạn 2010–2018 (giai đoạn 1 tại 14 tỉnh miền núi phía Bắc; giai đoạn 2 tại miền Trung và Tây Nguyên).
  • Bùi Tiến Diệu, Hoàng Nhật Đức và cộng sự: Phát triển các mô hình học máy lai ghép như Neural-Fuzzy tối ưu hóa và mô hình tích hợp thống kê Bayes (BayGmmKda) trong phân vùng không gian lũ quét không phụ thuộc trạm thủy văn.

Khoảng trống luận án lựa chọn giải quyết

Luận án tập trung giải quyết ba vấn đề kỹ thuật trọng tâm:

  1. Xây dựng quy trình tự động phát hiện và thành lập bản đồ hiện trạng lũ quét từ ảnh Radar SAR Sentinel-1A đa thời gian thông qua kỹ thuật phát hiện thay đổi (Change Detection), khắc phục tình trạng thiếu hụt dữ liệu vết lũ.
  2. Tích hợp thuật toán tối ưu hóa bầy đom đóm (Firefly Algorithm - FA) và Levenberg–Marquardt (LM) vào mạng nơ-ron nhân tạo (ANN) nhằm khắc phục hiện tượng rơi vào cực tiểu cục bộ và cải thiện tốc độ hội tụ.
  3. Ứng dụng thuật toán máy học cực đại (ELM) kết hợp tối ưu hóa bầy đàn (PSO), cùng mô hình học tập kết hợp Ensemble Learning (tích hợp FURIA, thuật toán di truyền GA và các kỹ thuật Bagging, AdaBoost, LogitBoost trên nền cây quyết định C4.5) để nâng cao độ chính xác phân loại không gian nguy cơ lũ quét.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết và khái niệm

  • Khung quản lý rủi ro thiên tai Sendai (UNISDR): Luận án vận dụng cấu trúc 3 chặng của Khung Sendai gồm: (1) Khu vực đánh giá và công tác chuẩn bị; (2) Phân tích rủi ro kỹ thuật; (3) Sử dụng kết quả và xây dựng quyết định.
  • Khái niệm lũ quét: Tổng hợp theo định nghĩa của WMO, AMS và các nhà khoa học Việt Nam, lũ quét được xác định là hiện tượng ngập lụt cục bộ xuất hiện đột ngột trong thời gian ngắn (dưới 6 giờ), dòng chảy xiết mang hàm lượng bùn đất, đá cao, có sức tàn phá lớn tại các lưu vực sườn dốc miền núi.
  • Giả thuyết thống kê phân loại nhị phân ("on-off"): Toàn bộ các pixel trong khu vực được phân thành 2 lớp: lớp lũ quét ($C_2 = +1$) và lớp không lũ quét ($C_1 = -1$). Xác suất đầu ra thuộc lớp lũ quét đóng vai trò là chỉ số nguy cơ lũ quét.

Phương pháp nghiên cứu thực tế

  1. Phương pháp Viễn thám:
    • Dữ liệu: 4 cảnh ảnh vệ tinh Sentinel-1A chế độ Interferometric Wide-Swath (IW), định dạng Level-1 GRD, phân cực kép VV, độ phân giải không gian 10 m, chụp vào các ngày 23/07/2017, 30/07/2017 (trước lũ quét) và 04/08/2017, 10/10/2017 (sau lũ quét).
    • Quy trình tiền xử lý trên phần mềm SNAP Toolbox: Lọc bỏ nhiễu nhiệt (Thermal Noise Removal); hiệu chỉnh bức xạ đưa về hệ số tán xạ ngược Sigma Nought ($\sigma^0$); lọc nhiễu đốm bằng bộ lọc Lee; nắn chỉnh hình học bằng mô hình số độ cao SRTM DEM về phép chiếu UTM Zone 48N - WGS84; tạo ảnh tổ hợp màu đa thời gian và khoanh vùng vết lũ quét; kiểm chứng thực địa bằng máy đo GPS cầm tay.
  2. Phương pháp GIS và Địa thống kê:
    • Xây dựng hệ thống cơ sở dữ liệu gồm 12 bản đồ thành phần: Độ cao (DEM), Độ dốc (Slope), Hướng dốc (Aspect), Độ cong địa hình (Curvature), Hình thái địa mạo (Toposhape), Chỉ số độ ẩm địa hình (TWI), Chỉ số năng lượng dòng chảy (SPI), Mật độ sông suối (Stream Density), Chỉ số thực vật (NDVI), Lượng mưa (Rainfall), Loại đất (Soil type), Thạch học (Lithology).
    • Kiểm tra đa cộng tuyến bằng hai chỉ số: Hệ số dung sai (Tolerance - TOL) và Hệ số phóng đại phương sai (Variance Inflation Factor - VIF).
  3. Phương pháp mô hình hóa học máy và tối ưu hóa:
    • Mô hình FA-LM-ANN: Mạng nơ-ron truyền thẳng nhiều lớp với hàm kích hoạt Log-sigmoid; thuật toán Firefly Algorithm (FA) tìm kiếm không gian toàn cục để khởi tạo không gian trọng số tối ưu; thuật toán Levenberg–Marquardt (LM) tối ưu hóa cục bộ dựa trên ma trận Jacobian để tăng tốc độ hội tụ và giảm thiểu hàm mất mát MSE.
    • Mô hình PSO-ELM: Mạng Extreme Learning Machine (ELM) với cấu trúc một lớp ẩn, các trọng số đầu ra được tính giải tích qua ma trận nghịch đảo Moore–Penrose; thuật toán tối ưu hóa bầy đàn (PSO) tối ưu hóa các vector trọng số đầu vào và độ lệch nhằm hạn chế việc gán ngẫu nhiên ban đầu.
    • Mô hình Ensemble Learning: Thuật toán mờ không thứ tự FURIA kết hợp giải thuật di truyền GA để tinh lọc luật mờ và lựa chọn thuộc tính; tích hợp kỹ thuật lấy mẫu lặp Bootstrap (Bagging), tăng cường thích ứng (AdaBoost) và tăng cường Logistic (LogitBoost) trên nền cây quyết định C4.5.
  4. Phương pháp đánh giá độ chính xác:
    • Phân chia tập dữ liệu thành tập huấn luyện (Training Set - 70%) và tập kiểm tra (Testing Set - 30%), kết hợp kiểm tra chéo (10-fold cross-validation và 5-fold cross-validation).
    • Các chỉ số đo lường hiệu suất: Sai số toàn phương trung bình (MSE), Căn bậc hai sai số toàn phương trung bình (RMSE), Sai số tuyệt đối trung bình (MAE), Sai số trung bình (ME).
    • Đánh giá phân loại: Ma trận nhầm lẫn (True Positive - TP, False Positive - FP, True Negative - TN, False Negative - FN), độ chính xác phân loại tổng thể, hệ số Kappa, đường cong đặc trưng máy thu (ROC) và diện tích dưới đường cong (AUC).
    • Kiểm định thống kê phi tham số: Kiểm nghiệm Wilcoxon signed-rank test để đánh giá mức độ khác biệt có ý nghĩa thống kê giữa các mô hình.
Yếu tố đầu vào / Bản đồ thành phần Nguồn dữ liệu gốc Ý nghĩa trong phân vùng lũ quét
Ảnh Radar Sentinel-1A (GRD) ESA Copernicus Open Access Hub Chiết xuất vết lũ quét và thành lập bản đồ hiện trạng
Độ cao (DEM), Độ dốc, Hướng dốc Dữ liệu SRTM DEM (độ phân giải 30 m) Xác định năng lượng thế năng, hướng và vận tốc dòng chảy
Độ cong (Curvature), Toposhape Xử lý phái sinh từ DEM trên GIS Phân tích mức độ tích tụ, phân tán dòng chảy mặt
Chỉ số TWI, Chỉ số SPI Tính toán thủy văn từ DEM Đánh giá độ ẩm bão hòa và sức cắt xói mòn của dòng nước
Mật độ sông suối (Stream Density) Bản đồ thủy văn và mạng lưới thủy đồ Mức độ tập trung nước và thoát lũ lưu vực
Chỉ số thực vật (NDVI) Xử lý từ ảnh Landsat / Sentinel-2 Khả năng giữ nước, cản trở dòng chảy của thảm phủ
Lượng mưa (Rainfall) Số liệu trạm khí tượng và trích xuất không gian Tác nhân kích hoạt trực tiếp hình thành lưu lượng đỉnh lũ
Loại đất (Soil), Thạch học (Lithology) Bản đồ thổ nhưỡng, bản đồ địa chất Khả năng ngấm, mức độ phong hóa và cung cấp bở rời cho bùn đá

Nội dung chính theo từng chương

Mở đầu

Phần Mở đầu thiết lập tính cấp thiết của đề tài trong bối cảnh các hiểm họa lũ quét liên tiếp gây thiệt hại nghiêm trọng tại vùng núi phía Bắc Việt Nam. Trình bày tổng quan mục tiêu nghiên cứu, đối tượng, phạm vi địa lý (huyện Bắc Hà và Bảo Yên, tỉnh Lào Cai), phạm vi khoa học, phương pháp luận nghiên cứu tiếp cận, 3 luận điểm bảo vệ, ý nghĩa khoa học và giá trị thực tiễn của luận án.

Chương 1: Tổng quan về nghiên cứu lũ quét

Chương 1 cung cấp bức tranh toàn diện về cơ sở lý thuyết và lịch sử nghiên cứu lũ quét:

  • Phân tích khung quy trình đánh giá rủi ro thảm họa thiên tai của UNISDR (3 chặng chính).
  • Chuẩn hóa các khái niệm, định nghĩa về lũ quét trên thế giới và tại Việt Nam.
  • Đánh giá ưu nhược điểm của 3 nhóm mô hình lũ quét trên thế giới: mô hình thống kê, mô hình mô phỏng mưa - dòng chảy (HEC-HMS, SWAT, MIKE) và mô hình phân loại "on-off" theo nhận dạng mẫu nhị phân.
  • Tổng kết tình hình nghiên cứu lũ quét tại Việt Nam từ các giai đoạn điều tra thống kê truyền thống, mô hình đa chỉ tiêu AHP, mô hình chuyên gia, đến các ứng dụng bước đầu của mạng trí tuệ nhân tạo.
  • Xác lập 2 điểm mới phát triển của luận án: (1) Ứng dụng ảnh Radar Sentinel-1A đa thời gian để lập bản đồ hiện trạng lũ quét; (2) Phát triển các mô hình học máy tích hợp tối ưu hóa metaheuristic mới (FA-LM-ANN, PSO-ELM, FURIA-GA-Ensemble).

Chương 2: Khu vực nghiên cứu và cơ sở khoa học

Chương 2 mô tả chi tiết điều kiện tự nhiên - kinh tế - xã hội của vùng nghiên cứu và luận giải cơ sở toán học của các thuật toán:

  • Khu vực nghiên cứu: Phân tích đặc điểm địa hình hiểm trở, hệ thống sông ngòi (sông Hồng, sông Chảy), chế độ khí hậu nhiệt đới gió mùa vùng núi (lượng mưa năm tại trạm Bắc Hà đạt 1.843,7 mm, tập trung hơn 80% vào mùa mưa từ tháng 4 đến tháng 9), dân số 145.208 người với tập quán canh tác sườn dốc dễ bị tổn thương.
  • Quy trình chiết xuất hiện trạng lũ từ ảnh Sentinel-1A: Cơ sở vật lý của độ tán xạ ngược sóng radar C-band phân cực VV đối với bề mặt nước và đất ướt sau lũ; các bước tiền xử lý ảnh trên SNAP (lọc nhiệt, hiệu chuẩn Sigma nought, lọc Lee, nắn hình học SRTM DEM).
  • Cơ sở toán học của mô hình FA-LM-ANN: Thiết lập kiến trúc mạng nơ-ron 3 lớp, công thức hàm biến đổi, cơ chế tối ưu hàm mất mát MSE bằng thuật toán Levenberg–Marquardt qua ma trận Jacobian, và thuật toán tìm kiếm vị trí tối ưu của cá thể đom đóm theo cường độ ánh sáng $I(r)$ và khoảng cách không gian đa chiều $r_{ij}$.
  • Cơ sở toán học của mô hình PSO-ELM: Thiết lập phương trình mạng ELM với ma trận nghịch đảo Moore–Penrose tính trọng số lớp ra; phương trình cập nhật vận tốc $V_i$ và vị trí $X_i$ của bầy hạt PSO để tối ưu hóa trọng số đầu vào $IW$ và độ lệch $BA$.
  • Cơ sở toán học của mô hình Ensemble Learning: Quy tắc luật mờ của thuật toán FURIA, toán tử mờ T-Norm, cơ chế chọn lọc - đột biến - lai ghép của thuật toán di truyền GA, và nguyên lý tạo tập con Bootstrap của kỹ thuật Bagging, kết hợp AdaBoost và LogitBoost trên nền cây quyết định C4.5.

Chương 3: Thực nghiệm và kết quả

Chương 3 trình bày toàn bộ quy trình thực nghiệm, xử lý dữ liệu và kết quả đạt được:

  • Xây dựng CSDL và bản đồ hiện trạng: Chiết xuất thành công các vị trí vết lũ từ ảnh Sentinel-1A và đối soát với dữ liệu kiểm chứng GPS thực địa; xây dựng 12 lớp bản đồ thành phần không gian dạng raster.
  • Phân tích đa cộng tuyến: Xác định các chỉ số TOL và VIF cho 12 bản đồ thành phần, chứng minh không có hiện tượng đa cộng tuyến nghiêm trọng giữa các biến đầu vào, đảm bảo tính độc lập để đưa vào mô hình hóa.
  • Thực nghiệm mô hình FA-LM-ANN:
    • Phân tích quá trình hội tụ tối ưu của thuật toán lai ghép FA-LM so với LM-ANN truyền thống.
    • Đánh giá hiệu suất dự báo của mô hình FA-LM-ANN so với các mô hình LM-ANN, FA-ANN, SVM và cây quyết định (CT).
    • Thành lập bản đồ phân vùng nguy cơ lũ quét tương ứng cho khu vực nghiên cứu.
  • Thực nghiệm mô hình PSO-ELM:
    • Đánh giá sự biến thiên giá trị tối ưu theo số vòng lặp trong quá trình huấn luyện; độ khớp của PSO-ELM trên tập huấn luyện và hiệu suất dự báo trên tập kiểm tra.
    • So sánh diện tích dưới đường cong AUC của PSO-ELM với các thuật toán cơ sở; thành lập bản đồ phân vùng nguy cơ lũ quét từ mô hình PSO-ELM và phân tích đặc trưng các lớp nguy cơ.
  • Thực nghiệm mô hình Ensemble Learning:
    • Đánh giá tầm quan trọng của các bản đồ thành phần theo phương pháp tiếp cận học máy.
    • Huấn luyện và kiểm định 3 mô hình tích hợp: FURIA-GA-LogitBoost, FURIA-GA-AdaBoost, FURIA-GA-Bagging.
    • So sánh độ chính xác phân loại, chỉ số ROC-AUC, Recall, Precision, tỷ lệ phần trăm mật độ tương đối lũ quét trên tập huấn luyện và tập kiểm tra.
    • Thực hiện kiểm nghiệm thống kê phi tham số Wilcoxon trên tập dữ liệu kiểm tra để xác nhận tính vượt trội và khác biệt có ý nghĩa thống kê của các mô hình đề xuất; xuất bản đồ phân vùng nguy cơ lũ quét tương ứng.

Kết luận và kiến nghị

Tổng kết các kết quả khoa học chính đã đạt được, khẳng định tính đúng đắn của các luận điểm bảo vệ; kiến nghị áp dụng bản đồ phân vùng nguy cơ lũ quét vào quy hoạch sử dụng đất, xây dựng công trình phòng chống thiên tai và phát triển các hệ thống cảnh báo sớm tại địa phương.


Kết quả và những đóng góp mới

Luận điểm bảo vệ của luận án

  1. Luận điểm 1: Sử dụng ảnh Sentinel-1 đa thời gian và kỹ thuật phát hiện thay đổi (Change Detection) cho phép xác định chính xác các khu vực bị ảnh hưởng bởi lũ quét ngoài thực địa mà không phụ thuộc vào điều kiện thời tiết hay mây mù.
  2. Luận điểm 2: Việc sử dụng tổ hợp 12 bản đồ thành phần (độ cao, độ dốc, hướng dốc, độ cong địa hình, hình thái địa mạo, TWI, SPI, NDVI, lượng mưa, loại đất, thạch học, mật độ sông suối) từ các nguồn dữ liệu địa không gian đa dạng cho phép thực hiện mô hình hóa và dự báo tin cậy các khu vực có nguy cơ xảy ra lũ quét.
  3. Luận điểm 3: Sự tích hợp của các thuật toán máy học và tối ưu hóa metaheuristic (FA-LM-ANN, PSO-ELM, FURIA-GA-Ensemble) cho phép xây dựng mô hình phân vùng nguy cơ lũ quét với độ chính xác cao vượt trội so với các mô hình đơn lẻ truyền thống.

Những đóng góp mới về mặt khoa học và thực tiễn

  • Về mặt khoa học:
    • Thiết lập thành công khung phương pháp luận khai thác ảnh vệ tinh Radar SAR Sentinel-1A (phân cực VV, chế độ IW) để tự động phát hiện vết lũ và xây dựng bản đồ hiện trạng lũ quét tại các vùng núi nhiệt đới có mây phủ dày đặc.
    • Phát triển mô hình tích hợp mới FA-LM-ANN, kết hợp thuật toán tối ưu bầy đom đóm (FA) tìm kiếm toàn cục với thuật toán Levenberg–Marquardt (LM) tìm kiếm nội bộ, giúp tự động tìm kiếm, cập nhật và tối ưu hóa trọng số của mạng nơ-ron nhân tạo cho bài toán dự báo lũ quét.
    • Phát triển mô hình tích hợp mới PSO-ELM, kết hợp giải thuật tối ưu hóa bầy đàn (PSO) với máy học cực đại (ELM), giải quyết triệt để nhược điểm gán trọng số ngẫu nhiên ban đầu của ELM.
    • Phát triển hệ mô hình Ensemble Learning mới kết hợp thuật toán di truyền (GA), hệ luật mờ không thứ tự (FURIA) với các kỹ thuật tổ hợp Bagging, AdaBoost và LogitBoost trên nền cây quyết định C4.5.
  • Về mặt thực tiễn:
    • Cung cấp bộ bản đồ phân vùng nguy cơ lũ quét chi tiết với độ chính xác cao cho hai huyện trọng điểm Bắc Hà và Bảo Yên (tỉnh Lào Cai).
    • Cung cấp một tài liệu hướng dẫn quy trình công nghệ hoàn chỉnh từ xử lý ảnh viễn thám, tích hợp GIS đến mô hình hóa AI, có khả năng chuyển giao và áp dụng tương tự cho các vùng miền núi khác tại Việt Nam mà không yêu cầu số liệu trạm quan trắc thủy văn dày đặc.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế

  • Mô hình tiếp cận dựa trên giả thuyết phân loại nhị phân không gian ("on-off") tập trung vào việc xác định vị trí và xác suất không gian có nguy cơ xảy ra lũ quét, chưa giải quyết được bài toán dự báo thời điểm chính xác khi nào lũ quét sẽ xảy ra trong tương lai theo thời gian thực.
  • Dữ liệu thực nghiệm của luận án mới chỉ tập trung trên phạm vi địa bàn hai huyện Bắc Hà và Bảo Yên thuộc tỉnh Lào Cai do các giới hạn về nguồn lực khảo sát thực địa và kinh phí nghiên cứu.

Hướng nghiên cứu tiếp

  • Mở rộng phạm vi thử nghiệm mô hình hóa sang các lưu vực sông suối và các tỉnh miền núi khác có điều kiện địa chất, địa mạo và khí hậu khác biệt (như khu vực Tây Bắc, miền Trung và Tây Nguyên).
  • Nghiên cứu tích hợp chuỗi số liệu dự báo mưa thời gian thực độ phân giải cao và các mô hình viễn thám khí tượng nhằm nâng cấp mô hình phân vùng tĩnh thành hệ thống cảnh báo sớm lũ quét động theo cả không gian và thời gian.

Giá trị tham khảo

Luận án là tài liệu chuyên khảo có giá trị học thuật và ứng dụng thực tiễn cao cho nhiều nhóm đối tượng:

  • Nghiên cứu sinh, học viên cao học và giảng viên: Thuộc các chuyên ngành Kỹ thuật Trắc địa - Bản đồ, Viễn thám và GIS, Địa lý tự nhiên, Thủy văn, Quản lý tài nguyên môi trường và Trí tuệ nhân tạo ứng dụng. Luận án cung cấp chi tiết cơ sở toán học, phương trình giải tích và quy trình tích hợp các thuật toán Metaheuristic (FA, PSO, GA) vào mạng nơ-ron và máy học.
  • Cán bộ quản lý và cơ quan hoạch định chính sách: Ban chỉ huy Phòng thủ dân sự, Phòng chống thiên tai và Tìm kiếm cứu nạn các cấp, Sở Tài nguyên và Môi trường, Sở Nông nghiệp và Phát triển nông thôn các tỉnh miền núi. Kết quả phân vùng nguy cơ là căn cứ khoa học trực tiếp phục vụ quy hoạch bố trí dân cư an toàn, lập kế hoạch sử dụng đất và xây dựng các công trình giảm thiểu rủi ro thiên tai.
  • Các phần đáng tham khảo nhất:
    1. Quy trình tiền xử lý ảnh Radar SAR Sentinel-1A trên phần mềm SNAP Toolbox để chiết xuất hiện trạng lũ (Mục 2.2).
    2. Phương pháp tích hợp thuật toán lai ghép FA-LM cho mạng ANN và PSO cho mậng ELM (Mục 2.3 và 2.4).
    3. Bộ tiêu chí 12 bản đồ thành phần không gian và kỹ thuật kiểm tra đa cộng tuyến TOL - VIF phục vụ mô hình hóa hiểm họa lũ quét (Chương 3).

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

1. Khu vực nghiên cứu thực nghiệm của luận án có những đặc điểm tự nhiên nào khiến lũ quét thường xuyên xảy ra?

Khu vực nghiên cứu gồm hai huyện Bắc Hà và Bảo Yên (tỉnh Lào Cai) có diện tích 1.510,4 km², độ cao biến thiên từ 38,9 m đến 1.878,69 m (trung bình 538,1 m). Địa hình đặc trưng miền núi hiểm trở với 85,4% diện tích có độ dốc từ 10° đến 40°. Vùng có hai hệ thống sông lớn chảy qua là sông Hồng (dài 28,7 km qua khu vực) và sông Chảy (dài 91,6 km, độ dốc lớn, dòng chảy xiết). Chế độ mưa tập trung cao độ từ tháng 4 đến tháng 9, chiếm hơn 80% tổng lượng mưa năm (lượng mưa năm tại trạm Bắc Hà đạt 1.843,7 mm), trong đó các tháng 6, 7, 8 chiếm trên 50% lượng mưa cả năm. Địa hình dốc kết hợp mưa lớn cục bộ là nguyên nhân chính hình thành các trận lũ quét dữ dội.

2. Vì sao luận án lựa chọn ảnh Radar SAR Sentinel-1A thay vì ảnh viễn thám quang học để thành lập bản đồ hiện trạng lũ quét?

Lũ quét tại miền núi thường xảy ra trong điều kiện thời tiết mưa bão cực đoan, bầu trời bị bao phủ bởi mây mù dày đặc, khiến ảnh viễn thám quang học (như Landsat hay Sentinel-2) bị che khuất và không thể ghi nhận được bề mặt đất. Ảnh Radar SAR Sentinel-1A sử dụng sóng viba C-band (bước sóng 3,5 cm, tần số 4–8 GHz) có khả năng xuyên mây và hoạt động không phụ thuộc vào điều kiện ánh sáng ngày/đêm. Với độ phân giải không gian 10 m, chu kỳ chụp lặp 12 ngày và định dạng phân cực kép VV, ảnh Sentinel-1A cho phép phát hiện rõ nét sự thay đổi hệ số tán xạ ngược của bề mặt đất ướt và vết lũ sau thiên tai.

3. Nguyên lý hoạt động của mô hình tích hợp FA-LM-ANN trong luận án là gì?

Mô hình FA-LM-ANN là sự kết hợp giữa mạng nơ-ron nhân tạo nhiều lớp (ANN) với hai thuật toán tối ưu:

  • Thuật toán bầy đom đóm (FA): Đóng vai trò là công cụ tìm kiếm và tối ưu hóa toàn cục, dựa trên quy luật tương tác và độ sáng giữa các cá thể đom đóm trong không gian đa chiều để xác định vùng trọng số tiềm năng ban đầu, giúp mạng tránh rơi vào bẫy cực tiểu cục bộ.
  • Thuật toán Levenberg–Marquardt (LM): Đóng vai trò là công cụ tìm kiếm nội bộ cục bộ, sử dụng ma trận sai số Jacobian để cập nhật và hiệu chỉnh trọng số với tốc độ hội tụ nhanh và độ ổn định cao.

4. Thuật toán PSO đã khắc phục nhược điểm cốt lõi nào của máy học cực đại (ELM)?

Mô hình máy học cực đại ELM truyền thống có ưu điểm là tốc độ huấn luyện rất nhanh do trọng số đầu ra được tính toán bằng phương pháp giải tích thông qua ma trận nghịch đảo tổng quát Moore–Penrose. Tuy nhiên, nhược điểm lớn của ELM là các vector trọng số đầu vào và độ lệch (bias) của lớp ẩn lại được khởi tạo hoàn toàn ngẫu nhiên, dẫn đến việc không tối ưu và có thể làm giảm hiệu suất dự báo trong nhiều trường hợp. Thuật toán tối ưu hóa bầy đàn (PSO) được tích hợp để mô phỏng trí thông minh tập thể của bầy chim, tự động tìm kiếm và tối ưu hóa toàn bộ các trọng số đầu vào và độ lệch trước khi tính toán trọng số đầu ra.

5. Những chỉ số thống kê nào được luận án sử dụng để kiểm chứng và so sánh độ chính xác của các mô hình?

Luận án sử dụng một hệ thống chỉ số thống kê toàn diện gồm:

  • Chỉ số sai số số học: Căn bậc hai sai số toàn phương trung bình (RMSE), Sai số bình phương trung bình (MSE), Sai số tuyệt đối trung bình (MAE), Sai số trung bình (ME).
  • Chỉ số phân loại: Ma trận nhầm lẫn (TP, TN, FP, FN), Tỷ lệ phân loại chính xác tổng thể, Hệ số tương đồng Kappa.
  • Đường cong và diện tích: Đường cong đặc tính máy thu (ROC) và Diện tích dưới đường cong ROC (AUC) trên cả tập huấn luyện (70%) và tập kiểm tra (30%), kết hợp phương pháp kiểm tra chéo 10-fold và 5-fold cross-validation.
  • Kiểm định phi tham số: Kiểm nghiệm Wilcoxon signed-rank test để chứng minh sự khác biệt về mặt thống kê giữa mô hình đề xuất và các mô hình cơ sở.

Kết luận

Luận án tiến sĩ của NCS Ngô Thị Phương Thảo đã giải quyết một bài toán khoa học và thực tiễn có ý nghĩa cấp bách trong công tác phòng chống thiên tai tại Việt Nam thông qua việc phát triển các mô hình trí tuệ nhân tạo tích hợp mới (FA-LM-ANN, PSO-ELM, FURIA-GA-Ensemble). Công trình đã thiết lập thành công quy trình công nghệ hoàn chỉnh từ khâu khai thác ảnh vệ tinh Radar SAR Sentinel-1A đa thời gian để thành lập bản đồ hiện trạng vết lũ, xử lý 12 lớp dữ liệu địa không gian GIS, đến khâu mô hình hóa học máy tối ưu hóa metaheuristic đạt độ chính xác cao. Kết quả phân vùng nguy cơ lũ quét tại huyện Bắc Hà và Bảo Yên (Lào Cai) là cơ sở khoa học tin cậy phục vụ quy hoạch không gian, quản lý rủi ro thiên tai và bảo vệ an toàn dân sinh ở các vùng miền núi.