Tổng quan về luận án

Nghiên cứu về tính chất luyện kim của tài nguyên khoáng sản sắt đóng vai trò then chốt trong việc tối ưu hóa hiệu suất nhiệt động học và động học của các dây chuyền sản xuất gang lò cao (Blast Furnace - BF) cũng như công nghệ hoàn nguyên trực tiếp sản xuất sắt xốp (Direct Reduced Iron - DRI). Tại Việt Nam, tỉnh Cao Bằng sở hữu nguồn tài nguyên quặng sắt phong phú với tổng trữ lượng và tài nguyên dự báo vượt trên 44 triệu tấn phân bố tại 15 mỏ và biểu hiện khoáng sản (như Nà Rụa, Ngườm Cháng, Nà Lũng, Tà Phình, Nà Mè - Hào Lịch). Trong đó, mỏ sắt Nà Rụa là mỏ lớn nhất với trữ lượng đạt trên 18 triệu tấn, chiếm khoảng 40% tổng trữ lượng quặng sắt toàn tỉnh, đặc trưng bởi hàm lượng sắt cao dao động từ 55% đến 70% Fe. Tuy nhiên, phần lớn quặng sắt Nà Rụa khai thác dạng quặng vụn, đặc sít và có tính chất luyện kim phức tạp chưa từng được khảo sát toàn diện, dẫn đến hàng loạt sự cố công nghệ tại 06 cơ sở luyện gang lò cao và nhà máy sản xuất sắt xốp MIREX trên địa bàn tỉnh (điển hình là hiện tượng sụt liệu, tiêu hao nhiên liệu cao, độ hoàn nguyên không đồng đều và tạo xỉ dính bết).

Luận án tiến sĩ kỹ thuật vật liệu của nghiên cứu sinh Nguyễn Thanh Hoàn với đề tài "Nghiên cứu tính chất luyện kim và khả năng ứng dụng quặng sắt Cao Bằng" (Mã số: 9520309, Trường Đại học Bách khoa Hà Nội, dưới sự hướng dẫn khoa học của PGS. Nguyễn Sơn Lâm và PGS. Bùi Anh Hòa) là công trình khoa học tiên phong tại Việt Nam khảo cứu một cách có hệ thống, chuẩn mực và chuyên sâu về bản chất vật lý, hóa học, đặc tính biến mềm và động học hoàn nguyên của quặng sắt Nà Rụa.

flowchart TD
    A["Quặng sắt Nà Rụa Cao Bằng (Fe: 55-70%)"] --> B["Nghiên cứu chế tạo Quặng vê viên"]
    B --> C["Chế độ nung ô-xy hóa P500 (500°C)"]
    B --> D["Chế độ nung ô-xy hóa P1200 (1200°C)"]
    C --> E["Khảo sát động học hoàn nguyên với CO/Ar (60%, 80%, 100% CO)"]
    D --> E
    C --> F["Khảo sát động học hoàn nguyên với Các-bon rắn (7%, 10%, 15%, 20% C)"]
    D --> F
    E --> G["Mô hình hóa động học (Crank-Ginstling-Brounshtein & Shrinking Core)"]
    F --> G
    G --> H["Xác định Hằng số tốc độ k & Năng lượng hoạt hóa Ea"]
    H --> I["Định hướng công nghệ: Lò cao BF & Sản xuất sắt xốp DRI (MIREX)"]

Khoảng trống nghiên cứu (research gap) lớn nhất trong ngành luyện kim trong nước trước công trình này là sự thiếu hụt các dữ liệu thực nghiệm định lượng về:

  • Động học chuyển pha tuần tự $Fe_2O_3 \rightarrow Fe_3O_4 \rightarrow FeO \rightarrow Fe$ của quặng viên Nà Rụa trong môi trường khí khử CO và các-bon rắn ở nhiệt độ cao (900 – 1100 °C);
  • Cơ chế khuếch tán nội tại qua lớp sản phẩm xốp so với trở kháng phản ứng bề mặt;
  • Tương quan giữa chế độ nung ô-xy hóa quặng viên (nung ở 500 °C - P500 và nung ở 1200 °C - P1200) đến thể tích lỗ xốp, diện tích bề mặt riêng ($S_{BET}$), bán kính lỗ xốp, độ bền cơ học và tính hoàn nguyên.

