Tổng quan về luận án

Nghiên cứu về phản ứng địa chấn của nền đất cục bộ đóng vai trò then chốt trong kỹ thuật địa chấn công trình và quản lý rủi ro thiên tai đô thị. Luận án tiến sĩ chuyên ngành Vật lí Địa cầu (Mã số: 9440130.06) của nghiên cứu sinh Giang Kiên Trung với đề tài "Đánh giá ảnh hưởng của điều kiện địa chất địa phương tới rung động nền do động đất gây ra trên tuyến Tây Đông thành phố Hà Nội" (Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội, 2022, người hướng dẫn: TS. Phạm Đình Nguyên và TS. Nguyễn Đức Vinh) đại diện cho công trình nghiên cứu tiên phong, định lượng hóa có hệ thống các biến đổi của dao động địa chấn dưới tác động của cấu trúc trầm tích Đệ tứ tại thủ đô Hà Nội.

Vùng đô thị Hà Nội tọa lạc trên đới đứt gãy hoạt động Sông Hồng và Sông Chay, có tiềm năng phát sinh động đất với độ lớn $M = 6{,}1 \div 6{,}5$, tương ứng cấp chấn động cực đại cấp VIII theo thang MSK-64 trên nền đá gốc loại A (Nguyễn Đình Xuyên, 2004; Nguyễn Ngọc Thủy và Nguyễn Sinh Minh, 2004). Tuy nhiên, hầu hết diện tích thành phố lại được bao phủ bởi các lớp trầm tích sông hồ, đầm lầy biến tiến - biến lùi Đệ tứ ($amQ$) dày từ vài chục đến hơn một trăm mét với cơ tính biến đổi phức tạp. Khoảng trống nghiên cứu (research gap) cốt lõi tồn tại trước công trình này là sự thiếu hụt các mô hình truyền sóng động lực học 1D định lượng, phân giải cao chạy xuyên suốt không gian địa chất từ vùng đồi gò chuyển tiếp phía Tây sang vùng trầm tích bồi tích dày phía Đông. Phần lớn các công trình trước đây (Nguyễn Hồng Phương, 2002, 2006; Trịnh Việt Bắc et al., 2011; Lê Tử Sơn et al., 2007) chỉ dừng lại ở vi phân vùng theo phương pháp độ cứng địa chấn Medvedev thực nghiệm hoặc khảo sát cục bộ từng điểm đơn lẻ tại khu vực nội thành cũ và phía Tây.

Luận án giải quyết ba câu hỏi nghiên cứu (Research Questions - RQ) và kiểm chứng ba giả thuyết khoa học (Hypotheses - H):

  • RQ1: Cấu trúc địa tầng và phân bố vận tốc sóng cắt ($V_s$) biến đổi như thế nào dọc theo tuyến trắc cảnh Tây - Đông cắt ngang qua các cấu trúc kiến tạo và thủy văn chính của Hà Nội?
  • RQ2: Sự biến thiên về bề dày và độ xốp của các tầng trầm tích Đệ tứ ảnh hưởng định lượng ra sao tới hệ số khuếch đại phổ gia tốc ($RSRA$), gia tốc nền cực đại ($PGA$), và chu kỳ trội ($T_0$) khi chịu tác động của kịch bản động đất nguồn gần và nguồn xa?
  • RQ3: Cơ chế vật lý nào chi phối hiện tượng tiêu tán hoặc khuếch đại bất thường tại các vị trí chứa tầng sét bùn dày trên 30 m?
  • H1: Tồn tại quy luật tăng tuyến tính của gia tốc nền cực đại ($PGA$) và chu kỳ dao động trội ($T_0$) từ Tây sang Đông tỷ lệ thuận với sự gia tăng bề dày tầng bồi tích yếu.
  • H2: Tầng trầm tích mềm yếu gần bề mặt có vận tốc $V_s < 80\text{ m/s}$ với chiều dày lớn hơn 30 m sẽ gây ra hiệu ứng suy giảm phi tuyến tính (non-linear deamplification) đối với biên độ gia tốc rung động mạnh trong khi vẫn dịch chuyển chu kỳ trội về miền tần số thấp.
  • H3: Tích hợp mô hình số 1D phi tuyến với hệ thống dữ liệu lỗ khoan chuẩn hóa cho phép dự báo phân vùng nguy hiểm địa chấn tiệm cận độ chính xác của các khảo sát địa chấn công trình quốc tế.

Khung lý thuyết của nghiên cứu được xây dựng trên cơ sở lý thuyết truyền sóng đàn hồi và dẻo đàn hồi qua môi trường nhiều lớp nằm ngang (Thomson, 1950; Haskell, 1953; Schnabel et al., 1972), kết hợp mô hình suy giảm chuyển động đất thế hệ mới NGA (Campbell & Bozorgnia, 2008) và kỹ thuật phân tích phản ứng nền phi tuyến trong phần mềm DEEPSOIL v7.0. Luận án khảo sát tập dữ liệu gồm hơn 324 lỗ khoan địa chất công trình và thủy văn sâu từ 20 m đến 130 m, lựa chọn 53 lỗ khoan trọng điểm để dựng mặt cắt trắc tuyến Tây - Đông dài hơn 30 km, đóng góp bức tranh toàn cảnh định lượng phục vụ trực tiếp thiết kế kháng chấn theo tiêu chuẩn TCVN 9386:2012 và quy hoạch bền vững đô thị mở rộng.

Literature Review và Positioning

