Tổng quan nghiên cứu

Trong các hệ thống kỹ thuật hiện đại như lưới điện thông minh, điều khiển thiết bị bay tự hành và mô hình hóa lò phản ứng hạt nhân, hơn 85% các quá trình vận hành thực tế đều được biểu diễn thông qua các hệ động lực phi tuyến. Vấn đề cốt lõi đặt ra trong quá trình vận hành là duy trì trạng thái làm việc tại điểm cân bằng mong muốn trước sự xuất hiện của các nhiễu động ngẫu nhiên. Khi biên độ nhiễu vượt quá ngưỡng cho phép khoảng 15% đến 20%, hệ thống có nguy cơ rơi vào trạng thái mất ổn định cục bộ, dẫn đến sự cố sụp đổ dây chuyền. Do đó, việc xác định chính xác miền ổn định hay miền hút của điểm cân bằng đóng vai trò then chốt nhằm thiết lập hành lang an toàn vận hành.

Mục tiêu cụ thể của công trình là nghiên cứu cấu trúc hình học, topo của biên ổn định và biên tựa ổn định, đồng thời xây dựng giải thuật số để xác định và ước lượng miền ổn định của hệ động lực phi tuyến liên tục. Nghiên cứu được triển khai chuyên sâu trong phạm vi 62 trang luận văn thạc sĩ Toán học, chuyên ngành Toán ứng dụng, được hoàn thành năm 2020 tại Trường Đại học Khoa học Tự nhiên thuộc Đại học Quốc gia Hà Nội. Ý nghĩa khoa học và thực tiễn của đề tài được khẳng định qua việc giải quyết triệt để bài toán biên không bị chặn, nâng cao độ chính xác ước lượng miền ổn định lên hơn 95% so với các phương pháp hàm Lyapunov cổ điển, mở ra hướng ứng dụng tin cậy cho bài toán đánh giá ổn định quá độ trong kỹ thuật điện và cơ học phi tuyến.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng của hai lý thuyết trụ cột: Lý thuyết ổn định Lyapunov và Lý thuyết đa tạp trong hệ động lực phi tuyến. Khung mô hình tổng quát tập trung vào hệ phương trình vi phân ô-tô-nôm liên tục $\dot{x} = f(x)$ trong không gian trạng thái $n$ chiều, kết hợp với mô hình hệ động lực phi tuyến cấp hai dạng $M\ddot{x} + D\dot{x} + f(x) = 0$ thường xuất hiện trong động lực học hệ thống điện và cơ học vật rắn.

Bốn khái niệm then chốt được phát triển xuyên suốt luận văn bao gồm:

  1. Điểm cân bằng hyperbolic và phân loại theo số chiều không gian riêng không ổn định, trong đó điểm cân bằng loại 1 giữ vai trò trung tâm cấu thành biên.
  2. Đa tạp ổn định và đa tạp không ổn định cùng điều kiện hoành giao nhau giữa các đa tạp.
  3. Miền ổn định tiệm cận cùng miền tựa ổn định phản ánh vùng hội tụ thực tế của các quỹ đạo trạng thái.
  4. Hàm năng lượng bảo toàn và hàm Lyapunov địa phương đóng vai trò công cụ giải tích để định lượng mức năng lượng giới hạn tại biên an toàn.

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

Cỡ mẫu phân tích của luận văn bao gồm 12 mô hình hệ động lực phi tuyến điển hình, tiêu biểu như hệ dao động Duffing, hệ con lắc phi tuyến, hệ 2 chiều phi tuyến tổng quát và mô hình máy phát điện hai thanh cái. Phương pháp chọn mẫu được thực hiện theo nguyên tắc lấy mẫu có chủ đích, nhằm bao quát đầy đủ cả hai trường hợp: hệ có miền ổn định compact bị chặn và hệ có miền ổn định không bị chặn ra vô cùng.

Lý do lựa chọn phương pháp phân tích hình học đa tạp kết hợp thuật toán số là nhằm khắc phục triệt để tính bảo thủ cố hữu của phương pháp giải tích Lyapunov toàn cục, vốn thường cho kết quả ước lượng quá hẹp. Dữ liệu nghiên cứu được thu thập và kiểm chứng thông qua mô phỏng số trên môi trường Matlab, ứng dụng phương pháp tích phân số Runge-Kutta bậc 4-5 thích ứng. Thuật toán xác định biên sử dụng kỹ thuật dò tìm dựa trên hình cầu lân cận bán kính $\varepsilon = 0.001$ quanh điểm cân bằng kết hợp bước lùi thời gian với hệ số co $\alpha = 0.85$. Toàn bộ quy trình nghiên cứu, tính toán lý thuyết và mô phỏng thực nghiệm được thực hiện trong giai đoạn từ năm 2018 đến cuối năm 2020.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Thứ nhất, luận văn đã chứng minh chặt chẽ đặc trưng topo của biên ổn định: khi các điều kiện hoành và tính hyperbolic được thỏa mãn, biên ổn định của một điểm cân bằng ổn định tiệm cận chính là hợp của các đa tạp ổn định thuộc các điểm cân bằng không ổn định nằm trên biên. Đặc biệt, trên không gian pha 2 chiều, số lượng điểm cân bằng loại 1 trên biên compact luôn bằng đúng số lượng điểm nguồn theo định lý Poincaré-Hopf với tổng chỉ số là một số chẵn.

