Tổng quan nghiên cứu

Quá trình đô thị hóa nhanh chóng tại Thành phố Hồ Chí Minh thúc đẩy nhu cầu khai thác không gian ngầm với các công trình quy mô lớn như Vietcombank Tower (4 tầng hầm), Bitexco Financial Tower (3 tầng hầm), Saigon One Tower (5 tầng hầm) hay hố đào dự án Nhiêu Lộc – Thị Nghè sâu đến 40 m. Trong điều kiện mật độ xây dựng cao, việc thi công hố đào sâu luôn tiềm ẩn rủi ro chuyển vị tường vây và lún nứt công trình lân cận, trong đó sai số từ mô hình hóa địa kỹ thuật là nguyên nhân hàng đầu. Vấn đề cốt lõi xuất phát từ việc lấy mẫu đất cát nguyên dạng tại hiện trường vô cùng khó khăn, mẫu cát dễ bị xáo trộn làm sai lệch đáng kể thông số độ cứng và góc ma sát trong phòng thí nghiệm.

Nghiên cứu này được thực hiện nhằm giải quyết triệt để bài toán xác định module biến dạng thực tế của đất cát bằng phương pháp phân tích ngược dựa trên số liệu quan trắc hiện trường qua toàn bộ các giai đoạn thi công. Phạm vi nghiên cứu tập trung vào môi trường địa chất hố đào sâu và vùng ảnh hưởng mở rộng gấp 2 lần bề rộng hố đào tại Thành phố Hồ Chí Minh, điển hình qua hai công trình trọng điểm là Vietcombank Tower và LIM Tower trong giai đoạn từ tháng 1 đến tháng 6 năm 2015.

Ý nghĩa khoa học và thực tiễn của công trình thể hiện ở việc làm giàu cơ sở dữ liệu địa kỹ thuật địa phương, thiết lập các mối tương quan thực nghiệm chính xác giữa module biến dạng và chiều sâu. Giải pháp này giúp nâng cao độ chính xác khi mô phỏng số, giảm thiểu sai số dự báo chuyển vị tường vây xuống dưới 10%, kiểm soát độ lún nền đất xung quanh và bảo đảm an toàn tuyệt đối cho các công trình hạ tầng đô thị lân cậ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 đồng thời hai mô hình ứng xử cơ học đất kinh điển và nâng cao để mô phỏng tương tác đất - kết cấu trong bài toán dỡ tải hố đào sâu:

Mô hình đàn dẻo lý tưởng Mohr-Coulomb: Sử dụng 5 thông số cơ bản gồm module biến dạng $E$, hệ số Poisson $\nu$, lực dính hiệu quả $c'$, góc ma sát trong $\phi'$ và góc giãn nở $\psi$. Mô hình này đơn giản, phổ biến nhưng xem độ cứng là hằng số tuyến tính trong từng phân lớp đất.

Mô hình Hardening Soil: Mô hình phi tuyến bậc cao mô phỏng sát thực tế ứng xử của đất chịu tải và dỡ tải thông qua đường cong hyperbol Duncan-Chang. Mô hình sử dụng các thông số độ cứng phụ thuộc ứng suất gồm module cát tuyến ba trục $E_{50}^{ref}$, module dỡ tải và tái gia tải $E_{ur}^{ref} = 3E_{oed}^{ref}$, module nén cố kết không nở hông $E_{oed}^{ref}$, số mũ phụ thuộc ứng suất $m = 0.5$, hệ số Poisson dỡ tải $\nu_{ur} = 0.2$ và hệ số phá hoại $R_f = 0.9$.

Khung lý thuyết tính toán mềm: Tích hợp lý thuyết Mạng thần kinh nhân tạo lan truyền tiến với thuật toán lan truyền ngược sai số kết hợp Thuật toán di truyền để giải quyết bài toán tối ưu hóa đa biến trong phân tích ngược địa kỹ thuật.

Các khái niệm trọng tâm bao gồm: Module biến dạng của đất cát, chuyển vị ngang tường vây (wall deflection), độ lún đất nền xung quanh hố đào và hàm sai số căn quân phương $R_h$.

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

Nghiên cứu kết hợp phương pháp phần tử hữu hạn với kỹ thuật trí tuệ nhân tạo thông qua chu trình tính toán ngược đa bước:

Nguồn dữ liệu và cỡ mẫu: Sử dụng chuỗi số liệu quan trắc độ sâu chuyển vị ngang tường vây thực tế thu thập từ hệ thống ống đo nghiêng Inclinometer qua 6 giai đoạn thi công đào đất tại công trình Vietcombank Tower, tương ứng với 140 điểm đo chuyển vị chuẩn hóa theo phương đứng.

