Tổng quan về luận án

Hoạt động công nghiệp khai thác và chế biến khoáng sản đóng vai trò đòn bẩy then chốt đối với tăng trưởng kinh tế và công nghiệp hóa, nhưng đồng thời là tác nhân nhân sinh mạnh mẽ nhất làm biến đổi sâu sắc địa chất môi trường. Tại khu vực Bắc Trung Bộ, đặc biệt là hai tỉnh Nghệ An và Hà Tĩnh, sự tập trung mật độ cao các điểm mỏ khai thác đá hoa trắng, đá xây dựng thông thường, sa khoáng thiếc (khu vực Quỳ Hợp) và quặng sắt, sa khoáng titan ven biển (khu vực Thạch Hà - Cẩm Xuyên, Thạch Khê) đã phá vỡ thế cân bằng sinh thái địa chất vốn được tích tụ qua hàng chục triệu năm. Luận án tiến sĩ kỹ thuật môi trường của Nghiên cứu sinh Nguyễn Thị Hòa với đề tài "Nghiên cứu đánh giá tai biến môi trường liên quan đến hoạt động khai thác khoáng sản khu vực Nghệ An - Hà Tĩnh và xây dựng giải pháp giảm thiểu" (thuộc Đề tài cấp Bộ mã số B2014-02-212) đặt nền móng tiên phong trong việc tiếp cận định lượng hóa nguy cơ tai biến địa chất - môi trường mỏ.

Khoảng trống nghiên cứu (research gap) cốt lõi được nhận diện rõ nét: Các công trình trước đây tại Việt Nam phần lớn tiếp cận ở quy mô vùng mang tính mô tả định tính đơn lẻ, phân tích các thông số ô nhiễm riêng biệt hoặc chỉ dừng lại ở các đánh giá tác động môi trường (ĐTM) tĩnh, thiếu tính liên kết động lực học không gian - thời gian. Chưa có một hệ phương pháp tích hợp đa quy mô có khả năng định lượng hóa nguy cơ tích hợp giữa biến động cơ học (trượt lở moong khai thác, vỡ đập quặng đuôi, lũ bùn đá) và quá trình lan truyền địa hóa chất ô nhiễm (kim loại nặng, dòng thải axit mỏ, phóng xạ tự nhiên) từ bãi thải vào hệ thống thủy văn - thổ nhưỡng.

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

  • Câu hỏi nghiên cứu 1 (RQ1): Những yếu tố tự nhiên nền và yếu tố kích hoạt nhân sinh nào đóng vai trò kiểm soát chính đối với các dạng tai biến trượt lở, lũ bùn đá và ô nhiễm môi trường do khai khoáng tại Nghệ An - Hà Tĩnh?
  • Câu hỏi nghiên cứu 2 (RQ2): Mức độ chính xác và tính tương thích của các mô hình toán thống kê không gian (Bayes WoE, CF, LR, ANN) và mô hình phần tử hữu hạn (CTRAN/W) trong việc mô phỏng lan truyền chất ô nhiễm và dự báo phân vùng nguy cơ tai biến môi trường đạt ngưỡng độ tin cậy nào?
  • Câu hỏi nghiên cứu 3 (RQ3): Làm thế nào để xây dựng một khung giải pháp giảm thiểu mang tính đồng bộ từ can thiệp công trình, công nghệ nổ mìn phi điện, kiểm soát quặng đuôi đến cơ chế quản trị rủi ro cấp chính sách?

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

  • Giả thuyết H1: Không gian phân bố của các vụ trượt lở sườn tầng khai thác và bãi thải mỏ có mối tương quan phi tuyến chặt chẽ với mật độ nứt nẻ/đứt gãy kiến tạo, độ dốc địa hình và chỉ số thực vật chuẩn hóa (NDVI) bị biến động do bóc phủ moong mỏ.
  • Giả thuyết H2: Mô hình mạng nơ-ron nhân tạo (ANN) và thống kê xác suất Bayes (WoE) mang lại độ chính xác dự báo (AUC/ROC) vượt trội so với các phương pháp đánh giá chuyên gia truyền thống khi tích hợp cơ sở dữ liệu viễn thám đa thời gian.
  • Giả thuyết H3: Quá trình di chuyển và lan truyền chất ô nhiễm từ bãi thải quặng đuôi (Bản Cô, Châu Hồng) vào tầng nước ngầm và lưu vực nước mặt tuân theo phương trình thấm - tán xạ phi bão hòa chịu sự chi phối trực tiếp của đường cong đặc trưng đất - nước (SWCC) và lượng mưa cực đoan 24 giờ.

