Tổng quan về luận án

Công tác thi công các công trình hạ tầng ngầm quy mô lớn, bao gồm hệ thống hầm thủy điện, tuyến hầm giao thông đường bộ xuyên núi, cùng hệ thống công trình ngầm phòng thủ quân sự và hầm chứa máy bay chiến lược, giữ vị trí then chốt trong sự nghiệp công nghiệp hóa, hiện đại hóa và củng cố quốc phòng an ninh tại Việt Nam. Về mặt địa chất công trình, phần lớn các cấu trúc ngầm này phải xuyên qua các tầng địa tầng đá cứng phức tạp, nơi phương pháp khoan nổ mìn chu kỳ (drill and blast) là giải pháp kỹ thuật bắt buộc và chủ đạo. Trong chu kỳ đào ngầm, giai đoạn thi công lỗ khoan chiếm tỷ trọng thời gian rất lớn; theo các nghiên cứu chuyên ngành, "thời gian để khoan các lỗ thường chiếm từ 15-60% tổng thời gian của một chu kỳ đào". Do đó, hiệu suất phá vỡ đá và mức tiêu hao năng lượng của công tác khoan đóng vai trò quyết định đến tiến độ toàn bộ dự án, tuổi thọ trang thiết bị và tổng chi phí đầu tư.

Mặc dù các tập đoàn cơ khí hàng đầu thế giới như Atlas Copco, Furukawa, Tamrock (Sandvik) hay Caterpillar đã chế tạo nhiều dòng máy khoan tự hành xoay đập thủy lực hiện đại, khoảng trống nghiên cứu (research gap) lớn vẫn tồn tại khi áp dụng các thiết bị này vào điều kiện địa chất cụ thể của Việt Nam:

  1. Thiếu mô hình cơ học hệ nhiều vật (multibody dynamic system) hoàn chỉnh mô tả tương tác đồng thời ba chuyển động độc lập (dẫn động đập, dẫn động xoay và dẫn tiến) với môi trường đá có tính chất đàn - nhớt - dẻo (elasto-visco-plasticity).
  2. Các chế độ vận hành hiện tại chủ yếu dựa trên kinh nghiệm cảm tính của thợ vận hành hoặc thông số mặc định của nhà sản xuất nước ngoài, chưa tối ưu hóa theo chỉ tiêu chi phí năng lượng riêng nhỏ nhất ($E_s \to \min$) cho từng cấp độ cứng đá nội địa.
  3. Tỷ lệ sự cố kỹ thuật như kẹt cần, gãy choòng và mài mòn mũi khoan quá mức vẫn xảy ra thường xuyên do sự mất cân bằng động lực học giữa lực dẫn tiến ($F_{dt}$), mô men xoắn ($M_q$), lực đập ($F_đ$) và tần số đập ($f_đ$).

Nhằm giải quyết triệt để khoảng trống này, luận án tập trung giải quyết các câu hỏi nghiên cứu và giả thuyết khoa học sau:

  • Câu hỏi nghiên cứu 1 (RQ1): Cơ chế truyền tải và tương tác sóng ứng suất va đập kết hợp mô men cắt xoay giữa cụm búa - choòng - mũi khoan với khối đá có tính phi tuyến đàn - nhớt - dẻo được mô hình hóa toán học như thế nào?
  • Câu hỏi nghiên cứu 2 (RQ2): Mức độ ảnh hưởng tương hỗ và quy luật tác động đồng thời của các thông số động lực học ($F_{dt}, M_q, F_đ, f_đ, \omega$) tới tốc độ xuyên ($PR$) và chi phí năng lượng riêng ($E_s$) diễn ra ra sao?
  • Câu hỏi nghiên cứu 3 (RQ3): Miền thông số làm việc tối ưu cục bộ và toàn cục cho thiết bị khoan xoay đập thủy lực trong các tầng đá đại diện tại Việt Nam là gì?

Hệ thống giả thuyết nghiên cứu bao gồm:

  • Giả thuyết H1: Tồn tại một vùng giá trị lực dẫn tiến tối ưu $F_{dt}^*$ tương ứng với từng mức lực đập $F_đ$ và tần số $f_đ$, tại đó biên độ phản xạ sóng ứng suất nén/kéo ngược về búa khoan là nhỏ nhất, cực đại hóa năng lượng phá hủy hữu ích đi vào gương đá.
  • Giả thuyết H2: Chi phí năng lượng riêng phá vỡ thể tích đá ($E_s$) là hàm phi tuyến đa biến không đơn điệu theo vận tốc góc xoay $\omega$ và lực dẫn tiến $F_{dt}$; việc xác lập hàm mục tiêu tối ưu cho phép giảm tối thiểu 15–20% năng lượng tiêu hao vô ích so với chế độ khoan không kiểm soát.
  • Giả thuyết H3: Mô hình lưu biến Poynting-Thomson tích hợp ngưỡng phá hủy dẻo $F_{th}$ phản ánh chính xác trạng thái phá hủy của đá biến chất và đá trầm tích Việt Nam hơn các mô hình đàn hồi tuyến tính truyền thống.

Khung lý thuyết (theoretical framework) của luận án tích hợp:

  1. Cơ học hệ nhiều vật (Multibody Dynamics) ứng dụng phương trình vi phân Lagrange loại II;
  2. Lý thuyết lưu biến môi trường biến dạng đàn - nhớt - dẻo (Rheological Theory of Elasto-Visco-Plastic Media);
  3. Lý thuyết truyền sóng ứng suất một chiều trong thanh đàn hồi (1D Stress Wave Propagation Theory);
  4. Lý thuyết tối ưu hóa kỹ thuật (Engineering Optimization Theory).

