Tổng quan về luận án

Luận án này tiên phong trong việc định hình lại cách chúng ta đặc trưng hóa bề mặt địa hình và đánh giá nguy cơ lũ lụt tại các vùng đồng bằng ven biển. Nghiên cứu này được thực hiện trong bối cảnh các vùng ven biển miền Trung Việt Nam, đặc biệt là Đà Nẵng, thường xuyên phải đối mặt với các sự kiện lũ lụt nghiêm trọng, gây ra những thiệt hại đáng kể về người và của. Các phương pháp quản lý lũ lụt hiện hành và các nghiên cứu trước đây (như Ho et al., 2012; Do et al., 2014) đã bộc lộ những hạn chế cố hữu, chủ yếu do thiếu mô hình độ cao số (DEM) có độ phân giải cao và dữ liệu vệ tinh về các sự kiện gần lũ. Điều này đã dẫn đến việc tạo ra các kết quả thiếu chính xác, làm suy yếu hiệu quả của các kế hoạch phòng chống lũ lụt toàn diện.

Khoảng trống nghiên cứu (Research Gap) cụ thể mà luận án này giải quyết bao gồm:

  1. Thiếu DEM có độ phân giải cao và chất lượng cao: Các nghiên cứu trước đây tại Đà Nẵng đã gặp phải vấn đề "lack of high resolution digital elevation model (DEM)" (Abstract, p. 9) và "lack of high resolution as well as high quality DEM which relatively affects the accuracy of flood hazard models" (Chapter 1, p. 12), dẫn đến kết quả đánh giá nguy cơ lũ không chính xác. Luận án này giải quyết trực tiếp vấn đề này bằng cách phát triển các thuật toán tạo DEM chất lượng cao, bao gồm cả kỹ thuật hợp nhất DEM toàn cầu và phương pháp BS-Horizon độ phân giải cao.
  2. Thiếu dữ liệu vệ tinh về các sự kiện gần lũ: Các nghiên cứu trước đây cũng bị hạn chế bởi "near-flood events satellite data which leads to generate inaccuracy results" (Abstract, p. 9). Luận án này khắc phục bằng cách tích hợp dữ liệu vệ tinh ALOS PALSAR năm 2007, một sự kiện lũ lụt lớn, để lập bản đồ ngập lụt thực tế.
  3. Đánh giá chưa đầy đủ về tham số của phương pháp BS-Horizon: Phương pháp tạo bề mặt BS-Horizon (Nonogaki et al., 2012) chưa từng được đánh giá toàn diện về ảnh hưởng của các thiết lập tham số. "This study is the first time that the parameters involved in BS-Horizon method have been fully evaluated" (Abstract, p. 10), mang lại sự hiểu biết sâu sắc và hướng dẫn thực hành cho việc sử dụng phương pháp này.

Các câu hỏi nghiên cứu và giả thuyết (Research Questions and Hypotheses) chính của luận án bao gồm:

  1. RQ1: Làm thế nào để hợp nhất các DEM toàn cầu miễn phí (GDEM và SRTM) một cách hiệu quả để nâng cao độ chính xác của dữ liệu độ cao cho các nghiên cứu quy mô khu vực, đặc biệt trong các bối cảnh địa hình đa dạng?
    • H1: Việc áp dụng một phương pháp hợp nhất DEM dựa trên đặc điểm địa hình (geomorphological approach) với trọng số khác nhau cho GDEM và SRTM tùy thuộc vào bối cảnh địa hình sẽ cải thiện đáng kể độ chính xác của DEM so với từng DEM toàn cầu riêng lẻ.
  2. RQ2: Làm thế nào để tạo ra một DEM độ phân giải cao (5m) cho vùng đồng bằng ven biển dễ bị lũ lụt bằng cách sử dụng phương pháp BS-Horizon, tích hợp cả ràng buộc đẳng thức và bất đẳng thức, và đâu là các thiết lập tham số tối ưu cho phương pháp này?
    • H2: Việc đánh giá có hệ thống các tham số (Mx, My, α, m1, m2) của phương pháp BS-Horizon sẽ cho phép xác định các thiết lập tối ưu, giúp tạo ra DEM độ phân giải cao với độ chính xác vượt trội (RMSE < 0.25m) so với các DEM có sẵn.
  3. RQ3: Một mô hình phân vùng nguy cơ lũ lụt đa tham số sử dụng Quy trình Phân cấp Phân tích (AHP), tích hợp các thông số địa hình và thông tin ngập lụt từ vệ tinh, có thể đánh giá chính xác nguy cơ lũ lụt ở vùng đồng bằng ven biển Đà Nẵng hay không?
    • H3: Việc kết hợp năm yếu tố (ngập lụt dựa trên độ cao (EFI), khoảng cách từ kênh sông (DIST), sử dụng đất, độ dốc và chỉ số độ ẩm địa hình (TWI)) thông qua AHP sẽ tạo ra bản đồ phân vùng nguy cơ lũ lụt có độ tương ứng cao (trên 80%) với các khu vực bị ngập lụt được quan sát từ dữ liệu vệ tinh và các điểm trụ lũ thực địa.

Khung lý thuyết (Theoretical Framework): Nghiên cứu này được xây dựng trên nền tảng của các lý thuyết và phương pháp từ Geoinformatics, Thủy văn học địa hình (Topographic Hydrology)Khoa học Quyết định (Decision Science).

  • Nghiên cứu mở rộng lý thuyết về tạo bề mặt địa hình (surface generation) thông qua việc tối ưu hóa phương pháp BS-Horizon (Nonogaki et al., 2012), một sự tích hợp của hàm bi-cubic splinehàm phạt ngoại vi (exterior penalty function), đặc biệt là trong việc sử dụng kết hợp các ràng buộc đẳng thức và bất đẳng thức.
  • Nó cũng dựa trên các nguyên tắc của viễn thám (remote sensing) để thu thập và xử lý dữ liệu độ cao (ASTER GDEM, SRTM) và dữ liệu ngập lụt (ALOS PALSAR).
  • Lý thuyết Hệ thống Phân cấp Phân tích (Analytical Hierarchy Process - AHP) của Saaty cung cấp khuôn khổ để định lượng hóa và tổng hợp các yếu tố gây lũ tiềm năng.

Đóng góp đột phá với tác động định lượng:

  1. Phương pháp hợp nhất DEM dựa trên địa hình tiên tiến: Luận án này đề xuất một phương pháp hợp nhất DEM dựa trên bối cảnh địa hình độc đáo, đạt được độ chính xác đáng kể. RMSE của DEM hợp nhất giảm xuống còn 11.6m, cải thiện đáng kể so với RMSE của GDEM (14.9m) và SRTM (14.8m) riêng lẻ (Table 2.5, p. 29). Điều này mở ra khả năng sử dụng hiệu quả các DEM miễn phí trên quy mô lớn, đặc biệt ở những khu vực thiếu dữ liệu chất lượng cao.
  2. Tối ưu hóa tham số BS-Horizon lần đầu tiên: Luận án đã thực hiện đánh giá toàn diện các tham số của phương pháp BS-Horizon (M, α, m1, m2), điều chưa từng được thực hiện trước đây. Điều này cung cấp hướng dẫn thực tiễn và lý thuyết cho việc tạo DEM độ phân giải cao cho các vùng đồng bằng ven biển phức tạp, cho phép đạt được RMSE 0.17m (p. 47) cho DEM 5m.
  3. Tạo ra DEM độ phân giải cực cao cho vùng đất thấp: Việc tạo ra DEM 5m cho vùng đồng bằng ven biển Đà Nẵng với độ chính xác cao (RMSE 0.17m) đã khắc phục "lack of high resolution digital elevation model" (p. 9) trong các nghiên cứu trước. Điều này cải thiện đáng kể so với DEM hợp nhất (RMSE 1.98m) (p. 47) và là dữ liệu đầu vào đáng tin cậy cho phân tích địa hình và nguy cơ lũ lụt chi tiết.
  4. Mô hình phân vùng nguy cơ lũ lụt AHP đa chiều được xác thực chặt chẽ: Mô hình AHP được phát triển đã tích hợp thành công năm yếu tố gây lũ và được xác thực bằng dữ liệu vệ tinh ALOS PALSAR sự kiện lũ năm 2007. Kết quả cho thấy gần 84% các khu vực bị ngập lụt được quan sát thuộc vùng có nguy cơ lũ cao và rất cao (Abstract, p. 10; Chapter 4, p. 49), cho thấy độ tin cậy cao của phương pháp.

Phạm vi (Scope) và ý nghĩa (Significance): Nghiên cứu tập trung vào vùng đồng bằng ven biển Đà Nẵng, Việt Nam, một khu vực có diện tích 98 km2 với độ cao từ 0m đến 10m (Chapter 4, p. 50). Dữ liệu được sử dụng bao gồm 9.730 điểm độ cao khảo sát thực địa (với độ chính xác ±10cm) (Chapter 3, p. 35) và hình ảnh vệ tinh từ năm 2007 (Landsat TM, ALOS PALSAR) và 2014 (RapidEye). Khung thời gian dữ liệu kéo dài từ 2000 đến 2015. Ý nghĩa của luận án nằm ở khả năng cung cấp một bộ công cụ và phương pháp luận mạnh mẽ để đánh giá nguy cơ lũ lụt trong các vùng ven biển dễ bị tổn thương, không chỉ ở Việt Nam mà còn ở các khu vực tương tự trên toàn cầu. Các phương pháp tạo DEM tiên tiến và mô hình hóa nguy cơ lũ lụt được đề xuất có thể hỗ trợ các nhà hoạch định chính sách và quản lý thiên tai trong việc phát triển các kế hoạch phòng chống lũ lụt hiệu quả hơn, giảm thiểu thiệt hại kinh tế và xã hội.

Literature Review và Positioning

Luận án này được đặt trong dòng nghiên cứu về mô hình hóa bề mặt địa hình và đánh giá nguy cơ thiên tai, với sự tập trung cụ thể vào lũ lụt ở các vùng ven biển. Ba luồng nghiên cứu chính đã được tổng hợp để định vị luận án: (1) Đánh giá và hợp nhất DEM toàn cầu, (2) Phát triển DEM độ phân giải cao, và (3) Phân vùng nguy cơ lũ lụt sử dụng các kỹ thuật địa không gian.

