Tổng quan nghiên cứu

Sự gia tăng tần suất các hiểm họa địa chấn tại khu vực Đông Nam Á, tiêu biểu như thảm họa động đất sóng thần Sumatra năm 2004 hay chuỗi động đất kích thích liên tiếp tại khu vực hồ chứa thủy điện Sông Tranh 2 (tỉnh Quảng Nam), đã đặt ra yêu cầu cấp thiết về việc nâng cao năng lực giám sát và đánh giá tai biến địa chất. Trước đây, mạng lưới quan trắc địa chấn Việt Nam chỉ có khoảng 3 đến 5 trạm ghi cơ học trên giấy ảnh, sau đó đã được đầu tư nâng cấp lên gần 30 trạm ghi số hiện đại trên toàn quốc, riêng khu vực công trình thủy điện Sông Tranh 2 được thiết lập bổ sung mạng lưới quan trắc cục bộ với 15 trạm địa chấn. Dù vậy, phần lớn các mô hình cấu trúc vận tốc vỏ Trái đất trước năm 2014 chỉ tập trung vào khu vực miền Bắc, trong khi khu vực miền Trung với chiều dài bờ biển lớn và cấu trúc kiến tạo phân dị phức tạp lại chưa có một mô hình chuẩn riêng biệt.

Vấn đề cốt lõi của nghiên cứu là xây dựng mô hình vận tốc truyền sóng địa chấn và hệ thống tốc đồ thời gian truyền sóng (biểu đồ thời khoảng) chuyên biệt cho vỏ Trái đất khu vực miền Trung Việt Nam. Mục tiêu cụ thể là tính toán các đường cong thời gian truyền sóng cho các pha sóng dọc và sóng ngang, làm chuẩn số liệu phục vụ quy toán và xác định chính xác các tham số cơ bản của chấn tiêu động đất (tọa độ chấn tâm, độ sâu chấn tiêu, thời điểm phát sinh chấn động). Nghiên cứu được triển khai trong phạm vi không gian bao phủ toàn bộ miền Trung và các vùng lân cận với bán kính chấn tâm mở rộng đến 800 km, dựa trên chuỗi số liệu của 275 sự kiện động đất ghi nhận từ các trạm trong nước và quốc tế. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn khi giúp giảm thiểu sai số định vị chấn tiêu xuống mức tối đa, cung cấp luận cứ khoa học chính xác cho công tác phân vùng động đất, dự báo hiểm họa và thiết kế kháng chấn các công trình trọng điểm quốc gia.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng lý thuyết truyền sóng đàn hồi trong môi trường địa chất phân lớp dị hướng và mô hình hóa toán học quá trình lan truyền sóng địa chấn. Các khái niệm nền tảng bao gồm: các pha sóng thể tích (sóng dọc P, sóng ngang S), sóng khối truyền trực tiếp trong tầng trên vỏ Trái đất (Pg, Sg), sóng phản xạ từ mặt ranh giới Conrad phân chia lớp Granit - Bazan (P*, S*), và sóng khúc xạ đầu truyền dọc theo mặt ranh giới Moho phân chia vỏ Trái đất với manti trên (Pn, Sn).

Khung lý thuyết vận dụng mô hình vận tốc chuẩn quốc tế IASPEI91 và mô hình vận tốc trung bình toàn cầu của Jeffreys-Bullen (1940, 1958) làm hệ quy chiếu ban đầu để tính toán sai số thời gian truyền sóng lý thuyết. Đồng thời, nghiên cứu kế thừa các đánh giá về tỷ số vận tốc sóng dọc trên sóng ngang ($K = V_P/V_S$) trong môi trường vỏ Trái đất Việt Nam (dao động ổn định từ 1,67 đến 1,70 theo đồ thị Vadati và phân tích sóng trao đổi PS). Mối quan hệ động học giữa thời gian truyền sóng $T$, khoảng cách chấn tâm $X$, vận tốc truyền sóng $V$ và độ sâu chấn tiêu $H$ được biểu diễn tổng quát qua phương trình thời khoảng lý thuyết $T = T(V, H, X)$.

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

Nguồn dữ liệu của nghiên cứu bao gồm tập hợp các cặp giá trị khoảng cách chấn tâm và thời gian truyền sóng ($X_i, T_i$) trích xuất từ 275 sự kiện động đất riêng biệt. Cỡ mẫu này được chọn lọc từ mạng lưới trạm địa chấn quốc gia của Viện Vật lý Địa cầu kết hợp cùng các trạm quan sát địa chấn quốc tế trong khu vực Đông Nam Á nhằm bù đắp sự thưa thớt của các trạm đo nội địa tại miền Trung. Phương pháp chọn mẫu là chọn mẫu có chủ đích kết hợp tiêu chuẩn lọc dữ liệu khắt khe: chỉ lựa chọn các trận động đất có ghi nhận rõ ràng các pha sóng (Pn, P, Pg, Sn, S, Sg), khoảng cách chấn tâm $\Delta \le 800$ km và có đầy đủ tham số chấn tiêu tin cậy.

