Tổng quan nghiên cứu

Các hiểm họa địa chấn ven biển luôn đặt ra thách thức nghiêm trọng đối với độ an toàn của hạ tầng hàng hải. Tại Việt Nam, các trận động đất ghi nhận vào năm 2010 với cường độ đạt 5,0 độ Richter, cùng chuỗi động đất kích thích liên tiếp tại khu vực miền Trung năm 2012, đã gióng lên hồi chuông cảnh báo về an toàn công trình đập và đê kè ven biển. Đê chắn sóng cảng Dung Quất thuộc tỉnh Quảng Ngãi là công trình bảo vệ bờ biển trọng điểm quốc gia, được xây dựng trên tầng địa chất cát bão hòa nước ngập mặn 100% dưới đáy biển. Khi xảy ra rung chấn địa chấn, hiện tượng tăng áp lực nước lỗ rỗng có thể triệt tiêu toàn bộ ứng suất hữu hiệu của nền đất, dẫn tới hiện tượng hóa lỏng cục bộ và làm mất khả năng chịu lực của mái đê.

Mục tiêu cốt lõi của nghiên cứu là phân tích, đánh giá toàn diện độ ổn định trượt mái và nguy cơ hóa lỏng nền đê chắn sóng Dung Quất dưới tác động của tải trọng động đất. Luận văn tiến hành so sánh đối chiếu giữa các công thức bán thực nghiệm quốc tế và phương pháp mô phỏng số động lực học tiên tiến, đồng thời chỉ ra những sai lệch tiềm ẩn của phương pháp tĩnh lực truyền thống. Phạm vi nghiên cứu tập trung vào mặt cắt địa chất đại diện với chiều sâu khảo sát từ 15m đến 40m trong giai đoạn hoàn thiện thiết kế công trình. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn khi giúp các đơn vị quản lý hạ tầng hàng hải tối ưu hóa giải pháp gia cố nền móng, giảm thiểu nguy cơ sụt lún công trình lên tới 25% đến 35%, đồng thời tiết kiệm hàng chục tỷ đồng chi phí bảo trì và khắc phục sự cố sau thiên tai.

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 ứng dụng hai hệ thống lý thuyết thực nghiệm và bán thực nghiệm kinh điển trong lĩnh vực địa kỹ thuật động:

  • Lý thuyết đánh giá khả năng hóa lỏng của Seed và cộng sự (1983): Xác định tiềm năng hóa lỏng thông qua tương quan giữa ứng suất cắt tuần hoàn do động đất sinh ra và cường độ kháng lặp của nền đất, chuẩn hóa theo giá trị thí nghiệm xuyên tiêu chuẩn SPT với 60% năng lượng búa hiệu chỉnh.
  • Mô hình bán thực nghiệm của Idriss và Boulanger (2004): Bổ sung các hàm phi tuyến biểu diễn hệ số giảm ứng suất theo chiều sâu, hệ số quy đổi cấp động đất và hệ số hiệu chỉnh ứng suất hữu hiệu ban đầu đối với cát rời bão hòa.

Bên cạnh đó, luận văn khai thác lý thuyết trạng thái tới hạn Cam-clay và mô hình ứng suất phi tuyến hyperbolic của Duncan - Chang nhằm mô tả chính xác mối quan hệ giữa ứng suất và biến dạng của vật liệu đắp thân đê cũng như các lớp đất nền. Các khái niệm nền tảng được làm rõ bao gồm: hiện tượng hóa lỏng đất nền, tỷ số ứng suất tuần hoàn (CSR), tỷ số cường độ kháng cắt tuần hoàn (CRR), áp lực nước lỗ rỗng thặng dư và hệ số an toàn ổn định mái dốc (FS).

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

Nguồn dữ liệu đầu vào được thu thập từ hồ sơ khảo sát địa chất công trình và khí tượng hải văn thực tế tại dự án đê chắn sóng cảng Dung Quất. Cỡ mẫu nghiên cứu bao gồm số liệu từ 12 lỗ khoan khảo sát địa chất với độ sâu khoan từ 20m đến 42m dọc theo tim đê, kết hợp 85 mẫu thí nghiệm cơ lý đất nguyên dạng trong phòng thí nghiệm cơ học đất. Phương pháp chọn mẫu phân tầng có hệ thống được áp dụng nhằm bao quát toàn bộ các tầng trầm tích bùn sét, cát hạt mịn rời rạc và tầng cuội sỏi chịu lực bên dưới đáy biển.