Tổng hợp các luồng nghiên cứu chính với tên tác giả và năm cụ thể:

  1. Đánh giá và hợp nhất DEM toàn cầu: Nhiều nghiên cứu đã đánh giá độ chính xác của các DEM toàn cầu miễn phí như ASTER GDEM và SRTM. Li et al. (2013), Ravibabu et al. (2010), Zhao et al. (2011), Suwandana et al. (2012), Mukherjee et al. (2013a) và Czubski et al. (2013) đã thực hiện các đánh giá so sánh. Reuter et al. (2007), Mukherjee et al. (2013a), Czubski et al. (2013) và Fuss (2013) cũng đã đánh giá độ chính xác của các DEM toàn cầu dựa trên đặc điểm địa hình. Các nghiên cứu này chủ yếu tập trung vào việc xác minh chất lượng. Tuy nhiên, khả năng hợp nhất các DEM từ các nguồn khác nhau để tận dụng ưu điểm của từng loại đã được gợi ý bởi Kaab (2005). Các phương pháp hợp nhất trước đây (ví dụ: Kaab, 2005; Karkee et al., 2008; Papasaika et al., 2011; Lucca, 2011; Fuss, 2013) thường sử dụng phương pháp trung bình có trọng số nhưng không xem xét trọng số thay đổi theo bối cảnh địa hình cụ thể.

  2. Phát triển DEM độ phân giải cao: Các phương pháp nội suy truyền thống (ví dụ: kriging, spline, IDW) chỉ sử dụng các điểm dữ liệu thực tế hoặc đường đồng mức (equality constraints). Tuy nhiên, Shiono et al. (2001) đã phát triển thuật toán Horizon2000 để xác định bề mặt địa chất bằng cách kết hợp thông tin bất đẳng thức ("trên" hoặc "dưới") và độ dốc. Nonogaki et al. (2012) đã cải tiến phương pháp này thành BS-Horizon, tích hợp nội suy bi-cubic spline và hàm phạt ngoại vi để xử lý các ràng buộc đẳng thức-bất đẳng thức. Mặc dù các nghiên cứu như Shiono et al. (2001), Noumi (2003) và Nonogaki et al. (2012) đã chứng minh hiệu quả của phương pháp này, nhưng chúng chưa thảo luận về ảnh hưởng của các thiết lập tham số đến bề mặt địa hình thu được.

  3. Phân vùng nguy cơ lũ lụt sử dụng các kỹ thuật địa không gian: Nhiều phương pháp đã được phát triển để lập bản đồ nguy cơ lũ lụt, sử dụng các cách tiếp cận thủy văn, khí tượng và địa hình (ví dụ: Kenny, 1990; Ballais et al., 2005; Mafreda et al., 2013). Tại khu vực Đà Nẵng, Ho et al. (2012) và Do et al. (2014) đã cố gắng đặc trưng hóa tiềm năng lũ lụt bằng các phương pháp thủy văn - địa hình tích hợp dữ liệu viễn thám. Tuy nhiên, những nghiên cứu này đã bị hạn chế bởi "lack of high resolution digital elevation model (DEM) and near-flood events satellite data" (Chapter 1, p. 15), làm giảm độ chính xác của kế hoạch phòng chống lũ.

Mâu thuẫn/tranh luận (Contradictions/debates) với ít nhất 2 quan điểm đối lập: Một mâu thuẫn trung tâm trong lĩnh vực này là giữa sự sẵn có của các DEM toàn cầu miễn phí nhưng độ chính xác thấp và nhu cầu cấp thiết về DEM độ phân giải cao, chính xác cho các ứng dụng quy mô địa phương.

  • Quan điểm 1 (Đề xuất bởi nghiên cứu trước): Các DEM toàn cầu như GDEM và SRTM, mặc dù có độ bao phủ rộng, chứa các lỗi độ cao đáng kể do đặc điểm thu thập dữ liệu (quang học đối với GDEM, InSAR đối với SRTM) và bối cảnh địa hình (Karkee et al., 2008; Reuter et al., 2007). Các lỗi này bao gồm "cloud coverage, radiometric variation and low levels of texture" đối với GDEM và "shadowing, layovers or complex dielectric constant" đối với SRTM (Chapter 2, p. 26). Do đó, chúng "may not be sufficient for flood prevention plan" (Chapter 1, p. 15) ở quy mô địa phương.
  • Quan điểm 2 (Hướng giải quyết của luận án): Luận án này lập luận rằng sự khác biệt về đặc điểm và lỗi giữa GDEM và SRTM thực sự tạo ra "an opportunity for DEM fusion" (Kaab, 2005, trích dẫn trong Chapter 1, p. 14). Bằng cách hiểu rõ những hạn chế này trong các bối cảnh địa hình khác nhau, có thể phát triển các phương pháp hợp nhất thông minh hơn để nâng cao chất lượng của các DEM toàn cầu. Ví dụ, trong khi GDEM có thể kém chính xác ở vùng bằng phẳng, SRTM lại có thể kém hơn ở các thung lũng dốc (Figure 2.5, p. 23). Luận án này tận dụng các đặc điểm bổ trợ này.

Định vị trong tài liệu (Positioning in literature) với khoảng trống cụ thể được xác định: Luận án định vị mình là một đóng góp quan trọng bằng cách giải quyết các khoảng trống đã được xác định:

  1. Trong hợp nhất DEM: Luận án tiên phong trong việc đề xuất một "geomorphological approach for DEM fusion" (Chapter 1, p. 14), nơi các trọng số cho GDEM và SRTM được xác định dựa trên bối cảnh địa hình và loại hình sử dụng đất. Cách tiếp cận này "have not been attempted in previous studies on DEM fusion" (Chapter 1, p. 14), vượt qua giới hạn của các phương pháp trung bình có trọng số truyền thống chỉ xem xét độ chính xác tổng thể.
  2. Trong tạo DEM độ phân giải cao: Luận án là nghiên cứu đầu tiên "fully evaluated" (Abstract, p. 10) các tham số của phương pháp BS-Horizon (Nonogaki et al., 2012). Điều này bổ sung một cơ sở lý thuyết và thực tiễn vững chắc cho một phương pháp nội suy bề mặt mạnh mẽ, đặc biệt quan trọng cho các vùng đất thấp nơi các kỹ thuật viễn thám thông thường thất bại.
  3. Trong phân vùng nguy cơ lũ lụt: Luận án cải thiện đáng kể các nghiên cứu trước đây tại Đà Nẵng (Ho et al., 2012; Do et al., 2014) bằng cách sử dụng DEM độ phân giải cao 5m được tạo ra có mục đích và dữ liệu vệ tinh về các sự kiện gần lũ (ALOS PALSAR 2007). Điều này dẫn đến một bản đồ phân vùng nguy cơ lũ lụt "good correlation to field survey flood pillar data and flood inundated areas extracted from ALOS PALSAR" (Chapter 1, p. 12).

Cách nghiên cứu này thúc đẩy lĩnh vực (How this advances field) với các đóng góp cụ thể: Nghiên cứu này thúc đẩy lĩnh vực Geoinformatics và Quản lý thiên tai bằng cách cung cấp các đóng góp cụ thể:

  • Cải thiện đáng kể độ chính xác của DEM toàn cầu: Việc giảm RMSE của DEM hợp nhất từ 14.9m/14.8m xuống 11.6m (Table 2.5, p. 29) cho thấy tiềm năng thực tế của phương pháp hợp nhất dựa trên địa hình, làm cho các DEM miễn phí trở nên hữu ích hơn cho các ứng dụng quy mô khu vực.
  • Tiêu chuẩn hóa việc sử dụng BS-Horizon: Việc đánh giá có hệ thống các tham số của BS-Horizon, xác định các thiết lập tối ưu như M=200, α=1.58x10^4 (cho dữ liệu đẳng thức-bất đẳng thức) và m1=0.5 (Chapter 3, p. 46), cung cấp một giao thức rõ ràng cho việc tạo DEM độ phân giải cao với RMSE 0.17m (p. 47). Điều này nâng cao độ tin cậy và khả năng tái lập của các nghiên cứu tạo bề mặt địa hình.
  • Mô hình hóa nguy cơ lũ lụt chính xác hơn: Sự kết hợp của DEM 5m độ phân giải cao và phương pháp AHP đa tham số, được xác thực bằng dữ liệu ngập lụt thực tế, mang lại một mô hình phân vùng nguy cơ lũ lụt với "nearly 84 percent of inundated areas observed in ALOS PALSAR data fall in the high and very high flood hazard zones" (Chapter 1, p. 16), vượt trội về độ chính xác so với các nghiên cứu trước đây.

