Tổng quan luận án

Luận án tiến sĩ địa chất của tác giả Đào Minh Đức với đề tài "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)", thực hiện tại Học viện Khoa học và Công nghệ – Viện Hàn lâm Khoa học và Công nghệ Việt Nam dưới sự hướng dẫn của GS. Vũ Cao Minh và bảo vệ năm 2023, chuyên ngành Địa chất học (Mã số: 9 44 02 01).

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

Báo cáo của Ngân hàng Thế giới năm 2018 cảnh báo các hiện tượng thiên tai cực đoan tại Việt Nam có xu hướng gia tăng về tần suất và cường độ. Theo số liệu của đề án "Điều tra, đánh giá và phân vùng cảnh báo nguy cơ trượt lở đất đá các vùng miền núi Việt Nam", các khối trượt nhỏ và nông xuất hiện với số lượng nhiều nhưng tổng thể tích nhỏ hơn nhiều so với các khối trượt lớn (KTL). Mặc dù số lượng KTL không nhiều, sự xuất hiện của chúng gây ra sức tàn phá nghiêm trọng đối với khu dân cư, công trình hạ tầng giao thông, ruộng bậc thang và đường dây truyền tải điện tại khu vực Miền núi phía Bắc (MNPB). Các dẫn chứng thực tế điển hình gồm:

  • Mùa mưa 2007–2008, hiện tượng trượt khối lớn (TKL) xuất hiện tại thị trấn Cốc Pài (huyện Xín Mần, tỉnh Hà Giang), nơi trụ sở UBND huyện và các phòng ban nằm trên khối trượt có chiều dài 350–500 m, chiều rộng 150–200 m.
  • Cuối tháng 7 năm 2012, KTL tại chợ Tân Sơn (xã Nấm Dẩn, huyện Xín Mần) phá hủy một phần khu chợ và làm sập nhà của 5 hộ dân.
  • Nhiều KTL khác nằm rải rác tại các xã Bản Díu, Nấm Dẩn, Xín Mần, Chế Là, Quảng Nguyên, Nà Chì.

Khoảng trống nghiên cứu được xác định: Các nghiên cứu trước đây chủ yếu tập trung mô tả từng điểm trượt riêng lẻ hoặc đánh giá nguy cơ trượt lở chung mà chưa hệ thống hóa các đặc trưng hình thái, cấu trúc của KTL trên quy mô toàn vùng MNPB. Đồng thời, còn thiếu việc ứng dụng kết hợp các mô hình học máy với hệ thông tin địa lý (GIS) để đánh giá nguy cơ KTL diện rộng và thiếu các nghiên cứu quan trắc địa kỹ thuật dài hạn nhằm định lượng ngưỡng kích hoạt động lực học của KTL phục vụ cảnh báo sớm.

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

  1. Làm sáng tỏ đặ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.
  2. Làm sáng tỏ đặc điểm động lực học của quá trình phát triển KTL bằng các công nghệ quan trắc hiện trường tại huyện Xín Mần, tỉnh Hà Giang.

Luận điểm bảo vệ

  • Luận điểm 1: Khối trượt lớn ở khu vực MNPB được đặc trưng bởi thể tích lớn hơn 4.500 m³, 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.
  • Luận điểm 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.
  • Luận điểm 3: Dịch chuyển khối trượt lớn 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 ($r_u$) ở vị trí mặt trượt $\ge 0,53$.

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

  • Đối tượng nghiên cứu: Một số KTL ở khu vực MNPB và nghiên cứu chi tiết động lực học các KTL tại huyện Xín Mần, tỉnh Hà Giang. Loại trừ: dòng lũ bùn đá độc lập, các khối trượt hình thành do khai thác khoáng sản, xây dựng và vận hành công trình thủy điện – thủy lợi.
  • Phạm vi không gian: Toàn bộ khu vực MNPB gồm 15 tỉnh (Hà Giang, Cao Bằng, Lào Cai, Bắc Kạn, Lạng Sơn, Tuyên Quang, Yên Bái, Thái Nguyên, Phú Thọ, Bắc Giang, Lai Châu, Điện Biên, Sơn La, Hòa Bình, Quảng Ninh) với tổng diện tích 95.000 km²; nghiên cứu chi tiết tại huyện Xín Mần (17 xã, diện tích 582 km²).
  • Phạm vi thời gian: Số liệu điều tra, khảo sát, ảnh viễn thám và tài liệu thu thập từ năm 2002 đến năm 2021; số liệu quan trắc biến dạng và áp lực nước lỗ rỗng tại hiện trường được duy trì liên tục từ năm 2017 đến năm 2021.

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

