Tổng quan nghiên cứu

Thảm họa động đất và sóng thần luôn là hiểm họa địa chất khốc liệt trên phạm vi toàn cầu, điển hình như trận động đất sóng thần Sumatra năm 2004 gây thiệt hại hơn 220.000 sinh mạng, hay trận động đất Tōhoku năm 2011 tại Nhật Bản mạnh 9,0 độ Richter tàn phá nghiêm trọng kinh tế và môi trường. Tại Việt Nam, hoạt động địa chấn diễn ra phức tạp tại đới đứt gãy sông Hồng, sông Mã và hiện tượng động đất kích thích liên tiếp tại khu vực đập thủy điện Sông Tranh 2 thuộc tỉnh Quảng Nam năm 2012 đã đặt ra yêu cầu cấp thiết về công tác dự báo. Các phương pháp dự báo trung và dài hạn trước đây chủ yếu đưa ra xác suất trong khung thời gian từ vài chục năm đến hàng nghìn năm, chưa thể chỉ ra chính xác thời điểm xảy ra động đất và độ lớn chấn cấp trong tương lai gần.

Mục tiêu cốt lõi của nghiên cứu là xây dựng thuật toán mới, thiết lập sơ đồ khối và hoàn thiện phần mềm tính toán dự báo động đất theo mô hình thống kê xác suất. Đề tài tập trung giải quyết bài toán xác định kỳ vọng toán học và khoảng tin cậy của thời điểm phát sinh cùng độ lớn chấn cấp của trận động đất kế tiếp. Phạm vi nghiên cứu thực nghiệm bao gồm khu vực Tây Bắc Việt Nam từ vĩ độ 20,5 đến 23,0 độ Bắc, kinh độ 101,5 đến 105,5 độ Đông với 442 trận động đất độc lập giai đoạn 1970–2011, cùng danh mục địa chấn Đông Nam Á giai đoạn 1997–2007. Kết quả của luận văn đóng góp trực tiếp vào Đề tài khoa học độc lập cấp Nhà nước mã số VAST.09/11-12 của Viện Hàn lâm Khoa học và Công nghệ Việt Nam, cung cấp công cụ tự động hóa với độ chính xác kiểm nghiệm đạt trên 85% trong các thử nghiệm hồi quy lịch sử.

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 phát triển dựa trên mô hình chuỗi thời gian ngẫu nhiên của các sự kiện địa chấn, tiếp thu và cải biên phương pháp thống kê toán học do tác giả A. Grishin đề xuất. Bản chất của mô hình là xem tập hợp các trận động đất trong một tiểu vùng không gian xác định như một chuỗi thời gian ngẫu nhiên gồm các tham số thời điểm xuất hiện và độ lớn chấn cấp. Mô hình vận dụng quy luật phân phối xác suất chuẩn Gauss nhằm xác định khoảng tin cậy cho các đại lượng dự báo với xác suất bảo đảm được xác lập trước.

Bên cạnh đó, nghiên cứu kết hợp lý thuyết phân bố cực trị Gumbel loại I và loại III cùng định luật thực nghiệm Gutenberg-Richter biểu diễn mối quan hệ giữa tần suất xuất hiện động đất và chấn cấp theo phân phối Poisson. Ba khái niệm nền tảng được định nghĩa chặt chẽ trong mô hình bao gồm: sự kiện tựa là trận động đất cuối cùng được ghi nhận trong danh mục dữ liệu của tiểu vùng; trận động đất dự báo là sự kiện kế tiếp xảy ra sau sự kiện tựa cần tìm các giá trị kỳ vọng toán học; và kích thước tiểu vùng không gian được giới hạn theo bước lưới tọa độ nhằm đảm bảo tính đồng nhất về cấu trúc địa chất.

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

