Tổng quan về luận án

Khai thác than hầm lò tại Việt Nam đang bước vào giai đoạn chuyển dịch cấu trúc sâu sắc theo định hướng của Quyết định số 403/QĐ-TTg ngày 14 tháng 3 năm 2016 của Thủ tướng Chính phủ về việc phê duyệt điều chỉnh quy hoạch phát triển ngành than Việt Nam đến năm 2020, xét triển vọng đến năm 2030. Theo đó, sản lượng than thương phẩm toàn ngành được hoạch định đạt từ 47 - 50 triệu tấn vào năm 2020 và tăng lên 55 - 57 triệu tấn vào năm 2030. Để hiện thực hóa mục tiêu sản lượng trong bối cảnh các diện khai thác lộ thiên ngày càng thu hẹp và cạn kiệt, các mỏ than hầm lò vùng Quảng Ninh (nơi chiếm gần 80% tổng trữ lượng khai thác hầm lò cả nước với 2,172 tỷ tấn tài nguyên huy động) buộc phải đẩy nhanh tốc độ xuống sâu xuống các mức âm từ -150m, -300m đến -500m (điển hình như mỏ Khe Chàm II-IV) và tăng cường áp dụng cơ giới hóa đồng bộ. Sự gia tăng độ sâu khai thác cùng tốc độ bóc tách gương than kéo theo lượng khí metan ($CH_4$) thoát ra các đường lò chuẩn bị và lò chợ tăng đột biến, gia tăng hiểm họa cháy nổ thảm khốc đe dọa sinh mạng thợ mỏ và an ninh năng lượng quốc gia.

Khoảng trống nghiên cứu (research gap) trọng tâm nằm ở chỗ: các công trình đánh giá chế độ khí trước đây tại Việt Nam (Trần Tú Ba, Lê Văn Thao) chủ yếu mang tính cục bộ, phân tán tại một số diện khai thác có độ thoát khí cao, thiếu tính bao quát toàn bồn than và chưa thiết lập được mô hình toán học dự báo chính xác độ thoát khí khi khai thác xuống sâu cho từng vỉa than antraxit có cấu trúc kiến tạo phức tạp. Luận án tiến sĩ kỹ thuật "Nghiên cứu chế độ thoát khí metan khi khai thác xuống sâu trong mỏ than hầm lò vùng Quảng Ninh" của nghiên cứu sinh Nguyễn Văn Thịnh (Trường Đại học Mỏ - Địa chất Hà Nội, 2019; chuyên ngành Khai thác mỏ, mã số 9 52 06 03; dưới sự hướng dẫn khoa học của PGS.TS Đặng Vũ Chí và TS. Lê Văn Thao) đã giải quyết triệt để khoảng trống này.

Hệ thống câu hỏi và giả thuyết nghiên cứu của luận án được xác lập cụ thể:

  1. Câu hỏi nghiên cứu 1 (RQ1): Quy luật phân bố độ chứa khí metan tự nhiên ($W_{tn}$) biến thiên theo chiều sâu địa chất và cấu trúc kiến tạo địa tầng bể than Quảng Ninh như thế nào?
  2. Câu hỏi nghiên cứu 2 (RQ2): Mức độ tương quan định lượng giữa độ thoát khí metan tương đối ($M$), độ thoát khí tuyệt đối ($V_{CH4}$) với chiều sâu khai thác ($H$), sản lượng lò chợ ($A$) và sơ đồ công nghệ khai thác chia lớp là gì?
  3. Câu hỏi nghiên cứu 3 (RQ3): Mô hình dự báo quốc tế (mô hình mỏ thực nghiệm Barbara - Ba Lan) cần được hiệu chỉnh bằng hàm số thực nghiệm nào để đạt độ chính xác tối ưu trong điều kiện địa chất mỏ antraxit Quảng Ninh?
  4. Giả thuyết nghiên cứu 1 (H1): Độ chứa khí metan tự nhiên và độ thoát khí metan tỷ lệ thuận phi tuyến với độ sâu khai thác, chịu sự chi phối trực tiếp của đặc tính khoáng sàng dạng kín hoặc hở.
  5. Giả thuyết nghiên cứu 2 (H2): Trong cùng điều kiện địa chất - kỹ thuật mỏ, độ thoát khí metan lò chợ lớp vách luôn lớn hơn lò chợ lớp trụ khi áp dụng công nghệ khai thác chia lớp nghiêng.
  6. Giả thuyết nghiên cứu 3 (H3): Tồn tại một hàm quy chuẩn tuyến tính hiệu chỉnh mô hình Barbara có hệ số xác định $R^2 > 0,95$, phản ánh chính xác chế độ thoát khí thực tế tại Quảng Ninh.

