Tổng quan về luận án

Nghiên cứu địa chất và đánh giá độ tin cậy trong công tác thăm dò than đóng vai trò sống còn đối với hiệu quả kinh tế và an toàn kỹ thuật của các dự án khai thác mỏ hầm lò xuống sâu. Luận án tiến sĩ kỹ thuật địa chất của Nghiên cứu sinh Phạm Tuấn Anh (Trường Đại học Mỏ - Địa chất Hà Nội, 2023, mã số 9520501) với đề tài "Nghiên cứu ứng dụng một số mô hình toán - địa chất đánh giá mức độ tin cậy của công tác thăm dò than khu mỏ Khe Chàm, thành phố Cẩm Phả, tỉnh Quảng Ninh" đã giải quyết triệt để bài toán chuyển đổi phương pháp luận từ hình học mỏ truyền thống sang mô hình toán - địa chất và địa thống kê hiện đại.

                    ┌────────────────────────────────────────────────────────┐
                    │     BỂ THAN QUẢNG NINH: ĐỊA HÀO HÒN GAI (T3n-r hg2)    │
                    └──────────────────────────┬─────────────────────────────┘
                                               │
                    ┌──────────────────────────▼─────────────────────────────┐
                    │  KHU MỎ KHE CHÀM (16.2 km² - Giới hạn bởi F.A-A & F.BH)│
                    └──────────────────────────┬─────────────────────────────┘
                                               │
                 ┌─────────────────────────────┼─────────────────────────────┐
                 │                             │                             │
    ┌────────────▼────────────┐  ┌─────────────▼────────────┐  ┌─────────────▼────────────┐
    │     KHỐI TÂY NAM        │  │      KHỐI TRUNG TÂM      │  │      KHỐI ĐÔNG BẮC       │
    │  (Tây Nam đứt gãy F.E)  │  │   (Giữa F.E và đứt gãy   │  │ (Đông Bắc đứt gãy F.L)   │
    │  Nhóm mỏ thăm dò: III   │  │    F.L) - Nhóm mỏ: II    │  │  Nhóm mỏ thăm dò: III    │
    └─────────────────────────┘  └──────────────────────────┘  └──────────────────────────┘

Trong lịch sử thăm dò bể than Quảng Ninh, mạng lưới công trình chủ yếu được thiết kế theo các quy định định tính chung ngoài đồng bằng sông Hồng (Quyết định số 25/2007/QĐ-BTNMT). Tại khu mỏ Khe Chàm (diện tích 16,2 km²), dù mạng lưới thăm dò tỷ mỷ $500 \times (250 \div 300)\text{ m}$ về mặt lý thuyết đã đủ cơ sở cấp trữ lượng 122 và 121 để lập dự án đầu tư khai thác hầm lò (như mỏ Khe Chàm III công suất 2,5 triệu tấn/năm khai thác đến mức -500m; mỏ Khe Chàm II-IV công suất 3,5 triệu tấn/năm khai thác đến mức -700m), nhưng thực tế thi công luôn phát sinh sai lệch nghiêm trọng. "Khi triển khai dự án đầu tư khai thác, tài liệu địa chất vẫn cần phải bổ sung, điều chỉnh mặc dù các khối trữ lượng đã được thăm dò ở cấp cao (cấp 121, 122)." Sự xuất hiện của các đứt gãy bậc cao có biên độ dịch chuyển $5 \div 30\text{ m}$, hiện tượng uốn nếp cục bộ, vát mỏng vỉa đột ngột và ô cửa sổ không than đã buộc các đơn vị khai thác phải chia nhỏ lò chợ cơ giới hóa đồng bộ (Combie) thành các lò chợ thủ công, gây lãng phí đầu tư và gián đoạn sản xuất.

Research gap cốt lõi được xác định là: Thiếu vắng một khung phương pháp luận định lượng đa biến hóa để khảo sát đồng thời ba phương diện biến tính của vỉa than (biến hóa ngẫu nhiên, biến hóa có quy luật/xu thế và tương quan không gian khu vực hóa), dẫn đến việc phân loại nhóm mỏ thăm dò áp đặt cào bằng và lựa chọn thuật toán tính trữ lượng không phản ánh đúng cấu trúc hình học thực tế.

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

  1. Câu hỏi 1 (Q1): Mức độ phức tạp cấu trúc địa chất và hình thái vỉa than tại khu mỏ Khe Chàm có thực sự đồng nhất để áp dụng một nhóm mỏ thăm dò duy nhất (loại II) hay tồn tại các phân khối cấu trúc bậc cao mang đặc tính biến dị riêng biệt?
    • Giả thuyết 1 (H1): Khu mỏ Khe Chàm phân hóa thành các khối đồng nhất tương đối bậc cao thuộc các nhóm mỏ thăm dò khác nhau (Nhóm II và Nhóm III) dưới tác động của hệ thống đứt gãy kiến tạo chính ($F.E, F.L$).
  2. Câu hỏi 2 (Q2): Những thông số địa chất công nghiệp nào đóng vai trò kiểm soát chính xác độ tin cậy của mô hình trữ lượng khi khai thác xuống sâu?
    • Giả thuyết 2 (H2): Chiều dày vỉa (tự nhiên và riêng than) cùng góc dốc vỉa là các thông số chi phối lớn nhất, chịu sự chi phối mạnh của tính bất đẳng hướng không gian.
  3. Câu hỏi 3 (Q3): Việc ứng dụng phương pháp địa thống kê (Kriging) trên mô hình khối 3 chiều (3D Block Model) có giải quyết được các sai số cố hữu của phương pháp hình học mỏ truyền thống (Sêcăng/Cosêcăng) hay không?
    • Giả thuyết 3 (H3): Mô hình hóa Kriging thông dụng kết hợp phân tích hàm cấu trúc không gian giúp tối ưu hóa mạng lưới thăm dò và nâng cao độ tin cậy trữ lượng thương mại, giảm thiểu rủi ro thiết kế đường lò.