Lịch sử nghiên cứu hiệu ứng nền địa phương (local site effects) ghi nhận sự chuyển dịch từ các quan sát kinh nghiệm mô tả thiệt hại sang các phương pháp vi địa chấn và mô phỏng số động lực học. Tác động khốc liệt của điều kiện địa chất bề mặt đã được chứng minh qua thảm họa động đất Michoacan, Mexico năm 1985 ($M = 8{,}1$), nơi lớp trầm tích sét hồ ao xốp mềm tại trung tâm Mexico City cách tâm chấn hơn 400 km đã khuếch đại sóng địa chấn lên gấp nhiều lần ($PGA = 170\text{ cm/s}^2$ so với $35\text{ cm/s}^2$ tại nền đá ngoài rìa), phá hủy 317 tòa nhà cao tầng và làm hơn 10.000 người thương vong (Anderson et al., 1986; Semblat, 2009). Tương tự, trong trận động đất Chi-Chi tại Đài Loan năm 1999 ($M_w = 7{,}6$), bồn trũng trầm tích Đài Bắc cách xa 160 km vẫn chịu chấn động cấp VIII (MMI) với dao động chu kỳ dài bất thường do hiện tượng cộng hưởng sóng và bẫy sóng mặt (Shin et al., 2000; Wen, 2012).

Trong y văn học thuật quốc tế, tồn tại hai luồng quan điểm đối lập lớn về phương pháp luận đánh giá phản ứng nền:

  1. Trường phái thực nghiệm vi địa chấn (Nakamura, 1989; Field & Jacob, 1995; Apostolidis et al., 2004): Ủng hộ phương pháp tỷ số phổ Fourier giữa thành phần nằm ngang và thẳng đứng ($HVSR$ hay kỹ thuật Nakamura) từ sóng vi địa chấn do nhiễu môi trường, cho rằng phương pháp này đủ độ tin cậy để xác định tần số cộng hưởng cơ bản $f_0$ và hệ số khuếch đại nền với chi phí thấp. Như Nakamura (2019) đã bổ sung: "Phương pháp Nakamura không chỉ sử dụng tỷ số phổ H/V ở mặt đất để ước tính tần số cộng hưởng và hệ số khuếch đại của lớp bề mặt mà cũng sử dụng tỷ số phổ H/V ở đỉnh của tòa nhà để ước tính tần số cộng hưởng và hệ số khuếch đại của tòa nhà và mặt đất."
  2. Trường phái mô phỏng truyền sóng cơ học đất (Idriss & Seed, 1968; Dobry et al., 1976; Bard & Riepl-Thomas, 2000; Hashash et al., 2015): Lập luận rằng phương pháp vi địa chấn tuyến tính không phản ánh được tính phi tuyến của đất cát dưới kích động động đất mạnh, do mô-đun cắt $G$ suy giảm và hệ số cản $\beta$ tăng mạnh khi biến dạng trượt vượt quá $10^{-4}$. Họ khẳng định việc giải phương trình truyền sóng 1D/2D/3D (như mã tính SHAKE91 hoặc DEEPSOIL) kết hợp số liệu địa kỹ thuật lỗ khoan ($N_{SPT}$, $V_s$, $\gamma$) là bắt buộc đối với công trình xây dựng chịu chấn động mạnh.

Tại Việt Nam, các nghiên cứu của Nguyễn Đình Xuyên và nnk. (1990, 1996, 2004) tại Viện Vật lí Địa cầu đã đặt nền móng vi phân vùng địa chấn theo phương pháp độ cứng Medvedev (1962). Tiếp đó, Nguyễn Hồng Phương (2002, 2006) áp dụng tiêu chuẩn NEHRP (1997) và công cụ HAZUS; Lê Tử Sơn et al. (2003, 2007) khảo sát $HVSR$ tại Điện Biên; Nguyễn Tiến Hùng và Kuo-Liang Wen (2011) đo chu kỳ trội vi địa chấn tại Hà Nội ghi nhận xu hướng $T_0$ tăng từ Tây sang Đông và vượt $0{,}6\text{ s}$ dọc Sông Hồng; Thái Anh Tuấn và Nguyễn Đức Vinh (2012) áp dụng phương pháp lai (hybrid method) với kịch bản đứt gãy Sông Hồng ($M = 6{,}5$) tính ra biên độ gia tốc $A_{max}$ nằm ngang từ $150 \div 545\text{ cm/s}^2$; Nguyễn Công Thăng và Phạm Đình Nguyên (2012) chỉ ra rằng: "Điều kiện địa hình trên mặt đất là một yếu tố quan trọng có ảnh hưởng đến trường sóng địa chấn; Điều kiện địa hình trên mặt đất có xu hướng làm gia tăng biên độ dao động nền ứng với các thành phần nằm ngang tại các vị trí quanh đỉnh núi..."