Phương pháp mô phỏng số: Mô phỏng bài toán biến dạng phẳng bằng phần mềm Plaxis 2D, chia nhỏ khối đất thành các phân lớp đơn để xác định bộ 10 thông số module biến dạng chưa biết đại diện cho các tầng địa chất.

Phương pháp phân tích ngược AI: Thiết lập mạng nơ-ron nhân tạo gồm 1 lớp vào (10 nút), 1 lớp ẩn (8 tế bào thần kinh) và 1 lớp ra (140 tế bào thần kinh) sử dụng hàm truyền Sigmoid để học xấp xỉ kết quả từ phần tử hữu hạn. Sau đó, Thuật toán di truyền thực hiện tìm kiếm toàn cục bộ thông số module biến dạng tối ưu nhằm cực tiểu hóa hàm mục tiêu sai số $R^*$.

Lý do lựa chọn: Sự kết hợp giữa mạng nơ-ron và thuật toán di truyền loại bỏ nhược điểm bẫy cực trị địa phương của các phương pháp tối ưu hóa truyền thống, giúp tái lập chính xác ứng xử phi tuyến của đất cát mà không cần thực hiện thí nghiệm đóng băng mẫu tốn kém.

Timeline nghiên cứu: Toàn bộ quá trình thu thập số liệu, huấn luyện mô hình số và kiểm chứng thực nghiệm hoàn thành trong khoảng thời gian từ tháng 1/2015 đến tháng 6/2015.

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

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

Quá trình tối ưu hóa và phân tích ngược đã mang lại các phát hiện cốt lõi:

Thứ nhất, xác định thành công bộ thông số độ cứng tối ưu cho đất nền Thành phố Hồ Chí Minh. Đối với mô hình Hardening Soil, bộ thông số tối ưu số 35 cho thấy giá trị tham chiếu $E_{50}^{ref}$ của các lớp đất cát dao động ổn định trong khoảng từ 20.000 kN/m² đến 45.000 kN/m². Đối với mô hình Mohr-Coulomb, bộ thông số tối ưu số 34 phản ánh module biến dạng $E$ tăng mạnh theo chiều sâu.

Thứ hai, thiết lập quy luật biến thiên module biến dạng theo độ sâu. Nghiên cứu chứng minh rõ rệt: Trong mô hình Mohr-Coulomb, module biến dạng $E$ của đất cát tăng tỷ lệ thuận với căn bậc hai của chiều sâu $z$ dưới tác động của ứng suất hữu hiệu bản thân gia tăng, tuân theo công thức thực nghiệm $E = (10.200 - 15.400)z^{0.5}$ (kN/m²). Trong khi đó, với mô hình Hardening Soil, module tham chiếu $E_{50}^{ref}$ duy trì gần như dạng hằng số do tính chất phụ thuộc ứng suất đã được tự động chuẩn hóa qua hàm toán học nội tại.

Thứ ba, giảm thiểu vượt bậc sai số mô phỏng. Kết quả tính toán chuyển vị tường vây từ bộ thông số tối ưu giúp giảm sai số căn quân phương $R_h$ hơn 45% so với việc gán thông số giả định ban đầu từ chỉ số SPT.

Thứ tư, kiểm chứng thành công trên công trình độc lập LIM Tower. Khi áp dụng các công thức tương quan vào mô hình hố đào LIM Tower, độ lệch chuyển vị ngang lớn nhất giữa mô phỏng và quan trắc thực tế chỉ ở mức 8% đến 12%, trong khi chuyển vị lún mặt đất lân cận đạt độ khớp trên 90%.

Thảo luận kết quả

Sự vượt trội của mô hình Hardening Soil so với mô hình Mohr-Coulomb bắt nguồn từ khả năng tách biệt module biến dạng khi gia tải sơ cấp ($E_{50}$) và dỡ tải tái gia tải ($Eur$). Trong thi công hố đào sâu, việc đào đất thực chất là quá trình dỡ tải mạnh mẽ. Mô hình Mohr-Coulomb do sử dụng một module biến dạng không đổi nên thường đánh giá quá cao độ trồi đáy móng và dự báo chuyển vị ngang thiếu chính xác nếu không hiệu chỉnh module tăng theo chiều sâu.

So sánh với các công bố quốc tế, kết quả này hoàn toàn tương đồng với nghiên cứu của Khoiri và Ou tại Cao Hùng (Đài Loan) cũng như 17 phân tích phần tử hữu hạn tại Berlin (Đức), nhưng đã được hiệu chỉnh chính xác theo đặc điểm địa chất trầm tích hạt mịn đặc thù của lưu vực sông Sài Gòn.

