Tổng quan nghiên cứu

Thảm họa thiên tai ven biển luôn là mối đe dọa thường trực đối với các quốc gia có đường bờ biển dài. Lịch sử nhân loại đã chứng kiến những tổn thất khốc liệt với hơn 230.000 người thiệt mạng trong trận sóng thần Ấn Độ Dương năm 2004 và thiệt hại kinh tế vượt mốc 210 tỷ USD trong thảm họa động đất - sóng thần Tohoku năm 2011 tại Nhật Bản. Nằm ở rìa phía tây của Biển Đông với đường bờ biển trải dài hơn 3.260 km cùng hơn 28 tỉnh, thành phố tập trung đông đảo dân cư và cơ sở hạ tầng trọng điểm, Việt Nam đứng trước những nguy cơ tiềm tàng về hiểm họa sóng thần xuyên đại dương. Biển Đông là vùng biển nửa kín sở hữu cấu trúc địa chất phức tạp, đặc biệt là sự hiện diện của rãnh hút chìm Manila dài hơn 1.150 km với khả năng phát sinh các trận động đất có độ lớn từ 8.0 đến 9.0 độ Richter.

Nhận thức rõ thách thức này, nghiên cứu được thực hiện với mục tiêu trọng tâm là tính toán, mô phỏng quá trình thành tạo và lan truyền sóng thần từ các đới nguồn tiềm năng trên Biển Đông tới vùng bờ biển Việt Nam. Luận văn hướng đến việc xác định cụ thể hai thông số thủy động lực học cốt lõi: phân bố chiều cao sóng lớn nhất và thời gian lan truyền sóng từ nguồn tới bờ. Phạm vi không gian nghiên cứu bao phủ toàn bộ vùng trũng Biển Đông và dải bờ biển nước ta, với dữ liệu nguồn gồm 13 kịch bản động đất tiêu chuẩn do Bộ Tài nguyên và Môi trường đề xuất, được thực hiện và hoàn thành trong giai đoạn từ tháng 6/2013 đến tháng 12/2013. Kết quả nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc, cung cấp hệ cơ sở dữ liệu số trị đáng tin cậy phục vụ công tác xây dựng bản đồ cảnh báo nguy cơ sóng thần quốc gia, hỗ trợ quy hoạch các dải đê kè và công trình ven bờ với độ an toàn cao trình từ 5m đến 14m, trực tiếp góp phần bảo vệ tính mạng và tài sản cho hơn 20 triệu cư dân vùng 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 dựa trên sự kết hợp chặt chẽ giữa lý thuyết cơ học môi trường liên tục trong vỏ Trái Đất và lý thuyết thủy động lực học hải dương. Để mô tả quá trình khởi phát sóng thần, khung lý thuyết biến dạng đàn hồi do chuyển dịch phay đứt gãy của Mansinha và Smylie (1971) được áp dụng dựa trên công thức tích phân mặt Volterra tổng quát và các mở rộng của Okada (1985). Lý thuyết này thiết lập mối quan hệ định lượng giữa sự trượt tức thời của phay kiến tạo (bao gồm trượt bằng strike-slip và trượt nghiêng dip-slip) với chuyển vị thẳng đứng của đáy biển, tạo ra xung lực ban đầu làm dâng hạ mực nước đại dương.

Để mô phỏng quá trình lan truyền sóng trên quy mô đại dương, luận văn sử dụng hệ phương trình sóng nước nông phi tuyến hai chiều (Nonlinear Shallow Water Equations - NSWE) tích phân theo độ sâu. Khung lý thuyết này vận hành dựa trên định luật bảo toàn khối lượng và bảo toàn động lượng, trong đó gia tốc phương đứng của hạt nước được xem là rất nhỏ so với gia tốc trọng trường 9.81 m/s², cho phép áp dụng gần đúng phân bố áp suất thủy tĩnh. Bốn khái niệm chuyên ngành then chốt được tích hợp xuyên suốt mô hình gồm:

  • Biến dạng đồng chấn đáy biển (Coseismic seabed deformation): Sự chuyển dịch tức thời của lớp vỏ đáy biển dưới tác động của trường ứng suất kiến tạo.
  • Thông lượng xung lượng mặt cắt nước (Momentum flux): Tích phân vận tốc dòng chảy theo phương đứng nhân với tổng độ sâu nước động.
  • Hiệu ứng nước nông và biến đổi vận tốc sóng (Shoaling effect): Hiện tượng vận tốc truyền sóng suy giảm theo độ sâu địa hình đáy theo quy luật căn bậc hai của tích số gia tốc trọng trường và độ sâu tĩnh, dẫn đến sự dồn tụ năng lượng.
  • Chiều cao sóng lớn nhất ven bờ (Maximum tsunami run-up amplitude): Dao động mực nước cực đại trên mực nước tĩnh khi sóng tiến sát đới bờ.

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

