Tổng quan luận án

Điện khí hóa toàn bộ nền kinh tế và đời sống xã hội là nhiệm vụ trọng tâm trong chiến lược công nghiệp hóa tại Việt Nam. Trong bối cảnh nền kinh tế duy trì mức tăng trưởng trung bình cao từ 6% đến 8%/năm, nhu cầu tiêu thụ điện năng tăng nhanh chóng. Cùng với nhiệt điện, điện gió và điện mặt trời, thủy điện giữ vai trò nguồn năng lượng cơ bản. Sau khi các dự án điện hạt nhân tạm ngừng triển khai theo nghị quyết của Quốc hội do những bài học về an toàn từ các sự cố quốc tế tại Chernobyl (Liên Xô cũ) và Fukushima (Nhật Bản năm 2012), vai trò của các nhà máy thủy điện càng trở nên quan trọng.

Lịch sử phát triển ngành thủy điện Việt Nam đã ghi nhận nhiều công trình quy mô lớn: nhà máy thủy điện Thác Bà (1971–1972, công suất $3 \times 36 = 108\text{ MW}$), Đa Nhim (1964), Hòa Bình (1992, $8 \times 240 = 1920\text{ MW}$), Trị An ($4 \times 100 = 400\text{ MW}$), Sơn La (2005–2012, $6 \times 400 = 2400\text{ MW}$), và Lai Châu (2011–2016, $3 \times 400 = 1200\text{ MW}$). Bên cạnh các công trình lớn, các nhà máy thủy điện công suất vừa và nhỏ (như Đa Nhim, Vĩnh Sơn, Thác Mơ, Sông Hinh, Hàm Thuận, Sông Tranh, Rào Quán, Sông Con 2, Đắc Rinh, Ea Krong) phát triển rộng khắp nhằm khai thác tiềm năng sông suối tại các lưu vực địa phương. Nhờ khả năng khởi động và dừng nhanh, các nhà máy thủy điện vừa và nhỏ có giá trị trong việc phủ đỉnh công suất vào các giờ cao điểm.

Tuy nhiên, đặc điểm kỹ thuật cơ bản của các nhà máy thủy điện vừa và nhỏ là thường được xây dựng theo phương án tiết kiệm, không có hồ chứa lớn, không có bể điều tiết hoặc tháp điều áp. Do đó, năng lượng tích trong cột nước (bao gồm thế năng và động năng) thay đổi liên tục theo điều kiện thủy văn và thời tiết. Sự thay đổi này làm biến thiên các tham số động học của tổ hợp turbine – máy phát. Trong khi đó, phần lớn thiết bị điều khiển đi kèm tổ máy nhập khẩu từ nước ngoài sử dụng thuật toán PID kinh điển với các tham số cứng hóa, không có khả năng tự chỉnh thích nghi, dẫn đến việc không khai thác tối ưu nguồn năng lượng cột nước và suy giảm chất lượng ổn định tần số khi phụ tải lưới biến động. Mặt khác, các nghiên cứu và giải pháp công nghệ điều khiển của nhà sản xuất nước ngoài thường mang tính khép kín, chi phí bảo trì thay thế cao.

Mục tiêu nghiên cứu: Mục tiêu của luận án là xây dựng phương pháp và thuật toán tổng hợp điều khiển hệ cánh lái hướng cấp lưu lượng nước vào turbine ở nhà máy thủy điện vừa và nhỏ nhằm thích nghi với sự thay đổi của cột nước và sự thay đổi của tải để đạt hiệu quả kinh tế kỹ thuật.

Nội dung và nhiệm vụ nghiên cứu cụ thể:

  1. Xác định tham số mô hình mô tả các khâu trong hệ thống điều khiển tần số quay máy phát điện phụ thuộc vào năng lượng cột nước.
  2. Đo năng lượng cột nước bằng các phương tiện đo hiện đại (sử dụng cảm biến vi cơ quán tính MEMS và thuật toán lọc Kalman).
  3. Xây dựng thuật toán điều khiển ổn định tần số quay máy phát trong nhà máy thủy điện vừa và nhỏ trên cơ sở áp dụng các phương pháp điều khiển hiện đại.
  4. Mô phỏng và thực nghiệm để kiểm chứng các thuật toán đã đề xuất.

