Tổng quan về luận án

Hoạt động khai thác và chế biến khoáng sản đóng vai trò là động lực then chốt trong tiến trình công nghiệp hóa và hiện đại hóa kinh tế tại Việt Nam. Tuy nhiên, sự phát triển nóng của ngành công nghiệp khai khoáng tại khu vực Bắc Trung Bộ, đặc biệt là địa bàn Nghệ An - Hà Tĩnh, đã và đang tạo ra những xung đột gay gắt đối với tính bền vững của các hệ sinh thái và môi trường sống dân sinh. Quá trình bóc phủ bề mặt, nổ mìn, đào đắp quy mô lớn cùng với sự tích tụ khổng lồ của các bãi thải đất đá và hồ lắng quặng đuôi đã làm biến đổi hoàn toàn cân bằng địa mạo - địa động lực vốn được thiết lập qua hàng triệu năm. Hậu quả trực tiếp là sự xuất hiện dồn dập của các dạng tai biến môi trường phức hợp như trượt lở sườn tầng khai thác, sạt lở bãi thải, sụt lún địa hình, vỡ đập quặng đuôi, lũ bùn đá và ô nhiễm phát tán kim loại nặng độc hại trên diện rộng.

Khoảng trống nghiên cứu (research gap) cốt lõi trước thời điểm luận án được công bố nằm ở chỗ: các công trình trước đây đa phần chỉ tiếp cận vấn đề theo hướng điều tra hiện trạng đơn lẻ hoặc đánh giá tác động môi trường (ĐTM) mang tính định tính cục bộ, thiếu vắng các mô hình toán học định lượng tích hợp có khả năng liên kết đa biến giữa các yếu tố nội sinh (cấu trúc địa chất kiến tạo, thạch học, địa chất thủy văn - địa chất công trình) và yếu tố ngoại sinh (lượng mưa cực đoan, công nghệ nổ mìn, độ dốc đổ thải). Hơn nữa, công tác mô phỏng không gian lan truyền chất ô nhiễm từ bãi thải và dự báo rủi ro vỡ hồ chứa quặng đuôi chưa được lượng hóa theo chuỗi thời gian thực.

Luận án của nghiên cứu sinh Nguyễn Thị Hòa (2020) thuộc chuyên ngành Kỹ thuật Môi trường (Mã số: 9520320) tại Trường Đại học Bách khoa Hà Nội, dưới sự hướng dẫn khoa học của TS. Trịnh Thành và PGS. Nguyễn Phương, đã giải quyết triệt để khoảng trống này thông qua hệ thống câu hỏi và giả thuyết nghiên cứu tường minh:

  • Câu hỏi nghiên cứu 1 (RQ1): Các tác nhân tự nhiên và hoạt động nhân sinh nào đóng vai trò kiểm soát chính đối với nguy cơ phát sinh tai biến môi trường mỏ tại khu vực Nghệ An - Hà Tĩnh?
  • Câu hỏi nghiên cứu 2 (RQ2): Mức độ phân hóa không gian và xác suất xảy ra trượt lở sườn tầng, bãi thải và lũ bùn đá có thể được mô hình hóa định lượng với độ tin cậy bao nhiêu khi kết hợp GIS, viễn thám và học máy/thống kê xác suất?
  • Câu hỏi nghiên cứu 3 (RQ3): Quy luật di chuyển, phát tán và tích tụ của các chất gây ô nhiễm (đặc biệt là kim loại nặng As, Cr6+, Pb, Mn) từ bãi thải quặng đuôi ra môi trường đất và nguồn nước xung quanh diễn ra như thế nào theo thời gian mô phỏng từ 250 đến 1.500 ngày?
  • Giả thuyết nghiên cứu 1 (H1): Sự kết hợp giữa các mô hình thống kê xác suất (Trọng số bằng chứng - WoE, Hệ số tin cậy - CF) và thuật toán trí tuệ nhân tạo (Mạng nơ-ron nhân tạo - ANN, Hồi quy Logistic) trên nền tảng GIS sẽ tạo ra bản đồ phân vùng nguy cơ tai biến trượt lở với diện tích dưới đường cong ROC (AUC) vượt trội so với các phương pháp đánh giá chuyên gia đơn thuần.
  • Giả thuyết nghiên cứu 2 (H2): Dòng thấm bão hòa và không bão hòa trong thân bãi thải quặng đuôi đóng vai trò là cơ chế động lực quyết định quá trình phát tán ô nhiễm kim loại nặng ra môi trường xung quanh theo quy luật bình lưu - phân tán (advection-dispersion).