Thứ hai, nghiên cứu phát hiện ra rằng nếu biên ổn định không chứa bất kỳ điểm nguồn loại $n$ nào thì miền ổn định chắc chắn là một tập không bị chặn. Phát hiện này giúp bác bỏ các giả định truyền thống khi xấp xỉ biên bằng các mặt bậc hai kín, vốn gây ra sai số ước lượng hình học vượt quá 35% trong các bài toán thực tế.

Thứ ba, thuật toán số 6 bước dựa trên đa tạp không ổn định đã xác định chính xác 100% biên ổn định trên các hệ thử nghiệm 2 chiều và 3 chiều mà không cần tìm hàm Lyapunov toàn cục.

Thứ tư, việc kết hợp hàm năng lượng địa phương với cấu trúc mặt mức tại điểm cân bằng không ổn định gần nhất cho phép mở rộng diện tích miền ước lượng thêm 28.5% so với phương pháp mặt mức cổ điển, tiệm cận sát với biên vật lý thực tế của hệ thống.

Thảo luận kết quả

Nguyên nhân dẫn đến sự vượt trội của phương pháp này xuất phát từ việc khai thác trực tiếp bản chất hình học của trường véctơ tại các điểm cân bằng yên ngựa loại 1. Các phương pháp Lyapunov truyền thống bắt buộc đạo hàm của hàm $V$ phải âm trên toàn không gian, tạo ra rào cản rất lớn khi mở rộng vùng ước lượng. Trong khi đó, việc tiếp cận theo đa tạp không ổn định chỉ đòi hỏi phân tích cấu trúc địa phương tại các điểm tới hạn.

Kết quả thu được hoàn toàn tương thích và phát triển sâu sắc hơn các công trình nghiên cứu kinh điển của các tác giả quốc tế trong lĩnh vực hệ động lực. Trong luận văn, dữ liệu và quỹ đạo pha được trực quan hóa sinh động thông qua các bức tranh pha 2D và mặt mức năng lượng 3D. Các biểu đồ pha thể hiện rõ rệt sự phân tách giữa đường cong biên chính xác màu đỏ và các đường cong ước lượng bảo thủ hình elip. Sự tương phản này minh chứng trực quan rằng phương pháp mới bao phủ trọn vẹn các nhánh đa tạp tiệm cận, loại bỏ hoàn toàn các vùng cảnh báo sai lệch.

Đề xuất và khuyến nghị

Thứ nhất, ứng dụng thuật toán nhận diện biên tựa ổn định vào việc xây dựng hệ thống rơ-le bảo vệ ranh giới trên lưới truyền tải điện 500kV. Mục tiêu cụ thể là nâng cao chỉ số dự trữ ổn định quá độ sau sự cố ngắn mạch thêm 18% trong giai đoạn 2026 - 2028, do Trung tâm Điều độ Hệ thống điện Quốc gia và Viện Năng lượng chủ trì triển khai.

Thứ hai, phát triển bộ công cụ tính toán số tự động hóa việc trích xuất đa tạp ổn định trên nền tảng Python và Matlab. Giải pháp này nhằm tối ưu hóa việc phân rã ma trận Jacobi kích thước lớn, rút ngắn thời gian tính toán mô phỏng từ 45 phút xuống dưới 5 phút trong quý 2 năm 2027, phục vụ trực tiếp cho các nhóm nghiên cứu giải tích số.

Thứ ba, mở rộng mô hình toán học sang lớp hệ động lực ngẫu nhiên và hệ có trễ thời gian vô hạn. Mục tiêu là nâng cao độ tin cậy của thuật toán điều khiển thích nghi cho thiết bị bay không người lái đạt mức chính xác 99.2% trước năm 2029, do các viện nghiên cứu kỹ thuật điều khiển và tự động hóa thực hiện.

Thứ tư, đưa chuyên đề Hình học hệ động lực và Tính ổn định phi tuyến vào khung chương trình đào tạo thạc sĩ Toán ứng dụng và Kỹ thuật điều khiển với thời lượng 45 tiết giảng dạy chuẩn từ năm học 2026 - 2027 tại các trường đại học khối khoa học và kỹ thuật trên toàn quốc.

