Tổng quan về luận án

Công trình ngầm đô thị, đặc biệt là các tuyến đường sắt đô thị (Metro) và hệ thống công trình ngầm phục vụ quốc phòng - an ninh, giữ vai trò chiến lược trong cơ sở hạ tầng của các đô thị hiện đại. Tuy nhiên, tính phức tạp trong ứng xử cơ học của môi trường đất đá tương tác với kết cấu vỏ hầm dưới tác động của tải trọng động lực học biến đổi phi tuyến nhanh như động đất đòi hỏi các phương pháp phân tích tiên tiến. Tại Việt Nam, mạng lưới trạm quan trắc địa chấn ghi nhận dữ liệu thực tế về các trận động đất mạnh còn rất hạn chế, dẫn đến sự thiếu hụt nghiêm trọng nguồn dữ liệu giản đồ gia tốc nền phục vụ tính toán theo lịch sử thời gian (Time-history analysis).

Khoảng trống nghiên cứu (Research gap) then chốt mà luận án giải quyết chính là sự thiếu vắng công cụ toán - tin chuyên dụng và quy trình kỹ thuật chuẩn hóa để phát sinh giản đồ gia tốc nhân tạo đáp ứng tiêu chuẩn kháng chấn Việt Nam TCVN 9386-2012, tích hợp đồng thời các đặc trưng cơ học của nguồn chấn, đường truyền và tương tác đất - vỏ hầm có xét đến tính nửa cứng của khớp nối phân tố (Semi-rigid joints).

Luận án xác lập 03 câu hỏi nghiên cứu và giả thuyết khoa học cốt lõi:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để điều chỉnh các giản đồ gia tốc thực tế nhằm thỏa mãn điều kiện khớp phổ phản ứng mục tiêu theo TCVN 9386-2012 mà không làm biến dạng hàm năng lượng tích lũy Arias ($I_A$) và đặc tính không dừng của sóng địa chấn?
    • Giả thuyết 1 (H1): Ứng dụng biến đổi Wavelet với hàm cơ sở Cosine cải biên (Modified Cosine Wavelet) cho phép bổ sung các lượng hiệu chỉnh cục bộ trong miền thời gian, tối ưu hóa sai số phổ $\Delta S_a \to 0$ với tốc độ hội tụ nhanh hơn phương pháp điều chỉnh miền tần số truyền thống.
  • Câu hỏi nghiên cứu 2 (RQ2): Có thể mô phỏng ngẫu nhiên giản đồ gia tốc nền đại diện cho các kịch bản động đất cụ thể tại khu vực Hà Nội (như đới đứt gãy Sông Hồng - Sông Chảy) dựa trên các đặc trưng xác suất của gói Wavelet hay không?
    • Giả thuyết 2 (H2): Hệ phương trình hồi quy 13 tham số của gói Wavelet (Wavelet Packet Transform - WPT) kết hợp biến đổi ngược (iWPT) có khả năng tái lập các giản đồ gia tốc nhân tạo mang bản chất ngẫu nhiên của nguồn phát ($M_w$), tiêu cự ($R_{hyp}$), tâm cự ($R_{rup}$) và vận tốc truyền sóng cắt trung bình $V_{s,30}$.
  • Câu hỏi nghiên cứu 3 (RQ3): Mức độ ảnh hưởng của dạng kích thích gia tốc nhân tạo và tính chất liên kết nửa cứng giữa các phân tố vỏ hầm lắp ghép (Segmental lining) đến nội lực cực đại ($M_{max}, Q_{max}, N_{max}$) trong kết cấu công trình ngầm tại Hà Nội là bao nhiêu?
    • Giả thuyết 3 (H3): Phân tích động lực học trực tiếp trong Plaxis 2D với mô hình nền Hardening Soil (HS) và liên kết Jassen sẽ phản ánh sự phân bố lại ứng suất chính xác hơn, chứng minh rằng liên kết nửa cứng giúp giảm đáng kể mômen uốn cực đại so với giả thiết vỏ hầm cứng liên tục tuyệt đối.

Khung lý thuyết của nghiên cứu được xây dựng trên sự giao thoa giữa Lý thuyết Biến đổi Wavelet (Wavelet Theory), Cơ học Đất động lực học (Soil Dynamics) và Động lực học Công trình Ngầm (Dynamics of Underground Structures). Luận án khảo sát đối tượng thực tế là hầm ngầm thuộc tuyến đường sắt đô thị số 3 (Nhổn - Ga Hà Nội), vị trí trọng điểm tại quận Ba Đình (tọa độ $21{,}030^\circ\text{N}; 105{,}824^\circ\text{E}$), khai thác dữ liệu đới động đất Sông Hồng - Sông Chảy ($M_{w,max} = 6{,}5$; $R_{rup} = 15\text{ km}$; $R_{hyp} = 15\text{ km}$).

