Tổng quan về luận án

Sự bùng nổ của các dự án hạ tầng giao thông ngầm đô thị và xuyên núi tại Việt Nam cũng như trên thế giới đặt ra thách thức kỹ thuật chưa từng có: thi công các tuyến hầm mới bằng phương pháp khoan nổ mìn trong khoảng cách cực kỳ hạn chế với các đường hầm hiện hữu (như tuyến hầm song song tại đèo Cổ Mã cách hầm đường sắt cũ 47,0 m, hay mở rộng hầm lánh nạn đèo Hải Vân cách hầm chính chỉ 30,0 m). Khi nổ mìn tại gương hầm, chỉ một phần năng lượng được dùng để phá vỡ đất đá, phần còn lại truyền vào môi trường địa chất dưới dạng các xung sóng ứng suất động, tác động trực tiếp và đe dọa sự toàn vẹn của kết cấu vỏ chống bê tông liền khối lân cận. Luận án tiến sĩ của tác giả Đặng Văn Kiên 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" (Trường Đại học Mỏ - Địa chất, 2018) đã giải quyết trọn vẹn bài toán cơ học công trình ngầm phức tạp này thông qua tiếp cận tích hợp đa phương pháp: đo đạc thực nghiệm hiện trường, thí nghiệm cơ học động trong phòng và mô phỏng số ba chiều (3D FEM).

Khoảng trống nghiên cứu (research gap) cốt lõi được nhận diện là sự thiếu hụt các mô hình dự báo không gian 3D chính xác về trường vận tốc dao động phần tử đỉnh (Peak Particle Velocity - PPV) và mức độ tổn thương kết cấu bê tông cũ chịu tải trọng nổ vi sai. Các nghiên cứu trước đây (như Jong-Ho Shin et al., 2011; Hoàng Thị Hồng, 2004; Nguyễn Hữu Thế, 2014) chủ yếu dừng lại ở các bài toán giải tích phẳng 2D, áp dụng thông số cơ lý tĩnh cho bài toán động lực học, hoặc áp dụng cho công trình quân sự chịu nổ đơn độc trên mặt đất mà chưa tính toán tương tác động học 3D dọc trục đường hầm trong điều kiện khối đá nứt nẻ thực tế.

Luận án thiết lập 4 câu hỏi nghiên cứu (RQ) và 4 giả thuyết khoa học (H) tương ứng:

  • RQ1: Quy luật lan truyền sóng chấn động nổ mìn và suy giảm PPV trong môi trường đá nứt nẻ phụ thuộc như thế nào vào chỉ số phân cấp chất lượng khối đá (RMR) và tỷ lệ khoảng cách nạp thuốc ($D/\sqrt{Q}$)?
    H1: Tồn tại mối tương quan phi tuyến chặt chẽ giữa chỉ số RMR với các tham số suy giảm sóng địa chấn ($K$ và $\alpha$) trong công thức thực nghiệm của Chapot.
  • RQ2: Ứng xử cơ học động (dynamic stress-strain response) và độ bền của đá granit dưới tác dụng của tốc độ biến dạng cao biến đổi ra sao so với trạng thái tĩnh?
    H2: Mô đun đàn hồi động ($E_d$) và cường độ chịu nén động của đá granit tăng đáng kể dưới tải trọng xung động va đập trong thí nghiệm thanh đo áp lực Hopkinson phân tách (Split Hopkinson Pressure Bar - SHPB).
  • RQ3: Làm thế nào để loại trừ hoàn toàn hiện tượng phản xạ sóng tại biên mô hình số và lựa chọn chính xác hệ số giảm chấn khi giải bài toán động lực học đường hầm?
    H3: Tích hợp phần tử vô hạn (Infinite Element - IE) làm điều kiện biên không phản xạ (Non-Reflecting Boundary Condition - NRBC) kết hợp hệ số giảm chấn Rayleigh $\xi = 5.0%$ cho phép mô phỏng chính xác đáp ứng động học so với dữ liệu đo đạc thực tế.
  • RQ4: Mức độ hư hại và vùng phá hủy của vỏ bê tông hầm cũ được phân bố như thế nào dọc theo trục không gian 3 chiều khi gương hầm mới tiến triển?
    H4: Vận tốc dao động phần tử cực đại trên vỏ hầm cũ đạt đỉnh tại vị trí đối diện gương hầm mới ($D^- = 0\text{ m}$) và suy giảm nhanh theo quy luật hàm mũ dọc theo phương trục hầm, với ngưỡng phá hủy được định lượng chính xác qua chỉ số phá hủy nổ mìn (Blast Damage Index - BDI).

