Bối cảnh và vấn đề nghiên cứu

Trong các nhà máy công nghiệp hiện đại, hệ thống cung cấp năng lượng nhiệt và lò hơi (steam generator / boiler) đóng vai trò trung tâm trong việc đảm bảo nguồn hơi ổn định cho các dây chuyền sản xuất. Việc vận hành các thiết bị trao đổi nhiệt, buồng đốt công nghiệp (Fired Heater / Combustion Chamber) đòi hỏi sự kiểm soát nghiêm ngặt về mặt an toàn, hiệu suất chuyển hóa năng lượng và chất lượng hơi sinh ra. Tuy nhiên, quá trình trao đổi nhiệt và quá trình cháy trong lò hơi mang bản chất phi tuyến phức tạp, có sự liên kết chặt chẽ giữa các thông số nhiệt động lực học, lưu lượng, áp suất và thành phần hóa học của khí thải.

Đề tài "Tối ưu hóa quy trình công nghiệp với công nghệ mới" được thực hiện nhằm giải quyết bài toán mô hình hóa, tính toán truyền nhiệt và xây dựng cấu trúc điều khiển tự động cho hệ thống nhiệt công nghiệp, gắn liền với khảo sát thực tế tại doanh nghiệp (hợp tác kỹ thuật và thiết bị liên quan đến đối tác Thermax India Ltd). Đồ án tập trung vào các nhiệm vụ nghiên cứu cụ thể:

  1. Thiết lập cơ sở lý thuyết về truyền nhiệt bức xạ, cân bằng nhiệt động và phản ứng cháy của nhiên liệu trong buồng đốt công nghiệp.
  2. Xây dựng mô hình toán học mô tả động học biến thiên áp suất và dao động mức nước bao hơi dưới tác động của sự thay đổi tải tiêu thụ hơi và công suất nhiệt.
  3. Phân tích, tích hợp các giải pháp điều khiển tự động hóa quá trình công nghiệp sử dụng hệ thống điều khiển phân tán (DCS), bộ điều khiển logic khả trình (PLC) cùng các vòng lặp điều khiển phản hồi PID/PI.
  4. Đánh giá ứng dụng công nghệ mô phỏng hiện đại, chuyển dịch từ các công cụ tính toán truyền thống sang phần mềm mô phỏng chuyên dụng Aspen Hysys kết hợp cùng Matlab và Simulink Toolboxes để tối ưu hóa quy trình vận hành nhiệt.

Đối tượng nghiên cứu của đồ án là các thiết bị trao đổi nhiệt, buồng đốt công nghiệp, lò sinh hơi và các hệ thống phụ trợ (quạt cấp gió FD, quạt hút ID, bộ hâm nước Economiser, bộ quá nhiệt Superheater, bộ bay hơi Evaporator). Phạm vi nghiên cứu bao gồm việc phân tích lý thuyết nhiệt động, mô hình hóa động học quá trình sinh hơi và cấu trúc các mạch điều khiển tự động được triển khai trong môi trường sản xuất công nghiệp thực tế kéo dài qua các phân đoạn công nghệ đến trang 99 của văn bản.

Cơ sở lý thuyết và phương pháp