Khung lý thuyết của công trình tích hợp lý thuyết đánh giá rủi ro tai biến địa chất của Ủy ban Hỗn hợp Quốc tế JTC1 (2009), Khung đánh giá rủi ro biến đổi khí hậu AR5 của IPCC (2014), cùng các nguyên lý nhiệt động học thấm và truyền chất trong môi trường rỗng. Phạm vi không gian nghiên cứu bao quát các điểm mỏ trọng điểm tại Quỳ Hợp (Châu Hồng, Châu Tiến, Châu Lộc, Liên Hợp, bản Na Kỳ) và ven biển Hà Tĩnh (Thạch Khê, Thạch Đỉnh, Thạch Hải, Thạch Văn). Kết quả nghiên cứu tạo ra tác động định lượng đột phá: cung cấp bộ bản đồ số hóa tỷ lệ 1:25.000 đến 1:50.000 phân vùng 04 nhóm tai biến trọng yếu, thiết lập ngưỡng an toàn môi trường mỏ phục vụ trực tiếp cho thẩm định quy hoạch và cấp phép khai thác khoáng sản bền vững.


Literature Review và Positioning

Tổng quan y văn quốc tế ghi nhận sự chuyển dịch mô thức mạnh mẽ từ các phương pháp khảo sát thực địa địa mạo truyền thống sang mô hình hóa số và học máy không gian. Trong giai đoạn tiền 1990, các nghiên cứu kinh điển như phân loại trượt đất của Nemcok và cộng sự (1972) tại Tiệp Khắc, lập bản đồ địa mạo vùng Sendai tỷ lệ 1:200.000 của Hatano và cộng sự (1974), hay phân vùng trượt lở mỏ than Nam Xứ Wales của Conway (1980) chủ yếu dựa vào đo vẽ kiểm kê thực địa tốn nhiều thời gian và mang tính lịch sử cục bộ.

Từ sau năm 2000, sự phát triển vượt bậc của công nghệ GIS và viễn thám mở ra kỷ nguyên định lượng hóa xác suất:

  1. Trường phái thống kê xác suất và trọng số bằng chứng (Weight of Evidence - WoE, Certainty Factor - CF): Tiêu biểu với công trình của Chung và Fabbri (2001), Lan và cộng sự (2004) tại lưu vực sông Xiaojiang (Trung Quốc), hay Christos Chalkias và cộng sự (2014) tại bán đảo Peloponnese (Hy Lạp) khi thiết lập Chỉ số nhạy cảm trượt lở đất (LSI) kết hợp dữ liệu lượng mưa cực đoan.
  2. Trường phái mô hình hóa số địa kỹ thuật và động lực học dòng thấm (Geotechnical & Numerical Modeling): Sử dụng các mô hình phần tử hữu hạn (FEM), phân tích biến dạng không liên tục (DDA), mô hình SINMAP, CHASM, TRIGRS và bộ phần mềm GeoStudio (SLOPE/W, SEEP/W, CTRAN/W) được phát triển bởi các nhóm nghiên cứu như Kumar và Sanoujam (2006), Gilson (2008), GEO-Slope (2018).
  3. Trường phái mô hình đa tiêu chí không gian (Spatial Multi-Criteria Analysis - SMCA/AHP): Khởi xướng bởi Saaty (1980), Ayalew (2005), Castellanos và Van Westen (2007) nhằm xử lý các vùng có dữ liệu quan trắc lịch sử hạn chế bằng logic mờ và ma trận so sánh cặp.

Tại Việt Nam, các nghiên cứu môi trường mỏ trước đây thường tập trung vào hiện trạng ô nhiễm hóa học đơn lẻ của nước mặt hoặc khí bụi công nghiệp (như các nghiên cứu của Liên đoàn Địa chất Xạ - Hiếm, các đề tài ĐTM riêng lẻ cho từng mỏ đá, mỏ thiếc). Luận án của NCS. Nguyễn Thị Hòa định vị chính xác điểm giao thoa giữa 03 trường phái trên: tích hợp viễn thám đa thời kỳ (Landsat 5 TM, 7 ETM+, 8 OLI) để phát hiện biến động lớp phủ và khai thác trái phép; áp dụng so sánh chéo 04 thuật toán định lượng (Bayes WoE, CF, LR, ANN) để triệt tiêu tính chủ quan của chuyên gia; đồng thời sử dụng mô hình phần tử hữu hạn 2D CTRAN/W để mô phỏng chính xác đường biên lan truyền chất ô nhiễm từ quặng đuôi theo các mốc thời gian từ 250 ngày đến 2000 ngày.

