Tổng quan về luận án

Trong ngành công nghiệp sản xuất giấy hiện đại, khâu sấy là công đoạn tiêu tốn nhiều năng lượng nhất (chiếm tới 60–70% tổng năng lượng nhiệt toàn dây chuyền) và giữ vai trò quyết định đối với các chỉ tiêu cơ lý tính của thành phẩm. Theo tiêu chuẩn kỹ thuật TCVN 1867, độ ẩm của giấy thành phẩm bắt buộc phải duy trì ổn định nghiêm ngặt trong khoảng 5–7% (dung sai cho phép $\pm 1%$). Nếu độ ẩm vượt quá giới hạn, liên kết hydro giữa các sớ sợi xenlulo bị suy giảm nghiêm trọng, dẫn đến hiện tượng nhòe mực khi in ấn; ngược lại, nếu giấy quá khô, độ giòn tăng cao sẽ gây đứt băng giấy trên máy xeo tốc độ cao, làm ngắt quãng dây chuyền sản xuất và gia tăng đáng kể phế phẩm.

Về bản chất nhiệt động học, buồng sấy giấy là một đối tượng công nghiệp đa biến (MIMO – Multiple-Input Multiple-Output) cực kỳ phức tạp với đặc tính phi tuyến cao, quán tính lớn và tồn tại các tác động xen kênh sâu sắc giữa hai cơ chế sấy: sấy tiếp xúc (dùng nhiệt hơi bão hòa cấp vào các quả lô sấy kim loại) và sấy đối lưu (dùng gió nóng thổi trực tiếp lên bề mặt băng giấy). Mặc dù các dây chuyền xeo hiện đại đã bổ sung hệ thống khí động học buồng sấy, thực tế vận hành công nghiệp (điển hình tại Nhà máy Giấy Bãi Bằng hay các dây chuyền lắp đặt bởi Forbes Marshall) vẫn sử dụng cấu trúc các mạch vòng điều khiển phản hồi PID đơn biến (SISO) hoạt động hoàn toàn độc lập. Sự tách biệt này bỏ qua tương tác động học giữa nhiệt độ điểm sương ($T_d$), mặt phẳng áp suất không (Zero Level - NPP) và áp suất hơi bão hòa ($P_L$), dẫn đến việc chỉnh định thông số vận hành chủ yếu dựa vào kinh nghiệm thực nghiệm, gây dao động độ ẩm lớn khi có nhiễu và làm lãng phí nghiêm trọng nguồn năng lượng hơi bão hòa.

Luận án tiến sĩ kỹ thuật "Nghiên cứu nâng cao chất lượng điều khiển phản hồi đầu ra cho quá trình đa biến buồng sấy giấy" của Nghiên cứu sinh Trần Kim Quyên, dưới sự hướng dẫn khoa học của PGS. Bùi Quốc Khánh và PGS. Đoàn Quang Vinh tại Đại học Đà Nẵng (chuyên ngành Kỹ thuật điều khiển và Tự động hóa, Mã số: 62520216), đã giải quyết triệt để rào cản kỹ thuật này. Luận án đặt ra hệ thống câu hỏi nghiên cứu và giả thuyết khoa học cốt lõi:

  • Câu hỏi nghiên cứu 1 ($RQ_1$): Làm thế nào để mô hình hóa toán học chính xác động học truyền nhiệt và truyền khối của băng giấy chuyển động liên tục qua 4 giai đoạn sấy trên từng lô sấy khi chịu tác động đồng thời của sấy tiếp xúc và sấy đối lưu?
  • Câu hỏi nghiên cứu 2 ($RQ_2$): Đặc tính xen kênh giữa mạch vòng điều khiển độ ẩm ($\gamma$), mạch vòng cân bằng gió vào – ra ($ZL$) và mạch vòng nhiệt độ điểm sương ($T_d$) ảnh hưởng như thế nào đến độ ổn định của hệ thống khi xuất hiện nhiễu độ ẩm đầu vào ($\Delta \gamma_1$) và biến động môi trường?
  • Câu hỏi nghiên cứu 3 ($RQ_3$): Cấu trúc điều khiển đa biến tách kênh (Decoupling Control) và thuật toán điều khiển dự báo mô hình (MPC) theo phương pháp tối ưu hóa từng đoạn trong miền thời gian thực có thể triệt tiêu ảnh hưởng xen kênh và nâng cao chất lượng điều khiển đầu ra như thế nào so với các bộ điều khiển GPC truyền thống?

