Tổng quan về luận án

Sự phát triển vượt bậc của hạ tầng giao thông ngầm và khai thác khoáng sản đặt ra thách thức kỹ thuật phức tạp khi xuất hiện ngày càng nhiều dự án thi công hầm mới nằm song song hoặc mở rộng sát các công trình ngầm hiện hữu. Trong điều kiện địa chất đô thị và đồi núi phức tạp, phương pháp khoan nổ mìn truyền thống giải phóng nguồn năng lượng động học cực lớn, truyền sóng ứng suất trong khối đá và tác động tiêu cực đến độ ổn định kết cấu vỏ chống công trình lân cận. Luận án tiến sĩ kỹ thuật với đề tài "Nghiên cứu ảnh hưởng của chấn động nổ mìn khi thi công đường hầm đến kết cấu công trình ngầm lân cận" do nghiên cứu sinh Đặng Văn Kiên thực hiện dưới sự hướng dẫn khoa học của GS. Võ Trọng Hùng (Trường Đại học Mỏ - Địa chất Hà Nội, 2018, Chuyên ngành: Kỹ thuật xây dựng Công trình ngầm, Mã số: 9580204) đã giải quyết trọn vẹn bài toán động lực học công trình ngầm mang tính thời sự và cấp bách này.

Khoảng trống nghiên cứu (research gap) lớn nhất được xác định là sự thiếu hụt các mô hình dự báo động học 3 chiều (3D) tích hợp chính xác các thông số cơ học động của khối đá, trong khi phần lớn các nghiên cứu trước đây tại Việt Nam chỉ dừng lại ở các bài toán tĩnh, áp dụng phương pháp tải trọng tĩnh tương đương, hoặc mô phỏng các vụ nổ đơn lẻ trên mặt đất phục vụ công trình quân sự. Các mô hình hiện hữu hầu như bỏ qua sự suy giảm năng lượng sóng nổ theo thời gian, ảnh hưởng của hệ số giảm chấn ($\xi$), đặc tính liên kết tiếp xúc giữa vỏ hầm và khối đá ("lớp lát hàn" hay "lớp lát lồng"), cũng như mối tương quan định lượng giữa chất lượng khối đá theo hệ số đánh giá khối đá RMR (Rock Mass Rating) với các hệ số truyền sóng chấn động.

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

  • RQ1: Tải trọng động sinh ra do quá trình nổ mìn tại gương hầm biến đổi như thế nào trong môi trường khối đá rắn chắc và tác động ra sao đến ứng xử cơ học của vỏ chống bê tông liền khối của hầm cũ lân cận?
    • H1: Sóng chấn động nổ mìn truyền trong môi trường khối đá đàn hồi tuyến tính gây ra trường ứng suất kéo biến đổi có biên độ phụ thuộc nghịch đảo vào khoảng cách và tỷ lệ thuận với căn bậc hai của khối lượng thuốc nổ nạp lớn nhất cho một đợt nổ vi sai ($Q$).
  • RQ2: Các thông số động học của khối đá (mô đun đàn hồi động $E_d$, tốc độ biến dạng $\dot{\varepsilon}$, độ bền động) thu được từ thí nghiệm động chuyên dụng có sai khác như thế nào so với các thông số tĩnh truyền thống?
    • H2: Dưới tốc độ biến dạng cao ($\dot{\varepsilon} \ge 10^1 \text{ s}^{-1}$), mô đun đàn hồi động và độ bền của đá granit tăng đáng kể so với điều kiện tải tĩnh, đóng vai trò quyết định đến tính chính xác của mô hình mô phỏng số.
  • RQ3: Mối quan hệ định lượng giữa chỉ số chất lượng khối đá RMR với các thông số suy giảm sóng nổ ($K, \alpha$) trong công thức thực nghiệm suy biến diễn ra theo quy luật nào?
    • H3: Giá trị hệ số suy giảm sóng $\alpha$ và hệ số thực nghiệm $K$ tương quan phi tuyến chặt chẽ với RMR và vị trí tương đối của gương nổ mìn dọc theo trục đường hầm ($H$).
  • RQ4: Ứng dụng điều kiện biên không phản xạ (Non-Reflecting Boundary Condition - NRBC) và phần tử vô hạn (Infinite Elements - IE) trong mô hình phần tử hữu hạn 3D trên nền tảng Abaqus có thể triệt tiêu sóng phản xạ và mô tả chính xác cơ chế phá hủy của vỏ chống theo chỉ số phá hủy nổ mìn BDI (Blast Damage Index) hay không?
    • H4: Mô hình 3D tích hợp NRBC với hệ số giảm chấn tối ưu $\xi = 5{,}0%$ cho phép tái hiện chân thực vận tốc dao động phần tử đỉnh (PPV) thực tế, xác định chính xác vị trí tập trung ứng suất kéo nguy hiểm nhất trên vỏ hầm cũ.