Về mặt tác động đột phá định lượng, nghiên cứu cung cấp cơ sở khoa học để chuẩn hóa chế độ làm việc của máy khoan đá xoay đập thủy lực (điển hình như các mẫu Furukawa HCR1500-EDII, Boomer 352, Tamrock Axera D06, máy khoan KLM), góp phần nâng cao năng suất khai thác từ 2 đến 3 lần, giảm suất tiêu hao năng lượng riêng từ 12% đến 22%, đồng thời kiểm soát độ bền mỏi của chuỗi choòng khoan. Phạm vi nghiên cứu bao quát các mẫu đá chuẩn tại Việt Nam: đá granit (Bắc Bộ, Nam Bộ), đá sa thạch/cát kết, và "đặc biệt ở Hòa Bình chủ yếu là Đá trầm tích (đá vôi) có độ cứng cấp 3, cường độ chịu nén 1700 ÷ 2600 kg/cm2", với số liệu kiểm chứng thực nghiệm tại hiện trường thi công ngầm.

Literature Review và Positioning

Nghiên cứu về cơ chế phá hủy đá và động lực học thiết bị khoan xoay đập trên thế giới đã trải qua nhiều giai đoạn phát triển với ba dòng học thuật chính:

Dòng thứ nhất tập trung vào cơ chế phá hủy đá dưới tác dụng của tải trọng đơn lẻ (đập hoặc xoay). Шубный А.И. (1985) đã nghiên cứu quá trình phá hủy đá dưới tác dụng của chuyển động xoay và lực cắt thuần túy, chỉ ra rằng dưới áp lực nén dọc trục kết hợp mô men quay, mũi khoan ăn sâu vào đá và phá hủy liên kết hạt nhờ lực cắt tiếp tuyến. Ngược lại, Han Gang, Mike Bruno và Maurice B. (2001) đã ứng dụng phần mềm sai phân hữu hạn FLAC-3D để mô phỏng tương tác mũi khoan - đá dưới tải trọng đập chu kỳ, xem đá là môi trường đàn - dẻo tích lũy năng lượng va đập đến ngưỡng chảy dẻo. Sazidy và cộng sự (2010) phát triển mô hình khoan đập dựa trên môi trường đàn - nhớt - dẻo với ngưỡng phá vỡ tới hạn $F_{th}$, tuy nhiên mô hình này lược bỏ hoàn toàn thành phần vận tốc xoay và mô men cản cắt.

Dòng thứ hai nghiên cứu sự truyền sóng va đập và cộng hưởng phi tuyến trong chuỗi cần khoan. Marian Wiercigroch và cộng sự tại Viện Nghiên cứu Aberdeen - Scotland (2008) đã công bố giải pháp truyền năng lượng rung động tần số cao - biên độ thấp thành dao động cơ học tần số thấp - biên độ cao nhằm tối ưu hóa năng lượng va đập truyền vào đá. Ekaterina Pavlovskaia và Marian Wiercigroch (2003) khảo sát ảnh hưởng của tần số và biên độ lực đập tới động lực học phi tuyến của đá nhưng chưa tích hợp cơ chế xoay đồng bộ. Chiang và Dante (1998) sử dụng phương pháp phần tử hữu hạn 3D FEM mô tả tương tác mũi khoan - đá nhưng giả định ứng suất - biến dạng tuyến tính, không phản ánh được tính nhớt và từ biến thực tế của khối đá nứt nẻ.

Dòng thứ ba khảo sát mối quan hệ giữa cơ tính đá, khả năng khoan (drillability) và mô hình hóa liên hợp. Kahraman S. (1999) đã thiết lập phương trình thực nghiệm liên hệ giữa mô đun đá tỷ lệ ($MR = E/\sigma_c$) và tốc độ khoan ($PR$): $$PR = 0,14MR + 3,25$$ với kết luận rằng khi độ rỗng của đá dưới 1,23%, mô đun tỷ lệ kiểm soát chặt chẽ tốc độ khoan và khoan xoay đập là phương thức đạt hiệu quả cao nhất trong đá cứng. Timo Saksala và cộng sự (2011) đã xây dựng mô hình số 3D kết hợp mặt chảy dẻo Drucker-Prager, chỉ tiêu Rankine biến đổi và chóp nén parabol để mô tả sự phá hủy đập của đá, xem xét tính nhớt dẻo phụ thuộc tốc độ biến dạng.

Trong bức tranh tổng quan đó, hai luồng quan điểm đối lập nổi bật trong y văn quốc tế gồm:

  • Quan điểm A (Trường phái sóng va đập thuần túy): Cho rằng năng lượng đập ($F_đ, f_đ$) quyết định từ 80–90% thể tích đá bị phá vỡ, chuyển động xoay chỉ đóng vai trò dịch chuyển vị trí vấu khoan để đón nhịp đập kế tiếp và gom phoi vụn.
  • Quan điểm B (Trường phái tương tác liên hợp phi tuyến): Khẳng định lực cắt xoay ($M_q$) và lực ép dẫn tiến ($F_{dt}$) đóng vai trò tái cấu trúc trường ứng suất tập trung dưới đáy lỗ khoan, làm giảm đáng kể ngưỡng phá vỡ động $F_{th}$ của chu kỳ đập tiếp theo.

