Tổng quan về luận án

Nghiên cứu điều khiển các hệ thống đa biến (Multi-Input Multi-Output - MIMO) giữ vai trò then chốt trong kỹ thuật điều khiển quá trình hiện đại, đặc biệt tại các phân xưởng công nghiệp quy mô lớn như lọc hóa dầu, chế biến khí, nhiệt điện và công nghệ hóa chất. Bản chất của hệ MIMO là sự tồn tại của các tương tác chéo (cross-coupling interactions) phức tạp giữa các biến điều khiển và biến trạng thái ngõ ra: một sự biến thiên tại một ngõ vào bất kỳ đều gây ra dao động không mong muốn trên nhiều ngõ ra còn lại. Thực tế này làm cho việc phân tích, tổng hợp và tối ưu hóa bộ điều khiển trở thành bài toán đa chiều phi tuyến đầy thách thức. Trong bối cảnh khoa học quốc tế, dù các kỹ thuật điều khiển tiên tiến như điều khiển dự báo mô hình (Model Predictive Control - MPC) hay điều khiển thông minh (Intelligent Control) phát triển mạnh mẽ, chúng vẫn bộc lộ hạn chế về chi phí tính toán, độ phức tạp thuật toán và khó khăn khi triển khai trên các phần cứng điều khiển công nghiệp khả trình (PLC/DCS). Do đó, xu hướng tận dụng cấu trúc điều khiển đơn biến (Single-Input Single-Output - SISO) kinh điển như PID thông qua mạng phân ly (decoupling network) nhằm tách hệ đa biến thành các vòng lặp đơn biến độc lập đang là hướng tiếp cận có tính đột phá và ứng dụng cao.

Khoảng trống nghiên cứu (research gap) cốt lõi được xác định xuất phát từ những giới hạn nội tại của các phương pháp phân ly truyền thống. Phân ly lý tưởng (Ideal Decoupling) tạo ra ma trận hàm truyền đường chéo hoàn hảo nhưng cấu trúc ma trận phân ly vô cùng phức tạp, thường chứa các phần tử không thực hiện được về mặt vật lý (non-causal) và cực kỳ nhạy cảm với sai số mô hình (Morari et al., 1989; Gagnon et al., 1998). Ngược lại, phân ly nghịch (Inverted Decoupling) khó hiện thực hóa trong cấu hình hồi tiếp mở rộng (Garrido et al., 2011). Đối với phân ly đơn giản hóa (Simplified Decoupling), các công trình trước đây hầu như chỉ giới hạn ở cấu hình 2 ngõ vào - 2 ngõ ra ($2 \times 2$) và thiếu một phương pháp thiết kế giải tích tổng quát cho hệ $n \times n$. Hơn nữa, việc hiệu chỉnh thông số bộ điều khiển PI/PID sau phân ly thường dựa vào các phương pháp dò tìm kinh nghiệm hoặc giải thuật lặp như Ziegler-Nichols cải biên (McAvoy, 1983), Biggest Log Modulus Tuning (BLT - Luyben, 1986) với tiêu chí $\text{BLT} = 2n\text{ dB}$, hay Sequential Auto-Tuning (SAT - Shen & Yu, 1994). Các phương pháp này dễ dẫn đến đáp ứng có độ vọt lố cao, dao động kéo dài và suy giảm chất lượng khi quá trình có thời gian trễ (dead-time) lớn.

Luận án tiến sĩ "Nghiên cứu thiết kế bộ điều khiển phân ly dùng cho các hệ thống đa biến" của NCS. Lê Linh (chuyên ngành Kỹ thuật Cơ khí, mã số 9520103, Trường Đại học Sư phạm Kỹ thuật TP. Hồ Chí Minh) đã giải quyết triệt để các hạn chế trên thông qua hệ thống câu hỏi và giả thuyết nghiên cứu tường minh:

  • Câu hỏi nghiên cứu 1 (Q1): Làm thế nào để xây dựng một quy luật giải tích tổng quát cho ma trận phân ly đơn giản hóa trên hệ thống đa biến bậc $n \times n$ dựa trên ma trận độ lợi tương quan động (Dynamic Relative Gain Array - DRGA)?

  • Câu hỏi nghiên cứu 2 (Q2): Làm thế nào để thiết lập công thức tường minh xác định các tham số bộ điều khiển PID theo cấu trúc điều khiển mô hình nội (Internal Model Control - IMC) cho các đối tượng có thời gian trễ mà không cần qua quy trình dò tìm số phức tạp?

  • Câu hỏi nghiên cứu 3 (Q3): Kỹ thuật xấp xỉ bậc cao Padé [3/2] và cấu trúc bộ lọc nối tiếp (series filters) có khả năng triệt tiêu tương tác chéo và giảm thiểu dao động tín hiệu điều khiển trên mô hình thực nghiệm thời gian thực như thế nào?

  • Giả thuyết nghiên cứu 1 (H1): Việc khai thác tỷ lệ giữa các đặc tính vòng hở ban đầu và các phần tử đường chéo của DRGA cho phép tổng quát hóa ma trận phân ly đơn giản hóa cho hệ $n \times n$ với các phần tử đường chéo bằng đơn vị ($d_{ii}(s) = 1$).

  • Giả thuyết nghiên cứu 2 (H2): Ứng dụng xấp xỉ Padé [3/2] cho khâu trễ $e^{-\theta s}$ kết hợp cấu trúc IMC và phương pháp đồng nhất hệ số chuỗi Maclaurin sẽ cung cấp công thức tường minh tính toán $K_c, \tau_I, \tau_D$, nâng cao đáng kể biên độ ổn định của hệ thống.

  • Giả thuyết nghiên cứu 3 (H3): Tích hợp bộ lọc bậc thấp nối tiếp $F_i(s)$ sẽ tối ưu hóa chỉ số sai số tích phân tuyệt đối (IAE) và giảm tổng biến thiên tín hiệu điều khiển (Total Variation - TV) cả trong bài toán bám điểm đặt (servomechanism) và kháng nhiễu tải (regulator).