Khung lý thuyết của nghiên cứu được xây dựng trên nền tảng cơ học môi trường liên tục, lý thuyết truyền sóng đàn hồi của Navier-Stokes-Cauchy, mô hình tải trọng nổ áp lực xung suy giảm theo thời gian, kết hợp hệ thống chỉ số phá hủy BDI. Đóng góp đột phá của luận án nằm ở việc định lượng hóa chính xác giới hạn an toàn PPV, xác định lượng thuốc nổ nạp lớn nhất cho phép ($Q_{max}$) nhằm bảo vệ an toàn tuyệt đối cho vỏ hầm cũ. Phạm vi nghiên cứu bao quát dữ liệu quan trắc thực tế từ 54 vụ nổ tại dự án hầm đường bộ Croix-Rousse (Lyon, Cộng hòa Pháp) với chiều dài 1757,5 m, diện tích mặt cắt ngang 84,10 m² trong khối đá granit có độ bền nén đơn trục $UCS > 120\text{ MPa}$, kết hợp thực nghiệm động Split Hopkinson Pressure Bar (SHPB) tại phòng thí nghiệm chuyên sâu của Viện Khoa học Ứng dụng Quốc gia Lyon (INSA Lyon, Pháp) và các công trình hầm trọng điểm tại Việt Nam như hầm Cổ Mã (khoảng cách 47 m tới hầm đường sắt số 24) và hầm Hải Vân (khoảng cách 30 m giữa hầm chính và hầm lánh nạn).

Literature Review và Positioning

Lịch sử kỹ thuật công trình ngầm thế giới đã ghi nhận nhiều thảm kịch do sóng nổ chấn động gây sụp đổ công trình lân cận. Điển hình là thảm họa tại đường hầm lấy nước Port Huron (bang Michigan, Hoa Kỳ, 1971) thi công ở độ sâu 220 feet dưới đáy hồ Huron, khiến 21 công nhân thiệt mạng và 9 người bị thương do sóng nổ phá hủy hoàn toàn kết cấu vách hầm; hay sự cố sập đổ 123,0 m đường hầm đường sắt Jiuguaidao (Trung Quốc, 1980) với hơn 3000,0 m³ đất đá sạt lở do hoạt động nổ mìn đào mái dốc gần kề, gây thiệt hại hàng tỷ USD.

Trong dòng chảy học thuật quốc tế, các giải pháp giải tích đã được Pao Y. và Wersall C. (2008) tiếp cận thông qua biến đổi tích phân Laplace và Fourier cho bài toán đường hầm mặt cắt tròn chịu sóng nổ hình trụ. Tuy nhiên, các giải pháp giải tích thuần túy bộc lộ nhược điểm cố hữu khi không thể giải quyết hình học phức tạp (mặt cắt vòm móng ngựa, vòm tường thẳng) hay môi trường khối đá dị hướng. Ở nhánh mô hình vật lý thu nhỏ, Khosrow B. (1998) dùng mô hình hộp thép 1-G ($150 \times 100 \times 500\text{ mm}$) để đo áp lực nổ, nhưng gặp giới hạn lớn về hiệu ứng kích thước (scale effect).

Về hướng quan trắc thực nghiệm và mô phỏng số, Ansell A. (2004) cùng Ahmed A. & Ansell A. (2014) đã phát triển mô hình liên kết dầm - lò xo để đánh giá vỏ bê tông phun chịu tải trọng nổ mìn tập trung. Jong-Ho Shin et al. (2011) sử dụng mô hình 2D thiết lập ranh giới vùng bảo vệ công trình ngầm trong đá yếu dựa trên PPV và chuyển vị. Li et al. (2013) nghiên cứu tác động của nổ mìn đào hầm mới tới tuyến hầm đường sắt đang hoạt động tại Tân Cương (Xinjiang, Trung Quốc), qua đó chỉ ra rằng mức độ nứt vỡ khối đá tăng tuyến tính theo sự gia tăng của PPV và xác lập ngưỡng giới hạn vận tốc dao động $[PPV] = 0{,}22\text{ m/s}$ ($220\text{ mm/s}$) cho kết cấu chống giữ. Liang Q. (2012) và Xia et al. (2013) tiếp tục khai thác phần mềm LS-DYNA để khảo sát trường ứng suất động và sóng phản xạ.

Tại Việt Nam, các công trình của Hoàng Thị Hồng (nghiên cứu trường ứng suất quanh buồng nổ hình cầu/trụ theo mô hình lưu biến Kelvin và Poynting-Thomson), Nguyễn Xuân Mãn (xác định bán kính phá hủy của lượng nổ đơn độc trong môi trường đàn hồi), Lê Đình Tân (khảo sát động lực học kết cấu hầm dạng vòm tường thẳng chịu nổ bề mặt trên phần mềm SAP2000 bằng phương pháp tích phân trực tiếp Newmark kết hợp mô hình hộp san hô tại hầm Truông Khấp), hay Nguyễn Hữu Thế (mô phỏng nổ mìn trên đảo Sơn Ca - Trường Sa bằng Plaxis và Autodyn) chủ yếu tập trung vào kết cấu phòng thủ quân sự chịu tác động bom đạn nổ bề mặt.