Literature Review và Positioning

Nghiên cứu cơ học công trình ngầm chịu tác động động đất trên thế giới đã phát triển qua ba trường phái phân tích cơ bản:

  1. Trường phái phân tích tĩnh áp đặt chuyển vị biên (Impose Seismic Ground Deformation - ISGD): Tiêu biểu với các lời giải giải tích của Wang (1993), Penzien và Wu (1998), dựa trên giả thiết môi trường đàn hồi đồng nhất, xác định biến dạng trượt cực đại $\gamma_{max} = \frac{v_{max}}{V_{se}}$ (với $v_{max}$ là vận tốc cực đại mặt đất, $V_{se}$ là vận tốc truyền hiệu quả sóng cắt) rồi áp đặt chuyển vị cưỡng bức lên biên mô hình. Nhược điểm lớn của ISGD là chỉ xem xét đến gia tốc đỉnh (PGA), bỏ qua hoàn toàn diễn biến lịch sử thời gian và quán tính của hệ.
  2. Trường phái phản lực tương hỗ đàn hồi (Hyperstatic Reaction Method - HRM): Được khởi xướng bởi Duddeck & Erdmann (1985), Leca & Clough (1992), Oreste (2007) và được Đỗ Ngọc Anh (2014) phát triển cho vỏ hầm lắp ghép. HRM tách kết cấu khỏi môi trường, thay thế đất đá bằng các liên kết lò xo Winkler phi tuyến. Phương pháp này tính toán nhanh trên mô hình không gian nhưng đơn giản hóa cấu trúc địa tầng nhiều lớp.
  3. Trường phái phân tích động lực học lịch sử thời gian (Dynamic Time-History Analysis): Hashash et al. (2001) tổng hợp các phân tích phần tử hữu hạn (FEM) và sai phân hữu hạn (FDM) với các phần mềm chuyên dụng như FLAC3D, PLAXIS, SASSI. Tại Việt Nam, Mai Đức Minh (2010), Lê Bảo Quốc (2013), Trần Thị Thu Hằng (2016) đã ứng dụng mô hình số và thí nghiệm bàn rung để khảo sát phản ứng động của công trình ngầm, khẳng định tầm quan trọng của việc nạp sóng gia tốc chân thực tại nền đá gốc (bedrock).

Về lĩnh vực phát sinh giản đồ gia tốc nhân tạo, hai luồng quan điểm đối lập định hình bức tranh nghiên cứu:

  • Hiệu chỉnh trên miền tần số (Frequency-domain adjustment): Fahjan & Ozdemir (2008), Nghiêm Mạnh Hiến (2012) dùng biến đổi Fourier nhanh (FFT) điều chỉnh phổ biên độ Fourier và giữ nguyên pha. Phương pháp này dễ gây xáo trộn quy luật phân bố năng lượng theo thời gian của bản ghi thực.
  • Hiệu chỉnh trên miền thời gian (Time-domain adjustment): Khởi xướng bởi Lilhanand & Tseng (1987, 1988) trên phần mềm RspMatch và được Hancock (2006) hoàn thiện với hàm Cosine cải biên trên SeismoMatch. Hướng tiếp cận này bảo tồn cấu trúc pha và năng lượng tích lũy Arias $I_A(t)$, cho phép sóng nhân tạo duy trì tính chất vật lý của sóng địa chấn gốc.

Bên cạnh đó, hướng mô phỏng ngẫu nhiên phi dừng chứng kiến bước tiến từ mô hình phân tích độ lệch pha Fourier (Thrainsson & Kiremidjian, 2002), mô hình lọc nhiễu trắng (Rezaeian, 2008) đến mô hình gói Wavelet (Wavelet Packet Transform) của Yamamoto (2011).

Luận án định vị chính xác khoảng trống học thuật: Các phần mềm quốc tế như SeismoMatch hay RspMatch chỉ hỗ trợ các tiêu chuẩn Eurocode 8, AS 1170 mà không tích hợp phổ TCVN 9386-2012; trong khi mô hình của Yamamoto (2011) tạo sóng ngẫu nhiên theo hệ phương trình hồi quy NGA nhưng chưa có cơ chế cưỡng bức khớp phổ thiết kế. Nghiên cứu của tác giả Vũ Ngọc Anh đã giải quyết hoàn chỉnh sự ngắt quãng này bằng việc tự thiết lập hai chương trình PG01 và PG02, tạo nên quy trình khép kín từ tạo sóng đến mô phỏng số tương tác đất - vỏ hầm có khớp nối nửa cứng.

Đó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 lý thuyết biến đổi Wavelet trong kỹ thuật địa chấn công trình thông qua việc chứng minh tính khả thi của việc tích hợp hàm cơ sở Wavelet Cosine cải biên vào bài toán tối ưu hóa sai số phổ phản ứng đàn hồi theo tiêu chuẩn Việt Nam.