Khung lý thuyết của luận án được xây dựng dựa trên sự tích hợp của Lý thuyết Thống kê Địa chất và Xác suất Không gian (Spatial Probability Theory), Lý thuyết Biến dạng Liên tục/Phi liên tục Cơ học Đất đá (Rock/Soil Mechanics & Geodynamics), và Lý thuyết Lan truyền Ô nhiễm Thủy địa hóa Môi trường (Environmental Solute Transport Theory). Luận án khảo sát trên quy mô thực địa toàn diện với 354 cơ sở khai thác và chế biến khoáng sản trên toàn vùng Nghệ An - Hà Tĩnh, tập trung chuyên sâu vào tâm điểm mỏ thiếc sa khoáng Bản Cô, Suối Bắc, Châu Hồng (huyện Quỳ Hợp) và mỏ sắt Thạch Khê, mỏ titan ven biển (tỉnh Hà Tĩnh). Đóng góp đột phá của luận án là việc xây dựng thành công bộ cơ sở dữ liệu số hóa không gian đa tầng, phân tầng chi tiết 4 dạng tai biến trọng yếu và kiến tạo 3 nhóm giải pháp quản trị rủi ro có tính ứng dụng thực tiễn cao.

Literature Review và Positioning

Lịch sử nghiên cứu tai biến môi trường mỏ trên thế giới đã ghi nhận sự chuyển dịch sâu sắc từ phương pháp thống kê mô tả sang mô hình hóa không gian số và cơ học tính toán. Trong nghiên cứu trượt lở và tai biến địa mạo mỏ, các học giả quốc tế đã phát triển nhiều trường phái tiếp cận. Trường phái thứ nhất đại diện bởi Bonham-Carter (1994) và Carranza (2008) với việc ứng dụng mô hình Trọng số bằng chứng (Weight of Evidence - WoE) dựa trên định lý xác suất Bayes, cho phép tính toán khách quan mối quan hệ không gian giữa các điểm tai biến lịch sử và các lớp bản đồ nhân tố môi trường. Trường phái thứ hai xuất phát từ Shortliffe & Buchanan (1975) với mô hình Hệ số tin cậy (Certainty Factor - CF), giúp xử lý hiệu quả các dữ liệu địa mạo không đồng nhất và độ bất định của dữ liệu địa kỹ thuật.

Trong thập kỷ gần đây, làn sóng thứ ba bùng nổ với các công trình của Pradhan & Lee (2010) và Pourghasemi et al. (2012), áp dụng Mạng nơ-ron nhân tạo (ANN) và Hồi quy Logistic (LR) để giải quyết các mối quan hệ phi tuyến phức tạp trong cơ chế kích hoạt trượt lở tại các khu vực mỏ miền núi. Đối với thảm họa vỡ đập quặng đuôi và lũ bùn đá, sau hàng loạt thảm họa kinh hoàng như sự cố vỡ đập Baia Mare tại Romania năm 2000 hay thảm họa hồ bùn đỏ Bento Rodrigues (Mariana) và Brumadinho tại Brazil (Rico et al., 2008; Carson, 2015), các nghiên cứu của Liên minh Quốc tế về Tai biến Địa chất (JTC1) và Ủy ban Liên chính phủ về Biến đổi Khí hậu (IPCC AR5) đã thúc đẩy việc ứng dụng phân tích thứ bậc AHP của Saaty (1980) kết hợp đánh giá đa tiêu chí (MCDA) để thiết lập bản đồ nhạy cảm sinh thái mỏ.

Tại Việt Nam, các nghiên cứu trước đây của nhiều tác giả chuyên ngành Địa chất Môi trường và Địa kỹ thuật thường gặp phải sự chia rẽ trong quan điểm học thuật:

  • Quan điểm 1: Trường phái địa chất công trình truyền thống cho rằng các yếu tố nội sinh (thạch học, đứt gãy kiến tạo, độ bền cắt của đất đá) là nguyên nhân tiên quyết gây ra tai biến sườn dốc và bãi thải.
  • Quan điểm 2: Trường phái môi trường - khí tượng học lại khẳng định chế độ mưa cực đoan của khí hậu nhiệt đới gió mùa và các tác động nhân sinh tiêu cực (nổ mìn, mở rộng tầng khai thác quá độ dốc cho phép, bạt đồi không giật cấp) mới là yếu tố quyết định kích hoạt thảm họa tức thời.

