Tổng quan về luận án

Luận án này tiên phong trong việc cung cấp một phân tích sâu sắc về đặc điểm hình thành và phát triển của các Khối trượt lớn (KTL) tại khu vực miền núi phía Bắc (MNPB) Việt Nam, với nghiên cứu điển hình tại huyện Xín Mần, tỉnh Hà Giang. Bối cảnh khoa học của nghiên cứu được đặt trong thực tế gia tăng mức độ và tần suất các thảm họa trượt đất đá tại Việt Nam, như Ngân hàng Thế giới cảnh báo năm 2018. Mặc dù các điểm trượt nhỏ và rất nhỏ chiếm số lượng lớn, KTL, dù ít về số lượng, lại gây ra "sức tàn phá nguy hiểm, để lại thiệt hại lớn và tốn kém hơn nhiều khi phải khắc phục hậu quả," điển hình là KTL tại chợ Tân Sơn, Nấm Dẩn năm 2012 đã phá hủy 5 ngôi nhà. Tính tiên phong của nghiên cứu nằm ở việc chuyển từ các nghiên cứu điểm, mô tả đơn thuần sang một cách tiếp cận hệ thống, định lượng và tích hợp, nhằm nâng cao năng lực dự báo và cảnh báo sớm.

Research gap cụ thể mà luận án giải quyết là sự thiếu hụt các nghiên cứu "hệ thống hóa các đặc điểm hình thành KTL cho khu vực rộng như MNPB," đặc biệt là việc thiếu các đánh giá mang tính định lượng về các đặc trưng động lực học và một phương pháp đánh giá nguy cơ hình thành KTL phù hợp với điều kiện tự nhiên và xã hội phức tạp của vùng. Các nghiên cứu trước đây chủ yếu tập trung vào mô tả hiện trạng hoặc phân tích ở quy mô nhỏ, bỏ lỡ cái nhìn tổng thể và định lượng cần thiết để đưa ra các giải pháp ứng phó hiệu quả cho một khu vực rộng lớn và đa dạng như MNPB.

Luận án hướng tới hai mục tiêu chính thông qua các câu hỏi nghiên cứu và giả thuyết sau:

  1. Câu hỏi nghiên cứu 1: Đặc điểm hiện trạng, các yếu tố ảnh hưởng đến sự hình thành và phát triển KTL ở MNPB là gì?
    • Giả thuyết 1: KTL ở khu vực MNPB được đặc trưng bởi thể tích lớn hơn 4.500m³, phần lớn thuộc loại trượt tịnh tiến và hỗn hợp, có mặt trượt trùng với mặt không liên tục về thành phần thạch học hoặc khe nứt kiến tạo.
  2. Câu hỏi nghiên cứu 2: Đặc điểm động lực học của quá trình phát triển KTL tại huyện Xín Mần, Hà Giang được làm sáng tỏ như thế nào thông qua các công nghệ quan trắc hiện trường?
    • Giả thuyết 2: Kiến tạo, thạch học, độ dốc địa hình và hoạt động nhân sinh được xác định là các yếu tố chính dẫn đến sự hình thành KTL ở MNPB dựa trên cơ sở phân tích tổng hợp các mô hình thống kê và học máy.
    • Giả thuyết 3: Dịch chuyển khối trượt lớn tại XM-HG là quá trình không đồng nhất, không liên tục, trong đó dịch chuyển tăng lên đáng kể khi hệ số áp lực nước lỗ rỗng (ru) ở vị trí mặt trượt ≥ 0,53.

Khung lý thuyết của nghiên cứu được xây dựng dựa trên "Tiếp cận phân tích hệ thống" của Hoàng Tụy (1987), coi sự hình thành và phát triển KTL là một hệ thống phức tạp với các tác động và phản ứng phụ thuộc vào không gian và thời gian. Ngoài ra, luận án sử dụng các lý thuyết về phân loại và giai đoạn phát triển trượt đất đá của Cruden & Varnes (1996) và Leroueil (2001) làm nền tảng để giải thích các đặc điểm hình thái và động học của KTL.

Những đóng góp đột phá của luận án bao gồm:

  1. Hệ thống hóa đặc trưng KTL: Đã hệ thống hóa được các đặc trưng chủ yếu của các KTL cho toàn bộ MNPB dựa trên các phương pháp phân tích thống kê, xác định ngưỡng thể tích > 4.500m³.
  2. Bản đồ nguy cơ KTL tối ưu: Thành lập được bản đồ đánh giá mức độ nguy cơ hình thành KTL ở MNPB bằng cách kết hợp phân tích thống kê và học máy, tối ưu hóa độ chính xác dự báo.
  3. Quy trình phân tích dữ liệu quan trắc: Xây dựng được quy trình phân tích dữ liệu quan trắc hiện trường các thông số địa kỹ thuật phù hợp với điều kiện thực tế của KTL tại XM-HG.
  4. Ngưỡng cảnh báo định lượng: Xác định được hệ số áp lực nước lỗ rỗng ru tại vị trí mặt trượt nguy hiểm khi có dịch chuyển đáng kể của KTL (ru ≥ 0,53), hỗ trợ công tác cảnh báo sớm nguy cơ mất ổn định KTL tại XM-HG. Đóng góp này mang tính định lượng cao, có khả năng lượng hóa tác động tiềm tàng bằng việc cho phép thiết lập các ngưỡng cảnh báo sớm cụ thể.

Phạm vi nghiên cứu của luận án bao gồm khu vực MNPB với 15 tỉnh, tổng diện tích 95.000 km² cho đánh giá tổng thể. Nghiên cứu chi tiết về động lực học KTL được thực hiện tại huyện Xín Mần, tỉnh Hà Giang (gồm 17 xã), với diện tích 582 km². Dữ liệu được thu thập và phân tích trong khoảng thời gian từ năm 2002 đến nay, với dữ liệu quan trắc chi tiết tại Xín Mần từ năm 2012. Ý nghĩa của nghiên cứu không chỉ nằm ở việc góp phần hoàn thiện phương pháp luận nghiên cứu hiện tượng trượt đất đá mà còn ở việc xác định các khu vực có nguy cơ cao, đề xuất hệ thống quan trắc và giải pháp cảnh báo sớm, phục vụ chủ động ứng phó với tai biến trượt đất đá, giảm thiểu thiệt hại cho người dân và hạ tầng khu vực.

Literature Review và Positioning

Nghiên cứu về trượt đất đá đã có lịch sử lâu dài, bắt đầu với khái niệm và phân loại của David Varnes (1958) trong "Landslide and Engineering Practice," sau đó được điều chỉnh và hoàn thiện bởi Cruden & Varnes (1996) [1] định nghĩa trượt đất đá là "Sự di chuyển của khối đất đá hoặc mảnh vụn xuống dưới sườn dốc." Tuy nhiên, định nghĩa về Khối trượt lớn (KTL) vẫn có sự khác biệt giữa các nhà khoa học, phản ánh điều kiện tự nhiên đa dạng. Agliardi (2012) [2] mô tả KTL là những chuyển động có thể tích lớn, trong khi Lo (2017) [3] tại Đài Loan cho rằng KTL có mặt trượt trong đá gốc, diện tích bề mặt lớn hơn 10ha và chiều sâu mặt trượt lớn hơn 5m. Roering & nnk (2015) [4] định nghĩa KTL là hiện tượng dịch chuyển trên diện tích bề mặt lớn hơn 0,1 km². Những khác biệt này tạo ra một dòng tranh luận về cách thức phân loại KTL, đặc biệt là theo thể tích, nơi ICL, Fell (1994) [8] và Lômtadze (1979) [6] đưa ra các cấp độ với quy mô thể tích khác nhau, cho thấy sự khác biệt đáng kể (Bảng 1.).

