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à nơi có hoạt động kiến tạo phức tạp với mật độ động đất cao nhất cả nước. Trong giai đoạn từ năm 2005 đến năm 2012, mạng lưới quan trắc địa chấn dải rộng đã ghi nhận hơn 248 trận động đất có độ lớn cục bộ từ 3.0 độ Richter trở lên. Nhiều trận động đất trong số này đã gây rung chấn mạnh trên bề mặt đạt tới cấp VII theo thang chấn cấp quốc tế MSK-64, đe dọa trực tiếp đến an toàn dân sinh và các công trình hạ tầng trọng điểm.

Vấn đề cốt lõi đặt ra là việc xác định chính xác cơ cấu chấn tiêu của các trận động đất tại Việt Nam trước đây gặp nhiều hạn chế do mạng lưới trạm quan sát thưa thớt và phương pháp truyền thống chủ yếu dựa vào dấu sóng tới dọc ban đầu. Nhằm giải quyết thách thức này, nghiên cứu đặt mục tiêu ứng dụng phương pháp nghịch đảo ten-xơ moment địa chấn trong miền tần số thấp từ số liệu băng ghi địa chấn dải rộng của mạng lưới 25 trạm quan sát hiện đại. Phạm vi nghiên cứu tập trung vào chuỗi động đất điển hình xảy ra tại khu vực Tây Bắc trong hai năm 2009 và 2010, bao gồm cụm động đất Mường La – Bắc Yên, Quan Sơn và Sốp Cộp.

Ý nghĩa khoa học và thực tiễn của công trình thể hiện ở việc nâng cao độ chính xác trong xác định hình thái đứt gãy sinh chấn, cung cấp nghiệm cơ cấu chấn tiêu với hệ số tương quan dạng sóng đạt tới 80%. Kết quả này là cơ sở dữ liệu định lượng quan trọng giúp làm sáng tỏ bức tranh trường ứng suất kiến tạo hiện đại, phục vụ trực tiếp cho công tác đánh giá độ nguy hiểm địa chấn và thiết kế kháng chấn công trình.

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 đàn hồi và mô hình phát xạ sóng khối của các học giả chuyên ngành. Trọng tâm lý thuyết bao gồm:

Mô hình đứt gãy phẳng được xác định bởi ba tham số hình học cơ bản: góc phương vị đường phương, góc dốc mặt trượt và góc trượt của cánh treo so với cánh nằm. Bản chất quá trình giải phóng năng lượng tại chấn tiêu được mô hình hóa toán học thông qua ten-xơ moment địa chấn bậc hai với 9 thành phần lực tương đương, trong đó có 6 thành phần độc lập do tính chất đối xứng. Độ lớn động đất được lượng hóa qua moment địa chấn vô hướng theo quan hệ tích số giữa độ cứng môi trường, diện tích mặt trượt và khoảng dịch trượt trung bình.

Lý thuyết hàm Green trong môi trường đàn hồi phân lớp miêu tả phản ứng dịch chuyển của môi trường đất đá dưới tác động của một nguồn xung lực đơn vị. Băng ghi địa chấn lý thuyết được biểu diễn dưới dạng tích chập tuyến tính giữa hàm thời gian nguồn, hàm Green cấu trúc môi trường và hàm đáp ứng thiết bị ghi nhận.

Phép phân tích ten-xơ moment phân rã nghiệm thành phần đẳng hướng mô tả sự thay đổi thể tích và phần ten-xơ lệch. Phần ten-xơ lệch tiếp tục được phân tách thành nguồn ngẫu lực kép đại diện cho dịch trượt kiến tạo thuần túy và nguồn lưỡng cực véc-tơ tuyến tính bù đại diện cho các phá hủy phức tạp hoặc hiện tượng dị thường.

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

Nguồn dữ liệu nghiên cứu sử dụng tập số liệu dạng sóng ba thành phần ghi nhận từ mạng lưới 25 trạm địa chấn dải rộng được lắp đặt đồng bộ tại miền Bắc Việt Nam. Quy trình chọn mẫu tập trung vào các sự kiện động đất có độ lớn từ 2.8 đến 4.6 độ Richter, có tối thiểu từ 3 đến 7 trạm ghi nhận đạt tỷ số tín hiệu trên nhiễu cao và phân bố góc phương vị bao quanh chấn tâm đồng đều.