Đối tượng và phạm vi nghiên cứu:

  • Đối tượng nghiên cứu: Các vấn đề xử lý tín hiệu đo lường và điều khiển lưu lượng nước cấp cho turbine nhằm ổn định tần số quay turbine ở các nhà máy thủy điện công suất vừa và nhỏ.
  • Phạm vi nghiên cứu:
    • Thuật toán toán học nhận dạng tham số mô hình động học hệ thống điều khiển tổ máy phát điện dựa trên dữ liệu từ các thiết bị đo lường sẵn có.
    • Thuật toán xác định tham số tối ưu cho bộ điều khiển PID/PI áp dụng trên các thiết bị điều khiển tổ máy đang vận hành thực tế mà không cần can thiệp sâu vào phần cứng.
    • Thuật toán điều khiển tối ưu thích nghi dành cho thiết bị điều khiển tổ máy phục vụ mục tiêu tự chủ thiết kế, chế tạo phần mềm điều khiển trong nước.
    • Kiểm chứng các thuật toán đề xuất thông qua mô hình tính toán và mô phỏng trên máy tính.

Tổng quan tài liệu và vị trí của luận án

Tổng quan các công trình nghiên cứu trong nước và quốc tế cho thấy các hướng tiếp cận chính đối với bài toán điều khiển hệ thống năng lượng và điều tốc turbine thủy điện:

Tác giả / Tài liệu Năm công bố Hướng tiếp cận và nội dung chính Giới hạn / Khoảng trống nghiên cứu
Các luận án tiến sĩ trong nước [14, 15, 16, 17] Giai đoạn trước 2017 Tập trung vào luận chứng kinh tế kỹ thuật, biến đổi dạng năng lượng, ổn định lưới và tích hợp điện gió, điện mặt trời. Chưa đi sâu vào thiết kế thuật toán điều khiển tự thích nghi cho turbine thủy điện vừa và nhỏ.
Nguyễn Đắc Nam [18] (ĐHBK Hà Nội) 2017 Ứng dụng mạng mờ nơ-ron (Fuzzy Neural Network) để chỉnh định thông số bộ điều khiển PID cho hệ điều tốc turbine. Chưa xét đến bài toán biến động năng lượng cột nước đặc thù tại các NMTĐ vừa và nhỏ.
Công trình tại Nga [51] 2009 Xây dựng hệ thống điều khiển công suất nhà máy thủy điện phục vụ tái phân phối năng lượng khi tải lưới thay đổi lớn. Chưa giải quyết bài toán điều khiển tối ưu hoặc điều khiển thích nghi cho turbine.
Công trình tại Nga [50] 2014 Đề xuất giải pháp nhận dạng cấu trúc mô hình cánh lái hướng và tổng hợp cấu trúc hệ điều khiển turbine. Chưa giải quyết nhận dạng tham số mô hình khi độ cao cột nước thay đổi.
Setiawan, Ardyono Priyadi, Mauridhi Hery Purnomo [37] 2011 Sử dụng mạng nơ-ron nhận dạng sự thay đổi của phụ tải để hiệu chỉnh thông số bộ điều khiển thích nghi. Chưa đề cập đến việc tổng hợp lệnh điều khiển tối ưu thích nghi theo trạng thái cột nước.
Fraile-Ardanuy và cộng sự [38] 2006 Sử dụng mạng nơ-ron truyền thẳng đa lớp (ANN 2-5-1) nội suy tối ưu tốc độ từ đường cong đặc tuyến. Giả định mô hình động học hệ thống là cố định và đã biết trước.
Ardyono Priyadi, Mauridhi Hery Purnomo [39] 2008 Ứng dụng cấu trúc mạng nơ-ron tương đương PID với thuật toán Brandt-Lin cho tổ máy công suất nhỏ. Áp dụng cho tổ máy vận hành độc lập với giả thiết tải ít thay đổi.
Dianwet Qian và cộng sự [40] 2012 Ứng dụng điều khiển trượt (SMC) kết hợp giải thuật di truyền (GA) để tối ưu hóa mặt trượt. Chưa giải quyết tính thích nghi với sự biến thiên của năng lượng cột nước.
Ilyas Eker và cộng sự [41] 2003 Thiết kế bộ điều khiển đa tầng, đa biến cho turbine máy phát. Cấu trúc hệ thống đo lường phức tạp, khó triển khai trong thực tế.
Orelind [42] 1989 Thiết kế bộ điều khiển số dựa trên cực tiểu hóa hàm mục tiêu bậc hai tại các điểm vận hành cố định. Các hệ số xác định theo độ mở định trước, chưa tính đến mạch vòng điều khiển cánh lái hướng.
Batlle và cộng sự [43] 2007 Thiết kế bộ điều khiển dựa trên hàm Lyapunov cho mô hình phi tuyến khi không có tháp điều áp. Chưa xét đến các chỉ tiêu tối ưu hóa chất lượng điều khiển.
Dewi Jones, Sa’ad Mansor [44] 2008 Đề xuất bộ điều khiển dự báo dựa trên sai lệch tần số và công suất turbine. Chưa giải quyết tính tối ưu và khả năng thích nghi tham số.
Chen Jian và cộng sự [45] 2014 Thiết kế bộ điều khiển bền vững kết hợp khâu bù khi hằng số thời gian biến thiên trong khoảng xác định. Chưa tính đến sự biến thiên của hệ số khuếch đại mô hình.
Lansberry J.E., Wozniak L. [46] 1994 Sử dụng giải thuật di truyền (GA) tự chỉnh thông số bộ điều tốc theo hằng số thời gian và phụ tải. Chưa xét đến sự thay đổi đồng thời của hệ số khuếch đại đối tượng.

