Tổng quan về luận án

Khu vực Biển Đông giữ vị trí xung yếu trong cấu trúc địa động lực Đông Nam Á, chịu sự tương tác phức tạp giữa ba mảng thạch quyển kiến tạo quy mô toàn cầu: mảng Âu - Á ở phía Tây, mảng Thái Bình Dương - Biển Philippine ở phía Đông và mảng Ấn - Úc ở phía Nam. Đặc thù chuyển dịch tương đối với vận tốc biến thiên từ 2 đến 15 nm/năm đã tích tụ trường ứng suất kiến tạo khổng lồ dọc theo các đới hút chìm hoạt tính cao như rãnh sâu Manila (Manila Trench), tiềm ẩn nguy cơ bùng phát động đất cực lớn và sóng thần hủy diệt đe dọa trực tiếp đến hơn 3.260 km bờ biển Việt Nam cùng cơ sở hạ tầng kinh tế biển.

Nghiên cứu địa chấn học Biển Đông trước đây bộc lộ một khoảng trống học thuật (research gap) nghiêm trọng: hầu hết các công trình (Phạm Văn Thục, 1979; Nguyễn Hồng Phương, 1998; Trần Thị Mỹ Thành, 2002; Cao Đình Triều và nnk, 2009a) đều sử dụng các tập dữ liệu địa chấn lịch sử kết hợp số liệu quan trắc máy chưa qua chuẩn hóa thang độ lớn (magnitude), hoặc áp dụng các mô hình xác suất tiệm cận cổ điển (như phân bố Gumbel I, Gumbel III) với bước nhảy thời gian tùy tiện và không chỉ ra được mốc thời gian dự báo tường minh. Nhằm khắc phục triệt để các hạn chế mang tính cố hữu này, công trình tiến sĩ địa vật lý của tác giả Vũ Thị Hoãn (2018) đã tiên phong xây dựng một danh mục động đất độc lập, đồng nhất hoàn toàn về thang độ lớn mômen ($M_w$) và ứng dụng phân bố cực trị tổng quát (Generalized Extreme Value - GEV) nhằm lượng hóa rủi ro địa chấn cực đại.

Các câu hỏi nghiên cứu và giả thuyết khoa học được xác lập có hệ thống:

  • RQ1: Làm thế nào để loại bỏ triệt để hiện tượng bão hòa thang đo và thiết lập một danh mục địa chấn đồng nhất chuẩn $M_w$ từ các nguồn dữ liệu đa dạng trong giai đoạn 1900 - 2017?
  • RQ2: Cấu trúc phân bố không gian - thời gian, độ sâu chấn tiêu và cơ cấu chấn tiêu của các chuỗi động đất mạnh tại Biển Đông phản ánh bản chất địa động lực gì của hệ thống đứt gãy sinh chấn?
  • RQ3: Mức độ rủi ro phát sinh độ lớn động đất cực đại ($M_{max}$) tại bồn trũng Biển Đông và đới Manila Trench đạt giá trị bao nhiêu theo các mốc thời gian xác định?
  • H1: Tích hợp hồi quy bình phương tối thiểu đa phân đoạn giữa các thang $M_s, m_b, M_L$ với $M_w$ kết hợp thuật toán cửa sổ không gian - thời gian sẽ tạo ra danh mục động đất độc lập đạt mức đại diện địa chấn tin cậy ($M_c \le 4.7$).
  • H2: Quá trình giải phóng ứng suất trong vùng nguồn của các động đất $M_w \ge 7.0$ có mối liên hệ nhân quả trực tiếp với cơ chế nén ép nghịch chờm và trượt bằng của các đứt gãy hoạt động, cho phép định vị và kiểm chứng các đứt gãy giả định.
  • H3: Mô hình phân bố cực trị tổng quát (GEV) cho phép dự báo chính xác giá trị $M_{max}$ đạt ngưỡng tiệm cận 8.7 trong chu kỳ 100 năm với xác suất vượt ngưỡng $P = 0.8$.

Khung lý thuyết của công trình tích hợp chặt chẽ giữa Lý thuyết Địa động lực thạch quyển và Kiến tạo mảng (Le Pichon, 1968; Hayes, 1994), Lý thuyết Năng lượng và Mômen Địa chấn (Kanamori, 1977), Quy luật Lặp lại Địa chấn Thực nghiệm Gutenberg - Richter (1944, 1956) và Lý thuyết Thống kê Cực trị Hiện đại (Pisarenko et al., 2008, 2010).

Nghiên cứu được triển khai trên quy mô không gian toàn diện: vùng thu thập số liệu mở rộng ($\varphi = 11.2^\circ\text{S} - 35.5^\circ\text{N}; \lambda = 92.5^\circ\text{E} - 132^\circ\text{E}$) và vùng nghiên cứu trọng tâm Biển Đông ($\varphi = 5^\circ\text{S} - 26^\circ\text{N}; \lambda = 100^\circ\text{E} - 127^\circ\text{E}$). Tập dữ liệu sơ cấp khổng lồ gồm 377.784 sự kiện địa chấn quan trắc thuần túy bằng máy đo ($3.0 \le M \le 9.1$) kéo dài suốt 118 năm (1900 - 2017), tổng hợp từ Trung tâm Địa chấn Quốc tế (ISC: 283.485 trận), Cục Khảo sát Địa chất Hoa Kỳ (USGS: 86.644 trận) và Hệ thống Cảnh báo Sớm Đa rủi ro Khu vực RIMES (7.655 trận).