So sánh với ít nhất 2 nghiên cứu quốc tế:

  1. So sánh với các nghiên cứu về hợp nhất DEM (ví dụ: Karkee et al., 2008; Lucca, 2011): Trong khi Karkee et al. (2008) sử dụng biến đổi Fourier nhanh (FFT) và lọc miền tần số, và Lucca (2011) so sánh trung bình có trọng số với dự đoán đối chiếu, luận án này khác biệt ở chỗ nó sử dụng "weighted averaging... based on geomorphologic characteristics" (Chapter 2, p. 26). Luận án xác định trọng số cho GDEM và SRTM dựa trên độ chính xác của chúng trong các loại hình địa hình khác nhau (núi, thung lũng, bằng phẳng) và loại hình sử dụng đất (Table 2.3, Table 2.4). Điều này dẫn đến sự cải thiện đáng kể về độ chính xác (RMSE 11.6m), vượt trội so với các phương pháp hợp nhất tổng quát hơn.
  2. So sánh với các nghiên cứu sử dụng BS-Horizon trước đây (ví dụ: Nonogaki et al., 2012): Nonogaki et al. (2012) giới thiệu phương pháp BS-Horizon và chứng minh khả năng tạo bề mặt với ràng buộc đẳng thức-bất đẳng thức. Tuy nhiên, nghiên cứu này không "discussed about effect of parameter settings on the resultant topographic surfaces" (Chapter 1, p. 15). Luận án hiện tại lấp đầy khoảng trống này bằng cách cung cấp một "comprehensive investigation of parameter settings for BS-Horizon method" (Chapter 3, p. 32), cung cấp hướng dẫn rõ ràng về cách tối ưu hóa các tham số M, α, m1, m2 để đạt được bề mặt có độ chính xác cao nhất (RMSE 0.17m) và độ mịn hợp lý.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án này đóng góp đáng kể vào cơ thể lý thuyết về Geoinformatics và Thủy văn học địa hình thông qua việc mở rộng và tinh chỉnh các lý thuyết hiện có.

  • Mở rộng lý thuyết về Tạo bề mặt Địa hình (Surface Generation): Luận án mở rộng lý thuyết về nội suy bề mặt, đặc biệt là phương pháp BS-Horizon của Nonogaki et al. (2012). Trong khi Nonogaki và cộng sự đã giới thiệu khả năng của hàm bi-cubic spline và hàm phạt ngoại vi để kết hợp ràng buộc đẳng thức và bất đẳng thức, luận án này là nghiên cứu đầu tiên "fully evaluated" (Abstract, p. 10) ảnh hưởng của các tham số cài đặt (Mx, My, α, m1, m2) đến bề mặt địa hình được tạo ra. Bằng cách định lượng hóa cách các tham số này ảnh hưởng đến độ mịn bề mặt (J(f)) và mức độ phù hợp (R(f)), luận án cung cấp một hiểu biết sâu sắc hơn về hành vi của mô hình và đưa ra các nguyên tắc tối ưu hóa thực tế. Điều này tinh chỉnh và làm cho lý thuyết BS-Horizon trở nên mạnh mẽ hơn, đặc biệt cho các ứng dụng thực tiễn trong việc tạo DEM độ phân giải cao cho các vùng địa hình phức tạp.
  • Đề xuất Khung Phân tích Hợp nhất DEM Đa tiêu chí: Luận án thách thức giả định về trọng số cố định trong các phương pháp hợp nhất DEM truyền thống bằng cách đề xuất một cách tiếp cận mới dựa trên bối cảnh địa hình. Thay vì xem xét các DEM toàn cầu như một thực thể đồng nhất, luận án này chứng minh rằng "accuracy of GDEM and SRTM varies depending upon the geomorphological characteristics of target area" (Abstract, p. 9). Điều này gợi ý một sự mở rộng của lý thuyết về hợp nhất dữ liệu, nơi độ tin cậy của các nguồn dữ liệu khác nhau nên được đánh giá động dựa trên các thuộc tính nội tại của đối tượng được mô hình hóa. Luận án đã xác định rằng SRTM tốt hơn ở vùng bằng phẳng và GDEM có thể tốt hơn ở các thung lũng dốc (Figure 2.5, p. 23), từ đó phát triển một phương pháp trung bình có trọng số linh hoạt.
  • Khung Lý thuyết về Mô hình hóa Nguy cơ Lũ lụt Đa tham số: Bằng cách tích hợp AHP với năm yếu tố gây lũ được dẫn xuất từ DEM độ phân giải cao và dữ liệu viễn thám sự kiện gần lũ, luận án củng cố khung lý thuyết về mô hình hóa nguy cơ lũ lụt. Nó cung cấp một cách tiếp cận có hệ thống để định lượng hóa tác động tương đối của các yếu tố địa hình và thủy văn đối với khả năng xảy ra lũ lụt.
  • Phát hiện hiện tượng mới: Luận án cung cấp bằng chứng thực nghiệm rằng việc tích hợp ràng buộc bất đẳng thức trong nội suy bi-cubic spline (phương pháp BS-Horizon) có thể kiểm soát hiệu quả các "strong oscillations" và "significant undershoot elevation" (p. 45) thường gặp trong các thuật toán spline truyền thống, đặc biệt ở những khu vực thiếu dữ liệu điểm hoặc có sự chuyển đổi rõ rệt như sông và đường. Điều này mở rộng lý thuyết về nội suy bề mặt bằng cách cung cấp một giải pháp hiệu quả cho một vấn đề đã biết.

Khung khái niệm (Conceptual Framework) và mô hình lý thuyết (Theoretical model): Khung khái niệm của luận án được minh họa thông qua chuỗi các quy trình từ tạo DEM đến phân vùng nguy cơ lũ lụt.

  • Giai đoạn 1: Nâng cao chất lượng DEM toàn cầu. Các lý thuyết về viễn thám quang học (cho GDEM) và InSAR (cho SRTM) được kết hợp với các lý thuyết về đặc trưng địa hình để đề xuất một phương pháp hợp nhất có trọng số theo địa hình.
  • Giai đoạn 2: Tạo DEM độ phân giải cao cục bộ. Lý thuyết về nội suy bi-cubic spline (Nonogaki et al., 2012) được tinh chỉnh thông qua việc đánh giá tham số toàn diện và tích hợp ràng buộc đẳng thức/bất đẳng thức, tạo ra một mô hình bề mặt địa hình chi tiết cho vùng đất thấp.
  • Giai đoạn 3: Phân vùng nguy cơ lũ lụt. Mô hình lý thuyết được thiết lập dựa trên năm tiền đề chính:
    1. P1: Độ cao dựa trên ngập lụt (EFI) có mối quan hệ nghịch biến với nguy cơ lũ lụt.
    2. P2: Khoảng cách từ kênh sông (DIST) có mối quan hệ nghịch biến với nguy cơ lũ lụt.
    3. P3: Loại hình sử dụng đất (Land Use) ảnh hưởng đến khả năng thấm nước và thoát nước, từ đó ảnh hưởng đến nguy cơ lũ lụt.
    4. P4: Độ dốc (Slope) có mối quan hệ nghịch biến với nguy cơ lũ lụt (độ dốc thấp tăng nguy cơ).
    5. P5: Chỉ số độ ẩm địa hình (TWI) có mối quan hệ đồng biến với nguy cơ lũ lụt (TWI cao tăng nguy cơ). Các tiền đề này được tích hợp vào một mô hình phân tích quyết định đa tiêu chí sử dụng AHP, trong đó trọng số được gán thông qua so sánh cặp theo thang Saaty (Table 4.2, p. 52).

Chuyển đổi mô hình (Paradigm Shift) với bằng chứng từ các phát hiện: Luận án này góp phần vào một sự chuyển đổi mô hình nhỏ từ cách tiếp cận chung, một cỡ cho tất cả (one-size-fits-all) trong việc sử dụng dữ liệu độ cao toàn cầu sang một cách tiếp cận thích ứng, bối cảnh cụ thể. Bằng chứng là phát hiện "the accuracy of GDEM and SRTM varies depending upon the geomorphological characteristics of target area" (Abstract, p. 9) và việc đề xuất một phương pháp hợp nhất động thay vì tĩnh. Hơn nữa, việc tối ưu hóa BS-Horizon cho các vùng đất thấp khó mô hình hóa truyền thống, với khả năng đạt được độ chính xác 0.17m (p. 47), thể hiện một sự chuyển dịch sang việc phát triển các giải pháp địa không gian tùy chỉnh, chính xác cao cho các vấn đề cụ thể, phức tạp.

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp một cách độc đáo nhiều lý thuyết và phương pháp để giải quyết các vấn đề phức tạp trong tạo DEM và đánh giá nguy cơ lũ lụt.

  • Tích hợp các lý thuyết: Luận án tích hợp các nguyên lý từ Lý thuyết tạo bề mặt địa hình (Surface Generation Theory) thông qua phương pháp BS-Horizon (Nonogaki et al., 2012), các kỹ thuật Viễn thám (Remote Sensing) để thu thập và xử lý dữ liệu (ALOS PALSAR, Landsat TM, RapidEye), và Lý thuyết Phân tích Quyết định (Decision Analysis Theory) thông qua Quy trình Phân cấp Phân tích (AHP) của Saaty. Sự tích hợp này cho phép một cách tiếp cận toàn diện, từ dữ liệu thô đến sản phẩm cuối cùng là bản đồ nguy cơ lũ lụt có giá trị thực tiễn.
  • Cách tiếp cận phân tích mới lạ: Cách tiếp cận phân tích bao gồm một quy trình hai giai đoạn để tạo DEM, phản ánh độ phức tạp của các yêu cầu về dữ liệu:
    1. Hợp nhất DEM đa tiêu chí: Không giống như các phương pháp trước đây (Karkee et al., 2008; Lucca, 2011), luận án này sử dụng trung bình có trọng số động, nơi trọng số được điều chỉnh dựa trên phân loại dạng địa hình (núi, thung lũng, bằng phẳng) được dẫn xuất từ TPI và độ dốc (Dickson và Beier, 2006). Ví dụ, trong vùng bằng phẳng, SRTM được gán trọng số cao hơn GDEM vì "SRTM ... is better than GDEM" (p. 23).
    2. Tối ưu hóa tham số BS-Horizon độc đáo: Luận án thực hiện một phân tích độ nhạy toàn diện và chưa từng có về các tham số của BS-Horizon (Mx, My, α, m1, m2) để tối ưu hóa đồng thời độ mịn bề mặt (J(f)) và độ phù hợp với dữ liệu (R(f)). Việc thiết lập tiêu chí RMSE <= 0.25 để đánh giá bề mặt là một đóng góp cụ thể (p. 41).
  • Đóng góp khái niệm với các định nghĩa:
    • "Geomorphological Approach for DEM Fusion": Định nghĩa một phương pháp hợp nhất DEM không chỉ dựa trên độ chính xác tổng thể mà còn dựa trên hiệu suất thay đổi của các DEM thành phần trong các bối cảnh địa hình cụ thể.
    • "Equality-Inequality Constraints in Surface Generation": Khái niệm này được phát triển sâu sắc hơn thông qua việc chứng minh cách tích hợp dữ liệu "trên" hoặc "dưới" (từ bản đồ địa hình) cùng với dữ liệu điểm thực địa có thể kiểm soát hiệu quả các lỗi nội suy (undershoot/overshoot) ở các vùng địa hình phức tạp, nơi dữ liệu điểm thưa thớt (p. 45).
  • Điều kiện biên (Boundary Conditions) được nêu rõ ràng:
    • Phạm vi địa lý: Các phát hiện liên quan đến tạo DEM độ phân giải cao (sử dụng BS-Horizon) chủ yếu áp dụng cho "small lowland area" (Abstract, p. 9) và "coastal lowland" (Chapter 3, p. 32) với độ cao tương đối thấp (0-10m). Các phương pháp hợp nhất DEM toàn cầu được chứng minh hiệu quả cho "global or regional scale study" (Abstract, p. 9).
    • Chất lượng dữ liệu đầu vào: Độ chính xác của DEM được tạo ra phụ thuộc vào "high density field survey point elevation data" (Chapter 1, p. 12) và chất lượng của dữ liệu viễn thám (ALOS PALSAR, Landsat TM).
    • Điều kiện áp dụng AHP: Phương pháp AHP được đề xuất "demonstrates the effectiveness of applying AHP method in flood hazard assessment for the areas with larger extents" (Abstract, p. 10), nhưng việc xác định trọng số vẫn mang tính chủ quan ở một mức độ nào đó, dựa trên phán đoán chuyên gia.

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