Luận án của Đặng Văn Kiên định vị vững chắc trong khoảng trống học thuật: Là công trình tiên phong tại Việt Nam tích hợp dữ liệu chấn động thực địa quy mô lớn (dự án Croix-Rousse, Pháp), kết hợp thực nghiệm cơ học động SHPB xác lập chính xác thông số động ($E_d, \dot{\varepsilon}$), từ đó xây dựng mô hình phần tử hữu hạn 3D Abaqus với biên vô hạn không phản xạ NRBC, làm sáng tỏ cơ chế tương tác dọc trục hầm ($H$) và thiết lập hàm tương quan giữa RMR với hệ số suy giảm truyền sóng Chapot.

Đó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 căn bản lý thuyết đàn hồi động học (Elastodynamics) và lý thuyết truyền sóng địa chấn nổ mìn thông qua việc chứng minh rằng: Trong môi trường khối đá granit đồng nhất, sóng ứng suất nén sơ cấp (P-wave) và sóng cắt thứ cấp (S-wave) khi truyền đến bề mặt tiếp xúc giữa khối đá và vỏ bê tông liền khối ("lớp lát hàn") sẽ chuyển hóa thành trường ứng suất động phức hợp. Tại mặt biên trong tự do của vỏ hầm cũ, sóng nén phản xạ thành sóng kéo; do cường độ chịu kéo của bê tông chỉ bằng khoảng $8 - 10%$ cường độ chịu nén, thành phần ứng suất kéo chính là nguyên nhân cốt lõi gây ra nứt nẻ và bong tróc kết cấu.

Khung phân tích lý thuyết thiết lập mối liên kết đa tầng:

  1. Mô hình truyền sóng khối vi phân Navier-Cauchy: $$(\lambda + G)\nabla(\nabla \cdot \mathbf{u}) + G\nabla^2 \mathbf{u} = \rho \frac{\partial^2 \mathbf{u}}{\partial t^2}$$ Trong đó $\lambda, G$ là các hằng số Lamé, $\mathbf{u}$ là vector chuyển vị, $\rho$ là khối lượng thể tích của khối đá.
  2. Mô hình áp lực xung nổ trên biên hầm theo hàm thời gian: Áp lực nổ $P(t)$ tăng đột ngột đạt đỉnh cực đại $P_{max}$ trong thời gian cực ngắn ($t_r \approx 1 - 2\text{ ms}$) và suy giảm theo quy luật hàm mũ: $$P(t) = P_{max} \cdot e^{-\beta t}$$ mô tả chính xác quá trình giải phóng khí nén TNT và sóng kích nổ vi sai nhiều đợt.
  3. Chuyển đổi hệ hình (Paradigm shift): Chuyển dịch từ phương pháp tính toán tải trọng tĩnh tương đương sang phân tích lịch sử thời gian động lực học phi tuyến (Nonlinear Dynamic Time-History Analysis) với các thông số động thực nghiệm đo bằng SHPB.

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

Khung phân tích của luận án tích hợp 3 trụ cột lý thuyết chuyên sâu:

  • Lý thuyết đàn hồi động học môi trường liên tục (Continuum Dynamic Elasticity Theory): Xác định quy luật lan truyền vận tốc phần tử và ứng suất trong trường đàn hồi đẳng hướng.
  • Mô hình suy giảm sóng thực nghiệm mở rộng (Extended Empirical Blast Attenuation Model): Nâng cấp công thức cổ điển của Chapot, USBM (United States Bureau of Mines), Langefors và Ambraseys-Hendron bằng cách đưa hệ số chất lượng khối đá RMR vào hàm phụ thuộc của $K$ và $\alpha$: $$\text{PPV} = K \cdot \left(\frac{D}{\sqrt{Q}}\right)^{-\alpha} = f(\text{RMR}, H) \cdot \left(\frac{D}{\sqrt{Q}}\right)^{-\alpha(\text{RMR}, H)}$$
  • Lý thuyết cơ học phá hủy và chỉ số tổn thương BDI (Blast Damage Index Theory): Đánh giá định lượng trạng thái tới hạn của vỏ hầm thông qua tương quan giữa ứng suất kéo động sinh ra ($\sigma_t^d$) và độ bền kéo động của bê tông ($R_t^d$): $$\text{BDI} = \frac{\sigma_t^d}{R_t^d}$$
    • $\text{BDI} < 0{,}5$: Không gây phá hủy;
    • $0{,}5 \le \text{BDI} \le 1{,}0$: Xuất hiện vi khe nứt bề mặt;
    • $\text{BDI} > 1{,}0$: Phá hủy giòn và sụp đổ cục bộ kết cấu vỏ bê tông.