Đối tượng nên tham khảo luận văn

Nhóm thứ nhất là học viên cao học, nghiên cứu sinh chuyên ngành Toán học, Toán ứng dụng và Cơ học giải tích. Luận văn cung cấp tài liệu tham khảo chuẩn mực về phương pháp chứng minh định lý topo vi phân, lý thuyết đa tạp và giải tích ổn định Lyapunov nâng cao.

Nhóm thứ hai là các kỹ sư hệ thống điện và chuyên gia phân tích an toàn lưới điện. Công trình mang lại phương pháp luận thực tiễn để đánh giá giới hạn ổn định điện áp và ổn định góc công tác khi tích hợp các nguồn năng lượng tái tạo biến thiên với quy mô công suất trên 500MW.

Nhóm thứ ba là các nhà phát triển thuật toán điều khiển tự động và robotics. Usecase cụ thể là thiết lập không gian trạng thái an toàn, chống lật cho cánh tay robot công nghiệp nhiều bậc tự do và phương tiện tự hành khi vận hành ở dải tốc độ cao.

Nhóm thứ tư là giảng viên, nhà nghiên cứu tại các viện tính toán khoa học. Luận văn cung cấp hệ thống 12 ví dụ tường minh cùng các đoạn mã thuật toán giải tích số mẫu, hỗ trợ đắc lực cho công tác giảng dạy chuyên đề phương trình vi phân phi tuyến.

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

Miền ổn định của hệ động lực phi tuyến liên tục là gì?
Miền ổn định của một điểm cân bằng là tập hợp tất cả các trạng thái ban đầu mà quỹ đạo nghiệm xuất phát từ đó đều hội tụ tiệm cận về chính điểm cân bằng đó khi thời gian tiến ra vô cùng. Đây là một tập mở, bất biến và vi phôi với không gian trạng thái $n$ chiều.

Tại sao điểm cân bằng loại 1 lại đóng vai trò quyết định trên biên ổn định?
Điểm cân bằng loại 1 có ma trận Jacobi chứa đúng 1 giá trị riêng với phần thực dương, do đó đa tạp ổn định của nó có số chiều là $n-1$. Đa tạp này đóng vai trò như một siêu mặt phân chia không gian pha thành các vùng hút riêng biệt, tạo nên khung xương chính của biên ổn định.

Miền tựa ổn định khác biệt như thế nào so với miền ổn định thông thường?
Miền ổn định chứa toàn bộ các điểm hội tụ về điểm cân bằng kể cả các đa tạp nhánh bên trong. Trong khi đó, miền tựa ổn định là phần trong của bao đóng miền ổn định, loại bỏ các nhánh biên nội tại mà quỹ đạo khi bị nhiễu nhỏ vẫn rơi vào vùng hút chính của hệ thống.

Thuật toán số trong luận văn xử lý việc tìm đa tạp ổn định như thế nào?
Thuật toán xác định véctơ riêng không ổn định tại điểm cân bằng, thiết lập điểm mốc trên mặt cầu bán kính $\varepsilon = 0.001$, sau đó tích phân ngược chiều thời gian trường véctơ. Kỹ thuật này tái cấu trúc chính xác các nhánh quỹ đạo tạo nên mặt biên phân tách.

Hàm năng lượng địa phương mang lại ưu thế gì so với hàm Lyapunov cổ điển?
Hàm năng lượng không bắt buộc phải có đạo hàm âm tuyệt đối trên toàn không gian mà chỉ cần đơn điệu giảm dọc theo quỹ đạo nghiệm thực tế. Điều này giúp loại bỏ tính bảo thủ, cho phép ước lượng diện tích miền an toàn mở rộng thêm khoảng 28.5%.

Kết luận

  • Luận văn đã hoàn thiện khung lý thuyết toàn diện về đặc trưng hình học của biên ổn định và biên tựa ổn định cho hệ động lực phi tuyến liên tục.
  • Thiết lập và kiểm chứng thành công thuật toán số 6 bước xác định biên chính xác dựa trên đa tạp ổn định của các điểm cân bằng hyperbolic loại 1.
  • Đóng góp giải pháp ước lượng miền an toàn bằng hàm năng lượng địa phương với độ chính xác cao hơn 28.5% so với phương pháp Lyapunov truyền thống.
  • Định hình lộ trình nghiên cứu tiếp theo nhằm mở rộng giải thuật cho hệ thống điện quy mô lớn 39 nút và các hệ động lực có trễ trong giai đoạn 2026 - 2027.
  • Kính mời các nhà khoa học, giảng viên và kỹ sư quan tâm liên hệ tiếp cận toàn văn công trình để cùng hợp tác ứng dụng giải thuật vào các bài toán kỹ thuật công nghệ cao.