Nghiên cứu khai thác nguồn dữ liệu đầu vào chuẩn hóa bao gồm tập hợp 13 kịch bản động đất nguồn do Bộ Tài nguyên và Môi trường ban hành, kết hợp dữ liệu địa hình độ sâu đáy biển ETOPO2 với độ phân giải không gian 2 phút góc (khoảng 3.7 km) từ Cơ quan Quản lý Khí quyển và Đại dương Quốc gia Hoa Kỳ (NOAA) và dữ liệu từ Ủy ban Hải dương học Liên chính phủ (IOC/UNESCO).

Về cỡ mẫu và phương pháp chọn mẫu, toàn bộ 13 kịch bản được phân bổ toàn diện tại 4 vùng nguồn kiến tạo có hoạt động địa chấn mạnh nhất Biển Đông:

  • Đới hút chìm Manila (5 kịch bản): Độ lớn động đất từ 7.0 đến 9.0 độ Richter, chiều dài phay từ 35 km đến 646 km, độ trượt phay từ 1.6m đến 17.49m.
  • Đới hút chìm Ryukyu - Nam Đài Loan (5 kịch bản): Độ lớn từ 7.0 đến 9.0 độ Richter, chiều dài phay từ 35 km đến 646 km.
  • Đới đứt gãy Nam Trung Bộ (1 kịch bản): Độ lớn 7.0 độ Richter, chiều dài phay 50 km, độ sâu tiêu cự 10 km.
  • Đới đứt gãy Nam Hải Nam (2 kịch bản): Độ lớn 7.0 và 7.5 độ Richter, chiều dài phay lần lượt là 50 km và 89 km.

Phương pháp phân tích số trị được lựa chọn là phương pháp sai phân hữu hạn hiện kiểu nhảy cóc (Leap-frog finite difference scheme) trên lưới so le không gian và thời gian. Lý do lựa chọn sơ đồ Leap-frog là khả năng bảo toàn năng lượng sóng vượt trội, triệt tiêu sai số bậc nhất để đạt độ chính xác bậc hai O(Δx², Δt²), đồng thời kết hợp kỹ thuật sai phân ngược dòng (upwind scheme) cho các số hạng đối lưu phi tuyến để ngăn chặn hiện tượng mất ổn định số học và khuếch tán giả khi sóng truyền qua hàng nghìn kilômét đại dương. Toàn bộ chu trình tính toán, thiết lập thuật toán Fortran và xử lý kết quả được tiến hành liên tục trong thời gian 6 tháng từ tháng 6/2013 đến tháng 12/2013.

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

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

Phân tích số trị từ 13 kịch bản mô phỏng đã mang lại những kết quả quan trọng về trường sóng và thời gian truyền sóng tại bờ biển Việt Nam:

Thứ nhất, rãnh hút chìm Manila là vùng nguồn nguy hiểm nhất, có khả năng gây thiệt hại nặng nề cho toàn bộ dải duyên hải miền Trung. Tại kịch bản 4 (động đất Manila 8.5 độ Richter), sóng thần xuất hiện trên diện rộng với chiều cao sóng vượt quá 1m tại vùng biển từ Vũng Tàu đến Nam Định, vượt 2m tại khu vực từ Đà Nẵng đến Bắc Phan Thiết, và đặc biệt đạt đỉnh vượt quá 6m tại bờ biển tỉnh Quảng Ngãi.

Thứ hai, trong kịch bản cực đoan 5 (động đất Manila 9.0 độ Richter, chiều dài phay 646 km, trượt 17.49m), toàn bộ dải ven biển nước ta phải đối mặt với thảm họa nghiêm trọng. Chiều cao sóng tại Quảng Ngãi chạm ngưỡng trên 14m (tăng hơn 133% so với kịch bản 4), khu vực Quảng Trị đến Phan Thiết ghi nhận sóng cao trên 2m, và ngay cả khu vực vịnh kín như Hải Phòng cũng ghi nhận mực nước dâng cao hơn 1m.