Khung lý thuyết của công trình tích hợp chặt chẽ giữa Lý thuyết Điều khiển phân ly (Decoupling Control Theory), Lý thuyết Điều khiển mô hình nội (IMC Theory của Morari & Zafiriou), Lý thuyết Mảng độ lợi tương quan động (DRGA) và Kỹ thuật xấp xỉ hữu tỉ Padé. Phạm vi nghiên cứu bao quát từ kiểm chứng mô phỏng trên 3 hệ thống chuẩn quốc tế nổi tiếng bao gồm: Cột chưng cất Wood & Berry ($2 \times 2$), Cột chưng cất Vinante & Luyben ($2 \times 2$), Cột chưng cất Ogunnaike & Ray ($3 \times 3$) và hệ Wardle & Wood ($2 \times 2$), đến chế tạo và thực thi trên 2 mô hình thực nghiệm phần cứng: Hệ bồn nước đôi ($2 \times 2$) và Tháp chưng cất hỗn hợp Ethanol - Nước ($3 \times 3$). Nghiên cứu mang lại ý nghĩa khoa học và thực tiễn vượt bậc khi cung cấp công cụ toán học hoàn chỉnh giúp làm chủ công nghệ điều khiển các quá trình chiết tách công nghiệp.

Literature Review và Positioning

Nghiên cứu điều khiển đa biến trong nửa thế kỷ qua đã hình thành hai trường phái chính: điều khiển tập trung thuần túy (centralized control) và điều khiển phân ly kết hợp đa vòng lặp (decoupled decentralized multi-loop control).

Trường phái điều khiển tập trung thuần túy, tiêu biểu với công trình của Wang và các cộng sự (2000, 2002), tiếp cận bài toán trên miền tần số thông qua giải thuật đệ quy bình phương cực tiểu (Recursive Least Squares) nhằm xấp xỉ ma trận hàm truyền đạt. Dù đạt được thời gian đáp ứng nhanh trên mô hình lý thuyết, phương pháp tập trung thuần túy đòi hỏi cấu trúc thuật toán ma trận bậc cao, rất phức tạp khi cài đặt trên các thiết bị điều khiển biên (edge controllers) và thiếu tính linh hoạt khi cần chuyển đổi giữa chế độ vận hành bằng tay (Manual) và tự động (Auto) cho từng kênh riêng biệt.

Trường phái điều khiển phân ly kết hợp mạng đa vòng PI/PID được khởi xướng và phát triển mạnh mẽ nhờ tính trực quan và mức độ tương thích cao với hạ tầng công nghiệp. Khởi nguyên từ phương pháp hiệu chỉnh kinh điển của Ziegler & Nichols (1942, 1943), McAvoy (1983) đã phát triển quy luật điều khiển đa vòng lặp để giảm thiểu tương tác chéo. Tiếp đó, Luyben (1986) giới thiệu phương pháp Biggest Log Modulus Tuning (BLT), sử dụng hệ số tinh chỉnh $F \in [2, 5]$ để thỏa hiệp giữa độ ổn định và tốc độ đáp ứng, đưa biên độ vòng kín đạt đỉnh cộng hưởng cực đại $+2n\text{ dB}$ ($+4\text{ dB}$ cho hệ $2 \times 2$ và $+6\text{ dB}$ cho hệ $3 \times 3$). Tuy nhiên, BLT bộc lộ nhược điểm phản ứng chậm và đáp ứng kháng nhiễu kém. Nhằm khắc phục, Shen & Yu (1994) đề xuất phương pháp Sequential Auto-Tuning (SAT) dựa trên kỹ thuật rơ-le hồi tiếp (Åström & Hägglund, 1984), cho phép tinh chỉnh tuần tự từng vòng lặp. Dù SAT cải thiện tốc độ đáp ứng so với BLT, nó lại gây ra hiện tượng vọt lố lớn và dao động kéo dài khi hệ thống tồn tại thời gian trễ lớn hoặc có tính phi tuyến cao. Gần đây hơn, Truong & Lee (2010) đã đề xuất phương pháp thiết kế bộ điều khiển PI đa biến dựa trên khai triển chuỗi Maclaurin kết hợp chỉ số mảng độ lợi tương quan (RGA), mở ra hướng đi mới trong việc giải tích hóa tham số điều khiển.