Quy trình phân tích thực hiện thông qua chương trình nghịch đảo ten-xơ moment ISOLA phiên bản 4.0 kết hợp với bộ công cụ SEISAN. Dữ liệu vận tốc thô định dạng chuẩn được chuyển đổi sang dạng dịch chuyển và lọc qua dải tần số thấp từ 0.05 Hz đến 0.15 Hz để loại bỏ nhiễu tán xạ trường gần. Mô hình cấu trúc vỏ Trái Đất một chiều 6 lớp với độ sâu khảo sát đến 33 km và dải vận tốc sóng dọc từ 5.6 km/s đến 8.1 km/s được sử dụng để giải bài toán hàm Green. Phương pháp tìm kiếm dạng lưới theo không gian và độ sâu chấn tiêu kết hợp thuật toán bình phương tối thiểu được lựa chọn nhằm tìm ra nghiệm ten-xơ moment tối ưu nhất, hạn chế tối đa sự phụ thuộc chủ quan vào việc đọc dấu sóng ban đầu.

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

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

Phân tích nghịch đảo dạng sóng địa chấn cho chuỗi sự kiện động đất Tây Bắc đã mang lại các kết quả định lượng chi tiết:

Trận động đất chủ chấn Mường La – Bắc Yên xảy ra ngày 26 tháng 11 năm 2009 có độ lớn cục bộ 3.9 độ Richter và độ lớn moment đạt 3.6. Nghiệm nghịch đảo xác định chấn tiêu nằm ở độ sâu nông 5 km với cơ cấu trượt bằng chiếm ưu thế tuyệt đối. Mặt phẳng đứt gãy thứ nhất có góc phương vị 302 độ, góc dốc 80 độ cắm về phía Bắc và góc trượt 178 độ. Độ khớp dạng sóng giữa lý thuyết và thực tế đạt hệ số tương quan 80%, trong đó tỷ lệ ngẫu lực kép đạt mức 99.8% và thành phần lưỡng cực tuyến tính bù chỉ chiếm 0.2%.

Hai trận dư chấn đi kèm tại khu vực Bắc Yên gồm dư chấn thứ nhất có độ lớn 3.5 độ Richter và dư chấn thứ hai có độ lớn 2.8 độ Richter đều cho nghiệm cơ cấu chấn tiêu đồng nhất với trận chủ chấn. Mặt đứt gãy của dư chấn thứ nhất có góc phương vị 303 độ, góc dốc 84 độ và góc trượt âm 177 độ tại độ sâu 5 km, xác nhận sự tiếp diễn chuyển động trượt bằng trên cùng một đới nứt nẻ kiến tạo.

Trận động đất Quan Sơn – Thanh Hóa xảy ra ngày 19 tháng 9 năm 2010 và trận động đất Sốp Cộp – Sơn La xảy ra ngày 30 tháng 12 năm 2010 có độ sâu chấn tiêu dao động từ 5 km đến 12 km. Cả hai sự kiện đều thể hiện cơ cấu nguồn dạng trượt bằng kết hợp thành phần trượt thuận nhỏ, phù hợp với động học của các đứt gãy cấp khu vực. Khi so sánh nghiệm cơ cấu chấn tiêu trận Sốp Cộp tính bằng chương trình ISOLA với dữ liệu công bố của Cục Khảo sát Địa chất Hoa Kỳ và Trung tâm Địa chấn Quốc tế, sai lệch về góc phương vị và góc dốc đều nhỏ hơn 5 độ, chứng minh độ tin cậy vượt trội của thuật toán áp dụng.

Thảo luận kết quả

Cơ cấu chấn tiêu của chuỗi động đất xác định trục ứng suất nén ép cực đại có phương á kinh tuyến Bắc - Nam và trục ứng suất tách giãn có phương á vĩ tuyến Đông - Tây. Điều này giải thích rõ ràng cơ chế trượt bằng phương Tây Bắc - Đông Nam của các đới đứt gãy hoạt động mạnh như đứt gãy Sơn La và đứt gãy Sông Mã dưới tác động của trường ứng suất kiến tạo khu vực Đông Nam Á.