Điều kiện biên lý thuyết (Boundary Conditions) được xác lập nghiêm ngặt: Môi trường làm việc là lớp tiếp xúc cứng hoàn toàn ("lớp lát hàn" - bonded contact lining), thỏa mãn điều kiện mô đun biến dạng của khối đá $E_0 \ge E_{bê tông}$ và không xuất hiện hiện tượng trượt tách rời trên bề mặt tiếp xúc trong suốt quá trình truyền sóng.

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ chặt chẽ triết lý thực chứng (Positivism) kết hợp chủ nghĩa duy thực phê phán (Critical Realism). Thiết kế nghiên cứu đa phương pháp (Mixed-Method Triangulation) kết hợp 3 nhánh độc lập:

  1. Nhánh đo đạc thực nghiệm hiện trường (Field Empirical Monitoring): Quan sát 54 đợt nổ thực tế tại dự án hầm giao thông đô thị Croix-Rousse (Pháp).
  2. Nhánh thí nghiệm cơ học động trong phòng (Laboratory Dynamic Testing): Vận hành hệ thống thanh ép động Split Hopkinson Pressure Bar (SHPB) tại INSA Lyon.
  3. Nhánh mô phỏng số nâng cao (Advanced Numerical Simulation): Xây dựng mô hình 2D và 3D phần tử hữu hạn trên phần mềm Abaqus/Explicit CAE tích hợp biên vô hạn NRBC.

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu được chuẩn hóa cao độ:

  • Quan trắc hiện trường tại hầm Croix-Rousse:
    • Tuyến hầm đào mới có tiết diện hình móng ngựa (bán kính $R = 5{,}55\text{ m}$), hầm cũ song song dạng vòm tường thẳng ($R = 8{,}05\text{ m}$, chiều cao tường $1{,}0\text{ m}$), khoảng cách tim hai hầm là $42{,}6\text{ m}$. Khối đá vách là đá granit rắn chắc ($UCS > 120\text{ MPa}$, khối lượng thể tích $\rho = 2650\text{ kg/m}^3$).
    • Sử dụng hệ thống cảm biến đo địa chấn 3 thành phần (Geophone triaxial sensors: Cảm biến P, T, A và C100) lắp đặt dọc theo vách vỏ bê tông hầm cũ tại các vị trí lý trình PM100, PM220, PM340, PM1400.
    • Ghi nhận đồng thời 3 thành phần vector vận tốc dao động: phương thẳng đứng (Vertical - V), phương ngang vuông góc trục hầm (Horizontal - H), và phương dọc trục hầm (Longitudinal - L) ở dải tần số thấp ($f = 130\text{ Hz}$) cùng chuyển vị vi mô.
  • Thí nghiệm thanh nén động SHPB tại INSA Lyon:
    • Gia công các mẫu đá granit hình trụ tròn đạt chuẩn ISRM với kích thước đường kính $\Phi = 45{,}0\text{ mm}$, chiều cao $L = 100{,}0\text{ mm}$ (tỷ lệ $L/D = 2{,}22$).
    • Mẫu được kẹp giữa hai thanh thép hợp kim có độ đàn hồi cao: thanh tới (incident bar) và thanh truyền (transmitter bar).
    • Tải trọng động sinh ra bằng thanh đánh (striker bar) được phóng bằng khí nén với hai cấp áp lực: $0{,}30\text{ MPa}$ và $0{,}35\text{ MPa}$.
    • Ghi nhận biến dạng vi mô trên thanh tới ($\varepsilon_i, \varepsilon_r$) và thanh truyền ($\varepsilon_t$) thông qua hệ thống cảm biến điện trở dán (strain gauges) nối với mạch cầu Wheatstone khuếch đại tín hiệu điện áp cao tần.
  • Kiểm định độ tin cậy và giá trị:
    • Dữ liệu sóng động được lọc nhiễu dải thông (bandpass filter); hệ số biến thiên dữ liệu $CV < 5%$.
    • Hiệu chuẩn sai số giữa mô hình số Abaqus và số liệu đo đạc thực tế tại điểm kiểm chứng PM100 đạt độ khớp trên $92%$.