Nghiên cứu khảo sát tập dữ liệu khống chế trên 11 vỉa than công nghiệp thuộc phân hệ tầng Hòn Gai giữa ($T_3n-r hg_2$) chiếm trên 87% tổng tài nguyên mỏ, trong đó 4 vỉa than đại diện ($V_{14-5}, V_{14-1}, V_{13-1}, V_{10}$) được lập mô hình chi tiết từ 1.306 công trình khoan và khai đào sâu xuống mức -700m.

Literature Review và Positioning

Các công trình nghiên cứu cấu trúc bể than Quảng Ninh khởi đầu từ những khảo sát địa tầng của Zeiller R. (1903), Bourret R. (1922), Deprat J. (1915), cho đến công trình phân chia điệp Hòn Gai của Cheremnykh V.M. (1964) và bản đồ địa chất tỷ lệ 1:500.000 của Dovjikov A.E. (1965). Về mặt kiến tạo, Trần Văn Trị (1977, 1990) đã xác lập mô hình bể than Quảng Ninh gồm hai địa hào Bảo Đài và Hòn Gai phát triển trên móng Caledoni - Cathaysia và trũng chồng Mesozoi An Châu, chia cắt bởi các đứt gãy sâu phương á kinh tuyến và vĩ tuyến.

Trong lĩnh vực phương pháp luận thăm dò than và đánh giá trữ lượng, lịch sử tồn tại hai trường phái đối lập:

  • Trường phái hình học mỏ truyền thống (Classical Mining Geometry): Dựa trên lý thuyết của Sobolevsky P.K. (1932), Popoff C.C. (1966), Đặng Trần Bảng (1983), Trần Văn Miến (1995), Nguyễn Sỹ Quý (1992). Trường phái này sử dụng phương pháp đa giác, đường đẳng trị và phương pháp chiếu Sêcăng ($Q = \sum S_i \cdot \sec \alpha_i \cdot m_i \cdot D$) hoặc Cosêcăng. Giả định cơ bản là các thông số địa chất thay đổi tuyến tính hoặc phẳng giữa các công trình thăm dò. Điểm yếu cốt tử là bỏ qua tương quan không gian thực tế, dẫn đến sai số diện tích ($S$), sai số chiều dày ($m$), và góc dốc ($\alpha$) tăng cao tại các vùng uốn nếp phức tạp.
  • Trường phái địa thống kê và giải tích toán - địa chất (Spatial Geostatistics & Mathematical Geology): Phát triển từ Thuyết biến ngẫu nhiên khu vực hóa của Georges Matheron (1963, 1971), kết hợp phân tích bề mặt xu thế (Trend Surface Analysis) của Agterberg F.P. (1974) và Davis J.C. (1986). Trường phái này phân rã biến địa chất thành các thành phần: xu thế khu vực vĩ mô, cấu trúc tương quan ngẫu nhiên có định hướng và sai số ngẫu nhiên cục bộ (hiệu ứng khối).
                      TIẾN TRÌNH PHÁT TRIỂN PHƯƠNG PHÁP LUẬN ĐỊNH LƯỢNG
┌──────────────────────────────────────┐       ┌──────────────────────────────────────┐
│  HÌNH HỌC MỎ TRUYỀN THỐNG (1932-1990) │       │  ĐỊA THỐNG KÊ & MÔ HÌNH TOÁN (2023)  │
│  - Sobolevsky, Popoff, Đặng Trần Bảng │       │  - Phạm Tuấn Anh, Matheron, Surpac   │
├──────────────────────────────────────┤       ├──────────────────────────────────────┤
│ • Chiếu phẳng Secăng/Cosecăng        │ ────► │ • Hàm cấu trúc không gian Variogram  │
│ • Giả định biến đổi tuyến tính       │       │ • Phân tích Trend xu thế bậc 2       │
│ • Không đo lường bất đẳng hướng      │       │ • Mô hình khối 3D Ordinary Kriging   │
│ • Nguy cơ sai lệch lớn khi xuống sâu │       │ • Lượng hóa sai số và bán kính khống │
└──────────────────────────────────────┘       └──────────────────────────────────────┘

Trên thế giới, việc tích hợp địa thống kê vào thăm dò than đã được chuẩn hóa tại nhiều bể than lớn:

  • Nghiên cứu tại Bể than Karoo (Nam Phi) của Snyman C.P. & Botha G.A. (1993) và Armstrong et al. đã chứng minh việc sử dụng Ordinary Kriging trên vỉa than Permian phân tán giúp giảm thiểu phương sai ước lượng tới 34% so với phương pháp cổ điển.
  • Nghiên cứu tại Bể than Appalachian và Powder River (Hoa Kỳ) theo chuẩn USGS Circular 891 (Olea R.A. & Luppens J.A., 2011; Hohn M.E., 1999) đã sử dụng thực nghiệm Variogram đa hướng để tái xác lập khoảng cách mạng lưới khoan thăm dò phân cấp trữ lượng tin cậy (Measured/Indicated Resource) thay cho các cự ly cố định chủ quan.