Khung lý thuyết của luận án tích hợp lý thuyết hấp phụ phân tử khí trong cấu trúc xốp than đá của Ianôscôi, lý thuyết phân vùng miền khí biến chất sâu theo mô hình địa chất động lực công trình, và phương pháp cân bằng vật chất động lực học mỏ thực nghiệm Barbara (Ba Lan). Đóng góp đột phá của luận án là xác lập hàm hồi quy hiệu chỉnh độ thoát khí $y = 0,8651x + 0,00246$ với hệ số tương quan thực nghiệm vượt trội $R^2 = 0,9896$, giải quyết bài toán dự báo cho quy mô 30 mỏ than hầm lò thuộc ba vùng than trọng điểm Uông Bí - Đông Triều, Hòn Gai, Cẩm Phả trong khung thời gian dữ liệu kéo dài từ 2007 đến 2018.

Literature Review và Positioning

Tổng quan y văn quốc tế ghi nhận ba trường phái học thuật chủ đạo trong nghiên cứu và dự báo độ thoát khí metan mỏ hầm lò:

  • Trường phái Tây Âu (Shulza, Wintera và Stuffkena): Xây dựng công thức giải tích xác định độ thoát khí tương đối dựa trên độ chứa khí tự nhiên của vỉa khai thác ($W_{tn}$) và hệ số ảnh hưởng thoát khí từ các vỉa lân cận ($\eta_{lc}$) thông qua chiều dày lớp kẹp ($m_{lc}$) theo phương trình: $$M = W_{tn} + \sum \frac{m_{lc}}{m_{kt}} \cdot \eta_{lc} \cdot W_{lc} - W_{lc}^t$$
  • Trường phái Đông Âu và Liên Xô cũ: Tiếp cận theo phương pháp thống kê kinh nghiệm dựa trên thông số "Bậc giàu khí" ($L$), biểu thị mức gia tăng chiều sâu địa chất (m) để độ thoát khí tương đối tăng thêm $1\text{ m}^3/\text{tấn.ngày-đêm}$: $$M_p = M_{p1} + \frac{H - H_1}{L} \quad \text{với} \quad L = \frac{H_2 - H_1}{M_{p2} - M_{p1}}$$
  • Trường phái Ba Lan (Mỏ thực nghiệm Barbara): Thiết lập phương trình cân bằng động học phân rã nguồn khí lò chợ thành 4 thành phần riêng biệt gồm khí thoát từ than khai thác phá vỡ ($V_{kt}$), khí thoát từ mặt lộ gương than ($V_g$), khí từ vách trên ($V_{tr}$) và khí từ trụ dưới ($V_{di}$): $$V_{CH4} = V_{kt} + V_g + V_{tr} + V_{di} = \frac{L_{lc} \cdot m_{kh} \cdot \gamma_v \cdot W_v \cdot p}{1440} + V_g + V_{tr} + V_{di}$$

Trong y văn tồn tại cuộc tranh luận học thuật sâu sắc giữa hai quan điểm đối nghịch:

  1. Quan điểm 1 (Địa chất tĩnh học): Cho rằng độ chứa khí tự nhiên trong vỉa than ($W_{tn}$) là hằng số quyết định hoàn toàn lượng khí thoát ra đường lò, không phụ thuộc đáng kể vào các thao tác công nghệ khai thác cục bộ.
  2. Quan điểm 2 (Động lực học công nghệ mỏ): Khẳng định chu kỳ khai thác, tiến độ nổ mìn gương than, tốc độ hạ trần phá hỏa và cấu trúc thông gió (thông gió chữ U nghịch so với thông gió thuận) mới là yếu tố tạo nên các đỉnh xung kích phát thải metan đột biến.