Khung lý thuyết của nghiên cứu được xây dựng trên nền tảng của Lý thuyết đàn hồi động học (Dynamic Theory of Elasticity), Lý thuyết lan truyền sóng ứng suất trong môi trường liên tục (Continuum Stress Wave Propagation), Cơ học phá hủy động (Dynamic Fracture Mechanics) và Tiêu chuẩn đánh giá mức độ tổn hại chấn động quốc tế (DIN 4150-3, AS 2187, GB 6722-2003).

Đóng góp đột phá của công trình được định lượng cụ thể: Xây dựng thành công chuỗi công thức kinh nghiệm dự báo PPV theo chỉ số RMR cho các phương dao động; xác lập hệ số giảm chấn tối ưu $\xi = 5.0%$ cho mô hình số phần tử hữu hạn; và thiết lập bản đồ phân vùng phá hủy vỏ hầm cũ bằng chỉ số BDI dọc theo trục không gian 3D. Quy mô nghiên cứu bao gồm việc xử lý dữ liệu đo đạc từ 54 vụ nổ thực nghiệm tại dự án hầm Croix-Rousse (Lyon, Pháp) với chiều dài 1757,5 m, diện tích mặt cắt $84,10\text{ m}^2$, tiến hành hàng chục mẫu thử nghiệm SHPB trên đá granit đường kính 45,0 mm tại INSA Lyon, và xây dựng mô hình số 3D quy mô lớn trên phần mềm Abaqus CAE.

Literature Review và Positioning

Nghiên cứu cơ chế lan truyền sóng nổ và tác động động lực học của nó lên kết cấu công trình ngầm đã trải qua nhiều giai đoạn phát triển với những luồng quan điểm học thuật đối trọng:

  [Trường phái Giải tích Cổ điển]
  - Pao & Mow (1973), Wersall (2008), Li et al. (2013)
  - Mô hình hóa giải tích dạng sóng đàn hồi 2D (Toán tử Laplace-Fourier)
  - Giới hạn: Chỉ giải quyết được hình học tròn lý tưởng, môi trường đẳng hướng đồng nhất
  [Trường phái Thực nghiệm & Bán kinh nghiệm]
  - Chapot (1981), Sadovsky (1974), USBM, Ansell (2004), Ahmed & Ansell (2012)
  - Hàm suy giảm PPV theo tỷ lệ nạp thuốc (Scaled Distance D/√Q)
  - Tranh luận đối lập: Tác động nổ đơn độc vs. Nổ vi sai nhiều lỗ mìn
  [Trường phái Mô hình hóa Số học Tiên tiến]
  - Lu et al. (2005), Jong-Ho Shin et al. (2011), Xia et al. (2013), Liang (2012)
  - Mô phỏng 2D/3D FEM/DEM, đánh giá sóng phản xạ biên & tải trọng nổ theo thời gian
  - Hạn chế: Thường dùng thông số tĩnh, bỏ qua tương tác không gian 3D của gương hầm
  [Vị thế của Luận án Đặng Văn Kiên (2018)]
  - Cầu nối tích hợp: SHPB xác định thông số ĐỘNG + Mô hình 3D Abaqus với Biên vô hạn NRBC
  - Hiệu chỉnh bằng 54 vụ nổ thực nghiệm tại Croix-Rousse + Liên kết tường minh RMR vào PPV

Trường phái giải tích cổ điển (Pao Y.H. & Mow C.C., 1973; Wersall C., 2008; Li et al., 2013) sử dụng toán tử tích phân Laplace và chuỗi Fourier để tìm nghiệm dạng đóng cho đường hầm mặt cắt tròn chịu sóng nổ phẳng. Dù mang tính chặt chẽ về toán học, trường phái này bộc lộ điểm nghẽn nghiêm trọng khi không thể áp dụng cho các đường hầm có tiết diện móng ngựa/vòm tường thẳng phức tạp, hoặc khối đá có tính dị hướng và nứt nẻ tự nhiên.