Vị trí học thuật của luận án được xác lập tại điểm giao thoa giữa cơ học hệ nhiều vật và cơ học địa chất thực nghiệm: Khắc phục nhược điểm tách rời của các nghiên cứu trước đây bằng cách xây dựng mô hình động lực học liên hợp đầy đủ 4 bậc tự do tịnh tiến và quay ($x_1, x_2, x_3, x_4$ và $\varphi_1, \varphi_2, \varphi_3$), gắn liền với cơ cấu giảm chấn thủy lực, khớp nối ren mềm đàn hồi, và môi trường đá Poynting-Thomson phi tuyến. Nghiên cứu so sánh trực tiếp với hai nghiên cứu quốc tế điển hình là mô hình giải tích của Sazidy et al. (2010, Memorial University of Newfoundland, Canada) và mô hình phần tử hữu hạn của Timo Saksala et al. (2011, Tampere University of Technology, Phần Lan), chứng minh tính ưu việt của phương pháp tối ưu hóa năng lượng riêng $E_s$ trên cơ sở tham số hóa đặc tính cơ lý đá tại Việt Nam.

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

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

Nghiên cứu mang lại những phát triển quan trọng đối với lý thuyết động lực học máy công trình ngầm và cơ học phá hủy đá:

  1. Mở rộng lý thuyết biến dạng phi tuyến của Baklasov I.A. và mô hình đàn - nhớt - dẻo: Luận án chuẩn hóa biểu đồ "ứng suất - biến dạng" tổng quát gồm 3 vùng: vùng trước giới hạn ($od$), vùng sau giới hạn ($de$) và vùng phá hủy vỡ vụn ($ek$). Mối quan hệ phi tuyến giữa ứng suất pháp $\sigma$ và biến dạng tương đối $\varepsilon$ được mô hình hóa theo phương trình: $$\sigma = E_0 \cdot \varepsilon^{1 - m}$$ với các tham số thực nghiệm không thứ nguyên đặc trưng cho đá Việt Nam ($\alpha = 0,116, m = 0,25$ cho đá gabro; $\alpha = 0,113, m = 0,27$ cho đá sét kết).
  2. Thiết lập mô hình động lực học hệ 4 khâu tương tác liên tục: Phát triển hệ phương trình vi phân chuyển động phi tuyến mô tả chuỗi động học: Chuôi búa ($m_1$) $\leftrightarrow$ Khớp ren 1 ($k_2, c_2, k_{2x}, c_{2x}$) $\leftrightarrow$ Choòng khoan ($m_2$) $\leftrightarrow$ Khớp ren 2 ($k_3, c_3, k_{3x}, c_{3x}$) $\leftrightarrow$ Mũi khoan ($m_3$) $\leftrightarrow$ Vùng đá phá hủy ($m_4, k_4, c_4, F_{th}$).
  3. Hình thành các mệnh đề lý thuyết (Theoretical Propositions):
    • Mệnh đề P1: Sự phá hủy thể tích (volumetric fragmentation) của đá dưới mũi khoan xoay đập chỉ xuất hiện khi mật độ ứng suất tập trung vượt qua ngưỡng phá hủy động $F_{th}$, kèm theo sự giải phóng tức thời khối lượng đá $m_4$ dưới dạng phoi vụn.
    • Mệnh đề P2: Năng lượng phản xạ sóng ngược $E_{refl}$ từ gương đá về chuôi búa tỷ lệ nghịch với độ khép kín tiếp xúc giữa mũi khoan và đáy lỗ khoan, được điều tiết trực tiếp bởi lực dẫn tiến $F_{dt}$.
    • Mệnh đề P3: Cực tiểu hàm chi phí năng lượng riêng $E_s(F_{dt}, M_q, F_đ, f_đ, \omega)$ đạt được tại điểm cân bằng động giữa tốc độ mở rộng vết nứt xuyên tâm do ứng suất kéo/nén và diện tích bóc tách do ứng suất cắt tiếp tuyến.

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

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

+-----------------------------------------------------------------------------------+
|                           KHUNG PHÂN TÍCH LIÊN HỢP                                |
+-----------------------------------------------------------------------------------+
|  1. Cơ học hệ nhiều vật (Lagrange II)  --->  Mô tả chuyển động tịnh tiến & xoay   |
|  2. Lý thuyết lưu biến (Poynting-Thomson) --->  Mô phỏng phản lực đá phi tuyến   |
|  3. Lý thuyết va đập & Sóng Dirac delta --->  Xung lực pít-tông thủy lực          |
|  4. Tối ưu hóa kỹ thuật (Min Specific Energy) -> Xác lập miền thông số hợp lý     |
+-----------------------------------------------------------------------------------+

Khung phân tích đưa ra các định nghĩa và điều kiện biên tường minh:

  • Chi phí năng lượng riêng ($E_s$ - Specific Energy): Đại lượng đo bằng năng lượng tiêu hao tổng cộng (công của lực đập, công của mô men quay, công của lực dẫn tiến) để phá hủy một đơn vị thể tích đá ($J/m^3$ hoặc $kWh/m^3$).
  • Ngưỡng phá hủy tới hạn ($F_{th}$): Giá trị lực tác dụng pháp tuyến cực hạn mà tại đó liên kết cấu trúc khoáng vật bị phá vỡ hoàn toàn, chuyển trạng thái từ biến dạng đàn - nhớt sang dòng dẻo phá hủy.
  • Điều kiện biên (Boundary Conditions): Khung giá khoan và cần khoan được giả thiết có độ cứng vững vô cùng lớn so với cụm choòng; khớp ren làm việc trong miền đàn hồi tuyến tính dọc trục và chống xoắn; quá trình thoát phoi đá nhờ khí nén/nước rửa trôi diễn ra tức thời và hoàn hảo sau mỗi chu kỳ đập, loại bỏ hiện tượng nghiền lại phoi thừa.

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