Luận án định vị nghiên cứu tại điểm giao thoa giữa đặc thù thạch học than antraxit biến chất cao vùng Quảng Ninh (chiếm 96,19% trữ lượng, cấu trúc nứt nẻ vi kiến tạo dày đặc) và các mô hình tính toán động lực học. So sánh đối chuẩn với hai bể than lớn trên thế giới cho thấy sự khác biệt bản chất:

  • So với bể than Thượng Silesia (Ba Lan): Than bitum có độ sâu >1000m, vỉa thoải đều, khoáng sàng kín kinh điển không chứa khí ở tầng nông <600m nhưng tăng vọt đạt 11 - 12 $\text{m}^3/\text{tấn khối cháy}$ ở độ sâu >700m. Ngược lại, bể than Quảng Ninh là vỉa nghiêng đến dốc đứng, phân cắt bởi hệ thống đứt gãy kiến tạo phức tạp (tại mỏ Hà Lầm phát hiện 129 đứt gãy, Mạo Khê phát hiện 88 đứt gãy trong quá trình khai thác), tạo ra chế độ thoát khí không đồng nhất theo không gian.
  • So với bể than Đôn-bát và Karaganda (Nga/Kazakhstan): Cường độ thoát khí mặt gương than antraxit Quảng Ninh suy giảm theo quy luật hàm mũ nhanh hơn sau khi bóc tách gương than, tuy nhiên áp lực thoát khí tích tụ tại khoảng không gian đã khai thác (phá hỏa) lại duy trì ở mức cao và kéo dài phạm vi ảnh hưởng tới 100 - 250m dọc theo đường lò thông gió gió thải.

Đó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 nhiệt động lực học hấp phụ khí của Ianôscôi và Langmuir bằng cách chứng minh cấu trúc phân bố ba trạng thái của metan trong than antraxit Việt Nam: dạng tự do trong lỗ rỗng (chiếm 5 - 22% ở áp suất 50 bar; 36% ở 100 bar; 65% ở 800 bar), dạng hấp phụ vật lý bề mặt hạt vi mô than, và dạng hấp phụ xâm nhập "dung dịch rắn" liên kết hóa học thuận nghịch.

Mô hình lý thuyết của luận án thiết lập 3 mệnh đề khoa học (propositions) cốt lõi:

  • Mệnh đề P1: Tại khoáng sàng dạng kín, độ chứa khí tự nhiên biến thiên theo chiều sâu tuân theo quy luật hàm mũ tăng trưởng $y = a \cdot e^{bx}$, trong khi ở khoáng sàng hở có lớp đất phủ chắn khí, độ chứa khí đạt giá trị cực đại cục bộ dưới lớp phủ trước khi giảm nhẹ và tái tăng trưởng ở đới biến chất sâu.
  • Mệnh đề P2: Tồn tại sự phân hóa rõ rệt về chế độ thoát khí theo phân lớp công nghệ: lò chợ khai thác lớp vách giải phóng năng lượng liên kết khí của toàn bộ vỉa than mạnh hơn từ 25 - 40% so với lò chợ khai thác lớp trụ do tác động của sự dỡ tải áp lực mỏ và nứt nẻ sơ cấp của đá vách.
  • Mệnh đề P3: Nồng độ metan trong luồng gió thải biến thiên có chu kỳ đồng pha với biểu đồ sản xuất: đạt đỉnh cao nhất khi tiến hành khoan nổ mìn và phá hỏa (tăng từ 0,2% đến 0,33% nồng độ thể tích tại mỏ Brzeszcze và Mạo Khê) và giảm xuống mức thấp nhất trong giai đoạn xúc tải than tĩnh.

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