Các hướng nghiên cứu trên thế giới

  1. Khái niệm và phân loại: David Varnes (1958) và Cruden & Varnes (1996) xác lập định nghĩa nền tảng về trượt lở đất đá. Agliardi (2012) mô tả KTL có thể tích lớn, mặt trượt sâu, vận tốc chậm kéo dài. Lo (2017) nghiên cứu tại Đài Loan xác định KTL có mặt trượt trong đá gốc, diện tích > 10 ha, sâu > 5 m. Roering và cộng sự (2015) xác định diện tích KTL > 0,1 km². Lômtadze (1979) định nghĩa khối trượt gắn với sự phá hủy sườn dốc; Postoev (2010) phân chia thành 3 loại cơ chế trượt (trượt sâu do nén ép, trượt nông, trượt dạng dòng). Hiệp hội Quốc tế về Trượt đất đá (ICL) và Fell (1994) phân cấp thể tích trượt lở thành 7 cấp.
  2. Hình thái và kích thước: Jaboyedoff và cộng sự (2020) đề xuất công thức ước lượng thể tích khối trượt $V = \frac{1}{6}\pi \times D \times W \times L$. Heim (1932) đặt nền móng nghiên cứu góc Fahrböschung tương quan với thể tích trượt.
  3. Cơ chế và giai đoạn phát triển: Lômtadze (1979) phân chia 3 giai đoạn (chuẩn bị trượt, thành tạo trượt thực thụ, ổn định lại); Leroueil (2001) làm rõ 4 giai đoạn (tiền mất ổn định, bắt đầu mất ổn định, sau mất ổn định, kích hoạt lại).
  4. Phương pháp khảo sát và mô phỏng: Barry Voight (1978) khảo sát KTL Bắc Mỹ; Strozzi và cộng sự (2013) dùng 205 bộ ảnh bay chụp qua 50 năm tại Osco (Thụy Sỹ); Chen (2016) dùng UAV nghiên cứu quá trình mở rộng KTL Tsaoling (Đài Loan) giai đoạn 1931–2002. Phân tích mô hình số ổn định mái dốc và động lực học phát triển qua các công trình của Crosta và cộng sự (2002), Doug Stead (2013), Oldrich Hungr (1995), Moriwaki và cộng sự (2004) với mô hình dàn mưa nhân tạo, phần mềm LS-Rapid và hệ thống Anura3D.
  5. Quan trắc và cảnh báo sớm: 21 hệ thống quan trắc KTL tại châu Âu được tổng hợp (ghi đo lượng mưa, áp lực nước lỗ rỗng, dịch chuyển ngang sâu, mực nước ngầm); hệ thống cảm biến không dây truyền qua GPRS tại núi Longjingwan (Quý Châu, Trung Quốc); công nghệ radar giao thoa mặt đất (T-InSAR) theo dõi KTL Subiaco và Calatabiano tại Italia.