Nghiên cứu phát triển mô hình toán học giải phương trình vi phân dao động của hệ một bậc tự do (1BTD) chịu kích thích bởi lượng bù gia tốc Wavelet $\delta a(t)$: $$\delta Y_i = \int_0^{t_i} \delta a(t) h_i(t_i - t) dt = \sum_{j=1}^{n_{CK}} b_j \int_0^{t_i} \psi_j(t) h_i(t_i - t) dt$$ với $h_i(t) = -\frac{\omega_i}{\sqrt{1-\xi^2}} e^{-\xi\omega_i t} \sin\left(\omega_i \sqrt{1-\xi^2} t\right)$ là phản ứng xung gia tốc của hệ có chu kỳ $T_i$, tần số góc $\omega_i$, và hệ số cản $\xi$.

Hệ phương trình đại số tuyến tính xác định véc-tơ trọng số ${b}$ được thiết lập chặt chẽ: $${\Delta S_a} = [S] {b} \implies {b} = [S]^{-1} {\Delta S_a}$$ trong đó $s_{ij} = \int_0^{t_i} \psi_j(t) h_i(t_i - t) dt$ là ma trận độ nhạy phản ứng gia tốc đơn vị.

Mô hình lý thuyết đề xuất 03 mệnh đề khoa học:

  • Mệnh đề 1 (Proposition 1): Phép hiệu chỉnh sơ bộ theo điều kiện bình phương tối thiểu khoảng lệch phổ $\sum (\Delta S_a(i))^2 \to \min$ với hệ số tỷ lệ $k_{hc}$ đóng vai trò điều kiện biên tối ưu hóa, giảm số vòng lặp Wavelet đi ít nhất $40%$.
  • Mệnh đề 2 (Proposition 2): Biến đổi ngược gói Wavelet (iWPT) với 13 tham số thống kê hồi quy từ năng lượng Arias $I_A$, thời gian duy trì $t_{5-95}$, trọng tâm thời gian $E_t$, tần số $E_f$ cho phép tái hiện tính chất ngẫu nhiên bất định của nguồn chấn mà phương pháp xác định (deterministic) không thể đạt được.
  • Mệnh đề 3 (Proposition 3): Sự suy giảm độ cứng chống uốn tại các khớp nối nửa cứng giữa các tấm vỏ hầm (Segmental lining joints) đóng vai trò bộ tiêu tán năng lượng biến dạng cục bộ, làm thay đổi dạng dao động riêng của kết cấu hầm trong môi trường đất yếu.

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

Khung phân tích của luận án tích hợp đa tầng 03 trụ cột lý thuyết: (1) Lý thuyết Xử lý Tín hiệu Số Phi dừng (Mallat 1999, Hancock 2006); (2) Lý thuyết Nguồn chấn & Năng lượng Địa chấn (Hanks & Kanamori 1979, Kramer 1996); và (3) Lý thuyết Động lực học Môi trường Đất nền Đa pha Phi tuyến (Schanz et al., Mô hình Hardening Soil).

   +-------------------------------------------------------------------------+
   |                  ĐẶC TRƯNG ĐỊA CHẤN KHU VỰC HÀ NỘI                       |
   |   (TCVN 9386-2012: agR, Phổ Sa_TK | Đới Sông Hồng - Sông Chảy: Mw, Rhyp) |
   +------------------------------------+------------------------------------+
                                        |
                 +----------------------+----------------------+
                 |                                             |
                 v                                             v
   +---------------------------+                 +---------------------------+
   |     CHƯƠNG TRÌNH PG01     |                 |     CHƯƠNG TRÌNH PG02     |
   | (Khớp phổ - Sóng nhân tạo)|                 | (Hồi quy WPT - Ngẫu nhiên)|
   |  - Sóng gốc: Điện Biên,   |                 |  - 13 tham số gói Wavelet |
   |    Campano Lucano, Lancang|                 |  - Biến đổi ngược iWPT    |
   |  - Wavelet Cosine cải biên|                 |  - Khớp đới đứt gãy Hà Nội|
   +-------------+-------------+                 +-------------+-------------+
                 |                                             |
                 +----------------------+----------------------+
                                        |
                                        v
                 +---------------------------------------------+
                 |        GIẢN ĐỒ GIA TỐC NỀN ĐÁ GỐC           |
                 |     (a(t), PGA, PGV, PGD, IA, t5-95, aRMS)  |
                 +----------------------+----------------------+
                                        |
                                        v
                 +---------------------------------------------+
                 |    MÔ HÌNH PHẦN TỬ HỮU HẠN PLAXIS 2D        |
                 |  - Nền đất: Hardening Soil (HS)             |
                 |  - Kết cấu: Vỏ hầm lắp ghép Metro Tuyến 3   |
                 |  - Liên kết nửa cứng: Mô hình Jassen        |
                 |  - Biên hấp thụ sóng: Viscous boundaries    |
                 +----------------------+----------------------+
                                        |
                                        v
                 +---------------------------------------------+
                 |    PHÂN TÍCH ỨNG XỬ ĐỘNG LỰC HỌC VỎ HẦM     |
                 |      (Nội lực cực đại Mmax, Qmax, Nmax      |
                 |       Phân tích hồi quy đa biến & Đánh giá) |
                 +---------------------------------------------+