Trên cơ sở đó, các giả thuyết nghiên cứu được xác lập:

  • Giả thuyết 1 ($H_1$): Việc tích hợp phương trình nhiệt động Clausius-Clapeyron, phương trình áp suất hơi Antoine và công thức Magnus vào mô hình cân bằng năng lượng - khối lượng phân đoạn sẽ mô tả chính xác tốc độ bay hơi nước ($q_{bh}$) trên bề mặt băng giấy.
  • Giả thuyết 2 ($H_2$): Thiết kế ma trận bù Feed-forward đa biến sẽ phân tách hoàn toàn các kênh ghép nối giữa biến điều khiển áp suất hơi ($\Delta P^$) và lưu lượng gió ($W_{a1}^, W_{a2}^*$).
  • Giả thuyết 3 ($H_3$): Ứng dụng điều khiển MPC tối ưu hóa từng đoạn trong miền thời gian sẽ tối ưu hóa quỹ đạo phản hồi, triệt tiêu ảnh hưởng của nhiễu độ ẩm đầu vào lên tới $\Delta \gamma_1 = 25%$ và biến thiên điểm đặt $\Delta T_d = \pm 10%$, giúp tiết kiệm đến 40% lượng hơi tiêu thụ và cho phép tăng tốc độ máy xeo lên 1,2 lần.

Phạm vi nghiên cứu bao quát toàn bộ động học buồng sấy trong dây chuyền xeo công nghiệp với định lượng giấy từ 52 đến 80 $g/m^2$, kết hợp số liệu thực nghiệm thu thập liên tục trong 4 năm tại Tổng công ty Giấy Bãi Bằng.

Literature Review và Positioning

Lịch sử nghiên cứu về công nghệ sấy và điều khiển quá trình sấy giấy trên thế giới được định hình qua ba dòng nghiên cứu chính:

Dòng nghiên cứu thứ nhất tập trung vào động học truyền nhiệt sấy tiếp xúc truyền thống qua vỏ lô sấy. Các công trình kinh điển của Karlsson (2000), Wilhelmsson (1995) và Stenström (1999) đã xây dựng các mô hình truyền nhiệt một chiều từ hơi bão hòa ngưng tụ bên trong lòng lô qua lớp vỏ gang, lớp nước ngưng và chăn sấy truyền vào băng giấy. Slätteke (2006) tại Đại học Lund (Thụy Điển) đã phát triển mô hình hàm truyền đạt dạng IPZ (Integrator-Pole-Zero) để mô tả quan hệ giữa áp suất hơi cấp và độ ẩm đầu ra. Tuy nhiên, các mô hình này hoàn toàn xem môi trường không khí xung quanh là hằng số, bỏ qua tác động động học của hệ thống thông gió buồng sấy.

Dòng nghiên cứu thứ hai khai thác cơ chế sấy đối lưu và khí động học buồng sấy kín. Poirier & Mujumdar (1991) cùng các nghiên cứu kỹ thuật của tập đoàn Valmet và Forbes Marshall (2013) chứng minh rằng việc thổi gió nóng trực tiếp lên bề mặt giấy giúp phá vỡ lớp biên ẩm, gia tăng hệ số chuyển khối $K$. Mặc dù vậy, các tài liệu kỹ thuật này chỉ dừng lại ở các giải pháp thiết kế cơ khí và triển khai các mạch vòng điều khiển PID đơn biến độc lập cho quạt cấp gió và quạt hút khí thải, dẫn đến hiện tượng dao động áp suất buồng sấy và ngưng tụ hơi nước ngược trở lại băng giấy.

Dòng nghiên cứu thứ ba tập trung vào lý thuyết điều khiển đa biến và điều khiển dự báo mô hình (MPC) ứng dụng trong công nghiệp giấy. Clarke et al. (1987) đề xuất thuật toán Generalized Predictive Control (GPC), tiếp sau đó là Dynamic Matrix Control (DMC) của Cutler & Ramaker (1980) và Model Algorithmic Control (MAC). Camacho & Bordons (2004) đã hệ thống hóa MPC cho các hệ thống có ràng buộc. Tuy nhiên, việc áp dụng MPC vào buồng sấy giấy thường gặp thách thức lớn do tính phi tuyến mạnh của điểm sương và khối lượng tính toán trực tuyến quá lớn trên nền tảng phần cứng DCS công nghiệp.