Để giải quyết triệt để khoảng trống này, luận án thiết lập 3 câu hỏi nghiên cứu và giả thuyết tương ứng:

  1. Câu hỏi nghiên cứu 1 (RQ1): Chế độ nhiệt độ nung ô-xy hóa (500 °C so với 1200 °C) ảnh hưởng như thế nào đến sự chuyển biến cấu trúc pha, mật độ lỗ xốp vi mô và độ bền nén của quặng viên Nà Rụa?
    Giả thuyết 1 (H1): Quặng viên nung ở 1200 °C (P1200) tái kết tinh mạng ê-ma-tít ($\alpha\text{-}Fe_2O_3$) bền vững kèm theo sự nóng chảy pha silicat tạo cầu liên kết, làm giảm diện tích bề mặt riêng nhưng nâng cao vượt bậc độ bền nén cơ học so với mẫu nung 500 °C (P500).
  2. Câu hỏi nghiên cứu 2 (RQ2): Cơ chế nào kiểm soát tốc độ phản ứng hoàn nguyên quặng viên Nà Rụa bằng hỗn hợp khí CO-Ar và các-bon rắn trong khoảng nhiệt độ 900 – 1100 °C?
    Giả thuyết 2 (H2): Do quặng sắt Nà Rụa có độ xốp ban đầu thấp, quá trình hoàn nguyên quặng viên P1200 tuân theo mô hình khuếch tán qua lớp sản phẩm Crank–Ginstling–Brounshtein thay vì phản ứng bề mặt thuần túy.
  3. Câu hỏi nghiên cứu 3 (RQ3): Đặc tính biến mềm và giãn nở nhiệt của quặng sắt Nà Rụa chi phối như thế nào đến khả năng lưu thông khí động học trong cột vật liệu lò cao và lò quay sản xuất sắt xốp?
    Giả thuyết 3 (H3): Hàm lượng các-bon phối liệu làm thay đổi đáng kể khoảng biến mềm nhiệt độ cao, quyết định vị trí vùng dẻo (cohesive zone) trong lò cao.

Khung lý thuyết của nghiên cứu tích hợp lý thuyết nhiệt động học luyện kim (giản đồ Richardson-Jeffes, giản đồ cân bằng Baur-Glaessner, giản đồ pha Fe-O Darken-Gurry) và lý thuyết động học phản ứng dị thể khí - rắn. Nghiên cứu thực hiện trên mẫu quặng sắt nguyên khai Nà Rụa, các mẫu quặng vê viên P500, P1200, các phối liệu quặng - than (M1: 0% than, M2: 7% than, M3: 10% than, M4: 15% than, M5: 20% than) với giải nhiệt độ thử nghiệm lên đến 1300 °C, đóng góp dữ liệu định lượng chính xác phục vụ thiết kế và vận hành lò luyện kim công nghiệp.

Literature Review và Positioning

Nghiên cứu về quá trình hoàn nguyên ô-xít sắt đã trải qua hơn một thế kỷ phát triển với các trường phái lý thuyết kinh điển:

  • Trường phái nhiệt động học cân bằng pha: Đặt nền móng bởi Baur và Glaessner (1903), Richardson và Jeffes (1948, 1949), Darken và Gurry (1945, 1946). Các tác giả đã xây dựng hệ thống giản đồ thế ô-xy, xác định áp suất riêng phần ô-xy cân bằng $pO_2$ và tỷ lệ khí cân bằng $CO/CO_2$, $H_2/H_2O$ qua các phản ứng chuyển hóa ô-xít sắt.
  • Trường phái động học phản ứng dị thể: Phát triển bởi Crank (1956), Ginstling và Brounshtein (1950), Szekely, Evans và Sohn (1976), nghiên cứu sâu về các phương trình truyền khối, khuếch tán nội hạt và mô hình lõi co không phản ứng (Shrinking Core Model - SCM).
  • Trường phái công nghệ hoàn nguyên trực tiếp: Khởi xướng bởi Eketorp (1989) với các sơ đồ hoàn nguyên pha lỏng và pha rắn nhằm loại bỏ bước luyện gang gián tiếp qua lò cao.

Trong y văn thế giới tồn tại hai luồng quan điểm tranh luận lớn:

  1. Cơ chế kiểm soát tốc độ hoàn nguyên quặng viên đặc chắc: Nascimento và cộng sự (1997) [63] khi khảo sát quặng hematite Brazil cho rằng phản ứng biên pha tại ranh giới ô-xít/kim loại tuân theo mô hình $1 - (1-f)^{1/3} = kt$ là giai đoạn khống chế tốc độ chính. Ngược lại, các công trình của Turkdogan (1980) và Szekely (1976) chứng minh rằng đối với quặng viên nung kết có độ đặc sít cao, trở kháng khuếch tán của khí khử qua lớp vỏ sắt kim loại xốp ($1 - \frac{2}{3}f - (1-f)^{2/3} = kt$) mới là yếu tố quyết định tốc độ tổng thể.
  2. Hiệu ứng chuyển pha mạng tinh thể đến độ hoàn nguyên: Một số nghiên cứu quốc tế chỉ ra rằng quá trình chuyển biến từ $\alpha\text{-}Fe_2O_3$ (mạng lục giác xếp chặt) sang $Fe_3O_4$ (mạng lập phương spinel) làm tăng thể tích riêng dẫn đến nứt nẻ tế vi, qua đó thúc đẩy tốc độ hoàn nguyên nhanh hơn so với việc hoàn nguyên trực tiếp từ quặng manhetit nguyên sinh.

