Tổng quan nghiên cứu

Thảm họa sụp đổ hố đào sâu 33.5m tại dự án tàu điện ngầm MRT Circle Line trên đường Nicoll Highway (Singapore) xảy ra vào lúc 15:30 ngày 20/04/2004 là một trong những sự cố địa kỹ thuật nghiêm trọng nhất trong lịch sử xây dựng công trình ngầm thế giới. Sự cố này đã làm sụt lún hoàn toàn một đoạn đường cao tốc 6 làn xe, khiến 4 người thiệt mạng và gây thiệt hại kinh tế ước tính hàng trăm triệu USD. Đề tài luận văn thạc sĩ chuyên ngành Kỹ thuật Xây dựng Công trình ngầm tại Trường Đại học Bách khoa – ĐHQG TP.HCM, hoàn thành vào tháng 06/2016, tập trung giải quyết bài toán phân tích chuyên sâu cơ chế mất ổn định của hố đào quy mô lớn này nhằm ngăn ngừa các thảm họa tương tự trong tương lai.

Mục tiêu trọng tâm của nghiên cứu là tái hiện và mô phỏng chính xác toàn bộ quá trình thi công đào đất nhiều giai đoạn bằng phần mềm phần tử hữu hạn Plaxis 2D, đồng thời làm rõ sự tương tác phức tạp giữa kết cấu tường vây bê tông cốt thép dày hơn 1.0m với nền đất sét biển mềm yếu (Singapore Marine Clay). Luận văn đối chiếu các kết quả tính toán mô phỏng số với số liệu quan trắc thực địa trước khi công trình sụp đổ, bao gồm chuyển vị ngang cực đại của tường vây đạt 325mm và độ lún bề mặt đất xung quanh lên tới 150mm. Thông qua việc phân tích lực dọc trong 9 tầng thanh chống và biến dạng uốn của tường chắn, nghiên cứu xác định chính xác các điểm giới hạn phá hoại dẻo trong đất và hư hỏng cục bộ của hệ giằng chống thép.

Kết quả nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc đối với ngành kỹ thuật công trình ngầm đô thị. Luận văn cung cấp khung đánh giá rủi ro định lượng, chứng minh rằng việc áp dụng sai mô hình ứng xử của đất có thể làm sai lệch dự báo áp lực ngang lên tới 25% - 30%, đồng thời đưa ra các chuẩn mực an toàn mới để nâng hệ số ổn định tổng thể cho các dự án ngầm đô thị tại Việt Nam từ mức quy chuẩn 1.2 lên mức an toàn khuyến nghị trên 1.5.

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 cơ học đất cổ điển kết hợp với cơ học đất hiện đại và phương pháp giải tích tiên tiến:

  • Lý thuyết áp lực đất Coulomb (1773) và Mohr – Rankine (1915): Được áp dụng để tính toán áp lực đất chủ động và bị động tác dụng lên tường vây trong điều kiện thoát nước lẫn không thoát nước. Công thức xác định hệ số áp lực đất chủ động $K_a = \tan^2(45^\circ - \phi'/2)$ và bị động $K_p = \tan^2(45^\circ + \phi'/2)$ được sử dụng để thiết lập đường bao ứng suất giới hạn.
  • Lý thuyết ổn định bùng nền đáy hố đào của Terzaghi (1943) và Bjerrum – Eide (1956): Đánh giá nguy cơ trồi đáy hố móng trong tầng sét yếu bão hòa nước với yêu cầu hệ số an toàn chống bùng nền tối thiểu $FS \ge 1.2$.
  • Lý thuyết áp lực đất biểu kiến của Peck (1969): Dùng để ước tính sơ bộ tải trọng phân bố lên hệ thanh chống nhiều tầng (strut load) theo quy tắc diện tích truyền tải một nửa nhịp và phương pháp gối tựa đơn.
  • Mô hình vật liệu phi tuyến Hardening Soil (HS) và Mohr – Coulomb (MC): Mô hình HS mô tả độ cứng phụ thuộc vào ứng suất theo quy luật hyperbol, phân biệt rõ mô-đun biến dạng tải ban đầu $E_{50}^{ref}$, mô-đun nén cố kết $E_{oed}^{ref}$ và mô-đun gia tải/dỡ tải $E_{ur}^{ref}$, giúp nắm bắt chính xác hiện tượng biến cứng dẻo của đất sét biển.