So sánh với các nghiên cứu quốc tế điển hình:

  • So với mô hình trượt lở bán đảo Peloponnese của Chalkias và cộng sự (2014): Nghiên cứu này không chỉ dừng lại ở bản đồ LSI tự nhiên mà tích hợp trực tiếp các tham số kích hoạt do nhân sinh (nổ mìn, cắt tầng bờ moong, bãi thải quặng đuôi không đầm nén).
  • So với mô hình phát tán chất ô nhiễm mỏ tại Châu Phi (Dự án UNESCO/IUGS, 2004) hay thung lũng Tully (Mỹ, Pair & Kappel, 2002): Luận án đã làm sáng tỏ cơ chế bục nước, thẩm thấu chất ô nhiễm trong điều kiện nhiệt đới gió mùa có lượng mưa 24 giờ vượt 100mm, phản ánh đặc thù địa chất karst hóa tại Quỳ Hợp và vùng cát biển ven bờ Thạch Khê.

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

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

Luận án đóng góp quan trọng vào lý thuyết Địa chất môi trường mỏ (Mining Environmental Geology) và Lý thuyết Tai biến tự nhiên tích hợp nhân sinh (Anthropogenic Geo-hazard Theory):

  1. Mở rộng mô hình lý thuyết Đánh giá Tai biến Môi trường: Luận án chuẩn hóa khái niệm: "Tai biến môi trường là biểu hiện về điều kiện, hoàn cảnh, hiện tượng, vụ việc hoặc quá trình, được xuất hiện, diễn biến trong thiên nhiên, trong xã hội, có tiềm năng gây hại, gây nguy hiểm đe dọa đối với an toàn sức khỏe, tính mạng con người, tài sản kinh tế, tài sản văn hóa - xã hội... hoặc phá vỡ tính ổn định, an toàn mang tính hệ thống môi trường tự nhiên, môi trường văn hóa xã hội và môi trường nhân sinh".
  2. Khung lý thuyết tương tác Yếu tố Nền - Yếu tố Kích thích (Baseline vs. Triggering Factors): Luận án chứng minh trên phương diện toán học rằng các yếu tố nền (đặc điểm thạch học, mức độ phong hóa, mật độ đứt gãy kiến tạo) đóng vai trò điều kiện cần, xác định ngưỡng chịu tải của sườn dốc; trong khi các yếu tố kích thích nhân sinh (góc dốc cắt tầng vượt thiết kế, rung chấn nổ mìn, sự gia tải của đống thải quặng đuôi kết hợp mưa bão >75-100mm/24h) là điều kiện đủ kích hoạt tai biến tức thời.
  3. Mô hình hóa dịch chuyển chất ô nhiễm qua môi trường rỗng không bão hòa: Mở rộng lý thuyết lan truyền chất ô nhiễm truyền thống bằng cách kết hợp cơ chế thấm không bão hòa với hàm lực hút dính ma trận (matric suction) và đường cong đặc trưng đất - nước (Soil-Water Characteristic Curve - SWCC) đối với bùn thải quặng sa khoáng.

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

Khung phân tích của luận án tích hợp đa tầng dữ liệu:

  • Tích hợp liên ngành: Kết hợp thạch học địa chất, địa chất thủy văn - địa chất công trình, viễn thám vi sai đa phổ và thủy lực dòng thấm.
  • Tiếp cận tổng hợp 3 bước: Tiếp cận thừa kế (phân tích lịch sử tai biến trượt lở, bục nước, vỡ đập quặng đuôi) $\rightarrow$ Tiếp cận phát sinh (lượng hóa cơ chế tác động của độ dốc, lượng mưa, nứt nẻ qua GIS) $\rightarrow$ Tiếp cận tổng hợp (chồng ghép lớp thông tin có trọng số qua AHP, Bayes WoE, ANN và CTRAN/W).
  • Điều kiện biên xác định (Boundary Conditions): Mô hình được thiết lập với các ngưỡng lượng mưa kích hoạt cụ thể (lượng mưa 24h cực đại, chu kỳ số ngày mưa liên tục >75mm và >100mm), góc dốc tầng khai thác từ $30^\circ$ đến $75^\circ$, và hệ số thấm $k$ biến thiên theo độ chặt của bãi thải mỏ.

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái Thực chứng luận (Positivism) kết hợp Hiện thực phê phán (Critical Realism), sử dụng phương pháp hỗn hợp (Mixed-Methods) với trọng tâm là định lượng hóa không gian:

  • Quy mô đa cấp: Cấp vùng (toàn tỉnh Nghệ An - Hà Tĩnh, tỷ lệ 1:100.000 đến 1:50.000) và Cấp điểm mỏ chi tiết (Quỳ Hợp, Thạch Khê, tỷ lệ 1:25.000 đến 1:10.000).
  • Mẫu khảo sát thực địa toàn diện: Khảo sát 100% các điểm mỏ đang hoạt động, mỏ đóng cửa, mỏ chưa cấp phép và các xưởng chế biến khoáng sản tại 2 cụm trọng điểm; tiến hành đo đạc 94 điểm suất liều bức xạ gamma môi trường, 96 điểm nồng độ radon trong đất/không khí/nước tại mỏ granit Nam Giới.

