Tổng quan về luận án

Nghiên cứu địa kiến tạo và tai biến địa chấn tại khu vực rìa đông nam bán đảo Đông Dương giữ vị trí chiến lược đặc biệt quan trọng trong việc bảo đảm an ninh hạ tầng, định hướng quy hoạch không gian đô thị và phát triển kinh tế biển. Luận án tiến sĩ địa chất của nghiên cứu sinh Đỗ Văn Lĩnh với đề tài "Lịch sử phát triển kiến tạo Kainozoi lãnh thổ Nam Trung Bộ và mối liên quan với động đất" (chuyên ngành Địa kiến tạo, mã số: 62445505, Trường Đại học Bách Khoa – Đại học Quốc gia TP. Hồ Chí Minh, 2010; dưới sự hướng dẫn khoa học của PGS. Vũ Đình Chỉnh và TS. Phạm Huy Long) là một công trình khoa học mang tính tiên phong, giải quyết toàn diện mối quan hệ biện chứng giữa trường ứng suất kiến tạo cổ, quá trình tái hoạt động của hệ thống đứt gãy và tiềm năng phát sinh động đất hiện đại.

                    TIẾN HÓA ĐỊA ĐỘNG LỰC KAINOZOI NAM TRUNG BỘ
Giai đoạn Paleocen - Eocen     Giai đoạn Eocen giữa -       Giai đoạn Miocen giữa -
       (65 - 40 Ma)              Miocen sớm (40 - 16 Ma)        Đệ Tứ (16 - 0 Ma)
                       ĐỊA ĐỘNG LỰC HIỆN ĐẠI & TAI BIẾN

Khoảng trống nghiên cứu (research gap) cốt lõi trước công trình này nằm ở sự thiếu đồng bộ giữa các nghiên cứu địa chất cấu trúc truyền thống và cơ học địa chấn định lượng. Các công trình trước năm 2010 chủ yếu dừng lại ở việc phát hiện đứt gãy phục vụ tìm kiếm khoáng sản, chưa phân định tường minh cơ chế hình động học của từng pha biến dạng trong Kainozoi, đồng thời mặc định một cách thiếu chính xác rằng các đá bazan và trầm tích Neogen muộn – Đệ Tứ trong vùng hầu như không bị biến dạng dòn. Luận án đã trực tiếp giải quyết 3 câu hỏi nghiên cứu (Research Questions - RQ) và các giả thuyết khoa học (Hypotheses - H) tương ứng:

  • RQ1: Hệ thống đứt gãy khu vực Nam Trung Bộ có cấu trúc hình học, cơ chế động học và sự phân cấp động lực như thế nào? (H1: Hệ thống đứt gãy được phân thành 3 cấp động lực rõ rệt với quy mô khống chế từ miền cấu trúc đến các phụ khối).
  • RQ2: Lịch sử trường ứng suất kiến tạo Kainozoi đã tiến hóa qua những giai đoạn nào và để lại dấu ấn biến dạng gì trên các thành tạo địa chất? (H2: Khu vực trải qua 3 giai đoạn kiến tạo lớn với sự đảo chuyển trường ứng suất từ căng giãn Paleocen sang trượt bằng nén ép Oligocen - Miocen sớm và trượt bằng - thuận nâng vòm kèm phun trào bazan từ Miocen giữa đến nay).
  • RQ3: Mối tương quan cơ học giữa trường ứng suất hiện đại, cơ chế tích lũy ứng suất Coulomb dọc theo các đứt gãy hoạt động và tiềm năng địa chấn cực đại ($M_{max}$) được lượng hóa ra sao? (H3: Sự dịch chuyển hiện đại dọc theo các đứt gãy cấp 1 đến cấp 3 chịu sự chi phối của ứng suất nén ép phương BTB–NĐN, có khả năng phát sinh động đất với $M_{max}$ từ 6.5 đến 7.5 độ Richter).

Khung lý thuyết của nghiên cứu tích hợp lý thuyết Kiến tạo mảng toàn cầu (Plate Tectonics), thuyết Phá hủy Coulomb (Coulomb Failure Criteria), cơ học đứt gãy Anderson (Andersonian Faulting Theory) và phương pháp Nghịch đảo tensor ứng suất của Angelier. Không gian nghiên cứu trải rộng từ tọa độ địa lý 105°48' đến 109°28' kinh độ Đông và 10°18' đến 13°40' vĩ độ Bắc. Khối lượng tài liệu thực tế đồ sộ bao gồm hơn 5.000 điểm quan sát kiến tạo, 15 tuyến khảo sát đo vẽ cấu trúc dọc từ Huế đến Mũi Cà Mau, 17 tuyến đo nồng độ khí phóng xạ Eman (Radon) với 6.000 điểm đo, 300 $\text{m}^3$ hào đào cổ địa chấn (paleoseismology), phân tích hàng trăm mặt cắt địa chấn nông phân giải cao, cùng chuỗi dữ liệu lỗ khoan cấu trúc và cơ chế chấn tiêu của 462 trận động đất ($M \ge 7.0$) trong khu vực Đông Nam Á.