Trong lý thuyết phân ly, cuộc tranh luận khoa học sâu sắc diễn ra giữa ba cấu trúc: Phân ly lý tưởng, Phân ly nghịch và Phân ly đơn giản hóa:

  • Gagnon và cộng sự (1998) cùng Garrido và các cộng sự (2011, 2013) khi tiến hành so sánh chuyên sâu trên các mô hình tháp chưng cất đã chỉ ra rằng: Phân ly lý tưởng tuy triệt tiêu hoàn toàn các kênh tương tác ngoài đường chéo nhưng ma trận phân ly $D(s)$ lại quá phức tạp, chứa các hàm truyền không nhân quả và cực kỳ nhạy cảm trước sai số mô hình.
  • Morari và các cộng sự (1985, 1989), thông qua định lý độ lợi nhỏ (Small Gain Theorem) và phân tích giá trị kỳ dị cấu trúc (Structured Singular Value - SSV), đã chứng minh rằng phân ly lý tưởng có độ bền vững (robust stability) thấp hơn đáng kể so với phân ly đơn giản hóa.
  • Shinskey (1988) khẳng định cấu trúc phân ly đơn giản hóa với các phần tử đường chéo bằng đơn vị ($D_{ii} = 1$) là cấu trúc tối ưu nhất để hiện thực hóa trong công nghiệp, vừa đảm bảo chuyển mạch trơn tru Manual/Auto, vừa duy trì độ bền vững cao trước các bất định tham số.

Vị trí của luận án được xác lập tại giao điểm đột phá: kế thừa ưu thế bền vững của phân ly đơn giản hóa theo quan điểm của Morari và Shinskey, khắc phục giới hạn chiều không gian $2 \times 2$ để tổng quát hóa cho hệ $n \times n$ qua DRGA, đồng thời giải quyết triệt để bài toán trễ pha của Luyben và Shen & Yu bằng kỹ thuật xấp xỉ Padé [3/2] và bộ lọc nối tiếp trong khung điều khiển IMC. Nghiên cứu tạo ra bước nhảy vọt so với các công trình quốc tế khi chuyển đổi từ quy trình dò tìm tham số lặp (BLT, SAT) sang công thức tính toán giải tích đơn bước tường minh.

Đó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 ba lý thuyết nền tảng trong kỹ thuật điều khiển tự động: Lý thuyết Điều khiển mô hình nội (IMC của Morari & Zafiriou), Lý thuyết Phân ly tương quan động (DRGA) và Lý thuyết Xấp xỉ hàm hữu tỉ (Rational Approximation Theory).

Mô hình khái niệm được xây dựng trên nguyên lý chuyển đổi tương đương: hệ thống đa biến phức tạp $P(s)$ kích thước $n \times n$ với các liên kết chéo mạnh được ghép nối tiếp với ma trận phân ly đơn giản hóa $D(s)$ để tạo thành một đối tượng tổng hợp $Q(s) = P(s)D(s)$. Ma trận $Q(s)$ có dạng gần đường chéo, cho phép phân rã hoàn toàn thành $n$ vòng lặp đơn biến độc lập $q_{ii}(s)$. Trên từng kênh đơn biến, bộ điều khiển IMC-PID kết hợp bộ lọc $F_i(s)$ được tổng hợp giải tích:

$$\text{Mô hình lý thuyết: } u(s) = D(s) \cdot m(s) = D(s) \cdot \text{diag}[G_{c1}(s)F_1(s), \dots, G_{cn}(s)F_n(s)] \cdot e(s)$$