So với các công trình quốc tế như nghiên cứu của Pineau và cộng sự (2006) về động học hoàn nguyên quặng manhetit hạt mịn bằng hỗn hợp khí $H_2\text{-}CO$, hay nghiên cứu của Spreitzer và Schenk (2019) về sự suy giảm tính hoàn nguyên do tạo pha fayalite ($2FeO\cdot SiO_2$), luận án của Nguyễn Thanh Hoàn định vị rõ nét đối tượng nghiên cứu là loại quặng sắt skarn Nà Rụa có hàm lượng sắt cao nhưng cấu trúc hạt nguyên sinh đặc sít, chứa hàm lượng tạp chất $SiO_2, Al_2O_3$ nhất định. Luận án đã làm sáng tỏ việc lựa chọn chế độ nung ô-xy hóa chuyển manhetit thành ematit trước khi hoàn nguyên giúp tăng cường tính luyện kim của quặng viên Cao Bằng một cách vượt bậc.

Đó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 và hoàn thiện lý thuyết động học hoàn nguyên dị thể áp dụng riêng cho quặng sắt đặc hữu Việt Nam:

  1. Kiểm chứng và chuẩn hóa mô hình động học: Luận án đã chứng minh một cách thuyết phục rằng quá trình hoàn nguyên quặng viên Nà Rụa nung ở 1200 °C (mẫu P1200) bằng khí CO và các-bon rắn tuân theo mô hình Crank–Ginstling–Brounshtein [62] với phương trình động học: $$1 - \frac{2}{3}f - (1-f)^{2/3} = kt$$ trong đó $f$ là mức độ hoàn nguyên, $k$ là hằng số tốc độ phản ứng ($s^{-1}$), $t$ là thời gian. Điều này bác bỏ giả thiết cho rằng phản ứng biên pha là yếu tố kiểm soát duy nhất đối với quặng viên nung kết sâu.
  2. Xác định các thông số động học nhiệt động định lượng: Luận án xác lập chính xác năng lượng hoạt hóa ($E_a$) và hằng số tốc độ ($k$) trong từng điều kiện nồng độ khí CO (60%, 80%, 100%) và các-bon rắn, cung cấp bộ hằng số thực nghiệm quan trọng cho cơ sở dữ liệu nhiệt động học vật liệu luyện kim.
  3. Minh chứng quy luật chuyển pha topochemical: Sử dụng phương pháp nhiễu xạ rơnghen (XRD), luận án khẳng định cơ chế chuyển biến cấu trúc tuần tự xảy ra theo 3 giai đoạn rõ rệt ở nhiệt độ $\ge 570\text{ }^\circ\text{C}$: $$Fe_2O_3 \rightarrow Fe_3O_4 \rightarrow Fe_{1-x}O \text{ (vu-tít)} \rightarrow Fe$$ đặc biệt xác định pha vu-tít phi hợp thức ($Fe_{1-x}O$) là giai đoạn có trở kháng nhiệt động và động học lớn nhất trong toàn bộ chuỗi phản ứng.
graph LR
    subgraph "Chuỗi chuyển pha Topochemical (T >= 570°C)"
        H["Fe2O3 (Ê-ma-tít)<br/>Mạng lục giác AB"] -->|"Tách Oxy (1/6 lk)"| M["Fe3O4 (Ma-nhê-tít)<br/>Mạng lập phương Spinel ABC"]
        M -->|"Tách Oxy (1/3 lk)"| W["Fe1-xO (Vu-tít)<br/>Mạng lập phương NaCl"]
        W -->|"Giai đoạn trở kháng cao nhất (1/2 lk)"| Fe["Fe (Sắt kim loại)<br/>Mạng lập phương tâm khối"]
    end

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

Khung phân tích của luận án tích hợp đa chiều các trụ cột lý thuyết:

  • Lý thuyết nhiệt động học chiết ly kim loại: Đánh giá thế hóa học của ô-xy, quan hệ phân ly ô-xít theo giản đồ Richardson-Jeffes và cân bằng phản ứng Boudouard: $$2CO \rightleftharpoons CO_2 + C$$
  • Lý thuyết tương thích mặt mạng tinh thể: Phân tích sự định hướng ưu tiên và sai lệch mạng khi tái cấu trúc từ lớp anion $O^{2-}$ (bán kính $0,136\text{ nm}$) sang cation $Fe^{2+}$ ($0,078\text{ nm}$) và $Fe^{3+}$ ($0,065\text{ nm}$).
  • Điều kiện biên nhiệt động học và truyền khối: Xác lập mối liên hệ giữa hệ số truyền khối $\beta = D/\delta$ (với $D$ là hệ số khuếch tán theo định luật Fick, $\delta$ là chiều dày lớp biên) và áp suất riêng phần khí $p_{CO}, p_{CO_2}$.

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 quan điểm thực chứng (positivism) với phương pháp định lượng thực nghiệm nghiêm ngặt. Cấu trúc nghiên cứu phân tầng từ cấu trúc tế vi hạt khoáng (micro-level), đặc tính hạt quặng vê viên (meso-level) đến các chỉ tiêu luyện kim vĩ mô (macro-level).