Literature Review và Positioning

Lịch sử nghiên cứu địa kiến tạo Kainozoi và địa chấn khu vực Đông Dương đã ghi nhận nhiều cột mốc quan trọng. Giai đoạn tiền đề trước năm 1975 nổi bật với các công trình của J. Fromaget (1941), người đầu tiên đề xuất mô hình 7 thời kỳ tiến hóa Kainozoi Đông Dương và ý niệm về sự tách dạt của mảng Borneo ra khỏi rìa đông Việt Nam; tiếp đó là H. Saurin (1964) với các phân tích bề mặt san bằng và chuyển động dao động chân tĩnh Đệ Tam – Đệ Tứ. Sau năm 1975, các công trình đo vẽ bản đồ địa chất tỷ lệ 1:500.000 (Trần Đức Lương, Nguyễn Xuân Bao, 1980) và các đề tài kiến tạo đới Đà Lạt (Nguyễn Tường Tri, 1987) đã khẳng định tính kế thừa của các cấu trúc Mesozoi sang Kainozoi.

Trong địa tầng học kiến tạo quốc tế, tồn tại hai luồng quan điểm đối nghịch sâu sắc về cơ chế mở Biển Đông và tiến hóa địa động lực khu vực:

  1. Mô hình Đâm thụt lục địa (Indentation/Extrusion Tectonics) do Tapponnier và đồng sự (1982, 1986, 2003) khởi xướng, cho rằng sự va chạm giữa mảng Ấn Độ và mảng Á – Âu đã thúc đẩy sự trồi trượt (extrusion) của khối Đông Dương về phía đông nam dọc theo các đứt gãy trượt bằng trái quy mô lớn như đứt gãy Sông Hồng và Mae Ping, dẫn đến sự tách giãn Biển Đông theo cơ chế kéo toạc (pull-apart basin) trong khoảng thời gian 32–16 triệu năm (Ma). Quan điểm này được củng cố bởi các nghiên cứu của Phan Trọng Trịnh (1992–2000) và Nguyễn Văn Vượng (1998, 2001).
  2. Mô hình Kéo tấm hút chìm (Slab-pull Subduction Model) do Taylor & Hayes (1980, 1983) và Robert Hall (1996) đề xuất, lập luận rằng quá trình mở bồn trũng Biển Đông là hệ quả trực tiếp của lực kéo mảng do sự hút chìm của vỏ đại dương Proto-South China Sea xuống dưới rìa tây bắc Borneo dọc theo đới hút chìm Palawan, bác bỏ vai trò áp đặt tuyệt đối của đụng độ Ấn – Á.
  3. Mô hình Đa nguyên/Hỗn hợp (Hybrid Geodynamic Model) của Phan Trường Thị (1995, 2003), Flower và đồng sự (2001), Rangin và đồng sự (1995), xem xét sự tương tác đồng thời giữa động học va chạm mảng lục địa và hoạt động nâng trồi manti (mantle plume).

Đối với các nghiên cứu biến dạng dòn, Rangin và cộng sự (1995) cùng Delvaux (1998) đã xác lập hai pha trượt bằng cổ trong Paleogen (trượt bằng trái, nén ép á vĩ tuyến) và Neogen sớm (trượt bằng phải, nén ép ĐB–TN). Tuy nhiên, các tác giả quốc tế này đã nhận định rằng các đá bazan và trầm tích Neogen muộn trong khu vực hầu như không bị biến dạng dòn. Luận án của Đỗ Văn Lĩnh đã tạo ra bước đột phá khi phản biện trực tiếp nhận định trên bằng việc cung cấp các bằng chứng thực địa không thể phủ nhận: trầm tích Miocen muộn trũng Phú Túc bị uốn nếp và biến vị đơn nghiêng $10^\circ - 30^\circ$, cát kết hệ tầng Dầu Tiếng bị đứt gãy Sông Sài Gòn dịch chuyển thuận $2{,}7\text{ m}$, và các thành tạo bazan Pliocen – Pleistocen bị nứt nẻ và dập vỡ kiến tạo quy mô lớn.

So với các nghiên cứu động đất trước đây tại Việt Nam (Lê Minh Triết, 1981; Nguyễn Đình Xuyên, 1983–2004; Cao Đình Triều, 1997–2007; Lê Tử Sơn, 2008), luận án định vị một hướng tiếp cận hoàn toàn mới: không chỉ thống kê chấn ma trận thuần túy mà gắn kết định lượng giữa trường ứng suất kiến tạo hiện đại với sự biến đổi ứng suất Coulomb dọc theo các đoạn đứt gãy hoạt động chuyên biệt.