Các khái niệm kỹ thuật then chốt bao gồm: ứng xử không thoát nước dạng A và B (Undrained A, Undrained B), cơ chế phá hoại cắt tổng thể, hiện tượng chân tường bị đá văng (kick-out), và hiệu ứng suy giảm độ cứng dọc trục của hệ giằng chống do biến dạng lệch tâm.

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

Nghiên cứu sử dụng phương pháp mô phỏng số phần tử hữu hạn 2D kết hợp đối chiếu giải tích và phân tích dữ liệu điều tra thực địa:

  • Nguồn dữ liệu: Bộ số liệu địa chất và hồ sơ thi công được thu thập từ báo cáo chính thức của Ủy ban Điều tra Sự cố Singapore (COI Report 2005) cùng các công trình nghiên cứu hiện trường của Whittle & Davies (2006). Cỡ mẫu khảo sát bao gồm chuỗi dữ liệu địa kỹ thuật từ 10 giếng khoan sâu đại diện, xác định chi tiết các chỉ tiêu cơ lý của các lớp đất sét biển trên (Upper Marine Clay) và đất sét biển dưới (Lower Marine Clay) với sức kháng cắt không thoát nước $s_u$ dao động từ 20 kPa đến 80 kPa.
  • Phương pháp phân tích: Sử dụng phần mềm Plaxis 2D với lưới phần tử hữu hạn tam giác 15 nút có độ mịn cao để mô phỏng chính xác trường biến dạng và gradient ứng suất. Phương pháp này được lựa chọn vì có khả năng mô phỏng liên tục trình tự thi công theo thời gian thực, sự tương tác phi tuyến giữa đất và kết cấu qua các phần tử tiếp xúc (interface elements), cũng như việc gia tải/dỡ tải qua từng phase đào.
  • Tiến trình nghiên cứu: Thực hiện từ ngày 11/01/2016 đến ngày 17/06/2016 tại Bộ môn Địa cơ – Nền móng, Trường Đại học Bách khoa TP.HCM. Quá trình mô phỏng bao gồm 9 giai đoạn đào đất tương ứng với 9 tầng thanh chống được lắp đặt và kích tải trước (pre-loading) từ độ sâu tự nhiên xuống đáy hố móng 33.5m.

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ố liệu mô phỏng từ phần mềm Plaxis 2D và so sánh với dữ liệu quan trắc thực tế tại thời điểm trước khi sụp đổ mang lại những phát hiện quan trọng:

  • Chuyển vị ngang của tường vây ($U_{x,\max}$): Mô hình Hardening Soil (HS) cho kết quả chuyển vị ngang cực đại là 283mm, bám sát giá trị quan trắc thực địa tại hiện trường là 325mm (độ sai lệch chỉ 12.9%). Trong khi đó, mô hình Mohr – Coulomb theo phương pháp Undrained B (MC-B) dự báo chuyển vị lên tới 355mm (vượt 9.2% so với thực tế), còn mô hình Mohr – Coulomb theo phương pháp Undrained A (MC-A) chỉ đạt 255mm (đánh giá thấp biến dạng tới 21.5%).
  • Mô-men uốn lớn nhất của tường vây ($M_{\max}$): Mô hình MC-B cho giá trị mô-men uốn lớn nhất đạt 4719 kNm/m tại khu vực tầng sét dưới. Mô hình HS cho giá trị $M_{\max} = 4480\text{ kNm/m}$, cao hơn 17.7% so với giá trị tính toán bằng mô hình MC-A ($M_{\max} = 3807\text{ kNm/m}$). Sự gia tăng đột biến mô-men này vượt quá khả năng chịu uốn đàn hồi định mức của tiết diện tường vây bê tông.
  • Độ lún bề mặt nền đất xung quanh ($U_{y,\max}$): Mô hình HS dự báo độ lún tối đa là 163mm và MC-B là 162mm, tương thích hoàn toàn với số liệu quan trắc ngày 04/04/2004 tại công trường ghi nhận độ lún là 150mm. Mô hình MC-A cho kết quả thấp hơn đáng kể với 123mm.
  • Quá tải nghiêm trọng hệ thanh chống (Strut System): Nội lực tính toán tại hệ thanh chống tầng 8 và tầng 9 (đặc biệt là thanh S335) tăng vọt vượt quá 30% khả năng chịu tải thiết kế khi đào qua lớp Jet Grouting (JGP), trực tiếp kích hoạt quá trình phá hoại mất ổn định dây chuyền.