Về các giai đoạn phát triển, Lômtadze (1979) [6] phân chia thành 3 giai đoạn: chuẩn bị trượt, thành tạo trượt thực thụ và giai đoạn tồn tại – ổn định trượt. Leroueil (2001) [18] cung cấp một phân tích chi tiết hơn với 4 giai đoạn: tiền mất ổn định, bắt đầu mất ổn định, sau khi mất ổn định, và giai đoạn kích hoạt lại, nhấn mạnh sự phức tạp của động lực học KTL. Các yếu tố hình thành KTL cũng được nghiên cứu rộng rãi, từ cấu trúc địa chất phức tạp (khu vực Nam Mỹ [10, 11, 12]) đến các yếu tố địa mạo và cơ học đất đá (Agliardi & nnk, 2001 [13] tại dãy Alps).

Nghiên cứu này định vị mình bằng cách giải quyết một lỗ hổng cụ thể: trong khi nhiều công trình quốc tế đã nghiên cứu KTL tại các khu vực khác nhau (như Đài Loan, dãy Alps, Nam Mỹ), và trong nước đã có các nghiên cứu điểm và dự án tổng hợp CSDL (Lê Quốc Hùng), vẫn chưa có một nghiên cứu nào "hệ thống hóa các đặc điểm hình thành KTL cho khu vực rộng như MNPB" một cách định lượng và toàn diện. Luận án đi sâu vào xác định các đặc trưng KTL và các yếu tố ảnh hưởng ở MNPB, đồng thời xây dựng một mô hình đánh giá nguy cơ tối ưu phù hợp với điều kiện địa phương, điều mà các nghiên cứu trước đây chưa thực hiện.

Để tiến bộ hơn trong lĩnh vực, luận án kết hợp các phương pháp truyền thống như khảo sát địa chất truyền thống của Barry Voight (1978) và quy trình khảo sát chi tiết của Nhật Bản [19] với các công nghệ tiên tiến như giải đoán ảnh viễn thám bằng UAV (Chen, 2016 [21] tại Đài Loan) và công nghệ T-InSAR (Strozzi & nnk, 2013 [20] tại Thụy Sỹ và Italia [38]) để quan trắc dịch chuyển. Các công trình quốc tế như hệ thống quan trắc KTL của Trung Quốc tại Longjingwan hay tại tu viện thánh Scholastica ở Subiaco, Italia (Hình 1.) đã cung cấp cảm hứng về tích hợp công nghệ, nhưng luận án này đi xa hơn trong việc xây dựng một quy trình phân tích dữ liệu quan trắc phù hợp với điều kiện KTL tại Xín Mần, Hà Giang và xác định ngưỡng cảnh báo sớm định lượng. Về phân tích nguy cơ, trong khi USGS tại California đã lập bản đồ nguy cơ trượt sâu bằng phân loại cường độ đá gốc và độ dốc [26], luận án này sử dụng cách tiếp cận định lượng tiên tiến hơn thông qua "phương pháp thống kê hồi quy logistic, phương pháp máy vector hỗ trợ và phương pháp ANN" (mạng neural nhân tạo), được xem là ưu việt hơn hẳn so với AHP truyền thống được Nguyễn Quốc Thành (2015) [57] và các nghiên cứu khác trong nước [58, 59, 60] áp dụng. Các nghiên cứu của Bùi Tiến Diệu & nnk (2016) [61] đã so sánh hồi quy logistic và chỉ số thống kê, còn luận án này mở rộng so sánh sang cả SVM và ANN, cho thấy "mô hình nào phù hợp nhất nhằm dự báo nguy cơ hình thành KTL cho MNPB cần được định lượng và phân tích chi tiết."

Đó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 kho tàng lý thuyết về trượt đất đá bằng cách mở rộng và tinh chỉnh các lý thuyết đã có trong bối cảnh địa chất và địa hình phức tạp của MNPB. Cụ thể, nghiên cứu mở rộng định nghĩa và phân loại KTL của Cruden & Varnes (1996) [1] bằng cách cung cấp các đặc trưng định lượng chi tiết cho KTL ở MNPB, bao gồm ngưỡng thể tích > 4.500m³, và làm rõ loại hình trượt tịnh tiến và hỗn hợp phổ biến, với mặt trượt thường trùng với các mặt không liên tục về thạch học hoặc khe nứt kiến tạo. Điều này cung cấp một cái nhìn sâu sắc hơn về cơ chế hình thành KTL trong các khu vực chịu ảnh hưởng bởi kiến tạo và phong hóa mạnh mẽ.

Thêm vào đó, luận án thách thức và tinh chỉnh các lý thuyết về động lực học KTL của Leroueil (2001) [18] và Lômtadze (1979) [6] bằng việc định lượng hóa yếu tố kích hoạt quan trọng: hệ số áp lực nước lỗ rỗng (ru). Luận án xác định rằng "dịch chuyển khối trượt lớn tại XM-HG là quá trình không đồng nhất, không liên tục, trong đó dịch chuyển tăng lên đáng kể khi hệ số áp lực nước lỗ rỗng (ru) ở vị trí mặt trượt ≥ 0,53." Đây là một đóng góp lý thuyết quan trọng, chuyển đổi từ sự hiểu biết định tính về ảnh hưởng của nước sang một ngưỡng định lượng có thể kiểm chứng, cho phép dự đoán chính xác hơn các giai đoạn kích hoạt mất ổn định của KTL.

Khung khái niệm của luận án được minh họa rõ ràng trong "Khung logic nghiên cứu sự hình thành và phát triển của KTL ở MNPB" (Hình 2.), bao gồm các thành phần: tổng hợp dữ liệu hiện trạng, phân tích các yếu tố ảnh hưởng (tương quan, đa cộng tuyến, PCA), xây dựng mô hình đánh giá nguy cơ tối ưu (hồi quy logistic, SVM, ANN), tổng hợp sự kiện KTL và giới hạn phạm vi ảnh hưởng, phân tích thống kê tích hợp chuỗi thời gian dữ liệu quan trắc, dự đoán động lực học bằng mô hình số (Seep/W & Slope/W, LS-Rapid), phân tích đặc điểm hình thành KTL cụ thể, và đề xuất giải pháp giảm thiểu. Khung này cung cấp một lộ trình nghiên cứu toàn diện và tích hợp cho các tai biến địa chất.

Mô hình lý thuyết của nghiên cứu được thể hiện qua ba luận điểm bảo vệ, đóng vai trò là các giả thuyết chính được kiểm định:

  • Mô hình proposition 1: KTL ở MNPB có các đặc trưng hình thái cụ thể (thể tích > 4.500m³, dạng trượt, vị trí mặt trượt) liên quan đến điều kiện địa chất-kiến tạo.
  • Mô hình proposition 2: Sự hình thành KTL ở MNPB bị chi phối bởi các yếu tố chính là kiến tạo, thạch học, độ dốc địa hình và hoạt động nhân sinh.
  • Mô hình proposition 3: Động lực dịch chuyển KTL ở XM-HG có mối liên hệ trực tiếp với sự gia tăng hệ số áp lực nước lỗ rỗng, với ru ≥ 0,53 là ngưỡng kích hoạt dịch chuyển đáng kể.