Literature Review và Positioning

Lịch sử nghiên cứu địa chấn học khu vực Đông Nam Á và Biển Đông được phân định qua ba dòng học thuật chính:

  1. Dòng nghiên cứu thang đo năng lượng và quan hệ tương quan magnitude: Khởi xướng từ thang độ lớn địa phương $M_L$ (Richter, 1935), sóng mặt $M_s$ (Gutenberg, 1945a; Vanek et al., 1962), sóng khối $m_b$ (Gutenberg, 1945b; Abe, 1981) đến lý thuyết mômen địa chấn $M_0$ và magnitude mômen $M_w$ (Kanamori, 1977). Các nhà nghiên cứu quốc tế đã xây dựng phương trình chuyển đổi magnitude cục bộ như Papazachos et al. (1997) tại biển Aegean ($M_L = 0.58 M_s + 2.14$), Scordilis (2006) trên bình diện toàn cầu ($M_w = 0.85 m_b + 1.03$), hay Sankar Kumar Nath et al. (2017) tại Nam Á áp dụng hồi quy trực giao đa phân đoạn. Tại Việt Nam, Ngô Thị Lư (1999, 2003) và Lê Tử Sơn (2008) đã bước đầu thiết lập các hàm tương quan $M_{sh}-M_b$ và $M_L-M_D$ cho vùng Đông Dương nhưng quy mô mẫu còn rất hạn chế (dưới 120 sự kiện).
  2. Dòng nghiên cứu phân vùng phát sinh động đất và cơ cấu chấn tiêu: Nghiên cứu kiến tạo trẻ của Saurin (1965), Bùi Công Quế và nnk (1995, 2010, 2014), Phan Trọng Trịnh và nnk (2010, 2011, 2012) qua số liệu đo lặp GPS và địa chấn dầu khí; cùng các công trình phân vùng nguồn phát sinh sóng thần của Nguyễn Hồng Phương và nnk (2012, 2014), Xu (2014), Nguyễn Văn Lương và Cao Đình Triều (2014).
  3. Dòng nghiên cứu dự báo độ lớn động đất cực đại $M_{max}$: Điển hình là trường phái tham số hình học đứt gãy (Wells & Coppersmith, 1994; Stirling et al., 2002) và trường phái xác suất thống kê cổ điển dựa trên phân bố Gumbel I, Gumbel III (Gumbel, 1958; Nguyễn Kim Lạp, 1986; Nguyễn Hồng Phương, 1998; Trần Thị Mỹ Thành, 2002; Đặng Thanh Hải, 2003).
                      ┌────────────────────────────────────────┐
                      │  KHÔNG GIAN LÝ THUYẾT & THỰC TIỄN      │
                      │  ĐỊA CHẤN HỌC BIỂN ĐÔNG (1900 - 2017)  │
                      └──────────────────┬─────────────────────┘
                                         │
     ┌───────────────────────────────────┼───────────────────────────────────┐
     ▼                                   ▼                                   ▼
