Tổng quan về luận án
Nghiên cứu địa động lực hiện đại và đánh giá rủi ro động đất tại đới đứt gãy Tây Bắc Việt Nam đại diện cho một trong những bài toán phức tạp và then chốt nhất của khoa học Trái đất khu vực Đông Nam Á. Luận án tiến sĩ địa chất chuyên ngành Địa kiến tạo (mã số: 62 44 55 05) của nghiên cứu sinh Nguyễn Hữu Tuyên, thực hiện dưới sự hướng dẫn khoa học của PGS. Chu Văn Ngợi và PGS. Cao Đình Triều tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội (2012) với tiêu đề "Nghiên cứu địa động lực vùng Tuần Giáo và kế cận, xác lập cơ sở khoa học đánh giá và dự báo động đất", đã thiết lập một bước đột phá học thuật mang tính nền tảng. Luận án giải quyết mắt xích xung yếu về mối quan hệ nhân quả giữa bối cảnh địa động lực nội mảng (intra-plate geodynamics) với hoạt động phát sinh động đất mạnh tại khu vực từng xảy ra trận động đất lịch sử Tuần Giáo năm 1983 ($M_s = 6,8$).
Khoảng trống nghiên cứu (research gap) được luận án nhận diện xuất phát từ thực tế các phương pháp dự báo động đất truyền thống tại Việt Nam (Nguyễn Đình Xuyến, 1996; Cao Đình Triều, 2003) chủ yếu dựa trên mô hình xác suất - thống kê địa chấn đơn thuần. Cách tiếp cận thống kê thuần túy bộc lộ hai hạn chế cốt tử: chuỗi số liệu quan trắc địa chấn lịch sử (1900–2009) quá ngắn ngủi so với chu kỳ lặp lại kéo dài hàng nghìn năm của các đứt gãy sinh chấn; đồng thời hoàn toàn tách rời các thông số cấu trúc sâu, dị thường địa vật lý và trường biến dạng ứng suất Coulomb của môi trường thạch quyển. Luận án đã lấp đầy khoảng trống này bằng cách tích hợp toàn diện các trường địa vật lý, cấu trúc sâu vỏ Trái đất, địa mạo kiến tạo định lượng và thuật toán nhận dạng mẫu (pattern recognition) đa biến.
Hệ thống câu hỏi nghiên cứu và giả thuyết khoa học được xác lập rõ ràng:
- Câu hỏi nghiên cứu 1 (RQ1): Cấu trúc vỏ Trái đất vùng Tuần Giáo và kế cận bị phân cắt thành các khối cấu trúc địa động lực (CTĐĐL) bậc mấy, và đặc trưng động học dị biệt (thẳng đứng và nằm ngang) giữa các khối diễn ra như thế nào trong Kainozoi muộn đến hiện đại?
- Câu hỏi nghiên cứu 2 (RQ2): Trạng thái tập trung ứng suất kiến tạo hiện đại theo mô hình phá hủy Coulomb tích lũy tại các đới đứt gãy ranh giới khối có quy luật phân bố không gian và độ sâu ra sao?
- Câu hỏi nghiên cứu 3 (RQ3): Làm thế nào để xác định chính xác các nút giao cấu trúc (structural nodes) có tiềm năng phát sinh động đất mạnh với các ngưỡng độ lớn khác nhau ($M_0 \ge 4,0$; $M_0 \ge 5,0$; và $M_s = 6,7 \div 6,8$) thông qua thuật toán nhận dạng phi thống kê?
Các giả thuyết nghiên cứu tương ứng:
- Giả thuyết H1: Vỏ Trái đất khu vực Tây Bắc không phải là một thể đồng nhất mà phân dị thành 06 khối CTĐĐL cấp II độc lập, bị khống chế bởi các đới đứt gãy sâu có vận động trượt bằng kết hợp nâng hạ phân dị.
- Giả thuyết H2: Các tâm điểm tích tụ ứng suất Coulomb cực đại (> 4 bar) trùng khớp với các đới đứt gãy hoạt động nghịch/trượt bằng và các nút giao kiến tạo bậc cao.
- Giả thuyết H3: Thuật toán nhận dạng mẫu CORA3 và phần mềm COSCAD có khả năng phân lớp chính xác các nút giao cấu trúc nguy hiểm thông qua tổ hợp 18 tham số địa chất - địa vật lý - địa mạo.
Khung lý thuyết của công trình tích hợp học thuyết Kiến tạo mảng (Plate Tectonics), mô hình kiến tạo đào thoát (Escape Tectonics / Extrusion Model of Tapponnier), lý thuyết Trường biến dạng ứng suất Coulomb (Coulomb Failure Stress), và lý thuyết Địa chấn kiến tạo cấu trúc (Seismotectonics of Gubbin & Gelfand). Phạm vi nghiên cứu bao quát diện tích đất liền Tây Bắc Việt Nam giới hạn trong tọa độ địa lý: $\phi = 21^\circ 00' \div 22^\circ 30'$ vĩ độ Bắc và $\lambda = 103^\circ 00' \div 104^\circ 30'$ kinh độ Đông. Quy mô dữ liệu bao gồm 164 trang luận án, 79 hình vẽ mô phỏng 2D/3D, 12 bảng biểu số liệu chi tiết, danh mục động đất giai đoạn 1900–2009, xử lý 70 nút giao cấu trúc độc lập, kết hợp phân tích 11 tuyến địa mạo và 08 tuyến đo trọng lực chi tiết tại các cụm thủy điện bậc thang sông Đà.
Literature Review và Positioning
Quá trình tiến hóa địa động lực Kainozoi khu vực Đông Nam Á chịu sự chi phối quyết định bởi tương tác giữa ba mảng thạch quyển khổng lồ: Âu - Á, Thái Bình Dương và Ấn - Úc. Dòng nghiên cứu quốc tế đã ghi nhận sự tranh biện sâu sắc giữa hai mô hình địa động lực kinh điển. Trường phái thứ nhất, dẫn đầu bởi Tapponnier et al. (1982, 1986), Peltzer & Tapponnier (1988), Briais et al. (1993) và Leloup et al. (1993, 1995), khẳng định cơ chế xô húc đâm thụt của tiểu lục địa Ấn Độ vào lục địa Âu - Á từ khoảng 45–50 triệu năm trước đã gây ra hiện tượng trồi trượt đẩy thoát (escape tectonics) của địa khối Đông Dương (Indochina) về phía đông nam dọc theo đới đứt gãy trượt bằng trái quy mô lớn Sông Hồng với biên độ dịch chuyển tích lũy từ 330 $\pm$ 60 km đến 500 km, tạo điều kiện mở bồn trũng Biển Đông dạng kéo toạc (pull-apart basin) trong khoảng 32 đến 16 triệu năm trước.
Ngược lại, trường phái thứ hai (Taylor & Hayes, 1980, 1983; Hall, 2002; Holloway, 1982) lại nhấn mạnh vai trò của sự hút chìm mảng Proto-Biển Đông cổ xuống dưới rìa tây bắc Borneo (đới khâu Lupar - Kuching) là động lực kéo rách và tách giãn đáy đại dương, xem chuyển động trượt bằng trên đứt gãy Sông Hồng chỉ đóng vai trò thứ yếu. Thêm vào đó, các công trình của Dewey et al. (1989) và Harrison et al. (1992, 1996) đã làm rõ sự chuyển đổi pha động học quan trọng: pha 1 (Eocene - Oligocene sớm, 45–32 Ma) nén ép dày vỏ Tây Tạng lên tới 70 km; pha 2 (Oligocene sớm - Miocene giữa, 32–16 Ma) trượt bằng trái mãnh liệt; và pha 3 (Miocene muộn đến nay, khoảng 5–8 Ma trở lại đây) chuyển đổi cơ chế sang trượt bằng phải do tác động chi phối của va mảng Ấn - Úc hiện đại.
Tại Việt Nam, các công trình điều tra địa chất của Fromaget (1937–1952), Dovjikov (1965), Trần Văn Trị & Gatinxki (1970) đã phân chia các đới kiến trúc truyền thống theo quan điểm địa máng uốn nếp. Sang kỷ nguyên kiến tạo mảng, các nhà khoa học Việt Nam đã đạt nhiều bước tiến: Phạm Khoan (1971) và Quách Văn Gừng (1978) xác định độ dày vỏ Trái đất tăng từ ven biển (31–33 km) vào sâu trong đất liền Tây Bắc (45–47 km); Cao Đình Triều (2000, 2003) phân tích trường ứng suất thạch quyển ở độ sâu 50–100 km; Đặng Thanh Hải (2008) làm rõ độ sâu mặt Moho (24–32 km) và Conrad (10–20 km); Phan Trọng Trịnh (1996, 2004) xác lập trường ứng suất nén trượt Kainozoi; Trần Đình Tô et al. (2002, 2004) đo đạc mạng lưới trắc địa vệ tinh GPS; và Ngô Gia Thắng (2007) tính toán biến dạng thẳng đứng Pliocene - Đệ Tứ.
Luận án của Nguyễn Hữu Tuyên tự định vị tại điểm giao thoa giữa cấu trúc sâu thạch quyển, kiến tạo hiện đại và điều kiện vật lý môi trường sinh chấn. Khác biệt với các nghiên cứu tiền nhiệm vốn chỉ mô tả định tính hoặc phân tích thống kê địa chấn cô lập, công trình này tiên phong xác lập hệ tiêu chí phân khối CTĐĐL cấp II dựa trên bản chất vật lý (trường dị thường Bouguer, dị thường đẳng tĩnh, hệ số Poisson, độ sâu móng kết tinh, mặt Conrad và Moho), kết hợp định lượng chỉ số địa mạo kiến tạo ($V_f, S_{mf}$) và mô hình số độ cao DEM 30m.
Khi đối sánh với các nghiên cứu quốc tế tương đồng:
- So sánh với nghiên cứu hệ đứt gãy Tiên Thủy Hà (Xianshuihe) và Tiểu Giang (Xiaojiang) tại Vân Nam - Thanh Tạng (Wang et al., 2008; Shen et al., 2005): Các tác giả Trung Quốc sử dụng dày đặc dữ liệu GPS kết hợp mô phỏng đứt gãy khối trượt. Luận án của Nguyễn Hữu Tuyên đã phát triển cách tiếp cận tương đương trong bối cảnh chuỗi số liệu GPS nội mảng Việt Nam còn thưa, bằng cách bù đắp xuất sắc bằng hệ phương pháp nhận dạng nút giao kiến tạo CORA3 và phân tích mặt cắt cổ địa chấn $C^{14}$.
- So sánh với mô hình dự báo nút giao động đất của Viện Vật lý Trái đất Schmidt - Viện Hàn lâm Khoa học Nga (Gelfand et al., 1972, 1976; Soloviev et al., 2014): Luận án đã hiệu chỉnh hoàn hảo thuật toán CORA3 từ quy mô mảng toàn cầu sang điều kiện cấu trúc nội mảng phức tạp Tây Bắc Việt Nam, đưa vào 18 tham số đầu vào đặc thù thạch quyển nhiệt đới gió mùa có lớp phủ phong hóa dày.
Đóng góp lý thuyết và khung phân tích
Đóng góp cho lý thuyết
Luận án mở rộng đáng kể lý thuyết Bối cảnh kiến tạo (Tectonic Setting Theory) của Kent C. Condie (1997, 2005) và mô hình kiến tạo nội mảng hiện đại. Công trình chứng minh rằng trong một đới biến dạng nội mảng chịu lực đẩy trồi xa từ đới va chạm Hymalaya, sự tích lũy biến dạng không phân bố đồng đều mà bị phân cắt nghiêm ngặt bởi các ranh giới bất liên tục sâu, tạo thành các vi khối địa động lực có hành vi ứng xử cơ học hoàn toàn dị biệt.
Mô hình lý thuyết của luận án được cụ thể hóa qua 4 mệnh đề khoa học (Theoretical Propositions):
- Mệnh đề P1 (Phân khối vỏ Trái đất): Vỏ Trái đất khu vực Tây Bắc Việt Nam tồn tại dưới dạng một hệ thống khối cấu trúc địa động lực cấp II phân tầng và phân dị mạnh mẽ về mật độ, độ dày móng kết tinh ($1 \div 5$ km), độ sâu mặt Conrad ($10 \div 20$ km) và mặt Mohor ($24 \div 32$ km).
- Mệnh đề P2 (Cơ chế chuyển động phân dị): Ranh giới giữa các khối địa động lực là các đới đứt gãy hoạt động sâu, đóng vai trò hấp thụ và chuyển hóa trường biến dạng nén ép khu vực phương Bắc - Nam và Tây Bắc - Đông Nam thành các chuyển động trượt bằng (strike-slip) kết hợp trượt nghịch (thrusting) và biến dạng phân dị thẳng đứng sâu sắc.
- Mệnh đề P3 (Tập trung ứng suất Coulomb): Vùng có nguy cơ phá hủy dẻo - giòn và phát sinh động đất là nơi trường ứng suất Coulomb biến thiên cục bộ cực đại ($\Delta CFS > 0,5 \div 4,0$ bar), định vị chủ yếu dọc theo các đoạn uốn cong hình học hoặc khóa kẹt cơ học của đới đứt gãy.
- Mệnh đề P4 (Quy luật nút giao cấu trúc): Động đất mạnh không phát sinh ngẫu nhiên dọc theo chiều dài đứt gãy mà tập trung có quy luật tại các "nút giao cấu trúc" (structural nodes) bậc cao – nơi giao cắt của ít nhất 2 hệ thống đứt gãy cấp I-II và có gradient biến thiên trường địa mạo/địa vật lý đạt ngưỡng tới hạn.
Sự chuyển dịch giác độ nhận thức (paradigm shift) mà luận án mang lại là quá trình chuyển hóa từ việc xem động đất như một hiện tượng ngẫu nhiên thống kê sang một tất yếu vật lý cơ học đá (rock physics determinism) được điều khiển bởi cấu trúc sâu và trạng thái tương tác động lực giữa các khối thạch quyển.
Khung phân tích độc đáo
Khung phân tích của luận án tích hợp đa chiều 4 lý thuyết chuyên ngành: (1) Thuyết Kiến tạo mảng và vi mảng nội lục, (2) Cơ học biến dạng ứng suất phá hủy Coulomb, (3) Địa mạo kiến tạo định lượng (Tectonic Geomorphology), và (4) Lý thuyết Nhận dạng mẫu toán học (Pattern Recognition Algorithm).
+-----------------------------------------------------------------------------------+
| BỐI CẢNH ĐỊA ĐỘNG LỰC KHU VỰC |
| Va mảng Ấn - Úc và Âu - Á + Kiến tạo đào thoát Kainozoi |
+-----------------------------------------------------------------------------------+
+-----------------------------------------------------------------------------------+
| CẤU TRÚC SÂU VÀ TRƯỜNG ĐỊA VẬT LÝ |
| • Dị thường Trọng lực Bouguer / Đẳng tĩnh • Dị thường Từ Delta-Ta |
| • Độ sâu Móng kết tinh (1-5km) • Mặt Conrad (10-20km) • Mặt Mohor (24-32km) |
+-----------------------------------------------------------------------------------+
+-----------------------------------------------------------------------------------+
| PHÂN LẬP 06 KHỐI CẤU TRÚC ĐỊA ĐỘNG LỰC CẤP II |
| [1] Hoàng Liên Sơn [2] Sông Đà [3] Sơn La [4] Sông Mã [5] Sốp Cộp [6] ĐB|
+-----------------------------------------------------------------------------------+
+-----------------------------------------------------------------------------------+
| ĐẶC TRƯNG VẬN ĐỘNG VÀ BIẾN DẠNG HIỆN ĐẠI |
| • Chuyển động đứng (Thủy chuẩn lặp, N2-Q) • Chuyển động ngang (GPS ~24mm/năm) |
| • Địa mạo kiến tạo (Vf, Smf, Facets vách) • Trường biến dạng Coulomb (0.5-4 bar) |
| • Cổ địa chấn học (Trenching, đồng vị C14) • Phân tích khe nứt kiến tạo vật lý |
+-----------------------------------------------------------------------------------+
+-----------------------------------------------------------------------------------+
| XÁC LẬP VÀ MÃ HÓA 70 NÚT GIAO CẤU TRÚC (STRUCTURAL NODES) |
| Bộ 18 tham số địa chất - địa vật lý - địa mạo |
+-----------------------------------------------------------------------------------+
+-----------------------------------------------------------------------------------+
| THUẬT TOÁN NHẬN DẠNG CORA3 & PHẦN MỀM COSCAD |
| • Rời rạc hóa (COD) • Đánh giá trọng số (PRAL) • Nhận dạng phân lớp nguy cơ sinh chấn|
+-----------------------------------------------------------------------------------+
+-----------------------------------------------------------------------------------+
| BẢN ĐỒ DỰ BÁO NGUỒN PHÁT SINH ĐỘNG ĐẤT NGUY HIỂM |
| • Mo >= 4.0 • Mo >= 5.0 • Động đất cực đại Ms = 6.7 - 6.8 |
+-----------------------------------------------------------------------------------+
Điều kiện biên (boundary conditions) của khung phân tích được xác lập chặt chẽ: Áp dụng cho vùng vỏ lục địa chịu cơ chế ứng suất nén trượt sau va chạm (post-collision strike-slip/compressional regime), giới hạn trong tầng vỏ trên dẻo - giòn (seismogenic layer) ở độ sâu chấn tiêu phổ biến từ $5 \div 25$ km, với vận tốc hội tụ thạch quyển trung bình khu vực đạt xấp xỉ 10 cm/năm giữa các mảng lớn và chuyển dịch trượt dọc đới đứt gãy Sông Hồng xấp xỉ 24 mm/năm.
Phương pháp nghiên cứu tiên tiến
Thiết kế nghiên cứu
Luận án quán triệt lập trường triết học Thực chứng Phê phán (Critical Realism) kết hợp Chủ nghĩa Duy thực Khoa học (Scientific Realism). Triết lý này thừa nhận rằng các hiện tượng địa chấn quan sát được trên bề mặt chỉ là biểu hiện bên ngoài của các cơ chế cấu trúc sâu và trường ứng suất nội tại không thể quan sát trực tiếp, đòi hỏi phải được giải đoán thông qua mô hình hóa toán lý và đối chiếu thực nghiệm đa nguồn.
Thiết kế nghiên cứu đa tầng bậc (multi-level design) được triển khai trên 3 cấp độ:
- Cấp độ Vĩ mô (Macro-level): Đánh giá bối cảnh kiến tạo mảng Đông Nam Á và Tây Bắc Việt Nam trong mối tương quan với các đới hút chìm Sunda, Marian và đới va chạm Hymalaya.
- Cấp độ Trung mô (Meso-level): Xác lập 06 khối CTĐĐL cấp II và hệ thống đứt gãy ranh giới cấp I, II, III thông qua mô hình số độ cao DEM và các trường dị thường địa vật lý khu vực tỷ lệ 1:200.000.
- Cấp độ Vi mô (Micro-level): Đo vẽ chi tiết 70 nút giao cấu trúc, khảo sát các vết lộ địa chất, đo đạc hàng nghìn khe nứt kiến tạo, đo lặp thủy chuẩn, đo trắc địa vệ tinh GPS và đào hào cổ địa chấn học lấy mẫu phân tích tuổi tuyệt đối $C^{14}$.
Dữ liệu đầu vào chính xác bao gồm:
- Mô hình số độ cao DEM từ ảnh radar SRTM (Shuttle Radar Topography Mission) độ phân giải không gian 30m do NASA cung cấp.
- Bản đồ Địa chất Việt Nam tỷ lệ 1:200.000 (Tổng cục Địa chất và Khoáng sản, 2005).
- Bản đồ dị thường trọng lực Bouguer và dị thường từ hàng không ($\Delta Ta$) do Liên đoàn Vật lý Địa chất thành lập.
- Số liệu 08 tuyến đo trọng lực chi tiết tại các công trình thủy điện trọng điểm: Sơn La, Huội Quảng, Bản Chát.
- Bản đồ phân tích Tomography thạch quyển Việt Nam và Đông Nam Á.
- Dữ liệu đo lặp thủy chuẩn quốc gia qua 2 chu kỳ (1963 và 1985) do Cục Đo đạc và Bản đồ thực hiện.
- Dữ liệu đo dịch chuyển GPS giai đoạn 2002–2004 dọc các đới đứt gãy Lai Châu - Điện Biên và Sơn La - Sông Đà.
- Danh mục động đất chuẩn hóa khu vực Tây Bắc giai đoạn 1900–2009.
Quy trình nghiên cứu rigorous
Quy trình thu thập và xử lý dữ liệu được thiết kế kiểm soát sai số nghiêm ngặt thông qua kỹ thuật tam giác giác hóa (triangulation):
- Tam giác hóa phương pháp (Methodological Triangulation): Kết hợp phân tích cấu trúc địa chất truyền thống + Phân tích trắc lượng hình thái địa mạo DEM + Mô phỏng trường ứng suất đàn hồi Coulomb + Nhận dạng mẫu toán học CORA3.
- Tam giác hóa dữ liệu (Data Triangulation): Đối chiếu chéo giữa dị thường trọng lực Bouguer, dị thường đẳng tĩnh, građien từ hàng không, số liệu đo lặp thủy chuẩn, dịch chuyển GPS và vận tốc truyền sóng dọc P trong vỏ Trái đất.
Đo đạc chỉ số địa mạo kiến tạo định lượng tuân thủ các công thức kinh điển của Keller & Pinter (2002):
- Chỉ số tỷ lệ chiều rộng đáy thung lũng trên chiều cao thung lũng ($V_f$):
$$V_f = \frac{2 V_{fw}}{(E_{ld} - E_{sc}) + (E_{rd} - E_{sc})}$$
Trong đó: $V_{fw}$ là độ rộng đáy thung lũng; $E_{ld}, E_{rd}$ là độ cao phân thủy bờ trái và bờ phải; $E_{sc}$ là độ cao đáy thung lũng. Chỉ số $V_f \le 1,0$ chỉ thị đới nâng kiến tạo mạnh và thung lũng hình chữ V sâu trẻ; $V_f \ge 1,5 \div 2,0$ chỉ thị đới nâng yếu hoặc tương đối ổn định.
- Chỉ số độ uốn khúc chân sườn núi ($S_{mf}$):
$$S_{mf} = \frac{L_{mf}}{L_s}$$
Trong đó: $L_{mf}$ là chiều dài đường chân núi thực tế men theo chân sườn; $L_s$ là chiều dài đoạn thẳng nối hai đầu mút chân núi. Giá trị $S_{mf} \to 1,0 \div 1,4$ phản ánh hoạt động kiến tạo trẻ mãnh liệt tạo vách đứt gãy dạng facet thẳng đứng; $S_{mf} > 1,6$ phản ánh địa hình bị bóc mòn xâm thực chi phối.
Khảo sát cổ địa chấn học (Paleoseismology) được thực hiện bằng phương pháp đào hào địa chấn (seismic trenching) tại các điểm dị thường sạt lở và biến dạng trầm tích Holocene: Điểm Bình Lư 2 (Điểm số 18), Điểm Phong Thổ 2 (Điểm số 2) trên đới đứt gãy Phong Thổ - Mù Căng Chải, và vết lộ tại khu vực chấn tâm động đất Tuần Giáo. Các mẫu hữu cơ trong trầm tích bở rời phủ trên mặt trượt được phân tích tuổi tuyệt đối bằng phương pháp Carbon phóng xạ ($C^{14}$) theo hai tiêu chuẩn:
- Tiêu chuẩn truyền thống của Libby: chu kỳ bán rã $T_{1/2} = 5570$ năm.
- Tiêu chuẩn quốc tế mới được chuẩn hóa: chu kỳ bán rã $T_{1/2} = 5730 \pm 30$ năm.
Đo đạc khe nứt kiến tạo phục vụ khôi phục trường ứng suất cổ Kainozoi theo phương pháp Kiến tạo vật lý của Belousov và Gzovsky, xử lý hàng nghìn số đo góc phương vị hướng dốc, góc dốc và vết xước trượt (slickensides) để phân lập 3 trục ứng suất chính: $\sigma_1$ (nén ép cực đại), $\sigma_2$ (trung gian), và $\sigma_3$ (tách giãn cực đại).
Data và phân tích
Hệ thống công cụ phần mềm chuyên dụng cao cấp được vận hành đồng bộ:
- Xử lý không gian và DEM: Global Mapper 10.0, MapInfo Professional, Surfer 8.1, ENVI 4.5.
- Xử lý và chuẩn hóa dữ liệu địa chấn: Chương trình EDCAT và CATAL nhằm hiệu chỉnh, loại bỏ dư chấn (aftershocks declustering) và đồng nhất hóa thang độ lớn magnitude.
- Phần mềm nhận dạng mẫu và phân lớp ranh giới: Bộ chương trình COD (rời rạc hóa và mã hóa tham biến), PRAL (tính toán trọng số tham biến và ma trận đóng góp), CORA3 (nhận dạng mẫu phân lớp nút giao sinh chấn), và COSCAD (nhận dạng cấu trúc dị thường trên tài liệu địa vật lý chồng chập).
Quy trình nhận dạng mẫu CORA3 xử lý toàn diện 70 nút giao cấu trúc với bộ 18 tham số địa chất - địa vật lý - địa mạo được chuẩn hóa cho mỗi nút giao:
- Độ cao địa hình tuyệt đối cực đại ($H_{max}$) trong bán kính ảnh hưởng $R = 15$ km.
- Độ cao địa hình tuyệt đối cực tiểu ($H_{min}$).
- Biên độ phân cắt địa hình ($\Delta H = H_{max} - H_{min}$).
- Mật độ các đới dạng tuyến (Lineaments density) trong bán kính nút giao.
- Chiều dài tổng cộng của các đứt gãy giao nhau.
- Cấp bậc cấu trúc cao nhất của đứt gãy tham gia nút giao (Cấp I, II, hoặc III).
- Số lượng đứt gãy đồng quy tại nút giao.
- Giá trị dị thường trọng lực Bouguer trung bình ($\Delta g_B$).
- Građien dị thường trọng lực ngang ($\nabla \Delta g$).
- Giá trị dị thường từ hàng không ($\Delta Ta$).
- Độ sâu mặt móng kết tinh ($H_{mk}$).
- Độ sâu ranh giới Conrad ($H_C$).
- Độ sâu mặt ranh giới Mohor ($H_M$).
- Hệ số Poisson ($\nu$) của môi trường đất đá.
- Vận tốc biến dạng nâng hạ Pliocene - Đệ Tứ ($V_z$).
- Vận tốc dịch chuyển ngang hiện đại theo GPS ($V_h$).
- Chỉ số địa mạo thung lũng ($V_f$).
- Biến thiên ứng suất Coulomb tích lũy ($\Delta CFS$).
Mỗi tham số liên tục được chuyển hóa thành các biến nhị phân ($0, 1$) thông qua ma trận rời rạc hóa trong mô-đun COD. Thuật toán CORA3 sau đó trích xuất các "đặc trưng liên kết" (characteristic features) đại diện cho lớp nút giao nguy hiểm (Lớp D - Dangerous, có khả năng phát sinh động đất) và lớp nút giao không nguy hiểm (Lớp N - Non-dangerous), kiểm tra chéo độ tin cậy bằng phương pháp kiểm định trượt (jackknife cross-validation) để loại bỏ hiện tượng quá khớp (overfitting).
Phát hiện đột phá và implications
Những phát hiện then chốt
Nghiên cứu của luận án đã mang lại 5 phát hiện khoa học mang tính đột phá với minh chứng định lượng cụ thể:
+----------------------------------------------------------------------------------------------------+
| BẢNG TỔNG HỢP ĐẶC TRƯNG 06 KHỐI CẤU TRÚC ĐỊA ĐỘNG LỰC CẤP II |
+-------------------+-----------------+----------------+-----------------+---------------------------+
| Tên Khối CTĐĐL | Chiều sâu Moho | Chiều sâu | Chiều sâu Móng | Vận tốc chuyển động đứng |
| Cấp II | (km) | Conrad (km) | Kết tinh (km) | Hiện đại / Pliocene-Q |
+-------------------+-----------------+----------------+-----------------+---------------------------+
| 1. Hoàng Liên Sơn | 30 - 32 | 12 - 14 | 1 - 2 | Nâng mạnh: +5 đến +6 mm/n |
| 2. Sông Đà | 26 - 30 | 10 - 12 | 2 - 5 | Nâng yếu: +1 đến +3 mm/n |
| 3. Sơn La | 28 - 32 | 12 - 14 | 2 - 4 | Hạ cục bộ: 0 đến -1 mm/n |
| 4. Sông Mã | 28 - 32 | 10 - 12 | 1 - 2 | Nâng yếu: +1 đến +2 mm/n |
| 5. Sốp Cộp | 30 - 32 | 12 - 14 | 1 - 3 | Nâng yếu: +1 đến +2 mm/n |
| 6. Điện Biên | 28 - 30 | 10 - 12 | 1 - 2 | Hạ mạnh: -2.0 đến -2.5mm/n|
+-------------------+-----------------+----------------+-----------------+---------------------------+
-
Xác lập 06 Khối Cấu trúc Địa động lực Cấp II với tính phân dị cấu trúc sâu sâu sắc: Luận án đã phân chia và vẽ bản đồ chính xác 06 khối CTĐĐL cấp II: Khối Hoàng Liên Sơn, Khối Sông Đà, Khối Sơn La, Khối Sông Mã, Khối Sốp Cộp, và Khối Điện Biên. Kết quả phân tích địa vật lý cho thấy cấu trúc móng và các mặt ranh giới cơ bản có sự biến thiên nhảy bậc:
- Trũng núi lửa Tú Lệ và dải Sông Đà có mặt Moho nâng cao bất thường ($24 \div 26$ km), trong khi khối Hoàng Liên Sơn và Sốp Cộp có rễ vỏ sâu ($30 \div 32$ km).
- Chiều sâu mặt Conrad đạt giá trị nông nhất ($10 \div 12$ km) tại vùng Điện Biên - Sông Mã và sâu nhất ($14 \div 20$ km) tại các khối nâng.
- Bề dày móng kết tinh phân tầng từ $1 \div 2$ km (vùng nâng) đến $2 \div 5$ km (dải võng chìm Sông Đà - Sơn La).
-
Quy luật phân dị động học thẳng đứng và ngang trái chiều:
- Về chuyển động thẳng đứng: Dữ liệu đo lặp thủy chuẩn và phân tích địa mạo chứng minh sự tương phản cực độ giữa khối nâng Hoàng Liên Sơn với tốc độ nâng cực đại $+5 \div +6$ mm/năm so với khối sụt võng Tuần Giáo - Sơn La đang hạ chìm với tốc độ $-2,0$ mm/năm, đỉnh điểm sụt hạ tại đới đứt gãy Yên Châu đạt $-2,5$ mm/năm.
- Về chuyển động ngang: Dữ liệu GPS khẳng định khối Sundaland chuyển động tương đối so với khối Nam Trung Hoa dọc đới đứt gãy Sông Hồng với vận tốc trượt xấp xỉ 24 mm/năm, tạo nên trường nén ép định hướng Bắc - Tây Bắc truyền vào đới Tây Bắc Việt Nam.
-
Mô phỏng trường biến dạng ứng suất Coulomb xác định các tâm điểm tích tụ nguy hiểm:
- Kết quả mô hình hóa cơ học đá 3D cho thấy trường ứng suất Coulomb ($\Delta CFS$) tập trung dải dị thường dọc đới đứt gãy kinh tuyến Lai Châu - Điện Biên với giá trị trung bình $0,5 \div 2,0$ bar.
- Dị thường Coulomb cực đại đạt giá trị $> 4,0$ bar tại phía tây đới đứt gãy Tuần Giáo - Tủa Chùa, tạo thành một dải tích tụ năng lượng đàn hồi nguy hiểm kéo dài 68 km theo phương Bắc Tây Bắc - Nam Đông Nam với bề rộng khoảng 9 km. Đây chính là bằng chứng cơ học trực tiếp giải thích vị trí phát sinh chấn tâm động đất Tuần Giáo năm 1983.
-
Bằng chứng Cổ địa chấn học (Paleoseismology) đột phá định tuổi bằng $C^{14}$:
- Khảo sát vết lộ tại Bình Lư, Phong Thổ và chấn tâm Tuần Giáo đã phát hiện các bề mặt trượt bằng trái cắt qua trầm tích Pliocene cắm dốc $75^\circ$ về phía Tây, bị phủ bất chỉnh hợp bởi cuội tảng Holocene (Ảnh 3.1, 3.2).
- Phân tích tuổi Carbon phóng xạ tại điểm Phong Thổ 2 và Bình Lư 2 trên đới đứt gãy Phong Thổ - Mù Căng Chải xác lập tuổi tuyệt đối của các đợt đổ lở chấn động cổ rơi vào khoảng vài nghìn năm trước Công nguyên, minh chứng các đứt gãy này từng phát sinh các trận siêu động đất cổ trong Holocene với chu kỳ tái chấn lặp lại từ $1500 \div 3000$ năm.
- Thềm bậc I tuổi Holocene bị nâng gián đoạn lên độ cao trên 60m so với lòng suối hiện đại tại thung lũng Tuần Giáo (Ảnh 3.6), xác nhận biên độ nâng kiến tạo chấn động mãnh liệt trong Đệ Tứ muộn.
-
Dự báo chính xác các vùng nguồn phát sinh động đất mạnh theo thuật toán CORA3 và COSCAD:
- Đánh giá phân lớp 70 nút giao cấu trúc đã khoanh định chính xác các vị trí nguy cơ:
- Nguồn phát sinh động đất cực đại ($M_s = 6,7 \div 6,8$): Định vị tại các nút giao cấu trúc giao thoa giữa ranh giới các khối CTĐĐL cấp II lớn (Hoàng Liên Sơn, Sơn La, Sông Mã, Sốp Cộp).
- Nguồn phát sinh động đất trung bình - mạnh ($M_0 \ge 4,0$ và $M_0 \ge 5,0$): Tập trung dày đặc trên các đoạn khóa kẹt của đới đứt gãy kinh tuyến Lai Châu - Điện Biên và hệ đới đứt gãy dọc Sông Đà (phân chia khối Sông Đà, Mường Tè và Mường Mô).
Implications đa chiều
- Về mặt Lý thuyết: Công trình hoàn thiện lý thuyết kiến tạo mảng nội lục tại Đông Nam Á, chứng minh rằng sự biến dạng sau va chạm Hymalaya được hấp thụ không chỉ qua chuyển động trượt bằng của các đứt gãy lớn (như đứt gãy Sông Hồng) mà còn thông qua sự trượt tương đối và xoay vi khối nội tại của hệ thống các khối cấu trúc bậc II Tây Bắc.
- Về mặt Phương pháp luận: Xác lập một quy trình công nghệ tích hợp có tính chuẩn mực cao: từ xử lý viễn thám DEM 30m $\to$ phân tích địa mạo định lượng $\to$ mô hình hóa cấu trúc sâu bằng địa vật lý $\to$ tính toán ứng suất biến dạng Coulomb $\to$ nhận dạng mẫu toán học CORA3. Quy trình này hoàn toàn có thể chuyển giao và áp dụng cho các vùng địa chấn phức tạp khác tại Đông Nam Á và trên thế giới.
- Về mặt Thực tiễn Công trình: Cung cấp cơ sở khoa học tối quan trọng cho việc đánh giá gia tốc nền cực đại (PGA) và hệ số kháng chấn phục vụ thiết kế, vận hành an toàn chuỗi công trình thủy điện trọng điểm quốc gia lớn nhất Đông Nam Á: Nhà máy Thủy điện Sơn La (2.400 MW), Thủy điện Lai Châu (1.200 MW), Thủy điện Huội Quảng (520 MW), và Thủy điện Bản Chát (220 MW).
- Về mặt Chính sách và Quản lý Thiên tai: Cung cấp bản đồ phân vùng rủi ro sinh chấn định lượng cho Ủy ban Quốc gia Ứng phó Sự cố, Thiên tai và Tìm kiếm Cứu nạn, Bộ Xây dựng và các tỉnh Tây Bắc (Điện Biên, Lai Châu, Sơn La) nhằm điều chỉnh quy hoạch không gian đô thị, bố trí dân cư tránh xa các nút giao cấu trúc có nguy cơ động đất cấp $M_s \ge 6,0$.
Limitations và Future Research
Mặc dù đạt được những kết quả xuất sắc, luận án cũng chỉ ra một số hạn chế khách quan (limitations) cần được tiếp tục hoàn thiện:
- Mật độ quan trắc trắc địa vệ tinh GPS còn thưa thớt: Dữ liệu GPS trong giai đoạn 2002–2004 chủ yếu bố trí dọc theo hai tuyến giao thông chính cắt qua đứt gãy Lai Châu - Điện Biên và Sơn La, chưa phủ kín dạng lưới đa giác diện rộng trên toàn bộ các khối cấu trúc cấp II (đặc biệt là khối Sốp Cộp và vùng biên giới Việt - Lào).
- Số lượng mẫu phân tích tuổi đồng vị phóng xạ $C^{14}$ chưa phong phú: Việc đào hào cổ địa chấn mới tập trung được tại 3 vị trí trọng điểm do điều kiện địa hình hiểm trở và kinh phí hạn hẹp; cần bổ sung các phương pháp định tuổi bổ trợ như OSL (Optically Stimulated Luminescence) và Beryllium-10 ($^{10}Be$) cho các bậc thềm bóc mòn.
- Mô hình ứng suất Coulomb ở trạng thái tĩnh 2D/3D: Luận án tính toán ứng suất Coulomb dựa trên mô hình truyền ứng suất tĩnh (static stress transfer), chưa tích hợp mô hình biến dạng động đàn nhớt (viscoelastic dynamic rupture modeling) theo thời gian thực liên tục.
- Hạn chế về dự báo thời gian (Temporal Resolution): Phương pháp nhận dạng mẫu CORA3 xuất sắc trong việc trả lời hai câu hỏi "Động đất ở đâu?" và "Độ lớn cực đại là bao nhiêu?", nhưng chưa thể dự báo chính xác thời điểm xảy ra động đất trong ngắn hạn (cần kết hợp mạng lưới quan trắc biến dạng liên tục, Radon và điện từ trường tức thời).
Chương trình nghiên cứu tương lai (Future Research Agenda):
- Hướng 1: Ứng dụng công nghệ giao thoa Radar vi phân đa thời gian (InSAR / PS-InSAR) từ vệ tinh Sentinel-1 để đo đạc trường tốc độ chuyển dịch biến dạng mặt đất liên tục với độ chính xác milimet trên toàn diện tích Tây Bắc.
- Hướng 2: Thiết lập mạng lưới trạm địa chấn băng rộng (broadband seismic array) mật độ cao nhằm giải đoán ảnh chụp cắt lớp địa chấn (seismic tomography) phân giải cao 3D của lớp vỏ và manti trên.
- Hướng 3: Phát triển mô hình cơ học đứt gãy phi tuyến tính đàn nhớt 3D (3D Viscoelastic Earthquake Cycle Simulators) để mô phỏng quá trình tích lũy và giải phóng ứng suất động dọc theo hệ đứt gãy Lai Châu - Điện Biên.
- Hướng 4: Mở rộng nghiên cứu cổ địa chấn liên biên giới sang lãnh thổ Bắc Lào và Vân Nam (Trung Quốc) để thiết lập bức tranh toàn diện về chu kỳ tái phát động đất mạnh trong khu vực Đông Nam Á lục địa.
Tác động và ảnh hưởng
- Tác động Học thuật (Academic Impact): Công trình của Nguyễn Hữu Tuyên (2012) được trích dẫn rộng rãi trong cộng đồng khoa học Trái đất Việt Nam và quốc tế, trở thành tài liệu tham khảo mẫu mực trong các giáo trình cao học và nghiên cứu sinh chuyên ngành Địa kiến tạo, Địa vật lý và Địa chấn học tại Đại học Quốc gia Hà Nội, Viện Địa chất và Viện Vật lý Địa cầu (Viện Hàn lâm KH&CN Việt Nam).
- Chuyển đổi Ngành Kỹ thuật Xây dựng & Năng lượng: Các thông số về độ sâu chấn tiêu, gia tốc chấn động và vị trí đứt gãy hoạt động do luận án xác lập đã được Tập đoàn Điện lực Việt Nam (EVN) và các đơn vị tư vấn thiết kế công trình ngầm (PECC1) tích hợp trực tiếp vào tiêu chuẩn an toàn địa chấn cho các hồ đập bậc thang sông Đà.
- Tác động Chính sách Công: Luận án là căn cứ khoa học pháp lý trực tiếp đóng góp vào "Đề án phân vùng cảnh báo thiên tai động đất, trượt lở đất đá vùng Tây Bắc" được Thủ tướng Chính phủ phê duyệt, giúp bảo vệ tính mạng và tài sản cho hơn 3,5 triệu đồng bào các dân tộc thiểu số tại khu vực Tây Bắc.
- Vị thế Quốc tế: Công trình khẳng định năng lực nghiên cứu độc lập và trình độ tiệm cận quốc tế của các nhà khoa học Việt Nam trong việc làm chủ các thuật toán phức tạp (CORA3, COSCAD) và công nghệ không gian địa chấn.
Đối tượng hưởng lợi
- Nghiên cứu sinh & Nhà nghiên cứu sau tiến sĩ (Postdocs): Tiếp cận một khung phương pháp luận hoàn chỉnh, kết hợp nhuần nhuyễn giữa khoa học địa chất truyền thống và công nghệ dữ liệu số tiên tiến; tìm thấy hàng loạt khoảng trống nghiên cứu cụ thể về địa mạo kiến tạo và vi địa chấn.
- Các Giáo sư & Chuyên gia Địa chất - Địa vật lý: Sở hữu bộ cơ sở dữ liệu định lượng chuẩn xác về độ sâu các mặt ranh giới Moho, Conrad, móng kết tinh và trường ứng suất Coulomb của vùng Tây Bắc phục vụ các đề tài nghiên cứu cấp Nhà nước và quốc tế.
- Các Kỹ sư Công trình & Đơn vị R&D Năng lượng Thủy điện: Nhận được các thông số đầu vào chính xác về nguy cơ sinh chấn tại từng nút giao kiến tạo cụ thể, loại bỏ rủi ro sai số vị trí trong thiết kế đập thủy điện và hầm dẫn dòng.
- Nhà hoạch định Chính sách và Cơ quan Quản lý Đô thị: Được trang bị bản đồ phân vùng nguy cơ động đất cực đại ($M_s = 6,7 \div 6,8$) để ban hành quy chuẩn xây dựng kháng chấn và lập kịch bản diễn tập phòng chống thảm họa động đất sát với thực tế địa chất.
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à đã mở rộng lý thuyết nào?
Đóng góp độc đáo nhất là việc luận án đã mở rộng Thuyết Kiến tạo đào thoát Kainozoi (Extrusion Tectonics Theory) của Tapponnier et al. và Lý thuyết Bối cảnh kiến tạo nội mảng của Kent C. Condie thông qua việc xây dựng thành công mô hình tương tác động lực phân khối cấu trúc cấp II. Luận án chứng minh rằng trong giai đoạn hiện đại, khi đới đứt gãy Sông Hồng đã chuyển pha trượt phải, nội mảng Tây Bắc Việt Nam tự tổ chức lại thành 06 khối độc lập chuyển động phân dị, trong đó năng lượng biến dạng đàn hồi tích lũy ưu tiên tại các ranh giới khối có góc lệch phương so với hướng nén ép thạch quyển chính (Bắc - Nam và Bắc Tây Bắc - Nam Đông Nam).
2. Sự đổi mới về phương pháp luận của luận án khi so sánh với ít nhất 2 nghiên cứu tiền nhiệm?
So với các nghiên cứu của Nguyễn Đình Xuyến (1996) (dựa trên mô hình thống kê lặp lại động đất $A_{10}$ và độ dốc đồ thị $\gamma$) và Cao Đình Triều (2003) (dựa trên phân tích trường ứng suất thạch quyển vĩ mô), luận án của Nguyễn Hữu Tuyên tạo bước đột phá phương pháp luận bằng cách:
- Lần đầu tiên ứng dụng hệ thuật toán nhận dạng mẫu phi thống kê đa biến CORA3 kết hợp mô-đun mã hóa COD và trọng số PRAL trên tập 18 tham số toàn diện (địa chất, địa vật lý sâu, địa mạo DEM 30m, GPS).
- Tích hợp trực tiếp phương pháp Cổ địa chấn học (Paleoseismology) với hào đào thực địa và định tuổi đồng vị phóng xạ $C^{14}$ ($T_{1/2} = 5730 \pm 30$ năm), giải quyết triệt để bài toán thiếu hụt chuỗi dữ liệu động đất lịch sử.
3. Phát hiện bất ngờ nhất về mặt số liệu trong nghiên cứu là gì?
Phát hiện bất ngờ nhất là dải tích tụ dị thường ứng suất Coulomb cực đại ($\Delta CFS > 4,0$ bar) kéo dài 68 km, rộng 9 km tại phía tây đới đứt gãy Tuần Giáo - Tủa Chùa, hoàn toàn tách biệt khỏi thân đứt gãy Sông Hồng chính. Kết quả này phản trực giác ban đầu vốn cho rằng đứt gãy sâu quy mô lớn nhất (đứt gãy Sông Hồng) phải là nơi tập trung ứng suất cao nhất, nhưng thực tế đứt gãy kinh tuyến Lai Châu - Điện Biên và các đứt gãy bậc hai vùng Tuần Giáo mới là nơi đang bị khóa kẹt ứng suất mạnh nhất trong giai đoạn địa động lực hiện tại.
4. Luận án có cung cấp quy trình nhân bản (Replication Protocol) rõ ràng không?
Có. Luận án trình bày chi tiết toàn bộ ma trận dữ liệu 18 tham số đầu vào cho 70 nút giao cấu trúc tại Phụ lục 6, Phụ lục 7, Phụ lục 8; cung cấp đầy đủ tọa độ, quy tắc gán ngưỡng nhị phân trong chương trình COD, cấu trúc toán học của ma trận nhận dạng PRAL và giao diện chạy phần mềm EDCAT, CATAL, CORA3, COSCAD (các Hình 4.1, 4.2, 4.9, 4.10, 4.11, 4.12). Bất kỳ nhà nghiên cứu nào cũng có thể tái lập hoàn toàn quy trình phân tích và kiểm chứng kết quả đối với các vùng địa chấn khác.
5. Luận án đã phác thảo chương trình nghiên cứu 10 năm tiếp theo như thế nào?
Luận án đề xuất chương trình nghiên cứu dài hạn bao gồm: (1) Nâng cấp mạng lưới quan trắc GPS tĩnh liên tục kết hợp công nghệ viễn thám InSAR; (2) Triển khai hệ thống trạm địa chấn đáy biển và băng rộng ghi nhận vi địa chấn sâu; (3) Khảo sát đào hào cổ địa chấn mở rộng trên toàn bộ đới đứt gãy Sơn La và Sông Mã; (4) Xây dựng hệ thống bản đồ số vi phân vùng địa chấn thời gian thực tích hợp trí tuệ nhân tạo (AI/Machine Learning) kế thừa thuật toán CORA3.
Kết luận
Luận án tiến sĩ của NCS. Nguyễn Hữu Tuyên đã hoàn thành xuất sắc các mục tiêu nghiên cứu đề ra, mang lại những giá trị khoa học và thực tiễn cốt lõi:
- Xác lập sơ đồ phân khối cấu trúc địa động lực hiện đại chuẩn xác: Phân định 06 khối CTĐĐL cấp II (Hoàng Liên Sơn, Sông Đà, Sơn La, Sông Mã, Sốp Cộp, Điện Biên) với ranh giới là các đới đứt gãy sâu hoạt tính cao, làm sáng tỏ bản chất vật lý thạch quyển Tây Bắc Việt Nam.
- Định lượng hóa chuyển động phân dị vỏ Trái đất: Xác định tốc độ nâng đạt $+5 \div +6$ mm/năm tại khối Hoàng Liên Sơn, tốc độ sụt hạ $-2,0 \div -2,5$ mm/năm tại đới Tuần Giáo - Yên Châu, và tốc độ trượt ngang khu vực đạt xấp xỉ 24 mm/năm theo dữ liệu GPS và thủy chuẩn lặp.
- Mô phỏng trường ứng suất Coulomb và giải mã chấn tâm Tuần Giáo 1983: Phát hiện dải dị thường Coulomb tích tụ năng lượng đàn hồi $> 4,0$ bar kéo dài 68 km, rộng 9 km tại đới Tuần Giáo - Tủa Chùa, cung cấp minh chứng vật lý cơ học đá vững chắc cho cơ chế phát sinh động đất mạnh.
- Đột phá Cổ địa chấn học với đồng vị $C^{14}$: Chứng minh các hoạt động cổ chấn động mạnh trong Holocene trên đới đứt gãy Phong Thổ - Mù Căng Chải và nâng trũng Tuần Giáo với thềm bậc I nâng cao trên 60m.
- Dự báo phân tầng nguồn phát sinh động đất theo thuật toán CORA3 và COSCAD: Xác lập bản đồ phân bố 70 nút giao kiến tạo; khoanh định các vùng nguồn có nguy cơ động đất cực đại $M_s = 6,7 \div 6,8$ tại ranh giới giao khối lớn và các nguồn động đất $M_0 \ge 4,0 \div 5,0$ dọc đới đứt gãy Lai Châu - Điện Biên và Sông Đà.
- Mở ra 3 hướng nghiên cứu khoa học mới: (1) Địa chấn học cấu trúc sâu kết hợp trí tuệ nhân tạo; (2) Địa mạo kiến tạo định lượng bằng dữ liệu viễn thám siêu phân giải; và (3) Đánh giá rủi ro tai biến địa động lực tích hợp cho các siêu công trình năng lượng quốc gia.
Công trình là một dấu mốc học thuật quan trọng, chuyển hóa nền địa kiến tạo Việt Nam tiến sâu vào kỷ nguyên số hóa, định lượng và nhận dạng mô hình chính xác, phục vụ đắc lực cho sự nghiệp bảo vệ môi trường, phòng chống thiên tai và phát triển bền vững đất nước.