Các hướng nghiên cứu trong nước

  • Giai đoạn trước 1990: Các khảo sát sơ khai của Nguyễn Thanh, Hồ Chất phục vụ tái định cư và khai thác mỏ.
  • Nghiên cứu hiện trạng và đới đứt gãy: Các công trình nghiên cứu nứt trượt tại đồi Khau Cả, đồi Khí Tượng (Sơn La), đồi Ông Tượng (Hòa Bình). Đề án 15 tỉnh MNPB do TS. Lê Quốc Hùng chủ trì xây dựng cơ sở dữ liệu (CSDL) trượt lở; đề tài VT/UD-03/13-15 ứng dụng công nghệ viễn thám Tây Bắc; đề tài ĐTĐL.CN-23/17 tại Quảng Nam khảo sát trượt lở gây sóng trên sông Trường.
  • Đánh giá nguy cơ và phân vùng: Viện Địa chất (Nguyễn Quốc Thành, 2015) lập hệ thống bản đồ cảnh báo trượt lở; ứng dụng phương pháp AHP; Bùi Tiến Diệu và cộng sự (2016) so sánh mô hình chỉ số thống kê và hồi quy logistic tại Hòa Bình, tiếp nối bởi các nghiên cứu máy vector hỗ trợ (SVM), mạng neural nhân tạo (ANN) tại lưu vực thủy điện Sơn La và Quốc lộ 32.
  • Phân tích ổn định và quan trắc: Ứng dụng mô hình cân bằng giới hạn và phần tử hữu hạn (GeoStudio, Plaxis); Đỗ Minh Đức và cộng sự (2013, 2018) nghiên cứu ảnh hưởng của mưa và ngưỡng mưa gây trượt; Đoàn Huy Lợi và cộng sự (2017) mô phỏng trượt lở dạng dòng. Quan trắc hiện trường được triển khai tại Đồi Ông Tượng (2005), trạm Cốc Pài (đề tài KC.33/06-10), chương trình hợp tác Na Uy VINOGEO-SRV 07/056 tại Xín Mần và Quang Bình.

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

Luận án giải quyết khoảng trống về việc thiếu tiêu chí phân loại định lượng ranh giới thể tích cho KTL tại MNPB, chưa có sự kết hợp đối sánh giữa các thuật toán học máy nâng cao trong thành lập bản đồ nhạy cảm KTL toàn vùng, và thiếu các phân tích định lượng chuỗi số liệu quan trắc sâu đa chỉ tiêu dài hạn (Inclinometer, Piezometer, T-InSAR) để xác lập giá trị giới hạn áp lực nước lỗ rỗng kích hoạt chuyển động trượt.

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

Khung lý thuyết và cách tiếp cận

  1. Tiếp cận phân tích hệ thống: Dựa trên cơ sở lý thuyết hệ thống của Hoàng Tụy (1987), quá trình hình thành và phát triển KTL được mô tả bằng toán tử chuyển trạng thái: $$x(t_s) = G(x(t_0), u(.), t_0, t_s)$$ Trong đó, tập tác động $U$ gồm 3 nhóm:
    • $u_1(x)$: Biến đổi rất chậm theo không gian (thạch học, kiến tạo, vỏ phong hóa, địa hình).
    • $u_2(x, t)$: Biến đổi chậm theo không gian và thời gian (thủy văn, khí hậu, sử dụng đất, nhân sinh).
    • $u_3(x, t)$: Biến đổi nhanh và rất nhanh (mưa, động đất).
  2. Tiếp cận tích hợp chuỗi thời gian: Tích hợp các sự kiện trượt độc lập, quá trình mất ổn định ngắn hạn và diễn biến dịch chuyển biến dạng chậm dài hạn trong cùng một thang quy chiếu động lực học.

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

  • Khảo sát thực địa và địa chất công trình: Khảo sát chi tiết 26 KTL bằng phiếu điều tra chuẩn hóa; khoan khảo sát tại 02 KTL (Tân Sơn và Díu Thượng) bằng máy khoan XJ-1A, nòng đôi T6 không dùng dung dịch, xếp mẫu vào khay định vị 5 m; phân loại vỏ phong hóa (VPH) theo 6 đới của Garet Hearn (1995).
  • Khảo sát UAV và mô hình số bề mặt: Bay chụp bằng thiết bị bay không người lái (UAV), nắn ảnh và xây dựng mô hình số bề mặt (DSM), mô hình số độ cao (DEM) bằng phần mềm Agisoft Metashape theo phương pháp tam giác không gian ảnh và bình sai đa tia.
  • Thống kê và học máy: Phân tích phân phối thể tích (biểu đồ Cullen & Frey); phân tích tương quan cặp đôi, phân tích thành phần chính (PCA); kiểm tra đa cộng tuyến qua hệ số phóng đại phương sai (VIF); lựa chọn biến bằng phương pháp trung bình hóa mô hình Bayes (BMA); xây dựng bản đồ nguy cơ bằng Hồi quy Logistic, SVM, ANN và mô hình tổng hợp xếp chồng (Stacking/Ensemble); đánh giá mô hình bằng đường cong đặc trưng máy thu (ROC), diện tích dưới đường cong (AUC) và hệ số Cohen's Kappa trên tập dữ liệu chia ngẫu nhiên (80% huấn luyện, 20% kiểm định).
  • Quan trắc địa kỹ thuật hiện trường: Đo dịch chuyển ngang sâu trong hố khoan bằng cảm biến Inclinometer (InCL); đo áp lực nước lỗ rỗng bằng cảm biến Piezometer (PIE); đo mưa tự động; đo dịch chuyển bề mặt bằng radar giao thoa mặt đất (T-InSAR).
  • Phần mềm mô hình số địa kỹ thuật: Mô hình Seep/W (mô phỏng trường thấm), Slope/W (tính hệ số ổn định FoS theo phương pháp cân bằng giới hạn Morgenstern-Price) trong bộ phần mềm GeoStudio; mô phỏng phạm vi trượt động học và vận tốc chuyển động bằng phần mềm LS-Rapid kết hợp số liệu thí nghiệm cắt xoay vòng trong phòng thí nghiệm.