Tại Việt Nam, trước công trình của Phạm Tuấn Anh (2023), các nghiên cứu của Nguyễn Phương (1994, 2017) và Đoàn Văn Kiển (2012) mới dừng ở việc phân loại nhóm mỏ vĩ mô hoặc ứng dụng cục bộ địa thống kê mà chưa khảo sát toàn diện bản chất biến hóa 3 chiều cho một tổ hợp vỉa than cụ thể có cấu tạo bất đối xứng, biến chất cao (antraxit) tại đới bờ Nam địa hào Hòn Gai. Luận án đã định vị chính xác khoảng trống này, đưa ra lời giải tích hợp toàn diện giữa cấu trúc kiến tạo, hàm xu thế đa thức bậc 2 và mô hình hàm cấu trúc không gian ba chiều.

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

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

  1. Mở rộng Thuyết biến ngẫu nhiên khu vực hóa (Theory of Regionalized Variables - Georges Matheron) vào đối tượng vỉa than nếp uốn biến dạng mạnh: Luận án chứng minh rằng thông số chiều dày vỉa than ($m$) và độ cao trụ vỉa ($Z$) không tồn tại như một biến ngẫu nhiên thuần túy, mà là hàm khu vực hóa $Z(x)$ cấu thành từ hai đại lượng: thành phần có quy luật thể hiện qua hàm xu thế $m(x)$ và thành phần biến động ngẫu nhiên dừng cục bộ $Y(x)$: $$Z(x) = m(x) + Y(x)$$
  2. Xác lập lý thuyết phân khối đồng nhất tương đối bậc cao trong địa chất mỏ: Bằng chứng thực nghiệm từ luận án đã bẻ gãy quan điểm phân hạng nhóm mỏ đồng nhất toàn mỏ. Cấu trúc khu mỏ Khe Chàm bị chia cắt bởi hai đứt gãy kiến tạo chính là đứt gãy $F.E$ và đứt gãy $F.L$, tạo ra 3 khối đồng nhất tương đối bậc cao có cơ chế biến dạng và mức độ biến thiên thông số hoàn toàn dị biệt:
    • Khối Trung tâm (nằm giữa đứt gãy $F.E$ và $F.L$): Cấu trúc nếp uốn thoải, vỉa duy trì liên tục $\rightarrow$ Xếp nhóm mỏ thăm dò loại II.
    • Khối Tây Nam (phía Tây Nam đứt gãy $F.E$) và Khối Đông Bắc (phía Đông Bắc đứt gãy $F.L$): Cấu trúc uốn nếp phức tạp, mật độ đứt gãy dày đặc, vỉa than bị phình thóp, vát nhọn $\rightarrow$ Bắt buộc xếp nhóm mỏ thăm dò loại III.
  3. Mô hình định lượng hóa các chỉ tiêu hình học - cấu tạo vỉa: Hệ thống hóa và công thức hóa các chỉ số kiểm soát hình thái: hệ số gián đoạn vỉa ($K_d$), hệ số phức tạp cấu tạo ($K_{pt}$), môđun chu tuyến ($\mu$), chỉ tiêu hình dạng vỉa ($\Phi$), cùng chỉ tiêu phân tán góc dốc ($\sigma_\alpha$).
           ┌─────────────────────────────────────────────────────────────┐
           │        KHUNG PHÂN TÍCH BA PHƯƠNG DIỆN BIẾN HÓA ĐỊA CHẤT     │
           └──────────────────────────────┬──────────────────────────────┘
                                          │
        ┌─────────────────────────────────┼─────────────────────────────────┐
        │                                 │                                 │
┌───────▼────────┐             ┌──────────▼──────────┐           ┌──────────▼──────────┐
│ PHƯƠNG DIỆN 1  │             │   PHƯƠNG DIỆN 2     │           │   PHƯƠNG DIỆN 3     │
│  Biến hóa ngẫu │             │   Biến hóa có quy   │           │   Tương quan cấu    │
│  nhiên cục bộ  │             │   luật (Vĩ mô)      │           │   trúc không gian   │
├────────────────┤             ├─────────────────────┤           ├─────────────────────┤
│• Toán thống kê │             │• Phân tích Trend    │           │• Hàm thực nghiệm    │
│  cổ điển       │             │  mặt bậc 2          │           │  $\gamma(h)$        │
│• $V, S, \bar{X}$│            │• Tách dị thường     │           │• Mô hình Spherical  │
│• Chỉ tiêu $K_d$│             │  dư $\Delta Z$      │           │• Bán kính ảnh hưởng │
│  $\mu, \Phi$   │             │• Xu thế uốn nếp     │           │  $a$, Nugget $C_0$  │
└────────────────┘             └─────────────────────┘           └─────────────────────┘

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

Khung phân tích của luận án tích hợp ba lý thuyết nền tảng: Lý thuyết kiến tạo phân khối (Structural Geology), Lý thuyết hồi quy bề mặt xu thế (Trend Surface Analysis), và Lý thuyết địa thống kê (Geostatistics).

Quy trình phân tích 3 phương diện biến hóa được thiết lập:

  • Phương diện 1 (Biến hóa ngẫu nhiên): Sử dụng các đặc trưng thống kê cổ điển (kỳ vọng toán $\bar{X}$, phương sai $S^2$, hệ số biến thiên $V%$) để đánh giá định lượng độ phân tán của chiều dày tự nhiên ($m_t$), chiều dày riêng than ($m_r$) và độ tro khô ($A^k$).
  • Phương diện 2 (Biến hóa có quy luật): Thiết lập mô hình hồi quy Trend bậc 2 trên tọa độ không gian $(X, Y)$ nhằm bóc tách bề mặt cấu trúc khu vực tổng thể của trụ vỉa và nhận diện các dị thường dư vị trí ($\Delta Z$): $$Z(X, Y) = b_0 + b_1X + b_2Y + b_3X^2 + b_4XY + b_5Y^2 + \varepsilon$$
  • Phương diện 3 (Tương quan cấu trúc không gian): Khảo sát hàm cấu trúc bán phương sai thực nghiệm $\gamma(h)$ theo các hướng đường phương ($0^\circ$) và hướng dốc ($90^\circ$), quy nạp về các mô hình lý thuyết (Spherical Model, Exponential Model) để xác lập bán kính ảnh hưởng cực đại ($a_{max}$), cực tiểu ($a_{min}$) và tỷ số bất đẳng hướng ($k = a_{max} / a_{min}$).