Data và phân tích

  • Đặc trưng cơ học động của mẫu đá granit (Mẫu điển hình N°38):
    • Khối lượng thể tích: $\rho = 2645\text{ kg/m}^3$; vận tốc truyền sóng dọc $C_p = 5120\text{ m/s}$.
    • Dưới áp lực phóng $0{,}30\text{ MPa}$: tốc độ biến dạng mẫu đạt $\dot{\varepsilon} \approx 65\text{ s}^{-1}$, ứng suất nén động đỉnh đạt $148{,}5\text{ MPa}$ (tăng $23{,}7%$ so với độ bền nén tĩnh $120\text{ MPa}$).
    • Dưới áp lực phóng $0{,}35\text{ MPa}$: tốc độ biến dạng đạt $\dot{\varepsilon} \approx 112\text{ s}^{-1}$, ứng suất nén động đỉnh đạt $172{,}3\text{ MPa}$.
    • Mô đun đàn hồi động xác lập từ thực nghiệm: $E_d = 62{,}5\text{ GPa}$ (so với mô đun tĩnh $E_s = 48{,}0\text{ GPa}$).
  • Thiết lập mô hình số 2D và 3D trong Abaqus:
    • Kích thước miền tính toán 3D: Chiều rộng $120\text{ m}$, chiều cao $80\text{ m}$, chiều dài dọc trục hầm $100\text{ m}$.
    • Rời rạc hóa bằng phần tử khối 8 nút C3D8R kết hợp lưới phần tử vô hạn CIN3D8 tại 4 mặt biên bao quanh và đáy nhằm triệt tiêu hoàn toàn sóng phản xạ (NRBC).
    • Kích thước lưới phần tử tại vùng tập trung ứng suất quanh vách hầm thỏa mãn điều kiện truyền sóng chính xác của Kuhlemeyer & Lysmer: $\Delta l \le \lambda_{min} / (8 \div 10) = 0{,}8\text{ m}$.
    • Khảo sát độ nhạy hệ số giảm chấn Rayleigh ($\xi$): chạy mô phỏng biến thiên $\xi$ từ $1{,}0%$ đến $10{,}0%$.

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

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

Luận án đã đạt được 5 phát hiện khoa học mang tính đột phá với đầy đủ bằng chứng thực nghiệm và mô phỏng số:

  1. Xác lập quy luật ảnh hưởng nghịch đảo của hệ số giảm chấn ($\xi$) đến PPV: Kết quả mô phỏng số 3D chỉ ra mối quan hệ tỷ lệ nghịch phi tuyến giữa $\xi$ và giá trị cực đại của PPV trong vỏ bê tông hầm cũ. Khi $\xi$ tăng từ $1{,}0%$ lên $5{,}0%$, biên độ PPV giảm mạnh $41{,}8%$; từ $5{,}0%$ đến $10{,}0%$, mức độ suy giảm chậm dần. Tại giá trị $\xi = 5{,}0%$, kết quả mô phỏng số khớp chính xác nhất với giá trị PPV ghi nhận từ cảm biến thực địa Geophone tại điểm A và B ($PPV_{thực nghiệm} = 18{,}6\text{ mm/s}$ so với $PPV_{mô hình} = 18{,}2\text{ mm/s}$, sai số dưới $2{,}2%$). Đây là bằng chứng thực nghiệm then chốt chứng minh $\xi = 5{,}0%$ là thông số giảm chấn chuẩn tắc cho môi trường đá granit nguyên khối.

  2. Khám phá phân bố không gian 3D của trường chấn động dọc trục hầm ($H$): Khác biệt hoàn toàn với giả định truyền thống trên mô hình 2D (cho rằng chấn động chỉ xảy ra đồng phẳng), mô hình 3D và số liệu thực địa chứng minh: Giá trị PPV cực đại luôn xuất hiện tại mặt cắt ngang của vỏ hầm cũ trùng đúng với mặt phẳng gương nổ mìn của hầm mới ($H = 0\text{ m}$). Khi khoảng cách tương đối dọc trục hầm $|H|$ tăng lên, PPV suy giảm theo hàm mũ bất đối xứng. Trong đó, miền vỏ hầm nằm phía sau gương nổ ($H < 0$, vùng đã đào) có tốc độ suy giảm sóng chấn động nhanh hơn $18{,}5%$ so với miền khối đá nguyên khối phía trước gương nổ ($H > 0$, vùng chưa đào), do hiệu ứng che chắn năng lượng (shadow zone) của khoảng rỗng đường hầm đã hình thành.

  3. Thiết lập hệ phương trình thực nghiệm liên kết trực tiếp chỉ số RMR với các hệ số Chapot: Dựa trên phân tích hồi quy phi tuyến từ 54 vụ nổ tại dự án Croix-Rousse, tác giả đã xây dựng thành công hàm thực nghiệm dự báo hệ số $K$ và số mũ suy giảm $\alpha$ theo chỉ số RMR ($55 \le \text{RMR} \le 85$) cho từng phân vùng không gian:

    • Khi vị trí quan sát nằm phía trước gương nổ ($H > 0$): $$\ln(K) = 0{,}0382 \cdot \text{RMR} + 4{,}125 \quad (R^2 = 0{,}912)$$ $$\alpha = -0{,}0145 \cdot \text{RMR} + 2{,}418 \quad (R^2 = 0{,}886)$$
    • Khi vị trí quan sát nằm trong phạm vi ảnh hưởng mạnh ($0 < H < 45\text{ m}$): $$K = 128{,}4 \cdot e^{0{,}028 \cdot \text{RMR}}$$ $$\alpha = 1{,}85 - 0{,}0085 \cdot \text{RMR}$$ Phát hiện này giải thích hiện tượng phản trực giác: Khối đá có chất lượng càng cao ($\text{RMR} > 75$, đá rắn chắc, ít nứt nẻ) thì hệ số truyền năng lượng sóng nổ $K$ càng lớn, độ suy giảm $\alpha$ càng nhỏ, dẫn đến mức độ lan truyền chấn động xa hơn và gây rủi ro chấn động cao hơn cho vỏ hầm cũ so với môi trường đá nứt nẻ yếu.

  4. Định lượng hóa cơ chế phá hủy vỏ bê tông theo chỉ số BDI: Phân tích trường ứng suất động cho thấy vị trí nguy hiểm nhất trên vỏ hầm cũ nằm tại phần vách hông (hông tường và vai vòm) quay trực diện về phía tâm nổ của hầm mới. Tại đây, ứng suất kéo chính $\sigma_1$ đạt cực đại $3{,}85\text{ MPa}$, vượt quá giới hạn kéo động của bê tông M300 ($R_t^d = 2{,}85\text{ MPa}$), tương ứng chỉ số $\text{BDI} = 1{,}35 > 1{,}0$, hình thành dải phá hủy kéo giòn dạng vi nứt dọc bề mặt trong của vỏ hầm.

  5. Xây dựng biểu đồ tương quan vận tốc dao động cho phép [PPV] theo các tiêu chuẩn quốc tế: So sánh thực nghiệm chỉ ra rằng tiêu chuẩn DIN 4150-3 (Đức) và SN 640 312:1978 (Thụy Sĩ) áp đặt ngưỡng an toàn khắt khe ở dải tần số thấp ($[PPV] \le 15 - 20\text{ mm/s}$), trong khi tiêu chuẩn GB 6722-2003 (Trung Quốc) và tiêu chuẩn Pháp cho phép ngưỡng dao động động lực học trong kết cấu ngầm kiên cố lên tới $50 - 100\text{ mm/s}$ tùy thuộc tần số dao động chi phối ($FCPV$).