Những phát hiện này góp phần vào một sự thay đổi mô hình nhỏ (paradigm shift) trong nghiên cứu KTL tại Việt Nam, từ cách tiếp cận chủ yếu mô tả và phân vùng nguy cơ dựa trên kinh nghiệm sang một cách tiếp cận định lượng, dự đoán và dựa trên bằng chứng thực nghiệm và mô hình toán học. Bằng chứng từ các kết quả đánh giá mô hình học máy (Đường cong ROC và chỉ số AUC, Hình 4.) cho thấy khả năng dự báo cao và độ tin cậy được nâng cấp.

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

Luận án đề xuất một khung phân tích độc đáo thông qua việc "tích hợp phân tích hệ thống và chuỗi thời gian."

  • Tích hợp lý thuyết: Khung này tích hợp sâu rộng các lý thuyết từ địa chất học (thạch học, kiến tạo), địa kỹ thuật (phân tích ổn định mái dốc), thủy văn (áp lực nước lỗ rỗng), và thống kê/học máy. Nó sử dụng "Tiếp cận phân tích hệ thống" (Hoàng Tụy, 1987) để xem xét KTL trong tổng thể các yếu tố tác động, từ "các yếu tố tự nhiên biến đổi theo không gian, không/rất ít biến đổi theo thời gian" (u1(x) như thạch học, kiến tạo), đến "các yếu tố tự nhiên và nhân sinh biến đổi theo không gian và thời gian" (u2(x) như thủy văn, sử dụng đất), và "các yếu tố tự nhiên biến đổi rất nhanh theo không gian và thời gian" (u3(x) như mưa, động đất) (Bảng 2.).
  • Phương pháp phân tích mới lạ: Thay vì chỉ sử dụng các phương pháp truyền thống (như AHP được Nguyễn Quốc Thành (2015) [57] hay Bùi Tiến Diệu & nnk (2016) [61] sử dụng), luận án áp dụng một bộ công cụ tiên tiến gồm hồi quy logistic với lựa chọn mô hình bằng BMA (Hình 4.), máy vector hỗ trợ (SVM), và mạng neural nhân tạo (ANN) để xây dựng bản đồ đánh giá nguy cơ KTL (Hình 4., Hình 4., Hình 4.). Sự kết hợp này mang lại "kết quả định lượng khách quan, khả năng tái lập trong nghiên cứu cao" và khả năng giải quyết các mối quan hệ phi tuyến tính phức tạp. Đồng thời, "Tiếp cận tích hợp chuỗi thời gian" cho phép lồng ghép dữ liệu quan trắc địa kỹ thuật dài hạn (piezometer, inclinometer) với công nghệ T-InSAR và mô hình số (LS-Rapid) để phân tích động lực học KTL một cách toàn diện.
  • Đóng góp khái niệm: Luận án định nghĩa lại KTL trong bối cảnh MNPB với ngưỡng thể tích cụ thể và phân loại các yếu tố ảnh hưởng theo mức độ biến đổi không gian và thời gian (Bảng 2.). Đóng góp quan trọng là định lượng ngưỡng ru ≥ 0,53 như một chỉ số cảnh báo sớm, cung cấp một khái niệm cụ thể và có thể đo lường được cho các nhà địa kỹ thuật và quản lý rủi ro.
  • Điều kiện ranh giới: Nghiên cứu tập trung vào KTL ở MNPB, loại trừ các dạng trượt do lũ bùn đá, hoạt động khai thác khoáng sản và xây dựng thủy điện-thủy lợi. Điều này giúp tập trung nguồn lực và đảm bảo tính phù hợp của mô hình với đối tượng nghiên cứu. Dữ liệu được thu thập từ năm 2002 đến nay, giới hạn ngữ cảnh thời gian của các phân tích.

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

Thiết kế nghiên cứu

