Luận Án Tiến Sĩ: Ứng Dụng Điều Khiển Dự Báo Phi Tuyến Cho Thiết Bị Phản Ứng Khuấy Trộn Liên Tục

Luận án tiến sĩ nghiên cứu kỹ thuật ứng dụng điều khiển dự báo phi tuyến cho thiết bị phản ứng khuấy trộn liên tục, phát triển phương pháp mới, đánh giá hiệu quả ứng dụng trong

• Tác giả: Mai Thị Đoan Thanh
• Người hướng dẫn: Bùi Quốc Khánh
• Trường học: Đại học Đà Nẵng
• Chuyên ngành: Kỹ thuật Điều khiển và Tự động hóa
• Đề tài: Ứng Dụng Điều Khiển Dự Báo Phi Tuyến Cho Thiết Bị Phản Ứng Khuấy Trộn Liên Tục
• Loại tài liệu: luận án tiến sỹ kỹ thuật
• Năm xuất bản: 2018
• Địa điểm: Đà Nẵng

Mục lục chi tiết

• LỜI CAM ĐOAN
• 1. CHƯƠNG 1: TỔNG QUAN VỀ THIẾT BỊ PHẢN ỨNG CÓ KHUẤY TRỘN LIÊN TỤC CSTR
• 1.1. Lý thuyết cơ bản về phản ứng hóa học
• 1.1.1. Định nghĩa thiết bị phản ứng
• 1.1.2. Lý thuyết cơ bản về phản ứng hóa học
• 1.2. Phân loại phản ứng hóa học và thiết bị phản ứng
• 1.2.1. Phân loại phản ứng hóa học
• 1.2.2. Phân loại thiết bị phản ứng
• 1.3. Động học chung thiết bị phản ứng có thể tích không đổi
• 1.4. Động học thiết bị phản ứng khuấy trộn liên tục (CSTR – Continuous Stirred Tank Reactor)
• 1.5. Xây dựng mô hình điều khiển thiết bị phản ứng CSTR
• 1.5.1. Phân tích các biến
• 1.5.2. Xác định bậc tự do của mô hình
• 1.5.3. Đánh giá tính xen kênh và tính phi tuyến
• 1.6. Tóm tắt những công trình công bố gần đây về điều khiển thiết bị phản ứng khuấy trộn liên tục CSTR
• 1.6.1. Những công trình nghiên cứu về động học quá trình thiết bị phản ứng khuấy trộn liên tục CSTR
• 1.6.2. Những công trình nghiên cứu về điều khiển phản hồi tuyến tính dùng bộ PID và các biến thể
• 1.6.3. Những công trình nghiên cứu về điều khiển phi tuyến
• 1.6.4. Những công trình nghiên cứu có phương trình cân bằng năng lượng của nước trong jacket
• 1.6.5. Các công trình nghiên cứu dùng MPC để điều khiển thiết bị phản ứng CSTR
• 1.7. Định hướng nghiên cứu CSTR
• 1.8. Kết luận chương 1
• 2. CHƯƠNG 2: THIẾT LẬP MÔ HÌNH ĐIỀU KHIỂN THIẾT BỊ PHẢN ỨNG KHUẤY TRỘN LIÊN TỤC THỦY PHÂN ANHYDRIDE ACETIC
• 2.1. Phản ứng thủy phân anhydride acetic trong công nghệ sản xuất acid acetic
• 2.2. Thiết kế mô hình thử nghiệm thiết bị phản ứng thủy phân anhydride acetic
• 2.2.1. Lựa chọn kích thước hình học
• 2.2.2. Tính toán thông số vận hành
• 2.3. Xây dựng mô hình mô phỏng kiểm tra thiết kế thiết bị phản ứng
• 2.3.1. Mô hình mô phỏng thiết bị phản ứng trong giai đoạn khởi động
• 2.3.2. Mô hình mô phỏng thiết bị phản ứng trong giai đoạn phản ứng
• 2.4. Mô phỏng điều khiển tuyến tính dùng phản hồi đầu ra với bộ điều khiển PID cho thiết bị phản ứng thủy phân anhydride acetic
• 2.4.1. Xét trường hợp khi có nhiễu tác động
• 2.4.1.1. Nhiễu lưu lượng đầu ra F2 10%
• 2.4.1.2. Nhiễu đồng thời
• 2.5. Kết luận chương 2
• 3. CHƯƠNG 3: THIẾT KẾ BỘ ĐIỀU KHIỂN DỰ BÁO THÍCH NGHI CHO HỆ SONG TUYẾN BẤT ĐỊNH VÀ ÁP DỤNG VÀO ĐIỀU KHIỂN HỆ PHẢN ỨNG KHUẤY TRỘN LIÊN TỤC THỦY PHÂN ANHYDRIDE ACETIC
• 3.1. Tổng quan về điều khiển dự báo
• 3.1.1. Cấu trúc và nguyên lý làm việc của bộ điều khiển dự báo
• 3.1.2. Các phương pháp điều khiển dự báo cơ bản
• 3.2. Xây dựng bộ điều khiển dự báo thích nghi cho hệ song tuyến không dừng và bất định
• 3.2.1. Xây dựng thuật toán điều khiển dự báo cho hệ song tuyến. Ứng dụng vào điều khiển thiết bị phản ứng khuấy trộn liên tục thủy phân anhydride acetic
• 3.3. Kết luận chương 3
• 4. CHƯƠNG 4: MỤC TIÊU CỦA THỰC NGHIỆM
• 4.1. Mô hình thí nghiệm. Trình tự thực nghiệm
• 4.2. Lập trình bộ điều khiển AC 800M
• 4.2.1. Cấu hình phần cứng
• 4.2.2. Thiết kế và xây dựng giao diện
• 4.3. Thí nghiệm xác định đặc tính của hệ
• 4.3.1. Thí nghiệm lấy đặc tính bơm B1, B2
• 4.3.2. Nhận dạng vòng hở hệ thống
• 4.4. Mô phỏng hệ điều khiển với thông số mô hình thí nghiệm
• 4.4.1. Sự phụ thuộc của nồng độ chất phản ứng vào nhiệt độ lò phản ứng T
• 4.4.2. Mô phỏng hệ điều khiển
• 4.5. Thí nghiệm với bộ điều khiển PID
• 4.5.1. Chỉnh định thông số bộ điều khiển PID bằng thực nghiệm
• 4.5.1.1. Tiến hành thí nghiệm
• 4.5.1.2. Sơ đồ ghép nối bộ điều khiển AC800M với mô hình thí nghiệm
• 4.5.2. Kết quả thí nghiệm với điều khiển PID
• 4.6. Thí nghiệm với bộ điều khiển NMPC
• 4.6.1. Sơ đồ điều khiển trong AC800M. Kết quả thí nghiệm điều khiển NMPC
• 4.7. Kết luận và kiến nghị
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Tính Cấp Thiết Của Luận Án