Nguồn dữ liệu phục vụ nghiên cứu được thu thập và chỉnh lý đồng bộ từ hệ thống mạng lưới trạm địa chấn quốc gia thuộc Viện Vật lý Địa cầu và Trung tâm Địa chấn Quốc tế ISC. Cỡ mẫu nghiên cứu tại vùng Tây Bắc bao gồm 1020 sự kiện địa chấn được ghi nhận từ năm 1970 đến năm 2011. Phương pháp chọn mẫu áp dụng thuật toán cửa sổ không gian và thời gian để tách lọc toàn bộ các nhóm tiền chấn và dư chấn, thu được 442 trận động đất độc lập có chấn cấp từ 3,0 độ Richter trở lên. Lý do lựa chọn phương pháp này là vì mô hình thống kê chỉ phát huy tính chính xác khi chuỗi dữ liệu đầu vào là các biến cố độc lập tuân theo phân phối ngẫu nhiên.

Quy trình phân tích được thực hiện qua các bước tính toán thống kê trượt trên từng khoảng thời gian chứa từ 10 đến 20 sự kiện để xác định kỳ vọng toán học và độ lệch chuẩn. Chương trình tính toán được lập trình bằng ngôn ngữ C# trên nền tảng .NET Framework với kiến trúc hướng đối tượng, cho phép nạp dữ liệu bảng 9 cột, tự động chia lưới không gian với bước nhảy 0,5 độ và xuất kết quả sang định dạng bảng tính cùng tập tin XML. Toàn bộ quy trình thử nghiệm và kiểm chứng thuật toán được hoàn thành trong giai đoạn 2011–2012.

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

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

Quá trình nghiên cứu và thử nghiệm phần mềm đã mang lại các phát hiện khoa học quan trọng:

Thứ nhất, thuật toán phân loại và tách lọc tiền chấn, dư chấn bằng cửa sổ không gian - thời gian đã xử lý thành công 1020 trận động đất tại vùng Tây Bắc giai đoạn 1970–2011, loại bỏ 578 sự kiện phụ thuộc (chiếm tỷ lệ 56,67%) để xác lập bộ dữ liệu chuẩn gồm 442 trận động đất độc lập.

Thứ hai, trong thử nghiệm lùi thời gian về mốc ngày 30/11/2001 với xác suất tin cậy lựa chọn 80% và ngưỡng tối thiểu 15 sự kiện trên mỗi tiểu vùng lưới 0,5 độ, chương trình đã tính toán dự báo 7 trận động đất. Kết quả đối chiếu với dữ liệu thực tế sau năm 2001 ghi nhận 6 trên 7 trận động đất (đạt độ chính xác 85,71%) xảy ra hoàn toàn trùng khớp trong khoảng thời gian và phạm vi chấn cấp dự báo.

Thứ ba, khi thiết lập thời gian tựa tại thời điểm nghiên cứu với danh mục đầy đủ đến năm 2011, chương trình đã dự báo được 9 vị trí phát sinh động đất tiềm năng tại Tây Bắc với độ lớn chấn cấp dao động từ 3,2 đến 5,1 độ Richter.

Thứ tư, khi so sánh hai tập danh mục lịch sử giai đoạn 1903–2011 và 1970–2011, kết quả chỉ ra rằng bộ dữ liệu trước năm 1970 cho sai số thời gian luôn mang giá trị âm với độ lệch tuyệt đối vượt quá 1000 ngày do chất lượng quan trắc thủ công thời kỳ đầu, trong khi bộ dữ liệu từ năm 1970 đến năm 2011 cho sai số cân bằng từ -11 ngày đến 471 ngày.

Thảo luận kết quả

Sự sai lệch giữa thời gian dự báo và thời gian xảy ra động đất thực tế phản ánh bản chất ngẫu nhiên của các quá trình tích tụ ứng suất kiến tạo trong thạch quyển. Khi tăng xác suất tin cậy mong muốn từ 70% lên 99%, khoảng thời gian dự báo bị nới rộng từ 1,037 lần lên 2,56 lần độ lệch chuẩn (tương đương khoảng thời gian dự báo tăng từ 150 ngày lên hơn 384 ngày), làm giảm đi tính ứng dụng thực tiễn của cảnh báo ngắn hạn. Dữ liệu nghiên cứu cho thấy sự tồn tại của mối tương quan tuyến tính thuận giữa độ chênh lệch biên thời gian tối đa với sai số thực tế, cho phép người dùng ước lượng trước độ tin cậy của kết quả dự báo ngay cả khi chưa có dữ liệu đối chứng trong tương lai.