Luận án của Giang Kiên Trung (2022) định vị chính xác khoảng trống học thuật: tích hợp số liệu địa chất công trình chi tiết với mô phỏng lan truyền sóng 1D phi tuyến dẻo đàn hồi qua chương trình DEEPSOIL, thiết lập chuỗi dữ liệu liên tục dọc mặt cắt Tây - Đông Hà Nội, cung cấp các bằng chứng định lượng chính xác về sự phân hóa trường gia tốc và làm sáng tỏ hiện tượng bất thường của tầng bùn yếu dày.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đã mở rộng và hoàn thiện lý thuyết truyền sóng địa chấn trong môi trường phân lớp đàn dẻo đối với trầm tích trẻ đồng bằng sông châu thổ:

  • Mở rộng lý thuyết phản ứng nền phi tuyến của Seed - Idriss: Luận án chứng minh trong điều kiện địa chất châu thổ Sông Hồng, khi tầng bùn sét bão hòa nước có vận tốc sóng ngang cực thấp ($V_s < 80\text{ m/s}$) và chiều dày vượt quá 30 m (tiêu biểu tại lỗ khoan LK50 Đê La Thành), hiện tượng suy giảm gia tốc đỉnh ($PGA$ giảm từ $0{,}25\text{ g}$ xuống $0{,}18\text{ g}$) xuất hiện đồng thời với sự dịch chuyển mạnh mẽ của chu kỳ trội $T_0$ lên vùng $1{,}0 \div 1{,}5\text{ s}$. Điều này thách thức quan niệm kinh nghiệm truyền thống cho rằng trầm tích yếu luôn khuếch đại biên độ gia tốc cực đại trong mọi điều kiện chấn động.
  • Xác lập mô hình lý thuyết quan hệ cấu trúc - vận tốc: Xây dựng khung chuyển đổi tối ưu giữa chỉ số xuyên tiêu chuẩn $N_{SPT}$ và vận tốc sóng ngang $V_s$, phù hợp với đặc tính cơ lý của các phân vị địa tầng Đệ tứ miền Bắc Việt Nam, bổ sung cho các công thức kinh nghiệm toàn cầu của Mayne (2001, 2007) và Abolfazl Eslami et al. (2019).
  • Hệ thống mệnh đề lý thuyết (Theoretical Propositions):
    • Mệnh đề 1: Hệ số khuếch đại gia tốc phổ ($RSRA$) tại đỉnh trầm tích tỷ lệ thuận với độ tương phản trở kháng âm học giữa tầng phủ Đệ tứ và tầng đá gốc Mesozoic ($Z = \rho_{bedrock} V_{s,bedrock} / \rho_{soil} V_{s,soil}$).
    • Mệnh đề 2: Chu kỳ dao động trội $T_0$ dọc theo tuyến địa chất Tây - Đông biến thiên theo hàm đơn điệu tăng tương ứng với độ dày tầng trầm tích $H$ theo nguyên lý gần đúng Rayleigh mở rộng ($T_0 \approx 4\sum H_i / V_{si}$).
    • Mệnh đề 3: Hiện tượng trễ pha và tiêu tán năng lượng do biến dạng dẻo trong tầng sét bùn dày $amQ$ đóng vai trò như một bộ lọc thông thấp (low-pass filter) tự nhiên đối với các thành phần tần số cao của sóng cắt $S$.
+-------------------------------------------------------------------------+
|                  KHUNG PHÂN TÍCH TỔNG THỂ CỦA LUẬN ÁN                   |
+-------------------------------------------------------------------------+
                                    |
                                    v
+-------------------------------------------------------------------------+
|              TẬP DỮ LIỆU ĐỊA CHẤT CÔNG TRÌNH VÀ THỦY VĂN                |
|  - Thu thập > 324 lỗ khoan tại Hà Nội; Lọc 53 lỗ khoan trắc tuyến       |
|  - Dữ liệu địa tầng: Phân lớp, mực nước ngầm, chỉ số N_SPT, độ dẻo (PI) |
+-------------------------------------------------------------------------+
                                    |
                                    v
+-------------------------------------------------------------------------+
|                CHUẨN HÓA THÔNG SỐ VẬT LÝ NỀN ĐẤT DỰNG 1D                |
|  - Tính toán vận tốc sóng ngang V_s từ N_SPT (theo Mayne & Sisman)       |
|  - Tính trọng lượng riêng gamma (theo Mayne, 2001)                      |
|  - Xây dựng mô-đun cắt động lực học G_max = rho * V_s^2                 |
|  - Xác lập đường cong suy giảm G/G_max và hệ số cản beta theo gamma_xy  |
+-------------------------------------------------------------------------+
                                    |
                                    v
+-------------------------------------------------------------------------+
|                   MÔ PHỎNG ĐỘNG LỰC HỌC VỚI DEEPSOIL                    |
|  - Băng gia tốc đầu vào: Imperial Valley-06 (khớp phổ Campbell 2008)    |
|  - Phân tích truyền sóng 1D phi tuyến (Non-linear time domain)          |
|  - Đánh giá phổ gia tốc mặt đất, hệ số khuếch đại RSRA, PGA, và T_0     |
+-------------------------------------------------------------------------+
                                    |
                                    v
+-------------------------------------------------------------------------+
|             KẾT QUẢ ĐỊNH LƯỢNG & ĐẶC TRƯNG HÓA TUYẾN TÂY - ĐÔNG         |
|  - Phía Tây (Đồng Mô, Quốc Oai): Trầm tích mỏng -> PGA cao, T_0 nhỏ     |
|  - Phía Đông (Long Biên, Thanh Trì): Trầm tích dày -> T_0 kéo dài       |
|  - Vị trí dị thường (LK50 Đê La Thành): Bùn dày > 30m -> PGA suy giảm   |
+-------------------------------------------------------------------------+

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp chặt chẽ ba trụ cột lý thuyết:

  1. Lý thuyết lan truyền sóng đàn hồi trong môi trường liên tục nhiều lớp (Wave Propagation Theory): Thiết lập phương trình vi phân chuyển động một chiều của sóng cắt lan truyền thẳng đứng từ nền đá gốc lên bề mặt tự do.
  2. Lý thuyết cơ học đất động lực học phi tuyến (Nonlinear Soil Dynamics): Áp dụng mô hình phần tử cấu tạo phụ thuộc biến dạng (strain-dependent modulus reduction and damping curves) của Darendeli (2001) và Hashash et al. (2016) được tích hợp trong DEEPSOIL.
  3. Lý thuyết mô hình chuyển động đất mạnh NGA (Next Generation Attenuation): Ứng dụng mô hình Campbell & Bozorgnia (2008) để kiểm soát các thông số nguồn chấn tiêu và xác lập băng gia tốc chuẩn đầu vào tại bề mặt đá gốc.