Đồ án xây dựng nền tảng nghiên cứu dựa trên các định luật nhiệt động lực học, cơ học chất lưu và lý thuyết điều khiển tự động công nghiệp:

  • Lý thuyết trao đổi nhiệt bức xạ trong buồng lửa: Tác giả sử dụng định luật Bouguer để mô tả độ đen của lớp khí bức xạ: $$a_g = 1 - e^{-k_g p s}$$ Trong đó $k_g$ là hệ số hấp thụ của khí, $p$ là áp suất thành phần, và $s$ [m] là chiều dày hiệu dụng của lớp khí. Mật độ dòng bức xạ nhiệt từ ngọn lửa truyền đến bề mặt vách hấp thụ nhiệt được tính theo phương trình: $$q_{bx} = c_o a_{bl} \left(\frac{T_l}{100}\right)^4 \times 10^{-3} \quad [\text{kW/m}^2]$$ Tiêu chuẩn Boltzmann ($Bo$) được áp dụng để xác định trường nhiệt độ khói đầu ra của buồng đốt: $$Bo = \frac{B_{tt} Q_{bx}}{c_o a_{bl} \psi F_v T_a^4 \times 10^{-11}}$$ Nhiệt độ khói ra khỏi buồng lửa $t''{bl}$ [$^\circ\text{C}$] được xác định qua biểu thức phụ thuộc vào tiêu chuẩn $Bo$, hệ số hấp thụ $a{bl}$, bề mặt vách buồng lửa $F_v$ và nhiệt dung riêng thể tích trung bình của dòng khói $(Vc)_{tb}$.
  • Cơ chế hóa học và cân bằng nhiệt quá trình cháy: Quá trình cháy được phân tích dựa trên thành phần làm việc của nhiên liệu gồm cacbon ($C^{ar}$ hoặc $C^{lv}$), hydro ($H^{ar}$), oxy ($O^{ar}$), nitơ ($N^{ar}$), lưu huỳnh ($S^{ar}$), tro ($A^{ar}$), độ ẩm ($M^{ar}$ / $W^{lv}$) và nhiệt trị làm việc $Q^{ar}$ ($Q^{lv}$). Lượng oxy lý thuyết $V_{O2}^0$ và không khí khô $V_{kko}$ [$m^3_{tc}/\text{kg}$ nhiên liệu] được tính toán từ các phản ứng cháy hóa học tiêu chuẩn.
  • Mô hình động học lò hơi: Động học tích lũy nhiệt và áp suất được xây dựng thông qua hệ phương trình vi phân phi tuyến bậc nhất, thiết lập mối quan hệ giữa biến thiên khối lượng nước - hơi và biến thiên năng lượng tích tụ trong kim loại lò.

Phương pháp nghiên cứu kết hợp giữa giải tích giải phương trình vi phân truyền nhiệt, lập bảng đối chuẩn thực nghiệm trên các dòng lò hơi công nghiệp đã lắp đặt, và ứng dụng công cụ mô phỏng số trên phần mềm Aspen Hysys cùng Matlab/Simulink. Nguồn dữ liệu bao gồm các công thức nhiệt kỹ thuật, số liệu thực nghiệm từ các hệ lò công nghiệp (KPS-6, DKB, BKZ, PK), và số liệu vận hành thực tế tại doanh nghiệp.

Nội dung chính theo từng chương

Chương 1: Tổng quan hệ thống nhiệt công nghiệp và lò hơi

Chương mở đầu trình bày sơ đồ công nghệ tổng thể (Process Flow Diagram - PFD) của hệ thống buồng đốt công nghiệp (Fired Heater / Combustion Chamber) và lò sinh hơi. Nội dung phân tích cấu tạo, vai trò chức năng của các phân hệ chính:

  • Buồng đốt (Fired Heater / Combustion Chamber): Không gian diễn ra phản ứng cháy giữa nhiên liệu (than, gas, nhiên liệu lỏng) và không khí, giải phóng nhiệt lượng ở dải nhiệt độ cao lên tới $2200^\circ\text{C}$.
  • Hệ thống bề mặt trao đổi nhiệt: Gồm bộ hâm nước (Economiser) tận dụng nhiệt thừa khói thải để gia nhiệt nước cấp; bộ bay hơi (Evaporator) chuyển hóa nước thành hỗn hợp hơi bão hòa; và bộ quá nhiệt (Superheater) nâng nhiệt độ hơi lên các cấp nhiệt $360^\circ\text{C}$, $500^\circ\text{C}$, $540^\circ\text{C}$ đáp ứng tiêu chuẩn vận hành tuabin hoặc cấp nhiệt công nghệ.
  • Hệ thống cấp gió và thải khói: Quạt gió cưỡng bức (Force Draught Fan - FD Fan) cấp không khí sơ cấp và thứ cấp; quạt hút khói (Induced Draught Fan - ID Fan) duy trì áp suất âm trong buồng đốt, đảm bảo lưu thông dòng khí khói qua các dãy ống đối lưu và thải qua ống khói.

Chương 2: Quá trình nhiệt động, trao đổi nhiệt bức xạ và cơ chế cháy nhiên liệu