Điều kiện biên: Khung phân tích được xác lập tối ưu cho các vỉa than antraxit trầm tích lục địa nếp uốn vũng vịnh thuộc hệ tầng Hòn Gai ($T_3n-r hg$), có chiều dày từ trung bình đến rất dày ($1,5\text{ m} \div 38\text{ m}$), góc dốc trung bình $25^\circ \div 27^\circ$, và chiều sâu khảo sát từ mặt địa hình xuống mức -700m.

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

┌────────────────────────────────────────────────────────────────────────────────────────┐
│                        QUY TRÌNH THỰC HIỆN PHƯƠNG PHÁP NGHIÊN CỨU                      │
└───────────────────────────────────────────┬────────────────────────────────────────────┘
                                            │
┌───────────────────────────────────────────▼────────────────────────────────────────────┐
│ 1. THIẾT LẬP CƠ SỞ DỮ LIỆU ĐỊA CHẤT ĐỒNG BỘ TRÊN GEOVIA SURPAC VÀ ACCESS               │
│    • 1.306 công trình khoan/khai đào (V14-5: 350; V14-1: 223; V13-1: 489; V10: 276)    │
└───────────────────────────────────────────┬────────────────────────────────────────────┘
                                            │
┌───────────────────────────────────────────▼────────────────────────────────────────────┐
│ 2. XỬ LÝ TOÁN THỐNG KÊ & PHÂN TÍCH TREND BẬC 2 TRÊN GOLDEN SOFTWARE SURFER             │
│    • Đánh giá V%, Kd, μ, Φ; Bóc tách xu thế khu vực m(x) và dị thường dư ΔZ             │
└───────────────────────────────────────────┬────────────────────────────────────────────┘
                                            │
┌───────────────────────────────────────────▼────────────────────────────────────────────┐
│ 3. KHẢO SÁT HÀM CẤU TRÚC KHÔNG GIAN BÁN PHƯƠNG SAI VARIOGRAM γ(h)                      │
│    • Hướng 0° (đường phương), Hướng 90° (hướng dốc), góc quét 22.5°                     │
│    • Xác lập hiệu ứng khối C0, thềm C+C0, bán kính ảnh hưởng không gian (a)            │
└───────────────────────────────────────────┬────────────────────────────────────────────┘
                                            │
┌───────────────────────────────────────────▼────────────────────────────────────────────┐
│ 4. MÔ HÌNH HÓA KHỐI 3D & NỘI SUY ORDINARY KRIGING                                      │
│    • Khởi tạo Vi khối (Block Model): 50x50x0.8n và 20x20x0.8n                          │
│    • Đánh giá trữ lượng tài nguyên & Đối sánh trực tiếp với phương pháp Sêcăng         │
└────────────────────────────────────────────────────────────────────────────────────────┘

Thiết kế nghiên cứu

Nghiên cứu đứng trên lập trường nhận thức luận thực chứng phản biện (Critical Realism) kết hợp triết lý thực chứng định lượng (Positivism). Thiết kế nghiên cứu đa cấp độ (Multi-level design) kết hợp giữa phương pháp mô hình hóa cấu trúc không gian 3 chiều và thuật toán phân tích dữ liệu thực tế:

  • Cấp độ vĩ mô: Cấu trúc địa hào Hòn Gai và các hệ thống đứt gãy ranh giới bậc III ($F.A-A$ và đứt gãy Bắc Huy $F.BH$).
  • Cấp độ trung mô: Các phân khối đồng nhất bậc cao (Khối Tây Nam, Khối Trung tâm, Khối Đông Bắc).
  • Cấp độ vi mô: Các vi khối tính toán (Block model cell: $50\times50\times0,8n\text{ m}$ và $20\times20\times0,8n\text{ m}$) của từng vỉa than riêng biệt.

Quy trình nghiên cứu rigorous

Quy trình xử lý số liệu tuân thủ các tiêu chuẩn kiểm soát nghiêm ngặt:

  1. Chuẩn hóa cơ sở dữ liệu: Số liệu từ 1.306 công trình khoan và tuyến hào qua các thời kỳ (1963–2015) được chuẩn hóa tọa độ không gian VN-2000 (kinh tuyến trục $105^\circ$, múi chiếu $6^\circ$), số hóa toàn bộ chiều dày thực, độ dốc, độ tro ($A^k$), tỷ lệ đới phá hủy ($P_p$) trên hệ thống cơ sở dữ liệu quan hệ Access và Surpac.
  2. Triangulation (Tam giác đạc kiểm chứng): Kết hợp kiểm chứng chéo giữa tài liệu địa vật lý lỗ khoan (đo carota gamma, mật độ), vết lộ địa chất thực tế mặt đất, số liệu quan trắc trong các lò ngầm mỏ Khe Chàm I, Khe Chàm III và moong lộ thiên Cao Sơn.
  3. Kiểm định độ tin cậy: Khảo sát hệ số tin cậy mẫu thông qua giá trị tương quan đa biến, đánh giá độ lệch chuẩn $\sigma$ và loại trừ các giá trị dị thường (outliers) trước khi đưa vào tính toán hàm cấu trúc.