Khung phân tích của luận án tích hợp ba trụ cột lý thuyết: (1) Lý thuyết cân bằng địa cơ học khối đá vỉa than dỡ tải, (2) Lý thuyết nhiệt động lực học khuếch tán khí xốp, và (3) Lý thuyết động học dòng khí thông gió mỏ hầm lò.

Cách tiếp cận phân tích mới nằm ở việc lượng hóa "Hệ số hiệu chỉnh địa phương hóa" ($k_{VN} = 0,8651$) kết hợp hằng số dịch chuyển tự do ($c = 0,00246$) để chuyển đổi mô hình tính toán Barbara từ than mỡ năng lượng (Bituminous coal) của Ba Lan sang than antraxit biến chất cao của Việt Nam. Điều kiện biên xác lập cho mô hình bao gồm: áp dụng cho các vỉa than có góc dốc từ thoải đến dốc ($0^\circ - 90^\circ$), chiều sâu khai thác từ mức lộ thiên đến -500m, áp suất khí vỉa than dưới 100 bar, và áp dụng cho các hệ thống khai thác cột dài theo phương hoặc chia lớp nghiêng điều khiển đá vách bằng phá hỏa toàn phần.

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

Thiết kế nghiên cứu

Luận án tuân thủ triết lý nghiên cứu thực chứng (positivism) kết hợp trường phái hiện thực phê phán (critical realism), dựa trên việc đo đạc thực nghiệm khách quan các thông số hóa lý, thống kê mỏ quy mô lớn và xử lý hồi quy toán học chính xác.

Nghiên cứu áp dụng thiết kế đa tầng (multi-level research design) gồm ba cấp độ:

  • Cấp độ 1 (Macro-level): Phân tích địa chất toàn vùng bể than Quảng Ninh (diện tích 300 $\text{km}^2$, chiều dài 130 km từ Phả Lại đến Vạn Hoa) phân chia thành 3 tiểu vùng: Uông Bí - Đông Triều, Hòn Gai, Cẩm Phả.
  • Cấp độ 2 (Meso-level): Khảo sát 3 mỏ than hầm lò đại diện điển hình cho ba vùng cấu trúc địa chất gồm mỏ Mạo Khê (khu Tây), mỏ Hà Lầm (khu Trung tâm), mỏ Khe Chàm 1 (khu Đông).
  • Cấp độ 3 (Micro-level): Phân tích chi tiết 15 lò chợ cụ thể tại các tầng khai thác khác nhau (ví dụ: lò chợ 6ĐMR mức +30/+67, -25/+30, -80/-25, -150/-80 mỏ Hà Lầm; lò chợ vỉa 10 khu III mức +40/+24, +20/-18, -20/-46, -46/-65, -70/-100 mỏ Mạo Khê; các lò chợ 13.1, 13.2, 13.2-4 mức -100/-125 mỏ Khe Chàm 1).

Quy trình nghiên cứu rigorous

Quy trình thu thập dữ liệu áp dụng nguyên tắc tam giác giác đạc (triangulation) nghiêm ngặt:

  • Tam giác đạc dữ liệu: Kết hợp chuỗi số liệu quan trắc khí liên tục từ hệ thống đo tự động MMI-01 và đầu đo khí cố định với kết quả đo nhanh bằng thiết bị đo khí cầm tay chuyên dụng tại gương lò và luồng gió thải.
  • Tam giác đạc phương pháp: Đối chiếu phương pháp cân bằng thông gió mỏ với phương pháp giải hấp khí mẫu lõi khoan thăm dò địa chất trong phòng thí nghiệm.
  • Tính giá trị và độ tin cậy: Toàn bộ quy trình đo kiểm tuân thủ nghiêm ngặt Quy chuẩn kỹ thuật Quốc gia về an toàn trong khai thác than hầm lò QCVN 01:2011/BCT ban hành kèm Thông tư số 03/2011/TT-BCT. Dụng cụ đo được hiệu chuẩn định kỳ tại Trung tâm An toàn Mỏ (Viện Khoa học Công nghệ Mỏ - Vinacomin) với độ sai số cho phép $\pm 0,1% \text{ CH}_4$.

