Tổng quan nghiên cứu

Việt Nam sở hữu đường bờ biển dài hơn 3.000 km trải dọc từ Bắc vào Nam, đi qua hơn 20 tỉnh, thành phố ven biển cùng hệ thống gần 3.000 hòn đảo lớn nhỏ. Khu vực ven biển là nơi tập trung hàng triệu cư dân sinh sống, đồng thời là địa bàn trọng điểm của nhiều công trình kinh tế, khu công nghiệp và hạ tầng an sinh xã hội trọng yếu. Tuy nhiên, dải ven biển nước ta luôn tiềm ẩn rủi ro chịu tác động trực tiếp từ hiểm họa thiên tai biển, đặc biệt là sóng thần do động đất kiến tạo kích hoạt. Theo cơ sở dữ liệu địa chấn thế giới, tính từ năm 1627 đến nay, đã có 62 trận sóng thần được ghi nhận trong khu vực Biển Đông và vùng biển lân cận phía tây Philippines.

Vấn đề cốt lõi đặt ra là sự thiếu hụt các nghiên cứu định lượng chuyên sâu về cơ chế lan truyền sóng thần từ các nguồn động đất xa và gần bờ đến thềm lục địa Việt Nam. Mục tiêu cụ thể của luận văn là làm sáng tỏ bản chất vật lý của quá trình kích hoạt, quy luật lan truyền sóng thần và xây dựng các kịch bản mô phỏng số chi tiết về thời gian truyền sóng cùng biên độ tác động đến bờ biển và hải đảo Việt Nam. Nghiên cứu tập trung vào không gian Biển Đông với nguồn chấn động từ đới hút chìm Manila (độ lớn cực đại ước tính Mw 8,85) và các đới đứt gãy hoạt động nội vùng, trên cơ sở chuỗi số liệu địa chấn kéo dài từ năm 1627 đến năm 2011. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn trong việc cung cấp các tham số định lượng chính xác, cho phép xác định cửa sổ thời gian cảnh báo sớm từ 100 đến 320 phút, hỗ trợ đắc lực cho công tác quy hoạch kháng chấn đô thị và giảm thiểu rủi ro thảm họa thiên tai ven biể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 được xây dựng trên nền tảng lý thuyết đàn hồi môi trường phân lớp và thủy động lực học sóng phi nén trong trường trọng lực. Trọng tâm của khung lý thuyết bao gồm mô hình nguồn đứt gãy đàn hồi Okada (1985), cho phép chuyển đổi tenxơ mômen địa chấn thành trường dịch chuyển tĩnh thẳng đứng và dịch chuyển ngang của đáy biển. Đồng thời, mô hình suy biến năng lượng và thang độ lớn mômen Kanamori với phương trình tương quan thực nghiệm log M0 = 1,5M + 9,1 được áp dụng nhằm định lượng chính xác năng lượng giải phóng từ các trận động đất kiến tạo có mômen vô hướng M0 = μ A u.

Các khái niệm then chốt được vận dụng xuyên suốt gồm: tenxơ mômen địa chấn 9 thành phần biểu diễn cặp ngẫu lực, nguồn dịch trượt nghịch (dip-slip), nguồn trượt bằng (strike-slip), góc cắm đứt gãy, góc trượt rake, hàm kích hoạt sóng thần R(ω) và phép biến đổi tích phân Fourier - Bessel kết hợp hàm Hankel bậc không H0^(2). Sự kích hoạt sóng trong tầng nước được mô hình hóa thông qua phương trình dao động chất lỏng kết hợp điều kiện biên bề mặt tự do và đáy biển đàn hồi, tạo cơ sở toán học chặt chẽ để xác định trường áp suất động và độ nâng mặt nước biển.

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

Luận văn sử dụng nguồn dữ liệu địa chấn và hải dương học đa nguồn, bao gồm danh mục 62 sự kiện sóng thần từ Trung tâm Dữ liệu Địa vật lý Quốc gia Hoa Kỳ (NGDC/NOAA), số liệu khảo sát vết tích sóng thần lịch sử từ bộ chính sử Đại Nam Thực lục và các báo cáo khảo sát thực địa ven biển miền Trung. Cỡ mẫu nghiên cứu bao gồm 6 kịch bản nguồn phát sinh động đất tiêu biểu tương ứng với 6 đới cấu trúc đứt gãy nguy hiểm nhất khu vực Biển Đông, với dải độ lớn động đất cực đại được lựa chọn phân tích từ 7,0 đến 8,85 độ Richter.