Data và phân tích

Tập dữ liệu phân tích tập trung vào 4 vỉa than chủ lực thuộc tập giữa:

  • Vỉa 14-5: 350 công trình (318 lỗ khoan sâu, 32 vết lộ/hào đào); chiều dày tự nhiên biến thiên $0,24 \div 38,87\text{ m}$ (trung bình $6,06\text{ m}$); chiều dày riêng than $0,24 \div 27,37\text{ m}$ (trung bình $5,31\text{ m}$); số lớp đá kẹp $0 \div 9$ lớp (trung bình dày $0,59\text{ m}$); độ dốc $3^\circ \div 60^\circ$ (trung bình $27^\circ$).
  • Vỉa 14-1: 223 công trình; chiều dày tự nhiên $0,00 \div 8,67\text{ m}$ (trung bình $1,74\text{ m}$); chiều dày riêng than $0,00 \div 5,98\text{ m}$ (trung bình $1,54\text{ m}$); độ dốc $3^\circ \div 70^\circ$ (trung bình $25^\circ$).
  • Vỉa 13-1: 489 công trình; chiều dày tự nhiên $0,19 \div 28,90\text{ m}$ (trung bình $4,24\text{ m}$); chiều dày riêng than $0,19 \div 21,23\text{ m}$ (trung bình $3,47\text{ m}$); đá kẹp $0 \div 12$ lớp (dày trung bình $0,66\text{ m}$); độ dốc $3^\circ \div 74^\circ$ (trung bình $26^\circ$).
  • Vỉa 10: 276 công trình; chiều dày tự nhiên $0,28 \div 24,17\text{ m}$ (trung bình $4,88\text{ m}$); chiều dày riêng than $0,28 \div 20,00\text{ m}$ (trung bình $4,32\text{ m}$); đá kẹp $0 \div 8$ lớp (dày trung bình $0,52\text{ m}$); độ dốc $5^\circ \div 67^\circ$ (trung bình $26^\circ$).

Các công cụ phần mềm chuyên dụng được vận hành:

  • GEOVIA Surpac (Dassault Systèmes): Xây dựng CSDL địa chất mỏ, tính toán mô hình thực nghiệm Variogram, dựng bề mặt trụ/vách (Surface Modeling), khởi tạo mô hình khối 3D (3D Block Model) và thực hiện nội suy Ordinary Kriging.
  • Golden Software Surfer: Xây dựng ma trận mạng lưới ô vuông (Grid Data), thực hiện các phép toán đại số bản đồ (Grid Math), phân tích hàm xu thế Trend bậc 2 và thành lập bình đồ đẳng trị trụ vỉa, đẳng dày, đẳng độ dốc.

Hàm cấu trúc bán phương sai lý thuyết được thiết lập theo mô hình cầu (Spherical Model): $$\gamma(h) = \begin{cases} C_0 + C \left[ \dfrac{3}{2}\left(\dfrac{h}{a}\right) - \dfrac{1}{2}\left(\dfrac{h}{a}\right)^3 \right] & \text{khi } h \le a \ C_0 + C & \text{khi } h > a \end{cases}$$ Trong đó: $C_0$ là hiệu ứng khối (nugget effect), $C$ là phương sai cấu trúc, $C_0 + C$ là thềm (sill), $a$ là bán kính ảnh hưởng (range of influence), $h$ là khoảng cách dịch chuyển vector giữa các điểm khảo sát.

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

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

      PHÂN BỐ TRƯỢNG BÁN KÍNH ẢNH HƯỞNG (a) VÀ TỶ SỐ BẤT ĐẲNG HƯỚNG
┌────────────┬──────────────────────────────────────┬───────────────────────┐
│  VỈA THAN  │ BÁN KÍNH KHỐNG CHẾ HÀM CẤU TRÚC (a)  │ BẤT ĐẲNG HƯỚNG (k)    │
├────────────┼──────────────────────────────────────┼───────────────────────┤
│  Vỉa 14-5  │ • Hướng phương (0°): a = 180 - 220 m │ k = 1.38 - 1.47       │
│            │ • Hướng dốc (90°):   a = 130 - 150 m │                       │
├────────────┼──────────────────────────────────────┼───────────────────────┤
│  Vỉa 13-1  │ • Hướng phương (0°): a = 160 - 200 m │ k = 1.33 - 1.54       │
│            │ • Hướng dốc (90°):   a = 120 - 130 m │                       │
├────────────┼──────────────────────────────────────┼───────────────────────┤
│  Vỉa 10    │ • Hướng phương (0°): a = 200 - 240 m │ k = 1.43 - 1.60       │
│            │ • Hướng dốc (90°):   a = 140 - 150 m │                       │
└────────────┴──────────────────────────────────────┴───────────────────────┘
  1. Phát hiện tính phân dị kiến tạo và xếp hạng lại nhóm mỏ thăm dò: Trái với các báo cáo thăm dò giai đoạn trước xếp toàn bộ khu mỏ Khe Chàm vào nhóm mỏ loại II, kết quả tính toán định lượng các chỉ tiêu hình học mỏ ($\mu = 1,12 \div 1,45; \Phi = 0,65 \div 0,82; K_d = 0,05 \div 0,28$) đã khẳng định tính bất đồng nhất. Khối Trung tâm đủ tiêu chí nhóm II, nhưng Khối Tây Nam và Đông Bắc bắt buộc phải thuộc nhóm III do mức độ phức tạp kiến tạo cao và tính gián đoạn vỉa lớn.
  2. Xác định định lượng tính bất đẳng hướng không gian của vỉa than: Khảo sát mô hình Variogram trên các vỉa $V_{14-5}, V_{13-1}, V_{10}$ cho thấy rõ tính dị hướng không gian:
    • Theo hướng đường phương ($0^\circ$): Bán kính ảnh hưởng thực nghiệm đạt cực đại $a_{max} = 180 \div 240\text{ m}$.
    • Theo hướng dốc vỉa ($90^\circ$): Bán kính ảnh hưởng thực nghiệm bị thu hẹp đáng kể $a_{min} = 120 \div 150\text{ m}$.
    • Tỷ số bất đẳng hướng ($k = a_{max} / a_{min}$): Dao động từ $1,33 \div 1,60$. Phát hiện này bác bỏ hoàn toàn việc áp dụng mạng lưới thăm dò dạng ô vuông đẳng hướng trước đây.
  3. Vạch rõ bản chất sai số trữ lượng của phương pháp Sêcăng truyền thống: Đối sánh trữ lượng tính bằng phương pháp Kriging (mô hình khối $50\times50\times0,8n\text{ m}$ và $20\times20\times0,8n\text{ m}$) với phương pháp Sêcăng trong báo cáo năm 2015 cho thấy độ lệch tổng thể dao động từ $-3,45%$ đến $+7,12%$, cá biệt tại từng khối cục bộ bị đứt gãy phân cắt, độ sai lệch lên tới trên $15 \div 22%$. Nguyên nhân là phương pháp Sêcăng đã làm mượt giả tạo bề mặt trụ vỉa và bỏ qua hiện tượng vát nhọn cục bộ.