Data và phân tích

Phân tích dữ liệu sử dụng thuật toán hồi quy bình phương nhỏ nhất (Ordinary Least Squares - OLS) và giải tích hồi quy phi tuyến trên nền tảng phần mềm thống kê chuyên ngành để tối ưu hóa hàm tương quan giữa độ sâu ($H$), sản lượng ($A$) và độ thoát khí ($M, V_{CH4}$).

Các phân tích kiểm tra độ vững chắc (robustness checks) được thực hiện bằng cách đối chiếu chéo kết quả tính toán lý thuyết với chuỗi số liệu vận hành thực tế tại 15 diện lò chợ trong giai đoạn 2007 - 2015. Kết quả kiểm định hàm dự báo cho thấy hệ số xác định $R^2$ đạt giá trị 0,9896, sai số tuyệt đối trung bình (MAE) nhỏ hơn 5%, khẳng định tính tương thích cao của mô hình đối với điều kiện địa chất than antraxit Quảng Ninh.

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

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

  1. Quy luật tương quan đồng biến giữa độ sâu và độ thoát khí: Xác lập bằng chứng định lượng chứng minh độ chứa khí và độ thoát khí tăng đơn điệu theo chiều sâu khai thác. Cụ thể tại mỏ Hà Lầm, khi khai thác từ mức +30/+67 xuống mức -150/-80, độ thoát khí metan tuyệt đối tại lò chợ 6ĐMR tăng từ $0,65\text{ m}^3/\text{phút}$ lên $2,45\text{ m}^3/\text{phút}$; tại mỏ Mạo Khê vỉa 10 khu III, độ thoát khí tăng liên tục từ mức +40/+24 xuống mức -70/-100.
  2. Quy luật phân hóa độ thoát khí theo công nghệ chia lớp: Tại các vỉa than dày khai thác chia lớp nghiêng, độ thoát khí metan tại lò chợ khai thác lớp vách cao hơn từ 28% đến 35% so với lò chợ khai thác lớp trụ dưới trong cùng một khối trữ lượng, do tác động nứt nẻ và giải phóng khí hấp phụ sớm của lớp than nóc.
  3. Mối quan hệ phi tuyến giữa sản lượng khai thác và độ thoát khí: Phát hiện độ thoát khí tuyệt đối ($m^3/\text{phút}$) tỷ lệ thuận tuyến tính với sản lượng lò chợ ($A$), nhưng độ thoát khí tương đối ($m^3/\text{tấn}$) lại có xu hướng giảm nhẹ khi sản lượng tăng cao đột biến trong thời gian ngắn do hiệu ứng trễ khuếch tán khí từ than nguyên khối.
  4. Xây dựng hàm hồi quy chuẩn hóa mỏ Việt Nam: Thiết lập thành công phương trình dự báo thực nghiệm: $$y = 0,8651x + 0,00246 \quad (R^2 = 0,9896)$$ Trong đó $x$ là giá trị độ thoát khí tính toán theo mô hình mỏ thực nghiệm Barbara, và $y$ là giá trị độ thoát khí đo đạc thực tế tại các lò chợ mỏ than hầm lò vùng Quảng Ninh.
  5. Thành lập bản đồ phân vùng khí metan theo không gian và chiều sâu: Xác định ranh giới 4 miền khí biến chất cho toàn bồn than Quảng Ninh, chứng minh đới biến chất sâu bắt đầu xuất hiện đậm nét từ mức âm -100m trở xuống với hàm lượng $CH_4$ tự do trong hỗn hợp khí vỉa vượt trên 80%.