graph TD
    A["Thiết kế nghiên cứu đa tầng"] --> B["Cấp độ Vi mô (Micro-level)"]
    A --> C["Cấp độ Trung gian (Meso-level)"]
    A --> D["Cấp độ Vĩ mô (Macro-level)"]
    
    B --> B1["Mạng tinh thể, Anion O2- (0.136 nm), Cation Fe2+/Fe3+"]
    B --> B2["Phổ XRD Rietveld, SEM mặt vỡ, EDX Line-scan"]
    
    C --> C1["Diện tích bề mặt BET (0.199 m2/g)"]
    C --> C2["Thể tích lỗ xốp (0.001 cm3/g), Bán kính lỗ xốp (16.872 Å)"]
    C --> C3["Độ bền nén cơ học quặng viên"]
    
    D --> D1["Mức độ hoàn nguyên R (%) theo thời gian"]
    D --> D2["Kính hiển vi nhiệt độ cao: Giãn nở & Biến mềm (M1-M5)"]
    D --> D3["Động học Crank-Ginstling-Brounshtein & Năng lượng hoạt hóa Ea"]

Đối tượng nghiên cứu là mẫu quặng sắt lấy trực tiếp từ mỏ Nà Rụa, tỉnh Cao Bằng. Quặng được nghiền mịn đạt cấp hạt tiêu chuẩn luyện kim, phối trộn chất kết dính bentonite và vê viên trên đĩa tạo viên với đường kính danh định $10 \div 15\text{ mm}$.

Quy trình nghiên cứu rigorous

Quy trình thực nghiệm được chuẩn hóa qua các bước kiểm soát nghiêm ngặt:

  1. Chế tạo và thiêu kết mẫu quặng viên: Mẫu quặng vê viên sống được sấy khô và nung ô-xy hóa trong lò điện ở các mức nhiệt độ: 500 °C, 900 °C, 1000 °C, 1100 °C và 1200 °C với thời gian lưu nhiệt 2 giờ.
  2. Khảo sát tính chất nhiệt lý và biến mềm: Sử dụng thiết bị kính hiển vi nhiệt độ cao đo độ co rút diện tích và xác định các điểm nhiệt độ bắt đầu biến mềm ($t_1$), biến mềm rõ rệt ($t_2$) và chảy dẻo ($t_3$) của các mẫu quặng nguyên khai và quặng phối trộn 7%, 10%, 15%, 20% than.
  3. Thực nghiệm hoàn nguyên bằng khí và các-bon:
    • Hoàn nguyên bằng hỗn hợp khí $CO\text{-}Ar$: Thực hiện trong lò ống nằm ngang có kiểm soát lưu lượng khí chính xác ở các nồng độ 60% CO + 40% Ar, 80% CO + 20% Ar và 100% CO tại các mốc nhiệt độ 900 °C, 950 °C, 1000 °C, 1050 °C, 1100 °C trong các khoảng thời gian 20, 30, 40, 60 và 90 phút.
    • Hoàn nguyên bằng các-bon rắn: Mẫu quặng viên được chôn trong lớp than khử rắn, gia nhiệt đẳng nhiệt tại các mức nhiệt độ tương ứng, đo độ giảm khối lượng để xác định tỷ phần hoàn nguyên $f$ và mức độ hoàn nguyên $R%$.
sequenceDiagram
    autonumber
    actor Tech as Kỹ sư / Nhà nghiên cứu
    participant Pellet as Máy vê viên đĩa & Lò nung
    participant Phys as Máy đo Hấp phụ Autosorb iq
    participant React as Lò phản ứng hoàn nguyên
    participant Analytic as Hệ thống XRD / SEM-EDX
    
    Tech->>Pellet: Phối trộn quặng Nà Rụa + Bentonite, nung P500 & P1200 (2h)
    Pellet->>Phys: Cung cấp mẫu viên nung kết
    Phys-->>Tech: Trả kết quả: SBET = 0.199 m2/g, Vp = 0.001 cm3/g, r = 16.872 Å
    Tech->>React: Nạp mẫu P500 / P1200, cấp khí CO-Ar (60-100%) hoặc Than (7-20%)
    Note over React: Hoàn nguyên ở 900 - 1100°C theo các mốc thời gian 20 - 90 phút
    React->>Analytic: Chuyển mẫu sau hoàn nguyên để phân tích
    Analytic-->>Tech: Dữ liệu phổ XRD (pha Fe2O3, Fe3O4, FeO, Fe), ảnh SEM, phổ EDX

Data và phân tích

Luận án đã sử dụng tổ hợp các phương pháp và thiết bị phân tích hiện đại:

  • Phân tích diện tích bề mặt riêng và độ xốp: Thiết bị hấp phụ khí tự động Quantachrome Autosorb iq, đo diện tích bề mặt theo phương pháp BET, phân bố thể tích và bán kính lỗ xốp vi mô qua đường đẳng nhiệt giải hấp BJH.
  • Phân tích cấu trúc pha: Thiết bị nhiễu xạ tia X (XRD) với nguồn phát Cu-$K\alpha$, nhận diện định danh các pha khoáng vật qua ngân hàng dữ liệu PDF (Powder Diffraction File).
  • Phân tích tổ chức tế vi và thành phần nguyên tố vi vùng: Kính hiển vi điện tử quét (SEM) phân giải cao kết hợp phổ tán xạ năng lượng tia Rơnghen (EDX) dạng điểm và đường (line-scan).
  • Phân tích nhiệt vi sai và nhiệt trọng lượng: Thiết bị DTA-TG ghi nhận biến thiên entanpy và tổn hao khối lượng theo chương trình gia nhiệt liên tục đến 1300 °C.

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

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