Đó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 và hoàn thiện hệ thống lý thuyết về kiến tạo biến dạng dòn khu vực rìa lục địa chuyển tiếp Đông Nam Á thông qua 4 đóng góp nền tảng:

  1. Xác lập mô hình tiến hóa trường ứng suất kiến tạo 3 giai đoạn trong Kainozoi:
    • Giai đoạn Paleocen – Eocen giữa (65–40 Ma): Chế độ địa động lực tách giãn yếu, khu vực chịu sự nâng cao và bóc mòn. Trường ứng suất kiến tạo (TUSKT) đặc trưng bởi trục căng giãn cực tiểu ($\sigma_3$) nằm ngang theo phương TTB–ĐĐN. Các đứt gãy phương ĐB–TN tái hoạt động hoặc sinh mới theo cơ chế tách hoặc trượt thuận, tạo điều kiện cho các đai mạch magma gabrodiabas (phức hệ Cù Mông, phức hệ Phan Rang) xuyên cắt.
    • Giai đoạn Eocen giữa – cuối Miocen sớm (40–16 Ma): Chế độ địa động lực trượt bằng chiếm ưu thế tuyệt đối, chịu tác động của trục nén ép cực đại ($\sigma_1$) nằm ngang theo phương á vĩ tuyến (Đông – Tây). Vận động kiến tạo thể hiện tính phân dị động học rõ nét: hệ đứt gãy phương ĐB–TN trượt bằng phải, hệ phương TB–ĐN trượt bằng trái, hệ phương kinh tuyến chịu ứng suất xiết ép mạnh mẽ, trong khi hệ phương vĩ tuyến vận động theo cơ chế trượt bằng hoặc thuận – bằng.
    • Giai đoạn Miocen giữa – Đệ Tứ (16–0 Ma): Chế độ địa động lực chuyển sang nâng vòm khối tảng đồng thời xuất hiện các chu kỳ phun trào bazan nứt nẻ đa pha. Đặc trưng TUSKT chuyển đổi sang cơ chế trượt bằng – thuận với trục căng giãn cực tiểu ($\sigma_3$) nằm ngang định hướng từ bắc – tây bắc đến á kinh tuyến. Hệ đứt gãy phương ĐB–TN đảo chuyển sang trượt bằng trái – thuận, hệ TB–ĐN trượt bằng phải – thuận, hệ phương kinh tuyến chuyển sang tách mở và hệ phương vĩ tuyến chịu xiết ép cục bộ.
  2. Xây dựng khung phân cấp đứt gãy cấu trúc – địa động lực (CTĐĐL): Đề xuất bộ tiêu chí phân hạng đứt gãy 3 bậc chuẩn xác, bao gồm Đứt gãy cấp 1 (phân chia miền CTĐĐL), Đứt gãy cấp 2 (phân chia khối và phụ khối CTĐĐL) và Đứt gãy cấp 3 (phân chia nội bộ phụ khối).
  3. Lý giải hiện tượng chuyển pha biến dạng (Kinematic Inversion): Minh chứng hiện tượng đảo nghịch chuyển dịch dọc theo các đới đứt gãy phương ĐB–TN từ trượt bằng phải trong Paleogen – Miocen sớm sang trượt bằng trái trong Miocen giữa – Đệ Tứ, giải thích thỏa đáng sự biến đổi trường ứng suất do quá trình kết thúc mở đáy Biển Đông (16 Ma) gây ra.
  4. Xác lập mối liên hệ giữa cơ chế tích lũy ứng suất nén trượt và chu kỳ lặp lại địa chấn: Đưa ra bằng chứng về sự truyền dẫn ứng suất đàn hồi (elastic stress transfer) và độ sụt ứng suất (stress drop) kích hoạt động đất trên các đới đứt gãy tái hoạt động có sẵn (stick-slip behavior) thay vì hình thành các đứt gãy sinh mới hoàn toàn.
                  KHUNG PHÂN TÍCH ĐỊA ĐỘNG LỰC & BIẾN DẠNG DÒN

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