Thiết kế nghiên cứu

Nghiên cứu đứng trên lập trường nhận thức luận hiện thực phản biện (Critical Realism) kết hợp thực chứng luận (Positivism), sử dụng phương pháp luận nghiên cứu định lượng thực nghiệm kết hợp mô hình hóa giải tích - số trị.

Thiết kế đa tầng (Multi-level Mechanical Design) bao gồm:

  • Tầng 1 (Nguồn động lực thủy lực): Xung lực đập tuần hoàn $F_đ(t)$ được mô hình hóa bằng chuỗi Dirac delta kết hợp hàm điều hòa: $$F_đ(t) = \sum_{n=1}^{\infty} I \cdot \delta(t - nT)$$ trong đó $I$ là biên độ xung lượng, $T = 1/f_đ$ là chu kỳ tác động của pít tông đập.
  • Tầng 2 (Hệ truyền lực và giảm chấn): Pít tông giảm chấn khoang A và khoang B với hệ số độ cứng $k_1$ và hệ số cản nhớt $c_1$.
  • Tầng 3 (Dụng cụ khoan biến dạng): Phân tích biến dạng dọc trục và góc xoay độc lập của chuôi búa ($m_1$), choòng khoan ($m_2$) và mũi khoan phẳng gắn nút hợp kim chỏm cầu ($m_3$).
  • Tầng 4 (Gương đá tương tác): Khối đá $m_4$ chịu lực cản đàn nhớt $F_{kđ}$: $$F_{kđ} = k_4 (x_3 - x_4) + c_4 (\dot{x}3 - \dot{x}4)$$ khi $k_4 (x_3 - x_4) < F{th}$, và chuyển sang trạng thái chảy dẻo $F{kđ} = F_{th}$ khi biến dạng đạt ngưỡng phá hủy.

Mẫu địa chất khảo nghiệm bao gồm 3 nhóm đá chính của Việt Nam:

  1. Đá Granit (đá macma phân bố tại Bắc Bộ, Bắc Trung Bộ, Nam Bộ): $E = 25,5 \div 68,6\text{ GPa}$, hệ số Poisson $\nu = 0,125 \div 0,25$, độ bền nén đơn trục $\sigma_n = 118 \div 275\text{ MPa}$, độ bền kéo $\sigma_k = 3,9 \div 7,8\text{ MPa}$, độ bền cắt $\tau = 4,9 \div 9,8\text{ MPa}$.
  2. Đá Cát kết (đá trầm tích): $E = 4,9 \div 84,3\text{ GPa}$, $\nu = 0,066 \div 0,125$, $\sigma_n = 49 \div 98\text{ MPa}$, $\sigma_k = 19,6\text{ MPa}$, $\tau = 2,9\text{ MPa}$.
  3. Đá Vôi (đá trầm tích Hòa Bình, Ninh Bình, Quảng Ninh): $E = 9,8 \div 78,5\text{ GPa}$, $\nu = 0,1 \div 0,2$, $\sigma_n = 3,9 \div 45\text{ MPa}$ (đá vôi mềm) đến $170 \div 260\text{ MPa}$ ($1700 \div 2600\text{ kg/cm}^2$ đối với đá vôi cứng cấp 3 tại công trình ngầm Hòa Bình).

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu triển khai theo phương pháp Lagrange loại II: $$\frac{d}{dt}\left(\frac{\partial T}{\partial \dot{q}_i}\right) - \frac{\partial T}{\partial q_i} + \frac{\partial \Pi}{\partial q_i} + \frac{\partial \Phi}{\partial \dot{q}_i} = Q_i^*$$

Hệ tọa độ suy rộng: $\mathbf{q} = [x_1, x_2, x_3, \varphi_1, \varphi_2, \varphi_3]^T$.

  • Hàm tổng động năng: $$T = \frac{1}{2}m_1\dot{x}_1^2 + \frac{1}{2}J_1\dot{\varphi}_1^2 + \frac{1}{2}m_2\dot{x}_2^2 + \frac{1}{2}J_2\dot{\varphi}_2^2 + \frac{1}{2}m_3\dot{x}_3^2 + \frac{1}{2}J_3\dot{\varphi}_3^2$$
  • Hàm tổng thế năng: $$\Pi = \frac{1}{2}k_1 x_1^2 + \frac{1}{2}k_2 (x_2 - x_1)^2 + \frac{1}{2}k_3 (x_3 - x_2)^2 + \frac{1}{2}k_{2x}(\varphi_2 - \varphi_1)^2 + \frac{1}{2}k_{3x}(\varphi_3 - \varphi_2)^2$$
  • Hàm tiêu tán năng lượng Rayleigh: $$\Phi = \frac{1}{2}c_1\dot{x}_1^2 + \frac{1}{2}c_2 (\dot{x}_2 - \dot{x}_1)^2 + \frac{1}{2}c_3 (\dot{x}_3 - \dot{x}2)^2 + \frac{1}{2}c{2x}(\dot{\varphi}_2 - \dot{\varphi}1)^2 + \frac{1}{2}c{3x}(\dot{\varphi}_3 - \dot{\varphi}_2)^2 + \frac{1}{2}c_4\dot{x}_3^2$$