Nguồn tư liệu và dữ liệu

Nhóm dữ liệu Nguồn gốc và quy cách kỹ thuật Mục đích sử dụng
CSDL hiện trạng trượt lở Đề án 15 tỉnh MNPB (Viện KHĐC&KS), Báo cáo Tổng cục PCTT (2015–2021), đề tài VINOGEO-SRV, ĐTĐL Lập tập mẫu điểm trượt huấn luyện và kiểm định
Bản đồ địa chất 19 tờ bản đồ địa chất và khoáng sản MNPB tỷ lệ 1:200.000 (xuất bản 2005) Chiết xuất yếu tố thạch học, đứt gãy, ranh giới địa chất
Dữ liệu địa hình (DEM) 51 tấm ảnh vệ tinh ALOS PALSAR độ phân giải 12,5 $\times$ 12,5 m (NASA/JAXA, 2010) Trích xuất độ cao, độ dốc, hướng dốc, độ cong, TWI, độ sâu thung lũng
Ảnh viễn thám thực vật Ảnh vệ tinh Landsat 8 (USGS, chụp tháng 12/2019) Tính toán chỉ số khác biệt thực vật NDVI
Dữ liệu khí tượng Bản đồ đẳng mưa Atlas khí hậu; số liệu mưa giờ trạm KTTV Việt Bắc (từ 2010) Xác định lượng mưa năm và lượng mưa trận
Số liệu khảo sát & quan trắc 26 KTL đo vẽ trực tiếp; chuỗi đo Inclinometer, Piezometer, T-InSAR tại XM-HG (2017–2021) Phân tích cơ chế và kiểm chứng mô hình động lực học

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

Chương I: Tổng quan các nghiên cứu về khối trượt lớn trên thế giới và ở Việt Nam

Chương này tổng thuật toàn diện các khái niệm, khung phân loại và phương pháp tiếp cận nghiên cứu trượt lở đất đá. Về phân loại thể tích, tác giả tổng hợp thang phân cấp của ICL, Fell (1994) và Lômtadze, chỉ rõ sự khác biệt giữa các ngưỡng thể tích từ nhỏ ($<100\text{ m}^3$) đến đặc biệt lớn ($>10^8\text{ m}^3$). Về hình thái, chương làm rõ các bộ phận cấu trúc: mặt trượt, đỉnh, vách trượt chính/thứ cấp, thân trượt, chân trượt, phạm vi lắng đọng và góc Fahrböschung.

Chương tổng kết các yếu tố hình thành (nội sinh, thạch học, cấu trúc đứt gãy) và yếu tố kích hoạt (mưa, động đất, áp lực nước lỗ rỗng, hoạt động nhân sinh như cắt xẻ taluy chiếm 85% vị trí trượt ghi nhận). Chương phân tích 4 giai đoạn phát triển của Leroueil (2001) và điểm lại các công nghệ quan trắc hiện đại (Inclinometer, Piezometer, T-InSAR), các giải pháp giảm nhẹ công trình (thoát nước mặt, giếng khoan hạ áp) và phi công trình (công viên chứng tích thiên tai tại Đài Loan, Úc, Nhật Bản).