Phạm vi biên (boundary conditions) của khung phân tích được xác định rõ: áp dụng cho môi trường trầm tích Đệ tứ phân lớp nằm ngang trên nửa không gian đá gốc đàn hồi đẳng hướng, sóng tới là sóng khối cắt $S$ truyền thẳng đứng vuông góc với mặt phẳng phân lớp, bỏ qua tác động 3D phức tạp của địa hình bề mặt gồ ghề.

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

Thiết kế nghiên cứu

Luận án vận dụng quan điểm nhận thức luận thực chứng (positivism) với phương pháp tiếp cận định lượng suy diễn (quantitative deductive approach). Thiết kế nghiên cứu kết hợp đa tầng (multi-level design) bao gồm:

  • Tầng vĩ mô (Regional scale): Phân tích kiến tạo đới đứt gãy Sông Hồng, Sông Chay, Sông Lô, Sông Mã và thiết lập kịch bản nguồn động đất cực đại ($M = 6{,}5$, độ sâu chấn tiêu $h = 10\text{ km}$).
  • Tầng trung mô (Transect scale): Thiết lập mặt cắt trắc tuyến địa chất - địa vật lý Tây - Đông dài hơn 30 km đi qua các đơn vị hình thái địa mạo khác nhau của Hà Nội (từ Quốc Oai, Hoài Đức qua Cầu Giấy, Ba Đình, Đống Đa, Hoàn Kiếm sang Long Biên, Gia Lâm).
  • Tầng vi mô (Site-specific scale): Mô hình hóa cấu trúc cột địa tầng 1D chi tiết đến từng mét chiều sâu tại các lỗ khoan tiêu biểu.

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu tuân thủ các tiêu chuẩn kiểm chuẩn nghiêm ngặt:

  • Chiến lược chọn mẫu: Thu thập tổng số 324 lỗ khoan từ các nguồn lưu trữ địa chất công trình quốc gia, các dự án hạ tầng trọng điểm (Đường cao tốc Láng - Hòa Lạc, Cầu Nhật Tân, Cầu Thanh Trì, Tuyến đường sắt đô thị Nhổn - Ga Hà Nội, Vành đai 4, các khu đô thị Bắc Đại Kim, Cầu Bươu) và các báo cáo thăm dò nước dưới đất.
  • Tiêu chuẩn lựa chọn: 53 lỗ khoan nằm sát tuyến trắc cảnh Tây - Đông được chọn lọc với tiêu chí có chiều sâu khoan thăm dò tối thiểu chạm lớp sỏi cuội hoặc đá gốc ($H \ge 20 \div 130\text{ m}$), có mô tả địa tầng thạch học chi tiết và có đầy đủ chỉ số thí nghiệm xuyên tiêu chuẩn $N_{SPT}$.
  • Chuẩn hóa thông số vật lý:
    • Vận tốc sóng cắt $V_s$ được xác định thông qua tương quan thực nghiệm đã được kiểm chứng cho từng loại đất: $$\text{Đất dính (sét, sét pha): } V_s = 80{,}6 \cdot N_{SPT}^{0{,}331}$$ $$\text{Đất rời (cát, cát pha, sỏi cuội): } V_s = 97{,}2 \cdot N_{SPT}^{0{,}323}$$
    • Đơn vị trọng lượng riêng ($\gamma$, $\text{kN/m}^3$) được xác định theo mô hình tương quan của Mayne (2001): $$\gamma = 8{,}32 \cdot \log(V_s) - 1{,}61 \cdot \log(z)$$ với $z$ là độ sâu tính bằng mét, đảm bảo sự tương thích cao với số liệu thí nghiệm mẫu nguyên dạng trong phòng ($R^2 > 0{,}88$).
  • Độ tin cậy và kiểm chuẩn: Sử dụng phép tam giác giác hóa dữ liệu (data triangulation) đối chiếu giữa tài liệu khoan địa chất công trình, mặt cắt địa vật lý thăm dò điện lực/địa chấn nông và kết quả đo vi địa chấn $HVSR$ đã công bố của các tác giả tiền bối (Nguyễn Ánh Dương et al., 2017; Bùi Thị Nhung, 2018).

Data và phân tích

Mô hình hóa số được thực hiện bằng chương trình DEEPSOIL v7.0 (Board of Trustees of the University of Illinois). Các bước phân tích bao gồm:

  1. Xác lập tín hiệu kích động đầu vào: Sử dụng băng gia tốc Imperial Valley-06 (1979) thành phần nằm ngang ghi nhận trên nền đá gốc, được điều chỉnh tỷ lệ biên độ và phổ gia tốc để tương thích hoàn toàn với phổ phản ứng mục tiêu NGA Campbell & Bozorgnia (2008) cho kịch bản động đất cực đại đới Sông Hồng ($M = 6{,}5$, $R_{rup} = 10\text{ km}$, $PGA_{rock} \approx 0{,}20 \div 0{,}25\text{ g}$).
  2. Mô hình hóa động lực học 1D: Cột địa tầng được chia thành các lớp vi mô có chiều dày thỏa mãn điều kiện $h_i \le V_{si} / (4 \cdot f_{max})$ với $f_{max} = 25\text{ Hz}$ để ngăn ngừa sai số số học và phản xạ giả.
  3. Phương pháp phân tích: Thực hiện cả hai quy trình phân tích tuyến tính tương đương (Equivalent Linear - EQL) trong miền tần số và phân tích phi tuyến từng bước thời gian (Nonlinear - NL) với tiêu chuẩn bền dẻo Mohr-Coulomb kết hợp đường cong suy giảm mô-đun cắt $G/G_{max}$ và độ cản $\beta$ theo biến dạng trượt $\gamma_{xy}$.
  4. Kiểm tra tính vững (Robustness Checks): Tiến hành phân tích độ nhạy bằng cách thay đổi giả định cấu trúc địa tầng tại lỗ khoan dị thường LK50 (thay thế tầng bùn bằng tầng sét chặt, hoặc thay thế các tầng trên sỏi cuội bằng bùn mềm) để kiểm chứng cơ chế suy giảm gia tốc.