Thiết kế nghiên cứu

Luận án này áp dụng một thiết kế nghiên cứu phức hợp, kết hợp nhiều giai đoạn và phương pháp để đạt được các mục tiêu nghiên cứu.

  • Triết lý nghiên cứu (Research Philosophy): Nghiên cứu chủ yếu theo triết lý thực chứng (positivism). Nó tìm cách thiết lập các mối quan hệ nhân quả và định lượng giữa các yếu tố địa hình và nguy cơ lũ lụt. Sự nhấn mạnh vào độ chính xác có thể đo lường được (RMSE, R²), việc kiểm tra giả thuyết, và việc sử dụng các mô hình toán học (BS-Horizon, AHP) phản ánh một cách tiếp cận khách quan, dựa trên dữ liệu.
  • Phương pháp hỗn hợp (Mixed Methods) với lý do kết hợp cụ thể: Nghiên cứu sử dụng một phương pháp hỗn hợp, tích hợp dữ liệu định lượng (dữ liệu viễn thám vệ tinh như GDEM, SRTM, ALOS PALSAR, Landsat TM; dữ liệu điểm khảo sát thực địa có độ chính xác cao) với các kỹ thuật mô hình hóa (hợp nhất DEM, nội suy BS-Horizon, AHP). Lý do kết hợp là để khắc phục những hạn chế của từng phương pháp đơn lẻ: dữ liệu vệ tinh toàn cầu thiếu độ phân giải cục bộ, trong khi dữ liệu khảo sát thực địa tốn kém và không có sẵn trên diện rộng. Sự kết hợp này cho phép một cách tiếp cận toàn diện từ việc tạo DEM quy mô lớn đến mô hình hóa chi tiết quy mô nhỏ.
  • Thiết kế đa cấp (Multi-level design) với các cấp độ được xác định rõ ràng: Nghiên cứu áp dụng thiết kế đa cấp:
    1. Cấp độ toàn cầu/khu vực: Giai đoạn đầu tập trung vào hợp nhất các DEM toàn cầu (GDEM, SRTM) cho "the whole area of Danang City" (Chapter 1, p. 12) (khoảng 950 km2, Figure 2.1, p. 18) nhằm nâng cao chất lượng của dữ liệu độ cao miễn phí trên phạm vi rộng.
    2. Cấp độ cục bộ/chi tiết: Giai đoạn sau tập trung vào việc tạo DEM độ phân giải cao 5m cho "a small lowland area in the south of Danang City" (Chapter 2, p. 31) (98 km2, Figure 4.1, p. 50), nơi đòi hỏi sự đặc trưng chi tiết của bề mặt địa hình cho mục đích phân vùng nguy cơ lũ lụt.
  • Kích thước mẫu (Sample size) và tiêu chí lựa chọn chính xác:
    • Hợp nhất DEM toàn cầu: Sử dụng toàn bộ GDEM Version 2 và SRTM Version 4.1 phủ khắp Đà Nẵng (khoảng 950 km2), với độ phân giải ban đầu 30m (GDEM) và 90m (SRTM) (Chapter 2, p. 19). DEM tham chiếu được tạo từ bản đồ địa hình 1:10.000 với đường đồng mức 5m và dữ liệu điểm độ cao được khảo sát năm 2009.
    • Tạo DEM BS-Horizon: Sử dụng 9,730 điểm độ cao khảo sát thực địa có mật độ cao (khoảng cách trung bình 40m, độ chính xác ±10cm) làm ràng buộc đẳng thức cho vùng đồng bằng 98 km2 (Chapter 3, p. 35). Ngoài ra, thông tin độ cao liên quan đến đường đồng mức 0-10m từ bản đồ địa hình 1:25.000 năm 2010 được sử dụng làm ràng buộc bất đẳng thức.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu được thiết kế để đảm bảo tính nghiêm ngặt và độ tin cậy của các phát hiện.

  • Chiến lược lấy mẫu (Sampling strategy) với tiêu chí bao gồm/loại trừ:
    • Dữ liệu độ cao: Dữ liệu điểm khảo sát thực địa được thu thập có hệ thống qua lưới mẫu để đảm bảo mật độ cao và phân bố đều trong khu vực nghiên cứu, với tiêu chí bao gồm các điểm độ cao có độ chính xác ±10cm. Không có tiêu chí loại trừ rõ ràng nào được nêu, ngoài việc đảm bảo chất lượng dữ liệu.
    • Dữ liệu vệ tinh: Hình ảnh ALOS PALSAR năm 2007 được chọn vì nó ghi lại một "flood event in 2007" (Chapter 4, p. 49), cung cấp thông tin ngập lụt thực tế. Landsat TM năm 2007 được sử dụng cho phân loại sử dụng đất.
  • Giao thức thu thập dữ liệu (Data collection protocols) với các công cụ được mô tả:
    • Dữ liệu điểm độ cao: Thu thập thông qua khảo sát thực địa (field survey) sử dụng các công cụ đo đạc chính xác (không nêu cụ thể, nhưng độ chính xác ±10cm cho thấy công cụ chuyên nghiệp) (Chapter 3, p. 35).
    • Dữ liệu vệ tinh: ASTER GDEM Version 2 (từ quang học), SRTM Version 4.1 (từ InSAR) (Chapter 2, p. 19). ALOS PALSAR (radar băng L) và Landsat TM (quang học đa phổ).
  • Triangulation (dữ liệu/phương pháp/nhà nghiên cứu/lý thuyết):
    • Triangulation dữ liệu: Sử dụng nhiều nguồn dữ liệu độ cao (GDEM, SRTM, khảo sát thực địa, bản đồ địa hình) và dữ liệu vệ tinh khác nhau (ALOS PALSAR, Landsat TM) để đảm bảo độ tin cậy của DEM và bản đồ ngập lụt.
    • Triangulation phương pháp: Kết hợp các phương pháp khác nhau như hợp nhất DEM dựa trên địa hình, nội suy bi-cubic spline (BS-Horizon) và AHP để giải quyết cùng một vấn đề, tăng cường tính hợp lệ của các phát hiện.
    • Triangulation lý thuyết: Các kết quả được phân tích dưới ánh sáng của lý thuyết Geoinformatics, thủy văn học và khoa học quyết định, đảm bảo sự giải thích toàn diện và sâu sắc.
  • Tính hợp lệ (Validity) và độ tin cậy (Reliability):
    • Độ tin cậy của cấu trúc (Construct validity): Các tham số gây lũ (độ dốc, TWI, EFI, DIST, sử dụng đất) được lựa chọn dựa trên i) vai trò đã được thiết lập của chúng trong các nghiên cứu lũ lụt trước đây (ví dụ: Kenny, 1990; Ballais et al., 2005) và ii) khả năng dẫn xuất chúng từ DEM chất lượng cao.
    • Độ tin cậy bên trong (Internal validity): Quy trình nội suy BS-Horizon được kiểm tra chặt chẽ bằng cách đánh giá độ nhạy của các tham số (Mx, My, α, m1, m2) đến J(f) và R(f) (Chapter 3, p. 37-44), đảm bảo rằng các kết quả tạo bề mặt là ổn định và đáng tin cậy. Các phép kiểm tra độ chính xác (RMSE) và hệ số tương quan (R²) được sử dụng rộng rãi.
    • Độ tin cậy bên ngoài (External validity)/Tính tổng quát hóa: Các phương pháp được phát triển (hợp nhất DEM dựa trên địa hình, tối ưu hóa BS-Horizon, AHP đa tham số) có tiềm năng tổng quát hóa cho các vùng ven biển có đặc điểm địa hình và nguy cơ lũ lụt tương tự. Các điều kiện tổng quát hóa được nêu rõ (p. 61).
    • Độ tin cậy (Reliability): Độ tin cậy được đảm bảo thông qua việc báo cáo các giá trị thống kê định lượng như RMSE (Root Mean Square Error)R² (hệ số tương quan). Ví dụ, DEM hợp nhất có RMSE 11.6m và R² = 0.9986 (Chapter 2, p. 29). DEM BS-Horizon có RMSE 0.17m (Chapter 3, p. 47). Đối với AHP, tính nhất quán được đánh giá bằng Tỷ lệ nhất quán (CR), đảm bảo các trọng số được gán là hợp lý.

Data và phân tích

  • Đặc điểm mẫu (Sample characteristics) với thông tin nhân khẩu học/thống kê:
    • Khu vực nghiên cứu: Vùng đồng bằng ven biển Đà Nẵng, Việt Nam, diện tích 98 km2. Độ cao dao động từ 0m đến 10m (Chapter 4, p. 50).
    • Dữ liệu độ cao: 9,730 điểm độ cao khảo sát thực địa có độ chính xác ±10cm (Chapter 3, p. 35).
    • Dữ liệu vệ tinh: ALOS PALSAR (ảnh radar, 2007), Landsat TM (ảnh quang học, 2007), RapidEye (ảnh quang học, 2014) (Chapter 4, p. 50).
    • Dữ liệu lũ lụt: Dữ liệu trụ lũ khảo sát thực địa (tháng 3 năm 2015) và bản đồ ngập lụt được trích xuất từ ALOS PALSAR 2007.
  • Các kỹ thuật nâng cao (Advanced techniques) với phần mềm:
    • Hợp nhất DEM: Thuật toán trung bình có trọng số (weighted averaging) dựa trên đặc điểm địa hình (phân loại dạng địa hình từ TPI và độ dốc theo Dickson và Beier, 2006). Thuật toán khử nhiễu (denoising algorithm) của Sun et al. (2007) được áp dụng. Phần mềm: GRASS GIS (r.statistics, r.mapcalc, r.denoise, r.surf.rst) (Chapter 2, p. 21, 26, 28).
    • Tạo DEM độ phân giải cao: Phương pháp BS-Horizon (Nonogaki et al., 2012) sử dụng hàm bi-cubic spline và hàm phạt ngoại vi. Chương trình BS-Horizon Fortran và chương trình nội suy APP4OPT FORTRAN (Nonogaki et al., 2009) được sử dụng (Chapter 3, p. 32, 42).
    • Phân vùng nguy cơ lũ lụt: Quy trình Phân cấp Phân tích (AHP) của Saaty, bao gồm so sánh cặp (pair-wise comparison) và tính toán trọng số (Chapter 4, p. 51-52).
    • Xử lý DEM: Nội suy RST (Regularized Spline with Tension) trong GRASS GIS để tạo DEM tham chiếu và nội suy SRTM từ 90m lên 30m (Chapter 2, p. 20-21).
  • Kiểm tra độ vững chắc (Robustness checks) với các thông số kỹ thuật thay thế:
    • Độ vững chắc của phương pháp BS-Horizon được kiểm tra thông qua phân tích độ nhạy tham số rộng rãi (sensitivity analysis) đối với M (số lượng ô), α (tham số phạt), m1 và m2 (tham số làm mịn) (Chapter 3, p. 37-44). Ví dụ, các trường hợp M=50, 100, 200, 400 được so sánh; α thay đổi từ 1 đến 10^12. Các kết quả thống kê (J(f), R(f), Q(f)) được phân tích để tìm ra các thông số tối ưu, đảm bảo rằng các bề mặt được tạo ra ổn định và không nhạy cảm quá mức với các thay đổi nhỏ trong cài đặt.
  • Kích thước hiệu ứng (Effect sizes) và khoảng tin cậy (confidence intervals) được báo cáo:
    • Mặc dù p-values và confidence intervals không được báo cáo trực tiếp, các kích thước hiệu ứng được ngụ ý thông qua các chỉ số định lượng về cải thiện độ chính xác. Ví dụ, việc giảm RMSE từ 14.9m xuống 11.6m cho DEM hợp nhất (Chapter 2, Table 2.5) hoặc đạt được RMSE 0.17m cho DEM BS-Horizon (Chapter 3, p. 47) là những chỉ số mạnh mẽ về kích thước hiệu ứng của các phương pháp được đề xuất. Tỷ lệ gần 84% diện tích ngập lụt thực tế rơi vào vùng nguy cơ cao/rất cao (Abstract, p. 10) cũng là một thước đo kích thước hiệu ứng đáng kể cho mô hình AHP.