Khung phân tích xác lập rõ các điều kiện biên: Môi trường địa chất giả định theo sơ đồ biến dạng phẳng (Plane Strain); đáy mô hình đặt tại tầng đá gốc chịu kích thích sóng ngang cắt truyền thẳng đứng; mặt tiếp xúc giữa vỏ hầm và đất được gán hệ số tương tác trượt $R_{inter} = 0{,}7$.

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

Thiết kế nghiên cứu

Luận án theo đuổi thế giới quan thực chứng (Positivism) kết hợp phương pháp luận cơ học tính toán định lượng nghiêm ngặt (Quantitative Computational Mechanics). Thiết kế nghiên cứu tiếp cận đa cấp độ (Multi-level modeling):

  • Cấp độ 1 (Vĩ mô - Địa chấn học): Xác định phổ phản ứng gia tốc thiết kế mục tiêu $S_{a_TK}(T)$ theo TCVN 9386-2012 ứng với gia tốc đỉnh tham chiếu $a_{gR} = 0{,}1098\text{ g}$ tại quận Ba Đình, hệ số tầm quan trọng $\gamma_I = 1{,}0$, loại nền đất $D$, các chu kỳ đặc trưng $T_B = 0{,}2\text{ s}$, $T_C = 0{,}8\text{ s}$, $T_D = 2{,}0\text{ s}$, hệ số khuếch đại đất nền $S = 1{,}35$.
  • Cấp độ 2 (Trung gian - Xử lý tín hiệu sóng): Phát triển thuật toán và lập trình thành công hai bộ công cụ tính toán:
    1. Chương trình PG01 (Matlab): Tự động hóa thuật toán Hancock, tích hợp biến đổi Wavelet Cosine cải biên: $$\psi_j(t) = \cos\left(2\pi f_j (t - t_j) + \gamma_j\right) \exp\left(-\frac{(t - t_j)^2}{2 \sigma_j^2}\right)$$ với độ lệch chuẩn $\sigma_j$ và độ dịch pha $\gamma_j$ phụ thuộc tần số $f_j$.
    2. Chương trình PG02 (Matlab): Cải biên thuật toán Yamamoto (2011), sử dụng hệ phương trình hồi quy 13 biến số để phát sinh ma trận hệ số gói Wavelet $c_{ij,k}$ trên lưới thời gian - tần số, sau đó áp dụng phép biến đổi ngược iWPT: $$x(t) = \sum_{i=1}^{2^j} \sum_{k=1}^{2^{N-j}} c_{ij,k} \psi_{ij,k}(t)$$
  • Cấp độ 3 (Vi mô - Cơ học công trình ngầm): Mô hình hóa chi tiết kết cấu hầm tròn đường kính ngoài $D = 6{,}3\text{ m}$, chiều dày vỏ hầm $t = 0{,}3\text{ m}$, bê tông C45/55 ($E_c = 36\text{ GPa}$, trọng lượng riêng $\gamma_b = 25\text{ kN/m}^3$) tương tác với địa tầng trầm tích Đệ Tứ Hà Nội bằng phần mềm PLAXIS 2D.

Quy trình nghiên cứu rigorous

Quy trình phát sinh và kiểm định sóng nhân tạo được kiểm soát sai số tuyệt đối. Với PG01, tiêu chí dừng lặp dựa trên sai số khoảng lệch phổ cực đại: $$\text{MaxError} = \max_i \left| \frac{S_a(T_i) - S_{a_TK}(T_i)}{S_{a_TK}(T_i)} \right| \le 5%$$

Ba giản đồ gia tốc thực tế đầu vào được lựa chọn đại diện bao gồm:

  1. Trận động đất Điện Biên (Việt Nam, 19/02/2001, $M_s = 5{,}3$, trạm Điện Biên);
  2. Trận động đất Campano Lucano (Italy, 23/11/1980, trạm 321 Bagnoli Irpino);
  3. Trận động đất Lancang (Trung Quốc, 06/11/1988, băng ghi ca064).

Sau khi xử lý qua PG01, thu được ba giản đồ gia tốc nhân tạo chuẩn hóa: BaDinh_01A, BaDinh_02A, và BaDinh_03A. Luận án đã kiểm chứng chéo kết quả của PG01 với phần mềm thương phẩm quốc tế SeismoMatch, cho thấy sự trùng khớp hoàn toàn về phổ phản ứng và đường cong tích lũy năng lượng Husid: $$h(t) = \frac{I_A(t)}{I_A(t_0)} = \frac{\int_0^t a^2(t) dt}{\int_0^{t_0} a^2(t) dt}$$ Đồng thời, phương pháp bảo tồn chính xác khoảng thời gian duy trì dao động mạnh $t_{5-95} = t_{95} - t_5$.