Khoảng trống nghiên cứu được lựa chọn: Qua khảo sát, các công trình quốc tế và trong nước chưa giải quyết trọn vẹn bài toán điều khiển tự thích nghi đồng thời với sự thay đổi của hằng số thời gian $T$, hệ số khuếch đại $K$, áp lực thủy tĩnh $z_2$ và mômen tải $z_1$ phát sinh do biến thiên cột nước tại các nhà máy thủy điện vừa và nhỏ. Luận án lựa chọn giải quyết khoảng trống này bằng cách kết hợp nhận dạng tham số trực tuyến qua phương pháp bình phương tối thiểu, đánh giá biến thiên năng lượng cột nước qua cảm biến vi cơ quán tính (MEMS) và lọc Kalman, từ đó tổng hợp luật điều khiển tối ưu thích nghi.

Cơ sở lý thuyết và phương pháp nghiên cứu

Cơ sở lý thuyết và mô hình toán học

Năng lượng toàn phần chứa trong một đơn vị thể tích chất lỏng của dòng sông được xác định theo phương trình thủy khí động học Bernoulli: $$E = \left( h + \frac{p}{\gamma} + \frac{\alpha_i V^2}{2g} \right) w$$ Trong đó: $h$ là vị năng (chiều cao mức nước), $\frac{p}{\gamma}$ là áp năng ($p$ là áp suất dòng nước, $\gamma = 1000\text{ kg/m}^3$ là trọng lượng riêng của nước), $\frac{\alpha_i V^2}{2g}$ là động năng ($V$ là vận tốc dòng chảy, $g = 9{,}81\text{ m/s}^2$ là gia tốc trọng trường, $\alpha_i$ là hệ số hiệu chỉnh phân bố vận tốc trên mặt cắt), $w$ là thể tích dòng nước ($\text{m}^3$). Năng lượng trao đổi cho turbine phản kích $E_{12} = E_1 - E_2$ phụ thuộc vào độ chênh lệch cột áp trước và sau bánh xe công tác.

Công suất của turbine phản kích: $$N_T = \eta_T \gamma g Q H$$ Trong đó: $\eta_T$ là hiệu suất turbine ($\eta_T \le 1$), $H$ là độ chênh cột nước, $Q$ là lưu lượng dòng chảy ($\text{m}^3/\text{s}$). Lưu lượng $Q$ phụ thuộc tuyến tính vào góc mở cánh lái hướng $\alpha$: $$Q = K_{CLH} \alpha \quad (0 \le \alpha \le \alpha^*)$$ Hệ số $K_{CLH}$ phụ thuộc trực tiếp vào cao trình cột nước và vận tốc dòng chảy. Tần số góc chuẩn danh định của rotor máy phát được xác định bởi: $$\omega_0 = \frac{2\pi f_0}{p}$$ với $f_0 = 50\text{ Hz}$ và $p$ là số đôi cực từ của máy phát điện.

Mô hình động học tổ hợp turbine – máy phát được mô tả qua phương trình cân bằng mômen: $$J_T \frac{d\omega}{dt} + K_c \omega = M_d - M_c$$ Trong đó: $J_T$ là mômen quán tính của rotor tổ máy, $K_c$ là hệ số cản của turbine phụ thuộc áp suất thủy tĩnh ($K_c = \xi h$), $M_d = K_y \alpha$ là mômen động lực do dòng nước tạo ra ($K_y$ tỉ lệ thuận với năng lượng cột nước), $M_c$ là mômen cản phụ tải của lưới điện. Chuẩn hóa phương trình bằng cách chia cho $K_c$, ta thu được: $$T \frac{d\omega}{dt} + \omega = K \alpha - z_1$$ với các tham số phụ thuộc năng lượng cột nước: $$T = \frac{J_T}{K_c}, \quad K = \frac{K_y}{K_c}, \quad z_1 = \frac{M_c}{K_c}$$