Implications đa chiều

  • Về mặt lý thuyết: Luận án đã hoàn thiện cơ sở lý luận về động học phát thải khí metan trong than antraxit biến chất cao, cung cấp luận cứ khoa học để bác bỏ giả định áp dụng nguyên bản các công thức tính toán than mỡ Tây Âu vào điều kiện Việt Nam mà không qua hiệu chỉnh.
  • Về mặt phương pháp luận: Thiết lập khung quy trình 4 bước chuẩn hóa gồm: đo đạc mẫu vỉa $\rightarrow$ tính toán theo mô hình phân rã nguồn Barbara $\rightarrow$ hiệu chỉnh bằng hàm hồi quy OLS $\rightarrow$ lập bản đồ đẳng trị dự báo cho mức khai thác kế tiếp.
  • Về mặt ứng dụng thực tiễn: Luận án đề xuất giải pháp kỹ thuật đồng bộ kiểm soát cháy nổ khí metan:
    • Giải pháp chủ động: Áp dụng công nghệ khoan tháo khí metan trước và song hành trong quá trình khai thác bằng hệ thống máy khoan chuyên dụng WD-02EA, WD-02H, WDP-1C, WDP-2A, WDH-1 kết hợp máy bơm ép rửa áp lực cao WT-30-2PB tại các khám khoan chuyên trách dọc lò thông gió.
    • Giải pháp quản lý luồng khí: Xây dựng hai sơ đồ công nghệ xử lý khí: dẫn khí metan qua đường ống chuyên dụng xả trực tiếp lên mặt đất qua tháp xả an toàn hoặc pha loãng cưỡng bức và xả vào luồng gió thải của trạm quạt gió chính mỏ.
  • Về mặt chính sách: Cung cấp cơ sở khoa học để Bộ Công Thương và Tập đoàn Vinacomin điều chỉnh, bổ sung các tiêu chuẩn an toàn trong QCVN 01:2011/BCT, bắt buộc phân cấp mỏ theo mức độ nguy hiểm nổ khí khi khai thác xuống dưới mức -150m.

Limitations và Future Research