Luận án đã đạt được 5 phát hiện khoa học mang tính đột phá với đầy đủ bằng chứng thực nghiệm định lượng:

[!NOTE] Trích dẫn minh chứng thực nghiệm 1: "Các chỉ số về diện tích bề mặt, thể tích lỗ xốp và bán kính của các lỗ xốp của mẫu viên quặng sau khi được nung ở nhiệt độ 1200 °C được kiểm tra bằng máy đo Autosorb iq cho kết quả: Diện tích bề mặt = $0,199\text{ m}^2/\text{g}$, Thể tích lỗ xốp = $0,001\text{ cm}^3/\text{g}$, Bán kính lỗ xốp ($r$) = $16,872\text{ \AA}$."

[!NOTE] Trích dẫn minh chứng thực nghiệm 2: "Giản đồ phân tích rơnghen mẫu quặng viên P500 được hoàn nguyên bằng các-bon rắn tại các giá trị mức độ hoàn nguyên $R$ là (a) 13%, (b) 47%, (c) 64% và (d) 96%... khẳng định cơ chế hoàn nguyên là $\text{ê-ma-tít} \rightarrow \text{ma-nhê-tít} \rightarrow \text{vu-tít} \rightarrow \text{sắt nguyên chất}$."

[!NOTE] Trích dẫn minh chứng thực nghiệm 3: "Theo công bố của Nascimento [63], động học của phản ứng hoàn nguyên từ vu-tít về sắt tuân theo mô hình lõi co của phản ứng biên pha, tuân theo mô hình toán học $1 - (1 - f)^{1/3} = kt$. Quá trình hoàn nguyên tổng thể tuân theo mô hình động học Crank–Ginstling–Brounshtein [62] quá trình hoàn nguyên khuếch tán qua lớp sản phẩm, theo phương trình $1 - \frac{2}{3}f - (1 - f)^{2/3} = kt$."

  1. Quy luật tiến hóa vi cấu trúc theo nhiệt độ nung: Quặng viên nung ở 500 °C (P500) giữ cấu trúc xốp rời rạc, độ bền nén thấp. Khi nung nâng nhiệt lên 1200 °C (P1200), quá trình kết tinh lại của $\alpha\text{-}Fe_2O_3$ diễn ra hoàn toàn, tạo nên khung liên kết tinh thể vững chắc với độ bền nén vượt trội, mặc dù thể tích lỗ xốp giảm chỉ còn $0,001\text{ cm}^3/\text{g}$.
  2. Sự hình thành pha sắt kim loại phân tầng: Kết quả ảnh SEM và phổ EDX line-scan chứng minh quá trình tạo mầm sắt kim loại diễn ra từ ngoài vào trong. Tại mức độ hoàn nguyên $R = 39%$, pha sắt bắt đầu kết tụ ở viền hạt; đến $R = 60\div 64%$, lớp sắt tạo thành màng liên tục bao bọc lõi vu-tít; và ở $R \ge 95%$, mẫu chuyển biến thành cấu trúc sắt xốp đồng nhất với mạng lỗ rỗng phân bố đều.
  3. Sự biến thiên hằng số tốc độ phản ứng ($k$) và năng lượng hoạt hóa ($E_a$): Khi tăng hàm lượng CO từ 60% lên 100%, hằng số tốc độ phản ứng $k$ tăng tuyến tính rõ rệt. Năng lượng hoạt hóa $E_a$ tính toán từ mô hình Crank–Ginstling–Brounshtein cho thấy quá trình hoàn nguyên ở nhiệt độ cao bị chi phối mạnh bởi khuếch tán thể tích của khí qua lớp vỏ sắt đặc.
  4. Quy luật biến mềm của quặng sắt phối liệu than: Phép đo hiển vi nhiệt độ cao cho thấy quặng nguyên khai M1 có nhiệt độ biến mềm cao, nhưng khi phối trộn các-bon rắn (mẫu M2: 7% than đến M5: 20% than), sự tạo pha cục bộ làm khoảng nhiệt độ biến mềm mở rộng, thay đổi bản chất tạo xỉ ban đầu.
  5. Hiệu suất hoàn nguyên vượt bậc của các-bon rắn ở nhiệt độ cao: Ở nhiệt độ 1050 – 1100 °C trong thời gian 90 phút, mức độ hoàn nguyên của quặng viên P1200 và P500 khi có mặt các-bon rắn đều đạt trên 95%, khẳng định phản ứng khí hóa Boudouard ($C + CO_2 \rightarrow 2CO$) đóng vai trò động lực tái sinh khí khử cực mạnh tại chỗ.