Mô hình động học cơ cấu servo quay cánh lái hướng: $$J_c \frac{d^2\alpha}{dt^2} + K_\alpha \frac{d\alpha}{dt} = M_u - M_A$$ Trong đó: $J_c$ là mômen quán tính cụm cánh lái hướng, $K_\alpha$ là hệ số cản quay, $M_u = K_x U$ là mômen điều khiển từ xy lanh thủy lực tỉ lệ với điện áp điều khiển $U$, $M_A$ là mômen thủy tĩnh do áp lực cột nước cản trở chuyển động cánh lái. Chuẩn hóa phương trình ta có: $$T_\alpha \frac{d^2\alpha}{dt^2} + \frac{d\alpha}{dt} = K_u U - z_2$$ với $T_\alpha = \frac{J_c}{K_\alpha}$, $K_u = \frac{K_x}{K_\alpha}$, $z_2 = \frac{M_A}{K_\alpha}$.

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

  1. Phương pháp bình phương tối thiểu (Least Squares): Sử dụng để giải hệ phương trình đại số tuyến tính nhằm nhận dạng các tham số biến thiên $T, K, z_1$ của tổ hợp turbine – máy phát và $z_2$ của hệ cánh lái hướng từ các tập dữ liệu đo đạc rời rạc $(x_j, y_j)$. Hệ phương trình tối thiểu hóa sai số: $$\sum_{j=1}^M \left( y_j - \sum_{i=0}^N a_i x_j^i \right)^2 \to \min \implies \frac{\partial S}{\partial a_k} = 0$$
  2. Lý thuyết lọc Kalman (Kalman Filter): Sử dụng thuật toán lọc Kalman liên tục tuyến tính và lọc Kalman phi tuyến mở rộng (EKF) để xử lý tín hiệu từ các cảm biến vi cơ quán tính MEMS (con quay quán tính, cảm biến gia tốc), nhằm ước lượng trạng thái động học, xác định độ trôi cảm biến và phát hiện thời điểm biến đổi năng lượng cột nước.
  3. Đẳng thức Parseval trong miền tần số: Thiết lập hàm quan hệ giữa hệ số điều khiển và tích phân bình phương sai số bám tần số $\int_0^\infty (\omega(t) - \omega_0)^2 dt$ để tìm bộ tham số tối ưu $K_P, K_I$ cho các bộ điều khiển PID/PI hiện có trên tổ máy.
  4. Phương pháp Backstepping và Lý thuyết điều khiển tối ưu: Tổng hợp luật điều khiển phi tuyến và tối ưu cho cơ cấu van cấp nước để triệt tiêu sai lệch tần số góc $\omega - \omega_0$ dưới tác động của biến thiên phụ tải và áp lực dòng chảy.
  5. Phương pháp mô phỏng số: Xây dựng mô hình trên phần mềm MATLAB / Simulink để mô phỏng kiểm chứng các thuật toán nhận dạng và luật điều khiển thích nghi trong các điều kiện thay đổi tham số cột nước ($K, T$) và phụ tải ($P_e$).

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

Chương 1: Tổng quan vấn đề điều khiển tốc độ turbine trong nhà máy thủy điện vừa và nhỏ

Chương 1 trình bày tổng quan về hệ thống thủy điện Việt Nam, vai trò của các nhà máy thủy điện vừa và nhỏ trong việc phủ đỉnh biểu đồ phụ tải. Chương này phân tích chi tiết các dạng công trình thủy điện: nhà máy thủy điện kiểu đập, nhà máy thủy điện kiểu kênh dẫn (có hoặc không có hầm dẫn nước), và nhà máy thủy điện kiểu hỗn hợp. Tác giả phân tích nguyên lý chuyển hóa năng lượng của dòng nước qua hai loại turbine chính: turbine xung kích (Pelton, Turgo, Banki) và turbine phản kích (Francis, Kaplan cánh quạt, cánh quay, turbine hướng chéo, turbine dòng thẳng).

                                  +-----------------------+
                                  |     Cột nước H, Q     |
                                  +-----------+-----------+
                                              |
                                              v
+------------------+  u (V)   +---------------+-------+  α (độ)   +-------------------+  ω (rad/s)
| Bộ điều khiển    +--------->| Mạch vòng vị trí      +---------->| Turbine & Máy     +----------->
| tốc độ quay      |          | hệ cánh lái hướng     |           | phát điện         |
+--------^---------+          +-----------------------+           +---------+---------+
         |                                                                  |
         +-------------------------- Đo tốc độ quay ω ----------------------+