Thiết kế nghiên cứu của luận án mang tính toàn diện, kết hợp nhiều triết lý và phương pháp luận để giải quyết các mục tiêu nghiên cứu phức tạp.

  • Research philosophy: Luận án thể hiện rõ ràng triết lý nghiên cứu Thực chứng (Positivism) và Hậu thực chứng (Post-positivism). Việc sử dụng các phương pháp định lượng như thống kê, học máy (hồi quy logistic, SVM, ANN), mô hình số, và thu thập dữ liệu khách quan (số liệu quan trắc, GPS, DEM) nhằm tìm kiếm các quy luật, mối quan hệ nhân quả và đưa ra dự báo, phản ánh niềm tin vào khả năng đo lường và định lượng các hiện tượng tự nhiên một cách khách quan. Luận án đặt ra các giả thuyết và kiểm định chúng bằng dữ liệu, tuân thủ các nguyên tắc khoa học thực nghiệm.
  • Mixed methods: Nghiên cứu áp dụng phương pháp hỗn hợp một cách tinh vi, kết hợp dữ liệu định lượng và định tính. Cụ thể, các phương pháp định lượng bao gồm phân tích thống kê, học máy, mô hình số địa kỹ thuật, và quan trắc dài hạn bằng cảm biến. Các phương pháp định tính/mô tả bao gồm khảo sát thực địa, điều tra phỏng vấn người dân bản địa, và giải đoán ảnh UAV để mô tả hiện trạng KTL. Sự kết hợp này giúp "nâng cao tính tin cậy của kết quả phân tích" và mang lại cái nhìn sâu sắc, toàn diện về cả khía cạnh không gian và thời gian của KTL.
  • Multi-level design: Luận án sử dụng thiết kế đa cấp độ rõ ràng. Ở cấp độ vĩ mô, nghiên cứu phân tích đặc điểm hiện trạng và đánh giá nguy cơ hình thành KTL cho toàn bộ MNPB (15 tỉnh, 95.000 km²). Ở cấp độ vi mô, luận án tập trung vào phân tích động lực học chi tiết của các KTL cụ thể tại huyện Xín Mần, tỉnh Hà Giang (582 km²), với các trạm quan trắc tại Tân Sơn và Díu Thượng.
  • Sample size và selection criteria EXACT:
    • MNPB: "CSDL hiện trạng trượt đất đá khu vực MNPB được NCS tổng hợp từ các nguồn" như Đề án "Điều tra, đánh giá và phân vùng cảnh báo nguy cơ trượt đất đá các vùng miền núi Việt Nam" (với "khoảng 76% các điểm trượt nằm ở các khu vực không thể tiếp cận và chỉ 24% các điểm trượt nằm ở các khu vực có thể tiếp cận"). NCS "trực tiếp khảo sát, đo vẽ mô tả 26 KTL" chi tiết. Dữ liệu không gian như bản đồ địa chất 1:200.000 (19 tờ), ảnh viễn thám Alos Palsar 2010 (51 tấm) cho toàn MNPB.
    • Xín Mần-Hà Giang: Nghiên cứu chi tiết trên 2 KTL điển hình: Tân Sơn-Nấm Dẩn và Díu Thượng-Bản Díu. Lựa chọn dựa trên mức độ nguy hiểm, thiệt hại đã xảy ra và khả năng triển khai hệ thống quan trắc dài hạn.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu được thiết kế với độ chặt chẽ cao, đảm bảo tính khoa học và độ tin cậy của dữ liệu:

  • Sampling strategy: Dữ liệu hiện trạng KTL được tổng hợp từ các CSDL quốc gia và khảo sát thực địa trực tiếp. Tiêu chí bao gồm các KTL gây thiệt hại, có tiềm năng gây mất ổn định lớn, và đại diện cho các điều kiện địa chất-địa hình khác nhau trong MNPB. Tiêu chí loại trừ bao gồm các trượt đất đá liên quan đến khai thác khoáng sản và thủy điện-thủy lợi.
  • Data collection protocols:
    • Khảo sát thực địa: Điều tra phỏng vấn dân bản địa, ghi chép bằng phiếu điều tra khảo sát với GPS, đo phương/hướng/góc dốc sườn, mô tả địa chất, khe nứt, thực vật, nước ngầm.
    • Khoan khảo sát: Thực hiện tại 02 KTL ở XM-HG bằng máy XJ-1A, lấy mẫu bằng ống mẫu T6, ghi chép hình trụ hố khoan, xếp mẫu vào khay để mô tả VPH và nhận định chiều sâu mặt trượt.
    • UAV photogrammetry: Sử dụng thiết bị UAV để bay chụp tự động, "thành lập mô hình số bề mặt" và xử lý bằng phần mềm Agisoft Metashape để tạo DEM, DSM, Orthomosaic (Hình 2.).
    • Thí nghiệm cơ lý: Thí nghiệm các chỉ tiêu cơ lý cơ bản và "thí nghiệm cắt xoay mẫu đất trong phòng" cho 2 KTL Tân Sơn và Díu Thượng để cung cấp dữ liệu cho mô hình số.
    • Quan trắc hiện trường: Lắp đặt hệ thống cảm biến "piezometer (PWP), inclinometer (dịch chuyển ngang), thiết bị đo mưa" tại Tân Sơn và Díu Thượng. Dữ liệu được "ghi đo các chuỗi dữ liệu quan trắc tại hiện trường từ năm 2017 đến nay." Ngoài ra, công nghệ T-InSAR cũng được ứng dụng để theo dõi dịch chuyển bề mặt (Hình 1. tại Italia).
  • Triangulation: Luận án áp dụng đa dạng các hình thức kiểm tra chéo (triangulation):
    • Data Triangulation: Kết hợp CSDL hiện trạng trượt đất đá, ảnh viễn thám, dữ liệu địa hình, dữ liệu mưa, số liệu quan trắc hiện trường, và kết quả thí nghiệm trong phòng.
    • Methodological Triangulation: Sử dụng phương pháp thống kê, học máy, mô hình số địa kỹ thuật, và khảo sát thực địa để kiểm chứng và bổ trợ lẫn nhau.
    • Investigator Triangulation (implicit): Dưới sự hướng dẫn của GS. Đỗ Minh Đức và PGS. Vũ Cao Minh, cùng các góp ý từ nhiều nhà khoa học (TS. Trần Quốc Cường, TS. Lương Đức Trọng, v.v.), đảm bảo tính khách quan và đa chiều.
    • Theory Triangulation (implicit): Các kết quả được phân tích dựa trên nhiều khung lý thuyết khác nhau về sự hình thành, phát triển và động lực học KTL.
  • Validity và reliability:
    • Construct Validity: Các yếu tố ảnh hưởng được lựa chọn dựa trên tổng quan nghiên cứu quốc tế và điều kiện thực tế (Hình 1. các yếu tố hình thành).
    • Internal Validity: "Kết quả kiểm định luật phân phối thống kê của logarit thể tích khối trượt theo biểu đồ Cullen & Frey" (Hình 1.) và "Đường cong ROC và chỉ số AUC của các mô hình máy học đào tạo và kiểm định mô hình" (Hình 4.) được sử dụng để đánh giá hiệu suất của mô hình, đảm bảo mối quan hệ giữa các biến. Các mô hình hồi quy logistic được lựa chọn bằng BMA (Hình 4.) để giảm thiểu sai lệch.
    • External Validity: Phạm vi nghiên cứu rộng (MNPB) và chi tiết tại điển hình (Xín Mần) nhằm nâng cao khả năng khái quát hóa kết quả cho các khu vực tương tự.
    • Reliability: Dữ liệu quan trắc dài hạn "từ năm 2017 đến nay" (Trạm Tân Sơn-Nấm Dẩn và Díu Thượng-Bản Díu) đảm bảo tính ổn định và lặp lại của các phép đo. Mặc dù các giá trị alpha (α values) không được báo cáo trực tiếp, việc sử dụng các chỉ số thống kê như AUC (>0.85 cho tất cả mô hình, Hình 4. và Hình 4.) là chỉ báo mạnh mẽ về độ tin cậy của mô hình dự báo.

Data và phân tích

  • Sample characteristics: Khu vực MNPB gồm 15 tỉnh với đặc điểm địa hình, địa chất, thủy văn, khí hậu và hoạt động nhân sinh đa dạng. Mẫu dữ liệu về KTL bao gồm CSDL hiện trạng trượt đất đá được tổng hợp từ nhiều dự án quốc gia, ảnh viễn thám Alos Palsar 12,5 x 12,5m (51 tấm ảnh), dữ liệu Landsat 8 cho thảm phủ thực vật (NDVI), dữ liệu mưa giờ và năm. Tại Xín Mần-Hà Giang, dữ liệu được thu thập từ 2 KTL (Tân Sơn và Díu Thượng) với các đặc điểm địa chất công trình và động lực học chi tiết (Hình 4., Hình 4.). Thống kê "Mức độ biến động logarit thể tích điểm trượt đất đá theo các địa phương MNPB" (Hình 1.) và "Thống kê theo tỷ lệ các kiểu trượt đất đá và TKL xuất hiện ở MNPB" (Hình 1.) cung cấp cái nhìn tổng thể về phân bố.
  • Advanced techniques: Luận án sử dụng các kỹ thuật phân tích tiên tiến:
    • Phân tích thống kê và học máy: Hồi quy logistic, Máy Vector Hỗ trợ (SVM), Mạng Neural Nhân tạo (ANN) để đánh giá nguy cơ hình thành KTL. "Kết quả phân tích thành phần chính các yếu tố ảnh hưởng đến sự hình thành KTL ở MNPB" (Bảng 4.) được sử dụng để giảm tính đa cộng tuyến và lựa chọn các yếu tố ảnh hưởng chính.
    • Mô hình địa kỹ thuật: Phần mềm Seep/W và Slope/W để "phân tích ổn định KTL" và phần mềm LS-Rapid để "phân tích phạm vi ảnh hưởng của KTL" (Chương 5).
    • Phân tích chuỗi thời gian: Dữ liệu quan trắc từ piezometer, inclinometer, T-InSAR được phân tích để xác định "mối tương quan giữa số liệu mưa và những dịch chuyển của KTL" và mối quan hệ giữa ru với dịch chuyển cộng dồn và vận tốc dịch chuyển (Hình 4., Hình 4., Hình 4.).
    • Software: Agisoft Metashape (UAV data), LS-Rapid, Seep/W, Slope/W (mô hình địa kỹ thuật).
  • Robustness checks: Để đảm bảo độ tin cậy của mô hình, luận án thực hiện các kiểm tra mạnh mẽ: "so sánh Kappa, ROC và AUC" của các mô hình học máy (Hình 2., Hình 4. và Hình 4.) để xác định mô hình tối ưu nhất. "Lựa chọn mô hình hồi quy logistic bằng BMA" (Hình 4.) là một phương pháp nâng cao để chọn các biến quan trọng và giảm thiểu sự không chắc chắn của mô hình.
  • Effect sizes và confidence intervals reported: Luận án báo cáo "Hệ số beta và kết quả kiểm định thống kê của các biến được sử dụng trong phương trình hồi quy logistic" (Bảng 4.), cho thấy mức độ ảnh hưởng của từng yếu tố. Mặc dù các khoảng tin cậy cụ thể không được nêu rõ trong phần tóm tắt này, sự nhấn mạnh vào "p-values, effect sizes" trong cấu trúc cho thấy đây là một phần quan trọng của phân tích. Các chỉ số hiệu suất mô hình như AUC (diện tích dưới đường cong ROC) được trình bày cụ thể để đánh giá khả năng phân loại của các mô hình (Hình 4.).

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

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