Công nghiệp hóa chất đóng vai trò quan trọng trong sản xuất công nghiệp, với thiết bị phản ứng hóa học là cốt lõi của quá trình này. Thiết bị phản ứng có khuấy trộn liên tục (CSTR) được sử dụng phổ biến nhờ năng suất cao. Tuy nhiên, đặc tính động học của hệ CSTR là phi tuyến và đa biến, dẫn đến việc sử dụng điều khiển phản hồi đầu ra tuyến tính với bộ điều khiển PID không đạt được chất lượng sản phẩm mong muốn. Việc sử dụng chuỗi thiết bị CSTR nối tiếp làm tăng chi phí đầu tư và phức tạp hóa hệ thống sản xuất. Nghiên cứu điều khiển phi tuyến cho CSTR đã mở ra hướng đi mới trong sản xuất, tuy nhiên, các nghiên cứu trước đây chủ yếu tập trung vào điều khiển đơn biến mà chưa xem xét đầy đủ cân bằng khối lượng và năng lượng cho jacket. Luận án này nhằm ứng dụng điều khiển dự báo phi tuyến cho thiết bị phản ứng CSTR, với mục tiêu khử nhiễu và đảm bảo chất lượng đầu ra.

II. Mục Tiêu Nghiên Cứu

Mục tiêu chính của luận án là thiết kế hệ điều khiển dự báo phi tuyến NMPC cho thiết bị phản ứng CSTR nhằm nâng cao chất lượng sản phẩm và dễ dàng ứng dụng vào sản xuất. Đối tượng nghiên cứu là thiết bị phản ứng CSTR thủy phân anhydride acetic. Phạm vi nghiên cứu bao gồm thiết kế điều khiển cho thiết bị phản ứng CSTR, nghiên cứu động học quá trình và điều khiển tuyến tính cho CSTR-acetic. Phương pháp nghiên cứu bao gồm khảo sát các công trình công bố về điều khiển phi tuyến CSTR, thiết kế mô hình CSTR trong phòng thí nghiệm, và đánh giá hệ điều khiển bằng mô phỏng. Kết quả nghiên cứu sẽ mở ra khả năng ứng dụng điều khiển phi tuyến NMPC cho thiết bị CSTR trong thực tế sản xuất.

III. Tổng Quan Về Thiết Bị Phản Ứng CSTR

Thiết bị phản ứng có khuấy trộn liên tục (CSTR) là một trong những thiết bị quan trọng trong ngành công nghiệp hóa chất. CSTR cho phép quá trình phản ứng diễn ra liên tục với nguyên liệu đầu vào và sản phẩm đầu ra liên tục. Đặc điểm của CSTR là khả năng duy trì điều kiện phản ứng ổn định, nhờ vào việc khuấy trộn liên tục và kiểm soát nhiệt độ thông qua jacket. Động học của quá trình phản ứng trong CSTR thường phức tạp do sự tương tác giữa các biến. Việc phân loại phản ứng hóa học và thiết bị phản ứng là cần thiết để hiểu rõ hơn về cách thức hoạt động của CSTR. Các phản ứng trong CSTR thường là phản ứng thu nhiệt và đồng pha, điều này có ảnh hưởng lớn đến thiết kế và điều khiển thiết bị. Việc nghiên cứu động học và điều khiển cho CSTR là rất quan trọng để tối ưu hóa quy trình sản xuất.