Thảo luận kết quả

Các kết quả mô phỏng có thể được trực quan hóa qua biểu đồ trắc dọc chuyển vị tường vây theo chiều sâu và biểu đồ bao mô-men uốn qua 9 giai đoạn thi công đào đất:

+-------------------------------------------------------------------------+
|                  SO SÁNH CÁC CHỈ SỐ KỸ THUẬT QUAN TRỌNG                 |
+--------------------------+----------+----------+---------+--------------+
| Chỉ số phân tích         | Mô hình  | Mô hình  | Mô hình | Quan trắc    |
|                          | MC-A     | MC-B     | HS      | thực địa     |
+--------------------------+----------+----------+---------+--------------+
| Chuyển vị ngang cực đại  | 255 mm   | 355 mm   | 283 mm  | 325 mm       |
| Mô-men uốn lớn nhất      | 3807 kNm | 4719 kNm | 4480 kNm| Vượt thiết kế|
| Độ lún đất nền xung quanh| 123 mm   | 162 mm   | 163 mm  | 150 mm       |
| Đánh giá độ chính xác    | Thấp     | Khá      | Rất cao | Chuẩn đối chiếu|
+--------------------------+----------+----------+---------+--------------+

Nguyên nhân cốt lõi dẫn đến sự sụp đổ là sự kết hợp của hai yếu tố: đánh giá thấp áp lực đất trong điều kiện không thoát nước và sai sót nghiêm trọng trong thiết kế chi tiết cấu tạo liên kết thanh giằng. Việc sử dụng mô hình Mohr – Coulomb đơn giản (MC-A) trong giai đoạn thiết kế ban đầu đã bỏ qua hiện tượng biến dạng dẻo tích lũy của đất sét yếu khi dỡ tải, dẫn đến việc thiết kế thiếu khả năng chịu tải của hệ chống tới 20% - 25%.

Bên cạnh đó, việc nhà thầu thi công tự ý thay đổi cấu tạo liên kết gia cường đầu thanh chống từ thép sườn gia cường (stiffeners) dạng hàn đối xứng sang thép chữ C (C-channel) đã làm giảm hơn 50% khả năng chịu tải chống uốn dọc trục khi chịu tác động của độ lệch tâm. Khi đáy hố móng xuất hiện hiện tượng đẩy trồi làm hệ cọc phụ (king post) bị nâng lên, toàn bộ hệ giằng chống bị uốn lệch trục, gây mất ổn định cục bộ tại thanh chống S335 và dẫn đến sự sụp đổ dây chuyền toàn bộ công trình trong vòng chưa đầy 30 phút.

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