Luận án thẳng thắn chỉ ra 4 giới hạn nghiên cứu cụ thể:

  1. Giới hạn mẫu địa chất: Dữ liệu thực nghiệm tập trung chủ yếu vào than antraxit vùng Quảng Ninh, chưa bao quát các mỏ than nâu Đồng bằng Sông Hồng hoặc than mỡ vùng Thái Nguyên, Nông Sơn.
  2. Độ sâu quan trắc thực tế: Mức sâu tối đa đo đạc thực tế trong giai đoạn nghiên cứu đạt mức -150m (tại thời điểm 2018), các mức sâu từ -300m đến -500m chủ yếu dựa trên ngoại suy hàm toán học và mô hình hóa.
  3. Tác động của nước dưới đất: Chưa tích hợp sâu mô hình thủy văn động lực vào phương trình thoát khí trong điều kiện áp lực cột nước ngầm lớn khi xuống sâu.
  4. Biến dạng vi cấu trúc: Chưa sử dụng kính hiển vi điện tử quét (SEM) để định lượng kích thước lỗ xốp nano của hạt than trong quá trình giải hấp khí.

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

  • Hướng 1: Nghiên cứu ứng dụng công nghệ thu hồi khí metan từ vỉa than (Coalbed Methane - CBM) để phát điện thương mại, giảm phát thải khí nhà kính theo cam kết Net Zero.
  • Hướng 2: Xây dựng thuật toán trí tuệ nhân tạo (Machine Learning/Deep Learning) kết hợp mạng cảm biến IoT để dự báo đột biến thoát khí metan theo thời gian thực (Real-time dynamic forecasting).
  • Hướng 3: Khảo sát chế độ thoát khí trong điều kiện áp dụng cơ giới hóa đồng bộ hạng nặng (giàn chống tự hành kết hợp máy khấu than Combai công suất lớn).
  • Hướng 4: Nghiên cứu hiện tượng bùng nổ khí và than đá đồng thời (Coal and Gas Outburst) tại các đới đứt gãy kiến tạo sâu dưới mức -300m.

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

  • Tác động học thuật: Luận án là công trình nghiên cứu hệ thống và toàn diện bậc nhất về chế độ khí mỏ than antraxit tại Việt Nam, mở ra tiềm năng trích dẫn lớn trong các tạp chí thuộc danh mục ISI/Scopus chuyên ngành Mining Engineering và Geotechnical Safety.
  • Chuyển đổi công nghiệp mỏ: Được ứng dụng trực tiếp tại các công ty than nòng cốt thuộc Tập đoàn Công nghiệp Than - Khoáng sản Việt Nam (TKV) như Công ty Than Mạo Khê, Than Hà Lầm, Than Khe Chàm, Than Nam Mẫu, Than Vàng Danh. Giúp tối ưu hóa lưu lượng gió thông gió mỏ, giảm chi phí điện năng vận hành quạt gió chính từ 8 - 12% trong khi vẫn đảm bảo nồng độ metan an toàn $< 0,75 - 1,0%$ theo quy chuẩn.
  • Tác động an toàn xã hội: Ngăn ngừa nguy cơ xảy ra các thảm họa nổ khí thảm khốc (tương tự như các vụ tai nạn lịch sử tại Mạo Khê năm 1999 làm 19 người chết, Khe Chàm năm 2008 làm 11 người chết), bảo vệ an toàn tính mạng cho hàng vạn thợ mỏ hầm lò.
  • Giá trị quốc tế: Cung cấp bộ dữ liệu thực nghiệm quý giá về hành vi khí metan trong than antraxit có cấu trúc kiến tạo nứt nẻ mạnh, đóng góp vào cơ sở dữ liệu chung của Hiệp hội Mỏ Quốc tế về an toàn mỏ hầm lò sâu.

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

  • Nghiên cứu sinh và học giả chuyên ngành Mỏ: Tiếp cận phương pháp luận chuẩn hóa về đo đạc, mô hình hóa và hiệu chỉnh sai số thực nghiệm trong nghiên cứu thông gió và an toàn mỏ khí.
  • Kỹ sư trưởng và Cán bộ quản lý an toàn mỏ: Ứng dụng trực tiếp công thức dự báo $y = 0,8651x + 0,00246$ để tính toán nhu cầu lưu lượng gió sạch, lựa chọn công suất trạm quạt thông gió và lập hộ chiếu khoan tháo khí chi tiết cho từng lò chợ.
  • Lãnh đạo Tập đoàn Công nghiệp Than - Khoáng sản Việt Nam (TKV): Có căn cứ khoa học vững chắc để phê duyệt dự án đầu tư mở rộng mỏ xuống sâu, quy hoạch hệ thống giếng đứng và trạm khử khí metan tập trung.
  • Cơ quan quản lý Nhà nước (Cục Kỹ thuật An toàn và Môi trường Công nghiệp - Bộ Công Thương): Sử dụng kết quả nghiên cứu làm căn cứ sửa đổi, hoàn thiện Quy chuẩn kỹ thuật Quốc gia QCVN 01:2011/BCT và các văn bản quy phạm pháp luật về quản lý khí mỏ.

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ì? Trả lời: Đóng góp độc đáo nhất là mở rộng mô hình phân rã nguồn khí động học Barbara của Ba Lan bằng cách xác lập hàm hồi quy tuyến tính thực nghiệm $y = 0,8651x + 0,00246$ với $R^2 = 0,9896$, giải quyết thành công sự khác biệt căn bản về đặc tính thạch học và biến chất giữa than mỡ thoải đồng nhất Tây Âu với than antraxit dốc nứt nẻ mạnh của bể than Quảng Ninh.

  2. Điểm đổi mới phương pháp luận so với các công trình trước đây là gì? Trả lời: So với phương pháp giải tích thuần túy của Shulza - Wintera (chỉ tính tĩnh) và phương pháp thống kê Bậc giàu khí của Liên Xô (đòi hỏi chuỗi số liệu quá khứ cồng kềnh), luận án đã kết hợp phương pháp phân rã nguồn 4 thành phần với kỹ thuật quan trắc tự động tập trung MMI-01 và xử lý số liệu bằng giải tích hồi quy bình phương nhỏ nhất đa biến, cho phép dự báo chủ động cho các tầng khai thác chưa từng mở vỉa.

  3. Phát hiện nào mang tính bất ngờ nhất và có bằng chứng dữ liệu hỗ trợ ra sao? Trả lời: Phát hiện bất ngờ nhất là độ thoát khí tương đối ($m^3/\text{tấn}$) có xu hướng giảm nhẹ khi sản lượng lò chợ tăng đột biến trong chu kỳ khai thác ngắn hạn, mặc dù độ thoát khí tuyệt đối ($m^3/\text{phút}$) tăng mạnh. Dữ liệu thực nghiệm tại lò chợ vỉa XXXII mỏ than Brzeszcze và lò chợ vỉa 10 mỏ Mạo Khê chứng minh rằng thời gian ủ nhiệt giải hấp khí của khối than lớn hơn tốc độ bóc tách cơ học của máy khấu, dẫn đến hiện tượng trễ phát tán khí ra luồng gió.

  4. Luận án có cung cấp quy trình tái lập (replication protocol) đầy đủ không? Trả lời: Luận án cung cấp chi tiết toàn bộ quy trình tái lập thực nghiệm bao gồm: thông số kỹ thuật hệ thống máy khoan tháo khí WD-02EA, WD-02H; áp lực máy bơm rửa WT-30-2PB; sơ đồ bố trí lỗ khoan thăm dò tại khám khoan; các bước tính toán hệ số quy đổi năng lượng TOE ($1\text{ tấn FO} = 1\text{ TOE}$, $1\text{ tấn than antraxit} = 0,7\text{ TOE}$); và thuật toán hồi quy OLS tương thích với mọi bộ dữ liệu mỏ hầm lò.

  5. Chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào? Trả lời: Định hướng trọng tâm vào việc phát triển hệ thống điều khiển thông gió mỏ thông minh (Smart Mine Ventilation) ứng dụng AI/IoT; nghiên cứu công nghệ ép nứt vỉa than bằng thủy lực (Hydraulic Fracturing) để tháo kiệt khí metan trước khi khai thác; và chuyển hóa khí metan mỏ hầm lò (VAM/CMM) thành năng lượng sạch phục vụ mục tiêu phát triển mỏ xanh bền vững.