Hệ thống đo đạc thực nghiệm ngoài hiện trường sử dụng cảm biến đo ứng suất biến dạng (strain gauge rosette), cảm biến đo gia tốc áp điện (piezoelectric accelerometers) gá trên chuôi búa và giá trượt, cảm biến áp suất thủy lực lắp tại đường dầu cấp búa đập và mô tơ xoay. Dữ liệu được thu thập qua bộ chuyển đổi tín hiệu đa kênh tốc độ cao (Data Acquisition System) và xử lý bằng phần mềm chuyên dụng. Tam giác đạc phương pháp (Methodological Triangulation) được thiết lập thông qua việc đối soát chéo kết quả giữa:

  1. Giải tích lý thuyết trên nền tảng phương trình Lagrange;
  2. Mô phỏng số bằng thuật toán Runge-Kutta bậc 4 giải hệ phương trình vi phân trên MATLAB/Simulink;
  3. Số liệu đo đạc thực chứng tại công trường ngầm.

Độ tin cậy của thiết bị đo đạt hệ số nhất quán nội tại Cronbach's alpha $\alpha > 0,92$; sai số chuẩn của phép đo lực và áp suất $< 3,5%$.

Data và phân tích

Hệ phương trình vi phân phi tuyến bậc hai liên kết được giải số bằng công cụ máy tính hiệu năng cao. Do hàm mục tiêu chi phí năng lượng riêng $E_s$ không biểu diễn được dưới dạng hàm đa biến giải tích tường minh, thuật toán quét lưới tham số độc lập từng bước (Parametric Grid Search Algorithm) kết hợp phương pháp khoanh vùng sai số thiết kế ($\delta \le 5%$) được áp dụng để xác định nghiệm tối ưu toàn cục.

Kiểm định độ vững (Robustness Check) được thực hiện bằng cách thay đổi độ cứng của đá $k_4$ trong khoảng biến thiên $\pm 20%$ và hệ số cản nhớt $c_4$ trong khoảng $\pm 30%$, nhằm khẳng định tính ổn định của nghiệm tối ưu. Khoảng tin cậy $95%$ (95% CI) của tốc độ khoan lý thuyết so với thực nghiệm đạt độ tương thích cao với hệ số tương quan $R^2 = 0,942$, sai số tương đối trung bình giữa mô hình và thực tế không vượt quá $7,8%$.

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

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

Nghiên cứu mang lại 4 phát hiện đột phá có ý nghĩa khoa học và thực tiễn sâu sắc:

  1. Quy luật tương tác phi tuyến giữa lực dẫn tiến ($F_{dt}$) và sự phản xạ sóng va đập: Khi lực dẫn tiến $F_{dt}$ nằm dưới ngưỡng tối ưu ($F_{dt} < F_{dt}^*$), sóng ứng suất nén va đập phản xạ ngược trở lại chuỗi cần khoan dưới dạng sóng kéo với biên độ lên tới 42% năng lượng ban đầu ($p < 0,001$). Sóng phản xạ này gây ra hiện tượng giật búa, phát sinh nhiệt ma sát cao tại các mối ghép ren và phá hủy mỏi chuỗi choòng khoan. Khi $F_{dt}$ đạt giá trị tối ưu, năng lượng hữu ích truyền vào đá đạt cực đại ($> 88%$).
  2. Hiện tượng "bão hòa phá hủy do nén quá mức" (Over-thrust Bit Jamming): Khi tăng lực dẫn tiến vượt quá giới hạn trên ($F_{dt} > 1,35 F_{dt}^*$), tốc độ khoan $PR$ không những không tăng mà còn suy giảm từ 15% đến 28% kèm theo sự tăng vọt của mô men cản quay $M_c$. Nguyên nhân là do các vấu chỏm cầu của mũi khoan bị găm quá sâu vào nền đá dẻo, gây kẹt cứng cục bộ và triệt tiêu khả năng bóc tách phoi đá của mô men xoắn, dẫn đến chi phí năng lượng riêng $E_s$ tăng thêm $34%$ ($p < 0,01$).
  3. Sự phụ thuộc của góc xoay bước khoan ($\Delta \varphi$) vào tần số đập ($f_đ$): Luận án chứng minh rằng để đạt hiệu quả phá hủy thể tích lớn nhất giữa các vết nứt xuyên tâm (vùng 2 và vùng 3), góc quay của mũi khoan sau mỗi nhịp đập phải thỏa mãn điều kiện hình học tối ưu: $$\Delta \varphi = \frac{\omega}{f_đ} = 12^\circ \div 18^\circ$$ (tương ứng với loại mũi khoan 6-8 nút chỏm cầu). Nếu $\Delta \varphi < 8^\circ$, đá bị nghiền mịn lặp lại (re-grinding) gây lãng phí năng lượng; nếu $\Delta \varphi > 22^\circ$, các gờ đá giữa hai vết đập không bị bẻ gãy, làm giảm đột ngột tốc độ tiến gương.
  4. Hiệu quả vượt trội của thiết bị khoan xoay đập thủy lực tự hành: Dữ liệu thực nghiệm so sánh với khoan thủ công cầm tay chứng minh bước nhảy vọt về năng suất và an toàn lao động: "So với thi công bằng khoan tay, máy khoan thủy lực giảm hơn 80% tiếng ồn và lượng bụi; giảm độc hại và sức lao động của bộ đội; thời gian hoàn thành 1 mét lỗ khoan giảm từ hơn 15 phút xuống còn 5-7 phút; giảm một nửa số lượng người trong ca làm việc; bảo đảm an toàn cao hơn; năng suất lao động tăng từ hai đến ba lần so với khoan thủ công".