Implications đa chiều

  • Về mặt lý thuyết: Mở rộng lý thuyết tương tác sóng động - công trình ngầm ngầm hóa (Soil-Structure Dynamic Interaction), cung cấp khung giải tích - số hóa hoàn chỉnh cho bài toán hầm đôi song song thi công ngầm.
  • Về mặt phương pháp luận: Thiết lập quy trình mẫu (standard protocol) kết hợp thí nghiệm động SHPB, quan trắc địa chấn Geophone 3 trục và mô phỏng số Abaqus/Explicit với biên vô hạn NRBC, có thể chuyển giao áp dụng cho mọi loại hình công trình ngầm dân dụng và quốc phòng.
  • Về mặt thực tiễn kỹ thuật: Cung cấp công thức tính toán lượng thuốc nổ nạp lớn nhất cho một đợt nổ vi sai ($Q_{max}$) phục vụ trực tiếp công tác thiết kế hộ chiếu khoan nổ mìn: $$Q_{max} = D^2 \cdot \left(\frac{[\text{PPV}]}{K}\right)^{2/\alpha}$$ giúp các nhà thầu kiểm soát rung chấn tuyệt đối, không gây nứt vỡ công trình ngầm kề bên.
  • Về mặt chính sách và quy chuẩn: Cung cấp căn cứ khoa học chuẩn xác để Bộ Xây dựng và Bộ Giao thông Vận tải Việt Nam cập nhật quy phạm kỹ thuật an toàn nổ mìn cho công trình ngầm dân dụng, thay thế các tiêu chuẩn nổ mìn khai thác mỏ lộ thiên đã lỗi thời.

Limitations và Future Research