Chương 1 chỉ rõ sơ đồ chức năng và sơ đồ cấu trúc hệ thống điều khiển ổn định tốc độ quay tổ hợp turbine – máy phát điện. Tác giả phân tích hiện trạng công nghệ điều khiển tại Việt Nam: phần cứng servomotor và hệ thống cơ khí thủy lực đã làm chủ chế tạo trong nước với giá thành thấp hơn khoảng 40% so với nhập khẩu, nhưng phần mềm và thuật toán điều tốc vẫn phụ thuộc vào hệ thống PID cứng hóa của nước ngoài. Qua việc phân tích các công trình nghiên cứu trong nước và quốc tế, chương 1 xác định rõ mục tiêu, đối tượng, phạm vi nghiên cứu và các nhiệm vụ kỹ thuật cần giải quyết của luận án.

Chương 2: Xây dựng thuật toán nhận dạng tham số mô hình động học mô tả tổ hợp turbine - máy phát của nhà máy thủy điện vừa và nhỏ

Chương 2 tập trung thiết lập các thuật toán nhận dạng tham số mô hình dựa trên hai công cụ toán học: phương pháp bình phương tối thiểu và bộ lọc Kalman.

  • Mô hình hóa động học: Thiết lập phương trình vi phân mô tả mạch vòng vị trí cánh lái hướng và động học rotor turbine – máy phát, chỉ ra 4 tham số phụ thuộc vào năng lượng cột nước và phụ tải gồm: hằng số thời gian $T$, hệ số khuếch đại $K$, nhiễu phụ tải $z_1$, và áp lực thủy tĩnh $z_2$.
  • Thuật toán nhận dạng tham số bằng phương pháp bình phương tối thiểu: Sử dụng dữ liệu đo đạc thực tế gồm điện áp điều khiển $U$, góc mở cánh lái hướng $\alpha$, và tần số quay rotor máy phát $\omega$ để nhận dạng tham số mô hình turbine và hệ thống cánh lái hướng. Thuật toán rời rạc hóa các phương trình vi phân và giải hệ phương trình đại số tuyến tính để tìm nghiệm tối ưu cho $T, K, z_1, z_2$.
  • Đánh giá năng lượng cột nước bằng bộ lọc Kalman và cảm biến MEMS: Ứng dụng cảm biến vi cơ quán tính (MEMS) để đo tốc độ dòng chảy và áp lực, kết hợp thuật toán lọc Kalman ước lượng trạng thái thực, loại bỏ sai số do độ trôi con quay $\varepsilon_c$ và nhiễu trắng $\xi_m$. Kết quả ước lượng năng lượng cột nước là căn cứ xác định chính xác thời điểm hệ thống cần kích hoạt thuật toán nhận dạng lại tham số và hiệu chỉnh luật điều khiển cánh lái hướng.

Chương 3: Tổng hợp luật điều khiển ổn định tần số quay turbine trong nhà máy thủy điện vừa và nhỏ

Chương 3 áp dụng lý thuyết điều khiển hiện đại để giải quyết bài toán ổn định tần số máy phát điện khi công suất phát thay đổi theo yêu cầu của lưới:

  • Xác định tham số tối ưu cho bộ điều khiển PID/PI bằng công thức Parseval: Luận án thiết lập hàm tương quan giữa tham số tỉ lệ $K_P$, tích phân $K_I$ trong luật điều khiển PID với sai số bám tần số $\omega(t) - \omega_0$. Áp dụng tích phân Parseval trong miền phức Laplace $s$, tác giả giải bài toán cực tiểu hóa chỉ tiêu sai số: $$J = \int_0^\infty (\omega(t) - \omega_0)^2 dt = \frac{1}{2\pi j} \int_{-j\infty}^{j\infty} E(s) E(-s) ds \to \min$$ Giải pháp này cho phép chỉnh định tối ưu bộ thông số $K_P, K_I$ cho mạch vòng cánh lái hướng và mạch vòng turbine trên các thiết bị điều khiển PID sẵn có mà không cần thay đổi phần cứng.
  • Xây dựng thuật toán Backstepping và luật điều khiển tối ưu thích nghi: Thiết kế thuật toán Backstepping ổn định tần số quay turbine bằng cách phân rã hệ thống thành các khâu bậc thấp, từng bước xây dựng hàm Lyapunov để đảm bảo tính ổn định toàn cục.
  • Quan sát trạng thái và thuật toán điều khiển van cấp nước: Ứng dụng lọc Kalman xây dựng bộ quan sát trạng thái đầy đủ phục vụ thiết lập luật điều khiển tối ưu thích nghi cho thiết bị điều khiển turbine tự thiết kế.