Trường phái thực nghiệm hiện trường (Ansell, 2004; Ahmed & Ansell, 2012; Chapot, 1981) tập trung xác lập các tương quan bán thực nghiệm giữa PPV với khoảng cách $D$ và lượng thuốc nổ $Q$. Trong đó, tồn tại cuộc tranh luận học thuật sâu sắc: một nhóm tác giả (như Lu Y. et al., 2005; Saharan & Mitri, 2008) cho rằng đáp ứng động lực học của kết cấu phụ thuộc tuyệt đối vào tổng lượng thuốc nổ nạp trên toàn gương hầm; ngược lại, nhóm nghiên cứu do Xia et al. (2013) và Liang Q. (2012) dẫn đầu chứng minh rằng chấn động cực đại bị chi phối chủ yếu bởi khối lượng thuốc nổ lớn nhất kích nổ trong một đợt vi sai ($Q_{max}$, thông thường trong khoảng thời gian phân cách $\Delta t = 8,0\text{ ms}$) tại các lỗ mìn khoét đột trung tâm.

Tại Việt Nam, các nghiên cứu của Hoàng Thị Hồng (2004) tiếp cận theo mô hình lưu biến Kelvin và Poynting-Thomson với buồng nổ hình trụ; Nguyễn Xuân Mãn (2000) giải bài toán đàn dẻo xác định bán kính vùng phá hủy của khối thuốc đơn lẻ; Lê Đình Tân (2012) và Nguyễn Hữu Thế (2014) khảo sát kết cấu hầm quân sự chịu tải trọng nổ bom trên mặt đất hoặc trong nền đá san hô bằng phần mềm Plaxis và Autodyn/MineBlast. Điểm hạn chế chung của các nghiên cứu trong nước là việc áp dụng các thông số cơ lý tĩnh của đất đá vào phương trình dao động động học, sử dụng điều kiện biên gối tựa cố định gây phản xạ ngược sóng nổ vào miền tính toán, và chưa có nghiên cứu chuyên sâu nào đánh giá ảnh hưởng của việc thi công hầm đôi song song trong môi trường đô thị bằng phương pháp khoan nổ mìn.

Luận án của Đặng Văn Kiên định vị chính xác vị thế học thuật bằng cách kết nối trực tiếp dữ liệu đo đạc thực nghiệm từ dự án hầm đô thị Croix-Rousse (Lyon, Pháp) và thí nghiệm cơ học động SHPB tại phòng thí nghiệm INSA Lyon để xây dựng mô hình số 3D hoàn chỉnh. Công trình đặt trong mối tương quan so sánh đối sánh với hai nghiên cứu quốc tế điển hình:

  1. Nghiên cứu hầm tàu điện ngầm đa tầng tại Quảng Châu của Lu et al. (2005): Lu et al. sử dụng mô hình LS-DYNA nhưng bỏ qua liên kết tiếp xúc cứng giữa vỏ hầm và khối đá, dẫn đến việc ước tính sai lệch ứng suất kéo tiếp xúc; trong khi luận án đã phân định rõ điều kiện làm việc của "lớp lát hàn" (fully bonded composite interface) với giả thiết vỏ bê tông liền khối bám chặt và dao động đồng pha với khối đá granit rắn chắc.
  2. Nghiên cứu hầm đường sắt Xinjiang của Liang Q. et al. (2012): Liang et al. chỉ khảo sát giá trị PPV trên mặt cắt 2D ngang qua tâm nổ; luận án đã mở rộng toàn diện sang không gian 3D, làm sáng tỏ sự biến thiên của trường vận tốc dao động dọc theo trục đường hầm cũ khi gương hầm mới tịnh tiến từ khoảng cách $D^- = -30\text{ m}$ đến $D^+ = +30\text{ m}$.

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

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

Luận án đóng góp trực tiếp vào việc mở rộng Lý thuyết lan truyền sóng trong môi trường đàn dẻo nứt nẻ và Lý thuyết cơ học tiếp xúc công trình ngầm:

  1. Mở rộng định luật thực nghiệm Chapot - Sadovsky: Tác giả đã chứng minh và lượng hóa sự phụ thuộc của hai hệ số suy giảm sóng nổ $K$ và $\alpha$ vào chỉ số chất lượng khối đá RMR:

$$\ln(\text{PPV}) = \ln(K) - \alpha \ln\left(\frac{D}{\sqrt{Q}}\right)$$

Thay vì coi $K$ và $\alpha$ là các hằng số địa chất cố định, nghiên cứu thiết lập hàm phụ thuộc tường minh: $\ln(K) = f(\text{RMR})$ và $\alpha = g(\text{RMR})$, chỉ ra rằng khối đá có chất lượng càng cao (RMR lớn, ít nứt nẻ) thì khả năng truyền dẫn năng lượng sóng chấn động càng hoàn hảo, dẫn đến biên độ PPV tác động lên vỏ hầm cũ cao hơn rõ rệt so với môi trường đá phong hóa nứt nẻ mạnh. 2. Xác lập mô hình tương tác động học hệ "Khối đá - Vỏ hầm liền khối": Nghiên cứu mở rộng lý thuyết tính toán kết cấu ngầm chịu sóng nổ của Bulichev (1982) và Kartozia (1986) bằng việc xác định tường minh các điều kiện biên của trường hợp "lớp lát hàn": $E_0 > E$ (mô đun đàn hồi của đá lớn hơn vỏ bê tông) và biến dạng của vỏ chống tuân thủ hoàn toàn phương trình truyền sóng dọc phẳng:

$$\frac{\partial^2 u}{\partial t^2} = \frac{E}{\rho} \frac{\partial^2 u}{\partial x^2}$$

  1. Chuyển dịch mô thức (Paradigm Shift) trong phân tích động học công trình ngầm: Bẻ gãy quan niệm truyền thống coi thông số tĩnh có thể đại diện cho bài toán nổ mìn, luận án chứng minh rằng việc áp dụng mô đun đàn hồi động $E_d$ xác định từ sóng truyền siêu âm và thí nghiệm SHPB ($E_d \approx 1,2 \div 1,4 E_{tĩnh}$) là điều kiện tiên quyết để triệt tiêu sai số trong dự báo ứng suất phá hủy.

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

Khung phân tích của luận án là sự hợp nhất của ba trụ cột lý thuyết: (1) Lý thuyết truyền sóng đàn dẻo trong môi trường vô hạn; (2) Cơ học vật liệu dưới tốc độ biến dạng cao (High Strain-Rate Mechanics); và (3) Phương pháp phần tử hữu hạn động lượng nâng cao (Advanced Dynamic FEM).

Phương pháp tiếp cận giải quyết triệt để ranh giới biên (boundary conditions) trong bài toán động bằng cách tích hợp các phần tử vô hạn (Infinite Elements - IE) bao bọc toàn bộ biên ngoài của mô hình phần tử hữu hạn 3D. Giải pháp này triệt tiêu hoàn toàn hiện tượng phản xạ sóng nhân tạo tại biên mô hình trở lại kết cấu hầm, điều mà các nghiên cứu trước đây tại Việt Nam sử dụng biên gối tựa cố định hoặc biên trượt không thể xử lý được.

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái Thực chứng diễn dịch (Positivist-Deductive Paradigm) kết hợp chặt chẽ giữa nghiên cứu thực nghiệm vật lý (Empirical Physical Testing) và mô hình hóa toán học số trị (Numerical Simulation). Cấu trúc nghiên cứu đa cấp độ (Multi-level Design) bao gồm:

  • Cấp độ vĩ mô (Macro-scale): Quan trắc hiện trường toàn diện phản ứng động học của vỏ hầm giao thông Croix-Rousse dưới tác động của các đợt nổ vi sai bậc mili-giây.
  • Cấp độ vi mô/trung mô (Meso-scale): Đánh giá đặc tính động học, cơ chế phá hủy và mối quan hệ giữa tốc độ biến dạng ($\dot{\varepsilon}$) với ứng suất đỉnh thông qua hệ thống thanh đo áp lực Split Hopkinson Pressure Bar (SHPB).
  • Cấp độ mô hình hóa liên tục (Continuum-scale): Tái tạo không gian trường sóng và phân bố biến dạng bằng phương pháp phần tử hữu hạn (FEM) phi tuyến trong miền thời gian thực.

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu được kiểm soát nghiêm ngặt theo các tiêu chuẩn quốc tế:

  1. Giao thức quan trắc hiện trường: Sử dụng hệ thống cảm biến vận tốc địa chấn (Geophones) ba thành phần (Thẳng đứng - V, Ngang - H, Dọc trục - L) gắn cố định trên vỏ bê tông hầm cũ tại các lý trình PM 100, PM 220 $\div$ PM 340 và PM 1400. Dữ liệu được ghi nhận với tần số lấy mẫu cao, phân tích qua dải tần số thấp $f = 130\text{ Hz}$ để lọc nhiễu sóng âm và xác định giá trị PPV tức thời.
  2. Giao thức thí nghiệm SHPB tại INSA Lyon: Mẫu đá granit tự nhiên thu thập từ tuyến hầm được gia công hình trụ tiêu chuẩn có kích thước đường kính $D = 45,0\text{ mm}$, chiều cao $L = 100,0\text{ mm}$ (tỷ lệ $L/D \approx 2,2$). Thanh đánh (striker), thanh tới (incident bar) và thanh truyền (transmitter bar) được chế tạo từ hợp kim thép cường độ cao có mô đun đàn hồi $E_b = 190\text{ GPa}$, vận tốc sóng âm trong thanh $C_0 = 5188\text{ m/s}$. Áp lực khí nén đẩy thanh đánh được thí nghiệm ở hai cấp áp lực: $0,30\text{ MPa}$ và $0,35\text{ MPa}$. Cầu điện trở Wheatstone cùng cảm biến đo biến dạng (strain gauges) ghi nhận tức thời tín hiệu điện áp của sóng tới ($\varepsilon_i$), sóng phản xạ ($\varepsilon_r$) và sóng truyền qua ($\varepsilon_t$).
  3. Kỹ thuật tam giác đạc (Triangulation): Kết quả phản ứng cơ học được đối chiếu chéo đồng thời giữa: (a) Dữ liệu đo Geophone thực tế tại hiện trường; (b) Giá trị thu được từ thí nghiệm SHPB; và (c) Kết quả xuất ra từ mô hình số Abaqus. Độ tin cậy và tính hợp thức (construct validity) được chứng minh khi sai số tương đối giữa mô hình số và thực nghiệm đo đạc chỉ dao động dưới $8,5%$.