Lý do lựa chọn phương pháp phân tích: Tác giả áp dụng chu trình phân tích gần đúng liên tiếp theo phương pháp của Kondorskaia kết hợp giải thuật tối ưu hóa hàm mục tiêu sai số thời gian:

$$F = \sum_{i=1}^n (T_i - T_{pi}')^2$$

Hàm mục tiêu được cực tiểu hóa thông qua phương pháp Monte-Carlo và phương pháp bình phương cực tiểu với tiêu chuẩn kiểm định thống kê $\chi^2$. Phương pháp này được ưu tiên lựa chọn vì sở hữu mối quan hệ ngược ưu việt: trong quá trình hiệu chỉnh tọa độ chấn tiêu sẽ đồng thời chuẩn hóa thời gian truyền sóng, cho phép làm sáng tỏ các bất đồng nhất cấu trúc vỏ Trái đất ngay cả khi độ sâu chấn tiêu chưa được xác định trước. Timeline nghiên cứu được thực hiện đồng bộ trong giai đoạn 2012-2014, thuộc khuôn khổ Nghị định thư hợp tác khoa học công nghệ cấp Nhà nước giữa Viện Hàn lâm Khoa học và Công nghệ Việt Nam và Viện Hàn lâm Khoa học Liên bang Nga.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Nghiên cứu đã xây dựng thành công bộ 6 phương trình tốc đồ thời gian truyền sóng trung bình cho các pha sóng địa chấn tại khu vực miền Trung Việt Nam, đạt độ tương thích thống kê rất cao:

  • Sóng dọc P: $T_P = 0,125X + 3,143$ với hệ số tương quan $R = 0,999$
  • Sóng trực tiếp Pg: $T_{Pg} = 0,168X + 0,498$ với hệ số tương quan $R = 0,999$
  • Sóng khúc xạ Pn: $T_{Pn} = 0,124X + 4,412$ với hệ số tương quan $R = 0,999$
  • Sóng ngang S: $T_S = 0,219X + 5,800$ với hệ số tương quan $R = 0,997$
  • Sóng ngang trực tiếp Sg: $T_{Sg} = 0,288X + 0,501$ với hệ số tương quan $R = 0,999$
  • Sóng ngang khúc xạ Sn: $T_{Sn} = 0,214X + 7,533$ với hệ số tương quan $R = 1,000$

Hiệu quả lọc nhiễu qua quy trình 3 bước lặp đã cải thiện vượt bậc chất lượng số liệu thực nghiệm. Đối với pha sóng Pn, phương trình ban đầu trước xử lý có dạng $T_{PN} = 0,126X + 4,363$ với hệ số tương quan $R_1 = 0,982$. Sau 3 vòng lọc bỏ các điểm phân tán mạnh, hệ số tương quan đã tăng lên $R_2 = 0,999$, tương ứng mức tăng độ tin cậy thêm khoảng 1,73%. Phân tích độ lệch thời gian $\Delta T = T' - T$ cho thấy mức sai số từ vài chục giây ban đầu (có điểm cá biệt lệch tới $-54$ giây) đã được thu hẹp về phạm vi dưới 0,5 giây, giúp đồ thị hàm số đạt độ trơn tối ưu.

Mô hình lát cắt vận tốc 1D xây dựng cho khu vực miền Trung xác lập rõ các ranh giới thạch học chính: tầng trầm tích mặt có vận tốc sóng P khoảng 4,1 đến 4,6 km/s, tầng Granit trung bình đạt khoảng 6,0 đến 6,2 km/s (bề dày khoảng 14 km), tầng Bazan đạt 6,6 đến 7,0 km/s và độ sâu mặt Moho dao động trong khoảng 28 đến 31 km.

Thảo luận kết quả

Quy luật vật lý thể hiện rõ nét trên đồ thị thời khoảng khi toàn bộ nhánh tốc đồ của sóng dọc P luôn nằm thấp hơn nhánh tốc đồ sóng ngang S trên cùng trục khoảng cách chấn tâm $X$. Điều này chứng minh tính nhất quán cơ học: vận tốc truyền sóng dọc luôn lớn hơn sóng ngang ($V_P > V_S$), do đó thời gian lan truyền sóng dọc đến trạm ghi luôn ngắn hơn ($T_P < T_S$).

Khi đối sánh với các mô hình lát cắt vận tốc ở miền Bắc Việt Nam của các tác giả tiền nhiệm (như lát cắt 3 lớp và 4 lớp của Ngô Thị Lư với 196 trận động đất Tây Bắc, lát cắt tuyến Thái Nguyên – Hòa Bình của Đinh Văn Toàn), mô hình miền Trung thể hiện sự tương đồng về cấu trúc phân tầng nhưng có sự phân dị rõ rệt về chiều dày tầng trầm tích và độ sâu móng kết tinh do ảnh hưởng của các đới đứt gãy ven biển và cơ chế mở rộng bồn trũng Biển Đông.

Để trực quan hóa các phát hiện trên, toàn bộ hệ dữ liệu có thể được biểu diễn một cách khoa học qua:

  • Biểu đồ phân tán hai chiều tương quan giữa khoảng cách $X$ (km) và thời gian truyền $T$ (giây) cho từng pha sóng, thể hiện rõ các điểm loại trừ và đường hồi quy tuyến tính sau 3 lần lọc.
  • Bảng ma trận so sánh các thông số vận tốc lớp ($V_i$), chiều dày tầng ($h_i$) và độ sâu ranh giới ($H$) giữa mô hình thực nghiệm miền Trung với mô hình chuẩn IASPEI91 và mô hình miền Bắc.
  • Đồ thị lát cắt vận tốc 1D dạng bậc thang thể hiện sự biến thiên của vận tốc $V_P$ và $V_S$ theo chiều sâu thẳng đứng từ mặt đất đến độ sâu 40 km.

Đề xuất và khuyến nghị

  1. Tích hợp trực tiếp bộ phương trình tốc đồ chuẩn miền Trung vào phần mềm định vị chấn tiêu tự động: Viện Vật lý Địa cầu cần cập nhật ngay 6 phương trình tốc đồ trung bình ($P, Pg, Pn, S, Sg, Sn$) vào hệ thống phần mềm xử lý dữ liệu động đất quốc gia. Mục tiêu là rút ngắn thời gian tính toán thông số chấn tiêu xuống dưới 5 phút sau khi ghi nhận sự kiện và giảm sai số tọa độ chấn tâm xuống dưới 2 km, hoàn thành trong giai đoạn 2026-2027.

  2. Mở rộng và hiện đại hóa mạng lưới trạm quan trắc địa chấn khu vực miền Trung: Bộ Khoa học và Công nghệ phối hợp cùng Tập đoàn Điện lực Việt Nam nâng cấp số lượng trạm đo từ 15 trạm hiện hữu lên tối thiểu 25 trạm ghi số băng rộng dải động cao (24-bit). Ưu tiên lắp đặt tập trung tại các bậc thang thủy điện dọc sông Tranh, sông Bung và các đới đứt gãy Rào Quán - A Lưới trước năm 2028 nhằm thu thập chuỗi số liệu thời gian thực có độ phân giải cao.

  3. Thực hiện các tuyến thăm dò địa chấn sâu 2D và 3D ven biển: Đề xuất các viện chuyên ngành triển khai ít nhất 2 tuyến đo địa chấn dò sâu (Deep Seismic Sounding) cắt ngang qua đới kiến tạo miền Trung ra thềm lục địa nhằm chính xác hóa bề mặt phản xạ Moho và Conrad, hướng tới nâng độ chính xác của lát cắt vận tốc địa phương lên trên 95% trước năm 2030.

  4. Ứng dụng mô hình vận tốc trong đánh giá rủi ro động đất kích thích và vi phân vùng địa chấn: Các cơ quan quản lý an toàn đập và quy hoạch đô thị cần áp dụng mô hình vận tốc này để tính toán lại phổ gia tốc nền cực đại (PGA) cho các công trình trọng điểm tại Quảng Nam, Đà Nẵng, Quảng Ngãi, định kỳ rà soát và cập nhật bản đồ nguy hiểm động đất 3 năm một lần.

Đối tượng nên tham khảo luận văn

  1. Các nhà nghiên cứu, nghiên cứu sinh và học viên cao học ngành Địa vật lý, Vật lý địa cầu và Địa chất kiến tạo: Luận văn cung cấp phương pháp luận mẫu mực về thuật toán nghịch đảo chu trình lặp của Kondorskaia, kỹ thuật lọc nhiễu số liệu thời gian truyền sóng và quy trình xây dựng hàm mục tiêu bằng phương pháp Monte-Carlo.

  2. Kỹ sư tư vấn thiết kế công trình thủy điện, công trình ngầm và hạ tầng ven biển miền Trung: Tài liệu giúp các kỹ sư xây dựng nắm bắt chính xác vận tốc nền đá ($V_P, V_S$) và các tham số chấn tiêu địa phương, phục vụ tính toán kháng chấn chính xác cho thân đập, hầm dẫn nước và các công trình giao thông trọng điểm.

  3. Cơ quan quản lý nhà nước về phòng chống thiên tai và an toàn năng lượng: Cán bộ thuộc Cục Quản lý Đê điều và Phòng chống Thiên tai, Ban chỉ huy Phòng thủ dân sự các tỉnh miền Trung có thể sử dụng dữ liệu luận văn để xây dựng kịch bản ứng phó sự cố đập và quy hoạch không gian an toàn trước thảm họa động đất.

  4. Giảng viên và sinh viên các trường đại học khối Khoa học Tự nhiên, Mỏ - Địa chất, Bách khoa: Đây là nguồn tài liệu tham khảo giá trị phục vụ giảng dạy các học phần Địa chấn học ứng dụng, Xử lý số liệu địa vật lý và Cấu trúc sâu vỏ Trái đất.

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

Tại sao phải xây dựng mô hình vận tốc riêng cho miền Trung thay vì dùng mô hình toàn cầu IASPEI91? Mô hình toàn cầu IASPEI91 được tính toán dựa trên số liệu trung bình toàn cầu nên không phản ánh chính xác cấu trúc thạch học phân tầng và các đới đứt gãy đặc thù tại Việt Nam. Sử dụng mô hình riêng giúp loại bỏ sai số hệ thống về thời gian truyền sóng, nâng cao độ chính xác định vị chấn tiêu cục bộ.

Số liệu đầu vào của luận văn được thu thập từ nguồn nào và quy mô ra sao? Nghiên cứu sử dụng tập hợp dữ liệu thời gian truyền sóng từ 275 trận động đất riêng biệt, ghi nhận bởi các trạm địa chấn Việt Nam kết hợp với các trạm địa chấn quốc tế lân cận, đảm bảo xây dựng đường tốc đồ bao phủ khoảng cách chấn tâm đến 800 km.

Quy trình lọc nhiễu 3 bước mang lại hiệu quả định lượng như thế nào? Quy trình lọc lặp 3 bước giúp giảm thiểu sai số thời gian $\Delta T$ từ mức sai lệch ban đầu lên tới 54 giây xuống dưới 0,5 giây, đồng thời nâng hệ số tương quan $R$ của phương trình sóng Pn từ 0,982 lên mức 0,999, đảm bảo dữ liệu đạt độ tin cậy thống kê cao nhất.

Tỷ số vận tốc $V_P/V_S$ trong nghiên cứu có ý nghĩa gì đối với công tác dự báo động đất? Tỷ số $V_P/V_S$ (dao động trong khoảng 1,67 đến 1,73) là thông số động học phản ánh độ cứng và trạng thái ứng suất của môi trường thạch quyển. Sự biến động bất thường của tỷ số này theo thời gian là dấu hiệu quan trọng để cảnh báo sự tích tụ ứng suất kiến tạo trước các trận động đất lớn.

Kết quả của luận văn có áp dụng được cho khu vực thủy điện Sông Tranh 2 không? Hoàn toàn áp dụng được. Hệ thống phương trình tốc đồ và mô hình vận tốc của luận văn là công cụ cơ bản giúp định vị chính xác độ sâu chấn tiêu các trận động đất kích thích tại Sông Tranh 2, làm rõ nguồn gốc kích hoạt do thấm nước hồ chứa hay do đứt gãy tự nhiên.

Kết luận

  • Luận văn đã giải quyết trọn vẹn bài toán xây dựng mô hình vận tốc truyền sóng địa chấn và hệ thống tốc đồ thời gian truyền sóng đầu tiên dành riêng cho vỏ Trái đất khu vực miền Trung Việt Nam.
  • Bộ 6 phương trình tốc đồ trung bình cho các pha sóng P, Pg, Pn, S, Sg, Sn đạt độ chính xác cao với hệ số tương quan $R$ từ 0,997 đến 1,000 trên cơ sở xử lý chuẩn hóa 275 sự kiện động đất trong bán kính 800 km.
  • Đã xác lập quy trình xử lý và lọc nhiễu 3 bước tối ưu, làm giảm độ lệch thời gian truyền sóng thực tế xuống dưới 0,5 giây, tạo tiền đề nâng cao chất lượng quy toán tham số chấn tiêu.
  • Lộ trình giai đoạn 2026-2030 cần tiếp tục mở rộng mạng quan trắc lên 25 trạm dải rộng và triển khai các tuyến thăm dò địa chấn sâu 2D/3D nhằm chi tiết hóa cấu trúc dị thường vận tốc vỏ Trái đất.
  • Các đơn vị quản lý vận hành công trình năng lượng, viện nghiên cứu và trường đại học hãy chủ động khai thác, tích hợp các kết quả khoa học của luận văn vào hệ thống quan trắc địa chấn thực tế để tối ưu hóa công tác cảnh báo sớm và giảm nhẹ thiên tai.