Nội dung chương 2 đi sâu vào tính toán cân bằng năng lượng và các hiện tượng hóa lý diễn ra trong buồng đốt:

  • Nhiệt lượng hữu ích và cân bằng năng lượng: Cân bằng nhiệt buồng lửa $Q_{bl}$ [$kJ/kg$] được thiết lập dựa trên nhiệt trị nhiên liệu và các thành phần tổn thất nhiệt: $$Q_{bl} = (Q_{lv}^{tr} + Q_{kk} + Q_{kkng} - Q_{tth}) \frac{100 - q_3 - q_4 - q_6}{100}$$ Trong đó $q_3, q_4, q_6$ lần lượt là các tổn thất nhiệt do cháy không hoàn toàn về mặt hóa học, cháy không hoàn toàn về mặt cơ học (tro xỉ mang nhiệt và sót than chưa cháy), và tổn thất nhiệt tỏa ra môi trường xung quanh: $$q_3 = \frac{Q_3}{Q_{dv}} \times 100%; \quad q_4 = \frac{Q_4}{Q_{dv}} \times 100%; \quad q_6 = \frac{Q_6}{Q_{dv}} \times 100%$$
  • Nhiệt độ bắt cháy của nhiên liệu: Văn bản tổng hợp nhiệt độ tự bốc cháy của các loại nhiên liệu thông dụng: than bùn ($225^\circ\text{C}$), than nâu ($300 - 350^\circ\text{C}$), than hoa ($360^\circ\text{C}$), và khí methane $\text{CH}_4$ ($700^\circ\text{C}$).
  • Các phản ứng cháy và tỏa nhiệt lượng: $$C + O_2 \to CO_2 + 8.084 \quad [\text{kcal/kg Cacbon}]$$ $$2C + O_2 \to 2CO + 2.430 \quad [\text{kcal/kg Cacbon}]$$ $$CO + 0,5O_2 \to CO_2 + Q_{C2}$$ $$2H_2 + O_2 \to 2H_2O + 28.922 \quad [\text{kcal/kg } H_2]$$ $$S + O_2 \to SO_2 + 2.224 \quad [\text{kcal/kg S}]$$
  • Hệ số không khí thừa và kiểm soát khí thải: Phân tích sự phụ thuộc giữa nồng độ $\text{CO}_2$, $\text{O}_2$, $\text{CO}$ trong khói thải với lượng không khí thừa (Excess air từ 0% đến 100%). Chế độ cháy tối ưu đạt được khi duy trì $\text{CO}_2$ ở dải $14,5 - 15%$ và oxy thừa $\text{O}_2$ ở mức $2 - 3%$. Ngoài ra, chương nghiên cứu nhiệt độ điểm sương axit của $H_2SO_4$ và nguy cơ ăn mòn bề mặt truyền nhiệt do hợp chất vanadi $V_2O_5$ tại dải nhiệt độ $550 - 580^\circ\text{C}$ (đặc biệt khi vượt trên $550^\circ\text{C}$).

Chương 3: Xây dựng mô hình toán học và động học quá trình nhiệt trong lò hơi