Data và phân tích

Phân tích số trị được thực hiện trên phần mềm Abaqus/Explicit với quy trình kiểm soát lưới nghiêm ngặt. Để tránh hiện tượng tán sắc số và đảm bảo tái hiện chính xác bước sóng địa chấn ngắn nhất ($\lambda_{min}$), kích thước phần tử hữu hạn ($l_e$) được khống chế chặt chẽ theo điều kiện Kuhn (1979):

$$l_e \le \left(\frac{1}{8} \div \frac{1}{10}\right) \frac{C_s}{f_{max}}$$

Trong mô hình số 3D khảo sát tương tác giữa đường hầm mới (bán kính $R_1 = 5,55\text{ m}$) và đường hầm cũ (bán kính vòm $R_2 = 8,05\text{ m}$, chiều cao tường $1,0\text{ m}$, cách nhau $42,6\text{ m}$ theo trục), toàn bộ miền tính toán được rời rạc hóa với hàng trăm nghìn phần tử khối 8 nút C3D8R kết hợp lớp phần tử vô hạn CIN3D8. Tải trọng áp lực nổ tại gương hầm được lập trình theo hàm suy giảm thời gian thực tế của đợt nổ vi sai với áp lực đỉnh tính toán từ phương trình trạng thái JWL của thuốc nổ.

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

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

  1. Định lượng chính xác ảnh hưởng của hệ số giảm chấn $\xi$ đến giá trị PPV: Mô hình số chứng minh giá trị PPV xuất hiện trong vỏ bê tông hầm cũ tỷ lệ nghịch với hệ số giảm chấn của khối đá. Khi hệ số giảm chấn tăng từ $\xi = 2,0%$ lên $\xi = 10,0%$, biên độ đỉnh PPV ghi nhận tại các điểm quan sát A và B giảm hơn $40%$. Luận án chứng minh hệ số giảm chấn Rayleigh $\xi = 5,0%$ là giá trị tối ưu, phản ánh trung thực nhất sự tiêu tán năng lượng trong đá granit khối nứt nẻ nhẹ.
  2. Quy luật chấn động nghịch đảo theo chất lượng khối đá RMR: Phân tích dữ liệu từ 54 vụ nổ tại Croix-Rousse chỉ ra rằng trong khối đá có chỉ số RMR cao (đá rắn chắc, ít khe nứt), năng lượng sóng nổ ít bị suy hao do tán xạ và ma sát khe nứt, khiến giá trị PPV truyền tới vỏ hầm cũ cao hơn đáng kể so với vùng đá có RMR thấp. Đây là kết quả có tính cảnh báo sâu sắc cho các kỹ sư thiết kế hầm: môi trường địa chất "tốt" lại tiềm ẩn rủi ro chấn động nổ mìn phá hủy công trình lân cận cao hơn.
  3. Hiệu ứng gia cường động học của đá granit qua SHPB: Dưới tác dụng của tốc độ biến dạng cao ($\dot{\varepsilon} = 10^2 \div 10^3\text{ s}^{-1}$), mô đun đàn hồi động $E_d$ của đá granit đạt giá trị cao hơn $20% \div 35%$ so với mô đun tĩnh thông thường. Tốc độ biến dạng và phản ứng ứng suất trong mẫu N°38 biến thiên tức thời theo áp lực xung kích của thanh đánh ($0,30\text{ MPa}$ và $0,35\text{ MPa}$), thiết lập cơ sở dữ liệu vật liệu động học chính xác cho Abaqus.
  4. Quy luật phân bố chấn động 3D dọc trục hầm: Khảo sát khoảng cách từ vị trí gương hầm mới đến điểm quan sát trên vỏ hầm cũ theo phương dọc trục ($D^-$) cho thấy: Chấn động lớn nhất không xuất hiện phân tán mà tập trung cục bộ tại tiết diện đối diện trực tiếp với gương nổ ($D^- = 0\text{ m}$), sau đó suy giảm theo hàm mũ khi khoảng cách dọc trục vượt quá $\pm 15\text{ m}$ ($D^- > 15\text{ m}$ hoặc $D^- < -15\text{ m}$).
  5. Xác lập ranh giới phá hủy vỏ bê tông theo Chỉ số BDI: Áp dụng chỉ số phá hủy nổ mìn BDI kết hợp tiêu chuẩn phá hủy bê tông, nghiên cứu chỉ ra giá trị ngưỡng PPV giới hạn cho vỏ bê tông hầm cũ là $[PPV] = 0,22\text{ m/s}$ ($220\text{ mm/s}$). Vùng phá hủy nguy hiểm nhất do ứng suất kéo tập trung cục bộ tại vị trí chân tường vòm phía đối diện trực tiếp với hướng truyền sóng nổ.