Các kết quả tính toán có thể được trực quan hóa hiệu quả thông qua biểu đồ phân tán hồi quy biểu diễn mối quan hệ giữa sai số thời gian thực tế và biên độ dự báo, cũng như bản đồ đường đồng mức 2D thể hiện phân bố chấn cấp cực đại trên hệ tọa độ kinh vĩ tuyến. So với các thuật toán dự báo trung hạn quốc tế như CN hay M8 vốn chỉ đưa ra vùng nguy hiểm trong chu kỳ nhiều năm, mô hình thống kê được lập trình bằng C# trong luận văn có ưu điểm vượt trội khi lượng hóa rõ ràng ranh giới thời gian sớm nhất, thời gian muộn nhất và biên độ chấn cấp cho từng ô tọa độ cụ thể.

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

Dựa trên các kết quả đạt được, luận văn đưa ra 4 nhóm giải pháp mang tính ứng dụng thực tiễn cao:

Thứ nhất, nâng cấp mạng lưới quan trắc địa chấn quốc gia bằng việc lắp đặt bổ sung 25 trạm địa chấn số hóa băng rộng tại các khu vực hoạt động đứt gãy mạnh như Lai Châu, Điện Biên và Sơn La trong giai đoạn 2026–2030, do Viện Vật lý Địa cầu chủ trì nhằm tăng độ đầy đủ của danh mục động đất lên trên 95%.

Thứ hai, tích hợp thuật toán thống kê xác suất với các phương pháp vật lý kiến tạo, kết hợp 5 tham số cấu trúc vỏ Trái đất bao gồm bề dày thạch quyển, gradient vận tốc tân kiến tạo và mật độ dòng nhiệt để nâng cao độ chính xác dự báo chấn cấp cực đại vượt ngưỡng 5,5 độ Richter.

Thứ ba, hoàn thiện và nâng cấp phần mềm C# thành nền tảng dịch vụ web tích hợp hệ thông tin địa lý GIS, tự động tiếp nhận dữ liệu quan trắc thời gian thực và xuất báo cáo cảnh báo trong vòng 60 giây sau khi có sự kiện tựa mới phát sinh.

Thứ tư, Bộ Xây dựng và Tập đoàn Điện lực Việt Nam cần ứng dụng bản đồ phân vùng dự báo địa chấn từ mô hình này để rà soát, kiểm định tiêu chuẩn an toàn kháng chấn cho 100% các công trình thủy điện bậc thang và hồ chứa trọng điểm tại khu vực miền núi phía Bắc và miền Trung.

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

Nội dung luận văn mang lại giá trị tham khảo sâu sắc cho 4 nhóm đối tượng chuyên môn:

Nhóm thứ nhất là các nhà khoa học, nghiên cứu viên trong chuyên ngành Vật lý Địa cầu, Địa chấn học và Địa mạo kiến tạo, sử dụng luận văn làm tài liệu nghiên cứu về phương pháp xử lý chuỗi thời gian ngẫu nhiên và giải thuật tách lọc dữ liệu địa chấn.

Nhóm thứ hai là các kỹ sư phát triển phần mềm mô hình hóa và chuyên gia công nghệ thông tin địa lý, tham khảo kiến trúc hướng đối tượng trên ngôn ngữ C# và phương pháp tối ưu hóa thuật toán xử lý dữ liệu ma trận không gian đa chiều.

Nhóm thứ ba là cán bộ quản lý thuộc Cục Quản lý đê điều và Phòng chống thiên tai, cùng các kỹ sư quản lý an toàn hồ đập thủy điện, áp dụng công cụ để lập phương án ứng phó rủi ro động đất cục bộ và xây dựng kế hoạch sơ tán dân cư kịp thời.

Nhóm thứ tư là giảng viên, học viên cao học và sinh viên các ngành Khoa học Trái đất, Địa chất công trình và Toán tin ứng dụng tại các trường đại học khối khoa học tự nhiên và kỹ thuật, dùng làm học liệu chuyên đề về mô hình hóa tai biến thiên nhiên.

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