Trong bức tranh tổng quan đó, tồn tại một cuộc tranh luận học thuật sâu sắc giữa hai quan điểm đối nghịch: Một bên ủng hộ việc duy trì các bộ điều khiển PID phân tán phi tập trung nhờ tính đơn giản, dễ hiệu chỉnh tại hiện trường; bên còn lại chỉ ra rằng việc bỏ qua ghép kênh trong hệ sấy khí đối lưu khiến hệ thống không thể triệt tiêu hiện tượng "mưa trong buồng sấy" hoặc quá khô gây cháy giấy khi vận hành ở dải tốc độ cao.

Luận án của NCS. Trần Kim Quyên đã xác lập vị trí tiên phong bằng cách bắc nhịp cầu hoàn hảo giữa hai trường phái: Xây dựng một mô hình toán học giải tích nhiệt động học 4 pha toàn diện, đồng thời thiết kế cấu trúc điều khiển đa biến tách kênh kết hợp MPC tối ưu hóa từng đoạn trong miền thời gian thực. So với công trình của Slätteke (2006, Thụy Điển) vốn chỉ điều khiển áp suất hơi, và hệ thống của Forbes Marshall (2013, Anh Quốc) chỉ dùng PID độc lập cho quạt gió, nghiên cứu này giải quyết trọn vẹn bài toán điều khiển đồng thời cả 4 thông số trạng thái: độ ẩm giấy, nhiệt độ điểm sương, điểm áp suất không và nhiệt độ gió cấp.

Đó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 nhiệt động lực học quá trình sấy và lý thuyết điều khiển hệ thống đa biến phi tuyến thông qua các đóng góp cụ thể:

Thứ nhất, luận án đã mở rộng lý thuyết truyền chất của Treybal và định luật nhiệt động Clausius-Clapeyron trong môi trường sợi xốp chuyển động liên tục. Bằng cách kết hợp phương trình thực nghiệm áp suất hơi bão hòa Antoine:

$$\log_{10}(p_v^0) = A - \frac{B}{T_p + C}$$

với phương trình cân bằng Clausius-Clapeyron:

$$\Delta H_s = -\frac{R_g}{M_w} \frac{d(\ln p_v)}{d(1/T_p)}$$

tác giả đã thiết lập được công thức giải tích xác định chính xác năng lượng cần thiết để làm bay hơi nước từ bề mặt băng giấy ($\Delta H_{bh} = 2260\text{ kJ/kg}$) và nhiệt dung riêng quy đổi của băng giấy ẩm:

$$C_p = \frac{C_{p,f} + \gamma C_{p,w}}{1 + \gamma}$$

trong đó nhiệt dung của sợi giấy khô được xác định chính xác là $C_{p,f} = 1256\text{ J/(kg}\cdot\text{K)}$.

Thứ hai, luận án xây dựng hệ phương trình vi phân mô tả quy luật biến thiên vi mô của độ ẩm và nhiệt độ băng giấy qua 4 giai đoạn sấy kế tiếp nhau trên một quả lô sấy:

  1. Giai đoạn 2 (Gia nhiệt kín): Băng giấy được ép sát mặt lô nhờ chăn sấy ($w_{bh} \approx 0$), phương trình truyền nhiệt thuần túy: $$m_p (\gamma + 1) C_p \frac{dT_p}{dt} = w_p (1 + \gamma_1) C_p (T_{p1} - T_{p2}) + K_{T1} A_1 \eta (T_m - T_p)$$
  2. Giai đoạn 1 & 3 (Bay hơi tự do mặt thoáng): Giấy rời chăn sấy nhưng còn tiếp xúc mặt lô: $$m_p (\gamma + 1) C_p \frac{dT_p}{dt} = w_p (1 + \gamma_1) C_p (T_{p1} - T_{p2}) + K_{T1} A_1 \eta (T_m - T_p) - w_{bh1} \Delta H_{b1}$$
  3. Giai đoạn 4 (Sấy đối lưu cưỡng bức): Giấy nằm giữa hai khoảng lô, nhận nhiệt hoàn toàn từ vòi phun gió nóng: $$m_p (\gamma + 1) C_p \frac{dT_{p2}}{dt} = K_{T2} A_2 \eta (T_G - T_p) + w_p (1 + \gamma_1) C_p (T_{p1} - T_{p2}) + w_{bh} C_p T_{p2} - w_{bh} \Delta H_b$$