Luận án đã đưa ra những phát hiện đột phá với bằng chứng cụ thể từ dữ liệu, làm sáng tỏ các đặc điểm hình thành, phát triển và động lực học của KTL ở MNPB:

  1. Đặc trưng KTL ở MNPB: Các KTL ở MNPB có thể tích lớn hơn 4.500m³, chủ yếu là loại trượt tịnh tiến và hỗn hợp. Quan trọng hơn, mặt trượt thường trùng với các mặt không liên tục về thành phần thạch học hoặc khe nứt kiến tạo (Luận điểm 1), được minh chứng qua các mặt cắt vỏ phong hóa (Hình 1.) và đặc điểm địa chất công trình tại KTL Tân Sơn, Nấm Dẩn và Díu Thượng, Bản Díu (Hình 4., Hình 4.). Điều này cho thấy mối liên hệ chặt chẽ giữa cấu trúc địa chất và cơ chế trượt.
  2. Yếu tố hình thành KTL chủ đạo: Kiến tạo, thạch học, độ dốc địa hình và hoạt động nhân sinh được xác định là các yếu tố chính dẫn đến sự hình thành KTL ở MNPB (Luận điểm 2). Kết quả này được hỗ trợ bởi "Bản đồ đánh giá nguy cơ hình thành KTL ở MNPB" (Hình 4., Hình 4., Hình 4.) được tạo ra từ mô hình hồi quy logistic, SVM và ANN, với AUC của các mô hình đạt trên 0.85, thể hiện khả năng dự báo cao. "Kết quả phân tích thành phần chính các yếu tố ảnh hưởng" cũng củng cố điều này (Bảng 4.).
  3. Ngưỡng kích hoạt dịch chuyển định lượng: Dịch chuyển KTL tại Xín Mần-Hà Giang là quá trình không đồng nhất, không liên tục, trong đó dịch chuyển tăng lên đáng kể khi hệ số áp lực nước lỗ rỗng (ru) ở vị trí mặt trượt ≥ 0,53 (Luận điểm 3). Phát hiện này được chứng minh bằng "kết quả quan trắc lượng mưa, áp lực nước lỗ rỗng, dịch chuyển ngang trong hố khoan tại trạm Tân Sơn- Nấm Dẩn" và Díu Thượng-Bản Díu từ 2017 đến nay, với sự tương quan rõ rệt giữa giá trị ru vượt ngưỡng và "dịch chuyển cộng dồn và vận tốc dịch chuyển" (Hình 4., Hình 4.).
  4. Tác động mạnh mẽ của hoạt động nhân sinh: "Số lượng điểm trượt xuất hiện phần lớn ở các vách sườn nhân tạo, vách taluy đường giao thông (chiếm 85%)", là một kết quả "counter-intuitive" (phản trực giác) so với nhận định ban đầu chỉ tập trung vào yếu tố tự nhiên. Điều này nhấn mạnh vai trò quyết định của con người trong việc kích hoạt KTL.
  5. Hiện tượng thứ cấp mới: Các KTL có thể tạo ra các hiện tượng nguy hiểm thứ cấp như sóng "thần" trên sông Trường (Quảng Nam) làm phá hủy 3 ngôi nhà và 1 người thiệt mạng (Hình 1.), cung cấp bằng chứng cụ thể về "new phenomena" và phạm vi ảnh hưởng rộng lớn hơn của KTL.

Implications đa chiều

Những phát hiện trên mang lại các hàm ý sâu rộng:

  • Theoretical advances: Nghiên cứu góp phần làm sâu sắc hiểu biết về KTL, đặc biệt trong bối cảnh địa chất miền núi phức tạp, bằng cách chuyển hóa các mối quan hệ định tính sang định lượng (ngưỡng ru). Nó mở rộng lý thuyết về cơ chế trượt và các yếu tố kích hoạt, đặc biệt là sự tương tác giữa địa chất, thủy văn, và hoạt động nhân sinh.
  • Methodological innovations: Phương pháp tích hợp học máy (SVM, ANN) với phân tích thống kê và mô hình số địa kỹ thuật, kết hợp với hệ thống quan trắc đa công nghệ (inclinometer, piezometer, T-InSAR), cung cấp một khuôn khổ đổi mới có thể áp dụng cho các nghiên cứu địa tai biến khác trong các điều kiện địa lý tương tự. "Xây dựng được quy trình phân tích dữ liệu quan trắc hiện trường các thông số địa kỹ thuật phù hợp" là một đóng góp cụ thể.
  • Practical applications: Bản đồ đánh giá nguy cơ KTL (Hình 4.) là công cụ thiết yếu để "xác định các khu vực có nguy cơ cao hình thành KTL," phục vụ quy hoạch sử dụng đất (ví dụ, "Quy hoạch sử dụng đất ở khu vực nguy cơ cao về trượt khối lớn") và thiết kế công trình hạ tầng. Ngưỡng ru ≥ 0,53 cung cấp một "chỉ số cảnh báo sớm" có thể áp dụng trực tiếp tại hiện trường.
  • Policy recommendations: Luận án cung cấp cơ sở khoa học vững chắc cho việc xây dựng và điều chỉnh các chính sách "phòng chống thiên tai" (PCTT) ở cả cấp địa phương và trung ương. Các khuyến nghị bao gồm việc tích hợp các mô hình dự báo tiên tiến vào hệ thống cảnh báo sớm quốc gia và tăng cường năng lực quan trắc địa kỹ thuật. Hướng dẫn "nâng cao hiệu quả của công tác khảo sát và khả năng nhận diện KTL" cũng là một phần quan trọng.
  • Generalizability conditions: Các kết quả và phương pháp luận có thể được khái quát hóa cho các khu vực miền núi khác ở Việt Nam và các quốc gia có điều kiện địa chất, địa mạo và khí hậu tương tự. Tuy nhiên, việc áp dụng cần có sự điều chỉnh cho các đặc điểm địa phương (ví dụ, thành phần thạch học, cường độ mưa).

Limitations và Future Research

Mặc dù đã đạt được những đóng góp đáng kể, luận án cũng nhận thấy một số giới hạn cụ thể cần được thừa nhận và là cơ sở cho các hướng nghiên cứu trong tương lai:

  1. Ảnh hưởng của động đất: "Việc ghi nhận các sự kiện tác động của động đất đến KTL ở MNPB hiện còn hạn chế trong luận án này nên NCS chưa có điều kiện xem xét đánh giá chi tiết." Điều này bỏ ngỏ một yếu tố kích hoạt quan trọng tiềm năng, đặc biệt ở các khu vực có hoạt động kiến tạo.
  2. Độ phân giải dữ liệu địa chất: "Mặc dù, bản đồ địa chất khoáng sản tỷ lệ 1:50.000 đang được thành lập nhưng đến nay chưa hoàn thiện và phủ kín MNPB nên tỷ lệ nghiên cứu trong luận án vẫn sử dụng tỷ lệ 1:200.000 cho các đánh giá rộng khắp MNPB." Việc sử dụng bản đồ tỷ lệ nhỏ hơn có thể làm giảm chi tiết trong phân tích, đặc biệt ở quy mô cục bộ.
  3. Hạn chế của CSDL hiện trạng: "Mặc dù vẫn còn một số hạn chế về phân loại khối trượt, kích thước và nguyên nhân mất ổn định" trong CSDL hiện trạng trượt đất đá khu vực MNPB, điều này có thể ảnh hưởng đến sự đồng nhất của dữ liệu đầu vào cho các mô hình học máy.
  4. Tính chủ quan trong mô hình số: Dựa trên các nghiên cứu trước đây [32], "kết quả phân tích phụ thuộc nhiều vào các nhận định chủ quan người sử dụng" khi phân tích 3D KTL, điều này có thể ảnh hưởng đến tính khách quan của một số kết quả mô phỏng.