Bảng đối sánh trữ lượng tài nguyên than giữa các phương pháp

Số liệu đối sánh chi tiết trữ lượng than (nghìn tấn) giữa mô hình khối Kriging ($50\times50\times0,8n$) và phương pháp Sêcăng truyền thống tại các vỉa than đại diện khu mỏ Khe Chàm:

Tên vỉa than Cấp trữ lượng / TN Trữ lượng theo Secăng (BC 2015) (nghìn tấn) Trữ lượng mô hình Kriging (Surpac) (nghìn tấn) Chênh lệch tuyệt đối ($\Delta Q$) (nghìn tấn) Tỷ lệ chênh lệch ($%$)
Vỉa 14-5 Cấp 122 + 121 22.450,6 21.675,8 -774,8 -3,45%
Vỉa 14-1 Cấp 122 + 121 4.820,3 5.163,5 +343,2 +7,12%
Vỉa 13-1 Cấp 122 + 121 38.910,2 37.860,0 -1.050,2 -2,70%
Vỉa 10 Cấp 122 + 121 28.340,5 29.110,4 +769,9 +2,72%
Tổng cộng Toàn bộ 4 vỉa 94.521,6 93.809,7 -711,9 -0,75%

Bảng mạng lưới thăm dò đề xuất tối ưu theo nhóm mỏ và cấp trữ lượng

Phân khối cấu trúc Nhóm mỏ xác lập Cấp trữ lượng / TN Mạng lưới hiện trạng trước đây (m) Mạng lưới đề xuất tối ưu theo Luận án (m) Hiệu quả kỹ thuật - kinh tế
Khối Trung tâm (giữa đứt gãy F.E và F.L) Nhóm II Cấp 121 (Trữ lượng tin cậy cao) $250 \times 150$ $200 \times 100 \div 120$ Định vị chính xác phình thóp, đảm bảo lò chợ Combie
Cấp 122 (Trữ lượng tin cậy) $500 \times (250 \div 300)$ $(300 \div 400) \times (150 \div 200)$ Giảm rủi ro sai số diện tích, khớp nối biến thiên dốc
Khối Tây Nam & Đông Bắc (ngoài F.E, F.L) Nhóm III Cấp 121 Không áp dụng (chưa khống chế) $(100 \div 150) \times (60 \div 80)$ Khống chế biên độ đứt gãy $5 \div 30\text{ m}$
Cấp 122 $500 \times 250$ (Không đạt độ tin cậy) $(150 \div 200) \times (80 \div 100)$ Ngăn chặn hiện tượng trượt vỉa và ô cửa sổ không than

Implications đa chiều

  • Ý nghĩa lý luận: Chuẩn hóa phương pháp luận đánh giá độ tin cậy thăm dò khoáng sản rắn uốn nếp tại Việt Nam, kết hợp hài hòa giữa toán thống kê, hàm hồi quy Trend và hàm bán phương sai Variogram.
  • Ý nghĩa phương pháp luận: Thiết lập quy trình công nghệ tích hợp phần mềm Surpac - Surfer - Access thành chuỗi xử lý khép kín từ số liệu thô lỗ khoan đến mô hình khối mỏ 3D.
  • Ý nghĩa thực tiễn sản xuất: Cung cấp thông số hình học mỏ chính xác phục vụ thiết kế các tuyến lò dọc vỉa, giếng nghiêng và định vị gương khai thác cơ giới hóa tại Công ty Than Khe Chàm và Công ty Than Hạ Long, loại bỏ tình trạng phải dừng lò chuyển diện đột xuất.
  • Hàm ý chính sách quản lý tài nguyên: Làm cơ sở khoa học để Hội đồng Đánh giá Trữ lượng Khoáng sản Quốc gia xem xét, cập nhật quy chế tính trữ lượng bằng phương pháp mô hình 3D Kriging thay thế dần cho phương pháp mặt cắt Sêcăng cổ điển.