Thứ ba, luận án đã chứng minh toán học định lý về tính ổn định của mặt phẳng áp suất không (Neutral Pressure Plane - NPP) trong buồng sấy kín, thông qua phương trình biến thiên áp suất tại khe đo:

$$\frac{dp}{dt} = \frac{R T}{M V} \left( W_{a1} + W_{bh} + a C_v \sqrt{\frac{2 |\Delta P|}{\nu}} - W_{a2} \right)$$

Chứng minh này đặt nền tảng lý thuyết vững chắc cho việc kiểm soát hiện tượng xâm nhập khí lạnh từ môi trường ngoài vào buồng sấy.

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

Khung phân tích của luận án tích hợp liên ngành ba trụ cột lý thuyết: Lý thuyết Nhiệt - Ẩm công nghiệp, Lý thuyết Điều khiển Đa biến Hiện đại và Lý thuyết Tối ưu hóa Ràng buộc Không gian Trạng thái.

Mô hình cấu trúc điều khiển buồng sấy được mô tả tổng quát dưới dạng ma trận hàm truyền đạt $3 \times 3$ hoặc mô hình không gian trạng thái liên tục:

$$\begin{bmatrix} \gamma(s) \ ZL(s) \ T_d(s) \end{bmatrix} = \begin{bmatrix} G_{11}(s) & G_{12}(s) & G_{13}(s) \ G_{21}(s) & G_{22}(s) & G_{23}(s) \ G_{31}(s) & G_{32}(s) & G_{33}(s) \end{bmatrix} \begin{bmatrix} \Delta P^(s) \ W_{a1}^(s) \ W_{a2}^*(s) \end{bmatrix} + \begin{bmatrix} G_{D1}(s) \ G_{D2}(s) \ G_{D3}(s) \end{bmatrix} \Delta \gamma_1(s)$$

Các điều kiện biên được xác lập tường minh: Nhiệt độ không khí trong buồng $0^\circ\text{C} < T < 160^\circ\text{C}$, độ ẩm tương đối $0 < RH < 100%$, và nhiệt độ điểm sương $T_d$ luôn phải duy trì thấp hơn nhiệt độ trung bình buồng sấy một khoảng an toàn từ $15^\circ\text{C}$ đến $20^\circ\text{C}$ để triệt tiêu hoàn toàn khả năng ngưng tụ giọt nước trong buồng sấy.

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ thế giới quan thực chứng (Positivism) kết hợp phương pháp mô hình hóa hộp xám (Grey-box Modeling), dung hòa hoàn hảo giữa mô hình toán học giải tích xây dựng từ định luật bảo toàn vật lý và việc nhận dạng tham số từ thực nghiệm sản xuất.

Thiết kế hệ thống được cấu trúc theo 3 cấp độ phân tầng:

  1. Cấp 1 (Cơ cấu chấp hành & Cảm biến): Bao gồm van điều áp tuyến tính $C_v$, quạt thổi biến tần $W_{a1}$, quạt hút $W_{a2}$, hệ thống Scanner đo độ ẩm online theo chiều dọc (MD) và ngang (CD) băng giấy, cùng hệ thống cảm biến nhiệt độ phân tầng đo điểm áp suất không.
  2. Cấp 2 (Điều khiển giải ghép & Tối ưu hóa): Bộ điều khiển đa biến Feed-forward kết hợp thuật toán MPC tuyến tính tối ưu hóa từng đoạn.
  3. Cấp 3 (Giám sát QCS): Hệ thống điều khiển chất lượng tổng thể dây chuyền xeo.

Quy trình nghiên cứu rigorous

Để thu thập số liệu phục vụ nhận dạng và kiểm chứng, tác giả đã thiết kế giao thức đo đạc thực nghiệm đặc biệt cho cảm biến NPP (Neutral Pressure Plane Sensor Array). Dàn 4 cảm biến nhiệt độ được bố trí chính xác tại khe gió cạnh buồng sấy: cảm biến $n$ đặt tại mặt trên của lô dưới (nhiệt độ $T_n$), cảm biến $n+1$ phía trên ($T_{n+1}$), cảm biến $n-1$ ($T_{n-1}$) và $n-2$ ($T_{n-2}$) phía dưới.

       Cửa thoát khí thải (Quạt hút Wa2)
       Cửa cấp gió nóng (Quạt thổi Wa1)