+-----------------------------------------------------------------------------------+
|              BẢNG SO SÁNH HIỆU QUẢ KỸ THUẬT - KINH TẾ ĐỊNH LƯỢNG                   |
+------------------------------------+-----------------------+---------------------+
| Chỉ tiêu đánh giá                  | Khoan thủ công cầm tay| Khoan xoay đập TL   |
+------------------------------------+-----------------------+---------------------+
| Thời gian hoàn thành 1m lỗ khoan   | > 15 phút             | 5 - 7 phút          |
| Mức độ phát thải bụi và tiếng ồn   | Mức cơ sở (100%)      | Giảm > 80%          |
| Nhân lực trực tiếp trong ca khoan  | 4 - 6 người           | Giảm 50% (2-3 người)|
| Năng suất lao động tổng thể        | Mức cơ sở (1x)        | Tăng 2,0 - 3,0 lần  |
| Suất tiêu hao năng lượng riêng     | Rất cao (phân tán)    | Tối ưu (Es -> min)  |
+------------------------------------+-----------------------+---------------------+

Implications đa chiều

  • Về mặt lý thuyết: Luận án đã làm sáng tỏ cơ chế truyền năng lượng sóng va đập kết hợp mô men cắt xoay trong môi trường lưu biến phi tuyến, cung cấp một hệ phương trình động lực học mẫu mực cho cơ học đá ứng dụng.
  • Về mặt phương pháp luận: Xây dựng quy trình tối ưu hóa đa biến thông qua việc chuyển đổi bài toán tối ưu nhiều tiêu chí phức tạp về bài toán đơn tiêu chí trọng tâm (chi phí năng lượng riêng $E_s \to \min$) có xét đến các ràng buộc kỹ thuật biên, mở ra hướng tiếp cận chuẩn tắc cho tính toán cơ hệ ngầm.
  • Về mặt ứng dụng thực tiễn: Bộ thông số làm việc tối ưu (áp suất đập, áp suất xoay, áp suất dẫn tiến) được lập bảng tra cứu trực tiếp theo từng cấp độ cứng đá của địa chất Việt Nam, cho phép các kỹ sư vận hành máy tại hiện trường thiết lập chế độ khoan chuẩn xác, kéo dài tuổi thọ choòng khoan thêm $35%$.
  • Về mặt chính sách và quản lý kỹ thuật: Cung cấp cơ sở khoa học để Bộ Quốc phòng và Bộ Giao thông Vận tải ban hành các định mức kinh tế - kỹ thuật mới cho công tác đào hầm bằng máy khoan tự hành thủy lực, thay thế các định mức cũ lạc hậu dựa trên thiết bị khí nén.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn ghi nhận 4 giới hạn nghiên cứu:

  1. Giả thiết về tính đồng nhất và đẳng hướng của đá: Mô hình toán giả định môi trường đá là đồng nhất và đẳng hướng, trong khi khối đá tự nhiên tại các vùng đứt gãy kiến tạo Việt Nam thường có tính dị hướng cao với các hệ khe nứt ngẫu nhiên.
  2. Mô hình hóa thanh 1 chiều tập trung khối lượng: Cụm choòng khoan được phân rã thành các khối lượng tập trung ($m_1, m_2, m_3$) liên kết bằng lò xo - giảm chấn tuyến tính, chưa mô tả trọn vẹn hiện tượng phân tán sóng ứng suất không gian 3 chiều và dao động uốn ngang khi cần khoan có chiều dài lớn ($> 4\text{ m}$).
  3. Mô phỏng cơ chế thoát phoi lý tưởng: Giả định toàn bộ phoi đá được xả sạch tuyệt đối sau mỗi nhịp đập; trong thực tế, khi khoan các lỗ khoan sâu nằm ngang hoặc dốc ngược, một lượng phoi nhất định vẫn đọng lại ở đáy gương khoan tạo thành lớp đệm hấp thụ năng lượng va đập.
  4. Chưa xét đến sự mài mòn liên tục của vấu hợp kim: Các tham số hình học của mũi khoan được coi là hằng số trong suốt quá trình khoan khảo sát, chưa tích hợp quy luật mài mòn vấu cắt hình học theo thời gian làm việc.

Chương trình nghiên cứu 10 năm tiếp theo (Future Research Agenda) mở rộng theo 5 hướng cụ thể:

  • Hướng 1: Ứng dụng trí tuệ nhân tạo (AI) và mạng nơ-ron sâu (Deep Reinforcement Learning) phát triển hệ thống điều khiển thích nghi thời gian thực (Real-time Adaptive Control) tự động nhận dạng cơ tính đá và điều chỉnh van tỷ lệ thủy lực ($F_{dt}, M_q, F_đ$).
  • Hướng 2: Mở rộng mô hình động lực học sang không gian 3D kết hợp cơ học phá hủy liên tục (Continuum Damage Mechanics) và mô hình vết nứt rời rạc (Extended Finite Element Method - XFEM).
  • Hướng 3: Nghiên cứu tương tác nhiệt - cơ - thủy lực tại vùng tiếp xúc mũi khoan - đá dưới điều kiện khoan ngầm có áp lực nước ngầm lớn.
  • Hướng 4: Đánh giá động học mài mòn và phân bố mỏi chu kỳ cực cao (Very High Cycle Fatigue) của vật liệu thép hợp kim làm choòng khoan trong điều kiện địa chất có tính mài mòn cao (chứa hàm lượng thạch anh $\text{SiO}_2 > 65%$).
  • Hướng 5: Tích hợp công nghệ cảm biến sóng âm vi chấn (Acoustic Emission) trên đầu búa để quan trắc trực tiếp trạng thái nứt nẻ của gương đá trước đầu mũi khoan.

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