Chương 4: Mô phỏng kiểm nghiệm

Chương 4 xây dựng mô hình mô phỏng trên nền tảng MATLAB / Simulink để kiểm chứng tính đúng đắn và hiệu quả của các thuật toán nhận dạng và điều khiển đã đề xuất.

+-----------------------------------------------------------------------------------+
|                            CẤU TRÚC HỆ THỐNG MÔ PHỎNG                             |
|                                                                                   |
|  +------------------------------------+    +-----------------------------------+  |
|  |     KHỐI NHẬN DẠNG THAM SỐ         |    |     KHỐI ĐIỀU KHIỂN TỐI ƯU        |  |
|  | Nhận dạng trực tuyến T(t) và K(t)  +--->| Hiệu chỉnh tín hiệu u(t) thích    |  |
|  | theo phương pháp bình phương       |    | ứng với biến thiên cột nước/tải   |  |
|  | tối thiểu                          |    +-----------------+-----------------+  |
|  +-----------------+------------------+                      |                    |
|                    ^                                         v                    |
|                    |               +-------------------------+---+                |
|                    +---------------+     ĐỐI TƯỢNG TURBINE -     |                |
|                                    |          MÁY PHÁT           |                |
|                                    +-----------------------------+                |
+-----------------------------------------------------------------------------------+

Các kịch bản mô phỏng được thực hiện bao gồm:

  1. Mô phỏng nhận dạng tham số: Kiểm chứng đồ thị nhận dạng hằng số thời gian $T$ và hệ số khuếch đại $K$ của mô hình động học khi đối tượng chịu tác động của nhiễu đo lường.
  2. Mô phỏng luật điều khiển tối ưu và so sánh: Đánh giá đáp ứng tần số quay turbine $\omega(t)$, góc quay cánh lái hướng $\alpha(t)$ và giá trị hàm chỉ tiêu chất lượng $J$. Kết quả cho thấy khi thực hiện đúng luật điều khiển tối ưu, hàm chỉ tiêu $J$ đạt giá trị nhỏ nhất; trong khi các trường hợp hệ số điều khiển $K_x$ tăng hoặc giảm sai lệch so với giá trị tối ưu, hàm $J$ đều tăng và chất lượng ổn định tần số suy giảm.
  3. Mô phỏng tính thích nghi khi tham số cột nước biến thiên: Khảo sát các trường hợp:
    • Hệ số $K$ và hằng số thời gian $T$ cùng tăng.
    • Hệ số $K$ và hằng số thời gian $T$ cùng giảm.
    • Hệ số $K$ tăng và $T$ giảm. Đồ thị mô phỏng ghi nhận thời điểm biến thiên tham số, sự điều chỉnh góc mở cánh lái hướng $\alpha$ và quá trình phục hồi tần số góc $\omega$ về giá trị danh định $\omega_0$.
  4. Mô phỏng đáp ứng khi tải tiêu thụ thay đổi: Khảo sát đáp ứng góc quay cánh lái hướng $\alpha$ và tần số quay $\omega$ tại các thời điểm tải yêu cầu $P_e$ giảm và tải yêu cầu $P_e$ tăng, chứng minh khả năng duy trì ổn định tần số của hệ thống.

Kết quả và những đóng góp mới