Chương 3 thiết lập hệ thống phương trình động học mô tả các quá trình quá độ trong bao hơi và giàn ống sinh hơi:

  • Phương trình vi phân cân bằng vật chất: $$\frac{d}{d\tau}(V_n \rho' + V_h \rho'') = D_{nc} - D \quad [\text{kg/s}]$$ Trong đó $V_n, V_h$ là thể tích phần nước và phần hơi trong bao hơi; $\rho', \rho''$ là khối lượng riêng của nước sôi và hơi bão hòa; $D_{nc}$ là lưu lượng nước cấp và $D$ là lưu lượng hơi tiêu thụ.
  • Phương trình vi phân cân bằng năng lượng: $$\frac{d}{d\tau}(V_n \rho' i' + V_h \rho'' i'' + G_{kl} C_{kl} t_{kl}) = D_{nc} i_{nc} - D i'' + Q \quad [\text{kW}]$$ Với $i_{nc}, i', i''$ lần lượt là enthalpy của nước cấp, nước sôi và hơi bão hòa; $G_{kl}, C_{kl}, t_{kl}$ là khối lượng, nhiệt dung riêng và nhiệt độ của kim loại vách tiếp nhiệt; $Q$ là nhiệt lượng buồng đốt cung cấp.
  • Mô hình biến thiên áp suất và nhiệt tích lũy: Biến thiên áp suất $\frac{dp}{d\tau}$ được biểu diễn qua phương trình: $$\frac{dp}{d\tau} = \frac{\Delta i_{hn} D_{nc} + Q - \varepsilon_2 D}{\varepsilon_3 V_n + \varepsilon_4 V_h + \varepsilon_5 G_{kl}} \quad [\text{MN/m}^2\text{s}]$$ Trong đó đại lượng tích tụ nhiệt của hệ thống là $Q_T = \varepsilon_3 V_n + \varepsilon_4 V_h + \varepsilon_5 G_{kl}$ [$\text{kJ/MN/m}^2$].
  • Động học mức nước bao hơi (Drum level dynamics): Tốc độ biến thiên mức nước $\frac{dh}{d\tau}$ [m/s] chịu ảnh hưởng của hiệu ứng bọt hơi ngập ("nước sôi giả"): $$\frac{dh}{d\tau} = a_1 Q + a_2 D_{nc} - a_3 D + \frac{1}{F_{bh}} \frac{dV_{ho}}{d\tau}$$ Biến thiên thể tích bọt hơi ngập $V_{ho}$ được mô tả theo phụ tải tương đối $D/D_0$ và tốc độ thay đổi áp suất $dP/dP_0$: $$\frac{dV_{ho}}{d\tau} = \frac{V_{ho}}{D} \frac{dD}{d\tau} - \frac{V_{ho}}{p} \frac{dp}{d\tau}$$

Dưới đây là bảng đối chuẩn thông số tích lũy động học của một số chủng loại lò hơi công nghiệp được tác giả tổng hợp:

Loại lò hơi Áp suất định mức $p$ [$\text{MPa}$] Sản lượng hơi định mức $D$ [t/h] Tốc độ tích lũy áp suất [$\text{MPa/s}$]
KPS-6 1,27 6,0 0,0042
DKB 1,27 6,5 0,0027
BKZ 4,32 75,0 0,0230
PK 10,80 230,0 0,0435
Lò trung gian 3,34 180,0 0,0104

Chương 4: Hệ thống điều khiển tự động và mô phỏng trên Aspen Hysys, Matlab/Simulink

Chương 4 hoàn thiện giải pháp kỹ thuật bằng việc thiết kế cấu trúc điều khiển và xây dựng mô hình mô phỏng:

  • Triển khai cấu trúc điều khiển phân tán DCS kết hợp PLC, bao gồm các trạm vận hành và trạm kỹ thuật (Engineering Work Station - EWS), hệ thống ghi nhận dữ liệu (Data Logger System - DLS), và hệ thống điều khiển toàn nhà máy (Plant Control System - PCS).
  • Cấu hình các bộ điều khiển tích phân tỷ lệ - vi phân (PID/PI) tương ứng với các biến quá trình (Process Variable - PV), giá trị đặt (Setpoint - SP), và tín hiệu đầu ra điều khiển (Output - OP):
    • Bộ điều khiển áp suất hơi (Pressure Integral Control - PIC).
    • Bộ điều khiển nhiệt độ hơi quá nhiệt (Temperature Integral Control - TIC) kết hợp van phun nước làm nguội (Quench Water).
    • Bộ điều khiển lưu lượng nước cấp (Flow Integral Control - FIC) theo thuật toán điều khiển mức nước bao hơi 3 phần tử (lưu lượng nước cấp, mức nước, lưu lượng hơi).
  • Sử dụng phần mềm Aspen Hysys kết hợp Matlab / Simulink Toolboxes để giải các phương trình nhiệt động phức tạp, thay thế phương pháp giải tích rời rạc thủ công, cho phép kiểm chứng đáp ứng quá độ của toàn hệ thống kéo dài qua chuỗi hình vẽ từ Hình 4-1 đến trang 99.

Thiết kế và triển khai hệ thống kỹ thuật

Hệ thống điều khiển và mô phỏng kỹ thuật được thiết kế với kiến trúc đồng bộ từ tầng đo lường cảm biến đến tầng giám sát điều khiển:

[Biến quá trình: PV] ---> [Bộ điều khiển: PIC / TIC / FIC (PID/PI)] ---> [Cơ cấu chấp hành: OP]
                                     |
                                     v
                        [Trạm kỹ thuật: DCS / PLC / PCS]
                                     |
                        [Mô phỏng: Aspen Hysys + Matlab/Simulink]
  • Kiến trúc công nghệ điều khiển: Hệ thống sử dụng mạng điều khiển phân tán DCS kết nối các trạm PLC chuyên trách. Trạm kỹ thuật EWS và hệ thống lưu trữ dữ liệu DLS thu thập liên tục các giá trị PV (nhiệt độ vách ống, nhiệt độ hơi, áp suất buồng đốt, mức nước bao hơi) và so sánh với giá trị đặt SP để xuất tín hiệu điều khiển OP tới các van điều khiển nước cấp, quạt gió FD Fan, quạt hút ID Fan.
  • Quy cách ống trao đổi nhiệt: Các thông số kết cấu được khảo sát và áp dụng bao gồm ống đường kính ngoài $51\text{ mm}$, độ dày vách $1,5\text{ mm}$, bước ống bố trí khoảng cách $120\text{ mm}$ và $38\text{ mm}$. Đường ống nhánh dẫn hơi có đường kính tiêu chuẩn $2,54\text{ cm}$ (1 inch), với các tuyến ống chính chạy xa trên $6\text{ m}$.
  • Môi trường phần mềm tích hợp: Việc tích hợp Aspen Hysys với Matlab/Simulink cho phép mô phỏng chính xác trạng thái nhiệt động bão hòa của nước - hơi, enthalpy hỗn hợp ở các dải nhiệt độ $360^\circ\text{C} - 540^\circ\text{C}$ và giải các phương trình vi phân biến thiên áp suất và mức nước trong thời gian thực.

Kết quả và đóng góp

  1. Thiết lập mô hình toán học hoàn chỉnh cho quá trình nhiệt: Đồ án đã mô hình hóa thành công mối quan hệ giữa các thông số động học của lò hơi thông qua hệ phương trình vi phân cân bằng vật chất (3-1), cân bằng năng lượng (3-2), và phương trình biến thiên áp suất (3-11).
  2. Xác định các thông số vận hành buồng đốt tối ưu: Phân tích định lượng quá trình cháy cho thấy chế độ vận hành tối ưu của buồng đốt cần duy trì nồng độ $\text{CO}_2$ trong khói thải ở mức $14,5 - 15%$ và nồng độ $\text{O}_2$ dư ở mức $2 - 3%$ nhằm cực tiểu hóa các tổn thất hóa học ($q_3$) và cơ học ($q_4$), đồng thời kiểm soát nhiệt độ tiếp nhiệt để tránh hiện tượng ăn mòn do $V_2O_5$ tại nhiệt độ trên $550^\circ\text{C}$ ($(550 - 580)^\circ\text{C}$).
  3. Hiện đại hóa phương pháp mô phỏng quy trình công nghiệp: Đóng góp chính của tác giả là việc ứng dụng thành công phần mềm chuyên dụng Aspen Hysys và Matlab/Simulink Toolbox thay thế cho các phương pháp tính toán gần đúng truyền thống, hỗ trợ công tác thiết kế các mạch vòng điều khiển PID (FIC, TIC, PIC) trên nền tảng DCS và PLC.

Hạn chế và hướng nghiên cứu tiếp

  • Hạn chế: Các số liệu thực nghiệm và hệ số tích lũy động học của mô hình được đối chuẩn chủ yếu trên các dòng lò hơi công nghiệp tiêu chuẩn (KPS-6, DKB, BKZ, PK) và các thiết bị đặc thù của Thermax India Ltd; chưa bao quát toàn bộ các cấu hình lò hơi tầng sôi tuần hoàn hoặc lò hơi siêu tới hạn phức tạp.
  • Hướng nghiên cứu tiếp: Mở rộng nghiên cứu thuật toán điều khiển thích nghi, điều khiển dự báo mô hình (MPC) trên nền tảng DCS hiện đại, tích hợp sâu hơn giữa mô hình dòng chảy đa pha trong ống trao đổi nhiệt với dữ liệu thời gian thực thu thập từ hệ thống DLS.

Giá trị tham khảo

Đồ án là tài liệu tham khảo hữu ích cho sinh viên và kỹ sư thuộc các ngành:

  • Kỹ thuật Nhiệt / Năng lượng: Tham khảo chi tiết phương pháp tính toán cân bằng nhiệt buồng lửa, truyền nhiệt bức xạ (định luật Bouguer, tiêu chuẩn Boltzmann), và cơ chế phản ứng cháy nhiên liệu rắn/khí.
  • Tự động hóa / Kỹ thuật Điều khiển: Tham khảo cấu trúc thiết kế vòng lặp điều khiển mức nước bao hơi 3 phần tử, điều khiển áp suất và nhiệt độ hơi (FIC, TIC, PIC), kiến trúc hệ thống DCS/PLC.
  • Kỹ thuật Hóa học / Quá trình: Tham khảo quy trình xây dựng mô hình công nghệ PFD và ứng dụng phần mềm Aspen Hysys liên kết với Matlab/Simulink để mô phỏng động thái hệ thống nhiệt.

Câu hỏi thường gặp

1. Tiêu chuẩn Boltzmann ($Bo$) trong đồ án được sử dụng nhằm mục đích gì?
Tiêu chuẩn Boltzmann ($Bo$) được thiết lập tại công thức (2-8) và (2-9) để mô tả tỷ số giữa nhiệt lượng hữu ích truyền trong buồng lửa $Q_{bl}$ và bức xạ nhiệt cực đại của môi trường cháy, qua đó tính toán chính xác nhiệt độ khói thải rời khỏi buồng lửa $t''_{bl}$ [$^\circ\text{C}$].

2. Điều kiện tối ưu của nồng độ khí thải trong buồng đốt được xác định như thế nào?
Theo kết quả phân tích cân bằng cháy ở Hình 2-6 và Hình 2-7, chế độ cháy hoàn toàn đạt hiệu quả cao nhất khi nồng độ $\text{CO}_2$ trong khói duy trì ở mức $14,5 - 15%$ và nồng độ $\text{O}_2$ thừa nằm trong khoảng $2 - 3%$.

3. Hiện tượng dao động mức nước bao hơi ("nước sôi giả") được giải thích qua phương trình nào?
Hiện tượng này được mô tả qua phương trình vi phân (3-16) và (3-17), thể hiện tốc độ thay đổi mức nước $\frac{dh}{d\tau}$ phụ thuộc trực tiếp vào biến thiên thể tích bọt hơi ngập dưới mặt nước $\frac{dV_{ho}}{d\tau}$ khi có sự thay đổi đột ngột về áp suất $\frac{dp}{d\tau}$ hoặc tải tiêu thụ hơi $\frac{dD}{d\tau}$.

4. Vì sao đồ án đề xuất ứng dụng Aspen Hysys kết hợp Matlab/Simulink?
Việc sử dụng Aspen Hysys cùng Matlab/Simulink Toolboxes cho phép tự động hóa quá trình tra cứu trạng thái nhiệt động (enthalpy, tỷ trọng pha) và giải đồng thời các hệ phương trình vi phân phi tuyến động học lò hơi, khắc phục các giới hạn tính toán rời rạc thủ công.

5. Hiện tượng ăn mòn nhiệt độ cao nào được cảnh báo trong quá trình truyền nhiệt?
Văn bản cảnh báo hiện tượng ăn mòn tro xỉ do hợp chất vanadi $V_2O_5$ xảy ra ở dải nhiệt độ $550 - 580^\circ\text{C}$ (đặc biệt khi nhiệt độ vượt trên $550^\circ\text{C}$), cùng với nguy cơ ăn mòn do điểm sương axit $H_2SO_4$ tại đuôi lò.

Kết luận

Đồ án tốt nghiệp "Tối ưu hóa quy trình công nghiệp với công nghệ mới" đã giải quyết trọn vẹn bài toán mô hình hóa nhiệt động lực học và thiết kế hệ thống điều khiển tự động cho lò hơi và buồng đốt công nghiệp. Tác giả đã kết hợp chặt chẽ giữa lý thuyết truyền nhiệt bức xạ, cơ chế phản ứng cháy nhiên liệu và hệ phương trình vi phân mô tả động thái áp suất, mức nước bao hơi. Việc tích hợp các giải pháp điều khiển DCS, PLC (với các vòng điều khiển PIC, TIC, FIC) cùng công cụ mô phỏng Aspen Hysys và Matlab/Simulink mang lại giá trị thực tiễn cao cho công tác thiết kế, vận hành và tối ưu hóa hệ thống năng lượng trong các nhà máy sản xuất hiện đại.