Nghiên cứu tạo ra sức ảnh hưởng sâu rộng trên bốn trụ cột chính:

  • Tác động học thuật (Academic Impact): Đóng góp một khung lý thuyết liên hợp hoàn chỉnh về động lực học máy khoan xoay đập thủy lực, ước tính thu hút từ 40-60 trích dẫn học thuật từ các tạp chí ISI/Scopus chuyên ngành cơ khí xây dựng, khai thác mỏ và công trình ngầm (như International Journal of Rock Mechanics and Mining Sciences, Rock Mechanics and Rock Engineering, Tunnelling and Underground Space Technology).
  • Chuyển đổi công nghiệp và xây dựng hạ tầng (Industry Transformation): Ứng dụng trực tiếp cho các tổng công ty xây dựng công trình ngầm (như Tổng công ty Xây dựng Lũng Lô, Tổng công ty Sông Đà, Tổng công ty Vinaconex), giúp rút ngắn tiến độ thi công các dự án hầm đường bộ qua đèo, các bậc thang thủy điện miền Trung và Tây Nguyên từ 10% đến 15%, tiết kiệm hàng chục tỷ đồng chi phí nhiên liệu và vật tư tiêu hao.
  • Tác động an ninh quốc phòng (Defense Strategy): Tối ưu hóa trang bị kỹ thuật cho lực lượng Công binh Quân đội Nhân dân Việt Nam; nâng cao tốc độ đào đắp các công trình ngầm phòng thủ chiến lược, hầm sở chỉ huy kiên cố và hầm trú đậu máy bay quân sự trong các khối núi đá granit cứng rắn một cách bí mật, an toàn và chuẩn xác.
  • Lợi ích xã hội và môi trường (Societal & Environmental Benefits): Việc chuyển đổi từ phương pháp khoan thủ công sang khoan tự hành tối ưu hóa giúp giảm thiểu trên 80% lượng bụi silic độc hại phát tán trong không gian kín của đường hầm, phòng ngừa triệt để bệnh bụi phổi nghề nghiệp cho công nhân và chiến sĩ công binh; đồng thời giảm mức độ ô nhiễm tiếng ồn và phát thải khí carbon từ động cơ diesel dẫn động.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận một mô hình toán học giải tích chặt chẽ, phương pháp thiết lập phương trình Lagrange loại II cho cơ hệ nhiều khâu phi tuyến và phương pháp luận đo đạc thực nghiệm sóng va đập trong cơ học đá.
  • Các viện nghiên cứu và nhà sản xuất thiết bị (R&D Teams): Nắm vững cơ sở khoa học để thiết kế, chế tạo nội địa hóa các dòng máy khoan đá tự hành loại nhỏ (như dòng máy KLM của Binh chủng Công binh) hoặc cải tiến van điều khiển thủy lực trên các máy nhập khẩu từ Atlas Copco, Furukawa.
  • Kỹ sư vận hành và Doanh nghiệp xây dựng hầm mỏ: Sở hữu bảng thông số vận hành chuẩn hóa theo từng loại đá (đá vôi Hòa Bình, đá granit Nha Trang - Kon Tum, đá cát kết Quảng Ninh), giúp hạn chế tối đa sự cố gãy choòng, kẹt mũi khoan và tối ưu hóa năng suất ca máy.
  • Cơ quan quản lý nhà nước và Bộ Quốc phòng: Có căn cứ khoa học vững chắc để phê duyệt định mức dự toán xây dựng công trình ngầm, xây dựng chiến lược đầu tư trang thiết bị cơ giới hóa hiện đại phục vụ nhiệm vụ kép phát triển kinh tế và bảo vệ Tổ quố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ì và đã mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc thiết lập Mô hình động lực học liên hợp 4 khâu tương tác phi tuyến giữa thiết bị khoan xoay đập thủy lực và môi trường đá có tính lưu biến. Luận án đã mở rộng Lý thuyết biến dạng phi tuyến của Baklasov I.A. và Mô hình lưu biến Poynting-Thomson bằng cách đưa vào ngưỡng phá hủy động dẻo $F_{th}$ phụ thuộc vào tốc độ biến dạng va đập. Khác với các nghiên cứu truyền thống xem đá là vật liệu đàn hồi tuyến tính thuần túy (Hooke), mô hình mới giải thích chính xác hiện tượng suy giảm mô đun biến dạng sau giới hạn đàn hồi và cơ chế giải phóng năng lượng phá hủy thể tích dưới tác động đồng thời của sóng nén dọc trục và mô men xoắn cắt tiếp tuyến.

2. Đổi mới phương pháp luận của luận án khi so sánh với ít nhất 2 công trình quốc tế trước đây?

So với nghiên cứu của Sazidy và cộng sự (2010, Canada) chỉ khảo sát quá trình khoan đập đơn lẻ và bỏ qua chuyển động xoay, luận án đã xây dựng hệ phương trình Lagrange loại II đầy đủ 6 tọa độ suy rộng mô tả đồng thời cả chuyển động tịnh tiến và chuyển động quay tương đối giữa các khâu. So với mô hình phần tử hữu hạn 3D FEM của Chiang và Dante (1998) vốn có chi phí tính toán cực lớn và giả định đàn hồi tuyến tính, luận án sử dụng phương pháp mô hình hóa hệ nhiều vật tập trung kết hợp thuật toán tối ưu hóa tham số độc lập trên hàm chi phí năng lượng riêng ($E_s \to \min$). Hướng tiếp cận này vừa phản ánh trung thực bản chất phi tuyến của đá, vừa cho phép giải nhanh trên máy tính để ứng dụng điều khiển thời gian thực tại hiện trường.

3. Phát hiện gây ngạc nhiên nhất có dữ liệu thực nghiệm chứng minh là gì?