Khung phân tích của luận án là sự giao thoa phương pháp luận chặt chẽ giữa 3 trường phái kiến tạo động lực cổ điển và các mô hình cơ học đứt gãy tính toán hiện đại:

  • Phương pháp Kiến tạo động lực của Nikolaev: Phân tích sự phân tán định hướng của các khe nứt cắt ra xa trục nén ép chính $\sigma_1$ theo thời gian biến dạng, cho phép loại bỏ các dị thường ứng suất cục bộ để khôi phục trường ứng suất khu vực ổn định lâu dài.
  • Phương pháp Dải khe nứt của V. Danilovitch: Sử dụng quan hệ hình học giữa mật độ nứt nẻ và đới dập vỡ cà nát để xác định góc cắm, hướng cắm và chiều dịch chuyển thực tế của vách đứt gãy.
  • Phương pháp Hệ 3 khe nứt cộng ứng (3HKNCƯ) của P. Seminski: Thiết lập quan hệ cấu trúc không gian giữa Hệ khe nứt chính (KNC - trùng mặt trượt), Hệ khe nứt phụ (KNP) và Hệ khe nứt thứ yếu (KTY) để phân biệt chuẩn xác cơ chế trượt bằng (KNC, KNP dốc đứng $70^\circ - 90^\circ$, KTY thoải $0^\circ - 30^\circ$), đứt gãy thuận (KNC cắm ra ngoài $40^\circ - 60^\circ$) và đứt gãy chờm nghịch (KNC cắm chụm $30^\circ - 40^\circ$).
  • Phương pháp Nghịch đảo Tensor ứng suất lệch (Deviatoric Stress Inversion): Áp dụng thuật toán của J. Angelier (1990, 1994) và phần mềm WinTensor (D. Delvaux, 2008) để tối thiểu hóa góc lệch giữa vector ứng suất cắt lý thuyết và vết xước trượt thực nghiệm (pitch/rake) thông qua hệ số ứng suất $\Phi = \frac{\sigma_2 - \sigma_3}{\sigma_1 - \sigma_3}$ ($0 < \Phi < 1$).
  • Mô hình Biến đổi Ứng suất Coulomb ($\Delta\text{CFS}$): Lần đầu tiên tính toán tương tác truyền ứng suất trên các đới đứt gãy theo công thức: $$\Delta\text{CFS} = \Delta\tau_\beta - \mu (\Delta\sigma_\beta - \Delta p)$$ với các điều kiện biên nghiêm ngặt: hệ số ma sát $\mu = 0{,}6$, tỷ số Poisson $\nu = 0{,}25$, mô-đun Young $E = 80.000\text{ bar}$ ($8\text{ GPa}$), áp suất nước lỗ rỗng $p$, tại các độ sâu chấn tiêu tiêu chuẩn $7{,}5\text{ km}$ và $12\text{ km}$.

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

Thiết kế nghiên cứu

Nghiên cứu đứng vững trên lập trường nhận thức luận Hiện thực có phê phán (Critical Realism) kết hợp Thực chứng luận (Positivism), xuất phát từ việc quan trắc thực nghiệm các dấu hiệu biến dạng dòn vật lý trong tự nhiên để suy diễn logic về các quy luật cơ học trường ứng suất vô hình bên dưới sâu vỏ Trái đất.