Các định đề lý thuyết chính được xác lập bao gồm:

  • Định đề 1 (Proposition 1): Mọi phần tử ngoài đường chéo $d_{ij}(s)$ của ma trận phân ly đơn giản hóa $n \times n$ đều được xác định tường minh thông qua định thức con và tỷ số hàm truyền thành phần từ ma trận DRGA, thỏa mãn triệt tiêu tương tác tĩnh và động tại dải tần số làm việc.

  • Định đề 2 (Proposition 2): Khâu trễ vô tỉ $e^{-\theta s}$ được thay thế bằng hàm hữu tỉ Padé [3/2] chuẩn tắc:

    $$e^{-\theta s} \approx \frac{1 - \frac{3}{5}\theta s + \frac{3}{20}\theta^2 s^2 - \frac{1}{60}\theta^3 s^3}{1 + \frac{2}{5}\theta s + \frac{1}{20}\theta^2 s^2}$$

    Kỹ thuật này bảo toàn đặc tính động học pha và triệt tiêu sai số cắt bậc thấp của xấp xỉ Padé [1/1] hoặc Taylor truyền thống.

  • Định đề 3 (Proposition 3): Áp dụng khai triển chuỗi Maclaurin cho ma trận độ nhạy $S(s) = [I + Q(s)G_c(s)]^{-1}$ và hàm truyền vòng kín $H(s)$, toàn bộ các tham số $K_c, \tau_I, \tau_D$ của bộ điều khiển PID được rút ra dưới dạng công thức đại số đóng phụ thuộc vào độ lợi quá trình $K$, hằng số thời gian $\tau$, thời gian trễ $\theta$ và tham số điều chỉnh vòng kín $\lambda$.

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

Khung phân tích của luận án tích hợp phương pháp đồng nhất hệ số (Method of Undetermined Coefficients) để xử lý tính nhân quả và ổn định của ma trận phân ly $D(s)$. Trong thực tế, các phần tử $d_{ij}(s) = -P_{ij}(s)/P_{ii}(s)$ thường xuất hiện khâu dự báo thời gian $e^{+\Delta \theta s}$ (nếu thời gian trễ của kênh chéo nhỏ hơn kênh chính), làm cho bộ phân ly không thể thực thi vật lý. Khung phân tích giải quyết vấn đề này bằng cách bổ sung độ trễ nhân tạo vào các kênh thích hợp hoặc xấp xỉ đồng nhất hệ số bậc động học tương đương, đảm bảo toàn bộ các phần tử $d_{ij}(s)$ đều là các hàm truyền hữu tỉ thực hiện được (causal transfer functions).

Điều kiện biên (boundary conditions) của khung phân tích được xác định rõ ràng:

  1. Áp dụng cho các quá trình đa biến có mô hình biểu diễn được dưới dạng hàm truyền FOPDT (bậc một có trễ), SOPTD (bậc hai có trễ) hoặc hệ có điểm không ngoài vòng tròn đơn vị (Non-minimum Phase - NMP).
  2. Ma trận độ lợi tĩnh $G_0$ không suy biến ($\det(G_0) \neq 0$).
  3. Tỷ số trễ trên hằng số thời gian $\theta/\tau \le 2.0$, đảm bảo tính hội tụ của khai triển chuỗi và độ bền vững của cấu trúc IMC.

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

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên nền tảng triết lý thực chứng thực nghiệm (Positivism/Deductive Engineering Science), kết hợp chặt chẽ giữa suy luận diễn dịch toán học thuần túy và kiểm chứng thực nghiệm nghiêm ngặt. Thiết kế nghiên cứu phân tầng 3 cấp độ (Multi-level design):

  • Cấp độ 1 (Modeling & Formulation): Giải tích hóa toán học ma trận phân ly $n \times n$ và thiết lập công thức thông số PID-IMC bằng kỹ thuật Padé [3/2].
  • Cấp độ 2 (Numerical Benchmarking): Mô phỏng số động học trên MATLAB/Simulink đối với các mô hình kinh điển thế giới để so sánh định lượng với các thuật toán BLT và SAT.
  • Cấp độ 3 (Hardware-in-the-Loop Implementation): Cài đặt thuật toán điều khiển thời gian thực trên 2 hệ thống phần cứng thực tế qua giao thức thu thập dữ liệu chuyên dụng.

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu thực nghiệm tuân thủ các chuẩn mực đo lường nghiêm ngặt:

  • Nhận dạng mô hình thực nghiệm (System Identification): Đối với các hệ thống vật lý chưa có phương trình giải tích, luận án áp dụng phương pháp nhận dạng hộp đen (Black-box identification) sử dụng đáp ứng nấc (step response) và tín hiệu kích thích chuỗi nhị phân giả ngẫu nhiên (PRBS). Dữ liệu ngõ vào - ngõ ra được thu thập qua card DAQ chuyên dụng của hãng National Instruments (NI PCI/USB DAQ).
  • Thuật toán xử lý và đồng bộ thời gian thực: Chương trình điều khiển được nhúng trực tiếp trong môi trường MATLAB/Simulink chạy ở chế độ Real-Time Windows Target (RTWT) với chu kỳ lấy mẫu cố định $T_s = 0.1\text{ s}$, đảm bảo không có hiện tượng trượt thời gian (timing jitter).
  • Tam giác hóa kiểm chứng (Triangulation): Kết quả được đối chiếu chéo giữa 3 phương thức: Phương trình giải tích toán học $\leftrightarrow$ Mô phỏng số Simulink $\leftrightarrow$ Đo đạc thực nghiệm trên phần cứng.
  • Độ tin cậy và độ giá trị (Validity & Reliability): Tính giá trị cấu trúc được bảo đảm qua việc sử dụng các chỉ số chuẩn hóa quốc tế: Sai số tích phân tuyệt đối (Integral Absolute Error - IAE) và Tổng biến thiên tín hiệu điều khiển (Total Variation - TV). Độ tin cậy được kiểm chứng qua việc lặp lại thí nghiệm 5 lần ở các chế độ tải khác nhau với độ lệch chuẩn sai số $< 2.5%$.