Limitations và Future Research

  1. Hạn chế về dữ liệu thăm dò xuống sâu dưới mức -500m: Dù diện tích nghiên cứu trải rộng 16,2 km², số lượng lỗ khoan xuyên sâu khống chế tập vỉa dưới (từ vỉa 1 đến vỉa 8 ở độ sâu dưới -500m đến -900m) còn thưa thớt (ví dụ lỗ khoan LK2525 gặp vỉa 1 ở độ sâu cos -933,84m). Mật độ công trình sâu chưa đủ dày để lập Variogram chi tiết cho tập vỉa dưới.
  2. Hiệu ứng khối ($C_0$) cao tại các đới đứt gãy cục bộ: Tại các vùng giao cắt kiến tạo phức tạp, biến sai ngẫu nhiên đo đạc chiều dày và độ tro tăng cao do hiện tượng vò nhàu, trượt tầng kiến tạo sinh biến dạng dẻo, làm giảm hệ số tương quan của mô hình cầu.
  3. Giới hạn biên công nghệ phần mềm: Luận án áp dụng mô hình vi khối trực giao ($50\times50\times0,8n\text{ m}$), chưa áp dụng mô hình lưới uốn dẻo theo lớp (Stratigraphic Gridding / Dynamic Block Modeling) để phản ánh hoàn hảo các vỉa than bị uốn nếp đứng góc dốc $>60^\circ \div 75^\circ$.

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Triển khai công nghệ mô phỏng địa thống kê điều kiện (Conditional Simulation - Sequential Gaussian Simulation) để lượng hóa rủi ro phân bố không gian và dải độ tin cậy $95%$ cho từng khối trữ lượng.
  • Mở rộng khảo sát địa thống kê 3 chiều cho toàn bộ các vỉa than sâu thuộc dải Cẩm Phả - Hòn Gai dưới mức -300m đến -1.200m.
  • Ứng dụng trí tuệ nhân tạo (Machine Learning/Deep Learning) kết hợp mạng Geostatistics để tự động nhận diện và phân tích đồng danh vỉa than qua các đứt gãy kiến tạo phức tạp.

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

               LAN TỎA GIÁ TRỊ VÀ TÁC ĐỘNG ĐA LĨNH VỰC CỦA CÔNG TRÌNH
┌───────────────────────────────────────┬───────────────────────────────────────┐
│     HỌC THUẬT & ĐÀO TẠO ĐỊA CHẤT      │     CÔNG NGHIỆP KHAI THÁC THAN        │
├───────────────────────────────────────┼───────────────────────────────────────┤
│ • Tiên phong chuẩn hóa Geostatistics  │ • Tiết kiệm hàng tỷ đồng khoan thăm dò│
│ • Giáo trình tham khảo sau đại học    │ • Triệt tiêu sự cố vỡ diện lò chợ     │
│ • Kết nối chuẩn quốc tế JORC/NI 43-101│ • Tối ưu thu hồi tài nguyên >87%      │
└───────────────────────────────────────┴───────────────────────────────────────┘
                                   │
                                   ▼
┌───────────────────────────────────────────────────────────────────────────────┐
│                    QUẢN TRỊ TÀI NGUYÊN & HOẠCH ĐỊNH CHÍNH SÁCH                │
├───────────────────────────────────────────────────────────────────────────────┤
│ • Cung cấp cơ sở chuyển đổi số toàn diện cho Tập đoàn TKV                     │
│ • Định hình khung quy chuẩn phê duyệt trữ lượng 3D cấp Quốc gia (HĐTLQG)      │
└───────────────────────────────────────────────────────────────────────────────┘
  • Tác động học thuật: Đặt nền móng vững chắc cho hướng nghiên cứu Toán - Tin Địa chất ứng dụng tại Việt Nam. Công trình là tài liệu tham khảo chuẩn mực trong đào tạo Tiến sĩ, Thạc sĩ ngành Kỹ thuật Địa chất, Địa kỹ thuật tại Trường Đại học Mỏ - Địa chất và các viện nghiên cứu chuyên ngành.
  • Tác động tới ngành công nghiệp mỏ: Giúp Tập đoàn Công nghiệp Than - Khoáng sản Việt Nam (TKV) và Tổng công ty Đông Bắc tối ưu hóa chi phí hàng triệu USD cho các chiến dịch khoan thăm dò bổ sung không cần thiết; đồng thời giảm thiểu tối đa các sự cố kỹ thuật đổ lò, bẹp giá chống khi gặp đứt gãy bất ngờ trong quá trình cơ giới hóa khai thác than hầm lò.
  • Tác động chính sách và xã hội: Đảm bảo an ninh năng lượng quốc gia thông qua việc đánh giá chính xác, minh bạch tài nguyên than sạch; nâng cao hệ số thu hồi than antraxit quý hiếm, giảm thiểu tổn thất tài nguyên và hạn chế sụt lún bề mặt ảnh hưởng đến môi trường sinh thái thành phố Cẩm Phả.

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

  • Nghiên cứu sinh và Giới học thuật Địa chất: Tiếp cận được phương pháp luận kết hợp chặt chẽ giữa toán địa chất, phân tích Trend và mô hình bán phương sai Variogram; giải quyết bài toán đồng danh và khống chế cấu trúc vỉa phức tạp.
  • Kỹ sư Trắc địa - Mỏ và Địa chất công trình: Nắm vững quy trình dựng mô hình khối 3D trên GEOVIA Surpac và Golden Software Surfer, nâng cao năng lực thiết kế mỏ hiện đại.
  • Lãnh đạo doanh nghiệp mỏ (TKV, Tổng công ty Đông Bắc): Có công cụ định lượng đáng tin cậy để ra quyết định đầu tư hàng nghìn tỷ đồng vào các dự án giếng đứng, giếng nghiêng xuống sâu một cách an toàn và hiệu quả.
  • Cơ quan quản lý Nhà nước (Bộ Tài nguyên và Môi trường, Hội đồng Đánh giá Trữ lượng): Có cơ sở khoa học xác thực để thẩm định các báo cáo kết quả thăm dò khoáng sản, phê duyệt trữ lượng minh bạch theo chuẩn mực hiện đại.

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 Thuyết biến ngẫu nhiên khu vực hóa (Theory of Regionalized Variables) của Georges Matheron vào điều kiện nếp uốn biến dạng kiến tạo phức tạp của bể than Hòn Gai. Luận án đã tích hợp thành công hàm xu thế Trend bậc 2 với hàm bán phương sai Variogram để bóc tách triệt để 3 phương diện biến hóa (ngẫu nhiên cục bộ, quy luật vĩ mô và tương quan không gian), chứng minh sự tồn tại của 3 khối đồng nhất tương đối bậc cao có nhóm mỏ thăm dò khác nhau trong cùng một mỏ.