Luận án thẳng thắn chỉ rõ 3 giới hạn nghiên cứu chính:

  1. Giả định môi trường đàn hồi tuyến tính đồng nhất: Khối đá trong mô hình số được giả thiết là môi trường đàn hồi đồng nhất và đẳng hướng, chưa mô phỏng chi tiết mạng lưới khe nứt ngẫu nhiên (Discrete Fracture Network - DFN) và các đới đứt gãy kiến tạo phức tạp.
  2. Điều kiện tiếp xúc liên kết cứng: Mô hình chỉ xét trường hợp "lớp lát hàn" (bê tông bám dính tuyệt đối vào đá), chưa đánh giá trường hợp "lớp lát lồng" có khe hở, lớp chống thấm đàn hồi hoặc hiện tượng trượt tách rời cục bộ.
  3. Phạm vi thí nghiệm vật liệu: Thực nghiệm SHPB tập trung chủ yếu trên mẫu đá granit rắn chắc ($UCS > 120\text{ MPa}$), chưa bao quát đầy đủ các thang đo trên các loại đá biến chất mềm yếu (đá phiến sét, đá vôi nứt nẻ mạnh).

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

  • Ứng dụng phương pháp phần tử rời rạc (DEM/3DEC) hoặc phương pháp liên tục - rời rạc kết hợp (FEM-DEM) để mô phỏng sự mở rộng vết nứt vi mô trong khối đá bất đẳng hướng.
  • Khảo sát ảnh hưởng tích lũy tổn thương mỏi (fatigue damage accumulation) của vỏ hầm dưới tác động của hàng nghìn chu kỳ nổ mìn liên tiếp trong suốt quá trình đào đường hầm.
  • Nghiên cứu cơ chế giảm chấn chủ động bằng các giải pháp kỹ thuật như tạo hào giảm chấn (presplitting/smooth blasting slot) nằm giữa hai đường hầm.
  • Tích hợp trí tuệ nhân tạo (Machine Learning/Deep Neural Networks) huấn luyện trên cơ sở dữ liệu lớn quan trắc sóng chấn động để dự báo PPV thời gian thực.

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

Luận án tạo ra tác động khoa học và thực tiễn sâu rộng:

  • Ảnh hưởng học thuật: Đặt nền móng cho các nghiên cứu chuyên sâu về động lực học công trình ngầm tại Việt Nam; cung cấp bộ dữ liệu thực nghiệm mẫu và phương pháp luận SHPB chuẩn xác phục vụ các chương trình đào tạo tiến sĩ, thạc sĩ ngành Xây dựng Công trình ngầm và Khai thác Mỏ.
  • Chuyển đổi công nghệ ngành xây dựng ngầm: Ứng dụng trực tiếp vào tối ưu hóa hộ chiếu nổ mìn tại các dự án giao thông trọng điểm quốc gia như: Dự án mở rộng hầm lánh nạn đường bộ qua Đèo Hải Vân (khoảng cách 30 m tới hầm chính), dự án hầm Cổ Mã (Đèo Cả), các dự án thủy điện tích năng và các tuyến metro ngầm đô thị tại Hà Nội và TP. Hồ Chí Minh.
  • Lợi ích kinh tế - xã hội định lượng: Ngăn ngừa hoàn toàn sự cố nứt vỡ kết cấu vỏ hầm lân cận, tiết kiệm hàng chục tỷ đồng chi phí gia cố sửa chữa, giảm thiểu rủi ro gián đoạn giao thông huyết mạch quốc gia, đồng thời giảm lượng tiêu hao thuốc nổ không cần thiết, góp phần bảo vệ môi trường địa chất đô thị.

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

  • Nghiên cứu sinh và nhà khoa học chuyên ngành Công trình ngầm/Cơ học đá: Tiếp cận quy trình nghiên cứu tích hợp hoàn chỉnh từ lý thuyết, thực nghiệm SHPB phòng thí nghiệm, đo đạc địa chấn thực địa đến mô phỏng số 3D cao cấp.
  • Kỹ sư tư vấn thiết kế và giám sát nổ mìn: Sở hữu công cụ và công thức kinh nghiệm chuẩn xác để tính toán lượng thuốc nổ cho phép $Q_{max}$, thiết kế hộ chiếu nổ mìn vi sai tối ưu an toàn.
  • Các Ban quản lý dự án và Nhà thầu thi công hạ tầng giao thông: Có cơ sở khoa học tin cậy để thẩm định biện pháp thi công, lắp đặt hệ thống cảnh báo chấn động Geophone tự động, đảm bảo an toàn tuyệt đối cho các công trình hiện hữu kề bên.
  • Cơ quan quản lý nhà nước (Bộ Giao thông Vận tải, Bộ Xây dựng, Bộ Công Thương): Sử dụng các phát hiện khoa học làm căn cứ hoàn thiện hệ thống tiêu chuẩn, quy chuẩn kỹ thuật quốc gia về đánh giá tác động chấn động nổ mìn trong công trình dân dụng.

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à việc mở rộng lý thuyết suy giảm sóng nổ Chapot thông qua việc tích hợp chỉ số chất lượng khối đá RMR và thông số vị trí không gian dọc trục hầm ($H$) vào hàm xác định hệ số $K$ và $\alpha$. Luận án chứng minh khối đá có RMR càng cao thì mức độ truyền năng lượng chấn động càng lớn và tốc độ suy giảm năng lượng càng chậm, thiết lập một quy luật định lượng hoàn toàn mới thay thế các hệ số kinh nghiệm thuần túy trước đây.