Quy tắc xác định vị trí mặt phẳng áp suất không dựa trên đạo hàm nhiệt độ gián đoạn:

  • Nếu $\Delta T_1 = T_n - T_{n-1}$ đạt cực đại $\rightarrow$ Điểm NP nằm chuẩn xác tại vị trí $n$.
  • Nếu $\Delta T_2 = T_{n+1} - T_n$ đạt cực đại $\rightarrow$ Điểm NP bị dịch lên trên (gió lạnh $W_{kk}$ tràn vào).
  • Nếu $\Delta T_3 = T_{n-1} - T_{n-2}$ đạt cực đại $\rightarrow$ Điểm NP bị tụt xuống dưới (thất thoát nhiệt ra ngoài).

Mọi thiết bị đo lường chuẩn hóa theo tiêu chuẩn quốc gia TCVN 1270 (định lượng), TCVN 3652 (độ chặt), TCVN 1867 (độ ẩm) và TCVN 1865 (độ trắng), đảm bảo độ tin cậy và tính lặp lại tuyệt đối của dữ liệu thực nghiệm.

Data và phân tích

Quá trình mô phỏng động học và kiểm chứng thuật toán điều khiển được thực hiện trên môi trường MATLAB/Simulink với bộ thông số công nghệ chính xác từ nhà máy:

Thông số công nghệ Giá trị thực nghiệm Đơn vị đo Ý nghĩa vật lý
Định lượng giấy chuẩn ($g$) $52 - 80$ $g/m^2$ Khối lượng giấy trên đơn vị diện tích
Độ ẩm giấy đầu vào khâu sấy ($\gamma_1$) $50$ $%$ Độ ẩm sau tổ hợp ép ướt
Độ ẩm giấy thành phẩm yêu cầu ($\gamma$) $5 - 7 \pm 1$ $%$ Tiêu chuẩn chất lượng TCVN 1867
Nhiệt dung riêng sợi giấy ($C_{p,f}$) $1256$ $J/(kg\cdot K)$ Nhiệt dung phần xơ sợi khô
Nhiệt ẩn hóa hơi nước ($\Delta H_{bh}$) $2260$ $kJ/kg$ Năng lượng bốc hơi ở áp suất khí quyển
Nhiệt độ buồng sấy trung bình ($T$) $80 - 110$ $^\circ\text{C}$ Nhiệt độ không khí làm việc
Nhiệt độ điểm sương đặt ($T_d^*$) $T - (15 \div 20)$ $^\circ\text{C}$ Ngưỡng an toàn chống đọng sương
Bước nhảy nhiễu ẩm đầu vào ($\Delta \gamma_1$) $+25$ $%$ Xung nhiễu thử nghiệm độ bền vững
Bước nhảy điểm đặt điểm sương ($\Delta T_d$) $\pm 10$ $%$ Thử nghiệm khả năng thích ứng

Các thử nghiệm Robustness checks được thiết lập khắt khe thông qua việc tác động đồng thời xung nhiễu độ ẩm đầu vào $\Delta \gamma_1 = 25%$ (mô phỏng sự cố tại máy ép keo) và thay đổi đột ngột điểm đặt nhiệt độ điểm sương $\Delta T_d = \pm 10%$.

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

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

Thứ nhất, nghiên cứu phát hiện và lượng hóa chính xác mức độ ghép kênh chéo giữa mạch vòng độ ẩm và mạch vòng điểm sương. Khi điểm đặt nhiệt độ điểm sương thay đổi $\Delta T_d = \pm 10%$, nếu sử dụng hệ thống PID đơn biến truyền thống, độ ẩm băng giấy $\gamma$ bị dao động lệch khỏi giá trị đặt tới $\pm 1.8%$, kéo dài thời gian quá độ lên hơn 120 giây. Khi áp dụng bộ điều khiển đa biến tách kênh, độ lệch này được triệt tiêu hoàn toàn về mức $0%$, đưa thời gian xác lập về dưới 15 giây.