Quy trình nghiên cứu rigorous

  1. Thu thập dữ liệu viễn thám đa thời gian:
    • Sử dụng ảnh vệ tinh Landsat 5 TM (giai đoạn 1990-1995), Landsat 7 ETM+ (2000-2005) và Landsat 8 OLI (2015-2017) thuộc các cảnh Path/Row 126/47 và 127/47 với độ phủ mây kiểm soát dưới 10%.
    • Nắn chỉnh hình học sai số $<1\text{ pixel}$; phân loại lớp phủ sử dụng thuật toán Phân loại xác suất cực đại (Maximum Likelihood Classification - MLC) và Ánh xạ góc phổ (Spectral Angle Mapper - SAM) để nhận diện các moong khoáng sản và bãi thải lộ thiên sáng màu, loang lổ.
  2. Thiết bị quan trắc hiện trường chuẩn quốc tế:
    • Định vị không gian: Máy định vị vệ tinh GPS 12 Garmin.
    • Quan trắc khí độc hại: Bộ thiết bị DSTOX-ST-PPC và GrayWolf (USA).
    • Đo nồng độ bụi: Bơm lấy mẫu bụi thể tích chuẩn BD Portable Air Sampling Pump Model 5959P (USA).
    • Đo ồn và chấn động: Máy đo mức âm tích phân Extech 407780, Extech 45160 (Đài Loan), thiết bị NL-04.
    • Quan trắc vi khí hậu: Đầu đo kỹ thuật số TESTO 615.
  3. Tiêu chuẩn kiểm chuẩn chất lượng môi trường:
    • Nước mặt: QCVN 08-MT:2015/BTNMT.
    • Nước ngầm: QCVN 09-MT:2015/BTNMT.
    • Nước thải công nghiệp: QCVN 40:2011/BTNMT.
    • Đất và chất thải rắn: QCVN 03-MT:2015/BTNMT.
    • Không khí xung quanh và tiếng ồn: QCVN 05:2013/BTNMT, QCVN 26:2010/BTNMT.

Data và phân tích