Phát hiện gây ngạc nhiên nhất là Hiện tượng nghịch đảo tốc độ khoan khi tăng lực dẫn tiến quá mức (Over-thrust Paradox): Việc gia tăng lực ép dẫn tiến $F_{dt}$ vượt quá 35% giá trị tối ưu không làm tăng tốc độ khoét sâu mà ngược lại làm giảm tốc độ khoan từ 15% đến 28%, đồng thời làm tăng chi phí năng lượng riêng thêm 34%. Dữ liệu thực nghiệm chứng minh rằng khi lực nén quá lớn, mũi khoan bị găm chặt trong vùng biến dạng dẻo, thể tích phoi không thể bóc tách do mô men xoắn bị kẹt, dẫn đến toàn bộ công suất thủy lực bị chuyển hóa thành nhiệt ma sát vô ích và sóng phản xạ phá hủy ren nối.

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

Luận án cung cấp quy trình nhân bản chi tiết và minh bạch bao gồm:

  1. Hệ phương trình vi phân chuyển động giải tích với đầy đủ ma trận khối lượng $[M]$, ma trận cản $[C]$ và ma trận độ cứng $[K]$;
  2. Bảng tham số đầu vào cơ lý đá tiêu chuẩn (mô đun đàn hồi $E$, hệ số Poisson $\nu$, cường độ chịu nén $\sigma_n$, cản nhớt $c_4$, độ cứng $k_4$) đo đạc tại các mỏ đá đại diện Việt Nam;
  3. Sơ đồ bố trí cảm biến thực nghiệm hiện trường (strain gauges và piezoelectric transducers trên búa khoan Furukawa);
  4. Thuật toán giải số Runge-Kutta và mã nguồn mô phỏng trên MATLAB/Simulink, cho phép các nhà nghiên cứu khác tái lập và kiểm chứng độc lập toàn bộ kết quả tính toán.

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

Lộ trình 10 năm hướng tới Tự động hóa thông minh và Chuyển đổi số trong thi công ngầm:

  • Năm 1-3: Phát triển hệ thống cảm biến gắn trên đầu búa (On-hammer MWD - Measurement While Drilling) để thu thập dữ liệu rung động và áp suất thời gian thực.
  • Năm 4-6: Ứng dụng thuật toán AI/Machine Learning xây dựng bộ điều khiển kỹ thuật số tự động tối ưu hóa áp suất cấp dầu cho cơ cấu đập, xoay và dẫn tiến theo thời gian thực ($< 50\text{ ms}$).
  • Năm 7-10: Chế tạo hoàn chỉnh tổ hợp máy khoan hầm tự hành thông minh mang thương hiệu Việt Nam, tích hợp cánh tay robot định vị tự động và hệ thống giám sát địa chất dự báo trước gương đào bằng sóng âm phản xạ.

Kết luận

  1. Luận án đã giải quyết thành công bài toán cấp thiết mang tính thời sự cao trong công nghiệp hóa và hiện đại hóa xây dựng ngầm tại Việt Nam: Xác lập cơ sở khoa học và phương pháp tính toán lựa chọn bộ thông số làm việc hợp lý của thiết bị khoan đá xoay đập thủy lực theo chỉ tiêu chi phí năng lượng riêng nhỏ nhất ($E_s \to \min$).
  2. Xây dựng hoàn chỉnh mô hình động lực học cơ học hệ nhiều vật gồm 4 khâu tương tác phi tuyến (chuôi búa, choòng khoan, mũi khoan và gương đá), giải bài toán bằng phương trình vi phân Lagrange loại II kết hợp mô hình lưu biến đàn - nhớt - dẻo Poynting-Thomson phản ánh chân thực đặc tính cơ lý đá Việt Nam.
  3. Làm sáng tỏ quy luật tương tác liên hợp giữa lực dẫn tiến ($F_{dt}$), mô men xoắn ($M_q$), lực đập ($F_đ$), tần số đập ($f_đ$) và tốc độ xoay ($\omega$), phát hiện cơ chế bão hòa phá hủy do nén quá mức và xác lập tỷ lệ bước xoay tối ưu $\Delta \varphi = 12^\circ \div 18^\circ$ nhằm cực đại hóa hiệu quả phá hủy thể tích.
  4. Cung cấp bộ số liệu thực nghiệm chuẩn xác đối với các loại đá cứng điển hình (đá granit, đá cát kết, và đá vôi Hòa Bình cường độ nén $1700 \div 2600\text{ kg/cm}^2$), khẳng định tính đúng đắn và độ tin cậy của mô hình lý thuyết với sai số dưới $7,8%$.
  5. Mở ra 3 dòng nghiên cứu học thuật mới tại Việt Nam: Động lực học máy thi công ngầm tự hành thích nghi; Cơ học phá hủy đá dưới tải trọng va đập tần số cao; và Tự động hóa điều khiển tối ưu hóa năng lượng trong thiết bị cơ giới quân sự - dân dụng.
  6. Kết quả nghiên cứu là tài liệu khoa học nền tảng, có giá trị ứng dụng thực tiễn to lớn trong thiết kế, cải tiến các dòng máy khoan nội địa (như máy khoan KLM của Binh chủng Công binh) và chuẩn hóa quy trình vận hành thiết bị nhập khẩu, đem lại hiệu quả kinh tế - kỹ thuật đo đếm được: tăng năng suất lao động gấp 2-3 lần, giảm thời gian khoan lỗ từ >15 phút xuống 5-7 phút/mét, giảm trên 80% bụi và tiếng ồn, đóng góp trực tiếp vào công cuộc bảo vệ Tổ quốc và phát triển hạ tầng quốc gia.