Thứ hai, kết quả thử nghiệm cho thấy đáp ứng vượt trội của bộ điều khiển MPC tối ưu hóa từng đoạn trong miền thời gian thực khi chịu xung nhiễu độ ẩm đầu vào cực đại $\Delta \gamma_1 = 25%$. Trong khi bộ điều khiển PID truyền thống làm giấy bị ẩm vượt ngưỡng $9.2%$ (nguy cơ đứt giấy) và bộ điều khiển GPC trong Matlab Toolbox vẫn để độ ẩm lệch $\pm 0.8%$, thì thuật toán MPC tối ưu hóa từng đoạn của luận án khống chế độ ẩm dao động không quá $\pm 0.15%$, đưa hệ thống trở lại trạng thái cân bằng trong vòng 8 giây.

  Độ ẩm giấy (%)
            0       20       40       60       80      100      120
           [Tác động nhiễu ẩm Δγ1 = +25%]

Thứ ba, kiểm soát chính xác mặt phẳng áp suất không (Zero Level) tại vị trí lô dưới giúp giảm lưu lượng gió lạnh lọt vào buồng sấy ($W_{kk}$) từ $1.85\text{ kg/s}$ xuống xấp xỉ $0.05\text{ kg/s}$, duy trì nhiệt độ buồng sấy ổn định mà không cần tăng lưu lượng sấy hơi.

Thứ tư, cơ chế sấy kết hợp giúp giảm căn bản suất tiêu hao hơi bão hòa cấp cho các cụm lô sấy. Số liệu định lượng chứng minh việc điều khiển tối ưu môi trường nhiệt ẩm buồng sấy giúp tiết kiệm tới 40% lượng hơi bão hòa cấp vào lô sấy và cho phép dây chuyền tăng tốc độ xeo lên 1,2 lần mà vẫn đảm bảo độ ẩm giấy đầu ra đạt chuẩn $5 - 7%$.

Implications đa chiều

Về mặt lý thuyết, nghiên cứu cung cấp một phương pháp luận hoàn chỉnh để mô hình hóa các hệ thống nhiệt động đa biến phi tuyến có chuyển pha chất lỏng - khí, đóng góp trực tiếp vào kho tàng lý thuyết điều khiển quá trình hóa học và công nghệ vật liệu xốp.

Về mặt phương pháp luận, thuật toán MPC tối ưu hóa từng đoạn trong miền thời gian thực chứng minh tính khả thi trong việc giải bài toán tối ưu hóa có ràng buộc phức tạp mà không đòi hỏi năng lực tính toán quá cao, mở ra hướng ứng dụng rộng rãi cho các hệ thống sấy công nghiệp khác như sấy nông sản, sấy gỗ, và sấy màng phim polymer.

Về mặt thực tiễn và chính sách, giải pháp của luận án mang lại lợi ích kép: Giúp các doanh nghiệp sản xuất giấy nâng cao năng suất sản xuất thêm 20%, giảm chi phí nhiên liệu đốt nồi hơi (than đá, sinh khối hoặc dầu FO), giảm phát thải khí nhà kính $CO_2$, đóng góp trực tiếp vào cam kết quốc gia về sử dụng năng lượng tiết kiệm và hiệu quả theo Nghị quyết phát triển kinh tế xanh của Chính phủ.

Limitations và Future Research

Dù đạt được những kết quả đột phá, luận án thẳng thắn chỉ ra 3 hạn chế nội tại:

  1. Mô hình hóa không gian 1 chiều (1D Model): Nghiên cứu chỉ tập trung phân tích động học độ ẩm dọc theo chiều chuyển động của băng giấy (Machine Direction - MD), giả thiết độ ẩm theo chiều ngang khổ giấy (Cross Direction - CD) là phân bố đồng đều. Trên thực tế các máy xeo khổ rộng ($> 4\text{ m}$), sự không đồng đều theo chiều ngang vẫn đòi hỏi các cơ cấu vòi phun vi chỉnh CD độc lập.
  2. Giả thiết lý tưởng hóa bề mặt truyền nhiệt: Mô hình chưa tính đến sự suy giảm hệ số truyền nhiệt $K_{T1}$ do hiện tượng đóng cặn cáu bẩn bên trong thành lô và sự lão hóa cơ học của chăn bạt sấy sau thời gian dài vận hành liên tục.
  3. Thực thi thuật toán: Thuật toán MPC tối ưu hóa từng đoạn mới được kiểm chứng mô phỏng thời gian thực (Real-time Simulation) và kiểm thử trên các đoạn mạch vòng thử nghiệm bán công nghiệp, chưa được nạp trực tiếp toàn phần vào firmware của hệ thống DCS phân tán thế hệ cũ.