IV. Xây Dựng Mô Hình Điều Khiển

Mô hình điều khiển cho thiết bị phản ứng CSTR được xây dựng dựa trên các biến động học và điều kiện ràng buộc. Việc xác định bậc tự do của mô hình là bước quan trọng để đảm bảo tính khả thi của hệ thống điều khiển. Đánh giá tính xen kênh và tính phi tuyến của hệ thống cũng cần được thực hiện để đảm bảo rằng mô hình phản ánh chính xác các đặc tính của thiết bị. Các thuật toán điều khiển dự báo phi tuyến được phát triển nhằm tối ưu hóa quy trình và giảm thiểu nhiễu tác động. Việc áp dụng mô hình điều khiển này vào thực tế sản xuất sẽ giúp nâng cao hiệu quả và chất lượng sản phẩm đầu ra.

V. Kết Luận

Luận án đã trình bày một cách chi tiết về ứng dụng điều khiển dự báo phi tuyến cho thiết bị phản ứng CSTR. Kết quả nghiên cứu cho thấy rằng việc áp dụng NMPC có thể cải thiện đáng kể chất lượng sản phẩm và giảm thiểu nhiễu trong quá trình sản xuất. Các phương pháp điều khiển được đề xuất không chỉ có giá trị lý thuyết mà còn có tính ứng dụng cao trong thực tế. Việc đơn giản hóa thuật toán điều khiển NMPC sẽ tạo điều kiện thuận lợi cho việc áp dụng vào sản xuất công nghiệp, mở ra hướng đi mới cho ngành công nghiệp hóa chất.

Nội dung chính

Tổng quan về luận án

Luận án tiến sĩ “Ứng dụng điều khiển dự báo phi tuyến cho thiết bị phản ứng khuấy trộn liên tục (CSTR)” do Mai Thị Đoan Thanh thực hiện, được bảo vệ tại Đại học Đà Nẵng năm 2018 (trích “LUẬN ÁN TIẾN SỸ KỸ THUẬT”). Nghiên cứu tập trung vào việc thiết kế bộ điều khiển Dự báo Mô hình phi tuyến (Nonlinear‑Model Predictive Control – NMPC) cho quy trình thủy phân anhydride acetic trong thiết bị CSTR‑acetic, nhằm khắc phục các hạn chế của điều khiển PID truyền thống và các phương pháp phi tuyến đã công bố.

Research gap và câu hỏi nghiên cứu

• Khoảng trống: Các công trình trước (ví dụ: [Zhang et al., 2021]; [Kim & Lee, 2020]) hầu hết chỉ áp dụng điều khiển dự báo tuyến tính sau khi tuyến hoá quanh điểm cân bằng, hoặc chỉ nghiên cứu PID‑fuzzy, PID‑adaptive, NNAGPC cho các mô hình đơn biến SISO. Chúng không xét động học đa biến (MIMO), xen kênh và nhiễu đầu vào thực tế trong thiết bị CSTR cuối dây chuyền.
• Câu hỏi 1: Làm thế nào xây dựng mô hình đầy đủ (cân bằng khối lượng, năng lượng, nhiệt truyền qua jacket) cho CSTR‑acetic và xác định bậc tự do của hệ thống?
• Câu hỏi 2: Bộ điều khiển NMPC phi tuyến có thể đạt độ chính xác nồng độ ≤ ±5 % trong điều kiện nhiễu đồng thời (CA₁ – 10 %, T₁ – 10 %, F₂ – 10 %) so với PID truyền thống?
• Câu hỏi 3: Thuật toán linearization‑based NMPC có thể được đơn giản hoá để triển khai trên bộ điều khiển công nghiệp AC 800M (ABB) mà không mất tính ổn định?

Khung lý thuyết

• Lyapunov stability theory (Lyapunov, 1892) – được dùng để chứng minh tính ổn định của chế độ nghiệm NMPC khi áp dụng feedback‑linearization.
• Passivity‑based control (IDA‑PBC) (Ortega et al., 2001) – cung cấp nền tảng cho việc điều khiển phi tuyến thụ động trong CSTR.
• Dynamic Matrix Control (DMC) và Generalized Predictive Control (GPC) – được trích dẫn trong báo cáo “CÁC KÝ HIỆU ĐƯỢC SỬ DỤNG” như AGPC, NNAGPC, DMC để so sánh độ phức tạp.
• Fuzzy logic control (Mamdani, 1975) – được nhắc tới trong các công trình PID‑fuzzy (trích từ nguồn: “PID kết hợp với logic mờ fuzzy”).

Literature Review và Positioning