Phát hiện đột phá và implications

Những phát hiện then chốt

Luận án đã đạt được các phát hiện khoa học mang tính bước ngoặt, được tổng hợp qua hai luận điểm mới cốt lõi của công trình:

"Hiệu ứng nền do động đất gây ra trên tuyến Tây Đông là rõ rệt khi động đất mạnh xảy ra, thể hiện qua xu thế chung của việc tăng biên độ gia tốc và chu kỳ trội của rung động nền tại đầu phía Đông - nơi bề dày các lớp trầm tích yếu lớn hơn - có giá trị lớn hơn khi so với phía Tây - nơi bề dày các lớp trầm tích yếu nhỏ hơn." (Giang Kiên Trung, 2022, tr. 14)

"Tuy nhiên, với trường hợp cá biệt, tại vị trí có lớp đất yếu (lớp có vận tốc <80 m/s) trên mặt với bề dày lớn hơn 30 m, khi động đất xảy ra biên độ gia tốc dao động nền lại có xu hướng giảm, mặc dù chu kỳ trội của rung động nền vẫn có xu hướng tăng lên." (Giang Kiên Trung, 2022, tr. 14)

Cụ thể, bốn phát hiện chi tiết gồm:

  1. Quy luật phân hóa không gian phổ gia tốc Tây - Đông:
    • Phần phía Tây thành phố (từ Đồng Mô, Quốc Oai qua Hoài Đức; gồm các lỗ khoan tiêu biểu CDS-14, VBPC CL-LK3, HKC-CS04) có bề dày trầm tích Đệ tứ mỏng ($H \approx 20 \div 35\text{ m}$), đá gốc nông hoặc lộ diện, vận tốc $V_s$ trung bình cao ($V_{s30} > 250\text{ m/s}$, nền loại B và C). Kết quả mô phỏng cho thấy phổ gia tốc phản ứng có đỉnh khuếch đại vượt trội ở dải chu kỳ ngắn ($T = 0{,}1 \div 0{,}3\text{ s}$), giá trị $PGA$ tại bề mặt tăng từ $1{,}2 \div 1{,}6$ lần so với đá gốc.
    • Phần phía Đông thành phố (khu vực trũng Sông Hồng, Long Biên, Gia Lâm, Thanh Trì; gồm các lỗ khoan LK5DT, LK3P16, M2) có bề dày trầm tích Đệ tứ tăng mạnh ($H = 60 \div 90\text{ m}$), bao gồm nhiều tầng cát pha, bùn sét bão hòa nước xen kẽ ($V_{s30} < 180\text{ m/s}$, nền loại D và E). Chu kỳ trội của rung động nền dịch chuyển rõ rệt về phía chu kỳ dài ($T_0 = 0{,}6 \div 1{,}2\text{ s}$), với hệ số khuếch đại phổ gia tốc $RSRA$ đạt giá trị từ $2{,}0 \div 3{,}2$.
  2. Phát hiện dị thường tiêu tán phi tuyến tại lỗ khoan LK50 (Đê La Thành): Tại vị trí lỗ khoan LK50, địa tầng ghi nhận lớp bùn sét màu xám đen, trạng thái chảy - dẻo chảy thuộc tướng bồi tích lòng sông cổ/đầm lầy Sông Hồng ($amQ$) có vận tốc sóng cắt cực kỳ thấp ($V_s \approx 65 \div 78\text{ m/s}$) với bề dày tích lũy liên tục vượt quá 30 m. Kết quả mô phỏng động lực học chỉ ra hiện tượng phi tuyến đặc biệt: biên độ gia tốc nền cực đại ($PGA$) tại mặt đất giảm xuống chỉ còn $0{,}18\text{ g}$ (so với gia tốc nền kích động đầu vào $0{,}25\text{ g}$ tại đá gốc, tương đương giảm $28%$). Nguyên nhân vật lý là do biến dạng trượt trong tầng bùn vượt qua giới hạn đàn hồi ($\gamma_{xy} > 10^{-2}$), kích hoạt độ cản nhớt và tiêu tán dẻo cực lớn ($\beta > 15 \div 20%$), làm suy hao năng lượng sóng tần số cao. Tuy nhiên, chu kỳ trội $T_0$ tại vị trí này vẫn tăng vọt lên mức $T_0 \approx 1{,}15\text{ s}$.
  3. Mặt cắt phân tầng vận tốc địa chấn hoàn chỉnh đầu tiên: Luận án đã thành lập thành công mặt cắt vận tốc sóng cắt $V_s$ liên tục theo phương vị Tây - Đông cắt qua lòng chảo Hà Nội. Cấu trúc gồm 4 tầng địa chấn - địa kỹ thuật chính:
    • Tầng 1 (Lớp phủ bề mặt): Sét, sét pha, bùn hữu cơ ($V_s = 65 \div 150\text{ m/s}$, dày $2 \div 35\text{ m}$).
    • Tầng 2 (Trầm tích giữa): Cát mịn, cát trung, sét dẻo cứng ($V_s = 160 \div 280\text{ m/s}$, dày $10 \div 30\text{ m}$).
    • Tầng 3 (Trầm tích cuội sỏi cuội kết bở rời): Cuội sỏi laterit, sạn sỏi bão hòa nước chặt xít ($V_s = 320 \div 550\text{ m/s}$, dày $10 \div 25\text{ m}$).
    • Tầng 4 (Đá gốc kỹ thuật / Bedrock): Đá biến chất, cát kết, đá vôi nứt nẻ vừa đến nguyên khối ($V_s \ge 800\text{ m/s}$).
  4. Mối tương quan giữa cấu trúc địa chất sông cổ và nguy cơ cộng hưởng công trình: Đối chiếu với sơ đồ hệ thống dòng chảy cổ Sông Hồng (Nguyễn Xuân Đạo, 1984), các khu vực có chu kỳ trội kéo dài ($T_0 > 0{,}8\text{ s}$) trùng khớp hoàn toàn với các đới lòng sông cổ bị vùi lấp. Điều này tạo ra nguy cơ cộng hưởng phá hủy rất cao đối với các công trình nhà cao tầng từ 10 đến 25 tầng ($T_{structure} \approx 0{,}1 \cdot N \approx 1{,}0 \div 2{,}5\text{ s}$) tại khu vực nội đô lịch sử và phía Đông.