Chương II: Cách tiếp cận và hệ phương pháp nghiên cứu

Chương II thiết lập cơ sở phương pháp luận gồm tiếp cận phân tích hệ thống (xây dựng toán tử chuyển trạng thái với các véc-tơ biến đổi theo thời gian $u_1, u_2, u_3$) và tiếp cận tích hợp chuỗi thời gian. Khung logic nghiên cứu 8 bước được cụ thể hóa thành quy trình khép kín:

  1. Tổng hợp dữ liệu hiện trạng trượt lở và các yếu tố ảnh hưởng.
  2. Phân tích tương quan, đa cộng tuyến (VIF) và thành phần chính (PCA).
  3. Xây dựng mô hình nhạy cảm KTL tối ưu bằng thống kê và học máy.
  4. Tổng hợp sự kiện và giới hạn phạm vi ảnh hưởng.
  5. Phân tích thống kê chuỗi thời gian dữ liệu quan trắc dài hạn.
  6. Dự đoán động lực học bằng mô hình số (Seep/W, Slope/W, LS-Rapid).
  7. Phân tích chi tiết đặc điểm hình thành từng KTL điển hình.
  8. Đề xuất giải pháp giảm thiểu thiệt hại và quy trình cảnh báo sớm.

Chương trình bày chi tiết kỹ thuật bay chụp UAV, quy trình xử lý bình sai ảnh trong Agisoft Metashape để thành lập DEM/DSM, kỹ thuật khoan máy XJ-1A và tiêu chuẩn mô tả 6 đới vỏ phong hóa theo Garet Hearn.

Chương III: Đặc điểm các khối trượt lớn ở miền núi phía Bắc

Chương III tập trung phân tích đặc trưng định lượng của các KTL tại 15 tỉnh MNPB:

  • Phân phối thống kê thể tích: Kiểm định phân phối logarit thể tích ($\text{Vlog}$) trên biểu đồ Cullen & Frey cho thấy quy luật phân phối xác suất rõ rệt. Kết quả xác lập ngưỡng thể tích để phân loại một khối trượt là KTL tại MNPB là $V > 4.500\text{ m}^3$.
  • Đặc điểm hình thái và hình học: Kết quả đo vẽ chi tiết 26 KTL điển hình xác định mối tương quan chặt chẽ giữa chiều sâu mặt trượt ($D_r$), chiều dài ($L_r$), chiều rộng ($W_r$) và chiều cao ($H_r$).
  • Cấu trúc mặt trượt và VPH: Mặt trượt của các KTL tại MNPB không tạo vòm cong đồng nhất trong đất mà hầu hết trùng với các mặt bất liên tục thạch học, mặt phân lớp hoặc đới khe nứt kiến tạo (được minh chứng qua các mặt cắt VPH tại Thác Khanh - Hòa Bình, Bắc Yên - Sơn La, Nậm Búng - Yên Bái, Móng Sến - Lào Cai, Cốc Pài - Hà Giang, Chế Cu Nha - Yên Bái).
  • Kiểu loại và diễn biến: Phân loại cho thấy KTL dạng tịnh tiến và dạng hỗn hợp chiếm ưu thế tuyệt đối. Diễn biến dịch chuyển mang tính trượt chậm, mở rộng qua nhiều chu kỳ mưa lũ (như ghi nhận biến dạng KTL Chế Cu Nha qua các năm 2010, 2016, 2020).

Chương IV: Đánh giá nguy cơ hình thành khối trượt lớn ở miền núi phía Bắc