Dữ liệu nghiên cứu được trực quan hóa sinh động thông qua biểu đồ biến thiên chuyển vị ngang dạng vòm của tường chắn theo độ sâu hố đào, biểu đồ đường cong võng lún mặt đất ngoài hố đào và đồ thị hội tụ sai số $R^*$ và $R_h$ qua các thế hệ lai ghép của thuật toán di truyền. Các bảng so sánh thông số đối chiếu chi tiết giữa giá trị thí nghiệm SPT, giá trị nén ba trục và kết quả phân tích ngược khẳng định độ tin cậy tuyệt đối của phương pháp.

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

Nhằm nâng cao chất lượng thiết kế và an toàn thi công hố đào sâu trong đô thị, luận văn đưa ra 4 nhóm khuyến nghị mang tính ứng dụng thực tiễn cao:

  1. Chuẩn hóa quy trình áp dụng mô hình Hardening Soil trong phân tích kết cấu ngầm: Các đơn vị tư vấn thiết kế cần chuyển đổi từ mô hình Mohr-Coulomb sang mô hình Hardening Soil khi mô phỏng hố đào sâu trên 10 m bằng Plaxis 2D/3D. Mục tiêu nhằm kiểm soát độ lệch chuyển vị ngang tường chắn dưới 10 mm, hoàn thành cập nhật quy trình trong vòng 6 tháng kể từ khi phê duyệt dự án.

  2. Ứng dụng công thức tương quan module biến dạng theo chiều sâu cho địa chất cát Thành phố Hồ Chí Minh: Các kỹ sư địa kỹ thuật nên áp dụng trực tiếp công thức thực nghiệm $E = (10.200 - 15.400)z^{0.5}$ (kN/m²) khi sử dụng mô hình Mohr-Coulomb để tính toán sơ bộ. Giải pháp này giúp cắt giảm 30% thời gian hiệu chỉnh mô hình mà vẫn đảm bảo độ chính xác trên 85%.

  3. Tích hợp công nghệ quan trắc tự động thời gian thực với phân tích ngược AI: Ban quản lý dự án và nhà thầu thi công cần lắp đặt hệ thống đo nghiêng tự động Inclinometer với tần suất ghi nhận dữ liệu 2 lần/ngày trong suốt 100% thời gian thi công hố móng. Dữ liệu này được kết nối trực tiếp với thuật toán mạng nơ-ron để liên tục cập nhật thông số đất nền, cảnh báo sớm biến dạng nguy hiểm trước 48 giờ.

  4. Xây dựng cơ sở dữ liệu số hóa địa kỹ thuật đô thị đa tầng: Sở Xây dựng phối hợp cùng các trường đại học chuyên ngành triển khai cổng dữ liệu số hóa thông số cơ lý đất nền từ hơn 50 dự án hố đào sâu trọng điểm. Kế hoạch triển khai giai đoạn 2025-2027 nhằm mục tiêu giảm thiểu 20% chi phí khảo sát bổ sung cho các công trình hạ tầng ngầm tương lai.

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

Nội dung và kết quả nghiên cứu mang lại giá trị thực tiễn sâu sắc cho 4 nhóm đối tượng chính:

Kỹ sư thiết kế địa kỹ thuật và công trình ngầm: Cung cấp bộ thông số đầu vào chuẩn xác ($E_{50}^{ref}$, $E_{ur}^{ref}$, $E_{oed}^{ref}$) cho phần mềm Plaxis, giúp tối ưu hóa chiều dày tường vây barrette từ 0,8 m đến 1,2 m và giảm từ 10% đến 15% khối lượng giằng chống shoring.

Nhà thầu thi công nền móng và hố đào sâu: Giúp các chỉ huy trưởng công trình nắm bắt chính xác quy luật chuyển vị của vách đất cát qua từng giai đoạn hạ mực nước ngầm và đào đất, từ đó lập biện pháp thi công Top-down hoặc Semi-Top-down an toàn tuyệt đối cho các hố đào sâu từ 3 đến 5 tầng hầm.

Chuyên gia tư vấn giám sát và đơn vị quan trắc: Cung cấp phương pháp luận khoa học để thẩm tra báo cáo quan trắc chuyển vị thành hố đào, đánh giá tính ổn định tổng thể của hố móng và đưa ra giải pháp ứng phó kịp thời khi độ lún công trình lân cận tiệm cận ngưỡng cảnh báo 15 mm.

Giảng viên, nghiên cứu sinh và học viên cao học: Tài liệu tham khảo mẫu mực về việc ứng dụng trí tuệ nhân tạo (Mạng nơ-ron kết hợp Thuật toán di truyền) giải bài toán phân tích ngược trong cơ học đất hiện đại, mở ra hướng nghiên cứu liên ngành giữa khoa học dữ liệu và kỹ thuật xây dựng.

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