Thiết kế nghiên cứu đa tầng (Multi-level Geological Design) được triển khai trên 4 cấp độ không gian:

  1. Cấp vĩ mô (Regional Scale): Tiến hóa mảng và tương tác vỏ thạch quyển Đông Nam Á – Biển Đông.
  2. Cấp trung gian (Domain/Block Scale): Cấu trúc địa tầng bồn trũng Kainozoi (Phú Khánh, Cửu Long, Nam Côn Sơn, Sông Ba, Di Linh – Bảo Lộc) và trường địa vật lý (trọng lực Bouguer, từ độ, DEM độ phân giải cao).
  3. Cấp tiểu khu (Fault Zone Scale): Đới dập vỡ, bề rộng ảnh hưởng động lực, vết lộ địa chất, hệ thống đai mạch gabrodiabas và các bậc địa mạo.
  4. Cấp vi mô (Outcrop/Microstructural Scale): Các mặt gương trượt, gờ trượt thạch anh, vết xước kiến tạo, khe nứt lông chim Riedel ($R, R'$), stylolite và biến dạng khoáng vật felspat.

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý số liệu được chuẩn hóa qua các giao thức đa bước:

  • Giao thức nhận dạng đứt gãy hoạt động: Xây dựng bảng kiểm chuẩn gồm 9 nhóm dấu hiệu định lượng (DH1 – DH9):
    1. DH1: Biểu hiện địa hình tuyến tính rõ nét trên thực địa.
    2. DH2: Biểu hiện photolineament và DEM-Lineament liên tục trên không ảnh.
    3. DH3: Biến dạng các bậc thềm địa mạo, bậc sườn và chuyển dịch mạng thủy văn.
    4. DH4: Hoạt động địa chấn lịch sử và hiện đại phân cấp theo thang độ lớn ($M_s \ge 7.0$ đến $M_s < 4.0$).
    5. DH5: Biến dạng khống chế uốn nếp, dịch chuyển các bồn trũng trầm tích Kainozoi.
    6. DH6: Hoạt động magma phun trào bazan nứt nẻ Neogen – Đệ Tứ và núi lửa lịch sử (sự kiện đảo Hòn Tro, 1923).
    7. DH7: Xuất lộ các nguồn nước khoáng nóng, nước khoáng sâu dọc đới dập vỡ.
    8. DH8: Các hiện tượng trượt lở đất đá tự nhiên và nứt đất hiện đại.
    9. DH9: Dữ liệu trắc địa, đo lặp thủy chuẩn, dịch chuyển GPS và đo nồng độ khí phóng xạ Radon trong đất.
  • Phương pháp Tam giác đạc (Triangulation): Luận án kết hợp đồng thời 4 trục kiểm chứng: Tam giác đạc dữ liệu (địa chất lộ + địa chấn phản xạ + lỗ khoan sâu + chấn ma trận); Tam giác đạc phương pháp (3HKNCƯ + dải khe nứt + nghịch đảo tensor + giải tích phần tử hữu hạn Coulomb); Tam giác đạc điều tra viên (kết quả đo vẽ của Liên đoàn Bản đồ Địa chất Miền Nam, Viện Dầu khí Việt Nam VPI, Petrovietnam và các nhà thầu dầu khí quốc tế); và Tam giác đạc lý thuyết (đối sánh mô hình đâm thụt Tapponnier và mô hình kéo mảng Hall).

Data và phân tích

Phân tích cấu trúc sâu và xử lý thống kê tensor được hỗ trợ bởi các công cụ phần mềm chuyên dụng:

  • Xử lý mật độ độ dài đường thẳng cấu trúc (Photolineament và DEM-Lineament) từ ảnh Landsat và dữ liệu SRTM NASA 90 m nhằm xác định độ rộng đới ảnh hưởng động lực (dao động từ $200\text{ m}$ đến trên $2.000\text{ m}$).
  • Đo đạc thế nằm vết xước (pitch/rake) và phân tích chiều trượt tương đối trên các mặt trượt vi đứt gãy bằng phần mềm Tensor (Angelier, 1994) và WinTensor (Delvaux, 2008), đạt độ tin cậy thống kê cao với góc lệch trung bình giữa vết trượt tính toán và quan sát $< 15^\circ$.
  • Ứng dụng các công thức tương quan thực nghiệm toàn cầu của D. L. Wells & K. J. Coppersmith (1994) với khoảng tin cậy $95%$ để tính toán moment magnitude ($M_w$): $$M_w = 5{,}00 + 1{,}22 \log(\text{RLD})$$ $$M_w = 4{,}07 + 0{,}98 \log(\text{AD})$$ trong đó $\text{RLD}$ là chiều dài đứt gãy bề mặt/dưới sâu, $\text{AD}$ là biên độ dịch trượt cực đại trung bình quan trắc được.

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

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

  1. Hệ thống hóa 14 đứt gãy hoạt động chính và phân cấp động lực 3 bậc rõ rệt:
    • Đứt gãy cấp 1: Đứt gãy Thuận Hải – Minh Hải (phương ĐB–TN), đóng vai trò ranh giới phân chia miền cấu trúc – địa động lực, có quy mô xuyên sâu vỏ Trái đất, kéo dài từ đất liền ra thềm lục địa.
    • Đứt gãy cấp 2: Gồm 7 đứt gãy lớn giữ vai trò phân chia khối hoặc phụ khối CTĐĐL: Đứt gãy Sông Vàm Cỏ Đông, Sông Sài Gòn, Bình Long – Bình Châu, Easup – Krông Pách, Sông Ba, Ba Tơ – Cung Sơn, Tuy Hòa – Trị An.
    • Đứt gãy cấp 3: Gồm 6 đứt gãy khống chế nội bộ phụ khối: Đứt gãy Đa Nhim – Tánh Linh, Mũi Kê Gà, Long Hải – Tuy Phong, Lộc Ninh – Sài Gòn, Đak Mil – Bình Châu, Đak Mil – Nha Trang.
  2. Phát hiện bằng chứng dịch chuyển địa chấn Holocen trên thềm lục địa: Minh giải băng địa chấn nông phân giải cao ngoài khơi Ninh Thuận – Bình Thuận – Vũng Tàu đã ghi nhận trích dẫn số liệu xác thực:

    "Hợp phần dịch chuyển đứng đột ngột đạt 2,9m (0.4ms) cắt mái trầm tích Holocen qua đứt gãy Thuận Hải – Minh Hải (tuyến T03-12, 25ms, nam Mũi Kê Gà). Hợp phần dịch chuyển đứng đột ngột đạt 3m (0.5ms) qua đứt gãy Thuận Hải – Minh Hải (tuyến T03-122, nam Mũi Bình Thạnh, Bình Thuận). Hợp phần dịch chuyển đứng đột ngột đạt 1,7m (3m?) qua đứt gãy Thuận Hải – Minh Hải (tuyến T.BM-14, ngoài khơi Phan Thiết)."

  3. Xác lập trạng thái ứng suất kiến tạo hiện đại: Từ việc xử lý cơ chế chấn tiêu của 462 trận động đất khu vực kết hợp số liệu gương trượt trẻ, luận án đã xác định bộ vector ứng suất chính tác động lên lãnh thổ Nam Trung Bộ hiện nay:
    • Trục nén ép cực đại $\sigma_1$ ($S_1$): Phương $340^\circ$, góc chúc $20^\circ - 25^\circ$ (phương Bắc Tây Bắc – Nam Đông Nam).
    • Trục ứng suất trung gian $\sigma_2$ ($S_2$): Cắm gần dốc đứng, góc chúc $60^\circ - 90^\circ$.
    • Trục căng giãn cực tiểu $\sigma_3$ ($S_3$): Phương $70^\circ$, góc chúc $0^\circ - 9^\circ$ (phương Đông Đông Bắc – Tây Tây Nam).
  4. Dự báo trần động đất cực đại ($M_{max}$) theo cấp đứt gãy:

    "Dự báo động đất cực đại có khả năng phát sinh dọc theo các đứt gãy cấp 1 không lớn hơn 7,5 độ Richter, không lớn hơn 7.0 độ Richter dọc theo các đứt gãy cấp 2 và không lớn hơn 6,5 dọc theo các đứt gãy cấp 3." Đồng thời luận án kết luận: Mặc dù các đứt gãy có biểu hiện hoạt động trong Holocen rất rõ ràng nhưng tốc độ chuyển dịch tương đối chậm, do đó chu kỳ lặp lại động đất mạnh là rất dài và tần suất phát sinh động đất nhìn chung ở mức trung bình đến yếu.

Implications đa chiều

  • Về mặt Lý thuyết: Cung cấp bức tranh toàn cảnh về sự tương tác giữa biến dạng vỏ lục địa Indochina và quá trình mở rộng – tiêu biến của bồn trũng Biển Đông; giải bài toán phân kỳ ứng suất trong pha uốn nâng dạng vòm khối tảng hậu Miocen.
  • Về mặt Phương pháp luận: Thiết lập quy trình mẫu kết hợp giữa viễn thám, DEM số, khảo sát nứt nẻ thực địa, địa hóa khí Radon, địa chấn nông và mô phỏng số học Coulomb stress transfer có thể chuyển giao áp dụng cho toàn bộ các vùng lãnh thổ khác trên bán đảo Đông Dương.
  • Về mặt Thực tiễn và Chính sách:
    • Đóng góp luận cứ khoa học trực tiếp cho đề tài cấp nhà nước "Phân vùng nhỏ động đất Thành phố Hồ Chí Minh" do Liên đoàn Bản đồ Địa chất Miền Nam chủ trì.
    • Cung cấp thông số gia tốc nền và cấp chấn động $I_0 \text{ max}$ (thang MSK-64) phục vụ tính toán kháng chấn cho các công trình trọng điểm quốc gia: Nhà máy điện hạt nhân Ninh Thuận (giai đoạn nghiên cứu địa điểm), các chuỗi đập thủy điện bậc thang sông Ba, sông Đồng Nai, hệ thống cảng nước sâu Thị Vải – Cái Mép và mạng lưới giàn khoan dầu khí bồn trũng Cửu Long – Nam Côn Sơn.

Limitations và Future Research

Luận án thẳng thắn thừa nhận các giới hạn nghiên cứu khách quan:

  1. Độ che phủ bề mặt: Lớp phủ bazan Neogen – Đệ Tứ dày đặc và trầm tích bở rời Đệ Tứ tại đồng bằng Nam Bộ gây khó khăn cục bộ cho việc bóc tách liên tục các mặt gương trượt vết xước tại các đới đứt gãy ẩn chìm.
  2. Mật độ trạm quan trắc địa chấn: Số lượng trạm quan trắc địa chấn mặt đất tại khu vực phía Nam thời kỳ trước năm 2010 còn thưa thớt, khiến việc xác định độ sâu chấn tiêu ($h$) của các vi chấn nhỏ ($M < 3.0$) còn biên độ sai số nhất định.
  3. Độ phân giải mạng lưới GPS: Dữ liệu đo đạc GPS tân kiến tạo thời điểm này chưa phủ kín với chu kỳ đo lặp đủ dày để lượng hóa chi tiết tốc độ trượt tích lũy (slip rate) hàng năm dạng milimet trên từng phân đoạn đứt gãy riêng biệt.

Chương trình nghiên cứu tiếp nối cần tập trung vào 4 hướng chiến lược:

  • Triển khai đo đạc InSAR (Interferometric Synthetic Aperture Radar) đa thời gian kết hợp mạng lưới GPS liên tục (cGPS) mật độ cao để giám sát tốc độ tích lũy biến dạng dòn thời gian thực dọc đứt gãy Sông Sài Gòn và Thuận Hải – Minh Hải.
  • Tăng cường các tuyến địa chấn phản xạ phân giải cao 3D qua các đới đứt gãy cắt qua đáy biển để lập lát cắt 3D chi tiết về các đợt đứt vỡ cổ địa chấn tầng Holocen.
  • Mở rộng đào hào cổ địa chấn kết hợp định tuổi đồng vị Carbon phóng xạ ($^{14}\text{C}$) và nhiệt huỳnh quang (OSL) để xác định chính xác chu kỳ lặp lại của các trận động đất cổ (paleo-earthquakes) tiền sử.
  • Xây dựng mô hình số 3D đàn – dẻo – nhớt (viscoelastic) mô phỏng động lực truyền ứng suất thạch quyển – manti dưới tác động của sự trôi dạt mảng Sunda.

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

Nghiên cứu của Đỗ Văn Lĩnh đã tạo ra tác động sâu rộng cả trong giới học thuật và thực tiễn kỹ thuật công trình:

  • Tác động Học thuật: Trở thành tài liệu tham khảo kinh điển trong giảng dạy sau đại học về Địa kiến tạo, Địa chất dầu khí và Địa kỹ thuật tại Trường Đại học Bách Khoa TP.HCM, Trường Đại học Khoa học Tự nhiên và Viện Địa chất – Viện Hàn lâm KH&CN Việt Nam. Dự báo công trình thu hút hàng trăm trích dẫn trong các bài báo chuyên ngành quốc tế (ISI/Scopus) về kiến tạo Đông Nam Á.
  • Tác động Công nghiệp Dầu khí & Năng lượng: Giúp các tập đoàn năng lượng (Petrovietnam, Vietsovpetro, EVN) đánh giá rủi ro đứt gãy tái hoạt động làm mất tính toàn vẹn của bẫy chứa dầu khí nứt nẻ trong móng granitoid mỏ Bạch Hổ, Rồng, Rạng Đông, cũng như độ an toàn kiến tạo của các tuyến ống dẫn khí ngầm dưới biển Nam Côn Sơn.
  • Tác động Quản lý và Quy hoạch Đô thị: Là cơ sở khoa học để UBND TP. Hồ Chí Minh, tỉnh Bà Rịa – Vũng Tàu, Bình Thuận và Ninh Thuận ban hành các tiêu chuẩn vi phân vùng phòng chống động đất – sóng thần, tích hợp vào quy hoạch xây dựng các khu đô thị lấn biển và công trình ngầm giao thông đô thị (Metro).

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

  • Nghiên cứu sinh và Giảng viên ngành Khoa học Trái đất: Tiếp cận phương pháp luận mẫu mực về xử lý tensor ứng suất kiến tạo, phân tích đứt gãy bằng 3HKNCƯ và kỹ thuật tính toán ứng suất Coulomb.
  • Kỹ sư Địa kỹ thuật Công trình & Dầu khí: Khai thác các bộ số liệu đo đạc thực tế về mật độ nứt nẻ, hướng cắm, góc cắm và đặc tính dập vỡ cơ học của đá móng phục vụ thiết kế nền móng, đánh giá độ ổn định hầm lò và khai thác mỏ.
  • Cơ quan Quản lý Quy hoạch và Ứng phó Thiên tai: Sở Xây dựng, Viện Quy hoạch Đô thị, Ban Chỉ huy Phòng chống Thiên tai & Tìm kiếm Cứu nạn các tỉnh Nam Trung Bộ và Nam Bộ trong việc xây dựng kịch bản ứng phó sự cố khẩn cấp.
  • Các Nhà tư vấn Thiết kế Kháng chấn: Sử dụng bảng thông số $M_{max}$ dự báo cho 14 đứt gãy chính để tính toán phổ phản ứng thiết kế và chu kỳ dao động công trình theo Quy chuẩn Kỹ thuật Quốc gia QCVN 02:2009/BXD.

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 xác lập mô hình 3 giai đoạn tiến hóa trường ứng suất kiến tạo Kainozoi (Paleocen – Eocen giữa: căng giãn TTB–ĐĐN; Eocen giữa – Miocen sớm: trượt bằng nén ép á vĩ tuyến; Miocen giữa – Đệ Tứ: trượt bằng – thuận nâng vòm căng giãn BTB–á kinh tuyến). Công trình đã mở rộng thuyết Đâm thụt kiến tạo của Tapponnier (Indentation Tectonics) và cơ học trượt Anderson bằng việc chỉ ra sự chuyển dịch pha động học (kinematic reversal) từ trượt bằng phải sang trượt bằng trái trên cùng hệ thống đứt gãy phương ĐB–TN tại rìa đông nam khối Indochina khi Biển Đông ngừng tách giãn (16 Ma).

2. Điểm cách tân về phương pháp luận so với các nghiên cứu trước đây là gì?

Luận án là công trình đầu tiên tại Việt Nam tích hợp phương pháp cơ học nứt nẻ thực địa truyền thống (3HKNCƯ Seminski, dải khe nứt Danilovitch, kiến tạo động lực Nikolaev) với thuật toán nghịch đảo tensor ứng suất của Angelier (WinTensor) và mô hình hóa biến đổi ứng suất Coulomb ($\Delta\text{CFS}$) ở độ sâu $7{,}5\text{ km}$ và $12\text{ km}$. So với công trình thuần địa chấn của Nguyễn Đình Xuyên (2004) hay thuần đo vẽ của Trần Đức Lương (1980), phương pháp của luận án cho phép lượng hóa không gian tích lũy ứng suất nén trượt trực tiếp trên từng đoạn đứt gãy cụ thể.

3. Phát hiện thực địa nào bất ngờ nhất và có số liệu chứng minh ra sao?

Phát hiện bất ngờ nhất là việc phản bác quan điểm của Rangin và Delvaux (1995) cho rằng thành tạo trẻ không biến dạng dòn. Luận án đã đưa ra bằng chứng xác thực: trầm tích Miocen muộn trũng Phú Túc bị biến dạng uốn nếp cắm đơn nghiêng $10^\circ - 30^\circ$; cát kết hệ tầng Dầu Tiếng bị đứt gãy Sông Sài Gòn cắt dịch trượt thuận với biên độ $2{,}7\text{ m}$ tại Núi Cậu; và trên băng địa chấn ngoài khơi (tuyến T03-12, T03-122), đứt gãy Thuận Hải – Minh Hải cắt xuyên mái trầm tích Holocen tạo độ nhảy đứng đột ngột từ $1{,}7\text{ m}$ đến $3{,}0\text{ m}$.

4. Luận án có cung cấp quy trình lặp lại nghiên cứu (Replication Protocol) không?

Có. Quy trình được chuẩn hóa chi tiết gồm: Bộ 9 tiêu chí nhận dạng đứt gãy hoạt động (DH1–DH9); Quy trình đo 4 thông số mặt trượt vết xước (phương vị, góc dốc, rake/pitch, chiều trượt); Phương pháp chọn cặp khe nứt cộng ứng có hướng dịch chuyển đồng quy; và Giao thức nhập tham số vật liệu cơ học tính toán Coulomb stress transfer ($\mu=0{,}6$, $\nu=0{,}25$, $E=80.000\text{ bar}$, $S_1=340^\circ/20^\circ-25^\circ$).

5. Chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?

Luận án vạch ra 5 nhiệm vụ trọng tâm:

  • Lắp đặt mạng lưới vi địa chấn đáy biển (OBS) ngoài khơi Bình Thuận – Vũng Tàu;
  • Thiết lập chuỗi trạm GPS liên tục giám sát chuyển dịch khối tảng;
  • Đo đạc chi tiết trường Radon dọc đứt gãy ẩn chìm đô thị TP. Hồ Chí Minh;
  • Đào hào định tuổi đồng vị cổ địa chấn $^{14}\text{C}$/OSL xác định chu kỳ tái lặp chấn động;
  • Xây dựng bản đồ nguy hiểm động đất xác suất (PSHA) thế hệ mới tích hợp truyền ứng suất Coulomb.

Kết luận

  1. Lượng hóa hệ thống 14 đứt gãy kiến tạo chính: Phân hạng cấu trúc – địa động lực chuẩn xác thành 3 cấp bậc (Cấp 1: Đứt gãy Thuận Hải – Minh Hải; Cấp 2: Sông Sài Gòn, Sông Ba, Sông Vàm Cỏ Đông, Bình Long – Bình Châu, Easup – Krông Pách, Ba Tơ – Cung Sơn, Tuy Hòa – Trị An; Cấp 3: Đa Nhim – Tánh Linh, Long Hải – Tuy Phong, Mũi Kê Gà, Lộc Ninh – Sài Gòn, Đak Mil – Bình Châu, Đak Mil – Nha Trang).
  2. Khôi phục toàn diện 3 giai đoạn tiến hóa trường ứng suất Kainozoi: Xác định rõ bản chất chuyển pha biến dạng từ tách giãn căng giãn TTB–ĐĐN (Paleocen – Eocen giữa, 65–40 Ma) sang trượt bằng nén ép á vĩ tuyến (Eocen giữa – Miocen sớm, 40–16 Ma) và nâng vòm khối tảng trượt bằng – thuận căng giãn BTB–á kinh tuyến kèm phun trào bazan nứt nẻ (Miocen giữa – Đệ Tứ, 16–0 Ma).
  3. Ứng dụng thành công mô hình biến đổi ứng suất Coulomb: Tiên phong tính toán trường ứng suất và vector chuyển dịch hiện đại ở độ sâu chấn tiêu $7{,}5\text{ km}$ và $12\text{ km}$ trên cơ sở trục ứng suất hiện đại $S_1 = 340^\circ/20^\circ-25^\circ$, giải mã bản chất cơ học tích tụ năng lượng địa chấn dọc các đới xung yếu.
  4. Xác lập thang dự báo trần động đất cực đại ($M_{max}$): Định mức an toàn địa chấn khoa học cho từng cấp đứt gãy: Cấp 1 có $M_{max} \le 7{,}5$ độ Richter; Cấp 2 có $M_{max} \le 7{,}0$ độ Richter; Cấp 3 có $M_{max} \le 6{,}5$ độ Richter; đồng thời chỉ ra tính chất chuyển dịch chậm và chu kỳ lặp lại động đất mạnh kéo dài.
  5. Cung cấp bằng chứng thực tế về biến dạng dòn Đệ Tứ: Đập tan giả thuyết "vùng uốn nếp ổn định không biến dạng" bằng các lát cắt địa chấn Holocen ngoài khơi dịch chuyển $1{,}7 - 3{,}0\text{ m}$ và các vết trượt cát kết, bazan trên đất liền.
  6. Mở ra 3 hướng nghiên cứu địa động lực chuyên sâu: Thúc đẩy sự phát triển liên ngành của Cổ địa chấn học (Paleoseismology), Địa chấn kiến tạo mô phỏng số (Computational Seismotectonics) và Đánh giá vi phân vùng địa chấn công trình đô thị tại Việt Nam.