Luận án đã đạt được các kết quả khoa học và đóng góp mới sau:

  1. Về mặt lý luận và học thuật:

    • Xây dựng thành công mô hình toán học tích hợp mô tả động học của tổ hợp turbine – máy phát và cơ cấu quay cánh lái hướng, trong đó tường minh hóa các tham số phụ thuộc vào năng lượng cột nước ($T, K, z_1, z_2$).
    • Thiết lập thuật toán nhận dạng tham số trực tuyến cho tổ máy thủy điện dựa trên phương pháp bình phương tối thiểu xử lý các tín hiệu đo điện áp điều khiển $U$, góc mở cánh lái $\alpha$ và tốc độ quay $\omega$.
    • Đề xuất giải pháp ứng dụng cảm biến vi cơ quán tính MEMS kết hợp thuật toán lọc Kalman để đánh giá biến thiên năng lượng cột nước, xác lập tiêu chí khoa học xác định thời điểm cần nhận dạng lại tham số và hiệu chỉnh luật điều khiển.
    • Ứng dụng công thức Parseval thiết lập phương pháp giải tích xác định bộ tham số tối ưu cho bộ điều khiển PID/PI điều chỉnh lưu lượng nước cấp vào turbine.
    • Xây dựng thuật toán Backstepping và thuật toán điều khiển tối ưu thích nghi trạng thái cho cơ cấu van cấp nước turbine.

  2. Về mặt thực tiễn và công nghệ:

    • Cung cấp cơ sở khoa học để hiệu chỉnh, tối ưu hóa các hệ thống điều tốc PID đang vận hành tại các nhà máy thủy điện vừa và nhỏ mà không cần thay thế phần cứng.
    • Cung cấp thuật toán và cấu trúc phần mềm phục vụ việc tự chủ thiết kế, chế tạo hoàn chỉnh thiết bị điều khiển tổ hợp turbine – máy phát trong nước, giúp giảm giá thành đầu tư (tiết kiệm khoảng 40% chi phí thiết bị điều khiển servomotor) và chủ động trong bảo trì, vận hành.

  3. Danh mục các công trình khoa học đã công bố của tác giả:

    • Công trình về thuật toán nhận dạng và lọc Kalman: Công bố trên Tạp chí Nghiên cứu Khoa học và Công nghệ Quân sự, số 47 (trang 70–75) và số 50.
    • Công trình về tổng hợp luật điều khiển tối ưu và ứng dụng công thức Parseval: Công bố trên Tạp chí Nghiên cứu Khoa học và Công nghệ Quân sự, số 53 (trang 38–48); Tạp chí Khoa học, Tập 46, số 3A.
    • Công trình về mô phỏng kiểm nghiệm hệ thống điều khiển thích nghi: Công bố trên Tạp chí Khoa học và Công nghệ Năng lượng (Trường Đại học Điện lực), số 14.

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

  • Hạn chế: Các kết quả kiểm chứng thuật toán nhận dạng và luật điều khiển tối ưu thích nghi trong luận án chủ yếu được thực hiện thông qua mô phỏng số trên phần mềm MATLAB / Simulink và mô hình toán học; chưa thực hiện thử nghiệm đầy đủ trên hệ thống tổ máy phát điện công nghiệp thực tế tại nhà máy thủy điện đang vận hành do điều kiện an toàn lưới điện.
  • Hướng nghiên cứu tiếp:
    • Triển khai thử nghiệm bán tự nhiên (Hardware-in-the-Loop - HIL) trên các hệ thống điều khiển thời gian thực (DSP/PLC công nghiệp).
    • Tích hợp hoàn thiện thuật toán vào phần cứng điều khiển nhúng chuyên dụng và tiến hành thử nghiệm thực tế tại một nhà máy thủy điện công suất vừa và nhỏ cụ thể.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên sâu cho các đối tượng sau:

  • Nghiên cứu sinh, học viên cao học và giảng viên ngành Kỹ thuật điều khiển và Tự động hóa, Kỹ thuật điện, Năng lượng tái tạo: Tham khảo phương pháp mô hình hóa hệ thống cơ điện – thủy lực phi tuyến, ứng dụng lý thuyết lọc Kalman, phương pháp bình phương tối thiểu và kỹ thuật Backstepping trong điều khiển thích nghi.
  • Kỹ sư vận hành và chuyên gia kỹ thuật tại các nhà máy thủy điện vừa và nhỏ: Ứng dụng quy trình tính toán tối ưu tham số PID theo công thức Parseval (Chương 3) để nâng cao chất lượng ổn định tần số tổ máy hiện hữu.
  • Các viện nghiên cứu và doanh nghiệp chế tạo thiết bị điện: Tham khảo cấu trúc thuật toán điều khiển và chương trình mô phỏng MATLAB/Simulink (Phần Phụ lục 1–4) phục vụ thiết kế, chế tạo thiết bị điều tốc turbine nội địa.

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

1. Tại sao sự biến thiên của năng lượng cột nước lại ảnh hưởng trực tiếp đến chất lượng điều khiển tần số ở nhà máy thủy điện vừa và nhỏ?

Ở các nhà máy thủy điện vừa và nhỏ không có hồ chứa lớn hoặc tháp điều áp, độ cao cột nước $h$ và lưu lượng dòng chảy thay đổi thường xuyên. Theo các phương trình động lực học, áp lực thủy tĩnh và động năng dòng nước quyết định trực tiếp đến hệ số cản turbine $K_c$ và hệ số mômen động lực $K_y$. Khi đó, hằng số thời gian $T = J_T/K_c$ và hệ số khuếch đại đối tượng $K = K_y/K_c$ bị biến thiên liên tục. Nếu bộ điều khiển sử dụng các thông số cố định, hệ thống sẽ không duy trì được chất lượng ổn định tần số danh định ($f_0 = 50\text{ Hz}$) khi phụ tải thay đổi.

