Tổng quan nghiên cứu

Khu vực miền Bắc Việt Nam nằm trong bối cảnh kiến tạo phức tạp do ảnh hưởng từ sự va húc giữa các mảng thạch quyển với hệ thống đới đứt gãy hoạt động kéo dài trên 1000 km như đứt gãy Sông Hồng, đứt gãy Lai Châu - Điện Biên, đứt gãy Sơn La và đứt gãy Mường La - Bắc Yên. Trong lịch sử quan trắc, khu vực Tây Bắc từng xảy ra các trận chấn động nghiêm trọng như động đất Tuần Giáo năm 1983 với độ lớn M = 6.7 gây rung động cấp VII tại Điện Biên cách tâm chấn 70 km, và trận động đất năm 2001 với độ lớn M = 5.3 gây chấn động cấp VII ở khoảng cách 20 km. Mặc dù hiểm họa địa chấn luôn hiện hữu, mạng lưới máy đo gia tốc nền phục vụ tính toán kháng chấn tại nước ta mới bắt đầu được triển khai từ năm 2000 với số lượng trạm rất thưa thớt, dẫn đến tình trạng thiếu hụt dữ liệu định lượng phục vụ công tác thiết kế công trình trọng điểm.

Mục tiêu cốt lõi của nghiên cứu là xây dựng quy trình chuyển đổi số liệu vận tốc dải rộng sang gia tốc nền thực tế, đồng thời tiến hành kiểm nghiệm tính tương thích của mô hình suy giảm dao động thế hệ mới Campbell (2008) đối với điều kiện địa chất khu vực Tây Bắc và phụ cận. Nghiên cứu thực hiện trên tập dữ liệu gồm 8 trận động đất có độ lớn Ml từ 4.0 trở lên trong giai đoạn 2009 - 2014 với tổng số 165 băng ghi địa chấn. Kết quả kiểm nghiệm đạt hệ số tương quan chéo trên 0.85, cung cấp luận cứ khoa học chính xác để chuẩn hóa tham số gia tốc đỉnh phục vụ quy hoạch hạ tầng và thiết kế kháng chấn các công trình năng lượng lớn.

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 vận dụng lý thuyết lan truyền sóng đàn hồi trong môi trường liên tục và lý thuyết tắt dần chấn động theo mô hình thế hệ mới (Next Generation Attenuation - NGA) do Trung tâm Nghiên cứu Kỹ thuật Địa chấn Thái Bình Dương phát triển. Động lực học công trình được mô tả thông qua mô hình phản ứng đàn hồi của hệ một bậc tự do (SDOF) dưới tác động của các kích động nền ngẫu nhiên.

Các khái niệm then chốt bao gồm:

  • Gia tốc nền cực đại (Peak Ground Acceleration - PGA): Đại lượng đo lường lực quán tính cực đại tác động trực tiếp lên kết cấu công trình.
  • Phổ phản ứng gia tốc (Spectral Acceleration - SA): Tham số biểu thị phản ứng động lực lớn nhất của kết cấu tại từng chu kỳ dao động riêng biệt, đặc biệt quan trọng trong dải tần số nguy hiểm từ 0 đến 10 Hz.
  • Vận tốc truyền sóng ngang trung bình trong 30 m đất mặt (VS30): Chỉ số phân loại điều kiện địa kỹ thuật nền móng theo tiêu chuẩn quốc tế.

Mô hình suy giảm chấn động Campbell (2008) tích hợp cấu trúc toán học phi tuyến phức hợp với hơn 20 hệ số hồi quy, phản ánh toàn diện các hàm thành phần: hàm độ lớn động đất, hàm suy giảm hình học theo khoảng cách ngắn nhất đến mặt trượt đứt gãy (RRUP), hàm cơ chế nguồn trượt bằng và hàm phi tuyến nền đất.

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