Mẫu quặng viên & Chế độ nung $S_{BET}$ ($\text{m}^2/\text{g}$) Thể tích lỗ xốp ($\text{cm}^3/\text{g}$) Bán kính lỗ xốp $r$ ($\text{\AA}$) Cơ chế động học chủ đạo Độ hoàn nguyên $R_{max}$ (1100°C, 90 min)
P500 (Nung 500 °C, 2h) Cao hơn ($>0,5$) Lớn ($>0,008$) Lớn ($>25$) Hỗn hợp: Phản ứng biên pha & Khuếch tán $96%$ (với Than) / $92%$ (với CO 100%)
P1200 (Nung 1200 °C, 2h) $0,199$ $0,001$ $16,872$ Crank–Ginstling–Brounshtein (Khuếch tán qua lớp vỏ Fe) $95%$ (với Than) / $89%$ (với CO 100%)

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp bằng chứng thực nghiệm chính xác xác lập ranh giới chuyển tiếp giữa mô hình phản ứng biên pha hạt nhân co và mô hình khuếch tán màng xỉ/kim loại Crank-Ginstling-Brounshtein trong điều kiện quặng manhetit biến tính nung kết.
  • Về mặt phương pháp luận: Xây dựng quy trình tích hợp giữa hiển vi nhiệt độ cao, hấp phụ đẳng nhiệt Autosorb iq, XRD Rietveld và xử lý động học thực nghiệm có thể áp dụng chuẩn hóa cho các mỏ quặng sắt khác tại Việt Nam.
  • Về ứng dụng thực tiễn:
    • Đối với sản xuất sắt xốp (DRI - Nhà máy MIREX): Khuyến cáo bắt buộc phải duy trì nhiệt độ vùng hoàn nguyên tối thiểu từ $1050\text{ }^\circ\text{C}$ đến $1100\text{ }^\circ\text{C}$, tỷ lệ than phối trộn tối ưu $15 \div 20%$ để đạt mức độ kim loại hóa $>95%$, triệt tiêu hiện tượng vón cục dính bết do pha wustite dư.
    • Đối với luyện gang lò cao (BF): Quặng sắt Nà Rụa cần được vê viên và nung ô-xy hóa ở $1200\text{ }^\circ\text{C}$ (P1200) trước khi nạp lò nhằm nâng cao độ bền cơ học nguội, giảm hiện tượng vỡ vụn sinh bụi và tối ưu hóa độ thông thoáng khí trong thân lò.
  • Về mặt chính sách: Cung cấp cơ sở khoa học xác đáng cho Sở Công Thương Cao Bằng và Bộ Công Thương trong việc thẩm định công nghệ, quy hoạch vùng nguyên liệu và cấm xuất khẩu quặng sắt thô chưa qua chế biến sâu.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án thẳng thắn ghi nhận 4 hạn chế nghiên cứu mang tính biên điều kiện:

  1. Thực nghiệm hoàn nguyên bằng khí $CO\text{-}Ar$ được tiến hành trong điều kiện đẳng nhiệt ở quy mô phòng thí nghiệm, chưa mô phỏng đầy đủ gradien nhiệt độ phi đẳng nhiệt và áp suất động học biến đổi dọc theo chiều cao cột liệu lò cao công nghiệp.
  2. Nguồn khí khử thực nghiệm chủ yếu là CO tinh khiết pha Ar; chưa khảo sát sâu ảnh hưởng cạnh tranh của khí Hydro ($H_2$) và hơi nước ($H_2O$) vốn luôn hiện diện trong khí than khí hóa công nghiệp.
  3. Luận án chưa đánh giá chi tiết ảnh hưởng của các hợp chất kiềm ($K_2O, Na_2O$) và kẽm ($Zn$) có trong quặng sắt skarn đến hiện tượng trương nở thể tích bất thường (catastrophic swelling) của viên quặng trong môi trường khử.
  4. Thử nghiệm biến mềm mới thực hiện thông qua đo quang học không tải; chưa tích hợp thiết bị đo độ biến mềm dưới tải trọng cơ học (RUL - Reduction Under Load) để mô phỏng chính xác vùng dẻo lò cao.