Chương IV thực hiện đánh giá không gian nguy cơ KTL trên toàn vùng MNPB:

  • Phân tích các yếu tố ảnh hưởng: Xem xét tổ hợp 12 biến đầu vào gồm thành phần thạch học, khoảng cách đến ranh giới thạch học, khoảng cách đến đứt gãy kiến tạo, cao độ, hướng dốc, độ dốc, độ cong ngang sườn, độ cong dọc sườn, độ sâu thung lũng, chỉ số ẩm địa hình (TWI), khoảng cách đến sông suối, lượng mưa năm/lượng mưa trận, chỉ số NDVI và khoảng cách đến đường giao thông.
  • Kiểm định thống kê: Phân tích tương quan cặp đôi, kiểm định VIF loại trừ hiện tượng đa cộng tuyến và dùng BMA chọn lọc mô hình hồi quy logistic tối ưu.
  • Đánh giá mô hình học máy: So sánh hiệu năng giữa Hồi quy Logistic, Máy vector hỗ trợ (SVM), Mạng neural nhân tạo (ANN) và mô hình xếp chồng tổng hợp (Stacking). Kết quả kiểm định qua diện tích dưới đường cong AUC trên tập dữ liệu kiểm tra độc lập (20%) khẳng định các mô hình học máy đạt độ tin cậy cao trong dự báo vùng nhạy cảm KTL.
  • Quy luật phân bố: KTL tập trung mật độ cao nhất tại các đới cấu trúc dập vỡ kiến tạo mạnh, sườn dốc $> 25^\circ$, vùng chuyển tiếp địa tầng thạch học có tính bất đồng nhất cao và các sườn dốc chịu tác động đào bạt taluy làm đường.

Chương V: Phân tích đặc điểm động lực học khối trượt lớn khu vực huyện Xín Mần, tỉnh Hà Giang

Chương V đi sâu phân tích động lực học thực địa tại hai trạm quan trắc Tân Sơn (xã Nấm Dẩn) và Díu Thượng (xã Bản Díu):

  • Phân tích chuỗi quan trắc địa kỹ thuật (2017–2021): Tích hợp số liệu mưa giờ, mực nước ngầm, áp lực nước lỗ rỗng (PWP) từ Piezometer và dịch chuyển ngang từ Inclinometer. Kết quả chỉ ra quá trình trượt diễn ra gián đoạn. Khi hệ số áp lực nước lỗ rỗng đạt giá trị tới hạn $r_u \ge 0,53$, vận tốc và biến dạng dịch chuyển ngang tăng vọt rõ rệt.
  • Kiểm chứng bằng T-InSAR: Số liệu đo biến dạng mặt đất bằng radar giao thoa (điểm P14) tương thích chặt chẽ với số liệu đo dịch chuyển sâu bằng Inclinometer trong giai đoạn từ tháng 8/2019 đến tháng 12/2019.
  • Mô hình số địa kỹ thuật: Kiểm toán hệ số ổn định mái dốc (Slope/W) tại các mức $r_u$ khác nhau xác nhận mái dốc chuyển sang trạng thái mất ổn định ($FoS < 1,0$) khi $r_u \ge 0,53$. Mô phỏng động lực học bằng phần mềm LS-Rapid (sử dụng thông số thí nghiệm cắt xoay) xác định chính xác phạm vi ảnh hưởng động học ($L_f$), vận tốc và bề dày tích tụ của KTL Tân Sơn và Díu Thượng.
  • Đề xuất giải pháp: Nhận diện các dấu hiệu mất ổn định ngoại vi; định hướng quy hoạch sử dụng đất; xây dựng quy trình 4 bước phân tích dữ liệu quan trắc phục vụ cảnh báo sớm KTL.

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

Đóng góp mới về mặt khoa học và lý luận

  1. Hệ thống hóa toàn diện các đặc trưng hình học, hình thái và cơ chế trượt của KTL tại khu vực MNPB; xác lập định lượng ngưỡng thể tích $V > 4.500\text{ m}^3$ làm tiêu chí nhận diện KTL cho toàn khu vực.
  2. Hoàn thiện phương pháp luận nghiên cứu KTL thông qua việc kết hợp phân tích hệ thống (toán tử chuyển trạng thái) với tiếp cận tích hợp chuỗi thời gian đa nguồn số liệu.
  3. Làm sáng tỏ cơ chế kích hoạt động lực học của KTL tại huyện Xín Mần: Dịch chuyển KTL là quá trình không liên tục, trong đó hệ số áp lực nước lỗ rỗng tại mặt trượt $r_u \ge 0,53$ là ngưỡng động lực quyết định sự gia tăng đột biến của vận tốc dịch chuyển và gây mất ổn định khối trượt.