Phần mềm trong luận văn có thể đưa ra thời điểm dự báo chính xác đến từng ngày giờ cụ thể hay không? Mô hình thống kê cung cấp thời điểm dự báo trung tâm đi kèm khoảng tin cậy với giới hạn sớm nhất và muộn nhất dựa trên độ lệch chuẩn Gauss. Động đất mang tính ngẫu nhiên nên kết quả dự báo luôn là một khoảng thời gian xác suất chứ không phải một thời điểm tuyệt đối cố định.

Tại sao bước lọc bỏ tiền chấn và dư chấn lại đóng vai trò quyết định trong phương pháp này? Mô hình thống kê toán học dựa trên giả định các sự kiện địa chấn diễn ra độc lập theo quy luật phân phối ngẫu nhiên. Việc loại bỏ tiền chấn và dư chấn giúp loại trừ các tương tác dây chuyền cục bộ, đưa danh mục 1020 sự kiện về 442 trận động đất độc lập chuẩn xác.

Xác suất tin cậy lựa chọn ảnh hưởng như thế nào đến giá trị thực tế của kết quả dự báo? Khi chọn xác suất tin cậy 80%, độ mở rộng thời gian là 1,282 lần độ lệch chuẩn, mang lại cửa sổ thời gian hẹp và hữu ích. Nếu chọn xác suất 99%, khoảng thời gian bị nới rộng lên 2,56 lần độ lệch chuẩn, khiến khung thời gian quá dài và giảm giá trị cảnh báo.

Lý do tác giả lựa chọn ngôn ngữ C# thay vì các ngôn ngữ tính toán chuyên dụng như Matlab là gì? Ngôn ngữ C# chạy trên nền tảng .NET Framework sở hữu cấu trúc lập trình hướng đối tượng chặt chẽ, khả năng kiểm soát lỗi dữ liệu đầu vào xuất sắc, giao diện người dùng trực quan và thuận lợi cho việc đóng gói, chuyển giao công nghệ cho các cơ quan nghiên cứu địa phương.

Chất lượng của dữ liệu địa chấn trước năm 1970 ảnh hưởng như thế nào đến độ chính xác mô hình? Số liệu trước năm 1970 chủ yếu thu thập qua ghi chép thủ công và mạng lưới trạm thưa thớt, dẫn đến thiếu sót nhiều sự kiện nhỏ. Khi đưa vào tính toán, độ chênh lệch thời gian dự báo luôn bị âm và sai số tuyệt đối lớn hơn rất nhiều so với giai đoạn 1970–2011.

Kết luận

Luận văn đã giải quyết thành công các nhiệm vụ khoa học và thực tiễn trọng tâm:

  • Xây dựng thành công thuật toán và sơ đồ khối dự báo động đất ngắn hạn theo mô hình thống kê chuỗi thời gian ngẫu nhiên.
  • Phát triển hoàn chỉnh chương trình phần mềm chuyên dụng bằng ngôn ngữ C#, có khả năng tự động xử lý danh mục địa chấn và xác định khoảng tin cậy theo phân phối Gauss.
  • Chuẩn hóa bộ dữ liệu 442 trận động đất độc lập tại khu vực Tây Bắc Việt Nam giai đoạn 1970–2011 sau khi loại trừ 56,67% sự kiện tiền chấn và dư chấn.
  • Kiểm chứng thực nghiệm đạt độ chính xác 85,71% trong thử nghiệm hồi quy tại mốc lịch sử năm 2001 và dự báo thành công 9 vị trí tiềm năng cho giai đoạn tiếp theo.
  • Xác lập mối liên hệ định lượng giữa khoảng tin cậy lý thuyết và sai số thực tế, tạo cơ sở ước tính chất lượng dự báo khi ứng dụng trong điều kiện thực tiễn.

Trong lộ trình giai đoạn 2026–2030, hướng phát triển tiếp theo của nghiên cứu là tích hợp tham số trường địa vật lý và dữ liệu quan trắc vệ tinh để tự động hóa hoàn toàn hệ thống cảnh báo sớm tai biến địa chấn quốc gia. Các đơn vị nghiên cứu địa chất và quản lý hạ tầng công trình trọng điểm được khuyến nghị khai thác và ứng dụng bộ công cụ này nhằm tối ưu hóa công tác phòng chống và giảm thiểu rủi ro thiên tai.