Implications đa chiều

  • Ý nghĩa học thuật: Bổ sung cơ sở dữ liệu cơ học động học cho khối đá granit chịu nổ; hoàn thiện khung lý thuyết dự báo sóng nổ không gian 3D kết hợp chỉ số phân loại địa kỹ thuật RMR.
  • Ý nghĩa phương pháp luận: Chuẩn hóa quy trình mô phỏng số động lực học công trình ngầm bằng cách tích hợp phần tử biên vô hạn (IE) và hiệu chuẩn thông số giảm chấn Rayleigh trong phần mềm Abaqus, mở ra phương pháp luận có thể chuyển giao cho các bài toán phân tích chấn động khác (như động đất, va đập, nổ công sự).
  • Ứng dụng thực tiễn công trình: Cung cấp công thức tính toán lượng thuốc nổ nạp lớn nhất cho phép cho một đợt nổ ($Q_{max}$) tương ứng với khoảng cách thi công $D$ và chỉ số RMR của từng đoạn tuyến, đảm bảo an toàn tuyệt đối cho các công trình lân cận tại các dự án hầm Hải Vân 2, Cổ Mã, Đèo Cả và các tuyến metro ngầm tại Hà Nội, TP. Hồ Chí Minh.
  • Khuyến nghị chính sách và quy chuẩn: Đề xuất bổ sung quy định về giá trị giới hạn PPV cho kết cấu hầm ngầm trong Quy chuẩn Kỹ thuật Quốc gia về an toàn nổ mìn tại Việt Nam, dựa trên sự hài hòa giữa các tiêu chuẩn quốc tế tiên tiến như DIN 4150-3 (Đức) và AS 2187 (Úc).

Limitations và Future Research