Data và phân tích

Địa tầng khu vực khảo sát tại Ba Đình được chia thành 6 lớp đất chính, khai báo bằng mô hình Hardening Soil (HS) để phản ánh tính chất phụ thuộc ứng suất và suy giảm độ cứng khi biến dạng:

  • Lớp 1: Đất lấp sét pha ($h = 2{,}5\text{ m}$, $\gamma = 18\text{ kN/m}^3$, $c = 12\text{ kPa}$, $\varphi = 15^\circ$);
  • Lớp 2: Sét dẻo mềm ($h = 4{,}5\text{ m}$, $E_{50}^{ref} = 8\text{ MPa}$, $E_{ur}^{ref} = 24\text{ MPa}$, $c = 15\text{ kPa}$);
  • Lớp 3: Bùn sét lẫn hữu cơ ($h = 7{,}0\text{ m}$, $E_{50}^{ref} = 4\text{ MPa}$, $E_{ur}^{ref} = 12\text{ MPa}$);
  • Lớp 4: Cát hạt mịn đến trung ($h = 12\text{ m}$, $E_{50}^{ref} = 35\text{ MPa}$, $\varphi = 32^\circ$);
  • Lớp 5: Cuội sỏi lẫn cát ($h = 15\text{ m}$, $E_{50}^{ref} = 80\text{ MPa}$, $E_{ur}^{ref} = 240\text{ MPa}$);
  • Lớp 6: Đá phiến sét/Đá gốc phong hóa ($E = 500\text{ MPa}$, $\nu = 0{,}25$).

Mô hình liên kết nửa cứng giữa các tấm vỏ hầm (LKNC) được thiết lập theo giải pháp liên kết phi tuyến Jassen: $$k_\theta = \frac{dM}{d\theta} = f(E_c, l_t, b, h, N)$$ trong đó $l_t$ là chiều cao vùng tiếp xúc chịu nén giữa các tấm phân tố, được mô hình hóa trong Plaxis 2D bằng phần tử bản trượt có độ cứng chống xoay biến thiên theo lực dọc ép $N$.

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

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

  • Phát hiện 1 (Khớp phổ tối ưu bảo toàn năng lượng): Chương trình PG01 tạo ra các giản đồ gia tốc BaDinh_01A, BaDinh_02A, BaDinh_03A đạt hệ số tương quan phổ $R^2 > 0{,}98$ so với phổ TCVN 9386-2012 chỉ sau 5 đến 8 vòng lặp. Cường độ Arias sau hiệu chỉnh duy trì ở mức thực tế ($I_A = 0{,}21 - 0{,}38\text{ m/s}$), triệt tiêu hoàn toàn hiện tượng trôi đường nền (baseline drift) và hiện tượng méo pha năng lượng thường gặp ở phương pháp biến đổi Fourier miền tần số.
  • Phát hiện 2 (Mô phỏng ngẫu nhiên theo đới đứt gãy đặc thù): Chương trình PG02 đã phát sinh thành công bộ 5 giản đồ gia tốc ngẫu nhiên (bd01-01a đến bd01-05a) tương thích với đới đứt gãy Sông Hồng - Sông Chảy ($M_w = 6{,}5$, $R_{rup} = 15\text{ km}$, $V_{s,30} = 250\text{ m/s}$). Gia tốc đỉnh đạt giá trị cực đại trung bình $\text{PGA} = 0{,}1098\text{ g}$, vận tốc đỉnh $\text{PGV} = 0{,}12 - 0{,}16\text{ m/s}$, khớp hoàn hảo với các kịch bản động đất nguy hiểm nhất tại nội đô Hà Nội.
  • Phát hiện 3 (Sự chênh lệch bản chất giữa phân tích Động và Tĩnh): So sánh kết quả tính toán giữa phương pháp Động lực học thời gian thực với phương pháp Tĩnh áp đặt chuyển vị biên (ISGD) và phương pháp HRM cho thấy:
    • Phương pháp ISGD đánh giá thấp mômen uốn cực đại $M_{max}$ trong vỏ hầm từ $18%$ đến $27%$ do không xét đến hiệu ứng khuếch đại dao động của các lớp bùn sét yếu phía trên.
    • Phương pháp HRM cho kết quả lực dọc ép $N_{max}$ lớn hơn thực tế $15%$, nhưng lại bỏ qua sự lệch pha giữa đỉnh sóng gia tốc và thời điểm đạt đỉnh của chuyển vị kết cấu.
  • Phát hiện 4 (Hiệu ứng giải tỏa ứng suất của liên kết nửa cứng): Khi xét đến liên kết nửa cứng Jassen giữa các phân tố vỏ hầm, mômen uốn lớn nhất xuất hiện trong vỏ hầm $M_{max}$ giảm từ $28{,}4%$ đến $34{,}6%$ so với mô hình vỏ hầm đúc liền khối liên tục (Monolithic continuous lining). Sự xoay mềm dẻo tại các chốt liên kết giúp giải phóng ứng suất uốn cục bộ, trong khi lực dọc $N$ chỉ biến thiên nhẹ ($< 6%$), đảm bảo độ kín khít chống thấm cho công trình ngầm.
  • Phát hiện 5 (Tương quan đa biến giữa đặc trưng dao động và nội lực): Phân tích hồi quy xác lập phương trình dự báo mômen uốn $M_{max}$ và lực cắt $Q_{max}$ theo các đại lượng động lực học của sóng gia tốc: $$\ln(M_{max}) = \beta_0 + \beta_1 \ln(\text{PGA}) + \beta_2 \ln(I_A) + \beta_3 \ln(t_{5-95}) + \beta_4 \ln(a_{RMS})$$ với mức ý nghĩa thống kê $p < 0{,}01$, trong đó năng lượng Arias $I_A$ và gia tốc hiệu dụng $a_{RMS}$ thể hiện mức độ chi phối nội lực uốn mạnh hơn đáng kể so với bản thân giá trị gia tốc đỉnh đơn lẻ PGA.