Dữ liệu nghiên cứu được biểu diễn trực quan qua biểu đồ cơ cấu chấn tiêu dạng quả bóng bãi biển hai màu thể hiện rõ các góc phần tư nén ép và tách giãn. Đồng thời, bảng tổng hợp tham số nguồn và đồ thị so sánh dạng sóng ba thành phần ghi nhận giữa đường lý thuyết màu đỏ và đường thực tế màu đen trên từng trạm quan trắc minh chứng cho sự hội tụ nghiệm với sai số bình phương tối thiểu rất thấp. Việc sử dụng trọn vẹn dạng sóng dải tần số thấp từ 0.05 Hz đến 0.15 Hz đã khắc phục triệt để nhược điểm của các phương pháp cổ điển như FOCMEC vốn hay gặp sai số khi tín hiệu sóng dọc pha đầu bị che khuất bởi nhiễu môi trường.

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

Dựa trên kết quả nghiên cứu thực nghiệm, các giải pháp mang tính hành động cụ thể được đề xuất nhằm nâng cao năng lực quan trắc và phòng chống rủi ro thiên tai:

Thứ nhất, Viện Vật lý Địa cầu cần chủ trì nâng cấp và mở rộng mạng lưới trạm địa chấn dải rộng từ 25 trạm hiện nay lên tối thiểu 40 trạm phủ kín các vùng đứt gãy hoạt động trọng điểm tại Tây Bắc trong giai đoạn từ năm 2026 đến năm 2030, nhằm giảm khoảng cách trung bình giữa các trạm xuống dưới 50 km.

Thứ hai, Trung tâm Báo tin động đất và Cảnh báo sóng thần cần chuẩn hóa quy trình tự động hóa nghịch đảo ten-xơ moment bằng chương trình ISOLA tích hợp vào hệ thống giám sát thời gian thực. Mục tiêu đặt ra là công bố cơ cấu chấn tiêu và độ sâu chấn tiêu chính xác trong vòng 15 phút sau khi ghi nhận các sự kiện có độ lớn từ 3.0 độ Richter trở lên.

Thứ ba, các cơ quan nghiên cứu chuyên ngành cần phối hợp triển khai đề tài xây dựng mô hình cấu trúc vận tốc vỏ Trái Đất ba chiều chi tiết cho từng phân vùng địa chất Tây Bắc trong vòng 3 năm tới, giúp giảm sai số tính toán hàm Green từ 20% hiện tại xuống dưới 10%.

Thứ tư, Bộ Xây dựng và Tập đoàn Điện lực Việt Nam cần cập nhật tham số cơ cấu chấn tiêu và gia tốc nền cực đại vào quy chuẩn xây dựng công trình ngầm, đập thủy điện bậc thang trên dòng sông Đà và sông Mã trước năm 2028, đảm bảo công trình chịu được rung động chấn động cấp VII đến cấp VIII theo thang MSK-64.

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

Công trình nghiên cứu mang lại giá trị tham khảo chuyên sâu và ứng dụng thực tế cho bốn nhóm đối tượng chính:

Nhóm nghiên cứu sinh và học viên cao học chuyên ngành Vật lý Địa cầu, Địa chất kiến tạo và Địa kỹ thuật: Sử dụng luận văn như tài liệu học thuật tiêu chuẩn về phương pháp giải bài toán nghịch đảo nguồn địa chấn, tính toán hàm Green và xử lý số liệu địa chấn dải rộng.

Chuyên viên phân tích dữ liệu địa chấn tại các viện nghiên cứu và trung tâm cảnh báo thiên tai: Ứng dụng quy trình xử lý tín hiệu và các bước thao tác trên phần mềm ISOLA để xác định nhanh các thông số cơ cấu chấn tiêu phục vụ báo cáo nghiệp vụ.

Kỹ sư thiết kế kết cấu công trình thủy điện, giao thông và dân dụng vùng núi phía Bắc: Khai thác số liệu về cơ cấu đứt gãy, độ sâu chấn tiêu và trường ứng suất nén ép để tính toán hệ số kháng chấn thực tế cho các dự án xây dựng hạ tầng.

Cơ quan quản lý nhà nước về phòng chống thiên tai cấp tỉnh và khu vực: Tham khảo bản đồ phân bố chấn tâm và đường đẳng chấn để xây dựng kịch bản ứng phó khẩn cấp và quy hoạch không gian phát triển kinh tế an toàn.

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