Quy trình mô phỏng số được thực hiện qua chuỗi 4 mô-đun liên kết chặt chẽ trong bộ phần mềm địa kỹ thuật chuyên dụng GeoStudio:

  • Mô-đun Seep/W: Phân tích trường thấm và xác lập đường bão hòa thủy tĩnh qua thân đê ứng với mực nước triều thiết kế.
  • Mô-đun Sigma/W: Mô phỏng trạng thái ứng suất tĩnh ban đầu và trường biến dạng tự nhiên của nền móng trước khi xuất hiện địa chấn.
  • Mô-đun Quake/W: Gán băng ghi gia tốc động đất theo thời gian với gia tốc nền cực đại từ 0,09g đến 0,22g, từ đó tính toán sự lan truyền sóng cắt, mức độ gia tăng áp lực nước lỗ rỗng và khoanh vùng đất bị hóa lỏng.
  • Mô-đun Slope/W: Tính toán hệ số an toàn trượt mái đê theo các phương pháp cân bằng giới hạn lát cắt Bishop và Morgenstern-Price tại từng bước thời gian dao động.

Lý do lựa chọn phương pháp phân tích động học tích hợp này là vì nó khắc phục hoàn toàn nhược điểm của phương pháp giải tích giản đơn, cho phép phản ánh chính xác sự thay đổi phi tuyến của độ cứng vật liệu và áp lực nước lỗ rỗng thặng dư theo thời gian thực của trận động đất kéo dài từ 15 đến 30 giây.

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

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

Quá trình tính toán theo công thức kinh nghiệm và mô phỏng số trên mô hình mặt cắt đê chắn sóng Dung Quất đã mang lại các phát hiện khoa học quan trọng:

  • Đánh giá theo phương pháp Seed (1983) và Idriss - Boulanger (2004) cho thấy tại tầng cát hạt mịn bão hòa ở độ sâu từ -15m đến -25m, hệ số an toàn chống hóa lỏng đều giảm xuống dưới 1,0 (dao động trong khoảng 0,72 đến 0,86), khẳng định tầng đất này bị hóa lỏng hoàn toàn dưới tác động của động đất cấp VII có gia tốc đỉnh 0,09g.
  • Kết quả mô phỏng động lực học bằng Quake/W chỉ ra áp lực nước lỗ rỗng thặng dư tích tụ cực nhanh trong 10 giây đầu của rung chấn, làm suy giảm từ 48% đến 62% cường độ kháng cắt không thoát nước của tầng nền cát yếu.
  • Phân tích ổn định mái dốc bằng Slope/W ghi nhận hệ số an toàn ổn định tổng thể của mái đê giảm từ 1,68 ở trạng thái tĩnh xuống mức nguy hiểm 1,06 khi chịu tác động của phổ gia tốc động đất thiết kế, và rơi xuống mức phá hoại 0,89 nếu không xử lý tầng đất hóa lỏng.
  • Khi áp dụng giải pháp thay thế lớp cát bề mặt bằng lớp đệm đá hộc dày 2,5m đến 3,5m, hệ số an toàn ổn định mái đê trong điều kiện động đất tăng vọt lên 1,34, tương đương mức cải thiện hơn 26% độ ổn định công trình.

Thảo luận kết quả

Nguyên nhân cốt lõi dẫn đến hiện tượng suy giảm ổn định đột ngột của đê chắn sóng là do cấu trúc hạt cát rời bão hòa bị phá vỡ dưới tác động của ứng suất cắt chu kỳ. Khi các hạt cát có xu hướng sắp xếp lại ở trạng thái chặt hơn nhưng nước trong lỗ rỗng không kịp thoát ra ngoài, toàn bộ áp lực truyền sang cột nước, biến đất thành một khối bùn lỏng mất hoàn toàn ma sát trong.

Khi so sánh với phương pháp giả tĩnh truyền thống, sự khác biệt xuất hiện rất rõ rệt. Phương pháp giả tĩnh xem tải trọng động đất như một lực quán tính không đổi tác dụng lên khối trượt, cho ra hệ số an toàn xấp xỉ 1,28. Kết quả này cao hơn từ 15% đến 20% so với kết quả phân tích động học thực tế theo thời gian. Sự chênh lệch nguy hiểm này bắt nguồn từ việc phương pháp giả tĩnh bỏ qua hoàn toàn hiện tượng suy giảm cường độ kháng cắt do áp lực nước lỗ rỗng sinh ra trong suốt quá trình rung lắc.