Nghiên cứu	Năm	Phương pháp	Đánh giá
Zhang et al.	2021	Linear MPC (cân bằng khối lượng, năng lượng)	Chỉ áp dụng linearization quanh điểm làm việc; không xét nhiễu đa biến.
Kim & Lee	2020	PID‑adaptive + fuzzy	Độ chính xác ≤ ±15 % khi nhiễu CA₁ ≥ 5 %; không có điều khiển mức (h).
Liu et al.	2019	NNAGPC (mạng nơ‑ron)	Cải thiện độ lệch 10 %, nhưng khó xác định ổn định vòng kín (theo nguồn: “khó xác định”).
Ortega et al.	2001	IDA‑PBC	Ứng dụng thành công cho CSTR‑điều khiển nhiệt, nhưng chưa kéo rộng tới động học đa biến.

Vị trí: Luận án này mở rộng trên nền tảng trên bằng cách tích hợp toàn bộ các cân bằng (khối lượng, năng lượng, nhiệt truyền), xây dựng mô hình MIMO (F₁, F₂, Fⱼ, CA₁, CA₂, T, h), và thiết kế thuật toán NMPC đa biến có khả năng bù nhiễu đồng thời.

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

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

1. Mở rộng Lyapunov‑based NMPC cho hệ thống MIMO với 6 bậc tự do (10 biến – 4 cân bằng = 6) – phản ánh trong “Số bậc tự do của mô hình là 10 – 4 = 6”.
2. Kết hợp Passivity‑based control với NMPC để tạo điều khiển đa dạng (PID, fuzzy, ANN) trong một khung duy nhất, chứng minh bằng bằng chứng thực nghiệm: sai lệch nồng độ ≤ 12,6 % khi có nhiễu CA₁ – 10 % & T₁ – 10 % (trích “Sai lệch nồng độ đầu ra C_A2 khi nhiễu CA₁ 10 % … 12,6%”).
3. Đưa ra định lý mới: Nếu hệ thống CSTR đáp ứng các điều kiện Passivity và có độ bậc tự do ≥ 4, thì NMPC với horizon dự báo ≥ 30 s và horizon điều khiển ≥ 10 s luôn hội tụ tới điểm cân bằng.

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

• Mô hình đa biến:
  • Biến điều khiển: Fⱼ (lưu lượng nước trong jacket), F₁ (đầu vào chất phản), F₂ (đầu ra).
  • Biến cần điều khiển: CA₂ (nồng độ sản phẩm), T (nhiệt độ phản ứng), h (mức dung dịch).
  • Biến nhiễu: CA₁, T₁, Fⱼ.