Nghiên cứu thẳng thắn thừa nhận 4 giới hạn khoa học nội tại:

  1. Giả thiết tính liên tục của môi trường đá: Mô hình số xem môi trường đất đá là đàn hồi tuyến tính, đồng nhất và đẳng hướng trong phạm vi cục bộ, chưa xét đến hệ thống khe nứt ngẫu nhiên 3D và các đới đứt gãy kiến tạo phức tạp bằng phương pháp phần tử rời rạc (DEM/DDA).
  2. Mô hình hóa áp lực nổ gián tiếp: Áp lực nổ được quy đổi thành xung áp lực tác dụng lên bề mặt biên đường hầm theo hàm thời gian, chưa mô phỏng trực tiếp tương tác khí động học dòng nổ đa pha và cơ chế nứt vỡ khí nén giãn nở (gas expansion fracture) từ bên trong từng lỗ mìn.
  3. Giả định tiếp xúc vỏ hầm: Xem mối liên kết giữa vỏ bê tông cũ và khối đá là "lớp lát hàn" bám dính tuyệt đối, chưa khảo sát các trường hợp vỏ hầm lắp ghép phân đoạn (segmental lining) có sự trượt tương đối hoặc tồn tại khoảng rỗng sau lưng vỏ hầm.
  4. Giới hạn loại đá thí nghiệm: Cơ sở dữ liệu động học SHPB mới tập trung chủ yếu trên mẫu đá granit đặc sít tại khu vực dự án nghiên cứu, chưa bao phủ các loại đá trầm tích hoặc biến chất có tính dị hướng cao (như đá phiến sét, đá vôi nứt nẻ).

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

  • Phát triển mô hình liên hợp phần tử hữu hạn - phần tử rời rạc (FEM-DEM/Particle Flow Code) để mô phỏng sự hình thành mạng lưới nứt vỡ vi mô trong khối đá bất liên tục.
  • Ứng dụng trí tuệ nhân tạo (Machine Learning/Deep Neural Networks) để tối ưu hóa mạng lưới lỗ nổ vi sai theo thời gian thực dựa trên tín hiệu quan trắc Geophone phản hồi.
  • Nghiên cứu cơ chế suy thoái mỏi động học (Dynamic Fatigue Degradation) của vỏ bê tông cốt thép hầm cũ dưới tác động của hàng nghìn chu kỳ nổ mìn lặp lại liên tục trong suốt quá trình đào hầm mới.
  • Mở rộng phân tích tương tác thủy địa cơ động (Hydro-Mechanical Dynamic Interaction) khi có sự hiện diện của áp lực nước ngầm cao trong khe nứt.

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

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

  • Tác động học thuật: Dự báo tạo tiền đề cho hàng loạt công bố khoa học chuyên ngành trên các tạp chí quốc tế uy tín (ISI/Scopus mảng Tunnelling and Underground Space Technology, Rock Mechanics and Rock Engineering); thiết lập khung tham chiếu phương pháp luận mô phỏng chấn động nổ mìn 3D cho các cơ sở đào tạo tiến sĩ ngành Kỹ thuật Xây dựng Công trình ngầm và Mỏ.
  • Chuyển dịch công nghiệp hạ tầng: Ứng dụng trực tiếp giải pháp kỹ thuật kiểm soát nổ mìn vi sai và công thức dự báo $Q_{max}$ tại các đại dự án hạ tầng giao thông huyết mạch của Việt Nam, giúp các nhà thầu thi công rút ngắn thời gian đào hầm, giảm thiểu chi phí gia cố bảo vệ công trình cũ, và triệt tiêu nguy cơ sập đổ sự cố có thể gây thiệt hại hàng trăm tỷ đồng.
  • Tác động an toàn và xã hội: Bảo vệ an toàn tuyệt đối cho các công trình dân dụng trên mặt đất và các công trình giao thông huyết mạch ngầm đang khai thác; loại trừ nguy cơ nứt nẻ nhà dân và công trình công cộng tại các khu vực đô thị có mật độ dân cư dày đặc khi thi công tuyến hầm mới bên dưới.

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

  • Nghiên cứu sinh và Giảng viên ngành Công trình ngầm: Tiếp cận khung phương pháp luận kết hợp chuẩn mực giữa đo đạc hiện trường, thí nghiệm động học SHPB cao cấp và kỹ thuật mô hình số Abaqus/Explicit 3D loại trừ phản xạ biên.
  • Tổng thầu và Kỹ sư nổ mìn công trình ngầm: Nhận chuyển giao trực tiếp bộ công thức tính toán lượng nạp thuốc an toàn $Q_{max}$, sơ đồ bố trí kíp nổ vi sai tối ưu và quy trình giám sát chấn động PPV tại hiện trường.
  • Cơ quan quản lý nhà nước và Ban quản lý dự án: Sở hữu cơ sở khoa học tin cậy để thẩm định biện pháp thi công nổ mìn, ban hành quy chuẩn an toàn chấn động kỹ thuật số cho các công trình ngầm cận kề trong tương lai.

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?

Đóng góp độc đáo nhất là việc mở rộng định luật thực nghiệm lan truyền sóng nổ Chapot - Sadovsky thông qua việc lượng hóa tường minh mối quan hệ giữa chất lượng khối đá theo chỉ số RMR với các thông số truyền sóng ($K, \alpha$). Nghiên cứu chứng minh quy luật nghịch lý thực nghiệm: Khối đá có RMR càng cao (chất lượng cơ lý tốt, ít khe nứt) thì khả năng bảo toàn năng lượng sóng nổ càng lớn, khiến vận tốc dao động phần tử đỉnh PPV tác động lên vỏ hầm cũ càng gia tăng.

2. Điểm đột phá về phương pháp luận so với các nghiên cứu trước đây là gì?