Dựa trên kết quả phân tích cơ chế phá hoại của công trình Nicoll Highway, luận văn đưa ra 4 nhóm giải pháp kỹ thuật cụ thể:

  1. Bắt buộc áp dụng mô hình đất phi tuyến nâng cao trong thiết kế hố đào sâu trên 20m: Các đơn vị tư vấn thiết kế phải chuyển đổi từ mô hình Mohr – Coulomb sang mô hình Hardening Soil (HS) hoặc Hardening Soil with Small-Strain Stiffness (HSS) khi phân tích đất dính bão hòa nước. Ứng xử cơ học phải được thiết lập theo phương pháp Undrained B để tính toán chính xác sự gia tăng áp lực nước lỗ rỗng và độ suy giảm cường độ chống cắt. Chủ thể thực hiện: Tư vấn thiết kế địa kỹ thuật. Thời hạn hoàn thành: Giai đoạn lập hồ sơ thiết kế kỹ thuật (1 - 3 tháng trước khi khởi công).
  2. Chuẩn hóa cấu tạo chi tiết liên kết và nẹp tăng cường hệ thanh chống thép: Nghiêm cấm tuyệt đối việc sử dụng các chi tiết liên kết đơn giản như thanh ốp chữ C mà không có sườn tăng cường chịu lực trực tiếp. Mọi liên kết tại nút giao giữa thanh chống (strut) và dầm văng (waler) phải được thiết kế với hệ sườn tăng cường hàn đối xứng, có khả năng chịu tải lệch trục tối thiểu bằng 15% lực dọc thiết kế. Chủ thể thực hiện: Đơn vị thiết kế kết cấu và nhà thầu chế tạo kết cấu thép. Thời hạn hoàn thành: Trước khi tiến hành gia công và lắp đặt tại hiện trường.
  3. Thiết lập hệ thống quan trắc tự động liên tục 24/7 với ngưỡng cảnh báo đa cấp độ: Triển khai các thiết bị đo nghiêng tự động (Inclinometer), cảm biến đo ứng suất thanh chống (Loadcell) và mốc đo lún bề mặt có khả năng truyền dữ liệu thời gian thực. Thiết lập ngưỡng báo động đỏ khi tốc độ dịch chuyển ngang của tường vây vượt quá 5mm/ngày hoặc độ lệch trục của thanh chống đạt quá 20mm. Chủ thể thực hiện: Nhà thầu quan trắc độc lập và Ban Quản lý dự án. Thời hạn thực hiện: Toàn bộ vòng đời thi công hố đào (kéo dài từ 12 đến 24 tháng).
  4. Kiểm soát nghiêm ngặt hiện tượng trồi cọc king post và thi công đáy móng: Quy định chiều sâu mỗi đợt đào đất không vượt quá 2.0m - 2.5m và phải hoàn thành việc lắp đặt, kích tải trước thanh chống trong vòng 24 giờ sau khi đào. Đối với các hố đào sâu trên 30m trong đất sét yếu, bắt buộc phải thi công lớp sàn gia cố đáy bằng cọc xi măng đất (Jet Grouting) với chiều dày tối thiểu 2.0m nhằm triệt tiêu nguy cơ bùng nền và đá chân tường. Chủ thể thực hiện: Chỉ huy trưởng công trường và Tư vấn giám sát trưởng. Thời hạn thực hiện: Hằng ngày trong suốt quá trình đào hố móng.

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

Luận văn là tài liệu chuyên khảo giá trị dành cho 4 nhóm đối tượng chính trong ngành xây dựng công trình ngầm:

  1. Kỹ sư thiết kế địa kỹ thuật và kết cấu ngầm: Nắm vững quy trình thiết lập thông số đầu vào cho phần mềm Plaxis 2D, hiểu rõ sự khác biệt bản chất giữa các mô hình đất (Mohr – Coulomb, Hardening Soil) và các phương pháp phân tích không thoát nước (Undrained A, B, C) để tránh các sai số chí mạng trong tính toán nội lực tường vây và hệ chống.
  2. Nhà thầu thi công và kỹ sư hiện trường: Hiểu rõ tầm quan trọng sống còn của việc tuân thủ tuyệt đối bản vẽ thiết kế biện pháp thi công, nhận diện các hiểm họa tiềm ẩn từ việc tự ý thay đổi chi tiết nẹp liên kết thép hoặc đào đất vượt quá chiều sâu cho phép mà chưa kịp lắp thanh chống.
  3. Chuyên viên tư vấn giám sát và quản lý dự án: Trang bị cơ sở khoa học để kiểm tra, đánh giá tính khả thi của biện pháp thi công hố đào sâu đô thị, đồng thời xây dựng quy trình giám sát số liệu quan trắc thực địa, thiết lập các kịch bản ứng phó sự cố khẩn cấp khi biến dạng vượt ngưỡng cho phép.
  4. Học viên cao học, nghiên cứu sinh và giảng viên chuyên ngành Địa kỹ thuật – Công trình ngầm: Sử dụng làm tài liệu tham khảo học thuật chuyên sâu về bài toán phân tích phá hoại công trình ngầm quy mô lớn, làm cơ sở để phát triển các đề tài nghiên cứu mở rộng về tương tác động lực học đất – kết cấu và mô hình hóa 3D không gian hố đào.

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

1. Nguyên nhân kỹ thuật chính gây ra sự cố sập hố đào Nicoll Highway 33.5m là gì?

Sự cố bắt nguồn từ việc mô hình tính toán thiết kế ban đầu đánh giá thấp áp lực đất không thoát nước, kết hợp với lỗi thay đổi chi tiết liên kết thép tăng cường từ dạng sườn hàn đối xứng sang thép chữ C. Điều này làm thanh chống S335 tại tầng 9 bị mất ổn định uốn dọc trục khi tường vây dịch chuyển ngang đạt 325mm, gây sập đổ dây chuyền.

2. Vì sao mô hình Hardening Soil (HS) lại chính xác hơn mô hình Mohr – Coulomb (MC) trong bài toán này?