Vị trí lỗ khoan Tọa độ & Khu vực Bề dày trầm tích ($H$, m) Phân loại đất (TCVN 9386) $V_{s30}$ trung bình (m/s) Chu kỳ trội $T_0$ (s) $PGA$ bề mặt mô phỏng (g) Đặc trưng phổ phản ứng
CDS-14 Ngọc Bài, Quốc Oai 34,2 B - C 275 0,22 0,32 Khuếch đại mạnh ở chu kỳ ngắn ($T < 0{,}3\text{ s}$)
HKC-CS04 Kênh Đồng Mô 20,0 B 310 0,18 0,35 Tần số cao chiếm ưu thế, biên độ đỉnh lớn
LK-P3 Hoàng Minh Giám 46,0 C - D 210 0,55 0,28 Phổ dẹt, dải cộng hưởng mở rộng $0{,}4 \div 0{,}7\text{ s}$
LK50 Đê La Thành 86,0 (Bùn > 30m) E (Đất đặc biệt) 125 1,15 0,18 (Suy giảm) Tiêu tán tần số cao, chu kỳ trội kéo dài
LK5DT Chương Dương Độ 80,0 D 165 0,85 0,29 Khuếch đại mạnh dải chu kỳ dài $0{,}7 \div 1{,}0\text{ s}$
LK3P16 Thanh Trì 68,0 D - E 150 0,92 0,27 Khuếch đại trầm tích sông, chu kỳ trội lớn
M2 Cổ Linh, Long Biên 55,0 D 170 0,80 0,30 Cộng hưởng rõ rệt tại dải tần số $1{,}2 \div 1{,}5\text{ Hz}$

Implications đa chiều

  • Ý nghĩa lý thuyết: Bổ sung cơ chế vật lý chi tiết về hiện tượng bão hòa biến dạng và tự tiêu tán động năng của tầng bùn dày trong mô hình lan truyền sóng 1D, giúp hoàn thiện các hướng dẫn phân loại nền đất đặc biệt trong cơ học địa chấn.
  • Đổi mới phương pháp luận: Thiết lập quy trình chuẩn hóa từ dữ liệu khảo sát địa kỹ thuật công trình thông thường ($N_{SPT}$) sang mô hình động lực học phi tuyến DEEPSOIL, có thể chuyển giao áp dụng trực tiếp cho các đô thị đồng bằng châu thổ khác tại Việt Nam như Hải Phòng, Nam Định, Cần Thơ, TP. Hồ Chí Minh.
  • Ứng dụng thực tiễn trong kỹ thuật xây dựng:
    • Cung cấp cơ sở trực tiếp để lựa chọn giải pháp móng: Đối với các công trình thấp tầng (< 5 tầng) tại phía Tây, cần chú trọng thiết kế kháng chấn chu kỳ ngắn ($T < 0{,}3\text{ s}$); đối với các công trình cao tầng (> 15 tầng) tại phía Đông và vùng đệm Sông Hồng, bắt buộc phải cắm cọc sâu qua tầng bùn sét $amQ$ vào tầng cuội sỏi chặt hoặc đá gốc ($V_s > 500\text{ m/s}$) để triệt tiêu hiệu ứng cộng hưởng dao động chu kỳ dài.
  • Khuyến nghị chính sách: Kiến nghị Bộ Xây dựng và Viện Khoa học Công nghệ Xây dựng (IBST) cập nhật phụ lục bản đồ phân vùng gia tốc nền và hệ số nền $S$ trong tiêu chuẩn quốc gia TCVN 9386:2012 cho khu vực Hà Nội mở rộng, phân tách rõ các phân vùng có nguy cơ suy giảm hoặc khuếch đại dị thường.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án thẳng thắn chỉ ra bốn giới hạn nghiên cứu (limitations):

  1. Giới hạn không gian của mô hình 1D: Mô hình truyền sóng một chiều chưa thể tích hợp đầy đủ các hiệu ứng tán xạ rìa bồn trũng (basin-edge effects), hội tụ năng lượng sóng bề mặt (Rayleigh và Love waves) phát sinh do sự bất đồng nhất hai bên mép bồn trầm tích như các mô hình 2D/3D phức tạp (Bard & Riepl-Thomas, 2000; Kawase, 1996).
  2. Nguồn dữ liệu băng gia tốc đầu vào: Do miền Bắc Việt Nam chưa có mạng lưới trạm quan trắc ghi nhận động đất mạnh dày đặc, nghiên cứu bắt buộc phải sử dụng băng gia tốc từ trận động đất quốc tế (Imperial Valley-06) đã qua chuẩn hóa phổ NGA thay vì số liệu thực đo tại chỗ.
  3. Độ sâu thăm dò của lỗ khoan: Một số lỗ khoan địa chất công trình chỉ đạt độ sâu $30 \div 50\text{ m}$, chưa chạm tới nền đá gốc tuyệt đối ($V_s \ge 800\text{ m/s}$), đòi hỏi phải ngoại suy vận tốc cho các tầng nằm dưới dựa trên tài liệu địa vật lý khu vực.
  4. Mật độ hố khoan phân bố chưa đồng đều: Tuyến nghiên cứu tập trung dọc các trục giao thông chính (Láng - Hòa Lạc, Đại lộ Thăng Long, các cầu qua Sông Hồng), một số khoảng cách giữa các hố khoan tại khu vực chuyển tiếp còn thưa.