Thứ ba, thời gian lan truyền sóng từ đới Manila đến đất liền Việt Nam diễn ra với tốc độ rất cao. Kết quả mô phỏng cho thấy sóng thần di chuyển với vận tốc trung bình khoảng 700 - 800 km/h ngoài khơi sâu, chỉ mất dưới 60 phút để đánh vào quần đảo Hoàng Sa và Trường Sa, khoảng 120 phút để tiến sát vùng bờ biển Bình Định - Phan Rang, và khoảng 180 phút để lan truyền đến dải ven biển Đà Nẵng - Quảng Ngãi.

Thứ tư, các nguồn chấn động nội tại như đới đứt gãy Nam Trung Bộ (kịch bản 11) và Nam Hải Nam (kịch bản 12 và 13) tạo ra chiều cao sóng ven bờ tương đối nhỏ, dao động dưới 1m tại Đà Nẵng và lân cận. Tuy nhiên, do khoảng cách từ nguồn đến đất liền ngắn (dưới 200 km), thời gian cảnh báo chỉ còn khoảng 120 phút dù vận tốc sóng bị suy giảm khi di chuyển qua thềm lục địa nông.

Thảo luận kết quả

Nguyên nhân cốt lõi khiến đới hút chìm Manila gây ra trường sóng tàn phá lớn nhất bắt nguồn từ đặc điểm hình học phay và hướng phát tán năng lượng. Phay Manila có góc dốc 15 độ, chiều sâu chấn tiêu nông từ 5 km đến 27 km, và phương vị phay gần như hướng Bắc - Nam (177 độ). Hướng giải phóng động năng chính vuông góc với trục phay, nghĩa là hướng thẳng trực diện về phía tây - đúng vào dải bờ biển miền Trung Việt Nam. Khi sóng di chuyển từ bồn trũng sâu trên 4.000m vào vùng nước nông ven bờ miền Trung, theo định luật Green, vận tốc sóng giảm mạnh nhưng năng lượng được bảo toàn dẫn đến biên độ sóng tăng vọt theo phương đứng.

Dữ liệu mô phỏng trong luận văn được minh họa rõ nét thông qua hai định dạng trực quan:

  • Biểu đồ chuỗi thời gian dao động mực nước (time-series) tại các trạm kiểm chứng ven biển trọng điểm (Đà Nẵng, Quảng Ngãi, Phú Yên), thể hiện rõ các pha dâng hạ ban đầu, chu kỳ sóng kéo dài từ 20 đến 40 phút và thời điểm xuất hiện đỉnh sóng cao nhất.
  • Bản đồ đẳng thời gian lan truyền (tsunami travel time map) với các đường đồng mức thời gian giãn cách 30 phút, phản ánh sự khúc xạ trường sóng qua địa hình phức tạp của quần đảo Hoàng Sa và bồn trũng Biển Đông.

So sánh với các công trình nghiên cứu của các tác giả trong và ngoài nước sử dụng mô hình MOST hoặc phương trình phi tuyến nước nông (như nghiên cứu của Vũ Thanh Ca năm 2008 và Trần Thị Mỹ Thành năm 2011), kết quả của luận văn thể hiện sự nhất quán cao về độ lớn chiều cao sóng cực đại (sai số dưới 8%) và thời gian truyền sóng. Điều này khẳng định tính đúng đắn, độ tin cậy của thuật toán sai phân Leap-frog được thiết lập trong mô hình.

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