Từ những hạn chế trên, 4 hướng nghiên cứu tiếp nối trong tương lai được vạch ra:

  • Mở rộng mô hình động học hệ thống sang hệ phương trình đạo hàm riêng 2 chiều (2D PDE) kết hợp đồng thời điều khiển phân bố ẩm MD và CD.
  • Ứng dụng mạng nơ-ron sâu tự thích nghi (Deep Reinforcement Learning MPC) để ước lượng trực tuyến sự biến thiên của hệ số truyền nhiệt $K_{T1}, K_{T2}$ theo thời gian thực.
  • Phát triển phần cứng nhúng biên (Edge-AI Embedded System) tích hợp trực tiếp thuật toán MPC giải trình tự tối ưu vào các bộ điều khiển khả trình PLC/DCS chuẩn công nghiệp.
  • Nghiên cứu tích hợp hệ thống thu hồi nhiệt thải nâng cao (Advanced Heat Recovery Unit - HRU) đa tầng nhằm tối ưu hóa nhiệt động học toàn bộ chu trình nhiệt nhà máy giấy.

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

Luận án đã tạo ra những tác động khoa học và công nghệ rõ nét:

Ảnh hưởng học thuật: Các bài báo khoa học trích xuất từ luận án được công bố trên các tạp chí chuyên ngành uy tín trong nước và quốc tế, trở thành tài liệu tham khảo nền tảng cho các nghiên cứu sinh chuyên ngành Tự động hóa và Công nghệ Hóa học.

Chuyển giao công nghệ công nghiệp: Giải pháp điều khiển tách kênh và cân bằng Zero Level cung cấp cẩm nang kỹ thuật thực tiễn cho các kỹ sư trưởng, kỹ sư điều khiển quá trình tại các nhà máy giấy lớn (Bãi Bằng, Tân Mai, An Hòa) trong việc cấu hình lại hệ thống điều khiển buồng sấy mà không cần đầu tư thay thế toàn bộ thiết bị cơ khí đắt tiền.

Tác động xã hội và môi trường: Việc giảm 40% tiêu thụ hơi bão hòa tương đương với việc tiết kiệm hàng nghìn tấn than tiêu chuẩn mỗi năm cho một nhà máy công suất trung bình, góp phần giảm thiểu bụi mịn và khí thải công nghiệp, nâng cao sức cạnh tranh của sản phẩm giấy nội địa trên thị trường quốc tế.

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

  • Nghiên cứu sinh & Nhà nghiên cứu trẻ: Tiếp cận một mô hình toán học mẫu mực kết hợp giữa phương trình vật lý nhiệt động học sâu sắc với kỹ thuật điều khiển hiện đại, giải quyết trọn vẹn bài toán hệ MIMO ghép kênh mạnh.
  • Giảng viên & Nhà khoa học cao cấp: Có thêm tài liệu giảng dạy chuyên sâu, hệ thống hóa lý thuyết điều khiển dự báo MPC và điều khiển tách kênh ứng dụng trong các đối tượng nhiệt công nghiệp phức tạp.
  • Kỹ sư R&D & Doanh nghiệp tự động hóa: Sở hữu thuật toán tối ưu hóa từng đoạn có thể lập trình trực tiếp lên PLC/SCADA để nâng cấp các dây chuyền sản xuất giấy hiện hữu.
  • Nhà quản lý & Hoạch định chính sách năng lượng: Có cơ sở khoa học xác đáng để xây dựng định mức tiêu hao năng lượng và tiêu chuẩn kỹ thuật quốc gia cho ngành công nghiệp sản xuất giấy và bột giấy.

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

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

Đóng góp độc đáo nhất là việc xây dựng thành công Mô hình động học 4 pha nhiệt - ẩm kết hợp trên từng lô sấy, tích hợp trực tiếp định luật Clausius-Clapeyron và phương trình Antoine vào phương trình vi phân cân bằng khối lượng và năng lượng của băng giấy. Luận án đã mở rộng lý thuyết sấy tiếp xúc cổ điển của Karlsson và Slätteke bằng cách bổ sung động học truyền khối đối lưu cưỡng bức chịu sự chi phối của nhiệt độ điểm sương buồng sấy ($T_d$).