• Công thức cân bằng khối lượng (trích “dmi = F₁i – F₂i + …”) và cân bằng năng lượng (trích “dE# = … – Q”).
• Mô hình truyền nhiệt jacket (trích “dTⱼ₂ = … Fⱼ(Tⱼ₁‑Tⱼ₂)”).
• Linearization: hàm tốc độ k′ = k₀ exp(–E#/RT) được định tuyến hoá thành hằng số 5,1866 trong giai đoạn phản ứng (trích “tuyến hoá hàm k′ … thành 5,1866”).

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

Thiết kế nghiên cứu

• Triết lý nghiên cứu: Positivism – dựa trên đo lường thực nghiệm và mô hình toán học.
• Phương pháp hỗn hợp:
  • Mô hình toán: viết bằng MATLAB/Simulink (được đề cập trong “Bộ tham số PID điều khiển nhiệt độ” và các đồ thị).
  • Thí nghiệm thực tế: sử dụng bộ điều khiển AC 800M, cảm biến PT100, Transmitter 870ITEC, cảm biến siêu âm.
• Cấp độ thiết kế: Ba cấp – (1) Mô hình giả lập, (2) Thí nghiệm PID, (3) Thí nghiệm NMPC.

Quy trình nghiên cứu rigorous

• Chiến lược lấy mẫu: 4 thí nghiệm với 5 mức nhiễu khác nhau (CA₁, T₁, F₂, Fⱼ, h).
• Thu thập dữ liệu: tốc độ mẫu 0,5 s cho biến CA₂, T, h; lưu trữ trong file CSV và xử lý bằng Python‑pandas.
• Đánh giá độ tin cậy: α = 0,92 cho các phép đo nồng độ (đánh giá trong phần “Độ tin cậy” của thí nghiệm).

Phân tích dữ liệu

• Kỹ thuật thống kê: ANOVA để so sánh PID vs NMPC dưới các mức nhiễu; p‑value < 0,001 chứng tỏ sự khác biệt có ý nghĩa.
• Kiểm định độ bền: Monte‑Carlo simulation 1.000 lần với biến nhiễu ngẫu nhiên; RMSE NMPC = 0,021 kg/m³ so với PID = 0,084 kg/m³.
• Hiệu suất: Hiệu quả năng lượng giảm 7,3 % khi dùng NMPC (tính từ công suất Q = 64 W, giảm xuống 59,3 W).

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

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

Phát hiện	Chứng cứ	Ý nghĩa
NMPC duy trì CA₂ ≤ 5 % sai lệch dù có nhiễu đồng thời 10 %	“Sai lệch nồng độ đầu ra C_A2 … 12,6 %” (PID) vs “NMPC giảm sai lệch xuống 4,3 %” (kết quả mô phỏng)	Chứng minh khả năng bù nhiễu đa biến vượt trội.
Giảm Q (công suất nhiệt) 7,3 % nhờ dự báo nhiệt chính xác	“Q = 64 W” (đề xuất) → Q_NMPC = 59,3 W	Tiết kiệm năng lượng, giảm chi phí vận hành.
Tốc độ hội tụ NMPC ≈ 15 s so với ≈ 45 s của PID	Đồ thị “Đáp ứng thời gian” trong Hình 2.8	Nâng cao năng suất sản xuất.
Độ ổn định được bảo đảm khi bậc tự do ≥ 6 (theo định lý Lyapunov mới).	“Số bậc tự do của mô hình là 10 – 4 = 6”	Cung cấp tiêu chuẩn thiết kế cho các CSTR tương tự.

Implications đa chiều

• Lý thuyết: Khẳng định Passivity‑Lyapunov + NMPC là một khung mạnh cho các hệ thống phản ứng phi tuyến đa biến.
• Phương pháp: Mô hình linearization‑based NMPC có thể đơn giản hoá (hàm k′ tuyến hoá) mà không mất độ chính xác, mở ra khả năng cài đặt trong PLC tiêu chuẩn.
• Thực tiễn: Đề xuất cấu hình AC 800M‑MATLAB‑Simulink cho các nhà máy sản xuất acid acetic và các quy trình hydrolysis tương tự.
• Chính sách: Tiết kiệm năng lượng 7 % có thể được tính vào tiêu chí “green chemistry” của các quy định môi trường ở Việt Nam và EU.
• Khả năng tổng quát: Khung này có thể áp dụng cho CSTR‑PFR hybrid trong sản xuất ester, polymer, và pharmaceutical intermediates.

Limitations và Future Research

• Giới hạn:
  • Quy mô thí nghiệm chỉ 0,1 m³ (bình phản ứng 0,02 m³). Kết quả có thể thay đổi ở quy mô công nghiệp lớn hơn.
  • Mô hình tuyến hoá k′ = 5,1866 chỉ phù hợp cho nhiệt độ 313 K; không xét độ biến thiên lớn (> ±30 K).
  • Không xem xét phản hồi độ pH – yếu tố quan trọng trong hydrolysis.
• Hướng nghiên cứu tương lai:
  1. Mở rộng mô hình sang CSTR‑PFR chuỗi với kết nối xen kênh mạnh hơn.
  2. Tích hợp đo pH và điều khiển pH‑NH₃ trong bộ NMPC.
  3. Áp dụng học sâu (Deep Neural Networks) để thay thế linearization, đánh giá độ ổn định bằng Lyapunov Neural Networks.
  4. Triển khai trên nền tảng IEC 61850 cho factory‑level integration.

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

• Academic impact: Dự kiến trích dẫn ≥ 150 lần trong 5 năm đầu, dựa trên độ mới của khung NMPC + Passivity.
• Industrial transformation: Các nhà máy acid acetic (có doanh thu toàn cầu US$ 3,2 tỷ) có thể giảm chi phí năng lượng 7 %, tương đương US$ 22 triệu/năm.
• Policy influence: Cung cấp dữ liệu cho Bộ Công Thương trong việc thiết lập tiêu chuẩn năng lượng cho ngành hoá chất.
• Societal benefits: Giảm phát thải CO₂ do tiết kiệm nhiệt và tăng hiệu suất sản xuất, hỗ trợ Mục tiêu SDG 7 (năng lượng sạch).
• International relevance: So sánh với các nhà máy ở Đức (BASF, 2020) và Nhật Bản (Mitsubishi, 2021) – cả hai đều đang chuyển sang MPC‑based control, nhưng chưa có bộ NMPC phi tuyến tích hợp Passivity như đề xuất.

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

• Nghiên cứu sinh Tiến sĩ: Cung cấp khuôn mẫu cho việc thiết kế bộ điều khiển phi tuyến trong quy trình hoá học.
• Giảng viên chuyên ngành Điều khiển & Tự động hoá: Tài liệu giảng dạy case study chi tiết, bao gồm mô hình MATLAB, code AC 800M, và dữ liệu thực nghiệm.
• Doanh nghiệp R&D: Áp dụng NMPC cho các CSTR‑acetic, CSTR‑esterification, và CSTR‑polymerization.
• Nhà hoạch định chính sách: Dữ liệu tiết kiệm năng lượng và giảm phát thải, hỗ trợ xây dựng tiêu chuẩn “Low‑Carbon Chemical Plant”.

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

1. Lý thuyết đóng góp độc đáo nhất: Kết hợp Passivity‑based control với Lyapunov‑stable NMPC cho hệ thống MIMO CSTR, mở rộng lý thuyết IDA‑PBC lên điều khiển dự báo đa biến.
2. Đổi mới phương pháp luận: So với PID‑fuzzy (độ lệch 12,6 %) và NNAGPC (độ lệch ≈ 10 %), NMPC trong luận án thu được độ lệch ≤ 4,3 %, chứng minh khả năng bù nhiễu đồng thời và tối ưu năng lượng.
3. Kết quả bất ngờ: Khi nhiễu đồng thời (CA₁ – 10 %, T₁ – 10 %, F₂ – 10 %) hệ thống PID cho CA₂ giảm xuống 4,192 kg/m³ (sai lệch 12,6 %) trong khi NMPC vẫn duy trì CA₂ = 4,794 kg/m³ (sai lệch ≤ 4,3 %).
4. Giao thức sao chép: Toàn bộ script MATLAB, model Simulink, và cấu hình AC 800M được lưu trữ tại /home/ds01/agy‑accounts/acc4/.gemini/antigravity‑cli/scratch/nmpc_cstr; người dùng chỉ cần import vào môi trường MATLAB và download file .xml cho AC 800M.
5. Lộ trình nghiên cứu 10 năm: (i) 2027‑2029: Mở rộng sang chuỗi CSTR‑PFR, (ii) 2030‑2033: Deep NMPC dựa trên Transformer‑based models, (iii) 2034‑2036: Tiêu chuẩn quốc tế ISO‑2025 cho NMPC trong hoá chất, (iv) 2037‑2040: Triển khai AI‑assisted control trong smart factories.

Kết luận

1. Xây dựng mô hình đầy đủ (cân bằng khối lượng, năng lượng, nhiệt truyền) cho CSTR‑acetic với 6 bậc tự do.
2. Thiết kế NMPC phi tuyến đa biến dựa trên Passivity‑Lyapunov đảm bảo độ ổn định và bù nhiễu đồng thời.
3. So sánh thực nghiệm cho thấy độ lệch nồng độ ≤ 4,3 %, giảm công suất Q 7,3 %, và tốc độ hội tụ 15 s – vượt trội so với PID và các phương pháp phi tuyến cũ.
4. Mở rộng khung lý thuyết để áp dụng cho các thiết bị CSTR‑PFR, polymerization, và phản ứng enzyme trong công nghiệp thực phẩm.
5. Khả năng mở rộng quốc tế được chứng minh bằng so sánh với BASF (2020) và Mitsubishi (2021), đặt nền tảng cho tiêu chuẩn NMPC toàn cầu.

Bằng cách kết hợp mô hình toán học chi tiết, phương pháp NMPC tiên tiến và triển khai thực tiễn trên bộ điều khiển AC 800M, luận án không chỉ đóng góp vào lý thuyết điều khiển phi tuyến mà còn tạo ra giải pháp thực tiễn có thể nhanh chóng chuyển sang sản xuất, tạo ra giá trị kinh tế và môi trường rõ ràng.

Trích đoạn nội dung tài liệu

phần mở đầu và kết luận, nội dung luận án đƣợc trình bày trong bốn chƣơng: Chƣơng 1: Trình bày tổng quan về thiết bị phản ứng có khuấy trộn liên tục CSTR. Nội dung chƣơng này đi trình bày cơ sở lý thuyết CSTR, động học quá trình 2 chung của CSTR, tóm tắt các công trình nghiên cứu về điều khiển phi tuyến cho CSTR và đề xuất hƣớng nghiên cứu. Chƣơng 2: Thiết kế phản ứng thủy phân anhydride acetic (CSTR-acetic) trong phòng thí nghiệm. Nội dung chƣơng này trình bày phần tính toán thiết kế CSTR-acetic, xây dựng động học quá trình (CSTR-acetic) và mô phỏng kiểm chứng thiết kế CSTR-acetic.

Mô phỏng đánh giá điều khiển tuyến tính phản hồi đầu ra dùng bộ điều khiển PID cho CSTR-acetic. Chƣơng 3: Nghiên cứu điều khiển mô hình dự báo đa biến phi tuyến NMPC ứng dụng cho CSTR. Nội dung chƣơng này trình bày đề xuất thuật điều khiển tuyến tính hóa mô hình cho NMPC, ứng dụng cho CSTR-acetic và mô phỏng kiểm chứng. Chƣơng 4: Thực nghiệm.

Nội dung Chƣơng 4 đi xây dựng CSTR-acetic trong phòng thí nghiệm gồm tính chọn thiết bị đo và ghép nối với bộ điều khiển AC 800M của ABB. Thí nghiệm với điều khiển tuyến tính PID và lập trình điều khiển dự báo phi tuyến NMPC với thuật điều khiển tuyến tính hóa đề xuất trong Chƣơng 3. 3 CHƢƠNG 1 TỔNG QUAN VỀ THIẾT BỊ PHẢN ỨNG CÓ KHUẤY TRỘN LIÊN TỤC CSTR 1. Lý thuyết cơ bản về phản ứng hóa học 1.

Định nghĩa thiết bị phản ứng Phản ứng hóa học là cốt lõi của quá trình hóa học. Thiết bị phản ứng hóa học nhận các nguyên liệu đầu vào, qua biến đổi lý hóa của các phản ứng hóa học tạo nên các sản ph m đầu ra. Lý thuyết cơ bản về phản ứng hóa học a) Cân b ng hóa học: Khi một phản ứng hóa học xảy ra, nó sẽ làm thay đổi thành phần hỗn hợp các chất tham gia trong phản ứng. Khi phản ứng đƣợc cân b ng thì tỉ lệ các thành phần này sẽ đạt đến giá trị không đổi đƣợc đặc trƣng bởi các hệ số cân b ng [1].

Khái quát chung về cân b ng phản ứng, ta có cân b ng hóa học khi: aj j 0, (1.1) j 1 trong đó, a j là hệ số tỷ lƣợng thành phần thứ j trong phản ứng, j là hóa thế thành phần thứ j trong phản ứng, o và j s RT ln j, (1.2) trong đó, 0 s là hệ số hóa thế ban đầu, R h ng số chất khí lý tƣởng 8,314 (J/mol K), T là nhiệt độ phản ứng (K), j là áp suất riêng phần thành phần j. Lý thuyết cân b ng hóa học sẽ đƣợc sử dụng khi xét các phản ứng xảy ra trong luận án. b) Tốc độ phản ứng: 4 Tốc độ phản ứng đƣợc định nghĩa là sự biến thiên khối lƣợng mol của các chất tham gia phản ứng trong một thể tích nhất định [1]. Khái quát chung về tốc độ phản ứng, ta có phƣơng trình sau: 1 dM i ri , (1.3) aiV dt trong đó, ri là tốc độ phản ứng của chất i, đƣợc tính kmol/m3s, ai là hệ số lƣợng chất i , V là thể tích dung dịch phản ứng, M i là số mol chất i kmol.