Hệ thống xử lý phân tích dữ liệu đa biến bao gồm:

  • Mô hình toán thống kê nhạy cảm trượt lở: Sử dụng Tỷ số tần suất (Frequency Ratio - FR), Trọng số bằng chứng Bayes (Weight of Evidence - WoE), Hệ số tin cậy (Certainty Factor - CF), Hồi quy Logistic (LR) và Mạng nơ-ron nhân tạo đa tầng truyền thẳng (ANN - Multi-layer Perceptron Back-propagation). Đánh giá độ chính xác qua ma trận sai số (Confusion Matrix) và đường cong đặc trưng hoạt động của bộ thu (ROC - AUC).
  • Mô hình đánh giá đa chỉ tiêu (MCDA - Saaty AHP): Thiết lập ma trận so sánh cặp đôi cho 12 chỉ tiêu tự nhiên và nhân sinh (độ dốc, độ cao DEM, hướng sườn, thạch học, mật độ đứt gãy lineament, điều kiện ĐCCT, ĐCTV, vỏ phong hóa, chỉ số thực vật NDVI, mật độ dòng chảy, loại đất, lượng mưa cực đại), tính toán chỉ số nhất quán ngẫu nhiên ($CR < 0.1$).
  • Mô hình phần tử hữu hạn GeoStudio 2012 (SEEP/W & CTRAN/W): Giải bài toán dòng thấm bão hòa - không bão hòa 2 chiều kết hợp phương trình vi phân truyền chất bình suy - phân tán (advection-dispersion equation), mô phỏng bước phát tán chất ô nhiễm tại bãi thải quặng thiếc Bản Cô và mỏ đá vôi Châu Hồng qua các mốc 250, 500, 750, 1000, 1250, 1500, 1750 và 2000 ngày.

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

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

      So Sánh Mức Độ Ô Nhiễm Nguồn Nước & Đất Tại Vùng Mỏ Ven Biển Thạch Hà
  1. Thực trạng suy thoái và ô nhiễm nguồn nước mặt - nước ngầm nghiêm trọng:
    • Tại lưu vực suối Chùa Tran (dài 8km, cấp nước 200 hộ dân và 820ha đất nông nghiệp các xã Thạch Đỉnh, Thạch Khê, Thạch Hải), nước có màu trắng đục, nồng độ pH từ 4,86 đến 5,94 (thấp hơn quy chuẩn QCVN 08-MT:2015 1,01-1,23 lần), COD dao động 19-47 mg/l (vượt quy chuẩn 1,9-4,7 lần), $\text{BOD}_5$ từ 13-26 mg/l (vượt 3,25-6,5 lần). Nước ngầm giếng sinh hoạt ô nhiễm đặc biệt nghiêm trọng với hàm lượng COD 28-56 mg/l, vượt quy chuẩn QCVN 09-MT:2015 từ 7 đến 14 lần do nước thải các mỏ đá Thạch Hải, Công ty 999 và mỏ sắt Thạch Khê xả trực tiếp.
    • Tại suối Đồng Bàn (xã Thạch Văn, cấp nước 50 hộ dân, 195ha đất lúa), nước mặt có pH 5,18-5,52; COD 22-58 mg/l (vượt 2,2-5,8 lần); $\text{BOD}_5$ 12-23 mg/l (vượt 3-5,8 lần). Nước ngầm nhiễm bẩn nặng nề với COD đạt ngưỡng 21-70 mg/l, vượt quy chuẩn từ 5,2 đến 17,5 lần từ nguồn xả thải của mỏ khai thác quặng titan sa khoáng.
  2. Ô nhiễm kim loại nặng trong đất vượt ngưỡng sinh thái:
    • Khu vực Thạch Khê (900ha đất canh tác ven chân núi Nam Giới): Hàm lượng Asen (As) dao động 15,5-29,7 ppm (vượt QCVN 03-MT:2015 từ 1,29 đến 2,48 lần), $\text{Cr}^{6+}$ đạt 368 ppm (vượt 1,47 lần).
    • Khu vực Thạch Văn (400ha đất canh tác): Hàm lượng As từ 15,0 đến 18,1 ppm (vượt 1,25-1,51 lần), $\text{Cr}^{6+}$ đạt 317 ppm (vượt 1,27 lần) do tích tụ quặng mịn và nước thải đuôi quặng titan.
  3. Hiện trạng bức xạ gamma và khí phóng xạ tự nhiên:
    • Kết quả đo 94 điểm gamma môi trường và 96 điểm Radon tại mỏ đá granit núi Nam Giới ghi nhận suất liều gamma từ 0,06 đến 0,43 $\mu\text{Sv/h}$, trong đó có 13 điểm ghi nhận suất liều $\ge 0,3\ \mu\text{Sv/h}$ (vượt tiêu chuẩn an toàn thứ cấp, tập trung cục bộ tại diện tích bóc đá gốc). Nồng độ khí Radon trong đất ghi nhận biến thiên 96,5-1620 $\text{Bq/m}^3$, trong không khí 0-54,4 $\text{Bq/m}^3$ và trong nước 32-182 $\text{Bq/m}^3$.
  4. Cơ chế lan truyền chất ô nhiễm theo mô hình phần tử hữu hạn CTRAN/W:
    • Kết quả mô phỏng động lực học phát tán chất ô nhiễm tại bãi thải quặng thiếc sa khoáng Bản Cô (xã Châu Thành) và mỏ đá vôi Châu Hồng (Quỳ Hợp) chứng minh: Quầng ô nhiễm kim loại nặng mở rộng bán kính thâm nhập sâu vào tầng chứa nước ngầm theo gradient thủy lực sau 500-1000 ngày, và đạt trạng thái lan truyền ổn định diện rộng sau 1500-2000 ngày, đe dọa trực tiếp các bồn trũng canh tác lúa hạ lưu.
  5. Độ chính xác vượt trội của mạng nơ-ron nhân tạo (ANN) trong dự báo trượt lở:
    • Kết quả đối chuẩn ma trận sai số khẳng định mô hình ANN đạt diện tích dưới đường cong AUC cao nhất ($>88%$), theo sau là mô hình Bayes WoE ($>83%$) và mô hình Hệ số tin cậy CF ($>81%$), chỉ ra các điểm nóng nguy cơ trượt tầng cao nhất tập trung tại các xã Châu Hồng, Châu Lộc và bản Na Kỳ (xã Liên Hợp).

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp bằng chứng thực nghiệm khẳng định sự tương tác phi tuyến giữa tải trọng bãi thải mỏ và trường ứng suất địa chất; củng cố khung lý thuyết lan truyền chất ô nhiễm đa môi trường.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình tích hợp liên thông giữa GIS, ảnh viễn thám đa thời kỳ phân giải cao và mô hình địa kỹ thuật phần tử hữu hạn, có khả năng chuyển giao áp dụng cho toàn bộ các tỉnh vùng mỏ Tây Bắc và Tây Nguyên.
  • Về mặt ứng dụng thực tiễn: Cung cấp bộ giải pháp công trình chi tiết gồm: thiết kế góc dốc bờ moong an toàn $\le 45^\circ-60^\circ$ theo cấu trúc thạch học, xây dựng hệ thống hào thoát nước đón đỉnh sườn bãi thải, ứng dụng sơ đồ nổ mìn vi sai phi điện (Non-electric delay blasting) giảm chấn động nứt nẻ đá vôi, và gia cố đập hồ lắng quặng đuôi bằng vải địa kỹ thuật kết hợp bẫy bùn sinh thái.
  • Về mặt chính sách quản lý nhà nước: Đề xuất cơ chế thẩm định báo cáo ĐTM bắt buộc phải tích hợp bản đồ số hóa dự báo lan truyền ô nhiễm và nguy cơ vỡ đập quặng đuôi; bổ sung quy định ký quỹ cải tạo phục hồi môi trường mỏ lũy tiến theo mức độ nhạy cảm sinh thái địa phương.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn chỉ ra 04 giới hạn khoa học mang tính khách quan:

  1. Độ phân giải của dữ liệu viễn thám: Nguồn ảnh Landsat (độ phân giải không gian 30m) tuy tối ưu cho phân tích chuỗi thời gian lịch sử nhưng chưa thể hiện được các vết nứt vi mô ($\le 1-2\text{m}$) trên từng bờ moong khai thác cục bộ.
  2. Dữ liệu quan trắc khí tượng cực đoan: Số lượng trạm quan trắc lượng mưa tự động tại vùng núi cao Quỳ Hợp còn thưa, việc nội suy trường mưa $24\text{h}$ dựa vào mô hình không gian có thể xuất hiện sai số cục bộ tại các thung lũng karst kín.
  3. Giới hạn không gian mô phỏng thấm 2D: Mô hình GeoStudio CTRAN/W được áp dụng trên các lát cắt địa chất 2D đại diện, chưa phản ánh đầy đủ hiệu ứng dòng thấm 3D phức tạp trong cấu trúc hang caxtơ ngầm phát triển bất đối xứng.
  4. Chuỗi số liệu quan trắc kim loại nặng theo mùa: Việc lấy mẫu môi trường nước và đất tập trung vào các đợt quan trắc định kỳ chính, chưa phản ánh liên tục 24/7 sự biến thiên nồng độ ô nhiễm tức thời trong các đợt lũ quét cực đoan.

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

  • Ứng dụng công nghệ ảnh viễn thám siêu cao (Sentinel-2 10m, ảnh UAV Flycam LiDAR) và kỹ thuật giao thoa sóng radar viễn thám (InSAR) để quan trắc biến dạng sườn tầng mỏ với độ chính xác milimét theo thời gian thực.
  • Mở rộng mô hình lan truyền chất ô nhiễm từ 2D sang không gian 3D ngầm đa pha (FEFLOW, MODFLOW-USG) tích hợp mạng lưới dòng chảy karst ngầm.
  • Phát triển hệ thống cảnh báo sớm tai biến thời gian thực (Real-time Early Warning System - EWS) dựa trên nền tảng cảm biến IoT kết hợp thuật toán học sâu (Deep Learning).

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

  • Tác động học thuật (Academic Impact): Công trình là tài liệu tham khảo mẫu mực trong đào tạo tiến sĩ và sau đại học các chuyên ngành Kỹ thuật Môi trường, Địa kỹ thuật công trình, Địa chất kinh tế tại Đại học Mỏ - Địa chất, Viện Khoa học & Công nghệ Môi trường (ĐHBK Hà Nội) và Đại học Quốc gia Hà Nội; mở ra hướng công bố trên các tạp chí ISI/Scopus Q1 về địa kỹ thuật và môi trường mỏ.
  • Tác động chuyển đổi ngành khai khoáng (Industry Transformation): Định hình lại quy chuẩn an toàn bãi thải và hồ lắng quặng đuôi cho các doanh nghiệp khoáng sản (Công ty CP Sắt Thạch Khê, Tổng Công ty Khoáng sản & Thương mại Hà Tĩnh, các hiệp hội đá hoa Quỳ Hợp), thúc đẩy chuyển dịch sang công nghệ khai thác mỏ xanh và kinh tế tuần hoàn chất thải mỏ.
  • Tác động quản trị chính sách (Policy Influence): Cung cấp luận cứ khoa học thực tiễn giúp Sở Tài nguyên & Môi trường Nghệ An và Hà Tĩnh rà soát, đình chỉ các cơ sở chế biến xả thải trái quy chuẩn; trực tiếp phục vụ công tác khoanh định các khu vực cấm và tạm thời cấm hoạt động khoáng sản theo Luật Khoáng sản.
  • Lợi ích kinh tế - xã hội: Bảo vệ an toàn nguồn nước sinh hoạt cho hơn 250 hộ dân hạ lưu, ngăn ngừa nguy cơ hoang mạc hóa và ô nhiễm kim loại nặng cho hơn 1.300 ha đất canh tác nông nghiệp tại các lưu vực suối Chùa Tran, Đồng Bàn và Quỳ Hợp.

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

  • Nghiên cứu sinh & Nhà khoa học: Tiếp cận khung phân tích phương pháp luận hoàn chỉnh tích hợp giữa toán học xác suất, mô hình phần tử hữu hạn và công nghệ viễn thám địa không gian.
  • Giảng viên & Viện nghiên cứu: Sử dụng bộ cơ sở dữ liệu thực địa chuẩn hóa gồm 36 bảng số liệu và 66 sơ đồ/bản đồ chuyên đề làm tài liệu giảng dạy chuyên sâu về Quản lý rủi ro môi trường.
  • Doanh nghiệp khai khoáng: Nắm bắt quy chuẩn kỹ thuật an toàn nổ mìn, góc dốc moong khai thác và quy trình gia cố đập quặng đuôi nhằm triệt tiêu rủi ro sự cố công trình gây thiệt hại kinh tế.
  • Cơ quan quản lý nhà nước (Bộ TN&MT, UBND cấp Tỉnh/Huyện): Sở hữu công cụ số hóa phục vụ công tác thẩm định báo cáo đánh giá tác động môi trường (ĐTM), cấp phép khai thác và kiểm tra thanh tra sau cấp phép.

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?