Mô hình Hardening Soil tính đến sự phụ thuộc của độ cứng đất vào mức độ ứng suất và phân biệt giữa độ cứng gia tải ban đầu ($E_{50}$) với độ cứng dỡ tải/tải lại ($E_{ur}$). Trong mô phỏng hố đào sâu 33.5m, mô hình HS cho chuyển vị tường 283mm và lún 163mm, sát nhất với thực tế quan trắc (325mm và 150mm), trong khi mô hình MC-A chỉ cho chuyển vị 255mm.

3. Sự khác nhau cốt lõi giữa phương pháp Undrained A và Undrained B trong Plaxis là gì?

Phương pháp Undrained A sử dụng các thông số sức kháng cắt hiệu quả ($c', \phi'$) và mô-đun đàn hồi hiệu quả, tự động tính toán áp lực nước lỗ rỗng thặng dư. Phương pháp Undrained B sử dụng trực tiếp sức kháng cắt không thoát nước ($s_u, \phi=0$), phản ánh sát thực tế trạng thái tới hạn trong đất sét biển bão hòa, cho kết quả chuyển vị lớn hơn khoảng 28% so với Undrained A.

4. Chi tiết nẹp bản cánh C-channel đã ảnh hưởng như thế nào đến khả năng chịu lực của thanh chống?

Các thử nghiệm kiểm chứng chỉ ra rằng việc sử dụng bản ốp chữ C làm tâm chịu lực của liên kết bị lệch khỏi trục tim của thanh chống thép hình H. Khi hố đào xuất hiện biến dạng ngang và đáy móng trồi nhẹ, độ lệch tâm này tạo ra mô-men xoắn phụ làm giảm hơn 50% khả năng chịu nén tới hạn của thanh chống so với thiết kế chuẩn ban đầu.

5. Biện pháp nào hiệu quả nhất để ngăn ngừa hiện tượng bùng nền và đá chân tường trong đất sét yếu?

Giải pháp hiệu quả nhất là kết hợp kéo dài chiều sâu ngàm của tường vây xuống tầng đất chịu lực tốt bên dưới và thi công lớp cọc gia cố đáy bằng công nghệ Jet Grouting (JGP) dày từ 1.5m đến 3.0m. Lớp sàn này đóng vai trò như một thanh chống ngầm dưới đáy, triệt tiêu hiện tượng chân tường bị đẩy ngang và nâng hệ số an toàn chống trồi đáy $FS$ lên trên 1.5.

Kết luận

Nghiên cứu cơ chế mất ổn định hố đào sâu 33.5m công trình Nicoll Highway – Singapore đã mang lại những đóng góp học thuật và thực tiễn cốt lõi:

  • Làm sáng tỏ toàn diện cơ chế phá hoại liên hoàn giữa nền đất sét biển mềm và kết cấu chống đỡ nhiều tầng dưới tác động của ứng suất không thoát nước.
  • Khẳng định tính vượt trội của mô hình Hardening Soil trong phần mềm Plaxis 2D khi tái hiện chính xác trường biến dạng với sai số dưới 13% so với số liệu quan trắc thực tế.
  • Chỉ ra bài học kinh nghiệm sâu sắc về việc kiểm soát độ cứng liên kết thanh chống, cảnh báo rủi ro phá hoại do biến dạng lệch tâm và hiện tượng trồi cọc king post.
  • Đề xuất phương pháp luận tính toán kiểm chứng kép kết hợp giữa mô phỏng số nâng cao và phân tích giải tích cổ điển Rankine – Terzaghi.
  • Xây dựng bộ chỉ dẫn kỹ thuật an toàn ứng dụng trực tiếp cho các dự án hố đào ngầm và tuyến đường sắt đô thị (Metro) đang triển khai tại Việt Nam.

Kế hoạch phát triển tiếp theo của nghiên cứu là mở rộng mô phỏng không gian 3D nhằm đánh giá hiệu ứng góc hố móng và ảnh hưởng của rung động động lực học trong thời gian từ 6 đến 12 tháng tới. Các kỹ sư, nhà thầu và cơ quan quản lý cần tích cực tham khảo và áp dụng nghiêm ngặt các bài học đắt giá từ thảm họa Nicoll Highway để đảm bảo an toàn tuyệt đối cho mọi công trình ngầm đô thị.