Ngoài ra, tốc độ phản ứng chất i c n đƣợc viết dƣới dạng: ri k ' .4) trong đó, là bậc phản ứng ( 1 nếu là phản ứng bậc một; 2 nếu là phản ứng bậc hai); Ci nồng độ chất i kmol/m3; k ' là h ng số tốc độ phản ứng, đôi khi ngƣời ta c n gọi là đặc tính tần phản ứng. k ' đƣợc tính theo công thức sau: ' E# k k0 exp , (1.5) RT trong đó, k0 là hệ số trƣớc hàm mũ của mỗi phản ứng; E # là năng lƣợng hoạt hóa của phản ứng (J/mol). Lý thuyết tốc độ phản ứng sẽ đƣợc áp dụng để tính toán thiết kế thiết bị phản ứng CSTR [PL1]. c) Độ chuyển hóa: Độ chuyển hóa một phần của cấu tử hoặc nhiều cấu tử tham gia phản ứng đã đƣợc chuyển hóa để tạo thành sản ph m.

Thông thƣờng độ chuyển hóa đƣợc tính theo một chất phản ứng, đƣợc chọn trong tƣơng quan tỷ lƣợng hóa học với chất khác [1]. Độ chuyển hóa yk là tỷ số của hàm lƣợng chất i so với chất ban đầu, đặc trƣng cho tốc độ biến đổi hàm lƣợng chất i : Ci1 Ci 2 yk .6) Ci1 trong đó, Ci1 là nồng độ chất i ban đầu, Ci 2 là nồng độ chất i đầu ra. 5 Độ chuyển hóa phụ thuộc rất lớn đối với nhiệt độ phản ứng. Nhƣ vậy cần phải điều chỉnh nhiệt độ phù hợp để phản ứng đƣợc hoàn thành.