2. Phương pháp bình phương tối thiểu được sử dụng trong luận án nhằm nhận dạng những thông số nào?

Phương pháp bình phương tối thiểu được sử dụng để nhận dạng trực tuyến 4 tham số biến thiên của mô hình động học:

  • Hằng số thời gian $T$ của tổ hợp turbine – máy phát.
  • Hệ số khuếch đại $K$ phản ánh mối quan hệ giữa góc mở cánh lái và mômen sinh ra.
  • Nhiễu mômen phụ tải quy đổi $z_1 = M_c/K_c$.
  • Áp lực mômen thủy tĩnh cản trở chuyển động quay cánh lái hướng $z_2 = M_A/K_\alpha$. Dữ liệu đầu vào cho quá trình nhận dạng là các đại lượng đo lường thực tế gồm điện áp điều khiển $U(t)$, góc mở cánh lái $\alpha(t)$, và tốc độ quay rotor $\omega(t)$.

3. Vai trò của thiết bị vi cơ quán tính (MEMS) và bộ lọc Kalman trong hệ thống điều khiển là gì?

Cảm biến MEMS (con quay quán tính và cảm biến gia tốc) được bố trí để thu thập thông tin về chuyển động và áp lực thủy động của dòng nước. Bộ lọc Kalman (tuyến tính liên tục và phi tuyến mở rộng) được áp dụng để khử nhiễu trắng $\xi_m$ và độ trôi cảm biến $\varepsilon_c$, từ đó đánh giá chính xác sự biến thiên năng lượng cột nước. Đây là tín hiệu chỉ thị để hệ thống quyết định kích hoạt thuật toán nhận dạng tham số và cập nhật luật điều khiển tối ưu thích nghi.

4. Việc ứng dụng công thức Parseval trong Chương 3 giải quyết vấn đề kỹ thuật nào?

Nhiều nhà máy thủy điện đang sử dụng bộ điều khiển PID nhập khẩu với các hệ số $K_P, K_I, K_D$ cố định chưa tối ưu. Việc can thiệp vào cấu trúc phần cứng của các hệ thống này rất phức tạp. Luận án ứng dụng công thức Parseval để tính tích phân bình phương sai số bám tần số trong miền phức, từ đó tìm ra nghiệm giải tích cho bộ hệ số tỉ lệ $K_P$ và tích phân $K_I$ tối ưu. Giải pháp này cho phép nâng cao hiệu suất khai thác năng lượng của các thiết bị hiện có mà không cần thay đổi phần cứng.

5. Các kịch bản mô phỏng trong Chương 4 đã kiểm chứng những nội dung gì?

Chương 4 thực hiện mô phỏng trên MATLAB/Simulink qua các kịch bản:

  • Kiểm chứng tốc độ hội tụ và độ chính xác của thuật toán nhận dạng tham số $T$ và $K$.
  • Đánh giá chỉ tiêu chất lượng $J$, góc mở cánh lái $\alpha$ và tốc độ $\omega$ giữa luật điều khiển tối ưu so với các trường hợp điều khiển không tối ưu ($K_x$ tăng hoặc giảm).
  • Khảo sát đáp ứng thích nghi của hệ thống khi tham số cột nước biến động ($K$ và $T$ cùng tăng; $K$ và $T$ cùng giảm; $K$ tăng và $T$ giảm).
  • Khảo sát đáp ứng điều khiển ổn định tần số khi phụ tải yêu cầu $P_e$ đột ngột tăng hoặc giảm.

Kết luận

Luận án tiến sĩ kỹ thuật của tác giả Đặng Tiến Trung đã giải quyết bài toán đo lường và điều khiển ổn định tần số quay tổ hợp turbine – máy phát trong nhà máy thủy điện công suất vừa và nhỏ trong điều kiện năng lượng cột nước và tải thay đổi. Bằng việc kết hợp phương pháp bình phương tối thiểu, kỹ thuật lọc Kalman với cảm biến MEMS, công thức Parseval và lý thuyết điều khiển tối ưu Backstepping, tác giả đã xây dựng hệ thống thuật toán nhận dạng và điều khiển thích nghi có cơ sở khoa học rõ ràng. Các kết quả mô phỏng trên MATLAB/Simulink đã khẳng định tính đúng đắn của thuật toán, đóng góp cơ sở lý luận và giải pháp thực tiễn quan trọng cho việc nâng cao chất lượng khai thác các tổ máy hiện có cũng như tự chủ thiết kế thiết bị điều khiển thủy điện tại Việt Nam.