Chương trình nghiên cứu tương lai (Future Research Agenda) mở ra 4 hướng tiếp nối:

  • Hướng 1: Phát triển mô hình mô phỏng số 3D đàn dẻo liên tục cho toàn bộ bồn trũng Đệ tứ Hà Nội, tích hợp ảnh hưởng của góc dốc nền đá gốc và cấu trúc đứt gãy sinh chấn Sông Hồng - Sông Chay.
  • Hướng 2: Triển khai mạng lưới quan trắc gia tốc mạnh nội đô theo thời gian thực (tương tự hệ thống P-Alert tại Đài Loan của Wu et al., 2019) để thu thập dữ liệu sóng động đất thực tế phục vụ hiệu chỉnh mô hình.
  • Hướng 3: Đánh giá khả năng hóa lỏng đất cát bão hòa nước ($amQ$) tại các quận ven Sông Hồng dưới tác động của các kịch bản động đất kích động mạnh.
  • Hướng 4: Nghiên cứu tương tác động lực học phi tuyến giữa kết cấu công trình ngầm (tuyến metro ngầm Nhổn - Ga Hà Nội) và môi trường đất nền phân lớp phức tạp.

Tác động và ảnh hưởng

  • Tác động học thuật: Công trình là tài liệu tham khảo mẫu mực trong đào tạo sau đại học ngành Vật lí Địa cầu, Địa kỹ thuật và Xây dựng công trình. Với 3 công trình công bố (1 bài báo quốc tế uy tín và 2 bài báo tạp chí chuyên ngành quốc gia), tiềm năng trích dẫn ước tính đạt hàng trăm lượt trong các nghiên cứu vi phân vùng chấn động và đánh giá rủi ro thiên tai đô thị tại Đông Nam Á.
  • Chuyển đổi công nghiệp xây dựng: Các tập đoàn tư vấn thiết kế, nhà thầu xây dựng hạ tầng ngầm và cao ốc tại Hà Nội có thể sử dụng trực tiếp bộ thông số phổ phản ứng $RSRA$, $V_{s30}$, và chiều sâu tầng bùn để tối ưu hóa chi phí xử lý nền móng, tiết kiệm hàng triệu USD chi phí kết cấu trong khi vẫn đảm bảo an toàn chịu chấn.
  • Ảnh hưởng chính sách đô thị: Cung cấp luận cứ khoa học định lượng vững chắc cho UBND Thành phố Hà Nội, Sở Quy hoạch - Kiến trúc, Sở Xây dựng trong việc phân vùng rủi ro thảm họa, lập quy hoạch các đô thị vệ tinh (Hòa Lạc, Xuân Mai, Sóc Sơn, Đông Anh) và hoàn thiện quy chuẩn xây dựng đô thị ngầm.
  • Giá trị xã hội và quốc tế: Giảm thiểu thiệt hại về người và tài sản trước nguy cơ thảm họa động đất kích hoạt từ đứt gãy Sông Hồng; đóng góp dữ liệu đối sánh quan trọng cho cộng đồng nghiên cứu địa chấn quốc tế về hành vi của trầm tích châu thổ sông nhiệt đới.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và học viên sau đại học: Tiếp cận phương pháp luận chuẩn mực về tích hợp cơ sở dữ liệu lỗ khoan địa chất công trình với công cụ mô phỏng số DEEPSOIL và lý thuyết truyền sóng địa chấn 1D.
  • Các nhà khoa học và chuyên gia địa chấn: Sử dụng bộ số liệu mặt cắt vận tốc sóng cắt $V_s$ Tây - Đông đã được số hóa làm nền tảng để phát triển các bài toán truyền sóng 2D/3D phân giải cao.
  • Kỹ sư kết cấu và địa kỹ thuật (R&D công nghiệp): Nhận được các thông số phổ thiết kế gia tốc thực tế tại các quận nội thành và ngoại thành, giúp thiết kế hệ thống móng cọc kháng chấn chính xác, tránh hiện tượng cộng hưởng phá hủy.
  • Cơ quan quản lý nhà nước và lập chính sách: Có cơ sở khoa học để rà soát quy chuẩn TCVN 9386, lập bản đồ rủi ro thiên tai chi tiết cấp phường/quận và xây dựng kế hoạch ứng phó khẩn cấp cho thủ đô Hà Nội.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?
Trả lời: Đóng góp độc đáo nhất là phát hiện và lượng hóa cơ chế suy giảm gia tốc nền cực đại ($PGA$) do hiệu ứng lọc phi tuyến tiêu tán năng lượng của tầng bùn sét dày > 30 m ($V_s < 80\text{ m/s}$) tại Hà Nội. Kết quả này mở rộng lý thuyết phản ứng nền phi tuyến của Seed & Idriss (1970) và mô hình Darendeli (2001) trong điều kiện trầm tích bão hòa nước châu thổ Đệ tứ, làm sáng tỏ rằng trầm tích mềm không chỉ đơn thuần khuếch đại biên độ sóng mà còn có thể làm tiêu hao biên độ đỉnh trong khi dịch chuyển năng lượng sang chu kỳ rất dài ($T > 1{,}0\text{ s}$).