Luận án của Nguyễn Thị Hòa định vị vững chắc trong dòng chảy học thuật bằng cách bắc cầu nối hai quan điểm trên, tích hợp đồng thời các biến số cơ học nền móng và các tham số khí hậu - nhân sinh vào một khung phân tích thống nhất. Khi đối chiếu với hai nghiên cứu quốc tế điển hình là công trình đánh giá lan truyền ô nhiễm mỏ kim loại của Johnson & Hallberg (2005) tại Anh và mô hình hóa trượt lở vùng mỏ của Lee & Choi (2012) tại Hàn Quốc, nghiên cứu của luận án vượt trội về tính thích ứng với điều kiện địa hình nhiệt đới gió mùa phân cắt sâu, thảm thực vật biến động mạnh và tính phức tạp của các mỏ sa khoáng quy mô vừa và nhỏ tại các quốc gia đang phát triển.

Đó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 Động lực học Tai biến Địa môi trường (Geo-environmental Hazard Dynamics) bằng cách làm sáng tỏ cơ chế tương tác đa quy mô giữa áp lực khai thác khoáng sản và tính dễ bị tổn thương sinh thái. Về mặt lý thuyết trượt lở, nghiên cứu đóng góp mô hình tích hợp phi tuyến giải quyết bài toán mâu thuẫn giữa độ bất định dữ liệu thực địa và tính chính xác thống kê. Bằng việc chứng minh mối liên hệ hữu cơ giữa mật độ đứt gãy kiến tạo, chiều dày tầng phong hóa sét - bột và chế độ mưa tích lũy trên 24 giờ (>75mm và >100mm), luận án đã khái quát hóa quy luật ngưỡng cân bằng giới hạn của sườn dốc nhân tạo tại các bãi đổ thải mỏ.

Đối với lý thuyết lan truyền chất ô nhiễm trong môi trường xốp (Porous Media Solute Transport), luận án đã chuẩn hóa các phương trình bình lưu - phân tán trong điều kiện đất không bão hòa: $$\frac{\partial (\theta C)}{\partial t} = \frac{\partial}{\partial x_i} \left( \theta D_{ij} \frac{\partial C}{\partial x_j} \right) - \frac{\partial}{\partial x_i} (q_i C) - \lambda \theta C$$ Trong đó $\theta$ là độ ẩm thể tích, $C$ là nồng độ ion kim loại nặng, $D_{ij}$ là tensor hệ số phân tán thủy lực, $q_i$ là vận tốc dòng thấm Darcy, và $\lambda$ là hệ số phân rã/hấp phụ. Luận án chứng minh rằng thành phần hạt mịn (bột và sét chiếm tới 60-70% trong chất thải mỏ thiếc sa khoáng Bản Cô) đóng vai trò là ma trận lưu giữ và đồng thời là nguồn phát tán chậm các kim loại độc hại (As, Cr6+, Pb) vào hệ thống thủy văn ngầm khi áp lực nước lỗ rỗng thay đổi theo chu kỳ mùa khô - mùa mưa.

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

Khung phân tích của luận án được vận hành trên cơ sở liên kết 3 trụ cột lý thuyết:

  1. Lý thuyết Phân tích Thứ bậc (AHP) của Saaty (1980) để cấu trúc hóa trọng số chuyên gia;
  2. Lý thuyết Thống kê Xác suất Không gian (Bayesian WoE & CF) để triệt tiêu tính chủ quan trong đánh giá nhạy cảm tai biến;
  3. Lý thuyết Phần tử Hữu hạn (Finite Element Method - FEM) để mô phỏng động học quá trình thấm và dịch chuyển ô nhiễm qua bộ phần mềm GeoStudio 2012 (mô-đun SEEP/W và CTRAN/W).
+-------------------------------------------------------------------------------+
|                       KHUNG PHÂN TÍCH TAI BIẾN TÍCH HỢP                       |
+-------------------------------------------------------------------------------+
                                        |
      +---------------------------------+---------------------------------+
      |                                                                   |
+-----v-----------------------------+           +-------------------------v-----+
| CƠ SỞ DỮ LIỆU ĐA NGUỒN (SPATIAL)  |           | THÔNG SỐ CƠ LÝ & THỦY ĐỊA HÓA |
| - DEM 30m, Ảnh Landsat 8 OLI      |           | - Thí nghiệm cơ học đất đá    |
| - Bản đồ Địa chất, Đứt gãy, ĐCCT  |           | - Nước mặt, Nước ngầm, Nước thải|
| - Mật độ sông suối, NDVI          |           | - Mẫu đất, Thành phần cấp hạt |
| - Dữ liệu mưa lớn (24h, >75mm)    |           | - Độ ẩm, Hệ số thấm k(theta)  |
+-----------------+-----------------+           +-----------------+-------------+
                  |                                               |
                  +-----------------------+-----------------------+
                                          |
                                          v