Implications đa chiều

  • Về mặt lý thuyết: Chứng minh tính ưu việt của toán tử Wavelet trong việc mô tả hiện tượng phi dừng của sóng động đất trong cơ học kết cấu ngầm, đặt nền tảng kết nối cơ học đứt gãy địa chấn với lý thuyết tương tác đất - kết cấu (SSI).
  • Về mặt phương pháp luận: Cung cấp bộ công cụ mã nguồn mở (PG01, PG02) độc lập, giải quyết sự phụ thuộc vào các phần mềm thương mại nước ngoài khi thực hiện thiết kế kháng chấn theo TCVN 9386-2012.
  • Về mặt thực tiễn kỹ thuật: Đưa ra chỉ dẫn tính toán trực tiếp cho Dự án Tuyến đường sắt đô thị số 3 Hà Nội (đoạn ngầm Kim Mã - Ga Hà Nội). Các kỹ sư kết cấu có thể tối ưu hóa cốt thép vỏ hầm nhờ lượng suy giảm mômen $30%$ từ hiệu ứng liên kết nửa cứng.
  • Về mặt chính sách và tiêu chuẩn: Cung cấp căn cứ thực nghiệm số để Bộ Xây dựng và Bộ Giao thông Vận tải bổ sung quy chuẩn tính toán động lực học cho hầm giao thông đô thị trong các phiên bản cập nhật của TCVN 9386.

Limitations và Future Research

Luận án thẳng thắn thừa nhận các giới hạn kỹ thuật:

  • Giới hạn mô hình không gian: Nghiên cứu tập trung vào bài toán biến dạng phẳng 2D, chưa xét đến hiệu ứng lan truyền sóng địa chấn xiên góc (Oblique wave propagation) dọc theo trục đường hầm 3D và hiện tượng uốn dọc thân hầm (Longitudinal bending).
  • Giới hạn cơ sở dữ liệu nền: Do thiếu mạng lưới đo đạc địa chấn nội địa dày đặc, các hệ số hồi quy trong PG02 vẫn kế thừa từ tập dữ liệu quốc tế NGA West của Yamamoto, cần được hiệu chỉnh thêm khi trạm quan trắc Việt Nam thu thập thêm dữ liệu thực tế.
  • Giới hạn vật liệu: Mô hình hóa chất lưu lỗ rỗng chưa xét đến quá trình phát sinh áp lực nước lỗ rỗng thặng dư gây hóa lỏng cát hạt mịn (Liquefaction) dưới rung động chu kỳ kéo dài.