Phát hiện đột phá và implications

Những phát hiện then chốt

Luận án này đã đưa ra những phát hiện đột phá với bằng chứng cụ thể từ dữ liệu, vượt qua những hạn chế của các nghiên cứu trước đây.

  1. Hợp nhất DEM dựa trên địa hình cải thiện đáng kể độ chính xác DEM toàn cầu: Phương pháp trung bình có trọng số dựa trên bối cảnh địa hình (núi, thung lũng, bằng phẳng) đã cải thiện hiệu quả chất lượng của các DEM toàn cầu. RMSE của DEM hợp nhất giảm xuống còn 11.6m, một sự cải thiện đáng kể so với GDEM ban đầu (RMSE 14.9m) và SRTM (RMSE 14.8m) (Chapter 2, Table 2.5, p. 29). Hơn nữa, "most artifacts were successfully eliminated" (Abstract, p. 9) và hệ số tương quan R² với DEM tham chiếu đạt 0.9986 (Chapter 2, p. 29).
  2. Tối ưu hóa các tham số BS-Horizon cho tạo DEM độ phân giải cao: Nghiên cứu đã xác định các thiết lập tham số tối ưu cho phương pháp BS-Horizon lần đầu tiên. Với dữ liệu đẳng thức-bất đẳng thức, các giá trị Mx = My = 200, α = 1.58×10^4 và m1 = 0.5 (Chapter 3, p. 46) được chứng minh là tạo ra DEM 5m độ phân giải cao với RMSE chỉ 0.17m (Chapter 3, p. 47). Phát hiện này là rất quan trọng vì nó cung cấp hướng dẫn thực tế để đạt được độ chính xác cao trong các vùng đồng bằng ven biển phức tạp, nơi các DEM toàn cầu thất bại.
  3. DEM độ phân giải cao vượt trội cho vùng đất thấp: DEM 5m được tạo ra bằng phương pháp BS-Horizon cho vùng đồng bằng ven biển Đà Nẵng đã thể hiện độ chính xác vượt trội (RMSE 0.17m) so với DEM hợp nhất (RMSE 1.98m) (Chapter 3, p. 47). Điều này xác nhận rằng các phương pháp chuyên biệt là cần thiết cho các vùng đất thấp nơi địa hình rất tinh tế và liên quan trực tiếp đến nguy cơ lũ lụt.
  4. Mô hình AHP đa tham số hiệu quả và được xác thực cho phân vùng nguy cơ lũ lụt: Mô hình AHP, tích hợp EFI, DIST, sử dụng đất, độ dốc và TWI, đã tạo ra bản đồ phân vùng nguy cơ lũ lụt có độ chính xác cao. "Nearly 84 percent of inundated areas observed in ALOS PALSAR data fall in the high and very high flood hazard zones" (Abstract, p. 10). Sự so sánh với dữ liệu trụ lũ thực địa cũng cho thấy "good correspondence with demarcated zones and estimated depth" (Abstract, p. 10).
  5. Hiệu quả của ràng buộc bất đẳng thức trong nội suy spline: Phát hiện cho thấy việc tích hợp ràng buộc bất đẳng thức ("trên" hoặc "dưới" một độ cao nhất định) có thể kiểm soát hiệu quả các "strong oscillations" và "significant undershoot elevation" (p. 45) thường xảy ra khi chỉ sử dụng ràng buộc đẳng thức, đặc biệt ở các khu vực không có dữ liệu điểm hoặc có sự chuyển đổi địa hình đột ngột (ví dụ: sông, đường). Điều này cải thiện độ tin cậy của bề mặt nội suy trong các khu vực không có dữ liệu.

Kết quả trái ngược/ngoại lệ (Counter-intuitive results) với giải thích lý thuyết: Mặc dù không có "kết quả trái ngược" lớn, nhưng có một sự sắc thái đáng chú ý: trong khi tăng tham số phạt α trong BS-Horizon làm giảm R(f) (độ phù hợp với ràng buộc), nó cũng làm tăng J(f) (mất độ mịn bề mặt) một cách hàm mũ, và do đó Q(f) cũng tăng theo cấp số nhân (Chapter 3, p. 38-39). Điều này nhấn mạnh tầm quan trọng của việc tìm kiếm sự cân bằng, thay vì chỉ đơn thuần tìm kiếm giá trị R(f) thấp nhất. Điều này giải thích tại sao không phải lúc nào mức độ phù hợp cao nhất cũng tạo ra bề mặt địa hình tốt nhất, vì nó có thể dẫn đến bề mặt quá lồi lõm với "anomalies with number of undershoot and overshoot elevation" (Chapter 3, p. 40).

Hiện tượng mới (New phenomena) với ví dụ cụ thể từ dữ liệu: Nghiên cứu làm sáng tỏ một hiện tượng đã được biết nhưng được mô hình hóa tốt hơn: "undershoot elevation" (độ cao thấp hơn thực tế) ở các vùng sông hoặc ranh giới giữa các loại sử dụng đất khác nhau khi chỉ sử dụng ràng buộc đẳng thức. Luận án chỉ ra rằng "These undershoots mainly locate on rivers where there are no field data or on the transition between two different land-use, mainly at the boundaries between paddy field and rivers, roads or other man-made structures" (p. 45). Việc giải quyết hiện tượng này bằng ràng buộc bất đẳng thức cung cấp một sự hiểu biết sâu sắc hơn về các hạn chế của nội suy spline trong các bối cảnh địa hình phức tạp, nơi có sự gián đoạn rõ rệt.

So sánh với các phát hiện nghiên cứu trước:

  • Các phát hiện về độ chính xác của GDEM và SRTM (ví dụ: lỗi cao ở vùng núi, thấp ở vùng bằng phẳng, Figure 2.8, p. 24) phù hợp với các đánh giá trước đây của Reuter et al. (2007) và Mukherjee et al. (2013a). Tuy nhiên, luận án này tiến xa hơn bằng cách sử dụng những hiểu biết này để xây dựng một phương pháp hợp nhất được tối ưu hóa.
  • Các nghiên cứu trước đây về phân vùng nguy cơ lũ lụt ở Đà Nẵng (Ho et al., 2012; Do et al., 2014) bị hạn chế bởi thiếu DEM độ phân giải cao và dữ liệu vệ tinh về các sự kiện gần lũ. Phát hiện của luận án này về độ chính xác 84% trong phân vùng nguy cơ lũ lụt (Abstract, p. 10) cho thấy một bước nhảy vọt đáng kể về độ tin cậy, trực tiếp khắc phục những hạn chế đó.

Implications đa chiều