Từ các hạn chế trên, 5 hướng nghiên cứu tiếp theo (future research agenda) được đề xuất:

  1. Nghiên cứu động học hoàn nguyên quặng viên Nà Rụa bằng nguồn khí giàu hydro ($H_2\text{-DRI}$) định hướng luyện kim xanh giảm phát thải carbon ($CO_2$).
  2. Thiết lập hệ thống đo độ biến mềm co ngót dưới tải trọng (RUL) với các chế độ tải trọng từ 0,5 đến $2,0\text{ kg/cm}^2$.
  3. Khảo sát động học tạo pha xỉ nhiệt độ cao giữa quặng viên Nà Rụa với các trợ dung đá vôi, đolomit địa phương để tối ưu hóa tỷ lệ bazơ $(CaO + MgO)/(SiO_2 + Al_2O_3)$.
  4. Mô phỏng động lực học chất lưu (CFD) tương tác khí - rắn 3D trong lò quay sản xuất sắt xốp dựa trên các hằng số động học $k, E_a$ đã thu được.
  5. Nghiên cứu công nghệ đóng bánh nguội quặng sắt - than composite (Carbon-Composite Agglomerates) sử dụng chất kết dính hữu cơ để hạ nhiệt độ hoàn nguyên.

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

  • Tác động học thuật: Luận án tạo lập một công trình tham chiếu tiêu chuẩn cho ngành kỹ thuật vật liệu luyện kim tại Việt Nam. Các bài báo khoa học công bố từ luận án trên các tạp chí chuyên ngành uy tín cung cấp bộ dữ liệu động học chuẩn cho các nhà nghiên cứu trong nước và quốc tế khi nghiên cứu về quặng sắt khu vực Đông Nam Á.
  • Tác động đến công nghiệp luyện kim: Giúp các doanh nghiệp luyện kim tại Cao Bằng (như Công ty Cổ phần Khoáng sản Xây dựng 30-4, Công ty Cổ phần Chế biến Khoáng sản Cao Sơn Hà, Nhà máy Sắt xốp MIREX) điều chỉnh chính xác chế độ nhiệt, thành phần phối liệu than - quặng, giúp giảm tiêu hao than cốc từ 5 - 8%, nâng cao năng suất lò cao và chất lượng sắt xốp thương phẩm.
  • Tác động chính sách và xã hội: Góp phần hiện thực hóa chiến lược phát triển công nghiệp khai khoáng chế biến sâu của tỉnh Cao Bằng, chuyển đổi từ xuất khẩu quặng thô sang sản xuất sản phẩm giá trị gia tăng cao (gang đúc, phôi thép), bảo đảm an ninh nguyên liệu và tạo việc làm bền vững cho lao động địa phương.
  • Ý nghĩa quốc tế: Đóng góp vào bức tranh khoa học toàn cầu về hành vi nhiệt động và động học của các mỏ quặng sắt skarn loại giàu magnetite khi ứng dụng các công nghệ hoàn nguyên trực tiếp thế hệ mới.

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

  • Nghiên cứu sinh và học viên cao học ngành Kỹ thuật Vật liệu/Luyện kim: Tiếp cận một cấu trúc nghiên cứu mẫu mực kết hợp chặt chẽ giữa mô hình hóa lý thuyết và thực nghiệm phân tích hiện đại.
  • Các nhà khoa học và giảng viên luyện kim: Có được nguồn dữ liệu định lượng đã kiểm chứng về hằng số tốc độ $k$, năng lượng hoạt hóa $E_a$, đặc tính giãn nở nhiệt và giản đồ chuyển pha để đưa vào bài giảng và công trình nghiên cứu chuyên sâu.
  • Kỹ sư vận hành và R&D tại các nhà máy luyện gang thép: Nắm vững thông số công nghệ cốt lõi: nhiệt độ nung quặng viên ($1200\text{ }^\circ\text{C}$), nhiệt độ hoàn nguyên sắt xốp ($1050 \div 1100\text{ }^\circ\text{C}$), nồng độ CO tối thiểu để tránh tích tụ pha wustite gây hại cho lò.
  • Các nhà quản lý tài nguyên và hoạch định chính sách: Sử dụng các luận cứ khoa học của luận án để quy hoạch định mức khai thác, kiểm soát tiêu chuẩn kỹ thuật nguyên liệu nạp lò tại các khu liên hợp luyện kim.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc xác định chính xác mô hình kiểm soát tốc độ hoàn nguyên quặng viên manhetit Cao Bằng nung kết sâu (P1200) tuân theo mô hình khuếch tán Crank–Ginstling–Brounshtein ($1 - \frac{2}{3}f - (1-f)^{2/3} = kt$). Luận án đã mở rộng lý thuyết truyền khối dị thể trong môi trường vật liệu xốp thấp ($V_p = 0,001\text{ cm}^3/\text{g}$), chỉ ra rằng sự tái kết tinh $\alpha\text{-}Fe_2O_3$ ở $1200\text{ }^\circ\text{C}$ tạo ra mạng tinh thể chặt chẽ khiến trở kháng khuếch tán qua lớp vỏ sắt kim loại mới sinh trở thành yếu tố quyết định tốc độ phản ứng ở nhiệt độ cao ($\ge 1000\text{ }^\circ\text{C}$).

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

So với các nghiên cứu trước đây tại Việt Nam (chủ yếu dừng lại ở phân tích hóa học tổng lượng và thử nghiệm nung tĩnh đơn giản) và nghiên cứu của Nascimento (1997) [63] (chỉ xét phản ứng biên pha trên mẫu hạt nhỏ), phương pháp luận của luận án có tính đột phá nhờ:

  • Tích hợp đồng thời kỹ thuật hấp phụ khí đẳng nhiệt Autosorb iq để định lượng diện tích bề mặt $S_{BET}$ và bán kính lỗ xốp vi mô;
  • Sử dụng kính hiển vi quang học nhiệt độ cao quan sát liên tục biến thiên diện tích mẫu để xác định chính xác các điểm biến mềm nhiệt độ cao ($t_1, t_2, t_3$);
  • Định lượng sự biến thiên pha tinh thể theo mức độ hoàn nguyên ($R = 13%, 47%, 64%, 96%$) qua phổ XRD Rietveld kết hợp phân tích vi vùng SEM-EDX.

3. Phát hiện bất ngờ nhất (counter-intuitive) có sự hỗ trợ của dữ liệu là gì?