2. Đổi mới phương pháp luận của luận án có gì vượt trội so với các nghiên cứu trước đây?

So với các nghiên cứu của Đặng Trần Bảng (1983) và Trần Văn Miến (1995) chỉ dựa trên hình học phẳng Sêcăng truyền thống, hay các nghiên cứu của Nguyễn Phương (1994) và Đoàn Văn Kiển (2012) mang tính phân loại định tính, luận án đã tiên phong thiết lập mô hình khối 3D trực giao kết hợp giải thuật nội suy Ordinary Kriging có định hướng bất đẳng hướng không gian thực nghiệm ($a_{max} = 180 \div 240\text{ m}$, $a_{min} = 120 \div 150\text{ m}$), giúp lượng hóa sai số trữ lượng đến từng vi khối khai thác.

3. Phát hiện bất ngờ và phản trực giác nhất từ dữ liệu thực nghiệm là gì?

Phát hiện phản trực giác nhất là: Mặc dù khu mỏ Khe Chàm đã được thăm dò tỷ mỷ và phê duyệt toàn bộ thuộc nhóm mỏ loại II với mạng lưới khoan $500 \times (250 \div 300)\text{ m}$, thực tế cấu trúc kiến tạo đã phân cắt mỏ thành các khối dị biệt. Hai khối rìa (Tây Nam và Đông Bắc) có hệ số gián đoạn vỉa cao và phân tán góc dốc mạnh, bắt buộc phải xếp vào nhóm mỏ loại III với mạng lưới thăm dò khống chế dày hơn gấp 1,5–2 lần mới đảm bảo độ tin cậy trữ lượng thương mại.

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

Hoàn toàn có. Luận án đã chuẩn hóa cấu trúc cơ sở dữ liệu quan hệ Access, quy trình 4 bước dựng mô hình biến sai trong Surpac, các tham số hồi quy ma trận trong Surfer, và thuật toán tính trữ lượng vi khối. Toàn bộ quy trình này có thể tái lập và áp dụng nguyên vẹn cho các khu mỏ than liền kề như Khe Tam, Cao Sơn, Bàng Nâu, Hà Ráng hoặc các mỏ than tương đồng trên toàn quốc.

5. Định hướng lộ trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?

Lộ trình 10 năm tới bao gồm:

  • Mở rộng mạng lưới mô hình hóa địa thống kê xuống các tầng than sâu dưới mức -1.000m toàn bể than Quảng Ninh.
  • Tích hợp mô phỏng địa thống kê đa biến liên tục (Sequential Gaussian Co-simulation) kết hợp giữa chiều dày vỉa, độ tro và độ bền cơ lý đá vách trụ.
  • Xây dựng hệ thống tự động hóa cập nhật dữ liệu khai thác thực tế (Real-time Model Updating) phục vụ mỏ than thông minh (Smart Mining 4.0).

Kết luận

  1. Xác lập thành công hệ phương pháp luận toán - địa chất hiện đại: Kết hợp đồng bộ toán thống kê cổ điển, hàm xu thế Trend bậc 2 và hàm cấu trúc bán phương sai Variogram để đánh giá toàn diện 3 phương diện biến hóa của vỉa than.
  2. Đột phá trong nhận thức cấu trúc kiến tạo mỏ Khe Chàm: Luận chứng khoa học chia khu mỏ 16,2 km² thành 3 khối đồng nhất tương đối bậc cao, xếp Khối Trung tâm vào nhóm mỏ thăm dò II và hai khối rìa Tây Nam - Đông Bắc vào nhóm mỏ thăm dò III.
  3. Định lượng hóa tính bất đẳng hướng không gian: Khám phá tỷ số dị hướng $k = 1,33 \div 1,60$ giữa hướng phương ($a_{max} = 180 \div 240\text{ m}$) và hướng dốc ($a_{min} = 120 \div 150\text{ m}$), làm tiền đề tái thiết kế mạng lưới công trình thăm dò hợp lý.
  4. Chứng minh độ tin cậy vượt trội của mô hình khối 3D Kriging: Loại bỏ các sai lệch nội suy của phương pháp hình học Sêcăng truyền thống, kiểm soát độ lệch trữ lượng ở mức an toàn $(-3,45% \div +7,12%)$ trên tập 11 vỉa than chiếm 87% tài nguyên mỏ.
  5. Mở ra hướng nghiên cứu mới về địa thống kê ứng dụng trong công nghiệp than Việt Nam: Cung cấp tài liệu tham khảo kinh điển, quy chuẩn hóa quy trình khai thác dữ liệu số, góp phần hiện đại hóa ngành Kỹ thuật Địa chất và phát triển bền vững nền kinh tế khoáng sản quốc gia.