2. Điểm cải tiến phương pháp luận so với các nghiên cứu quốc tế tiền nhiệm?

So với nghiên cứu của Slätteke (2006, Thụy Điển - chỉ xét hệ đơn biến áp suất hơi) và hệ thống của Forbes Marshall (2013 - dùng các mạch vòng PID độc lập), luận án đã thiết lập Ma trận tách kênh Feed-forward giải ghép hoàn toàn kết hợp MPC tối ưu hóa từng đoạn trong miền thời gian thực, cho phép xử lý đồng thời 3 biến vào, 3 biến ra có ràng buộc và triệt tiêu tương tác chéo trong buồng sấy kín.

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

Phát hiện bất ngờ nhất là Nhiệt độ điểm sương ($T_d$) có tác động xuyên kênh tức thời và mạnh hơn cả nhiệt độ gió cấp đối với tốc độ bay hơi nước ($q_{bh}$). Khi $T_d$ tăng tiến sát nhiệt độ buồng sấy (chênh lệch $< 15^\circ\text{C}$), áp suất riêng phần của hơi nước trong không khí tăng vọt làm triệt tiêu động lực chuyển khối, khiến băng giấy hoàn toàn không thể thoát ẩm dù công suất cấp nhiệt hơi ở mức tối đa.

4. Giao thức tái lập thực nghiệm (Replication Protocol) có được cung cấp đầy đủ không?

Luận án cung cấp tường minh 100% hệ phương trình vi phân toán học (từ phương trình 1.4 đến 2.15), bảng thông số vật lý của giấy ($C_p, g, \Delta H_{bh}$), sơ đồ kết nối khối MATLAB/Simulink và cấu hình dàn 4 cảm biến nhiệt độ phân tầng xác định mặt phẳng áp suất không NPP, cho phép bất kỳ phòng thí nghiệm nào cũng có thể tái lập hoàn toàn kết quả nghiên cứu.

5. Lộ trình nghiên cứu 10 năm tiếp theo được định hình ra sao?

Lộ trình 10 năm bao gồm: Phát triển mô hình điều khiển phân tán 2 chiều đồng thời MD-CD; Chuyển đổi thuật toán MPC sang kiến trúc mạng nơ-ron thích nghi sâu nhúng trên chip chuyên dụng (Edge-AI); Ứng dụng công nghệ bản sao số (Digital Twin) để chẩn đoán hư hỏng bạt sấy và cáu cặn lô sấy thời gian thực.

Kết luận

Luận án tiến sĩ của tác giả Trần Kim Quyên là một công trình nghiên cứu khoa học công phu, nghiêm túc và mẫu mực, giải quyết trọn vẹn một bài toán công nghệ bức thiết của ngành công nghiệp giấy với 5 đóng góp cốt lõi:

  1. Thiết lập hoàn chỉnh mô hình toán học động học 4 pha truyền nhiệt và truyền ẩm trên băng giấy chuyển động liên tục, tích hợp chính xác các quy luật nhiệt động học Clausius-Clapeyron, Antoine và Magnus.
  2. Xác lập mô hình động học và thuật toán định vị mặt phẳng áp suất không (Neutral Pressure Plane - NPP) tại khe buồng sấy, triệt tiêu sự xâm nhập của không khí lạnh ngoại vi.
  3. Thiết kế thành công bộ điều khiển đa biến tách kênh Feed-forward, phân rã hoàn toàn tính xen kênh giữa mạch vòng độ ẩm, nhiệt độ điểm sương và áp suất buồng sấy.
  4. Ứng dụng xuất sắc thuật toán điều khiển dự báo MPC tối ưu hóa từng đoạn trong miền thời gian thực, khống chế độ lệch ẩm dưới $\pm 0.15%$ trước các xung nhiễu khắc nghiệt lên tới $25%$.
  5. Mang lại hiệu quả kinh tế - kỹ thuật mang tính đột phá với việc tiết kiệm đến 40% năng lượng hơi bão hòa và tăng 20% tốc độ vận hành dây chuyền xeo giấy.

Công trình không chỉ xác lập một chuẩn mực mới trong điều khiển tự động hóa quá trình sấy công nghiệp mà còn mở ra những hướng nghiên cứu chuyên sâu về tối ưu hóa năng lượng cho nền sản xuất thông minh và bền vững.