Luận án đã mở rộng Lý thuyết Đánh giá Rủi ro Tai biến Tự nhiên - Nhân sinh (Coupled Human-Natural Hazards Theory) của IPCC AR5 và JTC1 bằng cách định lượng hóa mối quan hệ phi tuyến giữa các yếu tố nền địa chất (vỏ phong hóa, đứt gãy kiến tạo) với các tác nhân kích hoạt nhân sinh mỏ (góc dốc cắt tầng, độ rung chấn nổ mìn, tải trọng đống thải quặng đuôi). Đột phá nằm ở việc tích hợp thành công phương trình thấm - tán xạ phi bão hòa qua hàm lực hút dính ma trận đất đá thải vào bài toán dự báo tai biến môi trường diện rộng.

2. Tính đổi mới trong phương pháp nghiên cứu khi so sánh với ít nhất 2 nghiên cứu quốc tế trước đây?

So với nghiên cứu của Lan và cộng sự (2004) tại Trung Quốc chỉ dùng mô hình Bayes WoE tĩnh cho sạt lở tự nhiên và nghiên cứu của Chalkias và cộng sự (2014) tại Hy Lạp dùng chỉ số LSI dựa trên dữ liệu địa hình - lượng mưa, luận án đổi mới vượt bậc bằng cách: (1) Kết hợp so sánh chéo đồng thời 04 mô hình (WoE, CF, LR, ANN) để loại bỏ sai số chủ quan; (2) Tích hợp chuỗi ảnh viễn thám đa thời gian (1990-2017) nhận diện động thái bóc phủ moong mỏ trái phép; (3) Sử dụng mô hình số phần tử hữu hạn GeoStudio 2012 CTRAN/W để mô phỏng động học lan truyền chất ô nhiễm từ 250 đến 2000 ngày.