┌─────────────────────────┐ ┌─────────────────────────┐ ┌─────────────────────────┐
│ DÒNG 1: THANG ĐỘ LỚN    │ │ DÒNG 2: PHÂN VÙNG NGUỒN │ │ DÒNG 3: DỰ BÁO Mmax     │
│ & CHUẨN HÓA DỮ LIỆU     │ │ & ĐỊA ĐỘNG LỰC HỌC      │ │ & XÁC SUẤT CỰC TRỊ      │
├─────────────────────────┤ ├─────────────────────────┤ ├─────────────────────────┤
│• Richter (1935): ML     │ │• Bùi Công Quế (2010):   │ │• Wells & Coppersmith    │
│• Vanek et al.(1962): Ms │ │  37 đứt gãy hoạt động   │  │  (1994): Hình học rãnh │
│• Kanamori (1977): Mw    │ │• Phan Trọng Trịnh(2011):│ │• Nguyễn Hồng Phương     │
│• Sankar Nath et al.(2017│ │  GPS IGS05 biến dạng    │  │  (1998): Hợp lý cực đại│
│  Hồi quy đa phân đoạn   │ │• Xu (2014): Manila rãnh │ │• Gumbel I, III cổ điển  │
└────────────┬────────────┘ └────────────┬────────────┘ └────────────┬────────────┘
             │                           │                           │
             └───────────────────────────┼───────────────────────────┘
                                         │
                                         ▼
                     ┌───────────────────────────────────────┐
                     │ POSITIONING & RESEARCH GAPS CHÍNH:    │
                     │ 1. Trộn lẫn số liệu lịch sử & máy     │
                     │ 2. Bão hòa magnitude chưa chuẩn hóa Mw│
                     │ 3. Khuyết thiếu mốc thời gian dự báo  │
                     └───────────────────┬───────────────────┘
                                         │
                                         ▼
                     ┌───────────────────────────────────────┐
                     │ ĐÓNG GÓP TIÊN PHONG CỦA VŨ THỊ HOÃN:  │
                     │ • Đồng nhất hóa chuẩn Mw (1900-2017)  │
                     │ • Phân tách tiền/dư chấn không-thời   │
                     │ • Áp dụng GEV: Mmax = 8.7 (P=0.8/100n)│
                     └───────────────────────────────────────┘

Trong bức tranh tổng quan đó, tồn tại hai luồng tranh luận học thuật gay gắt:

  • Tranh luận 1: Phương pháp địa chất kiến tạo vs. Phương pháp xác suất thống kê. Phái kiến tạo (như Wells & Coppersmith, 1994) giả định mặt đứt gãy hình học cố định (chiều dài gấp đôi chiều rộng) để tính $M_{max}$. Tuy nhiên, tại môi trường biển sâu, dữ liệu thực địa ngầm khó phản ánh đầy đủ thông số hình học đứt gãy, dẫn đến sai số mô hình rất lớn. Ngược lại, phái xác suất (Nguyễn Hồng Phương, 1998) giải quyết bài toán qua lịch sử chấn động nhưng lại vấp phải nhược điểm nghiêm trọng khi áp dụng cho các vùng có dữ liệu thưa thớt (dưới 10 sự kiện/tiểu vùng), làm vi phạm điều kiện hội tụ thống kê.
  • Tranh luận 2: Vấn đề bão hòa thang đo và lựa chọn chu kỳ bước nhảy thời gian. Thang $M_L, m_b, M_s$ lần lượt bão hòa ở ngưỡng 6.8, 7.0 và 8.0, khiến việc đánh giá các siêu động đất bị đánh giá thấp đáng kể. Đồng thời, việc chọn bước nhảy 6 tháng (Nguyễn Kim Lạp, 1986) hay 2 năm (Đặng Thanh Hải, 2003) khi áp dụng hàm Gumbel hoàn toàn mang tính chủ quan, thiếu cơ sở toán lý vững chắc.

Công trình của Vũ Thị Hoãn đã định vị chính xác vị thế học thuật: thiết lập bước tiến đột phá so với các công trình quốc tế bằng cách giải quyết triệt để vấn đề bão hòa thang đo thông qua việc chuyển đổi toàn bộ 377.784 sự kiện về thang $M_w$, ứng dụng phương pháp cửa sổ không gian - thời gian để lọc sạch các cụm tiền chấn - dư chấn, và áp dụng mô hình phân bố cực trị tổng quát (GEV) của Pisarenko et al. (2008, 2010) cho toàn bộ bồn trũng Biển Đông. Kết quả này vượt trội so với nghiên cứu của Dietrich Stromeyer et al. (2004) tại Thụy Sĩ hay Sankar Kumar Nath et al. (2017) tại Nam Á nhờ khả năng đưa ra dự báo $M_{max}$ gắn liền với xác suất tin cậy và mốc thời gian cụ thể.

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

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

Nghiên cứu đóng góp trực tiếp vào việc mở rộng và kiểm chứng ba trụ cột lý thuyết địa chấn học hiện đại:

  • Mở rộng Lý thuyết Năng lượng và Mômen Địa chấn của Hiroo Kanamori (1977): Khẳng định tính ưu việt của thang $M_w$ trong việc phản ánh bản chất vật lý vùng nguồn thông qua tenxơ mômen địa chấn $M_0 = \mu D S$ (với $\mu$ là độ rắn chắc thạch quyển, $D$ là độ dịch trượt đứt gãy, $S$ là diện tích mặt đứt gãy). Nghiên cứu đã chứng minh rằng các hàm hồi quy phi tuyến đa phân đoạn giải quyết triệt để độ lệch hệ thống do hiệu ứng bão hòa sóng mặt ($M_s \ge 8.0$) và sóng khối ($m_b \ge 7.0$).
  • Bổ chỉnh Quy luật Phân bố Gutenberg - Richter (1944): Tái lập phương trình lặp lại tần số - độ lớn $\log_{10} N = a - b M_w$ trên tập dữ liệu Biển Đông đã khử cụm (declustered), chứng minh hệ số độ dốc địa chấn ($b$-value) phản ánh trực tiếp trạng thái ứng suất vi sai và cơ chế biến dạng vỏ Trái Đất trong khu vực.
  • Phát triển Lý thuyết Thống kê Cực trị Tiệm cận trong Khoa học Trái Đất: Kiểm chứng sự hội tụ của ba phân bố cực trị cổ điển (Gumbel, Fréchet, Weibull) về hàm phân bố cực trị tổng quát (GEV) do Pisarenko et al. (2008) đề xuất, xác lập mô hình toán học vững chắc cho phép xác định độ lớn cực đại $M_{max}$ không phụ thuộc vào các giả định hình học đứt gãy chủ quan.

Mô hình lý thuyết được lượng hóa thông qua hệ thống mệnh đề khoa học:

  • Proposition P1: Tồn tại mối quan hệ tương quan phi tuyến giữa magnitude sóng mặt ($M_s$), sóng khối ($m_b$), sóng địa phương ($M_L$) với magnitude mômen ($M_w$), trong đó dạng hàm bậc hai hoặc hàm phân đoạn bậc nhất tối thiểu hóa tổng bình phương sai số ($SSE$).
  • Proposition P2: Đám mây dư chấn của các trận động đất mạnh ($M_w \ge 7.0$) phát triển dọc theo phương cấu trúc của đứt gãy sinh chấn chính, đóng vai trò chỉ thị động học để nhận diện các đứt gãy hoạt động ngầm dưới đáy biển.
  • Proposition P3: Hàm phân bố cực trị tổng quát $G(x; \xi, \mu, \sigma)$ mô tả chính xác xác suất xuất hiện giá trị magnitude cực đại hàng năm tại đới hút chìm Manila và bồn trũng Biển Đông, với tham số hình dạng $\xi < 0$ chỉ thị giới hạn trên hữu hạn của năng lượng địa chấn tích lũy ($M_{max}$).

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

Khung phân tích của luận án tích hợp liên ngành ba hệ thống lý thuyết: (1) Lý thuyết Kiến tạo mảng và Đới hút chìm sâu (Subduction Tectonics), (2) Lý thuyết Tương quan Thống kê Địa chấn học thực nghiệm, và (3) Lý thuyết Phân bố Cực trị Tổng quát (GEV).

Quy trình phân tích triển khai tuần tự theo chuỗi thuật toán nghiêm ngặt:

$$\text{Tập số liệu ISC/USGS/RIMES} \xrightarrow{\text{Đối chiếu ISC}} \text{Hiệu chỉnh tọa độ, thời gian} \xrightarrow{\text{Hồi quy OLS}} \text{Quy chuẩn } M_w$$

$$\xrightarrow{\text{Cửa sổ Không - Thời gian}} \text{Tách TC - DC} \xrightarrow{\text{GEV Modeling}} M_{max} \text{ (Kèm } P \text{ và } T\text{)}$$

Điều kiện biên (boundary conditions) của mô hình phân tích được xác lập rõ ràng: giới hạn không gian bao trùm $\varphi = 5^\circ\text{S} - 26^\circ\text{N}$ và $\lambda = 100^\circ\text{E} - 127^\circ\text{E}$; ngưỡng lọc magnitude $M_w \ge 3.0$; và giả định tính độc lập ngẫu nhiên hoàn toàn (Poisson process) của các sự kiện địa chấn sau khi đã loại bỏ triệt để các chuỗi tiền chấn - dư chấn.

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

Thiết kế nghiên cứu

Công trình được xây dựng trên nền tảng triết lý thực chứng (positivism) với phương pháp tiếp cận định lượng chuyên sâu (rigorous quantitative methodology). Thiết kế nghiên cứu đa tầng (multi-level design) được cấu trúc qua ba giai đoạn:

  1. Giai đoạn 1: Chuẩn hóa và làm sạch tập dữ liệu không gian - thời gian lớn (Big Spatio-Temporal Seismic Data).
  2. Giai đoạn 2: Khử cụm địa chấn (seismic declustering) bằng thuật toán cửa sổ động học để thiết lập danh mục biến cố độc lập.
  3. Giai đoạn 3: Phân tích cơ cấu chấn tiêu (focal mechanism) 3D kết hợp ước lượng tham số cực trị GEV phi tuyến.

Quy trình nghiên cứu rigorous

Tác giả đã áp dụng chiến lược thu thập và xử lý dữ liệu qua các giao thức nghiêm ngặt:

  • Giao thức đối chiếu và khử trùng lặp (Cross-checking & De-duplication Protocol): Khi có sự sai lệch thông số giữa các trung tâm quan trắc quốc tế, hệ thống ưu tiên tuyệt đối dữ liệu của Trung tâm Địa chấn Quốc tế (ISC) – cơ quan chuẩn mực toàn cầu về hiệu chỉnh vị trí chấn tâm và độ sâu chấn tiêu.
  • Giao thức lọc chuỗi tiền chấn - dư chấn (Declustering Protocol): Ứng dụng phương pháp cửa sổ không gian - thời gian (Space-Time Window Method) với bán kính không gian $L(M)$ và khoảng thời gian $T(M)$ mở rộng theo hàm phi tuyến của magnitude kích động chính (mainshock). Kỹ thuật này phân tách toàn bộ dữ liệu thành hai tập chuyên biệt: Danh mục các trận động đất độc lập (phục vụ mô hình thống kê Poisson và tính $M_{max}$) và Danh mục các nhóm tiền chấn - dư chấn (phục vụ giải đoán động học đứt gãy).
  • Giao thức kiểm định độ tin cậy và mức đại diện (Completeness & Reliability): Độ hoàn thiện của danh mục ($M_c$) được kiểm tra thông qua biểu đồ độ lặp lại Gutenberg - Richter, xác định giá trị ngưỡng $M_w = 4.7$, đảm bảo dữ liệu từ ngưỡng này trở lên không bị bỏ sót bởi mạng lưới trạm quan trắc.

Data và phân tích

Phép phân tích hồi quy bình phương tối thiểu (Ordinary Least Squares - OLS) được thực thi để xây dựng các phương trình chuyển đổi magnitude thực nghiệm. Tác giả đã thử nghiệm cả dạng hàm tuyến tính bậc nhất và phi tuyến bậc hai, phân tích tổng bình phương sai số ($SSE$) để chọn ra mô hình tối ưu:

$$\text{Hàm tương quan } M_w - M_s: \quad M_w = f(M_s)$$

$$\text{Hàm tương quan } M_w - M_b: \quad M_w = f(M_b)$$

$$\text{Hàm tương quan } M_w - M_L: \quad M_w = f(M_L) \quad (\text{phân đoạn theo ngưỡng } M_L)$$

Mô hình phân bố cực trị tổng quát (GEV) thống nhất hàm phân bố xác suất tích lũy của giá trị cực đại hàng năm $X = \max(M_w)$ theo dạng:

$$G(x; \xi, \mu, \sigma) = \exp \left{ - \left[ 1 + \xi \left( \frac{x - \mu}{\sigma} \right) \right]^{-1/\xi} \right}$$

Trong đó:

  • $\mu \in \mathbb{R}$ là tham số vị trí (location parameter),
  • $\sigma > 0$ là tham số quy mô (scale parameter),
  • $\xi \in \mathbb{R}$ là tham số hình dạng (shape parameter). Khi $\xi < 0$, hàm phân bố thuộc dạng Weibull (Type III), chỉ thị sự tồn tại của chặn trên hữu hạn $M_{max} = \mu - \sigma/\xi$.

Toàn bộ các thuật toán tính toán thống kê, hồi quy phi tuyến, xử lý cửa sổ không - thời gian và vẽ mặt cắt 3D chấn tiêu được thực hiện tự động hóa bằng hệ thống phần mềm địa chấn chuyên dụng kết hợp các script tính toán toán học.

Thông số nghiên cứu Giá trị định lượng Nguồn kiểm chứng / Ý nghĩa
Tổng số sự kiện thu thập 377.784 trận động đất Tổng hợp ISC (283.485), USGS (86.644), RIMES (7.655)
Khoảng thời gian quan trắc 1900 – 2017 (118 năm) 100% dữ liệu ghi nhận bằng máy địa chấn
Tọa độ vùng nghiên cứu $\varphi = 5^\circ\text{S} - 26^\circ\text{N}$; $\lambda = 100^\circ\text{E} - 127^\circ\text{E}$ Bao trùm toàn bộ cấu trúc kiến tạo Biển Đông
Ngưỡng đại diện ($M_c$) $M_w = 4.7$ Xác định qua biểu đồ Gutenberg - Richter
Tỷ lệ độ sâu chấn tiêu $> 80%$ xảy ra ở độ sâu $h < 75\text{ km}$ Năng lượng tích tụ tập trung trong lớp vỏ Trái Đất
Độ lớn cực đại ($M_{max}$) $M_w = 8.7$ Tính theo GEV trong 100 năm với xác suất $P = 0.8$
Chu kỳ lặp lại $M_w \ge 7.5$ 3 – 5 năm / lần Đánh giá nguy cơ địa chấn quy mô toàn vùng

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

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

Công trình đã mang lại 5 phát hiện khoa học mang tính đột phá cho ngành Địa chấn học và Địa động lực biển:

  1. Xác lập danh mục động đất đồng nhất chuẩn $M_w$ đầu tiên cho Biển Đông (1900 - 2017): Khử bỏ hoàn toàn sai số ghép nối giữa các cơ quan quan sát quốc tế, chuẩn hóa 377.784 sự kiện về thang $M_w$ ($3.0 \le M_w \le 8.5$), đạt ngưỡng hoàn thiện $M_c = 4.7$.

  2. Làm sáng tỏ trường phân bố không gian và độ sâu chấn tiêu: Đại bộ phận chấn tâm động đất Biển Đông tập trung dày đặc dọc theo ranh giới mảng phía Đông (đới rãnh sâu Manila, Đài Loan, Philippines). Đáng chú ý, hơn 80% tổng số sự kiện tập trung ở độ sâu nông ($h < 75\text{ km}$), chứng minh hoạt động giải phóng ứng suất diễn ra chủ yếu trong tầng vỏ giòn (brittle crust) của Trái Đất – điều kiện tối ưu để truyền sóng xung kích gây chấn động mạnh và sóng thần.

  3. Định lượng chu kỳ tái lập địa chấn cực mạnh: Phân tích chuỗi thời gian chỉ ra rằng chu kỳ lặp lại của các trận động đất có độ lớn $M_w \ge 7.5$ trong toàn khu vực Biển Đông và lân cận là 3 - 5 năm/lần, khẳng định đây là khu vực có hoạt tính địa chấn vô cùng tích cực chứ không hề ở trạng thái "ngủ yên".

  4. Khám phá và kiểm chứng hệ thống đứt gãy hoạt động ngầm qua giải đoán 6 siêu động đất: Phân tích cơ cấu chấn tiêu và mặt cắt 3D không gian - thời gian của đám mây dư chấn đi kèm 6 trận động đất mạnh gồm:

    • Trận số 1 (21/01/2007, Molucca, $M_w 7.5$),
    • Trận số 2 (12/09/2007, Nam Sumatra, $M_w 8.4$),
    • Trận số 3 (11/02/2009, Talaud, $M_w 7.2$),
    • Trận số 4 (31/08/2012, Sulangan, $M_w 7.6$),
    • Trận số 5 (15/10/2013, Balilihan, $M_w 7.1$),
    • Trận số 6 (15/11/2014, Balilihan, $M_w 7.1$).

    Kết quả nghiên cứu đã chứng minh sự tồn tại thực tế của một đứt gãy hoạt động ngầm trước đây chỉ mang tính giả định, đồng thời phác thảo chính xác phương vị và góc cắm của 3 đoạn đứt gãy sinh chấn giả định mới trong khu vực.

  5. Xác lập giá trị động đất cực đại $M_{max} = 8.7$ theo mô hình GEV: Kết quả tính toán bằng hàm phân bố cực trị tổng quát cho thấy giá trị $M_{max}$ của khu vực Biển Đông đạt ngưỡng $M_w = 8.7$ trong khoảng thời gian 100 năm tính từ năm 2017 với mức xác suất tin cậy $P = 0.8$. Kết quả này tương đồng với các tính toán của Nguyễn Hồng Phương và nnk (2012) đối với vùng Bắc Manila Trench ($M_{max} = 8.93$, chu kỳ 2658 năm), nhưng vượt trội ở chỗ xác lập được khung thời gian cụ thể (100 năm) và loại bỏ được tính bất định của các biến hình học trung gian.

                      ┌────────────────────────────────────────┐
                      │    HỆ THỐNG PHÁT HIỆN ĐỘT PHÁ          │
                      │  LUẬN ÁN TIẾN SĨ - VŨ THỊ HOÃN (2018) │
                      └──────────────────┬─────────────────────┘
                                         │
     ┌───────────────────────────────────┼───────────────────────────────────┐
     ▼                                   ▼                                   ▼
┌─────────────────────────┐ ┌─────────────────────────┐ ┌─────────────────────────┐
│ 1. CHUẨN HÓA DANH MỤC   │ │ 2. BẢN ĐỒ KHÔNG - THỜI  │ │ 3. DỰ BÁO TOÁN HỌC GEV  │
├─────────────────────────┤ ├─────────────────────────┤ ├─────────────────────────┤
│• 377.784 sự kiện (100%  │ │• >80% chấn tiêu h < 75km│ │• Phân bố Weibull (ξ < 0)│
│  quan trắc bằng máy)    │ │  (tầng vỏ Trái Đất)     │ │• Mmax = 8.7             │
│• Đồng nhất chuẩn Mw     │ │• Chu kỳ lặp Mw ≥ 7.5:   │ │• Khung thời gian: 100   │
│  (Kanamori, 1977)       │ │  3 - 5 năm / lần        │ │  năm kể từ 2017         │
│• Mức đại diện Mc = 4.7  │ │• Giải mã 6 siêu động đất│ │• Xác suất tin cậy:      │
│• Phân tách sạch TC - DC │ │• Minh chứng 1 đứt gãy   │ │  P = 0.8 (80%)          │
│  qua Space-Time Window  │ │  thật & 3 đoạn giả định │ │• Khắc phục Gumbel cũ    │
└─────────────────────────┘ └─────────────────────────┘ └─────────────────────────┘

Implications đa chiều

  • Về mặt lý thuyết (Theoretical Implications): Khẳng định tính hiệu quả vượt trội của mô hình GEV trong phân tích cực trị địa vật lý; cung cấp bức tranh hoàn chỉnh về động học mảng và tương tác giữa đới hút chìm Manila với cấu trúc giãn đáy Biển Đông.
  • Về mặt phương pháp luận (Methodological Innovations): Thiết lập quy trình chuẩn mực (pipeline) từ thu thập số liệu đa nguồn, xử lý lọc nhiễu, hồi quy quy chuẩn magnitude, khử cụm không - thời gian đến dự báo xác suất cực trị, có khả năng chuyển giao áp dụng trực tiếp cho các vùng biển phức tạp khác trên thế giới.
  • Về mặt ứng dụng thực tiễn và chính sách (Practical & Policy Applications):
    • Tham số đầu vào thiết yếu cho mô hình sóng thần: Giá trị $M_{max} = 8.7$ tại rãnh Manila là thông số cốt lõi để các cơ quan khí tượng thủy văn và cứu nạn cứu hộ quốc gia cập nhật kịch bản ngập lụt sóng thần cho dải ven biển miền Trung từ Đà Nẵng đến Bình Thuận.
    • Quy chuẩn thiết kế công trình biển: Cung cấp cơ sở khoa học định lượng để Bộ Xây dựng, Bộ Giao thông Vận tải và Tập đoàn Dầu khí Việt Nam tái tính toán tải trọng kháng chấn cho các công trình giàn khoan dầu khí ngoài khơi, cảng biển nước sâu, tuyến cáp quang ngầm xuyên biển và các đô thị ven biển.

Limitations và Future Research

Nghiên cứu thẳng thắn thừa nhận các giới hạn nội tại mang tính khách quan:

  • Độ phân giải của mạng lưới trạm quan trắc: Các trạm ghi nhận địa chấn phân bố không đồng đều, tập trung chủ yếu trên đất liền và các hải đảo lớn, khiến độ chính xác trong việc xác định tọa độ chấn tâm và độ sâu chấn tiêu ($h$) ngoài biển khơi còn tồn tại dung sai nhất định.
  • Giới hạn thời gian quan trắc máy (118 năm): So với chu kỳ tích lũy và giải phóng năng lượng kiến tạo kéo dài hàng nghìn năm của các đới hút chìm lớn, chuỗi số liệu 1900 - 2017 vẫn là một lát cắt thời gian tương đối ngắn.
  • Phương pháp hồi quy đơn giản: Luận án sử dụng phương pháp bình phương tối thiểu thông thường (OLS) thay vì hồi quy trực giao chuẩn (Orthogonal Standard Regression) – vốn được một số tác giả (Deniz & Yucemen, 2010; Das et al., 2011) đánh giá là có khả năng giảm thiểu sai số đồng thời trên cả hai trục biến số.

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

  1. Tích hợp chuỗi số liệu quan trắc đáy biển (Ocean Bottom Seismometer - OBS) nhằm nâng cao độ chính xác định vị chấn tiêu các vi địa chấn ngầm.
  2. Ứng dụng kỹ thuật hồi quy trực giao phi tuyến và thuật toán máy học (Machine Learning) để tối ưu hóa việc phân tách cụm tiền chấn - dư chấn.
  3. Liên kết mô hình địa chấn GEV với dữ liệu biến dạng viễn thám vệ tinh InSAR và mạng lưới GNSS liên tục độ chính xác cao.
  4. Mô phỏng động lực học vết nứt đứt gãy 3D kết hợp lan truyền sóng thần thời gian thực (Real-time Tsunami Inundation Modeling) dựa trên kịch bản $M_{max} = 8.7$.

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

  • Tác động học thuật (Academic Impact): Công trình trở thành tài liệu tham khảo chuẩn mực về địa chấn học Biển Đông, đóng góp bộ cơ sở dữ liệu đồng nhất $M_w$ chất lượng cao phục vụ các nghiên cứu phân vùng rủi ro thiên tai và vật lý thạch quyển.
  • Chuyển đổi ngành công nghiệp (Industry Transformation): Cung cấp khung tham chiếu an toàn địa chấn cho các nhà thầu khai thác dầu khí quốc tế (Vietsovpetro, PVN, ENI) và các đơn vị thiết kế kết cấu công trình biển, giúp tối ưu hóa chi phí gia cố kháng chấn hàng triệu USD nhưng vẫn đảm bảo tuyệt đối an toàn vận hành.
  • Định hình chính sách công (Policy Influence): Là luận cứ khoa học then chốt để Ban Chỉ đạo Quốc gia về Phòng chống Thiên tai xây dựng Bản đồ phân vùng cảnh báo nguy hiểm sóng thần quốc gia và quy hoạch không gian biển bền vững.
  • Ý nghĩa quốc tế (Global Relevance): Đóng góp vào cơ sở dữ liệu giám sát địa chấn Thái Bình Dương của RIMES và UNESCO/IOC, củng cố hợp tác an ninh phi truyền thống giữa các quốc gia ven Biển Đông.

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

  • Nghiên cứu sinh và Giới học thuật Địa vật lý / Địa chất: Tiếp cận danh mục động đất 377.784 sự kiện đã được làm sạch và chuẩn hóa $M_w$, tiết kiệm hàng nghìn giờ xử lý dữ liệu thô và kế thừa phương pháp luận GEV tiên tiến.
  • Doanh nghiệp Dầu khí và Năng lượng Tái tạo Biển: Sử dụng bộ tham số $M_{max} = 8.7$ và gia tốc nền thiết kế để đảm bảo tuổi thọ công trình ngoài khơi, phòng ngừa thảm họa tràn dầu và sụp đổ giàn khoan.
  • Cơ quan Hoạch định Chính sách và Cảnh báo Thiên tai: Nhận diện chính xác nguy cơ từ đới hút chìm Manila, hoàn thiện hệ thống trạm phao cảnh báo sóng thần và diễn tập sơ tán dân cư ven biển có căn cứ khoa học.
  • Cộng đồng Cư dân Duyên hải: Được bảo vệ an toàn tính mạng và tài sản thông qua các tiêu chuẩn xây dựng đô thị và hệ thống cảnh báo sớm rủi ro thiên tai đị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 mở rộng Lý thuyết Năng lượng và Mômen Địa chấn của Hiroo Kanamori (1977) kết hợp với Lý thuyết Phân bố Cực trị Tổng quát (GEV) của Pisarenko et al. (2008) trên một bồn trũng biển nửa kín phức tạp. Luận án đã chứng minh trên thực nghiệm quy mô lớn rằng thang độ lớn mômen ($M_w$) là thang đo duy nhất phản ánh chính xác cơ chế vật lý vùng nguồn mà không bị bão hòa, đồng thời xác lập mô hình cực trị Weibull ($\xi < 0$) có chặn trên hữu hạn cho động đất Biển Đông.

2. Tính đổi mới trong phương pháp luận so với các nghiên cứu quốc tế trước đây?

So với nghiên cứu của Papazachos et al. (1997) tại biển Aegean hay Sankar Kumar Nath et al. (2017) tại Nam Á, luận án có hai điểm đột phá: (1) Chuẩn hóa thuần túy 100% dữ liệu quan trắc máy (loại bỏ hoàn toàn sai số từ dữ liệu lịch sử ước lượng), và (2) Tích hợp phương pháp cửa sổ không gian - thời gian phi tuyến để phân tách triệt để tiền chấn/dư chấn trước khi áp dụng mô hình GEV, giải quyết triệt để vấn đề phụ thuộc thống kê của chuỗi biến cố địa chấn.

3. Phát hiện nào gây bất ngờ nhất và bằng chứng định lượng đi kèm?

Phát hiện bất ngờ nhất là việc minh chứng tính hoạt động thực tế của một đứt gãy kiến tạo ngầm trước đây chỉ tồn tại ở dạng giả định lý thuyết, đồng thời định vị thêm 3 đoạn đứt gãy giả định mới. Bằng chứng định lượng được xác lập qua việc phân tích phân bố không gian 3D và cơ cấu chấn tiêu của các đám mây dư chấn từ 6 trận động đất cực mạnh ($M_w \ge 7.1$), đặc biệt là trận Nam Sumatra 2007 ($M_w 8.4$) và Sulangan 2012 ($M_w 7.6$).

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

Có. Luận án mô tả chi tiết quy trình 4 bước hoàn chỉnh: từ tọa độ thu thập dữ liệu ($\varphi = 11.2^\circ\text{S} - 35.5^\circ\text{N}; \lambda = 92.5^\circ\text{E} - 132^\circ\text{E}$), nguồn cơ sở dữ liệu (ISC, USGS, RIMES), các phương trình hồi quy quy đổi $M_w = f(M_s, m_b, M_L)$, bảng tham số cửa sổ không gian - thời gian $L(M)$ và $T(M)$, đến thuật toán ước lượng cực trị GEV, cho phép các nhà khoa học độc lập tái lập và kiểm chứng toàn bộ kết quả.

5. Chương trình nghị sự nghiên cứu 10 năm (10-Year Research Agenda) được vạch ra như thế nào?

Chương trình nghiên cứu định hướng: (1) Thiết lập mạng lưới máy ghi địa chấn đáy biển (OBS) băng rộng dọc đới Manila Trench; (2) Tích hợp mô hình dịch chuyển vỏ Trái Đất từ mạng lưới GPS liên tục IGS05 với trường tích lũy ứng suất Coulomb 3D; (3) Nâng cấp mô hình GEV thành mô hình cực trị không dừng (Non-stationary GEV) tích hợp yếu tố biến đổi thời gian; (4) Xây dựng hệ thống cảnh báo sớm sóng thần thời gian thực kết nối vệ tinh phục vụ phòng thủ bờ biển.

Kết luận

  1. Thành lập thành công danh mục 377.784 sự kiện địa chấn (1900 - 2017) thuần túy quan trắc bằng máy, đồng nhất hoàn toàn về thang độ lớn mômen $M_w$ ($3.0 \le M_w \le 8.5$), đạt ngưỡng đại diện địa chấn $M_c = 4.7$.
  2. Phân tách triệt để các chuỗi tiền chấn - dư chấn bằng phương pháp cửa sổ không - thời gian, tạo lập đồng thời tập dữ liệu độc lập cho thống kê và tập dữ liệu cụm cho nghiên cứu động học đứt gãy.
  3. Xác định quy luật phân bố độ sâu chấn tiêu với hơn 80% động đất xảy ra trong lớp vỏ Trái Đất ($h < 75\text{ km}$) và chu kỳ tái lặp động đất mạnh $M_w \ge 7.5$ là 3 - 5 năm/lần.
  4. Minh chứng sự tồn tại thực tế của 1 đứt gãy hoạt động ngầm và xác lập cấu trúc hình học của 3 đoạn đứt gãy sinh chấn mới thông qua giải đoán cơ cấu chấn tiêu 3D của 6 siêu động đất.
  5. Ứng dụng xuất sắc mô hình phân bố cực trị tổng quát (GEV), xác định độ lớn động đất cực đại khu vực Biển Đông đạt $M_{max} = 8.7$ trong thời hạn 100 năm kể từ năm 2017 với độ tin cậy xác suất $P = 0.8$, khắc phục triệt để nhược điểm của các phương pháp xác suất truyền thống.
  6. Mở ra 3 hướng nghiên cứu mũi nhọn mới cho khoa học Trái Đất Việt Nam: địa chấn học đáy biển sâu OBS, mô phỏng liên kết tai biến động đất - sóng thần thời gian thực, và đánh giá ứng suất kiến tạo hiện đại dựa trên đo lặp GNSS/InSAR độ phân giải cao.