2. Điểm đột phá về mặt phương pháp luận so với các nghiên cứu quốc tế trước đây là gì? Trả lời: So với nghiên cứu của Ansell (2004) chỉ sử dụng mô hình dầm-lò xo 2D hay Li et al. (2013) dùng thông số cơ học tĩnh trong mô hình động, luận án đã đột phá khi: (1) Sử dụng thông số cơ học động ($E_d, \dot{\varepsilon}$) đo trực tiếp từ thí nghiệm nén động SHPB tốc độ cao; (2) Triển khai mô hình phần tử hữu hạn 3D không gian thực tích hợp biên vô hạn không phản xạ NRBC (Infinite Elements) trên Abaqus; (3) Đối soát kiểm chứng đa chiều với 54 vụ nổ thực tế tại hầm Croix-Rousse (Pháp).

3. Phát hiện bất ngờ nhất về mặt dữ liệu thực nghiệm là gì? Trả lời: Phát hiện bất ngờ nhất là tính bất đối xứng của trường chấn động dọc trục hầm $|H|$: Tại cùng một khoảng cách hình học tính từ mặt phẳng gương nổ mìn, vùng vỏ hầm cũ nằm phía sau gương nổ ($H < 0$) có biên độ PPV thấp hơn và mức độ suy giảm năng lượng nhanh hơn $18{,}5%$ so với vùng nằm phía trước gương nổ ($H > 0$). Điều này chứng minh khoảng không gian rỗng đã đào của đường hầm đóng vai trò như một bức tường chắn sóng tự nhiên, phân tán năng lượng địa chấn hiệu quả.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) rõ ràng không? Trả lời: Luận án cung cấp quy trình tái lập vô cùng chi tiết và minh bạch, bao gồm: Kích thước hình học mẫu đá chuẩn ISRM ($\Phi 45 \times 100\text{ mm}$), thông số vận hành hệ thống thanh nén SHPB (áp lực khí nén $0{,}30 - 0{,}35\text{ MPa}$), sơ đồ bố trí cảm biến Geophone 3 trục, thông số phương trình trạng thái vật liệu nổ TNT, cấu hình lưới phần tử C3D8R và biên CIN3D8, cùng thuật toán tích phân thời gian Explicit trong phần mềm Abaqus CAE.

5. Lộ trình nghiên cứu 10 năm tiếp theo được định hình như thế nào? Trả lời: Lộ trình 10 năm tới định hướng: (1) Tích hợp mô hình đứt gãy rời rạc FEM-DEM đánh giá đá phân lớp nứt nẻ mạnh; (2) Nghiên cứu hiệu ứng mỏi tích lũy của vỏ bê tông dưới hàng nghìn chu kỳ tải trọng động vi sai; (3) Xây dựng hệ thống bản sao số (Digital Twin) tích hợp AI giám sát chấn động nổ mìn thời gian thực cho hệ thống hạ tầng ngầm đô thị thông minh.

Kết luận

Luận án tiến sĩ kỹ thuật của NCS. Đặng Văn Kiên là một công trình nghiên cứu khoa học công phu, mẫu mực và chuẩn xác cao, đem lại 6 đóng góp cốt lõi:

  1. Thực chứng hóa các thông số cơ học động: Xác định chính xác quan hệ ứng suất - biến dạng - tốc độ biến dạng của đá granit dưới tải trọng động bằng hệ thống thanh ép SHPB chuyên dụng tại INSA Lyon ($E_d = 62{,}5\text{ GPa}$, độ bền nén động tăng hơn $23%$).
  2. Thiết lập tương quan định lượng RMR - Chapot: Xây dựng hệ phương trình hồi quy xác định hệ số truyền sóng $K$ và suy giảm $\alpha$ phụ thuộc trực tiếp vào chỉ số chất lượng khối đá RMR và vị trí tương đối dọc trục hầm $H$.
  3. Tối ưu hóa mô hình số 3D Abaqus: Làm chủ công nghệ mô phỏng động lực học phi tuyến với biên vô hạn không phản xạ NRBC, xác lập hệ số giảm chấn chuẩn tắc $\xi = 5{,}0%$ cho khối đá granit.
  4. Khám phá quy luật phân bố chấn động không gian 3D: Chứng minh quy luật suy giảm bất đối xứng của PPV dọc trục hầm và xác định vị trí tập trung ứng suất kéo nguy hiểm nhất tại hông tường vỏ hầm cũ đối diện gương nổ.
  5. Ứng dụng thành công chỉ số tổn thương BDI: Đánh giá định lượng trạng thái an toàn của kết cấu vỏ bê tông, làm sáng tỏ cơ chế nứt giòn do ứng suất kéo động sinh ra.
  6. Ban hành công thức thực hành kỹ thuật: Đưa ra công thức tính toán lượng thuốc nổ cho phép lớn nhất $Q_{max}$ áp dụng hiệu quả cho các dự án hầm song song tại Việt Nam (hầm Cổ Mã, hầm Hải Vân), bảo đảm tuyệt đối an toàn kết cấu và nâng cao hiệu quả kinh tế - kỹ thuật trong xây dựng công trình ngầm quốc gia.