+-------------------------------------------------------------------------------+
|                         HỆ PHƯƠNG PHÁP TÍNH TOÁN ĐỊNH LƯỢNG                   |
|  [AHP - Saaty]  -->  [WoE & CF Bayes]  -->  [ANN / Logistic]  -->  [CTRAN/W]  |
+-----------------------------------------+-------------------------------------+
                                          |
      +-----------------------------------+-----------------------------------+
      |                                   |                                   |
+-----v-------------------+     +---------v---------------+     +-------------v-+
|  BẢN ĐỒ NGUY CƠ TRƯỢT LỞ|     | BẢN ĐỒ NGUY CƠ LŨ BÙN ĐÁ|     | MÔ PHỎNG PHÁT |
| SƯỜN TẦNG & BÃI THẢI MỎ |     | DO VỠ HỒ QUẶNG ĐUÔI     |     | TÁN KIM LOẠI  |
| (ROC-AUC Kiểm định)     |     | (MCDA - Phân tích đa CT)|     | (250-1500 ngày|
+-------------------------+     +-------------------------+     +---------------+

Điều kiện biên và phạm vi hữu dụng của khung phân tích được xác định rõ: áp dụng tối ưu cho các phức hệ mỏ nằm trong vùng đồi núi có độ dốc từ 15° đến trên 45°, chịu tác động của chế độ lượng mưa nhiệt đới gió mùa tập trung cao độ, và có đặc tính quặng hóa đa kim hoặc sa khoáng gắn liền với đá vôi/đá hoa hóa và macma xâm nhập.

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

Thiết kế nghiên cứu

Luận án vận hành theo hệ hình Hiện thực Phê phán (Critical Realism) và Thực chứng Hiện đại (Post-positivism), kết hợp chặt chẽ giữa quan trắc thực nghiệm, phân tích viễn thám không gian và mô hình hóa toán học. Thiết kế nghiên cứu đa tầng bao gồm cấp vùng (toàn diện 2 tỉnh Nghệ An - Hà Tĩnh) và cấp mỏ chi tiết (mỏ thiếc sa khoáng Bản Cô - Quỳ Hợp, mỏ thiếc Suối Bắc, mỏ sắt Thạch Khê).

Quy mô mẫu khảo sát và điều tra thực tế đạt độ bao phủ toàn diện:

  • Tổng số cơ sở khai thác, chế biến được điều tra, lập danh mục cơ sở dữ liệu: 354 mỏ và xưởng chế biến (gồm 1 mỏ than đá 17.000 tấn/năm; 43 mỏ kim loại như sắt 584.000 tấn/năm, mangan 105.000 tấn/năm, vàng 80 kg/năm; 20 mỏ khoáng chất công nghiệp; 13 mỏ nguyên liệu xi măng; 242 mỏ vật liệu xây dựng thông thường; và 34 xưởng chế biến khoáng sản quy mô lớn).
  • Riêng tại huyện Quỳ Hợp: thống kê và định vị tọa độ chính xác 14 mỏ thiếc (tổng diện tích cấp phép 125.775,96 ha, sản lượng khai thác thực tế 795,37 tấn/năm), 33 mỏ đá vôi trắng (sản lượng 1.939,8 nghìn m3/năm), và 17 mỏ đá xây dựng thông thường (sản lượng 872,30 nghìn m3/năm).

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý mẫu môi trường tuân thủ nghiêm ngặt các Tiêu chuẩn Việt Nam (TCVN) và Quy chuẩn Kỹ thuật Quốc gia (QCVN):

  • Phân tích ảnh viễn thám: Sử dụng tư liệu ảnh vệ tinh Landsat 8 OLI (ảnh tổ hợp kênh 5-4-1 chụp ngày 10/02/2018 và ảnh đa thời kỳ), giải đoán các biến động bề mặt moong khai thác, vị trí bãi thải lộ thiên, cấu trúc đứt gãy kiến tạo (lineament) và lớp phủ thực vật thông qua chỉ số NDVI.
  • Quan trắc và phân tích môi trường: Lấy mẫu định kỳ và phân tích phòng thí nghiệm đạt chuẩn VILAS đối với mẫu nước thải (NT), nước mặt, nước ngầm, mẫu đất (Đ) tại các tầng mặt (0-10cm) và tầng sâu (10-30cm), mẫu trầm tích và mẫu bụi lơ lửng.
  • Thử nghiệm cơ lý đất đá: Xác định thành phần cấp hạt bằng phương pháp rây và ống hút cơ giới; xác định đường cong đặc trưng đất - nước (Soil-Water Characteristic Curve - SWCC) và hàm hệ số thấm không bão hòa $k(\theta)$ phục vụ bài toán địa kỹ thuật.
  • Quy trình tam giác hóa (Triangulation): Kết hợp dữ liệu khảo sát lộ trình thực địa, dữ liệu viễn thám độ phân giải cao và kết quả phỏng vấn chuyên gia địa chất - môi trường để hiệu chỉnh các ma trận đánh giá rủi ro.

Data và phân tích

Phân tích định lượng được hỗ trợ bởi các công cụ phần mềm chuyên dụng hàng đầu: ArcGIS 10.4, GeoStudio 2012 (mô-đun SEEP/W và CTRAN/W), SPSS và MATLAB.

Để tính toán trọng số và phân vùng nguy cơ trượt lở, luận án đã thiết lập ma trận so sánh cặp đôi theo AHP của Saaty với 14 lớp nhân tố thành phần: độ dốc địa hình, hướng sườn dốc, độ cao DEM, mật độ đứt gãy kiến tạo, đặc tính thạch học, phân loại đất, chiều dày tầng phong hóa, thành phần cơ giới đất, điều kiện địa chất công trình (ĐCCT), điều kiện địa chất thủy văn (ĐCTV), mật độ mạng lưới sông suối, chỉ số thực vật (NDVI), lượng mưa lớn nhất 24 giờ, và số ngày mưa lớn (>75mm và >100mm).

Kiểm định tính nhất quán của ma trận trọng số AHP đảm bảo Tỷ số nhất quán (Consistency Ratio - CR) đạt $CR < 0,10$, đáp ứng hoàn hảo tiêu chuẩn lý thuyết của Saaty: $$CI = \frac{\lambda_{\max} - n}{n - 1}, \quad CR = \frac{CI}{RI}$$ Độ chính xác của các bản đồ dự báo trượt lở xây dựng theo mô hình Trọng số bằng chứng (Bayes - WoE), Hệ số tin cậy (CF) và Mạng nơ-ron nhân tạo (ANN) được kiểm định chéo bằng ma trận sai số (Error Matrix) và phân tích đường cong đặc trưng hoạt động của bộ thu (ROC curve). Kết quả cho thấy mô hình ANN và WoE đạt hệ số diện tích dưới đường cong (AUC) trên 85%, khẳng định độ tin cậy vượt trội trong việc phân tách các vùng có nguy cơ trượt lở từ "Rất thấp", "Thấp", "Trung bình", "Cao" đến "Rất cao".

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

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

  1. Sự vượt ngưỡng nghiêm trọng của nồng độ kim loại nặng và độc chất trong đất và nước tại các khu vực mỏ: Kết quả phân tích mẫu đất tại bãi thải và khu vực phụ cận mỏ thiếc sa khoáng Bản Cô (xã Châu Thành, Quỳ Hợp) chứng minh hàm lượng Asen (As) ở lớp đất mặt (0-10cm) tại điểm Đ1-1 đạt tới 189,4 mg/kg, vượt 12,6 lần so với giới hạn tối đa cho phép của QCVN 03-MT:2015/BTNMT (quy chuẩn đất nông nghiệp là 15 mg/kg). Hàm lượng Crom hóa trị 6 ($Cr^{6+}$) đạt 315 mg/kg (vượt 2,1 lần quy chuẩn 150 mg/kg); Mangan (Mn) đạt 1.435 mg/kg; Chì (Pb) đạt 107 mg/kg (vượt quy chuẩn 70 mg/kg); Đồng (Cu) đạt 112 mg/kg (vượt quy chuẩn 100 mg/kg); Kẽm (Zn) đạt 210 mg/kg (vượt quy chuẩn 200 mg/kg).

  2. Ô nhiễm nghiêm trọng trong dòng nước thải mỏ chưa qua xử lý chuẩn: Mẫu nước thải mỏ (NT1-1) tại moong và hồ lắng ghi nhận nồng độ Tổng chất rắn lơ lửng (TSS) lên tới 105 mg/l (vượt 2,1 lần so với QCVN 40:2011/BTNMT, Cột B là 50 mg/l); nồng độ Mangan (Mn) hòa tan đạt 5,74 mg/l (vượt 11,48 lần so với quy chuẩn 0,5 mg/l); Asen đạt 0,36 mg/l (vượt 7,2 lần so với quy chuẩn 0,05 mg/l); Cadmi (Cd) đạt 0,025 mg/l; độ pH nước thải giảm xuống mức 5,4 thể hiện tính axit nhẹ kích hoạt sự hòa tan kim loại.

  3. Phát hiện động học lan truyền chất ô nhiễm theo chuỗi thời gian 1.500 ngày: Mô phỏng phần tử hữu hạn 2D bằng CTRAN/W trên mặt cắt ngang bãi thải mỏ thiếc Bản Cô đã chỉ ra bức tranh trực quan về quá trình dịch chuyển chất ô nhiễm qua môi trường xốp. Tại mốc thời gian 250 ngày, chất ô nhiễm tập trung cục bộ tại đáy bãi thải. Đến 500 và 750 ngày, nồng độ ô nhiễm thâm nhập sâu vào tầng chứa nước ngầm lỗ rỗng. Đến mốc 1.000, 1.250 và 1.500 ngày, quầng ô nhiễm (plume) As và Mn đã phát tán lan rộng ra ngoài ranh giới bãi thải hàng trăm mét theo phương dòng chảy ngầm hướng về các thung lũng canh tác lúa và sông suối tự nhiên.

  4. Quy luật kích hoạt tai biến lũ bùn đá do vỡ đập quặng đuôi: Phân tích cấu trúc đập quặng đuôi tại vùng Quỳ Hợp (tiêu biểu như khu vực bản Na Kỳ, Suối Bắc) bộc lộ sự mất an toàn nghiêm trọng. Hầu hết các bờ bao hồ lắng được đắp sơ sài bằng chính đất đá thải có góc dốc tự nhiên cao, thành phần cơ giới chứa 38-45% hạt sét và 24-28% hạt bột. Khi lượng mưa ngày vượt 100mm, hiện tượng bão hòa nước làm suy giảm nhanh chóng lực dính kết ($c$) và góc ma sát trong ($\varphi$), đẩy hệ số an toàn mái dốc ($F_s$) xuống dưới 1,0, dẫn đến nguy cơ vỡ đập dây chuyền và hình thành lũ bùn đá tàn phá hạ lưu.

Implications đa chiều

  • Ý nghĩa lý thuyết: Nghiên cứu đã chứng minh rằng các mô hình tĩnh truyền thống không thể phản ánh chính xác rủi ro tai biến mỏ nếu không tích hợp động lực học dòng thấm không bão hòa và biến động khí tượng thủy văn cực đoan. Sự kết hợp thành công giữa lý thuyết độ tin cậy địa chất và mô hình phần tử hữu hạn cung cấp một bước tiến học thuật quan trọng trong ngành Địa kỹ thuật Môi trường.
  • Ý nghĩa phương pháp luận: Khung quy trình tích hợp AHP - WoE - GIS - CTRAN/W có thể chuyển giao và ứng dụng trực tiếp cho tất cả các vùng mỏ kim loại và vật liệu xây dựng trên toàn lãnh thổ Việt Nam và các khu vực có điều kiện địa hình tương đồng ở Đông Nam Á.
  • Ứng dụng thực tiễn và chính sách: Luận án đề xuất 3 nhóm giải pháp đồng bộ:
    1. Giải pháp công trình: Bắt buộc áp dụng công nghệ nổ mìn vi sai phi điện để giảm chấn động địa chấn kích hoạt nứt nẻ sườn tầng; xây dựng hệ thống đê quai chắn bùn cát có tầng lọc ngược bằng địa kỹ thuật; cắt tầng giật cấp bãi thải với góc dốc sườn tầng $< 28^\circ$; gia cố chân bãi thải bằng rọ đá; xây dựng mương phân lũ tách nước mặt không cho chảy tràn qua bãi thải.
    2. Giải pháp phi công trình: Thành lập mạng lưới trạm quan trắc tự động áp lực nước lỗ rỗng và độ dịch chuyển sườn tầng; ứng dụng viễn thám cập nhật biến động khai thác trái phép định kỳ 6 tháng/lần; thiết lập ngưỡng cảnh báo sớm lũ bùn đá theo lượng mưa thực đo.
    3. Giải pháp cơ chế - chính sách: Tăng mức tiền ký quỹ phục hồi môi trường mỏ; kiên quyết đình chỉ và không gia hạn giấy phép đối với các mỏ không thực hiện hoàn thổ và trồng cây phục hồi môi trường sau khai thác (như thực trạng tại mỏ titan Thạch Văn và Cẩm Hòa - Hà Tĩnh).

Limitations và Future Research

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

  1. Độ phân giải của dữ liệu đầu vào: Mô hình số độ cao DEM sử dụng có độ phân giải 30m, do đó các biến đổi vi địa hình cục bộ ở quy mô sườn tầng nhỏ ($< 10m$) chưa được phản ánh chi tiết tuyệt đối;
  2. Số lượng mẫu phân tích theo chiều sâu: Các mẫu đất phân tích mới chỉ dừng lại ở hai tầng phát sinh mặt (0-10cm) và tầng sâu (10-30cm), chưa thực hiện được các giếng khoan địa chất sâu qua toàn bộ chiều dày đới phong hóa để quan trắc biến thiên thủy địa hóa liên tục 3 chiều (3D);
  3. Giả định đồng nhất trong mô hình hóa CTRAN/W: Mô hình lan truyền chất ô nhiễm 2D giả định môi trường xốp có tính đẳng hướng tương đối, trong khi thực tế cấu trúc đứt gãy và nứt nẻ của đá gốc có tính dị hướng cao.

Chương trình nghiên cứu tương lai cần mở rộng:

  • Ứng dụng công nghệ quét LiDAR và thiết bị bay không người lái (UAV) để xây dựng mô hình số độ cao DEM độ phân giải siêu cao (dưới 1m) phục vụ phân tích ổn định sườn tầng mỏ theo thời gian thực;
  • Phát triển mô hình mô phỏng 3D tích hợp (3D coupled Hydro-Mechanical-Chemical modeling) để lượng hóa sự tương tác giữa biến dạng cơ học đất đá và dòng dịch chuyển hóa chất độc hại;
  • Đánh giá định lượng rủi ro sức khỏe cộng đồng (Human Health Risk Assessment - HHRA) thông qua việc phân tích chuỗi thức ăn, sự tích lũy sinh học của kim loại nặng trong cây trồng nông nghiệp và nguồn nước ngầm sinh hoạt tại hạ lưu các mỏ khoáng sản.

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

Công trình nghiên cứu mang lại tác động sâu rộng trên nhiều bình diện:

  • Tác động học thuật: Cung cấp bộ dữ liệu thực chứng chuẩn xác và phương pháp luận mẫu mực cho các nghiên cứu sinh, học viên cao học chuyên ngành Kỹ thuật Môi trường, Địa chất Thủy văn - Địa chất Công trình; nâng cao vị thế khoa học của các công trình công bố trong nước và quốc tế (tiêu biểu như các bài báo đã công bố tại Hội nghị ERSD 2018 và Tuần lễ Nước Quốc tế VACI 2019).
  • Tác động đối với ngành công nghiệp khai khoáng: Giúp các doanh nghiệp khai khoáng (như Công ty CP Kim loại màu Nghệ Tĩnh, các xí nghiệp khai thác đá hoa trắng và đá xây dựng) định hình lại quy trình đổ thải an toàn, thiết kế đập lắng quặng đuôi đạt chuẩn kỹ thuật, giảm thiểu tối đa tổn thất tài sản do sự cố sạt lở.
  • Tác động chính sách và xã hội: Cung cấp luận cứ khoa học thực chứng để Sở Tài nguyên và Môi trường tỉnh Nghệ An và Hà Tĩnh rà soát, thẩm định báo cáo ĐTM, đề án cải tạo phục hồi môi trường, và quy hoạch các vùng cấm hoặc tạm thời cấm hoạt động khoáng sản nhằm bảo vệ an toàn tính mạng và sinh kế bền vững cho hàng chục ngàn hộ dân địa phươ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 khung phương pháp luận định lượng hoàn chỉnh, kết hợp viễn thám, GIS, học máy và mô hình phần tử hữu hạn để giải quyết các bài toán địa môi trường phức tạp.
  • Cơ quan Quản lý Nhà nước (Bộ TN&MT, UBND các tỉnh Nghệ An, Hà Tĩnh): Sở hữu cơ sở dữ liệu số hóa và bản đồ phân vùng rủi ro chi tiết phục vụ công tác thanh tra, cấp phép và quy hoạch khai thác khoáng sản bền vững.
  • Doanh nghiệp Khai thác Mỏ: Có được cẩm nang kỹ thuật thực tiễn về quy chuẩn góc dốc sườn tầng, quy trình nổ mìn vi sai phi điện, và công nghệ xử lý ổn định đập quặng đuôi nhằm tối ưu hóa chi phí vận hành và phòng ngừa sự cố thảm họa.
  • Cộng đồng Dân cư Bản địa: Hưởng lợi trực tiếp từ việc giảm thiểu ô nhiễm nguồn nước ngầm, nước mặt và ngăn ngừa thảm họa trượt lở, lũ bùn đá đe dọa khu dân cư hạ lưu 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ì? Luận án đã mở rộng Lý thuyết Địa cơ học Môi trường bằng việc thiết lập mô hình tích hợp phi tuyến giải mã cơ chế kích hoạt trượt lở và lũ bùn đá do vỡ đập thải quặng đuôi tại vùng nhiệt đới gió mùa, liên kết định lượng giữa hàm lượng hạt mịn (sét - bột chiếm 60-70%), biến động áp lực nước lỗ rỗng và ngưỡng mưa cực đoan (>100mm/24h).

  2. Đổi mới phương pháp luận của luận án so với các nghiên cứu trước đây? So với nghiên cứu của Bonham-Carter (1994) và Rico et al. (2008), luận án không chỉ dừng lại ở đánh giá thống kê không gian tĩnh mà đã tích hợp chuỗi công cụ liên hoàn: từ giải đoán viễn thám Landsat 8 OLI, tính toán trọng số thứ bậc AHP và xác suất Bayes (WoE, CF), đến mô phỏng động học 2D dòng thấm và lan truyền ô nhiễm theo thời gian bằng mô hình phần tử hữu hạn CTRAN/W (GeoStudio 2012).

  3. Phát hiện thực nghiệm gây bất ngờ và có giá trị cảnh báo cao nhất? Nồng độ Asen tích tụ trong lớp đất mặt tại bãi thải mỏ thiếc Bản Cô lên tới 189,4 mg/kg (vượt 12,6 lần quy chuẩn QCVN 03-MT:2015), và kết quả mô phỏng CTRAN/W chứng minh quầng ô nhiễm As, Mn thâm nhập sâu vào tầng chứa nước ngầm và lan truyền nhanh chóng ra ngoài phạm vi mỏ sau 1.500 ngày, đe dọa trực tiếp các diện tích đất lúa và nguồn nước sinh hoạt hạ lưu.

  4. Luận án có cung cấp quy trình tái lập nghiên cứu (Replication Protocol) không? Toàn bộ quy trình từ thu thập, chuẩn hóa dữ liệu GIS, các bước chạy mô hình AHP, tính chỉ số CR, cài đặt điều kiện biên và thông số địa kỹ thuật trong GeoStudio CTRAN/W đều được mô tả chi tiết, minh bạch tại Chương 2 và Chương 3, cho phép các nhà khoa học khác tái lập hoàn toàn quy trình trên các vùng mỏ khác.

  5. Chương trình nghiên cứu 10 năm tiếp theo được định hình ra sao? Trọng tâm 10 năm tới là phát triển hệ thống giám sát cảnh báo sớm tai biến mỏ thông minh dựa trên mạng lưới cảm biến IoT kết hợp trí tuệ nhân tạo (AI) phân tích ảnh vệ tinh radar khẩu độ tổng hợp (InSAR) để đo đạc dịch chuyển sườn dốc ở mức milimet theo thời gian thực.

Kết luận

Luận án của nghiên cứu sinh Nguyễn Thị Hòa là công trình khoa học công phu, giải quyết xuất sắc bài toán cấp thiết về an toàn môi trường mỏ tại khu vực Nghệ An - Hà Tĩnh với các đóng góp cốt lõi:

  1. Xác lập cơ sở dữ liệu số hóa không gian toàn diện về 354 cơ sở khai thác, chế biến khoáng sản và hiện trạng tai biến môi trường trên địa bàn hai tỉnh Nghệ An và Hà Tĩnh;
  2. Chứng minh thực nghiệm mức độ ô nhiễm vượt chuẩn cực kỳ nghiêm trọng của kim loại nặng trong đất (As đạt 189,4 mg/kg) và nước thải mỏ (Mn đạt 5,74 mg/l; TSS đạt 105 mg/l);
  3. Ứng dụng thành công hệ mô hình định lượng tích hợp AHP, WoE, CF, ANN và CTRAN/W, tạo bước đột phá trong việc mô hình hóa không gian nguy cơ trượt lở sườn tầng và động học phát tán chất ô nhiễm theo chuỗi thời gian 1.500 ngày;
  4. Khoanh định chính xác các phân vùng rủi ro tai biến môi trường trọng điểm (vùng Quỳ Hợp, Tương Dương, Thạch Khê);
  5. Đề xuất hệ thống 3 nhóm giải pháp giảm thiểu mang tính khả thi cao: công trình (nổ mìn vi sai phi điện, cắt tầng bãi thải $< 28^\circ$, đê lọc ngược), phi công trình (quan trắc cảnh báo sớm, viễn thám giám sát) và cơ chế chính sách hoàn thổ nghiêm ngặt;
  6. Mở ra hướng tiếp cận học thuật mới cho ngành Kỹ thuật Môi trường và Địa kỹ thuật tại Việt Nam, đóng góp thiết thực vào mục tiêu phát triển bền vững ngành công nghiệp khoáng sản quốc gia.