Phương pháp phân tích chủ đạo là thuật toán số hóa giải tích dựa trên hàm Green hai chiều (2D Green's Function approach) kết hợp giải tích tiệm cận Levshin. Lý do lựa chọn phương pháp hàm Green thay vì các mô hình sai phân hữu hạn thuần túy là khả năng tính toán nhanh vượt trội, biểu diễn chính xác phản ứng động lực học của cột nước dưới tác động của trường lực nguồn đa chiều và mô tả hiệu quả sự tương tác giữa sóng biển với đáy đại dương biến đổi mà không đòi hỏi tài nguyên tính toán quá lớn. Quá trình mô phỏng được thực hiện trên cấu trúc hai lớp gồm lớp chất lỏng hữu hạn trên nửa không gian rắn đàn hồi đồng nhất.

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

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

Thứ nhất, phân tích cơ chế kích hoạt chỉ ra rằng các đứt gãy chờm nghịch (dip-slip) có góc cắm 45 độ tạo ra độ nâng đáy biển và biên độ kích hoạt sóng thần vượt trội hoàn toàn so với đứt gãy trượt bằng (strike-slip). Tại khoảng cách chấn tâm 500 km, thành phần dịch chuyển tỏa tia theo phương ngang của phần tử nước lớn hơn khoảng 3 lần so với thành phần chuyển động thẳng đứng, hình thành quỹ đạo chuyển động dạng elip đặc trưng của sóng nước nông.

Thứ hai, kết quả mô phỏng kịch bản động đất cực đại Mw 8,85 tại đới hút chìm Manila cho thấy thời gian sóng thần lan truyền đến dải bờ biển Quảng Ngãi - Bình Định chỉ mất khoảng 100 phút (khoảng 1 giờ 40 phút), và sau 320 phút (khoảng 5 giờ 20 phút), sóng sẽ lan tỏa đến bờ biển Minh Hải đồng thời đi sâu vào khu vực Vịnh Bắc Bộ. Biên độ sóng ghi nhận mức suy giảm dần theo khoảng cách lan truyền nhưng có xu hướng gia tăng đột biến khi tiến vào vùng nước nông ven bờ.

Thứ ba, nghiên cứu đã định lượng hóa mức độ nguy hiểm của 5 đới đứt gãy nội vùng trong Biển Đông (gồm đứt gãy Bắc Hoàng Sa, Kinh tuyến 110 độ, Thuận Hải - Minh Hải, Cảnh Dương - Phú Quý và Palawan) với độ lớn động đất cực đại từ 7,0 đến 8,0 độ Richter. Mặc dù biên độ nguồn thấp hơn đới Manila khoảng 40% đến 60%, các nguồn cận bờ này có thời gian truyền sóng cực ngắn (dưới 30 phút), tạo ra nguy cơ sóng thần địa phương đặc biệt nghiêm trọng.

Thứ tư, khi phân tích cấu trúc địa tầng đáy biển, kết quả so sánh giữa 3 mô hình cấu trúc cho thấy lớp trầm tích bở rời có tác động không đáng kể đến phổ kích hoạt sóng, trong khi độ sâu chấn tiêu nằm trong khoảng 9 đến 14 km thuộc phần trên vỏ Trái Đất là tầng phát sinh sóng thần mạnh nhất. Mô hình đáy biển mấp mô (2D) làm tăng biên độ sóng tổng hợp và gây hiện tượng trễ sóng rõ rệt so với mô hình đáy phẳng (1D).

Thảo luận kết quả

Sự gia tăng biên độ sóng khi truyền từ biển sâu vào đới bờ thoải giải thích rõ hiện tượng bảo toàn thông lượng năng lượng: khi vận tốc truyền sóng giảm theo độ sâu lớp nước, năng lượng sóng bị dồn nén làm chiều cao sóng tăng vọt theo định lý Green. Kết quả này tương thích với các quan sát thực địa về các đợt sóng dâng lịch sử tại Nghệ An và Bình Thuận ghi nhận trong tài liệu cổ.

Trong nghiên cứu, dữ liệu lan truyền sóng được minh họa trực quan qua hệ thống sơ đồ đường đẳng thời truyền sóng từ chấn tâm Manila tỏa khắp Biển Đông và các biểu đồ chuỗi thời gian dao động triều tổng hợp. So sánh với các mô hình thủy động lực quốc tế trong cuộc diễn tập Pacific Wave 06, kết quả mô phỏng bằng hàm Green cho độ sai lệch thời gian truyền sóng dưới 5%, đồng thời phản ánh rõ nét sự tán sắc pha của chuỗi sóng tại các trạm quan sát giả định trên thềm lục địa.

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

Thứ nhất, khẩn trương thiết lập hệ thống cảnh báo sớm sóng thần đa cấp độ tại khu vực miền Trung, do Viện Vật lý Địa cầu chủ trì phối hợp cùng Tổng cục Khí tượng Thủy văn. Mục tiêu đặt ra là rút ngắn thời gian phát bản tin cảnh báo động đất gây sóng thần xuống dưới 15 phút sau khi ghi nhận chấn động từ đới Manila, hoàn thành lắp đặt hệ thống còi báo động và trạm quan trắc mực nước biển chuyên dụng giai đoạn 2026 - 2028.

Thứ hai, tích hợp các kịch bản ngập lụt sóng thần vào quy hoạch xây dựng và phát triển đô thị ven biển tại 28 tỉnh, thành phố trực thuộc Trung ương. Ủy ban nhân dân các địa phương phối hợp với Bộ Xây dựng rà soát, áp dụng tiêu chuẩn kháng chấn và tải trọng sóng xung kích cho 100% công trình công cộng, bệnh viện, trường học nằm trong bán kính 1.000 mét tính từ mép nước biển, thực hiện xuyên suốt giai đoạn 2026 - 2030.

Thứ ba, củng cố và nâng cấp hệ thống đê biển, kè chắn sóng tại các cửa sông và dải bờ xung yếu chịu tác động trực tiếp từ đới hút chìm Manila. Bộ Nông nghiệp và Phát triển nông thôn cần xây dựng quy chuẩn kỹ thuật đê kè có khả năng chống chịu mức sóng dâng cực đại từ 3 đến 5 mét, kết hợp phát triển vành đai rừng ngập mặn phòng hộ ven biển đạt độ che phủ trên 70% diện tích bãi bồi vào năm 2030.

Thứ tư, định kỳ tổ chức diễn tập sơ tán khẩn cấp cấp cộng đồng cho cư dân ven biển và ngư dân hoạt động trên các ngư trường trọng điểm. Ủy ban Quốc gia Ứng phó Sự cố, Thiên tai và Tìm kiếm Cứu nạn cần ban hành kịch bản sơ tán chi tiết theo mốc thời gian vàng (từ 60 đến 90 phút trước khi sóng ập bờ), bảo đảm 95% người dân trong vùng nguy cơ nắm vững lộ trình di chuyển lên các khu vực cao an toàn từ năm 2026 trở đi.

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

Thứ nhất, các nhà khoa học, nghiên cứu viên và học viên cao học thuộc các chuyên ngành Vật lý Địa cầu, Hải dương học, Địa chất công trình và Quản lý Tai biến Thiên tai. Luận văn cung cấp khung giải tích toán học chuẩn mực về hàm Green và mô hình biến dạng đàn hồi Okada, đóng vai trò là tài liệu tham khảo chất lượng cho các nghiên cứu mô phỏng sóng biển và động lực học chất lỏng địa vật lý.

Thứ hai, kỹ sư công trình biển, chuyên gia tư vấn thiết kế và các cơ quan quy hoạch hạ tầng ven biển. Tài liệu mang lại dữ liệu nền tảng về trường ứng suất sóng, độ dịch chuyển đáy biển và thời gian truyền sóng, phục vụ thiết kế đê kè cảng biển, công trình dầu khí ngoài khơi và các dự án năng lượng tái tạo ven bờ tại thềm lục địa Việt Nam.

Thứ ba, cán bộ quản lý thuộc Cục Quản lý Đê điều và Phòng, chống Thiên tai, Trung tâm Báo tin Động đất và Cảnh báo Sóng thần. Nhóm đối tượng này có thể sử dụng trực tiếp các kịch bản lan truyền sóng và mốc thời gian truyền 100 - 320 phút để hoàn thiện quy trình phản ứng nhanh và lập phương án ứng phó khẩn cấp quốc gia.

Thứ tư, các doanh nghiệp vận hành chuỗi cung ứng logistics cảng biển, khu nghỉ dưỡng ven biển và đơn vị khai thác vận tải biển. Nghiên cứu giúp các tổ chức này đánh giá chính xác rủi ro kinh doanh, xây dựng kế hoạch quản trị rủi ro thiên tai và bảo toàn tài sản trước các biến cố sóng thần cực đoan.

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

Bờ biển Việt Nam có thực sự đối mặt với nguy cơ sóng thần không? Hoàn toàn có. Lịch sử địa chấn ghi nhận 62 trận sóng thần trong khu vực Biển Đông và lân cận từ năm 1627. Dù được bao bọc bởi các vòng cung đảo ngăn sóng từ Thái Bình Dương, bờ biển nước ta vẫn chịu nguy cơ trực tiếp từ đới hút chìm Manila và các đới đứt gãy nội vùng có khả năng gây sóng cao từ 2 đến 5 mét.

Tại sao đới hút chìm Manila lại là nguồn phát sinh sóng thần nguy hiểm nhất đối với Việt Nam? Đới hút chìm Manila dài trên 1.000 km, sở hữu cơ chế trượt chờm nghịch mạnh mẽ với khả năng tích tụ năng lượng gây động đất cực đại Mw 8,85. Vị trí đối diện trực diện với miền Trung Việt Nam khiến năng lượng bức xạ sóng thần truyền thẳng qua biển sâu mà không bị suy giảm đáng kể bởi vật cản địa hình.

Phương pháp hàm Green có ưu điểm gì vượt trội so với mô hình số sai phân truyền thống? Phương pháp hàm Green kết hợp phép biến đổi giải tích giúp rút ngắn thời gian xử lý dữ liệu tính toán tới hàng chục lần so với phương pháp sai phân lưới hữu hạn truyền thống. Phương pháp này cho phép phân tách chính xác đáp ứng động lực học của cột nước từ nguồn điểm lực sang nguồn hữu hạn phức tạp.

Thời gian cảnh báo sóng thần lan truyền từ đới Manila đến bờ biển Việt Nam là bao lâu? Theo tính toán mô phỏng trong luận văn, sóng thần mất khoảng 100 phút để di chuyển từ trung tâm đới Manila đến bờ biển Quảng Ngãi - Bình Định. Đối với các khu vực xa hơn như bờ biển Minh Hải và Vịnh Bắc Bộ, tổng thời gian truyền sóng kéo dài khoảng 320 phút, mở ra khung thời gian quý giá cho công tác sơ tán.

Các nguồn động đất cục bộ trong Biển Đông nguy hiểm như thế nào so với nguồn Manila? Các đứt gãy nội vùng như Kinh tuyến 110 độ hay Thuận Hải - Minh Hải có độ lớn động đất tối đa khoảng 7,0 đến 7,5 Richter, tạo biên độ sóng nhỏ hơn nguồn Manila. Tuy nhiên, do khoảng cách rất gần bờ, thời gian sóng ập vào đất liền chỉ từ 15 đến 30 phút, khiến việc phản ứng và sơ tán trở nên cấp bách hơn nhiều.

Kết luận

Luận văn đã giải quyết trọn vẹn bài toán mô phỏng lan truyền sóng thần bằng phương pháp giải tích hàm Green với những đóng góp khoa học và thực tiễn then chốt:

  • Xây dựng thành công bộ công cụ giải tích và thuật toán mô phỏng 2D dựa trên hàm Green cho bài toán lan truyền sóng thần trong môi trường chất lỏng trên nền đàn hồi.
  • Định lượng chính xác trường dịch chuyển đáy biển và khẳng định cơ chế đứt gãy trượt nghịch tạo ra biên độ kích hoạt sóng thần nguy hiểm nhất.
  • Thiết lập 6 kịch bản lan truyền sóng chi tiết từ các đới nguồn trọng yếu, xác định mốc thời gian cảnh báo từ 100 phút đến 320 phút đối với nguồn Manila.
  • Đánh giá toàn diện ảnh hưởng của địa hình đáy biển mấp mô và độ sâu chấn tiêu vỏ Trái Đất (9 - 14 km) đến hiện tượng tán sắc và khuếch đại biên độ sóng.
  • Đóng góp luận cứ khoa học tin cậy cho việc xây dựng bản đồ cảnh báo nguy cơ sóng thần và quy hoạch hạ tầng kháng chấn ven biển Việt Nam.

Trong giai đoạn 2026 - 2030, các cơ quan quản lý và giới học thuật cần tiếp tục mở rộng mô hình hóa 3D cho các nguồn vùng hữu hạn phi đối xứng. Hãy chủ động tham khảo và ứng dụng các kết quả nghiên cứu này vào công tác lập quy hoạch phát triển bền vững kinh tế biển và nâng cao năng lực ứng phó thiên tai quốc gia ngay hôm nay.