Phát hiện bất ngờ nhất là nghịch lý giữa độ xốp và khả năng hoàn nguyên cực đại của mẫu nung kết $1200\text{ }^\circ\text{C}$ (P1200). Mặc dù mẫu P1200 có thể tích lỗ xốp ban đầu cực thấp ($0,001\text{ cm}^3/\text{g}$) so với mẫu P500, nhưng ở nhiệt độ cao ($1100\text{ }^\circ\text{C}$) với các-bon rắn, mẫu P1200 vẫn đạt mức độ hoàn nguyên vượt trội ($R = 95%$). Dữ liệu SEM chỉ ra rằng chính sự chênh lệch thể tích riêng giữa pha mẹ $\alpha\text{-}Fe_2O_3$ và pha con $Fe_3O_4$ trong quá trình chuyển biến mạng tinh thể đã tự sinh ra hệ thống vi nứt rạn (micro-cracks) cục bộ, tạo đường dẫn nội hạt cho khí khử CO khuếch tán sâu vào tâm viên quặng.

4. Luận án có cung cấp quy trình nhân bản (replication protocol) chi tiết không?

Có. Luận án cung cấp đầy đủ và chi tiết toàn bộ quy trình chuẩn bị mẫu và thực nghiệm: thành phần phối liệu quặng mịn Nà Rụa với bentonite; thông số vê viên đĩa; chế độ nung ô-xy hóa P500 và P1200 (nâng nhiệt, đẳng nhiệt 2h, làm nguội); lưu lượng khí $CO/Ar$ (tỷ lệ 60/40, 80/20, 100/0); tỷ lệ phối trộn than M1 - M5 (0%, 7%, 10%, 15%, 20%); và quy trình phân tích XRD/SEM-EDX/Autosorb iq, cho phép tái lập nghiên cứu với độ chính xác cao.

5. Chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?

Chương trình 10 năm tới tập trung chuyển dịch từ luyện kim hoàn nguyên truyền thống sang luyện kim phi carbon:

  • Giai đoạn 1 (1-3 năm): Nghiên cứu hoàn nguyên quặng sắt Cao Bằng bằng khí hydro tinh khiết ($H_2$) và khí tổng hợp syngas sinh khối;
  • Giai đoạn 2 (4-6 năm): Nghiên cứu thử nghiệm thiết bị RUL xác lập bản đồ vùng dẻo lò cao cho quặng phối liệu đa mỏ;
  • Giai đoạn 3 (7-10 năm): Chế tạo thử nghiệm sắt xốp chất lượng siêu cao phục vụ nấu luyện thép hợp kim đặc chủng và vật liệu bột từ tính nano tại Việt Nam.

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Thanh Hoàn là một công trình nghiên cứu khoa học xuất sắc, mẫu mực và toàn diện về tính chất luyện kim của quặng sắt Nà Rụa - Cao Bằng, mang lại 6 đóng góp cốt lõi:

  1. Xác lập hệ thống dữ liệu toàn diện: Cung cấp bức tranh khoa học hoàn chỉnh đầu tiên về đặc tính khoáng vật học, nhiệt lý, diện tích bề mặt vi mô ($S_{BET} = 0,199\text{ m}^2/\text{g}$), thể tích lỗ xốp ($0,001\text{ cm}^3/\text{g}$) và độ bền nén của quặng vê viên Nà Rụa.
  2. Minh chứng bản chất chuyển biến cấu trúc pha: Khẳng định quy luật chuyển pha topochemical 3 giai đoạn nghiêm ngặt: $$Fe_2O_3 \rightarrow Fe_3O_4 \rightarrow Fe_{1-x}O \rightarrow Fe$$ qua định lượng phổ XRD tại các mức độ hoàn nguyên từ $13%$ đến $96%$.
  3. Xác lập mô hình động học chính xác: Chứng minh quá trình hoàn nguyên quặng viên nung kết P1200 tuân theo mô hình Crank–Ginstling–Brounshtein ($1 - \frac{2}{3}f - (1-f)^{2/3} = kt$), giải quyết dứt điểm các tranh luận lý thuyết về cơ chế kiểm soát tốc độ đối với quặng đặc chắc skarn.
  4. Làm chủ quy luật biến mềm nhiệt độ cao: Khảo sát tường tận ảnh hưởng của hàm lượng các-bon phối liệu đến các điểm nhiệt độ bắt đầu biến mềm và chảy dẻo, cung cấp cơ sở kiểm soát khí động học lò cao.
  5. Định hình giải pháp công nghệ công nghiệp đột phá: Đưa ra thông số tối ưu cho quy trình sản xuất sắt xốp DRI tại lò quay MIREX ($1050 \div 1100\text{ }^\circ\text{C}$, than $15 \div 20%$) và quy trình sản xuất quặng viên nung kết $1200\text{ }^\circ\text{C}$ nạp lò cao BF.
  6. Mở ra các hướng nghiên cứu công nghệ xanh tương lai: Đặt nền móng khoa học vững chắc để phát triển các công nghệ hoàn nguyên trực tiếp bằng hydro ($H_2\text{-DRI}$), đóng góp thiết thực vào mục tiêu phát triển bền vững và giảm phát thải ròng của ngành công nghiệp luyện kim Việt Nam.