Chương trình nghiên cứu tương lai (Future Agenda) gồm 4 hướng ưu tiên:

  1. Nâng cấp chương trình PG01/PG02 để phát sinh đồng thời 3 thành phần gia tốc không gian $(a_x(t), a_y(t), a_z(t))$ có xét đến tính tương quan tương hỗ giữa các phương.
  2. Xây dựng mô hình số 3D toàn diện khảo sát các vị trí đặc biệt như giếng chuyển làn, ga ngầm kết nối với thân hầm tròn chịu tác động của sóng địa chấn lan truyền không gian.
  3. Thực hiện thí nghiệm mô hình vật lý trên bàn rung ba chiều (3D Shake Table Test) hoặc máy ly tâm địa kỹ thuật (Geotechnical Centrifuge) để kiểm chứng cơ chế xoay của khớp nối nửa cứng.
  4. Tích hợp thuật toán học máy (Machine Learning / Deep Learning) để tự động hóa việc nhận dạng và trích xuất các hệ số Wavelet từ các nguồn dữ liệu địa chấn lớn trong khu vực Đông Nam Á.

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

  • Ảnh hưởng học thuật: Luận án mở ra hướng nghiên cứu kết hợp giữa xử lý tín hiệu tiên tiến và cơ học công trình ngầm, cung cấp tài liệu tham khảo chuẩn mực cho các nghiên cứu sinh chuyên ngành Kỹ thuật Xây dựng Công trình Đặc biệt (Mã số: 9 58 02 06).
  • Chuyển đổi công nghiệp xây dựng ngầm: Cung cấp công cụ tính toán chính xác, giúp các tổng thầu, đơn vị tư vấn thiết kế (như TEDI, VNCC, các liên danh quốc tế) tiết kiệm từ $8 - 12%$ chi phí vật liệu vỏ hầm nhờ việc mô hình hóa chính xác liên kết nửa cứng thay vì tăng chiều dày bê tông một cách cơ học.
  • Tác động an ninh - quốc phòng: Ứng dụng trực tiếp vào quy trình thiết kế, đánh giá độ an toàn của hệ thống công trình ngầm phòng thủ, công trình sở chỉ huy quốc phòng chịu tải trọng nổ và chấn động ngầm.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận cơ sở toán học tường minh của phương pháp biến đổi Wavelet (CWT, DWT, WPT, iWPT), thuật toán Hancock, Yamamoto và quy trình mô phỏng tương tác đất - hầm trong Plaxis 2D.
  • Kỹ sư Tư vấn Thiết kế Công trình Ngầm: Sở hữu phương pháp luận và công cụ PG01/PG02 để tạo dữ liệu đầu vào tính toán kháng chấn hợp chuẩn TCVN 9386-2012 cho các dự án Metro tại Hà Nội và TP. Hồ Chí Minh.
  • Ban Quản lý Đường sắt Đô thị (MRB Hà Nội, MAUR TP.HCM): Có cơ sở khoa học độc lập để thẩm tra, phản biện các giải pháp thiết kế vỏ hầm do tư vấn quốc tế đề xuất, đảm bảo an toàn tuyệt đối và tiết kiệm ngân sách nhà nước.

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ì?

Đóng góp độc đáo nhất là việc chuẩn hóa và mở rộng thuật toán hiệu chỉnh Wavelet miền thời gian của Hancock (2006) trên nền tảng hàm cơ sở Cosine cải biên, gắn kết trực tiếp với biểu thức toán học của phổ phản ứng đàn hồi quy định trong TCVN 9386-2012. Bằng cách thiết lập ma trận phản ứng đơn vị $[S]$ thông qua tích phân Duhamel cho từng chu kỳ riêng biệt, luận án đã biến đổi bài toán khớp phổ phi tuyến phức tạp thành hệ phương trình tuyến tính đại số ${\Delta S_a} = [S]{b}$, cho phép tìm nghiệm chính xác của lượng bù gia tốc $\delta a(t)$ mà không làm thay đổi các đặc trưng không dừng của sóng địa chấn gốc.

2. Sự đổi mới về mặt phương pháp luận so với các nghiên cứu trước đây?

So với các nghiên cứu trong nước của Nghiêm Mạnh Hiến (2012) vốn chỉ dùng phép biến đổi Fourier miền tần số làm xáo trộn cấu trúc năng lượng, hoặc công trình của Đinh Văn Thuật (2006) dựa trên giả thiết pha ngẫu nhiên dừng, luận án đã tạo bước đột phá phương pháp luận khi:

  • Tự phát triển hai phần mềm độc lập PG01 và PG02, vượt qua rào cản phụ thuộc vào phần mềm thương mại đóng gói SeismoMatch;
  • Kết hợp quy trình phát sinh sóng địa chấn với bài toán động lực học số tương tác phi tuyến đất - kết cấu trong Plaxis 2D, thay vì dừng lại ở các phân tích tĩnh đơn giản (ISGD, HRM) như các nghiên cứu của Lê Bảo Quốc (2013) hay Mai Đức Minh (2010).

3. Phát hiện bất ngờ nhất từ kết quả mô phỏng số là gì?

Phát hiện bất ngờ và có giá trị nhất là cơ chế tự điều tiết ứng suất của vỏ hầm lắp ghép: Khớp nối nửa cứng (Semi-rigid joints) đóng vai trò như các "khớp dẻo đàn hồi", làm suy giảm tới $34{,}6%$ mômen uốn cục bộ $M_{max}$ dưới tác động của động đất nhân tạo so với kết cấu đúc liền khối. Điều này bác bỏ quan điểm thiết kế truyền thống vốn cho rằng cần tăng độ cứng uốn toàn khối của vỏ hầm để kháng chấn, mở ra hướng tiếp cận thiết kế "mềm hóa có kiểm soát" nhằm tối ưu hóa khả năng hấp thụ năng lượng chấn động.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?