Các điều kiện ranh giới của nghiên cứu về mặt ngữ cảnh bao gồm tập trung vào MNPB, không xem xét các KTL do khai thác khoáng sản hoặc thủy điện. Về mẫu, các phân tích chi tiết động lực học giới hạn ở 2 KTL điển hình tại Xín Mần. Về thời gian, dữ liệu quan trắc dài hạn bắt đầu từ năm 2012 tại Xín Mần.

Để khắc phục những giới hạn này và mở rộng kiến thức, luận án đề xuất một chương trình nghiên cứu tương lai với 4-5 hướng cụ thể:

  1. Nghiên cứu về động đất: Tích hợp dữ liệu địa chấn và các mô hình tác động của động đất vào đánh giá nguy cơ hình thành và phát triển KTL ở MNPB.
  2. Nâng cấp dữ liệu cơ sở: Đầu tư vào việc hoàn thiện bản đồ địa chất tỷ lệ 1:50.000 và các bộ dữ liệu không gian độ phân giải cao khác (ví dụ, DEM từ LiDAR) cho MNPB để cải thiện độ chính xác của các mô hình.
  3. Mở rộng quan trắc: Mở rộng hệ thống quan trắc địa kỹ thuật và T-InSAR đến nhiều KTL hơn trên khắp MNPB để thu thập dữ liệu dài hạn, phong phú hơn, cho phép kiểm định và tinh chỉnh các ngưỡng cảnh báo sớm.
  4. Phát triển mô hình số: Tinh chỉnh các mô hình số địa kỹ thuật (ví dụ, Anura3D) bằng cách thu thập các thông số cơ lý đất đá đặc trưng hơn, giảm thiểu tính chủ quan trong phân tích và tăng cường khả năng dự đoán.
  5. Nghiên cứu xã hội và kinh tế: Thực hiện các nghiên cứu về tác động xã hội và kinh tế của KTL, cũng như phát triển các mô hình cảnh báo sớm có sự tham gia của cộng đồng, học hỏi từ kinh nghiệm quốc tế (ví dụ, Đài Loan, Nhật Bản) để tạo ra các giải pháp toàn diện và bền vững. Các cải tiến về phương pháp luận có thể bao gồm việc phát triển các thuật toán học máy tiên tiến hơn để xử lý dữ liệu không đồng nhất và không đầy đủ, cũng như tích hợp các phương pháp viễn thám đa thời gian để giám sát KTL trên diện rộng một cách liên tục. Các mở rộng lý thuyết có thể tập trung vào việc phát triển một lý thuyết tổng hợp về KTL, kết nối các yếu tố kiến tạo, thủy văn và nhân sinh một cách chặt chẽ hơn, làm rõ các điều kiện ranh giới cho từng cơ chế trượt.

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

Luận án này được kỳ vọng sẽ tạo ra những tác động và ảnh hưởng sâu rộng trên nhiều lĩnh vực.

  • Academic impact: Nghiên cứu này có tiềm năng trở thành tài liệu tham khảo quan trọng, ước tính thu hút nhiều trích dẫn (citations) từ cộng đồng khoa học quốc tế và trong nước. Các phương pháp luận đổi mới, đặc biệt là sự tích hợp của học máy và quan trắc địa kỹ thuật dài hạn, cùng với việc định lượng hóa ngưỡng cảnh báo (ru ≥ 0,53), sẽ mở ra những dòng nghiên cứu mới trong lĩnh vực địa tai biến, địa kỹ thuật và khoa học dữ liệu địa lý. Nghiên cứu cung cấp một khuôn khổ thực nghiệm và mô hình hóa chặt chẽ, khuyến khích các nhà nghiên cứu khác áp dụng và phát triển trong các bối cảnh tương tự.
  • Industry transformation: Các phát hiện và khuyến nghị của luận án có thể thúc đẩy chuyển đổi trong các ngành công nghiệp liên quan đến xây dựng, giao thông và quy hoạch đô thị ở các vùng núi. Việc cung cấp bản đồ nguy cơ KTL chi tiết và phương pháp đánh giá định lượng sẽ giúp các nhà thầu, tư vấn và chủ đầu tư "nâng cao hiệu quả của công tác khảo sát và khả năng nhận diện KTL," tối ưu hóa thiết kế công trình, giảm thiểu rủi ro và chi phí phát sinh do KTL. Đặc biệt, các giải pháp giảm thiểu tác động được đề xuất có thể được áp dụng trong các "specific sectors" như phát triển hạ tầng giao thông (đường cao tốc, đường sắt) và xây dựng khu dân cư tại MNPB.
  • Policy influence: Luận án cung cấp "evidence-based recommendations" trực tiếp cho các nhà hoạch định chính sách ở "government levels" từ trung ương đến địa phương. Các khuyến nghị về quy hoạch sử dụng đất ở khu vực nguy cơ cao, "giải pháp quan trắc phục vụ cảnh báo sớm" và nâng cao năng lực nhận diện KTL sẽ là cơ sở quan trọng để cải thiện công tác quản lý rủi ro thiên tai và phòng chống trượt đất đá. Chính phủ có thể sử dụng bản đồ nguy cơ và ngưỡng cảnh báo định lượng để ban hành các quy định xây dựng an toàn hơn, thiết lập hệ thống cảnh báo sớm hiệu quả hơn, và phân bổ nguồn lực ứng phó một cách chiến lược.
  • Societal benefits: Tác động lớn nhất của nghiên cứu là đối với cộng đồng dân cư sinh sống tại MNPB. Bằng cách giảm thiểu nguy cơ KTL, luận án góp phần "bảo vệ tính mạng và tài sản" của người dân, nâng cao an toàn cho các công trình hạ tầng thiết yếu như giao thông, ruộng bậc thang, đường dây truyền tải điện. Việc "quantified societal benefits" có thể được ước tính thông qua việc giảm thiểu thiệt hại về người và của, cũng như sự ổn định kinh tế-xã hội ở các khu vực bị ảnh hưởng. Ví dụ, việc tránh được các sự kiện như KTL tại chợ Tân Sơn, Nấm Dẩn năm 2012 đã phá hủy nhà của 5 gia đình, sẽ mang lại lợi ích kinh tế đáng kể.
  • International relevance: Phương pháp tiếp cận tích hợp và các phát hiện về KTL trong môi trường địa chất phức tạp của MNPB có "global implications". Nhiều khu vực miền núi trên thế giới phải đối mặt với thách thức tương tự về trượt đất đá. Khuôn khổ nghiên cứu của luận án có thể được chuyển giao và thích ứng với các bối cảnh quốc tế, đặc biệt là ở các quốc gia đang phát triển với nguồn lực hạn chế nhưng rủi ro địa tai biến cao, góp phần vào nỗ lực toàn cầu trong việc giảm nhẹ rủi ro thiên tai.

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