Data và phân tích

Các mô hình toán học chuẩn và dữ liệu thực nghiệm được khảo sát chi tiết:

  1. Mô hình Cột chưng cất Wood & Berry (WB - $2 \times 2$): Tách hỗn hợp Methanol - Nước, ma trận hàm truyền FOPDT:

    $$P_{WB}(s) = \begin{bmatrix} \frac{12.8 e^{-s}}{16.7s + 1} & \frac{-18.9 e^{-3s}}{21s + 1} \ \frac{6.6 e^{-7s}}{10.9s + 1} & \frac{-19.4 e^{-3s}}{14.4s + 1} \end{bmatrix}$$

  2. Mô hình Cột chưng cất Vinante & Luyben (VL - $2 \times 2$):

    $$P_{VL}(s) = \begin{bmatrix} \frac{-2.2 e^{-s}}{7s + 1} & \frac{1.3 e^{-0.3s}}{7s + 1} \ \frac{-2.8 e^{-1.8s}}{9.5s + 1} & \frac{4.3 e^{-0.35s}}{9.2s + 1} \end{bmatrix}$$

  3. Mô hình Cột chưng cất Ogunnaike & Ray (OR - $3 \times 3$): Quá trình 3 ngõ vào, 3 ngõ ra phức tạp với các khâu trễ chéo dao động từ $1.0\text{ s}$ đến $3.5\text{ s}$.

  4. Mô hình thực nghiệm Bồn nước đôi ($2 \times 2$): Điều khiển mức nước $h_1, h_2 \in [0, 300]\text{ mm}$ thông qua điều khiển điện áp 2 bơm cấp liệu.

  5. Mô hình thực nghiệm Tháp chưng cất Ethanol - Nước ($3 \times 3$): Điều khiển 3 biến công nghệ chủ chốt: Nhiệt độ bồn đáy ($y_1 \in [70, 100]^\circ\text{C}$), Nhiệt độ dòng cấp liệu ($y_2 \in [40, 80]^\circ\text{C}$) và Lưu lượng dòng cấp liệu ($y_3 \in [0, 50]\text{ lít/giờ}$).

Các tiêu chuẩn đánh giá định lượng được tính toán trực tiếp:

$$\text{IAE} = \int_0^\infty |e(t)| dt = \int_0^\infty |r(t) - y(t)| dt$$

$$\text{TV} = \sum_{k=1}^\infty |u(k+1) - u(k)|$$

Kiểm tra độ bền vững (Robustness Checks) được thực hiện bằng cách đưa vào sai số mô hình $\pm 20%$ trên hằng số thời gian trễ $\theta$ và độ lợi tĩnh $K$, đồng thời tác động nhiễu bậc thang trên các kênh điều khiển để thẩm định khả năng hồi phục trạng thái xác lập.

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

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

Trích dẫn trực tiếp từ văn bản luận án xác nhận các luận điểm cốt lõi:

"Phân ly đơn giản hóa với ưu điểm về sự đơn giản của ma trận phân ly, thành phần đường chéo bằng đơn vị, và sự bền vững khi có sai số mô hình là sự lựa chọn thích hợp nhất cho các ứng dụng thực tế." (Lê Linh, 2020, tr. iv)

"Xấp xỉ bậc cao cho khâu trễ sử dụng công thức xấp xỉ Padé 3/2 nhằm cải thiện độ chính xác và đặc tính động của hệ thống." (Lê Linh, 2020, tr. vi)

"Áp dụng các bộ lọc nối tiếp với bộ điều khiển tương ứng cho các quá trình khác nhau nhằm cải thiện khả năng đáp ứng của hệ thống." (Lê Linh, 2020, tr. vii)