Nguồn dữ liệu của luận văn được tổng hợp từ hệ thống trạm địa chấn quốc gia do Viện Vật lý Địa cầu vận hành trong giai đoạn 2009 - 2014. Tổng dung lượng mẫu khảo sát gồm 165 băng ghi từ 8 trận động đất độ lớn Ml từ 4.0 đến 5.3, bao gồm 150 băng ghi vận tốc dải tần rộng từ cảm biến STS-2 và 15 băng ghi gia tốc 3 thành phần từ thiết bị EpiSensor ES-T có dải động học đạt 155 dB.

Phương pháp chọn mẫu áp dụng là phương pháp lấy mẫu mục tiêu kết hợp kiểm tra tỷ số tín hiệu trên nhiễu (Signal-to-Noise Ratio - SNR). Quy trình phân tích số liệu diễn ra qua các giai đoạn:

  1. Thuật toán biến đổi Fourier nhanh (FFT) trên môi trường MATLAB được sử dụng để khảo sát phổ biên độ tín hiệu với cửa sổ phân tích 0.2 giây trước sóng P nhằm loại bỏ các bản ghi nhiễu cao.
  2. Phần mềm chuyên dụng SAC (Seismic Analysis Code) thực hiện giải tích chập khử đáp ứng xung thiết bị đo và vi phân số tín hiệu vận tốc trong dải tần lọc thông dải 0.01 - 40 Hz để xuất ra chuỗi thời gian gia tốc nền.
  3. Sử dụng hàm tương quan chéo đa kênh để so sánh độ tương thích pha và biên độ giữa tín hiệu gia tốc chuyển đổi và tín hiệu đo thực.

Lý do lựa chọn phương pháp phân tích này xuất phát từ thực tiễn mạng lưới quan trắc vận tốc có mật độ dày và lịch sử vận hành liên tục từ năm 2005, trong khi số lượng trạm gia tốc kế trực tiếp còn quá hạn chế.

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

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

Quá trình phân tích thực nghiệm và kiểm nghiệm mô hình đã mang lại các kết quả định lượng cụ thể:

  • Thuật toán chuyển đổi tín hiệu từ băng ghi vận tốc sang băng ghi gia tốc thông qua phần mềm SAC đạt độ chuẩn xác vượt trội. Tại trạm địa chấn Sơn La đối với trận động đất ngày 19/07/2014 (độ lớn Ml = 4.4), hệ số tương quan chéo giữa băng ghi chuyển đổi và băng ghi thực tế đạt trên 0.90 ở cả ba thành phần không gian (Bắc - Nam, Đông - Tây, thẳng đứng).
  • Đường cong tắt dần chấn động Campbell (2008) thể hiện sự tương thích rất cao với đặc tính dao động nền thực tế tại Tây Bắc, với hơn 85% điểm số liệu thực nghiệm phân bố hoàn toàn bên trong giới hạn một độ lệch chuẩn (+/- 1 sigma) của mô hình.
  • Phổ gia tốc phản ứng SA tính toán theo mô hình Campbell (2008) cho trận động đất Điện Biên (Mw = 4.8) ở cự ly nguồn chấn RRUP = 15 km phản ánh sát thực tế phân bố năng lượng trong dải chu kỳ dao động kết cấu từ 0.1 đến 2.0 giây.
  • Khắc phục triệt để sai số của các nghiên cứu trước: mô hình Nguyễn Lê Minh (2012) vốn chỉ giới hạn ở các trận động đất nhỏ (độ lớn dưới 4.0) và mô hình Trần Việt Hùng (2012) có 90% dữ liệu nguồn từ Nhật Bản thường đánh giá thấp gia tốc đỉnh ở vùng gần đứt gãy dưới 30 km.

Thảo luận kết quả

Sự phù hợp của mô hình Campbell (2008) bắt nguồn từ tính tương đồng sâu sắc về cơ chế kiến tạo. Động đất tại miền Bắc Việt Nam và khu vực Vân Nam (Trung Quốc) chủ yếu hình thành từ các đứt gãy trượt bằng trong vỏ lục địa nông ở độ sâu dưới 20 km. Dữ liệu này đồng nhất với tập mẫu 2155 băng ghi từ 170 trận động đất trên toàn cầu mà dự án NGA sử dụng để xây dựng công thức hồi quy.

