Tổng quan nghiên cứu
Khu vực miền Bắc Việt Nam, đặc biệt là vùng Tây Bắc nằm về phía tây hệ thống đứt gãy Sông Hồng, là địa bàn có hoạt động kiến tạo hiện đại diễn ra vô cùng phức tạp. Theo các khảo sát địa chấn học, khu vực này chịu sự khống chế của hàng loạt hệ thống đứt gãy quy mô lớn như đứt gãy Sơn La, Sông Mã, Sông Cả và Lai Châu - Điện Biên, là nguồn gốc phát sinh nhiều trận động đất có độ lớn từ trung bình đến mạnh kèm theo chuỗi dư chấn kéo dài. Trước năm 2005, việc nghiên cứu chi tiết cơ cấu nguồn động đất tại Việt Nam gặp nhiều trở ngại do mật độ trạm quan sát còn thưa thớt và thiếu hụt nguồn số liệu địa chấn chất lượng cao.
Nhằm giải quyết bài toán cấp thiết này, nghiên cứu tập trung ứng dụng phương pháp nghịch đảo ten-xơ moment địa chấn trên nền tảng nguồn dữ liệu ghi nhận từ mạng lưới 25 trạm địa chấn dải rộng hiện đại được lắp đặt từ năm 2005 trong khuôn khổ hợp tác khoa học giữa Viện Vật lý Địa cầu và Viện các Khoa học Trái Đất Đài Loan. Mục tiêu cốt lõi của đề tài là xác định chính xác các tham số cơ cấu chấn tiêu cho các trận động đất tiêu biểu xảy ra trong giai đoạn 2009 - 2010 tại khu vực Tây Bắc, bao gồm trận động đất Bắc Yên năm 2009, trận Quan Sơn năm 2010 và trận Sốp Cộp năm 2010.
Kết quả nghiên cứu không chỉ làm sáng tỏ bức tranh kiến tạo địa tầng, cơ chế dịch trượt và hướng ứng suất kiến tạo khu vực, mà còn cung cấp cơ sở dữ liệu định lượng quan trọng phục vụ công tác đánh giá độ nguy hiểm địa chấn. Với việc chuẩn hóa quy trình phân tích dạng sóng dải rộng, nghiên cứu góp phần giảm thiểu sai số định vị nguồn chấn, phục vụ trực tiếp cho việc hoạch định an toàn cho hơn 10 công trình thủy điện trọng điểm bậc thang sông Đà.
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 nguồn địa chấn hiện đại và giải tích ten-xơ moment. Nguồn động đất kiến tạo được mô hình hóa dựa trên quá trình phá hủy môi trường đàn hồi trong lòng đất, giải phóng tích tụ năng lượng dưới dạng sóng khối P, sóng S và sóng mặt. Chuyển động dịch trượt trên mặt đứt gãy được biểu diễn thông qua ba tham số hình học cơ bản: góc phương vị (strike) dao động từ 0 đến 360 độ, góc dốc (dip) từ 0 đến 90 độ, và góc trượt (rake) từ âm 180 đến dương 180 độ. Về mặt lực học, dịch trượt đứt gãy tương đương với một hệ cặp ngẫu lực kép (double couple), đặc trưng bởi ten-xơ moment địa chấn đối xứng bậc hai gồm 9 thành phần trong không gian 3 chiều, trong đó có 6 thành phần độc lập.
Bên cạnh đó, lý thuyết hàm Green đóng vai trò cầu nối xác định phản ứng động lực học của môi trường đàn hồi nhiều lớp đối với một xung lực Dirac đơn vị theo thời gian và không gian. Phương pháp phân tích ten-xơ moment độ lệch còn cho phép bóc tách tỷ lệ giữa nguồn ngẫu lực kép thuần túy và thành phần lưỡng cực véctơ tuyến tính bù (CLVD), từ đó phản ánh chân thực bản chất vật lý của nguồn phát sinh chấn động.
Phương pháp nghiên cứu
Nguồn dữ liệu đầu vào sử dụng tập hợp băng ghi sóng địa chấn 3 thành phần ở định dạng chuẩn SAC thu thập từ 25 trạm dải rộng phân bố khắp miền Bắc Việt Nam. Cỡ mẫu nghiên cứu bao gồm 3 trận động đất chính có độ lớn từ M 3.0 trở lên kèm theo các dư chấn tiêu biểu xảy ra trong khung thời gian từ tháng 1 năm 2009 đến tháng 12 năm 2010. Phương pháp chọn mẫu ưu tiên các trạm địa chấn nằm trong bán kính dưới 300 km tính từ chấn tâm nhằm giảm thiểu tác động nhiễu của cấu trúc dị thường sâu, đảm bảo tỷ số tín hiệu trên nhiễu (SNR) đạt mức tối ưu và có tối thiểu từ 3 đến 7 trạm ghi nhận rõ ràng pha sóng.
Phương pháp phân tích trọng tâm là giải bài toán nghịch đảo ten-xơ moment trong miền tần số thấp bằng phần mềm ISOLA phiên bản 4.0 phát triển bởi các chuyên gia quốc tế năm 2009. Lý do lựa chọn ISOLA thay vì các chương trình truyền thống như FOCMEC hay PINV là bởi ISOLA sử dụng toàn bộ dạng sóng địa chấn dải rộng thay vì chỉ phụ thuộc vào dấu sóng P đầu tiên, hạn chế triệt để sai số do nhiễu môi trường. Để tính toán hàm Green chính xác, nghiên cứu tích hợp mô hình cấu trúc vận tốc vỏ Trái Đất 1D gồm 6 lớp đất đá chi tiết có độ sâu khảo sát đến 33 km tại khu vực miền Bắc Việt Nam, kết hợp thuật toán tìm kiếm lưới (grid-search) tối ưu hóa độ sâu và vị trí chấn tiêu.
Kết quả nghiên cứu và thảo luận
Những phát hiện chính
Quá trình nghịch đảo ten-xơ moment bằng phần mềm ISOLA đã mang lại những kết quả định lượng chuẩn xác về cơ cấu chấn tiêu cho chuỗi động đất tại Tây Bắc:
Thứ nhất, đối với trận động đất Mường La - Bắc Yên xảy ra năm 2009 cùng 2 dư chấn liên tiếp mang mã hiệu ML01 và ML02, kết quả tính toán cho thấy hệ số tương quan giữa băng ghi sóng thực tế và sóng lý thuyết tổng hợp đạt trên 82%. Mặt trượt chính của trận động đất thể hiện cơ chế trượt bằng kết hợp trượt bằng trái với góc dốc lớn xấp xỉ 70 độ, độ sâu chấn tiêu tối ưu nằm trong khoảng 8 km đến 12 km.
Thứ hai, trận động đất Quan Sơn - Thanh Hóa xảy ra năm 2010 được xác định có cơ chế trượt nghịch kết hợp trượt bằng, phản ánh trường ứng suất nén ép mạnh mẽ theo phương kinh tuyến với hệ số tương quan dạng sóng đạt 85% trên 6 trạm quan sát gần chấn tâm.
Thứ ba, trận động đất Sốp Cộp - Sơn La năm 2010 xác lập nghiệm cơ cấu chấn tiêu dạng trượt bằng phương tây bắc - đông nam, hoàn toàn trùng khớp với hướng phát triển của đứt gãy Sông Mã. So sánh nghiệm tính toán từ ISOLA với các thông số do Trung tâm Địa chấn Quốc tế (ISC) và Cục Khảo sát Địa chất Hoa Kỳ (USGS) công bố cho thấy độ sai lệch góc phương vị và góc dốc dưới 8%, khẳng định độ tin cậy vượt trội của mô hình vận tốc 6 lớp nội địa.
Thứ tư, phân tích thành phần ten-xơ moment độ lệch chỉ ra rằng thành phần ngẫu lực kép (double couple) chiếm tỷ lệ áp đảo từ 86% đến 94%, trong khi thành phần CLVD chỉ chiếm dưới 14%, chứng minh toàn bộ các sự kiện đều là động đất kiến tạo thuần túy do tích lũy ứng suất trượt cắt.
Thảo luận kết quả
Các kết quả tính toán có thể được trực quan hóa rõ nét thông qua biểu đồ quả cầu bãi biểu diễn các góc phần tư nén và giãn, đi kèm với bảng đối chiếu các tham số góc phương vị, góc dốc và góc trượt. Đồ thị dạng sóng địa chấn thể hiện sự khớp nối hoàn hảo giữa đường cong lý thuyết màu đỏ và đường tín hiệu thực tế màu đen trên cả ba thành phần chuyển vị thẳng đứng, hướng tâm và vuông góc.
Nguyên nhân dẫn đến sự trùng khớp cao này xuất phát từ việc lựa chọn dải tần số lọc sóng thích hợp từ 0.05 Hz đến 0.15 Hz, giúp loại bỏ các dao động vi địa chấn nhân tạo. Khi so sánh với các nghiên cứu dùng phương pháp dấu sóng P trước đây trong cùng khu vực, phương pháp nghịch đảo dạng sóng toàn phần giúp giảm độ bất định của góc dốc từ mức sai số 25 độ xuống dưới 5 độ. Ý nghĩa khoa học cốt lõi là đã chứng minh được nhánh phía tây đứt gãy Sông Hồng hiện vẫn đang trong giai đoạn tích tụ ứng suất kiến tạo tích cực, tạo cơ sở hiệu chỉnh các thông số nguồn cho bản đồ tai biến địa chấn khu vực.
Đề xuất và khuyến nghị
Dựa trên kết quả nghiên cứu thực nghiệm, luận văn đề xuất 4 nhóm giải pháp chiến lược nhằm nâng cao năng lực quan trắc và ứng phó rủi ro động đất:
Nâng cấp và mở rộng mạng lưới trạm quan sát: Đề nghị Viện Vật lý Địa cầu phối hợp cùng các cơ quan hữu quan đầu tư lắp đặt thêm từ 15 đến 20 trạm địa chấn dải rộng ba thành phần tại các khu vực trọng điểm Tây Bắc, đặc biệt dọc tuyến đứt gãy Sơn La và Điện Biên, hoàn thành trước năm 2028 nhằm hạ thấp ngưỡng ghi nhận động đất xuống mức M 2.0.
Tự động hóa quy trình xác định cơ cấu chấn tiêu: Xây dựng hệ thống tính toán nghịch đảo tự động kết hợp phần mềm ISOLA với hệ thống truyền dữ liệu thời gian thực, rút ngắn thời gian xử lý và công bố thông số nguồn động đất từ 120 phút xuống dưới 15 phút sau khi sự kiện xảy ra, phục vụ công tác cảnh báo khẩn cấp.
Hoàn thiện mô hình cấu trúc vận tốc 3D: Chuyển đổi từ mô hình vận tốc lát cắt 1D 6 lớp hiện tại sang mô hình sóng khối 3D chi tiết cho toàn bộ lãnh thổ miền Bắc nhằm giảm thiểu sai số hàm Green xuống dưới mức 3% trong các phân tích chấn tiêu phức tạp.
Tích hợp tham số nguồn vào quy chuẩn kháng chấn công trình: Khuyến nghị Bộ Xây dựng và Ban quản lý các công trình năng lượng áp dụng trực tiếp các tham số góc trượt và trường ứng suất động đất từ nghiên cứu này vào tính toán thiết kế kháng chấn cho các đập thủy điện lớn và hồ chứa bậc thang vùng Tây Bắc theo lộ trình 2026 - 2030.
Đối tượng nên tham khảo luận văn
Nhà nghiên cứu và chuyên gia Địa chấn học: Khai thác sâu phương pháp giải bài toán nghịch đảo ten-xơ moment dạng sóng dải rộng, thuật toán tìm kiếm lưới và cách thức xây dựng ma trận hàm Green cho môi trường đàn hồi nhiều lớp.
Kỹ sư thiết kế công trình thủy điện và hạ tầng giao thông: Sử dụng các tham số cơ cấu chấn tiêu, hướng ứng suất kiến tạo và độ sâu nguồn động đất thực tế để tính toán tải trọng động học, gia tốc nền cực đại phục vụ thiết kế kháng chấn cho cầu đường, hầm ngầm và đập dâng.
Cơ quan quản lý nhà nước về phòng chống thiên tai: Vận dụng số liệu định lượng về các đứt gãy hoạt động để cập nhật bản đồ phân vùng rủi ro động đất, xây dựng kịch bản ứng phó sự cố vỡ đập liên hoàn và quy hoạch dân cư an toàn tại các tỉnh miền núi phía Bắc.
Học viên cao học và sinh viên ngành Vật lý Địa cầu, Địa chất công trình: Sử dụng tài liệu như một cẩm nang mẫu mực về xử lý số liệu địa chấn dải rộng, cách vận hành phần mềm chuyên khảo ISOLA kết hợp bộ công cụ SEISAN trong nghiên cứu học thuật thực tế.
Câu hỏi thường gặp
Cơ cấu chấn tiêu động đất phản ánh những thông tin vật lý quan trọng nào? Cơ cấu chấn tiêu cung cấp thông tin chính xác về hướng mặt phẳng đứt gãy, phương dịch trượt, góc dốc và trục trường ứng suất kiến tạo nén ép hoặc giãn nở gây ra động đất. Nhờ đó, các nhà khoa học có thể xác định được nguyên nhân phát sinh chấn động là do hoạt động trượt bằng, trượt thuận hay trượt nghịch của vỏ Trái Đất.
Tại sao phương pháp nghịch đảo ten-xơ moment bằng phần mềm ISOLA lại vượt trội hơn các phương pháp truyền thống? Phương pháp truyền thống như FOCMEC phụ thuộc hoàn toàn vào việc đọc dấu sóng P đầu tiên vốn rất dễ sai lệch khi tín hiệu có nhiều tạp âm. Trong khi đó, ISOLA sử dụng toàn bộ biên độ và pha của dạng sóng dải rộng trong miền tần số thấp, kết hợp giải phương trình ma trận hàm Green giúp nâng cao độ chính xác lên trên 85%.
Điều kiện tối thiểu về dữ liệu để áp dụng thành công phần mềm ISOLA là gì? Để chạy chương trình ISOLA đạt độ tin cậy cao, sự kiện địa chấn cần có độ lớn tối thiểu từ M 3.0 trở lên, được ghi nhận đồng thời bởi ít nhất 3 trạm địa chấn dải rộng có tỷ số tín hiệu trên nhiễu tốt, phân bố góc phương vị bao quanh chấn tâm và nằm trong cự ly dưới 300 km.
Mô hình cấu trúc vỏ 1D 6 lớp của chuyên gia có vai trò gì trong bài toán nghịch đảo? Mô hình vận tốc 1D đóng vai trò là tham số môi trường đầu vào để tính toán hàm Green lý thuyết. Mô hình 6 lớp với độ sâu khảo sát 33 km phản ánh đúng đặc tính thạch học vỏ Trái Đất miền Bắc, giúp dạng sóng mô phỏng tổng hợp trùng khít với dạng sóng quan trắc thực tế tại các trạm.
Kết quả nghiên cứu này đóng góp gì cho sự an toàn của hệ thống thủy điện Tây Bắc? Nghiên cứu làm rõ cơ chế dịch trượt và độ sâu chấn tiêu của các đứt gãy chạy qua khu vực sông Đà và sông Mã. Những dữ liệu này giúp các kỹ sư hiệu chỉnh chính xác phổ phản ứng thiết kế và hệ số an toàn kháng chấn cho thân đập tại các nhà máy thủy điện lớn như Sơn La, Hòa Bình và Huổi Quảng.
Kết luận
- Xác định thành công cơ cấu chấn tiêu cho 3 trận động đất lớn và 2 dư chấn tại miền Bắc Việt Nam giai đoạn 2009 - 2010 bằng số liệu dải rộng chất lượng cao.
- Chứng minh tính hiệu quả vượt trội của phần mềm ISOLA kết hợp mô hình vận tốc 1D 6 lớp với độ tương quan dạng sóng đạt trên 80%.
- Làm sáng tỏ hoạt động kiến tạo trượt bằng kết hợp trượt nghịch của các hệ thống đứt gãy chính phía tây đứt gãy Sông Hồng.
- Khẳng định tỷ lệ thành phần ngẫu lực kép chiếm trên 86%, xác nhận bản chất nguồn gốc động đất thuần túy do tích lũy ứng suất trượt cắt kiến tạo.
- Đề xuất khung giải pháp hoàn thiện quy trình quan trắc địa chấn tự động và cập nhật tiêu chuẩn an toàn công trình hạ tầng trọng điểm.
Kế hoạch triển khai tiếp theo trong giai đoạn 2026 - 2028 tập trung vào việc mở rộng mạng lưới trạm thu thập dữ liệu thời gian thực và xây dựng mô hình cấu trúc vận tốc 3D cho toàn vùng. Hãy liên hệ và hợp tác cùng các chuyên gia địa vật lý để ứng dụng ngay các dữ liệu cơ cấu chấn tiêu chuẩn xác này vào công tác tính toán kháng chấn và quy hoạch an toàn thiên tai cho dự án của bạn.