Các phát hiện thực nghiệm và mô phỏng đột phá bao gồm:

  1. Vượt trội toàn diện về chỉ số chất lượng bám điểm đặt và kháng nhiễu: Trên mô hình Cột chưng cất Wood & Berry (WB), phương pháp đề xuất triệt tiêu hoàn toàn hiện tượng dao động cộng hưởng kéo dài vốn xuất hiện ở phương pháp McAvoy và Z-N. So với BLT ($F = 2.45, K_{c1} = 0.375, \tau_{I1} = 8.29$) và SAT ($\beta = 0.4, K_{c1} = 0.98, \tau_{I1} = 6.4$), bộ điều khiển đề xuất giúp giảm chỉ số IAE từ $28%$ đến $45%$, đồng thời loại bỏ hoàn toàn độ vọt lố ($M_p = 0%$) trên kênh $y_1$ và giảm thời gian xác lập $t_s$ xuống dưới $35\text{ phút}$ (so với $> 80\text{ phút}$ của BLT).
  2. Triệt tiêu tương tác động học đa biến (Dynamic Decoupling Efficiency): Khi có sự thay đổi giá trị đặt bước nhảy tại vòng lặp 1 ($\Delta r_1 = 1.0$), độ lệch quá độ của ngõ ra vòng lặp 2 ($\Delta y_2$) ở phương pháp đề xuất tiệm cận mức 0 ($\Delta y_2 < 0.02$), trong khi ở phương pháp BLT và SAT, ngõ ra $y_2$ bị dao động cưỡng bức với biên độ lệch lên tới $0.25 - 0.38$ đơn vị.
  3. Hiệu năng vượt trội của xấp xỉ Padé [3/2]: Việc sử dụng xấp xỉ Padé [3/2] thay thế cho Padé [1/1] giúp bù pha chính xác tại dải tần số cắt biên độ $\omega_c$, ngăn ngừa hiện tượng suy giảm biên dự trữ pha (Phase Margin), từ đó cho phép tăng độ lợi tỷ lệ $K_c$ mà không gây mất ổn định vòng kín.
  4. Giảm thiểu độ mệt mỏi cơ cấu chấp hành nhờ bộ lọc nối tiếp: Phân tích chỉ số tổng biến thiên tín hiệu điều khiển TV cho thấy việc tích hợp bộ lọc nối tiếp $F_i(s)$ giúp giảm chỉ số TV từ $35%$ đến $52%$ so với cấu hình PID tiêu chuẩn, loại bỏ hiện tượng giật cục (chattering) của van điều khiển và máy bơm định lượng.
  5. Độ chính xác và tính khả thi trên hệ thống thực nghiệm: Kết quả điều khiển thời gian thực trên Tháp chưng cất Ethanol - Nước ($3 \times 3$) cho thấy hệ thống bám chính xác giá trị đặt nhiệt độ bồn đáy ($85^\circ\text{C} \pm 0.5^\circ\text{C}$), nhiệt độ cấp liệu ($65^\circ\text{C} \pm 0.8^\circ\text{C}$) và lưu lượng cấp liệu ($30\text{ l/h} \pm 0.2\text{ l/h}$), giúp nồng độ cồn sản phẩm đỉnh đạt độ tinh khiết yêu cầu $> 85%$ thể tích.

Implications đa chiều

  • Về mặt lý thuyết: Hoàn thiện lý thuyết điều khiển phân ly giải tích cho hệ đa biến bậc cao $n \times n$; cung cấp cầu nối toán học chặt chẽ giữa cấu trúc điều khiển mô hình nội IMC và bộ điều khiển PID công nghiệp có khâu lọc.
  • Về mặt phương pháp luận: Thiết lập quy trình thiết kế bộ điều khiển đơn bước (one-step analytical design workflow), loại bỏ hoàn toàn nhu cầu chạy các thuật toán tối ưu hóa số phi tuyến phức tạp hay dò tìm thử - sai tốn kém thời gian.
  • Về mặt ứng dụng thực tiễn: Cho phép các kỹ sư tự động hóa tại các nhà máy triển khai trực tiếp các thuật toán điều khiển phân ly tiên tiến trên các hệ thống DCS (Yokogawa Centum VP, ABB Ability, Siemens PCS7) và PLC thông dụng bằng các khối hàm FBD/ST tiêu chuẩn mà không cần mua các gói phần mềm MPC đắt đỏ.
  • Về mặt chính sách và công nghệ: Đóng góp giải pháp kỹ thuật tự chủ, nâng cao năng lực nội địa hóa công nghệ điều khiển trong các ngành công nghiệp mũi nhọn như dầu khí, hóa chất và năng lượng tái tạo, giảm phụ thuộc vào các giải pháp trọn gói từ nước ngoài.

Limitations và Future Research

Nhằm duy trì tính khách quan khoa học, luận án thẳng thắn chỉ ra 4 giới hạn nghiên cứu:

  1. Quy mô ma trận thực nghiệm: Các kiểm chứng thực nghiệm phần cứng mới dừng lại ở cấu hình $3 \times 3$ (tháp chưng cất). Với các hệ thống siêu đa biến ($n \ge 5$), việc tính toán ma trận phân ly giải tích có thể gặp hiện tượng phình to bậc hàm truyền (transfer function order inflation).
  2. Tính phi tuyến tại vùng làm việc mở rộng: Bộ điều khiển được thiết kế dựa trên mô hình tuyến tính hóa quanh điểm làm việc danh định. Khi hệ thống dịch chuyển sang các chế độ vận hành phi tuyến cực đoan, chất lượng phân ly có thể bị suy giảm cục bộ.
  3. Hiện tượng bão hòa cơ cấu chấp hành: Thiết kế giải tích trong miền tần số chưa tích hợp trực tiếp ràng buộc bất đẳng thức về giới hạn hành trình và tốc độ của van (actuator saturation constraints) như trong bộ điều khiển MPC.
  4. Hệ thống có điểm không ngoài vòng tròn đơn vị (NMP) phức tạp: Khi các phần tử truyền chéo chứa các nghiệm zero không ổn định đồng thời với trễ lớn, ma trận phân ly đơn giản hóa đòi hỏi cấu trúc bù trễ nhân tạo làm chậm tốc độ đáp ứng của vòng kín.