Trong thực tế phân tích, toàn bộ dữ liệu dao động nền được biểu diễn trực quan qua các biểu đồ phân bố điểm gia tốc đỉnh (PGA) theo khoảng cách chấn tiêu và bảng đối chiếu sai số tương đối giữa các thành phần nằm ngang và thẳng đứng. Việc kiểm chứng thành công phương pháp chuyển đổi dữ liệu vận tốc sang gia tốc mở ra khả năng khai thác hơn 90% khối lượng băng ghi sẵn có từ mạng lưới địa chấn quốc gia, tạo nguồn đầu vào chuẩn hóa cho các bài toán phân tích rủi ro địa chấn mà không phải chờ đợi hàng thập kỷ để tích lũy dữ liệu gia tốc kế trực tiếp.

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

  • Tích hợp chuẩn hóa module giải tích chập và chuyển đổi tự động dữ liệu vận tốc sang gia tốc nền trong dải tần số 0.01 - 40 Hz vào hệ thống xử lý thời gian thực của Viện Vật lý Địa cầu trong vòng 12 tháng tới, nhằm nâng cao 100% dung lượng dữ liệu phục vụ nghiên cứu kháng chấn.
  • Ứng dụng ngay đường cong tắt dần Campbell (2008) vào công tác tính toán tải trọng động đất và hiệu chỉnh bản đồ phân vùng nguy hiểm địa chấn cho khu vực miền Bắc, hướng tới mục tiêu giảm sai số dự báo dao động nền xuống dưới 15% trong giai đoạn 2026 - 2028 dưới sự phối hợp của Bộ Xây dựng và Viện Hàn lâm Khoa học và Công nghệ Việt Nam.
  • Đầu tư lắp đặt bổ sung 10 trạm gia tốc kế số chuyên dụng EpiSensor ES-T tại các vị trí trọng yếu dọc theo đới đứt gãy Mường La - Bắc Yên và đới đứt gãy Sông Mã, đảm bảo khoảng cách trung bình giữa các trạm quan sát dưới 50 km trước năm 2030.
  • Xây dựng cơ sở dữ liệu địa kỹ thuật công trình chuẩn hóa tham số vận tốc sóng cắt VS30 cho toàn bộ các đô thị loại II trở lên và các bậc thang thủy điện vùng Tây Bắc, hoàn thành kiểm định an toàn đập với độ tin cậy đạt mức 95%.

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

  • Kỹ sư thiết kế kết cấu và chuyên gia tư vấn xây dựng công trình: Khai thác trực tiếp các tham số gia tốc đỉnh PGA và phổ gia tốc SA để tính toán tải trọng kháng chấn cho nhà cao tầng, cầu lớn và đập thủy điện trọng lực.
  • Nhà nghiên cứu và giảng viên chuyên ngành Địa vật lý, Địa kỹ thuật công trình: Tiếp cận phương pháp luận xử lý số liệu địa chấn tiên tiến, thuật toán lọc tín hiệu trên MATLAB và quy trình kiểm nghiệm mô hình suy giảm chấn động thế hệ mới.
  • Cơ quan quản lý nhà nước về phòng chống thiên tai và quy hoạch đô thị: Sử dụng các kết quả đánh giá thực nghiệm làm căn cứ sửa đổi quy chuẩn xây dựng kháng chấn và lập bản đồ vi phân vùng hiểm họa động đất cấp địa phương.
  • Doanh nghiệp thẩm định rủi ro và các tổ chức bảo hiểm công trình: Ứng dụng mô hình dự báo gia tốc nền để mô phỏng kịch bản thiệt hại kết cấu hạ tầng với chu kỳ lặp 475 năm, phục vụ định phí bảo hiểm chính xác và xây dựng phương án ứng phó khẩn cấp.

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