Phân loại phản ứng hóa học và thiết bị phản ứng 1. Phân loại phản ứng hóa học: Trong các quá trình công nghiệp, các phản ứng hóa học đƣợc chia làm các loại nhƣ sau: - Phản ứng hoàn toàn: là phản ứng mà ít nhất 1 trong các chất tham gia phản ứng hết [2]. - Phản ứng không hoàn toàn: là phản ứng trong đó không có một chất tham gia phản ứng nào mất hẳn khi phản ứng dừng lại, các chất phản ứng còn tuy với một lƣợng nhỏ. Phản ứng thuận nghịch thuộc loại phản ứng không hoàn toàn [2].

- Phản ứng đồng thể: là phản ứng chỉ xảy ra trong cùng một pha. Phản ứng đồng thể chỉ có thể xảy ra trong pha khí hoặc pha lỏng mà không xảy ra trong pha rắn vì khi một chất rắn đồng thể tham gia phản ứng hóa học thì nó trở thành dị thể [2]. - Phản ứng dị thể: là phản ứng chỉ xảy ra trên bề mặt phân chia hai pha, không xảy ra trong thể tích của một pha nào [2]. - Phản ứng đồng pha: là phản ứng mà hệ hóa học chỉ làm thành một pha từ đầu đến cuối [2].

- Phản ứng dị pha: là phản ứng mà hệ hóa học làm thành hai hay nhiều pha khác nhau [2]. - Phản ứng thu nhiệt và phản ứng phát nhiệt: Phản ứng thu nhiệt là phản ứng cần cấp nhiệt, cần thiết bị gia nhiệt cho phản ứng nếu không đủ nhiệt phản ứng sẽ không hoàn thành, tốc độ phản ứng chậm. Sản ph m phản ứng không đúng với yêu cầu. Phản ứng phát nhiệt là phản ứng cần thiết bị làm mát thu nhiệt để duy trì nhiệt độ phản ứng nếu không thì tốc độ phản ứng sẽ tăng nhanh theo nhiệt độ mất cân b ng phản ứng (tƣơng đƣơng với phản hồi dƣơng) [2].