Chương trình nghiên cứu tiếp theo (Future Research Agenda) bao gồm:

  • Mở rộng thuật toán phân ly thích nghi (Adaptive Simplified Decoupling) tự động cập nhật ma trận $D(s)$ theo thời gian thực dựa trên các kỹ thuật nhận dạng tham số trực tuyến.
  • Tích hợp khâu chống bão hòa tích phân đa biến (Multivariable Anti-Windup) vào cấu trúc IMC-PID sau phân ly.
  • Nghiên cứu mở rộng cấu trúc phân ly cho các hệ thống đa biến không vuông (Non-square MIMO systems, số ngõ vào $m \ne$ số ngõ ra $n$).
  • Ứng dụng trí tuệ nhân tạo (AI/Neural Networks) để bù sai số phi tuyến cho các thành phần ngoài đường chéo của ma trận phân ly.

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

  • Tác động học thuật: Công trình đóng góp nền tảng phương pháp luận vững chắc cho cộng đồng nghiên cứu điều khiển quá trình trong nước và quốc tế. Phương pháp giải tích hóa thông số IMC-PID với Padé [3/2] mở ra hướng nghiên cứu mới trong việc tối ưu hóa hệ thống có trễ, với tiềm năng trích dẫn cao trên các tạp chí chuyên ngành uy tín (IEEE Transactions on Control Systems Technology, Control Engineering Practice, ISA Transactions).
  • Chuyển đổi công nghiệp: Ứng dụng thành công trên tháp chưng cất Ethanol - Nước mở ra khả năng chuyển giao công nghệ trực tiếp cho các nhà máy sản xuất cồn sinh học, nhà máy chiết xuất dược liệu và các phân xưởng tách phân đoạn dầu mỏ, giúp tiết kiệm từ $8 - 15%$ năng lượng tiêu thụ tại nồi đun đáy tháp nhờ kiểm soát chính xác nhiệt độ và giảm dao động quá độ.
  • Lợi ích xã hội và môi trường: Nâng cao hiệu suất phân tách hóa chất giúp giảm thiểu phát thải độc hại ra môi trường, tối ưu hóa quá trình sử dụng tài nguyên và đảm bảo an toàn vận hành, phòng chống sự cố quá nhiệt, quá áp trong các thiết bị áp lực công nghiệp.

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

  • Nghiên cứu sinh và Học viên cao học: Tiếp cận khung phân tích giải tích hoàn chỉnh, phương pháp nhận dạng hộp đen và quy trình mô phỏng - thực nghiệm chuẩn mực cho hệ đa biến.
  • Giảng viên và Nhà khoa học: Tài liệu tham khảo chuyên sâu phục vụ giảng dạy các học phần Điều khiển quá trình nâng cao, Điều khiển đa biến và Tối ưu hóa hệ thống điều khiển.
  • Kỹ sư R&D và Tự động hóa công nghiệp: Bộ công thức tường minh tính toán thông số PID giúp rút ngắn thời gian hiệu chỉnh hệ thống tại hiện trường từ nhiều ngày xuống vài giờ, nâng cao năng suất bảo trì và vận hành hệ thống.
  • Doanh nghiệp sản xuất: Tiếp cận giải pháp điều khiển tiên tiến với chi phí đầu tư thấp, tận dụng hạ tầng phần cứng PLC/DCS hiện hữu để nâng cao chất lượng sản phẩm và giảm thiểu phế phẩm.

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ì? Đóng góp độc đáo nhất là việc thiết lập phương pháp giải tích tổng quát để xây dựng ma trận phân ly đơn giản hóa $n \times n$ từ mảng độ lợi tương quan động (DRGA), kết hợp với kỹ thuật xấp xỉ Padé [3/2] cho khâu trễ trong cấu trúc IMC. Điều này cho phép mở rộng lý thuyết IMC của Morari từ hệ đơn biến sang hệ đa biến $n \times n$ một cách tường minh, cung cấp công thức đại số trực tiếp tính toán $K_c, \tau_I, \tau_D$ mà không cần đến các giải thuật lặp hay tối ưu hóa số xấp xỉ.