Kết luận

  1. Xác lập quy luật biến thiên độ chứa khí metan: Khẳng định quy luật gia tăng độ chứa khí tự nhiên theo chiều sâu khai thác tại bể than Quảng Ninh, làm rõ sự khác biệt giữa khoáng sàng kín và khoáng sàng hở có lớp đất phủ không thấm khí.
  2. Khám phá cơ chế phân hóa thoát khí theo công nghệ khai thác: Chứng minh định lượng rằng trong hệ thống khai thác chia lớp nghiêng, lò chợ lớp vách có cường độ thoát khí metan lớn hơn lò chợ lớp trụ từ 28 - 35%.
  3. Thành lập phương trình dự báo hiệu chỉnh thực nghiệm chuẩn xác cao: Xây dựng hàm hồi quy tương quan $y = 0,8651x + 0,00246$ với hệ số tin cậy $R^2 = 0,9896$, cho phép dự báo chính xác độ thoát khí tuyệt đối và tương đối cho các mức khai thác xuống sâu đến -500m.
  4. Thành lập bản đồ phân vùng khí metan toàn bồn than: Phân chia chi tiết 4 miền khí biến chất và lập bản đồ dự báo nguy cơ khí metan cho 30 mỏ than hầm lò vùng Quảng Ninh.
  5. Đề xuất gói giải pháp công nghệ phòng ngừa cháy nổ đồng bộ: Hoàn thiện giải pháp khoan tháo khí chủ động bằng dàn thiết bị WD kết hợp bơm áp lực cao WT-30-2PB và giải pháp xả khí an toàn, đáp ứng triệt để các yêu cầu nghiêm ngặt của Quy chuẩn Quốc gia QCVN 01:2011/BCT.
  6. Mở ra 3 hướng nghiên cứu học thuật đột phá mới: Đặt nền móng cho các nghiên cứu chuyên sâu về khai thác khí CBM thương mại, dự báo động lực học thời gian thực bằng trí tuệ nhân tạo, và kiểm soát hiện tượng bùng nổ than - khí ở độ sâu siêu giới hạn dưới -500m.