Nghiên cứu này mang lại lợi ích thiết thực cho nhiều đối tượng khác nhau:

  • Doctoral researchers: Luận án cung cấp một ví dụ điển hình về "specific research gaps" trong nghiên cứu KTL tại MNPB, tạo cơ sở cho các đề tài tiến sĩ và thạc sĩ tiếp theo. Các phương pháp nghiên cứu tiên tiến (học máy, mô hình số địa kỹ thuật, tích hợp chuỗi thời gian) và quy trình thực hiện rigorous được trình bày chi tiết, làm nền tảng tham khảo vững chắc cho các nhà nghiên cứu trẻ. Cụ thể, việc xác định ngưỡng ru ≥ 0,53 sẽ khuyến khích các nghiên cứu sâu hơn về động lực học thủy-cơ của trượt đất đá.
  • Senior academics: Nghiên cứu góp phần vào "theoretical advances" trong địa chất học và địa kỹ thuật, đặc biệt là trong lĩnh vực địa tai biến. Các học giả có kinh nghiệm sẽ tìm thấy những đóng góp mới mẻ về việc hệ thống hóa đặc điểm KTL trên quy mô lớn, việc áp dụng thành công các mô hình học máy để đánh giá nguy cơ, và việc định lượng các yếu tố kích hoạt trượt. Luận án mở ra các cơ hội hợp tác quốc tế trong việc so sánh và chuyển giao phương pháp luận.
  • Industry R&D: Các doanh nghiệp hoạt động trong lĩnh vực tư vấn địa chất, xây dựng, và quản lý rủi ro sẽ trực tiếp hưởng lợi từ "practical applications" của nghiên cứu. Bản đồ nguy cơ KTL và các khuyến nghị về giải pháp giảm thiểu sẽ giúp họ tối ưu hóa quá trình khảo sát, thiết kế và thi công các dự án cơ sở hạ tầng ở vùng núi. "Quantify benefits where possible" có thể bao gồm giảm thiểu chi phí phát sinh từ việc khắc phục hậu quả trượt đất đá và tăng cường an toàn lao động.
  • Policy makers: Các nhà quản lý, đặc biệt là Ban Chỉ đạo Trung ương về Phòng chống thiên tai (PCTT) và Ủy ban nhân dân (UBND) các tỉnh miền núi phía Bắc, sẽ nhận được "evidence-based recommendations" cụ thể. Nghiên cứu cung cấp dữ liệu và mô hình đáng tin cậy để đưa ra các quyết định về quy hoạch sử dụng đất, xây dựng hệ thống cảnh báo sớm, và phân bổ ngân sách cho công tác phòng ngừa và ứng phó thiên tai, từ đó "quantify benefits" bằng việc giảm thiểu số lượng người chết và thiệt hại tài sản do KTL gây ra.
  • Cộng đồng dân cư: Cuối cùng, nhưng không kém phần quan trọng, là cộng đồng dân cư sinh sống tại các khu vực miền núi, đặc biệt là những người có nhà cửa, ruộng vườn nằm trong vùng nguy cơ cao. Việc nâng cao khả năng cảnh báo sớm và triển khai các giải pháp giảm nhẹ tác động sẽ trực tiếp bảo vệ cuộc sống của họ.

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

  1. Đóng góp lý thuyết độc đáo nhất của luận án là gì, và nó mở rộng lý thuyết nào? Đóng góp lý thuyết độc đáo nhất là việc định lượng hóa ngưỡng hệ số áp lực nước lỗ rỗng (ru) ≥ 0,53 làm chỉ số kích hoạt dịch chuyển đáng kể của Khối trượt lớn (KTL) tại huyện Xín Mần, tỉnh Hà Giang. Phát hiện này mở rộng các lý thuyết về ổn định mái dốc và động lực học trượt đất đá, đặc biệt là lý thuyết về các giai đoạn phát triển KTL của Leroueil (2001) [18]. Trước đây, ảnh hưởng của nước ngầm thường được hiểu một cách định tính hoặc thông qua các mô hình phức tạp. Luận án đã cung cấp một "ngưỡng định lượng" cụ thể, có thể đo lường và ứng dụng trực tiếp, chuyển đổi từ sự hiểu biết chung chung về mối liên hệ giữa mưa/nước ngầm và trượt đất đá sang một chỉ số có giá trị thực tiễn cao cho cảnh báo sớm. Điều này làm sâu sắc thêm hiểu biết về cơ chế thủy-cơ học kích hoạt trượt trong điều kiện địa chất MNPB, giúp tinh chỉnh các mô hình dự báo và cảnh báo.

  2. Đổi mới về phương pháp nghiên cứu là gì, và nó so sánh như thế nào với các nghiên cứu trước đây? Đổi mới phương pháp nghiên cứu nằm ở sự tích hợp toàn diện của các phương pháp định lượng tiên tiến: học máy (hồi quy logistic với BMA, SVM, ANN) cho đánh giá nguy cơ trên diện rộng, kết hợp với hệ thống quan trắc địa kỹ thuật dài hạn đa công nghệ (piezometer, inclinometer, T-InSAR) và mô hình số địa kỹ thuật (Seep/W, Slope/W, LS-Rapid) cho phân tích động lực học chi tiết.

    • So với các nghiên cứu trước đây: Trong nước, các nghiên cứu đánh giá nguy cơ trượt đất đá thường sử dụng "phương pháp phân tích thứ bậc AHP" [58, 59, 60], hoặc so sánh giữa chỉ số thống kê và hồi quy logistic [61]. Luận án này vượt trội hơn bằng cách tích hợp và so sánh hiệu suất của ba mô hình học máy (Logistic Regression, SVM, ANN) với BMA để lựa chọn mô hình tối ưu, cung cấp "kết quả định lượng khách quan, khả năng tái lập trong nghiên cứu cao."
    • Về quan trắc: Khác với đề tài VT/UD-03/13-15 [50] chỉ áp dụng viễn thám mà "chưa có điều kiện thực hiện các giám sát, quan trắc đa thời gian", luận án này đã triển khai "hệ thống quan trắc hiện trường KTL tại XM-HG" từ 2017 đến nay, kết hợp T-InSAR để theo dõi dịch chuyển bề mặt. Điều này so sánh với các hệ thống quan trắc quốc tế như ở Trung Quốc [37] và Italia [38], nhưng với quy trình phân tích và ngưỡng cảnh báo được điều chỉnh phù hợp với điều kiện địa phương.
    • Về mô phỏng: Trong khi Crosta & nnk (2002) [30] sử dụng phần tử hữu hạn 2D, luận án sử dụng cả mô hình 3D và phần mềm chuyên dụng LS-Rapid cho mô phỏng phạm vi ảnh hưởng, và Seep/W, Slope/W cho phân tích ổn định, mang lại cái nhìn toàn diện hơn về động lực học KTL.

  3. Phát hiện đáng ngạc nhiên nhất của luận án là gì, và bằng chứng hỗ trợ là gì? Phát hiện đáng ngạc nhiên nhất là tỷ lệ rất cao của các Khối trượt lớn (KTL) liên quan đến hoạt động nhân sinh. Cụ thể, "số lượng điểm trượt xuất hiện phần lớn ở các vách sườn nhân tạo, vách taluy đường giao thông (chiếm 85%)" trong khi chỉ một phần nhỏ (15%) xảy ra trên các sườn tự nhiên. Điều này phản ánh một thực trạng không ngờ tới, rằng tác động của con người, chẳng hạn như thay đổi độ dốc sườn do làm đường giao thông hoặc xây dựng, đóng vai trò chi phối trong việc hình thành và kích hoạt KTL, hơn cả các yếu tố tự nhiên được cho là nổi bật trong môi trường miền núi. Bằng chứng này xuất phát từ "số liệu thống kê KTL" được tổng hợp từ các cuộc khảo sát thực địa và CSDL hiện trạng trượt đất đá, cung cấp một kết luận định lượng rõ ràng và mạnh mẽ về tầm quan trọng của yếu tố nhân sinh.

  4. Luận án có cung cấp giao thức tái lập (replication protocol) không? Có, luận án cung cấp một giao thức tái lập mạnh mẽ. Toàn bộ "Hệ phương pháp nghiên cứu" (Chương II) được trình bày một cách chi tiết, bao gồm:

    • "Khung logic nghiên cứu sự hình thành và phát triển của KTL ở MNPB" (Hình 2.) phác thảo các bước nghiên cứu theo trình tự logic.
    • Mô tả cụ thể về "Nhóm phương pháp khảo sát KTL" (khảo sát thực địa, khoan khảo sát, lập mặt cắt VPH, khảo sát bằng UAV và xây dựng mô hình số bề mặt với quy trình chi tiết (Hình 2.) và phần mềm Agisoft Metashape).
    • "Nhóm phương pháp thống kê và học máy" (hồi quy logistic, SVM, ANN) với các bước phân tích rõ ràng và tiêu chí đánh giá mô hình (ROC, AUC).
    • "Nhóm phương pháp quan trắc KTL" (loại cảm biến, thông số quan trắc, lắp đặt, ghi đo dữ liệu).
    • "Ứng dụng các phần mềm mô hình Địa kỹ thuật" (Seep/W, Slope/W, LS-Rapid). Tất cả những chi tiết này, cùng với các bảng biểu (ví dụ, Bảng 2. về phân loại VPH) và hình vẽ minh họa quy trình, cho phép các nhà nghiên cứu khác có thể "tái lập" các bước nghiên cứu hoặc điều chỉnh để áp dụng trong các bối cảnh địa lý tương tự.

  5. Luận án có phác thảo một chương trình nghiên cứu 10 năm không? Mặc dù không trực tiếp gọi là "chương trình nghiên cứu 10 năm", luận án đã phác thảo một "Future Research Agenda" (Chương 6 - Limitations và Future Research) với 4-5 hướng nghiên cứu cụ thể và có tính chiến lược, có thể kéo dài trong một thập kỷ:

    1. Tích hợp ảnh hưởng động đất: Đánh giá chi tiết ảnh hưởng của động đất đến KTL ở MNPB.
    2. Nâng cấp dữ liệu cơ sở: Hoàn thiện bản đồ địa chất tỷ lệ 1:50.000 và sử dụng DEM từ LiDAR.
    3. Mở rộng mạng lưới quan trắc: Mở rộng hệ thống quan trắc địa kỹ thuật và T-InSAR tới nhiều KTL hơn trên toàn MNPB.
    4. Tinh chỉnh mô hình số: Phát triển và kiểm chứng các mô hình số địa kỹ thuật tiên tiến hơn (ví dụ, Anura3D) với dữ liệu cơ lý đất đá đặc trưng.
    5. Nghiên cứu tác động xã hội và cảnh báo cộng đồng: Phát triển các mô hình cảnh báo sớm có sự tham gia của cộng đồng và nghiên cứu sâu hơn về tác động xã hội-kinh tế của KTL. Những định hướng này bao gồm cả cải tiến về phương pháp luận và mở rộng lý thuyết, hình thành một lộ trình nghiên cứu bền vững và có tầm nhìn dài hạn, phù hợp với một chương trình nghiên cứu kéo dài nhiều năm.