3. Phát hiện bất ngờ nhất (counter-intuitive/surprising finding) có bằng chứng thực nghiệm là gì?

Phát hiện bất ngờ nhất là tại khu vực mỏ sắt Thạch Khê và mỏ sa khoáng titan Thạch Văn: dù địa hình ven biển tương đối bằng phẳng (độ dốc $<5^\circ$, năng lượng địa hình thấp), nhưng mức độ ô nhiễm nước ngầm và đất nông nghiệp lại đạt ngưỡng đặc biệt nghiêm trọng (COD nước ngầm vượt tới 17,5 lần, Asen vượt 2,48 lần, $\text{Cr}^{6+}$ đạt 368 ppm). Điều này phản ánh rằng trong địa hình đồng bằng cát ven biển, cơ chế thẩm thấu trọng lực và lan truyền hóa chất trong tầng chứa nước bở rời nguy hiểm hơn gấp nhiều lần so với các tai biến động lực hình thái sườn dốc vùng núi.

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

Có. Luận án thiết lập một giao thức chuẩn hóa 5 bước nghiêm ngặt: (1) Thu thập và chuẩn hóa cơ sở dữ liệu GIS địa không gian 12 lớp chỉ tiêu; (2) Xử lý ảnh viễn thám Landsat bằng thuật toán MLC và SAM với sai số nắn chỉnh $<1\text{ pixel}$; (3) Thiết lập ma trận trọng số chuyên gia Saaty AHP kết hợp học máy ANN; (4) Khảo sát thực địa đo đạc mẫu khí, bụi, nước, đất bằng thiết bị chuyên dụng chuẩn Mỹ (GrayWolf, 5959P, Extech); (5) Chạy mô phỏng phần tử hữu hạn CTRAN/W với các thông số đặc trưng đất nước SWCC đo đạc thực nghiệm.