2. Đổi mới phương pháp luận của luận án so với các công trình trước đây tại Việt Nam?
Trả lời: So với phương pháp độ cứng địa chấn kinh nghiệm của Medvedev áp dụng năm 1964, 1990 (Viện Vật lí Địa cầu, 1990) và phương pháp vi địa chấn tuyến tính thuần túy của Nguyễn Tiến Hùng & Kuo-Liang Wen (2011) hay Nguyễn Ánh Dương et al. (2017), luận án của Giang Kiên Trung (2022) là công trình đầu tiên tại Việt Nam xây dựng mặt cắt trắc tuyến $V_s$ liên tục từ Tây sang Đông và giải bài toán truyền sóng động lực học phi tuyến 1D qua phần mềm DEEPSOIL v7.0 với dữ liệu chuẩn hóa từ 53 lỗ khoan thực nghiệm sâu tới 130 m, đồng thời kiểm soát phổ đầu vào bằng mô hình NGA Campbell & Bozorgnia (2008).

3. Phát hiện gây bất ngờ nhất từ tập dữ liệu thực nghiệm là gì?
Trả lời: Phát hiện bất ngờ nhất xuất hiện tại lỗ khoan LK50 (Đê La Thành), nơi có lớp bùn sét đầm lầy cổ dày trên 30 m. Trong khi xu thế chung dọc tuyến Tây - Đông là gia tốc nền cực đại $PGA$ tăng dần về phía Đông, tại LK50 giá trị $PGA$ bề mặt lại giảm sút mạnh ($0{,}18\text{ g}$ so với $0{,}25\text{ g}$ tại nền đá gốc), đi kèm sự biến dạng dẻo lớn của tầng đất. Thử nghiệm thay đổi giả định địa tầng trong luận án đã chứng minh chính lớp bùn xốp dày là nguyên nhân trực tiếp triệt tiêu gia tốc tần số cao.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?
Trả lời: Hoàn toàn có. Luận án cung cấp chi tiết toàn bộ tọa độ, độ sâu, phân tầng địa chất của 324 lỗ khoan (Phụ lục A và B), công thức tính đổi tương quan $N_{SPT} \rightarrow V_s$, công thức xác định trọng lượng riêng $\gamma$ theo Mayne (2001), các thông số đường cong phi tuyến dẻo đàn hồi, và tập tin băng gia tốc kích động đầu vào Imperial Valley-06. Mọi nhà nghiên cứu độc lập đều có thể tái lập chính xác mô phỏng trên DEEPSOIL hoặc các phần mềm tương đương (SHAKE, STrata).

5. Chương trình nghiên cứu 10 năm tới được phác thảo ra sao?
Trả lời: Chương trình 10 năm bao gồm: (1) Mở rộng mặt cắt trắc tuyến theo hướng Bắc - Nam cắt qua đứt gãy Sông Chay và đứt gãy Sông Lô; (2) Nâng cấp mô hình số 1D lên không gian 3D đàn dẻo cho toàn bồn trũng Hà Nội; (3) Tích hợp quan trắc vi địa chấn mảng (SPAC array) đa kênh; (4) Xây dựng hệ thống bản đồ số vi phân vùng gia tốc phổ tương tác phục vụ cấp phép xây dựng thông minh (Smart City Seismic Resilience).

Kết luận

Luận án tiến sĩ của Giang Kiên Trung (2022) đã khẳng định vị thế một công trình khoa học mẫu mực, giải quyết trọn vẹn mục tiêu định lượng hóa hiệu ứng nền địa phương tại thủ đô Hà Nội. Năm đóng góp cụ thể nổi bật bao gồm:

  1. Xây dựng thành công cơ sở dữ liệu địa kỹ thuật động lực học chuẩn hóa từ hơn 324 lỗ khoan, thiết lập mặt cắt vận tốc sóng cắt địa chấn ($V_s$) liên tục đầu tiên theo phương vị Tây - Đông thành phố Hà Nội.
  2. Ứng dụng thành công chương trình mô phỏng số phi tuyến DEEPSOIL v7.0 kết hợp phổ kích động chuẩn NGA, lượng hóa chính xác các tham số dao động nền ($PGA$, $RSRA$, $T_0$) phục vụ thiết kế kháng chấn theo TCVN 9386:2012.
  3. Xác lập quy luật biến thiên không gian: khu vực phía Tây có biên độ gia tốc ngắn chiếm ưu thế ($T < 0{,}3\text{ s}$), khu vực phía Đông có chu kỳ trội kéo dài ($T_0 = 0{,}6 \div 1{,}2\text{ s}$) do bề dày trầm tích Đệ tứ tăng cao.
  4. Phát hiện và lý giải khoa học hiện tượng suy giảm gia tốc dị thường tại các vùng chứa tầng bùn sét yếu dày trên 30 m ($V_s < 80\text{ m/s}$) do cơ chế tiêu tán năng lượng phi tuyến.
  5. Cung cấp luận cứ khoa học thực tiễn vững chắc phục vụ công tác quy hoạch không gian ngầm, định hướng kết cấu nền móng nhà cao tầng và giảm thiểu hiểm họa địa chấn đô thị.

Công trình tạo ra bước tiến quan trọng về phương pháp luận trong vật lý địa cầu ứng dụng tại Việt Nam, mở ra ba hướng nghiên cứu chuyên sâu về mô phỏng 3D bồn trũng châu thổ, phân tích rủi ro hóa lỏng đất nền và thiết lập hệ thống quan trắc cảnh báo sớm động đất đô thị trong kỷ nguyên mới.