Dựa trên các phát hiện định lượng về độ cao sóng và thời gian truyền sóng, luận văn đề xuất 4 nhóm giải pháp chiến lược nhằm giảm nhẹ rủi ro thiên tai:

  • Nâng cấp và mở rộng mạng lưới trạm cảnh báo sớm sóng thần đa tầng: Tiến hành khảo sát, lắp đặt bổ sung 15 trạm cảnh báo âm học chuyên dụng kết hợp cảm biến đo áp suất đáy biển tại các tỉnh ven biển miền Trung (từ Quảng Trị đến Bình Thuận), kết nối đồng bộ với 2 trạm cảnh báo hiện hữu tại Mân Thái (16,094°B; 108,251°Đ) và Xuân Thiều (16,098°B; 108,140°Đ) thuộc thành phố Đà Nẵng. Mục tiêu: Rút ngắn thời gian phát tin cảnh báo tự động xuống dưới 5 phút sau khi ghi nhận chấn động địa chấn. Thời gian thực hiện: Giai đoạn 2026 - 2028. Chủ thể chủ trì: Viện Vật lý Địa cầu phối hợp cùng Tổng cục Khí tượng Thủy văn.

  • Xây dựng bản đồ phân vùng ngập lụt chi tiết và cắm mốc cảnh báo ngập sâu: Triển khai thành lập hệ thống bản đồ ngập lụt do sóng thần tỷ lệ 1:10.000 cho 100% các đô thị và khu kinh tế ven biển trọng điểm (Đà Nẵng, Hội An, Quy Nhơn, Tuy Hòa, Nha Trang), xác định rõ ranh giới ngập sâu từ 2m đến trên 14m ứng với kịch bản động đất Manila cấp 8.5 đến 9.0. Thời gian thực hiện: Giai đoạn 2026 - 2029. Chủ thể chủ trì: Bộ Tài nguyên và Môi trường phối hợp cùng Ủy ban nhân dân các tỉnh duyên hải miền Trung.

  • Ban hành quy chuẩn thiết kế và gia cố công trình hạ tầng ven bờ: Thiết lập quy chuẩn kỹ thuật xây dựng bắt buộc các công trình công cộng, khu tái định cư ven biển phải có cao trình sàn vượt mức ngập tính toán tối thiểu 5m, đồng thời quy hoạch các tòa nhà kiên cố nhiều tầng làm điểm sơ tán thẳng đứng tại chỗ; kết hợp trồng dải rừng phòng hộ ven biển rộng tối thiểu 200m để tiêu tán năng lượng sóng. Thời gian thực hiện: Ban hành quy chuẩn trong quý 4 năm 2027. Chủ thể chủ trì: Bộ Xây dựng phối hợp với Bộ Nông nghiệp và Phát triển nông thôn.

  • Tổ chức diễn tập ứng phó cộng đồng định kỳ và chuẩn hóa quy trình sơ tán khẩn cấp: Định kỳ hàng năm tổ chức các đợt diễn tập thực địa quy mô cấp tỉnh tại các địa bàn xung yếu, phấn đấu đạt tỷ lệ 95% người dân trong vùng nguy cơ nắm vững lộ trình di chuyển tới vị trí an toàn cao trên 15m trong khung thời gian vàng từ 15 đến 30 phút. Thời gian thực hiện: Triển khai thường niên từ năm 2026. Chủ thể chủ trì: Ban Chỉ huy Phòng chống thiên tai và Tìm kiếm cứu nạn các địa phương.

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