Đóng góp mới về mặt thực tiễn và ứng dụng

  1. Xây dựng bản đồ phân vùng nguy cơ hình thành KTL cho 15 tỉnh MNPB bằng việc tích hợp các mô hình học máy (Logistic, SVM, ANN, Stacking) trên nền tảng GIS, cung cấp công cụ trực quan phục vụ quy hoạch lãnh thổ.
  2. Thiết lập và duy trì thành công mạng lưới quan trắc hiện trường đa thông số (Inclinometer, Piezometer, đo mưa, T-InSAR) tại hai KTL trọng điểm Tân Sơn và Díu Thượng (Xín Mần, Hà Giang).
  3. Đề xuất quy trình phân tích dữ liệu quan trắc theo thời gian thực và mô hình số hóa phạm vi ảnh hưởng (LS-Rapid), phục vụ trực tiếp cho công tác cảnh báo sớm và đề xuất các giải pháp công trình, phi công trình nhằm giảm nhẹ thiệt hại do TKL tại miền núi phía Bắc.

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

Hạn chế của luận án

  • Dữ liệu bản đồ địa chất nền sử dụng cho đánh giá diện rộng toàn vùng MNPB là tỷ lệ 1:200.000 do hệ thống bản đồ địa chất tỷ lệ 1:50.000 chưa phủ kín toàn bộ 15 tỉnh.
  • Chưa phân tích định lượng chi tiết ảnh hưởng của yếu tố kích hoạt do động đất đến sự hình thành KTL tại MNPB do thiếu chuỗi số liệu quan trắc địa chấn vi mô đồng bộ tại các điểm trượt.
  • Số lượng trạm quan trắc địa kỹ thuật chuyên sâu mới được đầu tư triển khai tại 02 vị trí đại diện thuộc huyện Xín Mần, tỉnh Hà Giang.

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

  • Cập nhật và nâng cấp bản đồ nhạy cảm KTL lên tỷ lệ chi tiết (1:50.000 và 1:25.000) khi CSDL địa chất công trình và địa chất khoáng sản được hoàn thiện.
  • Mở rộng lắp đặt mạng lưới quan trắc địa kỹ thuật tự động tại các KTL có nguy cơ đe dọa trực tiếp các đô thị, trung tâm huyện lỵ và các hồ chứa thủy điện lớn ở MNPB.
  • Tích hợp thêm các thuật toán học sâu và chuỗi ảnh radar viễn thám vệ tinh (InSAR đa thời gian) để theo dõi biến dạng diện rộng liên tục.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên sâu phục vụ các đối tượng:

  • Nghiên cứu sinh, học viên cao học và giảng viên: Thuộc các chuyên ngành Địa chất học, Địa kỹ thuật, Địa chất công trình, Địa lý tài nguyên và môi trường, Kỹ thuật xây dựng công trình ngầm và mỏ. Phần phương pháp luận tiếp cận hệ thống, xử lý đa cộng tuyến (VIF, BMA) và ứng dụng học máy trong phân vùng nguy cơ trượt lở (Chương II, Chương IV) có giá trị học thuật cao.
  • Cơ quan quản lý và quy hoạch: Ban Chỉ huy Phòng chống thiên tai và Tìm kiếm cứu nạn các cấp, Sở Tài nguyên và Môi trường, Sở Nông nghiệp và Phát triển nông thôn, Sở Giao thông Vận tải 15 tỉnh MNPB. Phần bản đồ phân vùng nguy cơ (Chương IV) và các giải pháp quy hoạch, nhận diện dấu hiệu mất ổn định (Chương V) là căn cứ thực tiễn để định hướng phân bổ dân cư và xây dựng công trình giao thông an toàn.
  • Kỹ sư địa kỹ thuật và đơn vị tư vấn thiết kế: Tham khảo trực tiếp quy trình khoan nòng ba lấy mẫu VPH, quy trình bố trí lắp đặt cảm biến Piezometer, Inclinometer trong hố khoan và phương pháp mô phỏng ổn định bằng Seep/W, Slope/W kết hợp LS-Rapid (Chương II, Chương V).

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

1. Luận án xác định tiêu chí thể tích để phân loại một khối trượt là "khối trượt lớn" (KTL) ở miền núi phía Bắc là bao nhiêu?