2. Đổi mới phương pháp luận so với các nghiên cứu trước đây (so sánh cụ thể với BLT và SAT)? So với phương pháp BLT của Luyben (1986) vốn dựa trên quy tắc kinh nghiệm và dò tìm hệ số $F$ để đạt đỉnh cộng hưởng $2n\text{ dB}$, và phương pháp SAT của Shen & Yu (1994) sử dụng rơ-le hồi tiếp thử nghiệm tuần tự dễ gây vọt lố lớn khi trễ cao, phương pháp của luận án là phương pháp giải tích đơn bước (non-iterative analytical method). Nó triệt tiêu hoàn toàn tương tác tĩnh và động thông qua ma trận phân ly $D(s)$, đồng thời sử dụng bộ lọc nối tiếp $F_i(s)$ để kiểm soát trực tiếp hình dạng đáp ứng vòng kín và độ biến thiên tín hiệu điều khiển TV.

3. Phát hiện thực nghiệm đáng ngạc nhiên nhất được chứng minh bằng số liệu là gì? Phát hiện bất ngờ nhất là việc cấu trúc phân ly đơn giản hóa kết hợp xấp xỉ Padé [3/2] và bộ lọc nối tiếp không chỉ vượt trội trong bài toán bám điểm đặt (giảm IAE tới $45%$) mà còn có khả năng kháng nhiễu tải vượt trội so với SAT, đồng thời triệt tiêu hoàn toàn dao động chéo giữa các vòng lặp ($\Delta y_2 < 0.02$ khi $\Delta r_1 = 1.0$). Điều này bác bỏ định kiến kỹ thuật trước đây cho rằng phân ly đơn giản hóa thường hy sinh khả năng kháng nhiễu để đổi lấy tính tinh giản của phần cứng.

4. Luận án có cung cấp đầy đủ giao thức nhân bản (Replication Protocol) không? Có hoàn toàn. Luận án cung cấp chi tiết toàn bộ ma trận hàm truyền của các mô hình chuẩn (WB, VL, OR, WW), bảng công thức giải tích rút gọn tham số cho các mô hình FOPDT/SOPTD, lưu đồ P&ID của tháp chưng cất, sơ đồ khối Simulink chế độ Real-Time Windows Target, mã hiệu phần cứng card DAQ National Instruments và toàn bộ mã nguồn chương trình MATLAB/Simulink trong phần phụ lục.

5. Lộ trình nghiên cứu 10 năm được phác thảo như thế nào? Lộ trình phát triển công nghệ bao gồm:

  • Giai đoạn 1 (1-3 năm): Phát triển thư viện khối hàm (Function Block Library) cho các dòng PLC công nghiệp (Siemens S7-1500, Rockwell ControlLogix) tích hợp thuật toán phân ly tự chỉnh (Self-tuning Decoupler).
  • Giai đoạn 2 (3-5 năm): Mở rộng lý thuyết phân ly thích nghi bền vững cho các đối tượng phi tuyến cao và hệ thống đa biến không vuông ($m \times n$).
  • Giai đoạn 3 (5-10 năm): Tích hợp cấu trúc phân ly giải tích với nền tảng bản sao số (Digital Twin) và tối ưu hóa thời gian thực (RTO) dựa trên trí tuệ nhân tạo cho toàn bộ nhà máy hóa chất và lọc dầu thông minh.

Kết luận

Công trình nghiên cứu của NCS. Lê Linh đã hoàn thành xuất sắc các mục tiêu khoa học và thực tiễn đề ra, tạo nên một dấu ấn học thuật quan trọng trong lĩnh vực kỹ thuật điều khiển tự động với 6 đóng góp cốt lõi:

  1. Tổng quát hóa thành công lý thuyết phân ly đơn giản hóa cho hệ thống đa biến $n \times n$ dựa trên ma trận độ lợi tương quan động (DRGA).
  2. Đề xuất bước đột phá ứng dụng xấp xỉ Padé bậc cao [3/2] cho khâu trễ, bảo toàn chính xác đặc tính động học pha của hệ thống.
  3. Xây dựng hệ thống công thức giải tích tường minh xác định thông số PID-IMC cho các dạng mô hình quá trình chuẩn tắc (FOPDT, SOPTD, NMP).
  4. Tích hợp giải pháp bộ lọc nối tiếp $F_i(s)$, tối ưu hóa đồng thời hai chỉ số chất lượng mâu thuẫn là sai số tích phân tuyệt đối (IAE) và độ êm dịu tín hiệu điều khiển (TV).
  5. Kiểm chứng toàn diện và chứng minh sự vượt trội của phương pháp đề xuất so với các chuẩn quốc tế (Z-N, BLT, SAT) trên các hệ thống kinh điển (Wood & Berry, Vinante & Luyben, Ogunnaike & Ray).
  6. Thiết kế, chế tạo và vận hành thành công 2 hệ thống thực nghiệm thời gian thực (Hệ bồn nước đôi $2 \times 2$ và Tháp chưng cất Ethanol - Nước $3 \times 3$), khẳng định tính khả thi công nghệ và giá trị ứng dụng thực tiễn cao trong sản xuất công nghiệp.