Tại sao việc xác định module biến dạng của đất cát từ thí nghiệm trong phòng thường thiếu chính xác? Đất cát có cấu trúc hạt rời rạc, rất khó lấy mẫu nguyên dạng ngoài hiện trường nếu không dùng biện pháp đóng băng tốn kém. Quá trình lấy và vận chuyển mẫu làm xáo trộn độ chặt tự nhiên, khiến giá trị module đàn hồi trong phòng thí nghiệm thường sai lệch từ 30% đến 50% so với độ cứng thực tế.

Mô hình Hardening Soil có ưu điểm gì vượt trội so với mô hình Mohr-Coulomb khi tính toán hố đào? Mô hình Hardening Soil mô phỏng chính xác tính phi tuyến của đất và phân biệt rõ độ cứng khi gia tải sơ cấp với dỡ tải tái gia tải ($E_{ur} \approx 3E_{50}$). Điều này giúp phản ánh đúng hiện tượng giải phóng ứng suất khi đào sâu, triệt tiêu sai số trồi đáy hố móng vốn bị thổi phồng ở mô hình Mohr-Coulomb.

Mạng thần kinh nhân tạo kết hợp thuật toán di truyền đóng vai trò gì trong phân tích ngược? Mạng nơ-ron đóng vai trò là hàm xấp xỉ hóa kết quả mô phỏng phần tử hữu hạn cực nhanh, trong khi thuật toán di truyền thực hiện dò tìm toàn cục để tìm ra bộ 10 thông số module biến dạng tối ưu có sai số $R_h$ nhỏ nhất so với 140 điểm đo chuyển vị thực tế.

Công thức tương quan module biến dạng của đất cát theo chiều sâu có thể áp dụng cho mọi công trình tại TP.HCM không? Công thức $E = (10.200 - 15.400)z^{0.5}$ áp dụng rất tốt cho các tầng cát bồi tích hạt mịn đến hạt trung tại khu vực trung tâm Quận 1, Quận 2 và Quận 7. Tuy nhiên, với các khu vực có địa tầng bùn sét dày biến tính, kỹ sư cần kết hợp thêm chỉ số SPT hiện trường để tinh chỉnh.

Sai số dự báo chuyển vị tường vây được kiểm soát ở mức nào khi áp dụng phương pháp này? Phương pháp phân tích ngược AI giúp kiểm soát sai số chuyển vị ngang lớn nhất của tường vây trong phạm vi từ 8% đến 12% so với thực tế quan trắc, vượt trội hoàn toàn so với mức sai số 35% đến 60% khi sử dụng các bảng tra thông số kinh nghiệm truyền thống.

Kết luận

Nghiên cứu mang lại những đóng góp khoa học và thực tiễn mang tính đột phá cho chuyên ngành Địa kỹ thuật xây dựng:

  • Xác định thành công bộ thông số module biến dạng thực tế của đất cát tại hiện trường Thành phố Hồ Chí Minh thông qua phân tích ngược từ dữ liệu hố đào sâu Vietcombank Tower.
  • Thiết lập quy luật toán học thực nghiệm chuẩn xác giữa module biến dạng $E$ theo chiều sâu $z$ cho mô hình Mohr-Coulomb và khoảng tham chiếu tối ưu $E_{50}^{ref}$ cho mô hình Hardening Soil.
  • Xây dựng quy trình tích hợp thành công giữa Mạng thần kinh nhân tạo (10-8-140 nút), Thuật toán di truyền và Phương pháp phần tử hữu hạn Plaxis 2D trong việc giải bài toán địa kỹ thuật phức tạp.
  • Kiểm chứng độc lập xuất sắc tại công trình LIM Tower với độ chính xác dự báo chuyển vị đạt trên 88%, giảm thiểu tối đa nguy cơ mất ổn định hố đào và lún nứt công trình lân cận.
  • Đặt nền móng vững chắc cho việc ứng dụng công nghệ số và trí tuệ nhân tạo vào thiết kế không gian ngầm đô thị an toàn, bền vững.

Trong giai đoạn tiếp theo từ năm 2025 đến 2027, các kỹ sư và nhà nghiên cứu nên tiếp tục mở rộng áp dụng bộ thông số này cho các dự án nhà ga ngầm Metro Bến Thành – Suối Tiên và các hố đào sâu trên 5 tầng hầm. Hãy chủ động chuẩn hóa mô hình Hardening Soil và số hóa dữ liệu quan trắc ngay hôm nay để tối ưu hóa chi phí và đảm bảo an toàn tuyệt đối cho mọi công trình ngầm đô thị.