Để biểu diễn trực quan các kết quả này trong hồ sơ kỹ thuật, dữ liệu biến thiên hệ số an toàn được thể hiện tối ưu qua biểu đồ diễn biến hệ số an toàn theo bước thời gian dao động. Đồng thời, phạm vi hóa lỏng và mức độ tập trung ứng suất cắt thặng dư được minh họa sinh động thông qua các bảng mặt cắt không gian 2D với thang màu đẳng trị áp lực nước lỗ rỗng và mạng lưới cung trượt nguy hiểm nhất đi qua chân mái đê.

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

Dựa trên các kết quả mô phỏng và phân tích chuyên sâu, luận văn đưa ra 4 nhóm giải pháp kỹ thuật cụ thể nhằm đảm bảo an toàn tuyệt đối cho đê chắn sóng:

  • Thay thế lớp đất nền dễ hóa lỏng: Nhà thầu thi công tiến hành nạo vét toàn bộ lớp bùn cát yếu ở độ sâu từ 2,0m đến 4,0m dưới đáy đê, thay thế bằng vật liệu đá đổ hoặc cát thô cấp phối tốt có hệ số đầm chặt tối thiểu đạt K = 0,95, thực hiện dứt điểm trong 6 tháng đầu của giai đoạn xử lý nền móng.
  • Gia cố kết cấu chân khay và lăng cơ tiêu sóng: Ban quản lý dự án chỉ đạo đơn vị tư vấn thiết kế mở rộng kích thước lăng cơ đá đổ chân đê thêm 20% và lắp đặt khối phủ bê tông tiêu sóng trọng lượng từ 12 tấn đến 16 tấn nhằm triệt tiêu lực đẩy thủy động và cố định vị trí chân mái trượt trong thời gian thi công thân đê.
  • Thiết lập hệ thống quan trắc tự động thông minh: Chủ đầu tư triển khai lắp đặt hệ thống cảm biến đo áp lực nước lỗ rỗng điện tử kết hợp máy đo gia tốc 3 trục tại 5 vị trí mặt cắt xung yếu nhất của tuyến đê, thiết lập chế độ cảnh báo tự động 24/7 khi áp lực nước lỗ rỗng vượt quá 70% áp lực bản thân có hiệu của đất nền.
  • Cập nhật quy chuẩn thiết kế kháng chấn công trình thủy: Các cơ quan quản lý chuyên ngành cần ban hành quy định bắt buộc áp dụng phương pháp phân tích động học theo chuỗi mô-đun phần tử hữu hạn cho các dự án đê kè cấp I và cấp II trên phạm vi cả nước trong giai đoạn 2024 - 2026, từng bước thay thế phương pháp giả tĩnh đối với các công trình ven biển có nền cát bão hòa.

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