Công trình nghiên cứu mang giá trị thực tiễn và hàm lượng học thuật cao, là tài liệu tham khảo thiết thực cho 4 nhóm đối tượng chính:

  • Cơ quan quản lý nhà nước và hoạch định chính sách phòng chống thiên tai (Bộ Tài nguyên và Môi trường, Ban Chỉ đạo Quốc gia về Phòng chống thiên tai): Luận văn cung cấp luận cứ khoa học chuẩn xác để xây dựng kịch bản ứng phó quốc gia và phân bổ nguồn lực công trình phòng thủ bờ biển. Trường hợp sử dụng tiêu biểu: Ứng dụng số liệu 13 kịch bản để quy hoạch hành lang an toàn bờ biển và vận hành hệ thống cảnh báo sớm trực ban 24/7.

  • Đơn vị tư vấn thiết kế, nhà thầu xây dựng công trình thủy công và công trình biển: Các kỹ sư tiếp cận được dữ liệu tải trọng sóng động lực học, cao trình mực nước dâng cực đại từ 1m đến 14m để tính toán ổn định công trình. Trường hợp sử dụng tiêu biểu: Thiết kế đê kè biển, bến cảng nước sâu, nhà máy lọc hóa dầu và các khu nghỉ dưỡng ven biển miền Trung đáp ứng tiêu chuẩn kháng chấn và chống tràn sóng.

  • Các viện nghiên cứu, trường đại học và nghiên cứu sinh chuyên ngành Hải dương học, Địa vật lý, Kỹ thuật Công trình biển: Tài liệu cung cấp phương pháp luận toán học chi tiết về thuật toán sai phân Leap-frog và kỹ thuật mô phỏng nguồn đứt gãy. Trường hợp sử dụng tiêu biểu: Phát triển các mô hình số trị thế hệ mới mô phỏng tương tác sóng - dòng chảy và sóng tràn qua công trình.

  • Các tập đoàn bất động sản nghỉ dưỡng và doanh nghiệp khai thác kinh tế biển: Giúp các nhà đầu tư nhận diện chính xác mức độ rủi ro thiên tai tại từng tọa độ dự án ven bờ. Trường hợp sử dụng tiêu biểu: Đánh giá tác động môi trường, tối ưu hóa quy hoạch mặt bằng tổng thể và thiết lập phương án bảo hiểm tài sản trước các kịch bản thiên tai 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ơ xảy ra thảm họa sóng thần hay không? Bờ biển Việt Nam hoàn toàn có nguy cơ hứng chịu thảm họa sóng thần nguy hiểm, chủ yếu xuất phát từ đới hút chìm rãnh Manila dài hơn 1.150 km nằm ở phía đông Biển Đông. Kết quả tính toán của luận văn chứng minh các trận động đất có độ lớn từ 7.5 đến 9.0 độ Richter tại đới nguồn này đủ khả năng kích hoạt chuỗi sóng thần vượt đại dương, tạo ra sóng dâng cao từ 2m đến 14m tấn công dải ven biển nước ta.

  • Vùng bờ biển nào tại Việt Nam chịu mức độ rủi ro sóng thần cao nhất? Khu vực duyên hải miền Trung, đặc biệt là dải bờ biển từ thành phố Đà Nẵng qua tỉnh Quảng Ngãi đến tỉnh Bình Định, là nơi chịu rủi ro nặng nề nhất. Theo mô phỏng kịch bản động đất Manila cấp 9.0, khu vực Quảng Ngãi ghi nhận chiều cao sóng cực đại vượt quá 14m do nằm trực diện trên trục bức xạ năng lượng chính của đới đứt gãy kiến tạo.

  • Thời gian để chính quyền và người dân triển khai sơ tán khẩn cấp là bao lâu? Thời gian lan truyền sóng từ rãnh Manila đến quần đảo Hoàng Sa là dưới 60 phút, đến vùng biển Bình Định - Phan Rang khoảng 120 phút và đến Đà Nẵng khoảng 180 phút. Do đó, các địa phương ven biển miền Trung có quỹ thời gian vàng quý giá từ 2 đến 3 giờ để kích hoạt còi báo động, phát tin trên hệ thống truyền thông và sơ tán người dân đến các vị trí an toàn.

  • Ưu thế vượt trội của phương pháp sai phân hữu hạn Leap-frog trong mô phỏng sóng thần là gì? Phương pháp sai phân Leap-frog trên lưới so le cho phép xấp xỉ đạo hàm bậc hai với sai số nhỏ O(Δx², Δt²), đảm bảo độ ổn định cao và triệt tiêu sai số làm tròn khi mô phỏng trên miền tính toán hàng triệu km² của Biển Đông. Phương pháp này cân bằng hoàn hảo giữa độ chính xác động lực học và hiệu quả tối ưu thời gian tính toán của máy tính.

  • Hiện nay Việt Nam đã có những công trình thực nghiệm nào để phòng chống và cảnh báo sóng thần? Việt Nam đã đưa vào vận hành thử nghiệm 2 trạm cảnh báo sóng thần chuyên dụng đầu tiên tại thành phố Đà Nẵng, gồm trạm Mân Thái (tọa độ 16,094°B; 108,251°Đ) và trạm Xuân Thiều (tọa độ 16,098°B; 108,140°Đ). Song song đó, các cơ quan chức năng đã tổ chức thành công nhiều đợt diễn tập thực địa ứng phó sóng thần quy mô lớn tại Đà Nẵng, Quảng Ngãi và Bình Định.

Kết luận

  • Luận văn đã xây dựng thành công mô hình số trị mô phỏng sự thành tạo và lan truyền sóng thần trên Biển Đông bằng phương pháp sai phân hữu hạn Leap-frog kết hợp phương trình sóng nước nông phi tuyến.
  • Xác lập cơ sở khoa học khẳng định đới hút chìm Manila là vùng nguồn nguy hiểm nhất, có thể tạo ra sóng thần cao vượt quá 14m tại Quảng Ngãi và trên 2m tại dải bờ biển từ Quảng Trị đến Phan Thiết trong kịch bản động đất cấp 9.0 độ Richter.
  • Cung cấp bản đồ thời gian truyền sóng chi tiết, xác định khung thời gian vàng sơ tán khẩn cấp cho vùng duyên hải miền Trung dao động từ 2 đến 3 giờ sau chấn động.
  • Kế hoạch hành động giai đoạn 2026 - 2030 cần tập trung hoàn thiện hệ thống quan trắc cảnh báo sớm đa điểm và tích hợp bản đồ phân vùng ngập lụt tỷ lệ 1:10.000 vào chiến lược phát triển kinh tế biển bền vững.
  • Các nhà quản lý, đơn vị tư vấn kỹ thuật và doanh nghiệp ven biển cần chủ động khai thác các thông số định lượng từ luận văn để hoàn thiện phương án phòng chống thiên tai và bảo vệ an toàn công trình ngay hôm nay.