Các phát hiện của luận án có những hàm ý sâu rộng trên nhiều khía cạnh.

  • Tiến bộ lý thuyết (Theoretical advances) với đóng góp cho 2+ lý thuyết:
    • Lý thuyết tạo bề mặt địa hình: Nghiên cứu mở rộng lý thuyết BS-Horizon (Nonogaki et al., 2012) bằng cách cung cấp một khuôn khổ toàn diện để tối ưu hóa các tham số của nó. Điều này làm sâu sắc thêm sự hiểu biết về cách các hàm bi-cubic spline và hàm phạt ngoại vi tương tác trong các tình huống thực tế, đặc biệt với dữ liệu điểm thưa thớt và các ràng buộc về độ cao. Nó cung cấp một lý thuyết được tinh chỉnh để tạo ra các mô hình bề mặt địa hình chính xác cho các vùng đồng bằng.
    • Lý thuyết hợp nhất dữ liệu địa không gian: Bằng cách chứng minh tính ưu việt của phương pháp hợp nhất dựa trên địa hình so với các phương pháp tổng quát hơn, luận án đóng góp vào lý thuyết về hợp nhất dữ liệu, cho thấy rằng độ tin cậy của dữ liệu đầu vào không phải là cố định mà thay đổi theo bối cảnh địa lý và đặc tính của cảm biến.
  • Đổi mới phương pháp luận (Methodological innovations) áp dụng cho các bối cảnh khác:
    • Giao thức tối ưu hóa BS-Horizon: Quy trình có hệ thống để đánh giá và lựa chọn các tham số BS-Horizon (M, α, m1, m2) có thể được áp dụng như một phương pháp luận tiêu chuẩn cho các nghiên cứu tạo bề mặt địa hình khác sử dụng phương pháp spline, đặc biệt ở những khu vực có dữ liệu điểm mật độ cao và nhu cầu về DEM độ phân giải cực cao.
    • Khung hợp nhất DEM thích ứng: Phương pháp hợp nhất DEM dựa trên địa hình có thể được áp dụng cho các khu vực khác trên thế giới để nâng cao chất lượng của các DEM toàn cầu miễn phí, đặc biệt là ở các quốc gia đang phát triển nơi dữ liệu độ cao chất lượng cao còn hạn chế.
    • Mô hình AHP đa tham số: Cách tiếp cận AHP tích hợp nhiều yếu tố gây lũ có thể được điều chỉnh cho các loại hình thiên tai khác (ví dụ: sạt lở đất, cháy rừng) hoặc các vấn đề môi trường khác yêu cầu tích hợp các yếu tố địa lý đa dạng.
  • Ứng dụng thực tiễn (Practical applications) với các khuyến nghị cụ thể:
    • Kế hoạch phòng chống lũ lụt: DEM độ phân giải cao và bản đồ phân vùng nguy cơ lũ lụt chính xác là những công cụ quan trọng cho các nhà quy hoạch đô thị và quản lý thiên tai. Các bản đồ này có thể được sử dụng để xác định các khu vực có nguy cơ cao để ưu tiên các biện pháp giảm thiểu (ví dụ: xây dựng đê điều, sơ tán, quy hoạch sử dụng đất).
    • Quy hoạch sử dụng đất: Các nhà quy hoạch có thể sử dụng dữ liệu này để đưa ra các quyết định sáng suốt về quy hoạch phát triển đô thị, cơ sở hạ tầng và nông nghiệp tại các vùng ven biển, tránh xây dựng ở các khu vực có nguy cơ cao.
    • Hệ thống cảnh báo sớm: DEM chính xác là nền tảng cho việc phát triển các mô hình thủy văn thời gian thực, hỗ trợ các hệ thống cảnh báo sớm lũ lụt.
  • Khuyến nghị chính sách (Policy recommendations) với lộ trình thực hiện:
    • Đầu tư vào dữ liệu độ cao chất lượng cao: Các chính phủ và cơ quan địa phương nên ưu tiên đầu tư vào việc thu thập dữ liệu điểm độ cao mật độ cao cho các vùng ven biển dễ bị tổn thương, vì đây là nền tảng cho việc tạo DEM chính xác.
    • Áp dụng phương pháp luận tiên tiến: Các cơ quan quản lý thiên tai và quy hoạch nên áp dụng các phương pháp luận tiên tiến được phát triển trong luận án này (ví dụ: tối ưu hóa BS-Horizon, AHP đa tham số) để đánh giá và phân vùng nguy cơ lũ lụt, thay vì dựa vào các phương pháp cũ, kém chính xác hơn.
    • Lồng ghép dữ liệu vệ tinh sự kiện gần lũ: Khuyến khích sử dụng thường xuyên dữ liệu vệ tinh về các sự kiện lũ lụt trong quá khứ để xác thực và hiệu chỉnh các mô hình nguy cơ lũ lụt.
  • Điều kiện tổng quát hóa (Generalizability conditions) được nêu rõ: Các phương pháp tạo DEM và phân vùng nguy cơ lũ lụt được đề xuất có thể tổng quát hóa cho các vùng đồng bằng ven biển khác có đặc điểm địa hình, thủy văn và sẵn có dữ liệu tương tự. Đặc biệt, các khu vực thường xuyên bị lũ lụt, có dữ liệu điểm khảo sát chất lượng cao và tiềm năng sử dụng dữ liệu vệ tinh về các sự kiện lũ lụt sẽ hưởng lợi nhiều nhất. Tuy nhiên, việc điều chỉnh các tham số (ví dụ: cho BS-Horizon) có thể cần thiết tùy thuộc vào đặc điểm riêng của từng khu vực nghiên cứu.

Limitations và Future Research

Mặc dù luận án đã đạt được những tiến bộ đáng kể, nhưng nó cũng nhận thức rõ những hạn chế và điều kiện biên của mình.

  • 3-4 hạn chế cụ thể được thừa nhận:
    1. Tính sẵn có của dữ liệu điểm khảo sát: Phương pháp BS-Horizon độ phân giải cao phụ thuộc vào dữ liệu điểm khảo sát thực địa mật độ cao (9.730 điểm, p. 35). Sự sẵn có của dữ liệu như vậy có thể hạn chế khả năng áp dụng phương pháp này cho các khu vực rộng lớn hơn hoặc các vùng mà việc thu thập dữ liệu tốn kém.
    2. Tính chủ quan trong AHP: Mặc dù AHP là một phương pháp mạnh mẽ, việc gán trọng số thông qua so sánh cặp vẫn có thể bị ảnh hưởng bởi phán đoán chủ quan của các chuyên gia (Chapter 4, p. 51). Điều này có thể ảnh hưởng đến kết quả phân vùng nguy cơ lũ lụt.
    3. Dữ liệu sử dụng đất tĩnh: Dữ liệu sử dụng đất được lấy từ Landsat TM năm 2007 (Chapter 4, p. 50). Tuy nhiên, sử dụng đất có thể thay đổi nhanh chóng, đặc biệt ở các vùng ven biển đang phát triển, điều này có thể ảnh hưởng đến độ chính xác của các phân tích theo thời gian nếu không được cập nhật.
    4. Hạn chế của DEM hợp nhất ở vùng đất thấp: Mặc dù DEM hợp nhất cải thiện độ chính xác tổng thể, nó vẫn có những hạn chế về độ phân giải và độ chính xác ở các vùng đồng bằng thấp nhỏ ("fused DEM is insufficient... fused DEM presents limitations," Chapter 2, p. 31). RMSE 1.98m của nó ở vùng đất thấp vẫn cao hơn đáng kể so với DEM BS-Horizon (0.17m), cho thấy một phương pháp đơn lẻ không đủ cho mọi bối cảnh.
  • Điều kiện biên về bối cảnh/mẫu/thời gian:
    • Bối cảnh: Các kết quả tối ưu hóa tham số BS-Horizon và tính hiệu quả của mô hình AHP được thiết lập trong bối cảnh cụ thể của vùng đồng bằng ven biển Đà Nẵng, Việt Nam. Các đặc điểm địa hình độc đáo (độ cao thấp 0-10m, mạng lưới sông phức tạp, mùa mưa rõ rệt) có thể yêu cầu điều chỉnh khi áp dụng cho các bối cảnh khác.
    • Mẫu: Độ chính xác cao của DEM BS-Horizon phụ thuộc vào chất lượng và mật độ của mẫu điểm độ cao được khảo sát.
    • Thời gian: Dữ liệu vệ tinh và khảo sát thực địa được sử dụng là tại các thời điểm cụ thể (ví dụ: ALOS PALSAR 2007, khảo sát trụ lũ 2015). Sự thay đổi khí hậu và phát triển đô thị liên tục có thể làm thay đổi các yếu tố gây lũ, đòi hỏi cập nhật định kỳ.
  • Chương trình nghiên cứu trong tương lai (Future research agenda) với 4-5 hướng cụ thể:
    1. Tích hợp dữ liệu LiDAR: Điều tra tiềm năng tích hợp dữ liệu LiDAR (Light Detection And Ranging) để tạo DEM độ phân giải cực cao, đặc biệt cho các khu vực đô thị và cơ sở hạ tầng quan trọng, vượt qua giới hạn của dữ liệu khảo sát điểm (vốn tốn kém).
    2. Mô hình hóa lũ lụt động lực học: Kết hợp DEM độ phân giải cao được tạo ra với các mô hình thủy văn động lực học 2D hoặc 3D (ví dụ: HEC-RAS 2D, MIKE 21C) để mô phỏng hành vi lũ lụt theo thời gian, cung cấp thông tin chi tiết hơn về tốc độ dòng chảy và thời gian ngập lụt.
    3. Phân tích sự thay đổi sử dụng đất/độ cao: Nghiên cứu tác động của sự thay đổi sử dụng đất và thay đổi độ cao do con người (ví dụ: san lấp mặt bằng) đối với nguy cơ lũ lụt theo thời gian, sử dụng chuỗi thời gian dữ liệu vệ tinh và DEM.
    4. Kết hợp AHP với các phương pháp học máy: Khám phá sự kết hợp của AHP với các thuật toán học máy (ví dụ: Support Vector Machine, Random Forest) để cải thiện quá trình gán trọng số và tự động hóa việc phân vùng nguy cơ lũ lụt.
    5. Mô hình hóa nguy cơ tổng hợp (Compound Flood Risk): Mở rộng mô hình để bao gồm các yếu tố nguy cơ lũ lụt phức hợp như lũ do triều cường, nước dâng do bão và mưa lớn đồng thời, đặc biệt quan trọng đối với các vùng ven biển.
  • Những cải tiến về phương pháp luận được đề xuất:
    • Phát triển các thuật toán hiệu quả hơn để tự động hóa quá trình xác định các tham số tối ưu cho BS-Horizon, giảm thiểu sự cần thiết của phân tích độ nhạy thủ công.
    • Nghiên cứu các phương pháp định lượng hóa tính chủ quan trong AHP, chẳng hạn như sử dụng AHP mờ (Fuzzy AHP) hoặc tích hợp ý kiến từ nhiều chuyên gia hơn để tạo ra trọng số mạnh mẽ hơn.
  • Đề xuất mở rộng lý thuyết: Đề xuất một lý thuyết toàn diện hơn về việc tích hợp các ràng buộc bất đẳng thức trong nội suy bề mặt cho các bề mặt có sự gián đoạn rõ ràng (ví dụ: đê, kè, ranh giới thảm thực vật), không chỉ giới hạn ở các ràng buộc độ cao mà còn mở rộng sang các đặc tính địa hình khác.

Tác động và ảnh hưởng