Tại sao cần chuyển đổi dữ liệu từ trạm đo vận tốc sang gia tốc nền? Do số lượng trạm đo gia tốc chuyên dụng tại Việt Nam còn rất hạn chế với chỉ 15 băng ghi đạt chuẩn trong tập mẫu nghiên cứu, việc chuyển đổi 150 băng ghi vận tốc dải rộng sang gia tốc nền thông qua phần mềm SAC với hệ số tương quan trên 0.90 giúp khai thác tối đa nguồn dữ liệu địa chấn sẵn có phục vụ kỹ thuật kháng chấn.

Mô hình suy giảm Campbell (2008) có ưu điểm gì vượt trội so với các công thức cũ? Mô hình Campbell (2008) được xây dựng từ tập dữ liệu khổng lồ gồm 2155 băng ghi của 170 trận động đất, tích hợp hơn 20 tham số hồi quy phi tuyến tính bao gồm cơ chế đứt gãy trượt bằng và vận tốc sóng cắt VS30, giúp giảm thiểu độ phân tán và khắc phục lỗi dự báo sai lệch ở cự ly chấn tâm gần.

Dải tần số dao động nào của động đất gây phá hoại công trình lớn nhất? Dải tần số nguy hiểm nhất nằm trong khoảng từ 0 đến 10 Hz. Đây là khoảng tần số dao động tự nhiên phổ biến của hầu hết các công trình dân dụng và công nghiệp, rất dễ xảy ra hiện tượng cộng hưởng pha làm khuếch đại biên độ rung lắc, dẫn đến nứt gãy và phá hủy kết cấu chịu lực.

Tham số vận tốc sóng cắt VS30 đóng vai trò như thế nào trong mô hình? Tham số VS30 đo lường độ cứng của lớp đất mặt dày 30 m, phản ánh mức độ khuếch đại địa tầng cục bộ. Đất nền mềm có giá trị VS30 dưới 300 m/s có thể làm gia tăng biên độ gia tốc đỉnh lên hơn 2 lần so với nền đá gốc có VS30 trên 760 m/s.

Kết quả kiểm nghiệm của luận văn có thể mở rộng cho toàn bộ Việt Nam không? Kết quả kiểm nghiệm hoàn toàn phù hợp cho khu vực miền Bắc, đặc biệt là vùng Tây Bắc nơi chi phối bởi cơ chế đứt gãy trượt bằng nông. Đối với các vùng kiến tạo khác ở miền Trung hoặc thềm lục địa phía Nam, cần tiến hành kiểm nghiệm bổ sung với các tập dữ liệu địa chấn đặc thù trước khi áp dụng đại trà.

Kết luận

  • Luận văn đã thiết lập thành công quy trình xử lý và chuyển đổi tín hiệu sóng địa chấn từ băng ghi vận tốc dải rộng sang băng ghi gia tốc nền với hệ số tương quan đạt trên 0.90.
  • Khảo sát toàn diện 165 băng ghi từ 8 trận động đất độ lớn Ml từ 4.0 trở lên trong giai đoạn 2009 - 2014 tại miền Bắc Việt Nam.
  • Chứng minh mô hình suy giảm chấn động Campbell (2008) hoàn toàn phù hợp với điều kiện kiến tạo Tây Bắc với hơn 85% số liệu thực tế nằm trong khoảng sai số cho phép.
  • Đóng góp luận cứ khoa học quan trọng giúp giải quyết bài toán thiếu hụt dữ liệu gia tốc kế, phục vụ trực tiếp công tác thiết kế kháng chấn và đánh giá độ nguy hiểm động đất.
  • Lộ trình 2026 - 2030 cần tiếp tục mở rộng trạm quan trắc số hóa và cập nhật tiêu chuẩn kháng chấn quốc gia; các đơn vị tư vấn và chủ đầu tư cần chủ động ứng dụng mô hình Campbell (2008) vào thẩm định an toàn công trình ngay từ giai đoạn thiết kế cơ sở.