So với các nghiên cứu của Jong-Ho Shin et al. (2011) và Liang Q. (2012), luận án đã: (1) Thay thế hoàn toàn các thông số đàn dẻo tĩnh bằng các thông số động học ($E_d$, tốc độ biến dạng) trích xuất trực tiếp từ thí nghiệm va đập SHPB; (2) Giải quyết triệt để vấn đề phản xạ sóng tại biên trong mô hình FEM bằng việc tích hợp hệ thống phần tử vô hạn CIN3D8; và (3) Mở rộng bài toán từ mặt phẳng 2D sang không gian 3D toàn diện dọc theo trục hầm.

3. Phát hiện nào mang tính bất ngờ và phản trực giác nhất trong dữ liệu?

Phát hiện phản trực giác nhất là sự tập trung cục bộ cao độ của năng lượng sóng chấn động: Biên độ dao động PPV trên vỏ hầm cũ đạt cực đại tuyệt đối chỉ trong phạm vi hẹp $\pm 5\text{ m}$ tại mặt phẳng đối diện trực diện gương nổ ($D^- = 0\text{ m}$) và suy giảm đột ngột khi khoảng cách dọc trục vượt quá $15\text{ m}$. Điều này chứng minh rằng việc gia cố tăng cường vỏ hầm cũ chỉ cần tập trung cục bộ theo bước tiến của gương hầm thay vì phải xử lý tốn kém trên toàn bộ chiều dài tuyến hầm.

4. Giao thức tái lập nghiên cứu (Replication Protocol) có được cung cấp đầy đủ không?

Luận án cung cấp chi tiết toàn bộ thông số đầu vào: Kích thước hình học hầm Croix-Rousse ($R_1 = 5,55\text{ m}$, $R_2 = 8,05\text{ m}$, khoảng cách trục $42,6\text{ m}$); thông số thanh thử SHPB (đường kính $45\text{ mm}$, $E_b = 190\text{ GPa}$, $C_0 = 5188\text{ m/s}$); thông số vật liệu đá granit và bê tông; thuật toán ma trận giảm chấn Rayleigh ($\xi = 5,0%$); và phương trình áp lực tải trọng nổ theo thời gian, cho phép tái lập 100% mô hình trên Abaqus.

5. Lộ trình phát triển nghiên cứu 10 năm tới được phác thảo như thế nào?

Lộ trình hướng tới: (1) Số hóa bản đồ chấn động ngầm thời gian thực (Real-time Blast Digital Twin); (2) Phát triển các loại vật liệu vỏ chống composite hấp thụ sóng chấn động; và (3) Hoàn thiện hệ thống tự động hóa điều khiển bước nổ vi sai bằng trí tuệ nhân tạo thích ứng với sự thay đổi địa chất tức thời tại gương hầm.

Kết luận

Công trình nghiên cứu tiến sĩ của Đặng Văn Kiên xác lập 5 đóng góp học thuật và thực tiễn mang tính nền tảng:

  1. Thiết lập cơ sở dữ liệu động học hoàn chỉnh của đá granit tự nhiên dưới tải trọng xung kích tốc độ biến dạng cao thông qua hệ thống thử nghiệm chuyên sâu SHPB tại INSA Lyon.
  2. Xác lập hệ công thức kinh nghiệm mới dự báo chính xác vận tốc dao động phần tử đỉnh PPV theo ba phương không gian tích hợp chỉ số chất lượng khối đá RMR.
  3. Phát triển thành công mô hình số 3D FEM tối ưu trên Abaqus với kỹ thuật biên vô hạn NRBC khử phản xạ sóng và xác định chuẩn xác hệ số giảm chấn Rayleigh $\xi = 5,0%$.
  4. Làm sáng tỏ quy luật phân bố ứng suất - biến dạng động và ranh giới phá hủy nứt nẻ của vỏ bê tông hầm cũ theo chỉ số tổn thương nổ mìn BDI dọc theo trục không gian 3 chiều.
  5. Cung cấp quy trình công nghệ tính toán an toàn nổ mìn có khả năng ứng dụng thực tiễn cao cho các dự án xây dựng hầm song song phức tạp tại Việt Nam và trên thế giới.

Bằng việc kết hợp chặt chẽ giữa lý thuyết cơ học sóng đàn dẻo, thực nghiệm hiện trường đô thị và kỹ thuật mô phỏng số tiên tiến, luận án đã thúc đẩy một bước tiến quan trọng trong lĩnh vực Cơ học công trình ngầm và Công nghệ khai thác mỏ, mở ra hướng nghiên cứu liên ngành bền vững cho sự phát triển hạ tầng giao thông ngầm hiện đại.