Luận án này có tiềm năng tạo ra tác động và ảnh hưởng đáng kể trên nhiều lĩnh vực, từ học thuật đến xã hội.

  • Tác động học thuật (Academic impact) với ước tính số trích dẫn tiềm năng:
    • Nâng cao hiểu biết về tạo DEM: Các đóng góp về tối ưu hóa tham số BS-Horizon và phương pháp hợp nhất DEM dựa trên địa hình sẽ là tài liệu tham khảo quan trọng cho các nhà nghiên cứu trong lĩnh vực Geoinformatics và viễn thám. Ước tính có thể đạt 50-100 trích dẫn trong 5-10 năm tới, đặc biệt từ các nghiên cứu tập trung vào tạo DEM độ phân giải cao và mô hình hóa bề mặt địa hình.
    • Cải thiện mô hình hóa nguy cơ thiên tai: Cách tiếp cận AHP đa tham số được xác thực nghiêm ngặt sẽ cung cấp một phương pháp luận mạnh mẽ cho các nghiên cứu về nguy cơ lũ lụt và các thiên tai địa không gian khác, thu hút các trích dẫn từ các nhà khoa học môi trường, thủy văn và quản lý thiên tai.
    • Mở ra các hướng nghiên cứu mới: Việc chỉ ra hiệu quả của ràng buộc bất đẳng thức trong nội suy spline mở ra các cơ hội nghiên cứu về cách ứng dụng các loại ràng buộc khác để cải thiện mô hình bề mặt địa hình.
  • Chuyển đổi ngành công nghiệp (Industry transformation) với các lĩnh vực cụ thể:
    • Quy hoạch đô thị và phát triển cơ sở hạ tầng: Các công ty quy hoạch và xây dựng có thể sử dụng các DEM độ phân giải cao và bản đồ nguy cơ lũ lụt chính xác để đưa ra các quyết định sáng suốt hơn về vị trí và thiết kế các dự án phát triển, giảm thiểu rủi ro và tổn thất liên quan đến lũ lụt. Điều này có thể áp dụng cho các ngành bất động sản, xây dựng và bảo hiểm.
    • Phát triển phần mềm địa không gian: Các nhà phát triển phần mềm có thể tích hợp các thuật toán tối ưu hóa BS-Horizon hoặc các nguyên tắc hợp nhất DEM dựa trên địa hình vào các công cụ GIS thương mại, giúp mở rộng khả năng của chúng trong việc tạo ra các mô hình bề mặt địa hình chất lượng cao.
    • Ngành bảo hiểm: Các công ty bảo hiểm có thể sử dụng các bản đồ nguy cơ lũ lụt chi tiết để đánh giá rủi ro chính xác hơn, định giá phí bảo hiểm công bằng hơn và phát triển các sản phẩm bảo hiểm mới cho các khu vực có nguy cơ cao.
  • Ảnh hưởng chính sách (Policy influence) với các cấp chính phủ:
    • Cấp chính quyền địa phương (ví dụ: Sở Tài nguyên và Môi trường Đà Nẵng): Luận án cung cấp "specific recommendations" (p. 60) và dữ liệu đáng tin cậy để cập nhật các kế hoạch phòng chống lũ lụt của thành phố, quy hoạch sử dụng đất và chính sách giảm thiểu rủi ro. Việc tích hợp DEM 5m và bản đồ nguy cơ lũ lụt vào hệ thống thông tin địa lý của thành phố sẽ hỗ trợ ra quyết định dựa trên bằng chứng.
    • Cấp chính quyền quốc gia (ví dụ: Bộ Tài nguyên và Môi trường Việt Nam): Các phương pháp luận được phát triển có thể được nhân rộng và áp dụng cho các vùng ven biển dễ bị lũ lụt khác trên toàn quốc, hỗ trợ việc xây dựng các chiến lược quản lý thiên tai cấp quốc gia.
  • Lợi ích xã hội (Societal benefits) được định lượng hóa nếu có thể:
    • Giảm thiểu thiệt hại do lũ lụt: Bằng cách cung cấp thông tin chính xác hơn về các khu vực có nguy cơ cao, luận án góp phần vào việc giảm thiểu thiệt hại về người và của. Ví dụ, trong giai đoạn 1997-2009, Đà Nẵng đã mất 6803 tỷ VND do lũ lụt (Chapter 1, p. 13). Việc giảm một phần nhỏ của con số này (ví dụ: 5-10%) thông qua quy hoạch tốt hơn sẽ mang lại lợi ích kinh tế đáng kể.
    • Nâng cao an toàn cộng đồng: Các bản đồ nguy cơ lũ lụt chính xác sẽ hỗ trợ các hệ thống cảnh báo sớm và kế hoạch sơ tán, giúp cứu sống và bảo vệ tài sản của người dân.
    • Thúc đẩy phát triển bền vững: Bằng cách lồng ghép rủi ro lũ lụt vào quy hoạch phát triển, nghiên cứu này góp phần vào việc xây dựng các cộng đồng ven biển có khả năng chống chịu tốt hơn và bền vững hơn.
  • Sự liên quan quốc tế (International relevance) với hàm ý toàn cầu:
    • Các vùng đồng bằng ven biển trên toàn thế giới đang đối mặt với nguy cơ lũ lụt gia tăng do biến đổi khí hậu và nước biển dâng. Các phương pháp tạo DEM tiên tiến và mô hình hóa nguy cơ lũ lụt được phát triển trong luận án có thể được áp dụng rộng rãi cho các thành phố ven biển và vùng đất thấp khác ở các quốc gia đang phát triển (ví dụ: ở Đông Nam Á, châu Phi, Mỹ Latinh) nơi có những hạn chế tương tự về dữ liệu độ cao chất lượng cao. Luận án cung cấp một mô hình cho việc quản lý rủi ro lũ lụt ở những bối cảnh như vậy.

Đối tượng hưởng lợi

Luận án này mang lại lợi ích cụ thể cho nhiều đối tượng khác nhau trong giới học thuật, ngành công nghiệp và hoạch định chính sách.

  • Các nhà nghiên cứu tiến sĩ (Doctoral researchers):
    • Khoảng trống nghiên cứu cụ thể: Cung cấp một ví dụ xuất sắc về cách xác định và lấp đầy các khoảng trống nghiên cứu cụ thể trong Geoinformatics, đặc biệt là trong các lĩnh vực hợp nhất dữ liệu và tạo mô hình bề mặt. Việc làm rõ các hạn chế của các nghiên cứu trước (ví dụ: Ho et al., 2012; Do et al., 2014) và cách luận án này giải quyết chúng là một hướng dẫn có giá trị.
    • Phương pháp luận nghiêm ngặt: Cung cấp một khuôn mẫu chi tiết về việc tiến hành phân tích độ nhạy tham số (cho BS-Horizon), xác thực mô hình đa chiều (AHP) và tích hợp các loại dữ liệu đa dạng.
    • Sử dụng hiệu quả DEM miễn phí: Cho thấy cách các DEM toàn cầu miễn phí có thể được xử lý và nâng cao chất lượng để phù hợp với các mục đích nghiên cứu phức tạp hơn, có ý nghĩa đối với các nhà nghiên cứu có nguồn lực hạn chế.
  • Các học giả cao cấp (Senior academics):
    • Tiến bộ lý thuyết: Các đóng góp vào lý thuyết BS-Horizon (Nonogaki et al., 2012) thông qua việc tối ưu hóa tham số và sự mở rộng của lý thuyết hợp nhất dữ liệu địa không gian sẽ được các học giả trong lĩnh vực Geoinformatics đánh giá cao.
    • Khung phân tích mới: Khung phân tích độc đáo kết hợp các phương pháp khác nhau (DEM fusion, BS-Horizon, AHP) cung cấp một mô hình cho các nghiên cứu liên ngành phức tạp.
    • Dữ liệu thực nghiệm mới: Các phát hiện định lượng (ví dụ: RMSE 0.17m cho DEM 5m, 84% độ chính xác AHP) cung cấp bằng chứng thực nghiệm mạnh mẽ để thúc đẩy các cuộc thảo luận khoa học và lý thuyết trong các lĩnh vực tương ứng.
  • R&D ngành công nghiệp (Industry R&D):
    • Ứng dụng thực tiễn: Các phương pháp tạo DEM độ phân giải cao có thể được các công ty tư vấn về môi trường, kỹ thuật và xây dựng sử dụng để thực hiện các đánh giá rủi ro địa chất và thủy văn chính xác hơn.
    • Công cụ và quy trình cải tiến: Các giao thức tối ưu hóa BS-Horizon cung cấp một quy trình rõ ràng để các nhóm R&D có thể phát triển các công cụ phần mềm hoặc dịch vụ tốt hơn để tạo mô hình bề mặt địa hình chính xác cho các ứng dụng như quy hoạch hạ tầng, phát triển đô thị.
    • Phân tích rủi ro chính xác hơn: Các công ty bảo hiểm và quản lý rủi ro có thể sử dụng các bản đồ nguy cơ lũ lụt được tạo ra để định lượng rủi ro tốt hơn, dẫn đến các quyết định kinh doanh hiệu quả hơn.
  • Các nhà hoạch định chính sách (Policy makers):
    • Khuyến nghị dựa trên bằng chứng: Luận án cung cấp các bản đồ phân vùng nguy cơ lũ lụt chi tiết và các khuyến nghị chính sách cụ thể (ví dụ: về đầu tư dữ liệu, áp dụng phương pháp luận) (p. 60) mà các cơ quan chính phủ có thể sử dụng để phát triển các chính sách quản lý thiên tai và quy hoạch sử dụng đất hiệu quả hơn.
    • Hiểu biết rõ ràng về các khu vực dễ bị tổn thương: Các bản đồ được tạo ra xác định rõ ràng các khu vực có nguy cơ lũ lụt cao và rất cao (ví dụ: 84% diện tích ngập lụt thực tế rơi vào các khu vực này, Abstract, p. 10), cho phép các nhà hoạch định chính sách ưu tiên các nguồn lực và can thiệp vào những nơi cần thiết nhất.
    • Công cụ hỗ trợ ra quyết định: Luận án cung cấp một mô hình tổng hợp thông tin địa lý phức tạp thành một định dạng dễ hiểu, hữu ích cho việc ra quyết định ở các cấp chính quyền khác nhau. Định lượng lợi ích nơi có thể:
  • Giảm thiểu tổn thất kinh tế do lũ lụt, ví dụ, nếu bản đồ nguy cơ lũ lụt chính xác giúp giảm 10% thiệt hại ở Đà Nẵng, điều này có thể tiết kiệm hàng trăm tỷ VND mỗi năm dựa trên dữ liệu thiệt hại 6803 tỷ VND từ 1997-2009 (Chapter 1, p. 13).
  • Cải thiện độ chính xác của DEM từ ~15m xuống 0.17m (p. 29, 47) là một lợi ích định lượng rõ ràng cho bất kỳ ứng dụng nào phụ thuộc vào dữ liệu độ cao.

Câu hỏi chuyên sâu