Dựa trên phân tích phân phối thống kê logarit thể tích ($\text{Vlog}$) của CSDL hiện trạng trượt lở tại MNPB và kiểm định biểu đồ Cullen & Frey, luận án xác định các khối trượt có thể tích $V > 4.500\text{ m}^3$ được xếp vào nhóm khối trượt lớn (KTL).

2. Những yếu tố tự nhiên và nhân sinh nào đóng vai trò chính quyết định sự hình thành KTL ở MNPB?

Kết quả phân tích thống kê đa biến, thành phần chính (PCA) và mô hình học máy chỉ ra 4 yếu tố chi phối chính:

  • Kiến tạo (mật độ đứt gãy và đới dập vỡ kiến tạo);
  • Thạch học (tính bất đồng nhất và mặt tiếp xúc địa tầng);
  • Độ dốc địa hình (sườn dốc lớn);
  • Hoạt động nhân sinh (đào bạt taluy làm đường giao thông, thay đổi độ dốc sườn tự nhiên).

3. Ngưỡng áp lực nước lỗ rỗng kích hoạt dịch chuyển lớn tại KTL huyện Xín Mần được xác định là bao nhiêu?

Qua phân tích chuỗi dữ liệu quan trắc hiện trường từ cảm biến Piezometer và Inclinometer tại trạm Tân Sơn và Díu Thượng (2017–2021), luận án xác định khi hệ số áp lực nước lỗ rỗng tại mặt trượt đạt giá trị $r_u \ge 0,53$, vận tốc dịch chuyển ngang và độ dịch chuyển tích lũy của KTL tăng lên đột biến, dẫn đến trạng thái mất ổn định ($FoS < 1,0$).

4. Hệ thống quan trắc hiện trường tại huyện Xín Mần sử dụng các thiết bị và công nghệ nào?

Hệ thống quan trắc tại hai KTL Tân Sơn (xã Nấm Dẩn) và Díu Thượng (xã Bản Díu) bao gồm:

  • Cảm biến Piezometer (PIE) đặt trong hố khoan đo áp lực nước lỗ rỗng;
  • Cảm biến Inclinometer (InCL) đặt trong hố khoan đo dịch chuyển ngang sâu;
  • Thiết bị đo mưa tự động ghi nhận lượng mưa giờ;
  • Công nghệ radar giao thoa mặt đất khẩu độ tổng hợp (T-InSAR) theo dõi vi dịch chuyển bề mặt.

5. Luận án sử dụng các mô hình số nào để kiểm toán ổn định và dự báo phạm vi ảnh hưởng của KTL?

Luận án ứng dụng tổ hợp phần mềm địa kỹ thuật chuyên dụng:

  • Module Seep/W để mô phỏng trường thấm và áp lực nước lỗ rỗng trong thân trượt;
  • Module Slope/W để tính toán hệ số ổn định mái dốc ($FoS$) theo phương pháp cân bằng giới hạn Morgenstern-Price;
  • Phần mềm LS-Rapid kết hợp thông số thí nghiệm cắt xoay vòng để mô phỏng động lực học chuyển động, tính toán vận tốc, bề dày lắng đọng và phạm vi ảnh hưởng ($L_f$) của khối trượt.

Kết luận

Luận án của tác giả Đào Minh Đức là công trình nghiên cứu toàn diện về đặc điểm hình thành, quy luật phân bố và động lực học phát triển của các khối trượt lớn tại khu vực miền núi phía Bắc, lấy ví dụ thực nghiệm tại huyện Xín Mần, tỉnh Hà Giang. Công trình đã kết hợp thành công phương pháp phân tích hệ thống và chuỗi thời gian, ứng dụng các mô hình học máy tiên tiến trong đánh giá nguy cơ diện rộng và thiết lập mạng lưới quan trắc địa kỹ thuật đa chỉ tiêu dài hạn tại hiện trường. Kết quả xác lập ngưỡng thể tích $V > 4.500\text{ m}^3$ cho KTL vùng MNPB và phát hiện ngưỡng hệ số áp lực nước lỗ rỗng $r_u \ge 0,53$ kích hoạt dịch chuyển trượt là những đóng góp khoa học có giá trị lý luận và thực tiễn cao, phục vụ hiệu quả cho công tác cảnh báo sớm và giảm nhẹ rủi ro thiên tai trượt lở đất đá.