5. Kế hoạch nghiên cứu 10 năm (10-year research agenda) được vạch ra như thế nào?

Kế hoạch 10 năm tập trung vào: (1) Xây dựng mạng lưới trạm quan trắc địa kỹ thuật - hóa môi trường tự động IoT tại 100% các mỏ quặng đuôi trọng điểm; (2) Phát triển nền tảng bản đồ số Digital Twin mô phỏng tai biến địa môi trường mỏ thời gian thực cho toàn lãnh thổ Việt Nam; (3) Nghiên cứu sâu về độc học sinh thái học (Ecotoxicology) sự chuyển hóa kim loại nặng vào chuỗi thức ăn nông nghiệp tại các vùng mỏ hoàn nguyên.


Kết luận

  1. Thành lập cơ sở dữ liệu địa môi trường mỏ toàn diện: Xây dựng thành công hệ thống dữ liệu không gian tích hợp 12 lớp chỉ tiêu tự nhiên, công nghệ và hiện trạng khai thác cho 02 vùng trọng điểm Quỳ Hợp (Nghệ An) và Thạch Hà - Cẩm Xuyên (Hà Tĩnh).
  2. Lượng hóa nguy cơ và kiểm chuẩn mô hình khoa học: Chứng minh mô hình mạng nơ-ron nhân tạo (ANN) và mô hình thống kê xác suất Bayes (WoE) mang lại độ chính xác cao nhất (AUC $>83-88%$) trong việc phân vùng dự báo nhạy cảm trượt lở sườn tầng và lũ bùn đá do vỡ đập quặng đuôi.
  3. Làm sáng tỏ động học phát tán ô nhiễm: Xác lập thành công mô hình phần tử hữu hạn CTRAN/W dự báo bán kính lan truyền chất ô nhiễm từ bãi thải quặng thiếc và đá vôi qua các mốc thời gian từ 250 đến 2000 ngày, cung cấp bằng chứng thực nghiệm về mức độ ô nhiễm nước ngầm (COD vượt tới 17,5 lần) và đất (As vượt 2,48 lần).
  4. Hệ thống giải pháp giảm thiểu đa tầng: Đề xuất đồng bộ 03 nhóm giải pháp công trình (hào đón nước, đập lắng, nổ mìn vi sai phi điện), phi công trình (quan trắc GIS viễn thám tự động) và cơ chế chính sách ký quỹ cải tạo môi trường mỏ.
  5. Mở ra các nhánh nghiên cứu tiên phong: Đặt nền móng cho các hướng nghiên cứu ứng dụng AI trong địa kỹ thuật mỏ, mô phỏng 3D dòng chảy karst ngầm và giám sát biến dạng vi mô mỏ bằng công nghệ vệ tinh Radar InSAR tại Việt Nam.