Trả lời với CÁC CHI TIẾT CỤ THỂ:

  1. Đóng góp lý thuyết độc đáo nhất (tên lý thuyết được mở rộng): Đóng góp lý thuyết độc đáo nhất là việc mở rộng và tinh chỉnh lý thuyết về tạo bề mặt địa hình BS-Horizon (Nonogaki et al., 2012). Luận án là nghiên cứu đầu tiên thực hiện "fully evaluated" (Abstract, p. 10) ảnh hưởng của tất cả các tham số cài đặt (Mx, My, α, m1, m2) đến kết quả tạo bề mặt. Điều này cung cấp một sự hiểu biết sâu sắc chưa từng có về cách các hàm bi-cubic spline và hàm phạt ngoại vi tương tác để tạo ra bề mặt địa hình, đặc biệt là khi tích hợp ràng buộc đẳng thức và bất đẳng thức. Cụ thể, nó đã xác định được các vùng tham số tối ưu (ví dụ: Mx=My=200, α=1.58x10^4, m1=0.5) (Chapter 3, p. 46) để đạt được sự cân bằng giữa độ mịn bề mặt (J(f)) và độ phù hợp với dữ liệu (R(f)), điều mà Nonogaki et al. (2012) không đi sâu vào. Điều này nâng cao khả năng áp dụng thực tế và tính vững chắc của phương pháp BS-Horizon.

  2. Đổi mới phương pháp luận (so sánh với 2+ nghiên cứu trước): Đổi mới phương pháp luận nổi bật là phương pháp hợp nhất DEM dựa trên đặc điểm địa hình (geomorphological approach for DEM fusion). Không giống như các phương pháp hợp nhất DEM trước đây của Karkee et al. (2008) sử dụng FFT hoặc Lucca (2011) so sánh trung bình có trọng số tổng quát, luận án này sử dụng một thuật toán trung bình có trọng số linh hoạt hơn. Các trọng số cho GDEM và SRTM được xác định động dựa trên sự thay đổi về độ chính xác của chúng trong các loại hình địa hình cụ thể (núi, thung lũng, bằng phẳng) được phân loại bằng TPI và độ dốc (Dickson và Beier, 2006) (Chapter 2, p. 26). Ví dụ, phương trình Fused DEM = (GDEM + SRTM*2)/3 được sử dụng ở vùng bằng phẳng, trong khi Fused DEM = (GDEM*2 + SRTM)/3 được sử dụng ở thung lũng (Chapter 2, p. 27-28), phản ánh hiệu suất khác nhau của GDEM và SRTM (Figure 2.5, p. 23). Sự đổi mới này đã dẫn đến việc giảm RMSE của DEM hợp nhất xuống còn 11.6m, cải thiện đáng kể so với GDEM (14.9m) và SRTM (14.8m) riêng lẻ (Chapter 2, p. 29), một cấp độ chi tiết và độ chính xác mà các phương pháp hợp nhất tổng quát hơn không đạt được.

  3. Phát hiện đáng ngạc nhiên nhất (với dữ liệu hỗ trợ): Phát hiện đáng ngạc nhiên nhất là việc hàm phạt ngoại vi (exterior penalty function) làm tăng J(f) (độ mịn bề mặt) một cách hàm mũ khi α tăng, dẫn đến việc mất độ mịn và tạo ra các bề mặt lồi lõm (anomalies with undershoot/overshoot elevation), ngay cả khi R(f) (độ phù hợp) đang giảm (Chapter 3, p. 38-40). Điều này trái ngược với kỳ vọng trực giác rằng việc cải thiện độ phù hợp với ràng buộc sẽ luôn dẫn đến một mô hình tốt hơn. Dữ liệu hỗ trợ là các biểu đồ của J(f), R(f) và Q(f) trong Figure 3.4 và Figure 3.6, cho thấy Q(f) tăng theo cấp số nhân với α, và các ví dụ về DEM được tạo ra ở các giá trị α cao có "elevation ranges... from 773.0×108" hoặc "from 564.0×108" (Chapter 3, p. 39). Điều này nhấn mạnh rằng việc tối ưu hóa phải tìm kiếm sự cân bằng giữa độ phù hợp và độ mịn, chứ không chỉ đơn thuần giảm thiểu lỗi.

  4. Giao thức tái lập được cung cấp? Có, giao thức tái lập cho việc tạo DEM độ phân giải cao bằng phương pháp BS-Horizon đã được cung cấp một cách chi tiết. Luận án đã mô tả cặn kẽ các bước cần thiết để tái lập các phát hiện, bao gồm:

    • Dữ liệu đầu vào: Yêu cầu dữ liệu điểm độ cao khảo sát thực địa mật độ cao (9,730 điểm với độ chính xác ±10cm) làm ràng buộc đẳng thức, và thông tin độ cao từ bản đồ địa hình (0-10m) làm ràng buộc bất đẳng thức (Chapter 3, p. 35-36).
    • Phần mềm: Cần chương trình BS-Horizon Fortran và chương trình nội suy APP4OPT FORTRAN (Nonogaki et al., 2012; Nonogaki et al., 2009) (Chapter 3, p. 32, 42).
    • Thiết lập tham số tối ưu: Các giá trị tham số tối ưu đã được xác định rõ ràng: Mx = My = 200, α = 1.58×10^4 (khi sử dụng dữ liệu đẳng thức-bất đẳng thức), và m1 = 0.5 (Chapter 3, p. 46).
    • Tiêu chí đánh giá: Tiêu chí RMSE <= 0.25m (hoặc R(f) <= 0.0625m) được đề xuất để xác định bề mặt tối ưu (Chapter 3, p. 41). Những chi tiết này cho phép các nhà nghiên cứu khác tái tạo quy trình và kết quả tạo DEM 5m độ phân giải cao.
  5. Chương trình nghiên cứu 10 năm được phác thảo? Có, một chương trình nghiên cứu tương lai (future research agenda) kéo dài hơn 10 năm đã được phác thảo, bao gồm 4-5 hướng nghiên cứu cụ thể với các cải tiến về phương pháp luận và mở rộng lý thuyết (Chapter 6, p. 61). Các hướng này bao gồm:

    • Tích hợp dữ liệu LiDAR: Điều tra tiềm năng của dữ liệu LiDAR để tạo DEM độ phân giải cực cao, vượt qua giới hạn của dữ liệu điểm khảo sát hiện tại.
    • Mô hình hóa lũ lụt động lực học: Kết hợp DEM độ phân giải cao với các mô hình thủy văn động lực học 2D/3D (ví dụ: HEC-RAS 2D, MIKE 21C) để mô phỏng lũ lụt theo thời gian thực và dự báo.
    • Phân tích sự thay đổi sử dụng đất/độ cao: Nghiên cứu tác động lâu dài của sự thay đổi môi trường do con người đối với nguy cơ lũ lụt bằng cách sử dụng dữ liệu chuỗi thời gian.
    • Kết hợp AHP với học máy: Khám phá tích hợp AHP với các thuật toán học máy (ví dụ: SVM, Random Forest) để cải thiện và tự động hóa quá trình phân vùng nguy cơ lũ lụt.
    • Mô hình hóa nguy cơ tổng hợp: Mở rộng mô hình để bao gồm các yếu tố lũ lụt phức tạp hơn như nước dâng do bão và mưa lớn đồng thời. Các đề xuất này không chỉ giải quyết các hạn chế hiện tại mà còn mở ra những con đường mới cho nghiên cứu trong Geoinformatics, thủy văn học và quản lý thiên tai trong thập kỷ tới.

Kết luận

Luận án này đã tạo ra những đóng góp cụ thể và đáng kể trong lĩnh vực Geoinformatics và quản lý thiên tai, đặc biệt là trong bối cảnh các vùng đồng bằng ven biển dễ bị lũ lụt.

  1. Phát triển một phương pháp hợp nhất DEM toàn cầu dựa trên địa hình tiên tiến, đạt được RMSE 11.6m, cải thiện đáng kể so với GDEM (14.9m) và SRTM (14.8m) riêng lẻ (Chapter 2, p. 29). Điều này cung cấp một giải pháp hiệu quả để nâng cao chất lượng của dữ liệu độ cao miễn phí ở các khu vực thiếu dữ liệu chất lượng cao.
  2. Thực hiện đánh giá tham số toàn diện lần đầu tiên cho phương pháp tạo bề mặt BS-Horizon, xác định các thiết lập tối ưu (Mx=My=200, α=1.58x10^4, m1=0.5) (Chapter 3, p. 46) để tạo DEM 5m độ phân giải cao với RMSE chỉ 0.17m (Chapter 3, p. 47). Điều này thiết lập một giao thức rõ ràng cho việc sử dụng BS-Horizon trong các ứng dụng thực tế.
  3. Tạo ra một DEM độ phân giải cao (5m) và chính xác cho vùng đồng bằng ven biển, trực tiếp khắc phục hạn chế về dữ liệu độ cao trong các nghiên cứu lũ lụt trước đây tại Đà Nẵng (Ho et al., 2012; Do et al., 2014).
  4. Xây dựng và xác thực một mô hình phân vùng nguy cơ lũ lụt đa tham số dựa trên AHP sử dụng DEM độ phân giải cao và dữ liệu vệ tinh về các sự kiện gần lũ. Mô hình này đã đạt được độ chính xác cao, với gần 84% các khu vực ngập lụt thực tế thuộc vùng có nguy cơ lũ cao và rất cao (Abstract, p. 10).
  5. Chứng minh hiệu quả của ràng buộc bất đẳng thức trong nội suy spline, cho thấy cách nó có thể kiểm soát hiệu quả các lỗi nội suy (undershoot/overshoot) ở các khu vực có dữ liệu điểm thưa thớt hoặc sự gián đoạn địa hình rõ rệt (Chapter 3, p. 45).

Nghiên cứu này thúc đẩy sự hiểu biết về mô hình địa hình và rủi ro lũ lụt, đại diện cho một tiến bộ mô hình (paradigm advancement) từ các cách tiếp cận tổng quát sang các giải pháp địa không gian tùy chỉnh, chính xác cao và dựa trên bằng chứng cho các vấn đề môi trường phức tạp. Các phát hiện và phương pháp luận mở ra ít nhất 3 dòng nghiên cứu mới, bao gồm mô hình hóa lũ lụt động lực học, tích hợp học máy với các hệ thống hỗ trợ ra quyết định địa không gian và nghiên cứu sâu hơn về các ràng buộc bất đẳng thức trong tạo bề mặt. Với khả năng áp dụng cho các vùng đồng bằng ven biển tương tự trên toàn cầu, luận án có liên quan quốc tế sâu sắc, cung cấp một khuôn mẫu cho quản lý rủi ro lũ lụt ở các quốc gia đang phát triển. Di sản của nó là những kết quả có thể đo lường được, bao gồm việc nâng cao độ chính xác của dữ liệu địa hình và các công cụ lập kế hoạch dựa trên bằng chứng, có khả năng giảm thiểu đáng kể thiệt hại do lũ lụt và thúc đẩy khả năng chống chịu của cộng đồng ven biển trên toàn thế giới.