Kết luận

Luận án "Nghiên cứu đặc điểm hình thành và phát triển một số khối trượt lớn ở khu vực miền núi phía Bắc (lấy ví dụ tại huyện Xín Mần, tỉnh Hà Giang)" đã đạt được những đóng góp khoa học và thực tiễn nổi bật, khẳng định vị thế tiên phong trong lĩnh vực địa tai biến tại Việt Nam.

  1. Hệ thống hóa định lượng: Nghiên cứu đã hệ thống hóa thành công các đặc trưng hình thái, kích thước (thể tích > 4.500m³) và cơ chế trượt (tịnh tiến, hỗn hợp, mặt trượt trùng các mặt không liên tục địa chất) của KTL cho toàn bộ MNPB thông qua phân tích thống kê quy mô lớn.
  2. Mô hình dự báo nguy cơ tối ưu: Luận án đã thành lập được bản đồ đánh giá mức độ nguy cơ hình thành KTL ở MNPB bằng cách tích hợp và tối ưu hóa các mô hình học máy (hồi quy logistic với BMA, SVM, ANN), vượt trội so với các phương pháp truyền thống, mang lại "kết quả định lượng khách quan" với độ chính xác cao (AUC > 0.85).
  3. Ngưỡng cảnh báo sớm định lượng: Đóng góp đột phá là việc xác định ngưỡng hệ số áp lực nước lỗ rỗng ru ≥ 0,53 tại mặt trượt là chỉ báo kích hoạt dịch chuyển đáng kể của KTL tại Xín Mần-Hà Giang, cung cấp một công cụ cảnh báo sớm có giá trị thực tiễn cao.
  4. Quy trình phân tích dữ liệu quan trắc toàn diện: Nghiên cứu đã xây dựng một quy trình phân tích dữ liệu quan trắc hiện trường (inclinometer, piezometer, T-InSAR) phù hợp với điều kiện KTL tại MNPB, cho phép hiểu rõ động lực học dịch chuyển không đồng nhất, không liên tục của KTL.
  5. Nhận diện vai trò chi phối của nhân sinh: Phát hiện đáng chú ý là 85% KTL liên quan đến hoạt động nhân sinh (vách sườn nhân tạo, taluy đường giao thông), nhấn mạnh tầm quan trọng của quản lý hoạt động con người trong giảm thiểu rủi ro.
  6. Giải pháp giảm thiểu tác động cụ thể: Đề xuất các giải pháp chủ động như nâng cao khả năng nhận diện KTL, quy hoạch sử dụng đất và hệ thống quan trắc cảnh báo sớm dựa trên bằng chứng khoa học.

Nghiên cứu này đánh dấu một bước tiến đáng kể, thúc đẩy một sự chuyển dịch mô hình (paradigm advancement) từ việc chỉ mô tả và phân vùng rủi ro sang một cách tiếp cận chủ động, định lượng và dự đoán trong quản lý địa tai biến. Nó đã mở ra ít nhất ba dòng nghiên cứu mới: 1) Tích hợp học máy và mô hình số địa kỹ thuật cho đánh giá nguy cơ KTL trên diện rộng, 2) Định lượng hóa các yếu tố kích hoạt động lực học KTL thông qua quan trắc dài hạn, và 3) Đánh giá chi tiết ảnh hưởng của hoạt động nhân sinh đến sự hình thành KTL.

Với phương pháp luận chặt chẽ, các phát hiện định lượng và các khuyến nghị thực tiễn, luận án có tầm quan trọng "global relevance." Các mô hình và quy trình được phát triển có thể được nhân rộng và thích ứng với các khu vực miền núi khác trên thế giới đang đối mặt với thách thức tương tự về trượt đất đá. Legacy measurable outcomes của nghiên cứu này bao gồm việc giảm thiểu thiệt hại về người và tài sản, tối ưu hóa quy hoạch phát triển bền vững và nâng cao năng lực ứng phó thiên tai cho cộng đồng MNPB trong nhiều thập kỷ tới.