Hoàn toàn đầy đủ. Luận án cung cấp chi tiết:

  • Sơ đồ khối thuật toán của chương trình PG01 và PG02;
  • Toàn bộ công thức toán học giải ma trận $[S]$, tích phân Newmark, và biến đổi gói Wavelet ngược (iWPT);
  • Bảng tham số vật liệu chi tiết của 6 lớp địa tầng Hà Nội theo mô hình Hardening Soil ($E_{50}^{ref}, E_{oed}^{ref}, E_{ur}^{ref}, c, \varphi, \psi, R_{inter}$);
  • Tọa độ địa lý, thông số nguồn chấn đới Sông Hồng - Sông Chảy ($M_w = 6{,}5; R_{rup} = 15\text{ km}; V_{s,30} = 250\text{ m/s}$);
  • Các bước khai báo điều kiện biên nhớt (viscous boundary) và bước thời gian tích phân động lực học $\Delta t = 0{,}01\text{ s}$ trong PLAXIS 2D.

5. Định hướng nghiên cứu 10 năm tới được phác thảo ra sao?

Chương trình nghiên cứu 10 năm xác lập 3 trụ cột phát triển:

  • Giai đoạn 2022-2025: Hoàn thiện mô hình bài toán 3D phân tích hiệu ứng không gian của sóng địa chấn truyền không đồng bộ lên toàn tuyến đường hầm và nhà ga ngầm;
  • Giai đoạn 2026-2028: Thực hiện thí nghiệm bàn rung kích thước lớn kiểm chứng ứng xử của liên kết nửa cứng dưới tác động của giản đồ gia tốc phát sinh từ PG01/PG02;
  • Giai đoạn 2029-2032: Ứng dụng Trí tuệ nhân tạo (AI) và Bản sao số (Digital Twin) để xây dựng hệ thống quan trắc, cảnh báo sớm và tự động đánh giá độ an toàn kết cấu ngầm đô thị theo thời gian thực khi xảy ra động đất.

Kết luận

  1. Luận án đã làm chủ toàn diện cơ sở lý thuyết và xây dựng thành công chương trình PG01 trên nền tảng Matlab, cho phép phát sinh các giản đồ gia tốc nhân tạo trên nền đá gốc bằng phương pháp hiệu chỉnh Wavelet miền thời gian theo thuật toán Hancock, đáp ứng hoàn hảo điều kiện khớp phổ phản ứng đàn hồi theo tiêu chuẩn TCVN 9386-2012 với sai số $< 5%$.
  2. Nghiên cứu đã cải biên thành công thuật toán Yamamoto, xây dựng chương trình PG02 mô phỏng ngẫu nhiên phi dừng các giản đồ gia tốc nền thông qua hệ phương trình hồi quy 13 tham số gói Wavelet kết hợp biến đổi ngược iWPT, phản ánh chân thực đặc trưng của đới động đất Sông Hồng - Sông Chảy đối với khu vực trung tâm Hà Nội.
  3. Thiết lập mô hình phần tử hữu hạn động lực học 2D trên phần mềm PLAXIS với mô hình đất nền Hardening Soil (HS), phản ánh chính xác tương tác phức tạp giữa địa tầng trầm tích mềm yếu Hà Nội và kết cấu hầm ngầm đô thị dưới tác động của các kích thích động lực học thời gian thực.
  4. Khám phá và định lượng hóa chính xác hiệu ứng giảm mômen uốn cực đại ($28{,}4% - 34{,}6%$) của liên kết nửa cứng Jassen giữa các tấm vỏ hầm lắp ghép, cung cấp luận cứ khoa học then chốt để tối ưu hóa giải pháp kết cấu vỏ hầm Metro Tuyến số 3 Hà Nội.
  5. Xác lập hệ phương trình hồi quy đa biến biểu diễn mối quan hệ phụ thuộc của nội lực cực đại ($M_{max}, Q_{max}$) vào các chỉ số năng lượng địa chấn ($I_A, a_{RMS}, t_{5-95}$), chứng minh rằng phân tích kháng chấn cho công trình ngầm bắt buộc phải dựa trên toàn bộ diễn biến lịch sử thời gian của sóng gia tốc thay vì chỉ dựa vào gia tốc đỉnh PGA.
  6. Mở ra hướng nghiên cứu liên ngành chuyên sâu giữa Xử lý Tín hiệu Số Địa chấn và Kỹ thuật Công trình Ngầm Đô thị, để lại di sản khoa học ứng dụng trực tiếp cho sự nghiệp phát triển hạ tầng giao thông ngầm và công trình phòng thủ quốc phòng của Việt Nam.