Phương pháp nghịch đảo ten-xơ moment bằng ISOLA có ưu điểm gì vượt trội so với phương pháp xác định dấu sóng đầu tiên? Phương pháp nghịch đảo dạng sóng toàn phần tận dụng toàn bộ biên độ và pha sóng dải tần thấp thay vì chỉ dùng một điểm dấu sóng P đầu tiên. Nhờ đó, phương pháp đạt hệ số tương quan trên 80% ngay cả khi mạng trạm thưa thớt và biên độ sóng bị nhiễu.

Tại sao dải tần số từ 0.05 Hz đến 0.15 Hz lại được lựa chọn trong quá trình nghịch đảo số liệu địa chấn địa phương? Dải tần số thấp từ 0.05 Hz đến 0.15 Hz ít bị ảnh hưởng bởi tính bất đồng nhất quy mô nhỏ của môi trường vỏ Trái Đất. Việc lọc dải tần này giúp sóng địa chấn giữ được dạng cơ bản, tối ưu hóa độ chính xác khi đối sánh với hàm Green một chiều.

Tỷ lệ thành phần ngẫu lực kép đạt 99.8% trong trận động đất Bắc Yên phản ánh điều gì về mặt địa chất? Tỷ lệ ngẫu lực kép xấp xỉ 100% chứng minh nguồn phát sinh động đất hoàn toàn do quá trình tích lũy và giải phóng ứng suất trượt cắt kiến tạo thuần túy dọc theo mặt phẳng đứt gãy, loại trừ các nguyên nhân sụt lún rỗng hoặc hoạt động phun trào thể tích.

Mô hình vỏ Trái Đất 6 lớp của khu vực miền Bắc đóng vai trò như thế nào trong kết quả tính toán? Mô hình cấu trúc 6 lớp với độ sâu 33 km và dải vận tốc sóng dọc từ 5.6 km/s đến 8.1 km/s cung cấp các tham số truyền sóng chính xác để tính toán hàm Green. Mô hình này giúp loại bỏ sai số lan truyền sóng từ chấn tiêu tới trạm thu.

Kết quả cơ cấu chấn tiêu đóng góp gì cho an toàn của chuỗi thủy điện lớn tại khu vực Tây Bắc? Kết quả xác định chính xác phương vị đứt gãy và độ sâu chấn tiêu nông ở mức 5 km giúp các nhà quản lý đánh giá đúng mức độ rung chấn cực đại cấp VII MSK-64, từ đó kiểm định độ an toàn đập và cửa xả lũ định kỳ.

Kết luận

Luận văn đã giải quyết trọn vẹn bài toán xác định cơ cấu chấn tiêu động đất bằng số liệu địa chấn dải rộng tại miền Bắc Việt Nam với các đóng góp cốt lõi:

  • Thiết lập quy trình chuẩn xử lý và nghịch đảo ten-xơ moment dạng sóng dải tần thấp cho các trận động đất địa phương có độ lớn từ 2.8 độ Richter trở lên.
  • Khẳng định tính ưu việt của chương trình ISOLA với độ khớp dạng sóng lý thuyết và thực tế đạt hệ số tương quan 80%.
  • Xác định chính xác cơ cấu trượt bằng chiếm ưu thế của chuỗi động đất Mường La – Bắc Yên năm 2009 và các sự kiện tại Thanh Hóa, Sơn La năm 2010 ở độ sâu nông từ 5 km đến 12 km.
  • Cung cấp bằng chứng định lượng về trường ứng suất kiến tạo khu vực với hướng nén ép chủ đạo Bắc - Nam và tách giãn Đông - Tây.
  • Chứng minh độ tin cậy của mạng lưới 25 trạm địa chấn dải rộng quốc gia trong việc cung cấp dữ liệu phục vụ nghiên cứu địa chấn chuyên sâu.

Trong giai đoạn tiếp theo từ năm 2026 đến năm 2030, việc mở rộng tích hợp mô hình vận tốc không gian ba chiều và tự động hóa quy trình phân tích sẽ là bước đi tất yếu. Các tổ chức nghiên cứu và kỹ sư công trình hãy khai thác ngay nguồn cơ sở dữ liệu quý giá này để tích hợp vào các đồ án quy hoạch hạ tầng an toàn trước hiểm họa động đất.