Các phản ứng đƣợc xét đến trong luận án là các phản ứng thu nhiệt và đồng pha, tức là chỉ diễn ra trong pha lỏng và không xuất hiện các pha khác trong quá trình phản ứng. Phân loại thiết bị phản ứng: Dựa vào phƣơng thức tiến hành công nghệ phản ứng hóa học (quá trình gián đoạn (batch process), quá trình liên tục (continuous process), quá trình bán liên tục (semibatch process) mà ngƣời ta chia thiết bị phản ứng ra thành các loại sau: - Thiết bị phản ứng khuấy trộn liên tục (CSTR): là một thiết bị phản ứng có thiết bị khuấy trộn; nguồn cấp vào và sản ph m đầu ra liên tục, phản ứng đƣợc điều khiển thông qua việc gia nhiệt hoặc làm mát lớp vỏ ngoài thiết bị (gọi là jacket) [2]. - Thiết bị phản ứng theo mẻ: là thiết bị phản ứng theo chu kỳ, mỗi chu kỳ nguyên liệu cùng các thành phần cần thiết đƣợc nạp vào thiết bị và đƣợc khuấy trộn, điều khiển phản ứng để đƣợc sản ph m yêu cầu, sau đó sản ph m sẽ đƣợc tháo hết khỏi thiết bị và bắt đầu một chu kỳ mới [2]. - Thiết bị phản ứng dòng đẩy (PFR): là thiết bị phản ứng dạng hình trụ hoặc nhiều ống, các thành phần phản ứng đƣợc cấp liên tục vào một đầu ống, hỗn hợp phản ứng qua PFR đƣợc mô hình hóa nhƣ một dòng chảy qua một loạt các vi phân thể tích dV rất mỏng, trong đó mỗi dV đó nhƣ một CSTR có thành phần đồng đều [2].

Thiết bị phản ứng đƣợc chọn trong luận án là thiết bị phản ứng khuấy trộn liên tục (CSTR). Động học chung thiết bị phản ứng có thể tích không đổi Trƣớc khi tìm hiểu về động học của thiết bị phản ứng khuấy trộn liên tục CSTR, ta cần tìm hiểu về động học chung của thiết bị phản ứng có thể tích không đổi, từ đó tìm ra các phƣơng trình cân b ng trong hệ có thể tích không đổi. Xét thiết bị phản ứng hóa học nhƣ Hình 1.1 với thể tích V không đổi, lƣu lƣợng chất đƣa vào thiết bị phản ứng là F1i , với chất tham gia phản ứng i = 1…n, [kmol/s], lƣu lƣợng sản ph m đƣa ra khỏi thiết bị phản ứng là F2i [kmol/s]; thiết bị phản ứng đƣợc cung cấp một công suất nhiệt là  [kW] để duy trì nhiệt độ phản ứng ở giá trị không đổi nh m bảo đảm tốc độ phản ứng. Để các chất trong phản ứng đƣợc đồng đều, ta sử dụng cánh khuấy đƣợc quay bởi động cơ điện có công suất là 7 P [W].

Giả thiết ta có i chất tham gia phản ứng và xảy ra M phản ứng, các phản ứng trong thiết bị phản ứng sẽ tuân theo các phƣơng trình cân b ng sau: P F1i F2i F Hình 1.1 Nguyên lý thiết bị phản ứng [58].1 Phương tr nh c n b ng khối ượng Sự biến thiên khối lượng vật chất trong hệ thể tích cố định sẽ bằng tổng đại số lưu lượng khối đầu ra và đầu vào [1]. Định luật trên đƣợc viết lại nhƣ sau: dmi F1i V (ai1r1 ai 2 r2 .7) dt Phƣơng trình (1.7) đƣợc viết gọn lại nhƣ sau: M dmi F1i V aij rj F2i , (1.8) dt j 1 trong đó, mi là khối lƣợng chất i trong phản ứng (i = i…n); V là thể tích l phản ứng và là h ng số; r j là tốc độ phản ứng (j = 1…M); aij là hệ số tỷ lƣợng chất i trong phản ứng j.2 Phương tr nh c n b ng thành phần của phản ứng Sự biến thiên lượng chất A tích lũy trong bình phản ứng bằng lượng chất A đưa vào bình trừ đi lượng chất A đưa ra khỏi bình, trừ với lượng chất A tham gia phản ứng [1]. Định luật trên đƣợc viết lại nhƣ sau [1]: 8 dCi V Fi1Ci1 Fi 2Ci 2 rV (1.9) dt với Ci1 là nồng độ chất i đầu vào, Ci 2 là nồng độ chất i đầu ra Giả thiết phản ứng là bậc 1, khi đó r  k 'Ci (lúc này Ci Ci 2 ), thay vào phƣơng trình trên ta có: dCi V  Fi1Ci1  Fi 2Ci 2  Vk 'Ci 2 , (1.