Công trình nghiên cứu cung cấp khối lượng dữ liệu và phương pháp luận giá trị cao cho 4 nhóm đối tượng trọng tâm:

  • Kỹ sư thiết kế công trình cảng và đường thủy: Vận dụng trọn vẹn quy trình mô phỏng 4 bước trên GeoStudio để tính toán chính xác ổn định mái dốc và tối ưu hóa khối lượng vật liệu đắp đê chắn sóng.
  • Đơn vị tư vấn giám sát và nhà thầu thi công nền móng: Nắm bắt các tiêu chuẩn kiểm soát đầm nén vật liệu đệm đá và giải pháp kỹ thuật ngăn ngừa cát chảy, hóa lỏng nền đê trong điều kiện thi công ngập nước biển.
  • Giảng viên, học viên cao học và nghiên cứu sinh ngành Địa kỹ thuật xây dựng: Sử dụng luận văn như một tài liệu tham khảo học thuật chuẩn mực về việc ứng dụng lý thuyết cơ học đất tới hạn và đánh giá hóa lỏng theo các mô hình bán thực nghiệm quốc tế.
  • Cơ quan thẩm định dự án và ban quản lý cảng biển: Khai thác bộ số liệu đối chứng giữa phương pháp tĩnh và động để xây dựng khung đánh giá rủi ro thiên tai và thẩm định hồ sơ kỹ thuật cho các dự án hạ tầng hàng hải quy mô lớn.

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

  • Tại sao đê chắn sóng xây dựng trên nền cát bão hòa lại có nguy cơ mất ổn định rất cao khi xảy ra động đất? Hiện tượng rung chấn địa chấn diễn ra trong thời gian ngắn ngủi từ 10 đến 30 giây khiến nước trong lỗ rỗng của tầng cát bão hòa không kịp thoát ra ngoài. Áp lực nước lỗ rỗng tăng vọt xấp xỉ áp lực tổng, làm ứng suất hữu hiệu triệt tiêu về 0, biến nền đất thành khối bùn lỏng và gây trượt sạt mái đê.

  • Phương pháp của Idriss và Boulanger (2004) có ưu điểm vượt trội gì so với phương pháp Seed (1983)? Phương pháp Idriss - Boulanger (2004) tích hợp các hàm hiệu chỉnh phi tuyến hiện đại cho hệ số giảm ứng suất theo độ sâu, hệ số quy đổi cấp động đất và hệ số ảnh hưởng của hàm lượng hạt mịn, giúp kết quả đánh giá hóa lỏng đạt độ chính xác cao hơn 15% đến 20% khi ứng dụng tại hiện trường.

  • Vì sao không nên chỉ sử dụng phương pháp giả tĩnh để tính toán ổn định đê chắn sóng chịu động đất? Phương pháp giả tĩnh quy đổi động đất thành một lực quán tính tĩnh không đổi, hoàn toàn bỏ qua quá trình tích tụ áp lực nước lỗ rỗng thặng dư và sự suy giảm độ cứng của đất theo thời gian. Điều này dẫn tới việc ước tính hệ số an toàn cao hơn thực tế từ 15% đến 20%, tạo ra sự an toàn ảo rất nguy hiểm.

  • Bộ phần mềm GeoStudio kết nối các mô-đun như thế nào để giải quyết bài toán động lực học? Quy trình thực hiện tuần tự qua 4 mô-đun: Seep/W xác định đường bão hòa thấm, Sigma/W tính ứng suất tĩnh ban đầu, Quake/W phân tích sóng động đất và xác định vùng hóa lỏng, cuối cùng Slope/W nhập toàn bộ trường ứng suất động để tính toán hệ số an toàn trượt mái.

  • Giải pháp xử lý nền nào mang lại hiệu quả cao nhất cho đê chắn sóng Dung Quất? Biện pháp nạo vét thay thế tầng cát bão hòa mặt bằng lớp đệm đá hộc dày từ 2,5m đến 3,5m kết hợp đầm chặt đạt K = 0,95 là giải pháp tối ưu nhất, giúp tăng hệ số an toàn ổn định động từ mức nguy hiểm 1,06 lên mức an toàn 1,34.

Kết luận

  • Luận văn đã phân tích và đánh giá thành công nguy cơ hóa lỏng nền đê chắn sóng Dung Quất dưới tác động của động đất từ cấp VI đến cấp VIII theo thang MSK-64.
  • Xây dựng thành công quy trình mô phỏng số động lực học 4 bước khép kín bằng bộ công cụ GeoStudio, phản ánh trung thực tương tác ứng suất - biến dạng và áp lực nước lỗ rỗng theo thời gian.
  • Chứng minh rõ nét tính thiếu an toàn của phương pháp giả tĩnh khi ước tính hệ số an toàn cao hơn thực tế khoảng 15% đến 20% do bỏ qua biến thiên áp lực nước lỗ rỗng.
  • Xác định giải pháp kỹ thuật thay thế nền cát yếu bằng đệm đá hộc dày 2,5m đến 3,5m giúp nâng cao hệ số an toàn chống trượt động lực học lên trên 1,34.
  • Đóng góp nguồn tài liệu tham khảo khoa học chất lượng cao, phục vụ công tác thiết kế, thẩm định và thi công các công trình bảo vệ bờ biển trên nền địa chất phức tạp tại Việt Nam.

Kế hoạch tiếp theo đòi hỏi các đơn vị thiết kế và quản lý hạ tầng cảng biển cần nhanh chóng tích hợp phương pháp phân tích động học vào quy trình thẩm tra dự án ngay trong giai đoạn 2024 - 2025. Hãy chủ động ứng dụng công nghệ mô phỏng số tiên tiến và hệ thống quan trắc địa kỹ thuật hiện đại để nâng cao tuổi thọ và độ an toàn tuyệt đối cho mọi công trình hạ tầng ven biển trọng điểm.