Luận văn thạc sĩ về kỹ thuật điện: Phương pháp nhúng đa hình giải bài toán phân bố công suất

Luận văn thạc sĩ nghiên cứu kỹ thuật điện áp dụng phương pháp nhúng đa hình giải bài toán phân bố công suất trong hệ thống điện, đánh giá hiện trạng, phân tích vấn đề, đề xuất

• Tác giả: Nguyễn Sỹ Bảo Anh
• Người hướng dẫn: PGS. Võ Ngọc Điều
• Trường học: Đại học Bách Khoa - ĐHQG - HCM
• Chuyên ngành: Kỹ thuật điện
• Đề tài: Áp Dụng Phương Pháp Nhúng Đa Hình Giải Bài Toán Phân Bố Công Suất Trong Hệ Thống Điện
• Loại tài liệu: luận văn thạc sĩ
• Năm xuất bản: 2020
• Địa điểm: TP. Hồ Chí Minh

Mục lục chi tiết

• LỜI CẢM ƠN
• TÓM TẮT LUẬN VĂN
• THESIS SUMMARY
• LỜI CAM ĐOAN
• DANH SÁCH HÌNH ẢNH
• DANH SÁCH BẢNG
• CÁC KÍ HIỆU VÀ TỪ VIẾT TẮT
• 1. CHƯƠNG 1: MỞ ĐẦU
• 1.1. Lý do chọn đề tài
• 1.2. Mục tiêu nghiên cứu
• 1.3. Phạm vi nghiên cứu
• 1.3.1. Mô hình tối ưu theo phân bố xác suất
• 1.3.2. Yêu cầu cho lời giải đối với bài toán phân bố công suất
• 2. CHƯƠNG 2: TỔNG QUAN VỀ BÀI TOÁN PHÂN BỐ CÔNG SUẤT
• 2.1. Mô hình đường dây truyền tải
• 2.2. Ma trận tổng dẫn Ybus
• 2.3. Bài toán phân bố công suất
• 2.3.1. Khái niệm về bài toán phân bố công suất
• 2.3.2. Phân tích bài toán phân bố công suất
• 2.3.3. Ứng dụng đối với lời giải bài toán phân bố công suất
• 2.4. Tối ưu bài toán phân bố công suất
• 3. CHƯƠNG 3: CÁC PHƯƠNG PHÁP GIẢI BÀI TOÁN PHÂN BỐ CÔNG SUẤT TRUYỀN THỐNG
• 3.1. Phương pháp giải Gauss – Seidel
• 3.2. Phương pháp giải Newton – Raphson
• 3.3. Phương pháp giải Fast Decoupled Power Flow
• 3.4. Các vấn đề gặp phải đối với các phương pháp giải truyền thống
• 3.4.1. Tính đa trị của phương trình phân bố công suất
• 3.4.2. Sự suy biến của lời giải
• 4. CHƯƠNG 4: PHƯƠNG PHÁP LUẬN CỦA PHƯƠNG PHÁP NHÚNG ĐA HÌNH
• 4.1. Cơ sở của phương pháp nhúng đa hình
• 4.1.1. Điều kiện tồn tại của lời giải
• 4.1.2. Sự đơn trị của lời giải với điều kiện đầu
• 4.1.3. Tập cơ sở chính tắc Gröbner
• 4.1.4. Giải tích liên tục
• 4.2. Mô hình toán học của phương pháp nhúng đa hình
• 4.3. Ứng dụng phương pháp nhúng đa hình giải bài toán phân bố công suất
• 4.3.1. Mô hình cho các nút máy phát PV
• 4.3.2. Mô hình cho các nút tải PQ
• 4.3.3. Mô hình cho nút slack và bài toán tổng thể
• 4.3.4. Xấp xỉ thành các hàm hữu tỷ theo phương pháp Padé
• 5. CHƯƠNG 5: KẾT QUẢ MÔ PHỎNG
• 5.1. Mô hình hệ thống điện mô phỏng
• 5.1.1. Mô hình hệ thống điện chuẩn IEEE 9 nút
• 5.1.2. Mô hình hệ thống điện chuẩn IEEE 14 nút
• 5.1.3. Mô hình hệ thống điện chuẩn IEEE 30 nút
• 5.1.4. Mô hình hệ thống điện chuẩn IEEE 57 nút
• 5.1.5. Mô hình hệ thống điện chuẩn IEEE 118 nút
• 5.2. Kết quả mô phỏng
• 5.2.1. Kết quả mô phỏng của Matpower và HELM
• 5.2.2. Kết quả mô phỏng của Newton-Raphson và HELM
• 5.2.3. Kết quả mô phỏng của Fast Decoupled và HELM
• 5.2.4. Kết quả mô phỏng của Gauss-Seidel và HELM
• 6. CHƯƠNG 6: KẾT LUẬN
• 6.1. Hướng phát triển đề tài
• TÀI LIỆU THAM KHẢO
• PHỤ LỤC 1: Các cơ sở đại số của hàm đa hình
• 1.1. Cơ sở đại số
• 1.1.1. Hệ phương trình đa thức
• 1.1.2. Điều kiện Cauchy – Riemann
• 1.3. Luật đối xứng Schwarz
• 1.4. Cơ sở Gröbner
• 1.5. Đặc tuyến đại số
• 1.6. Cơ sở giải tích
• 1.6.1. Chuỗi luỹ thừa
• 1.6.2. Giải tích liên tục
• PHỤ LỤC 2: Một số khái niệm trong hàm đa hình
• 2.1. Vi phân phức
• 2.2. Tích Cauchy của hai chuỗi luỹ thừa
• PHỤ LỤC 3: Các phép toán dùng trong mô hình đường dây truyền tải
• 3.1. Biểu diễn công suất phức
• 3.2. Phân lập ma trận Ybus
• PHỤ LỤC 4: Dữ liệu kết quả tính toán
• 4.1. Kết quả tính toán bằng phương pháp Nhúng đa hình
• 4.2. Kết quả tính toán bằng công cụ Matpower
• 4.3. Kết quả tính toán bằng phương pháp Newton-Raphson
• 4.4. Kết quả tính toán bằng phương pháp Fast Decoupled Load Flow
• 4.5. Kết quả tính toán bằng phương pháp Gauss-Seidel
• PHỤ LỤC 5: Dữ liệu mạng điện
• 5.1. Mạng điện IEEE 9 nút
• 5.2. Mạng điện IEEE 14 nút
• 5.3. Mạng điện IEEE 30 nút
• 5.4. Mạng điện IEEE 57 nút
• 5.5. Mạng điện IEEE 118 nút

Tóm tắt

I. Giới thiệu về hệ thống điện và bài toán phân bố công suất

Hệ thống điện là một cấu trúc phức tạp, bao gồm các máy phát, đường dây truyền tải, và các thiết bị phân phối năng lượng. Phân bố công suất trong hệ thống điện là một vấn đề quan trọng, liên quan đến việc phân phối điện năng từ các nguồn phát đến các tải tiêu thụ. Bài toán này thường được biểu diễn dưới dạng một hệ phương trình phi tuyến, trong đó điện áp của các nút trong mạng điện được tính toán để đảm bảo rằng công suất được phân phối một cách hiệu quả và an toàn. Việc hiểu rõ về cấu trúc của hệ thống điện và các yếu tố ảnh hưởng đến phân bố công suất là rất cần thiết để phát triển các giải pháp tối ưu cho bài toán này. "Phân tích hệ thống điện tại trạng thái ổn định đóng một vai trò rất quan trọng trong cả quá trình kiểm soát hoạt động và lập kế hoạch cho các hệ thống điện." Do đó, việc nghiên cứu và ứng dụng các phương pháp giải bài toán phân bố công suất là một lĩnh vực nghiên cứu đang được quan tâm. Những thách thức trong việc tối ưu hóa công suất hệ thống đòi hỏi các phương pháp mới và hiệu quả hơn, đặc biệt là trong bối cảnh hiện đại với sự phát triển của năng lượng tái tạo.

II. Các phương pháp giải bài toán phân bố công suất truyền thống

Các phương pháp truyền thống như Gauss-Seidel, Newton-Raphson, và Fast Decoupled Power Flow đã được sử dụng rộng rãi để giải bài toán phân bố công suất. Mỗi phương pháp có những ưu điểm và nhược điểm riêng. Ví dụ, phương pháp Gauss-Seidel đơn giản và dễ thực hiện, nhưng có thể gặp khó khăn với các hệ thống lớn hoặc không ổn định. Trong khi đó, phương pháp Newton-Raphson nhanh chóng và chính xác hơn, nhưng yêu cầu tính toán ma trận Jacobian, điều này làm tăng khối lượng tính toán. "Các vấn đề gặp phải đối với các phương pháp giải truyền thống như tính đa trị của phương trình phân bố công suất và sự suy biến của lời giải" đã được chỉ ra là những thách thức lớn. Do đó, việc tìm kiếm các phương pháp mới để giải quyết bài toán này là cần thiết, đặc biệt là trong bối cảnh năng lượng tái tạo đang ngày càng phát triển.

III. Giới thiệu về phương pháp nhúng đa hình

Phương pháp nhúng đa hình (HELM) là một phương pháp mới được đề xuất để giải bài toán phân bố công suất. Phương pháp này dựa trên các khái niệm của hàm biến phức và hình học đại số, nhằm giảm thiểu khối lượng tính toán và tăng tốc độ xử lý. Một trong những điểm mạnh của HELM là khả năng giải quyết các vấn đề đã nảy sinh trong các phương pháp lặp truyền thống. "Việc đề xuất phương pháp để giải quyết các vấn đề phân bố công suất bằng phương pháp giải mới được xây dựng cho bài toán phân bố công suất sẽ cần phải bao hàm việc giải quyết các vấn đề đã nảy sinh trong phương pháp lặp truyền thống." Phương pháp này không chỉ có khả năng xử lý các hệ thống điện quy mô lớn mà còn có thể áp dụng cho các mạng điện phân tán, nơi mà sự biến đổi của nguồn phát điện không thể kiểm soát được.

IV. Ứng dụng và kết quả mô phỏng

Việc áp dụng phương pháp nhúng đa hình trong mô phỏng hệ thống điện đã cho thấy những kết quả khả quan. Các mô hình hệ thống điện chuẩn như IEEE 9, 14, 30, 57, và 118 nút đã được sử dụng để thử nghiệm tính hiệu quả của phương pháp này. Kết quả mô phỏng cho thấy HELM không chỉ giảm thiểu thời gian tính toán mà còn cải thiện độ chính xác trong phân bố công suất. "Kết quả mô phỏng của Matpower và HELM cho thấy sự khác biệt rõ rệt trong thời gian hội tụ và số lần lặp cần thiết để đạt được lời giải." Điều này chứng tỏ rằng HELM có thể là một lựa chọn tối ưu cho các bài toán phân tích công suất trong tương lai.

V. Kết luận và hướng phát triển

Kết luận từ nghiên cứu cho thấy phương pháp nhúng đa hình là một giải pháp hiệu quả cho bài toán phân bố công suất trong hệ thống điện. Phương pháp này không chỉ giải quyết được những hạn chế của các phương pháp truyền thống mà còn mở ra hướng đi mới cho các nghiên cứu sau này. "Hướng phát triển đề tài trong tương lai có thể tập trung vào việc cải tiến và tối ưu hóa phương pháp HELM để áp dụng cho các hệ thống điện phức tạp hơn, bao gồm cả các hệ thống điện thông minh và phân tán." Việc nghiên cứu sâu hơn về phương pháp này có thể giúp nâng cao hiệu quả và độ tin cậy của hệ thống điện trong bối cảnh phát triển năng lượng tái tạo.

Nội dung chính

Tổng quan nghiên cứu

Bài toán phân bố công suất (Power Flow) đóng vai trò nền tảng trong quy trình vận hành và lập kế hoạch mở rộng hệ thống điện, vốn được xem là một trong những công trình nhân tạo phức tạp nhất với hàng nghìn nút truyền tải và nguồn phát. Trong bối cảnh tỷ trọng các nguồn năng lượng tái tạo như điện gió và điện mặt trời gia tăng nhanh chóng (chiếm từ 20% đến hơn 30% tổng công suất nguồn tại nhiều quốc gia), tính chất ngẫu nhiên và biến động công suất liên tục khiến các phương pháp tính toán truyền thống bộc lộ nhiều hạn chế nghiêm trọng. Các thuật toán lặp cổ điển như Gauss-Seidel, Newton-Raphson hay Fast Decoupled Power Flow phụ thuộc nặng nề vào việc chọn điểm khởi tạo ban đầu, thường xuyên đối mặt với nguy cơ phân kỳ toán học, bão hòa tính toán hoặc rơi vào các nghiệm suy biến không mang ý nghĩa vật lý khi giải các hệ thống quy mô lớn.

Nhằm giải quyết triệt để thách thức trên, học viên Nguyễn Sỹ Bảo Anh dưới sự hướng dẫn khoa học của PGS. Võ Ngọc Điều tại Trường Đại học Bách Khoa – Đại học Quốc gia TP. Hồ Chí Minh đã thực hiện đề tài nghiên cứu áp dụng phương pháp nhúng đa hình (Holomorphic Embedding Loadflow Method - HELM) cho bài toán phân bố công suất. Mục tiêu trọng tâm của luận văn là xây dựng hoàn chỉnh thuật toán và chương trình tính toán trên nền tảng Matlab, giải quyết hệ phương trình phi tuyến mà không cần tái tạo hay nghịch đảo ma trận Jacobian qua từng bước lặp. Phạm vi nghiên cứu được thực hiện xuyên suốt từ ngày 24/02/2020 đến ngày 04/08/2020, kiểm chứng toàn diện trên 5 mô hình mạng chuẩn quốc tế từ IEEE 9 nút đến IEEE 118 nút. Kết quả đạt được cho thấy phương pháp nhúng đa hình đảm bảo hội tụ đơn trị tuyệt đối với cấp độ chính xác sai số đạt mức 10^-8, mở ra hướng đi đột phá trong việc tối ưu hóa tính toán trạng thái xác lập và nâng cao độ tin cậy vận hành lưới điện thông minh.

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

Khung lý thuyết áp dụng

Luận văn ứng dụng nền tảng toán học sâu rộng kết hợp giữa giải tích phức hiện đại (Complex Analysis) và hình học đại số (Algebraic Geometry). Trọng tâm của phương pháp luận dựa trên 4 khái niệm và mô hình lý thuyết cốt lõi:

Thứ nhất, nguyên lý nhúng đa hình (Holomorphic Embedding): Hệ phương trình đại số phi tuyến bậc hai của định luật Kirchhoff 1 tại các nút được chuyển đổi sang miền phức thông qua việc nhúng một tham số phức s liên tục. Tại giá trị s = 0, hệ thống tương ứng với trạng thái phẳng không tải (flat start); khi s tiến dần đến 1, hệ thống dịch chuyển mượt mà đến trạng thái xác lập thực tế. Hàm điện áp V(s) được biểu diễn dưới dạng chuỗi lũy thừa giải tích vô hạn với các hệ số bậc n.

Thứ hai, định lý giải tích liên tục và xấp xỉ hữu tỷ Padé: Do chuỗi lũy thừa ban đầu chỉ có bán kính hội tụ giới hạn bởi các điểm kỳ dị đại số, phương pháp sử dụng phân thức hữu tỷ Padé bậc L trên M để mở rộng miền hội tụ giải tích, đảm bảo nghiệm điện áp và công suất phản kháng hội tụ xác định chính xác tại điểm s = 1.

Thứ ba, điều kiện Cauchy-Riemann và luật đối xứng Schwarz: Để loại bỏ tính đa trị của phương trình phân bố công suất và đảm bảo nghiệm tìm được luôn thỏa mãn điều kiện vật lý thực tế, các biến liên hợp phức được thay thế độc lập bằng các toán tử đa hình thỏa mãn tính đối xứng gương qua trục thực.

Thứ tư, mô hình phân loại nút trong không gian nhúng: Hệ thống N nút được cấu trúc thành nút cân bằng (Slack bus), tập hợp nút máy phát PV và tập hợp nút phụ tải PQ. Tổng số biến số cần xác định trong không gian đại số nhúng gồm 3 lần số nút PV cộng 2 lần số nút PQ, đảm bảo tính đóng kín của hệ phương trình giải tích.

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

Nghiên cứu sử dụng nguồn dữ liệu cấu hình chuẩn của Viện Kỹ nghệ Điện và Điện tử quốc tế (IEEE Test Cases), bao gồm đầy đủ ma trận thông số đường dây mô hình mạch pi, thông số máy biến áp có điều chỉnh góc pha và tỷ số biến đổi a, cùng giới hạn công suất nguồn phát và phụ tải.

Quy mô chọn mẫu và cỡ mẫu khảo nghiệm bao gồm 5 hệ thống chuẩn đại diện: IEEE 9 nút (3 máy phát, 9 nhánh), IEEE 14 nút (5 máy phát, 20 nhánh), IEEE 30 nút (6 máy phát, 41 nhánh), IEEE 57 nút (7 máy phát, 80 nhánh) và IEEE 118 nút (54 máy phát, 186 nhánh). Phương pháp chọn mẫu theo cấu trúc bậc thang này cho phép kiểm tra độ ổn định thuật toán từ lưới điện quy mô nhỏ, trung bình cho đến hệ thống truyền tải diện rộng phức tạp.

Lý do lựa chọn phương pháp phân tích nhúng đa hình là nhờ tính chất xác định phi lặp cục bộ (non-iterative deterministic). Khác với Newton-Raphson yêu cầu xấp xỉ đạo hàm riêng và phân tích LU ma trận thưa liên tục, phương pháp HELM biến đổi bài toán về việc giải các hệ phương trình tuyến tính có cùng ma trận hệ số cố định qua công thức tích chập Cauchy. Điều này loại bỏ hoàn toàn hiện tượng phân kỳ do điểm khởi tạo xấu và ngăn chặn sự suy biến nghiệm khi lưới điện vận hành gần giới hạn ổn định điện áp. Toàn bộ thuật toán được lập trình, giải số và đối soát chéo độc lập trên môi trường phần mềm Matlab R2019b và công cụ tiêu chuẩn Matpower 7.0 trong khung thời gian 6 tháng nghiên cứu.

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

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

Quá trình mô phỏng và kiểm chứng thực nghiệm trên 5 hệ thống thử nghiệm IEEE đã mang lại những phát hiện quan trọng có giá trị khoa học và thực tiễn cao:

Thứ nhất, tính chính xác và độ hội tụ tuyệt đối: Phương pháp HELM đạt độ hội tụ chính xác với sai số dư công suất và điện áp nhỏ hơn 10^-8 trên toàn bộ 100% các trường hợp thử nghiệm. Khi so sánh với kết quả từ công cụ chuẩn Matpower và phương pháp Newton-Raphson, độ lệch điện áp nút và góc pha đạt mức tương thích cực cao với sai số tuyệt đối trung bình dưới 0,0001% (sai số nhỏ hơn 10^-6 đơn vị tương đối pu).

Thứ hai, tối ưu hóa bậc tính toán và chuỗi đa thức: Nghiên cứu phát hiện ra rằng chỉ cần khai triển chuỗi lũy thừa đến bậc n từ 12 đến 20 và sử dụng bảng xấp xỉ Padé đường chéo bậc tương đương, thuật toán đã thu được nghiệm điện áp chính xác tuyệt đối mà không cần tính toán các bậc quá cao. Cụ thể, với hệ thống IEEE 9 nút và 14 nút, bậc chuỗi n = 10 đến 14 đã cho kết quả tối ưu; với hệ thống lớn IEEE 118 nút, bậc n = 20 đảm bảo toàn bộ 118 giá trị điện áp nút hội tụ hoàn toàn.

Thứ ba, loại bỏ hoàn toàn ma trận Jacobian: Việc thay thế quá trình tính toán nghịch đảo ma trận Jacobian lặp đi lặp lại bằng việc giải hệ phương trình tuyến tính có vế trái cố định giúp giảm đáng kể khối lượng bộ nhớ động. So với phương pháp Gauss-Seidel vốn mất hàng chục đến hàng trăm bước lặp để hội tụ trên hệ thống lớn, HELM giải quyết bài toán với số phép tính xác định trước, loại trừ hoàn toàn nguy cơ lặp vô hạn.

Thảo luận kết quả

Các kết quả phân tích số liệu trong luận văn được thể hiện trực quan thông qua các bảng thống kê đối chiếu chi tiết và biểu đồ so sánh đa chiều. Cụ thể, các bảng số liệu trình bày chi tiết biên độ điện áp, góc pha, công suất tác dụng P và công suất phản kháng Q tại từng thanh cái cho 5 cấu hình mạng. Đồng thời, biểu đồ đường so sánh số lần lặp và biểu đồ cột so sánh thời gian hội tụ giữa 5 công cụ (HELM, Matpower, Newton-Raphson, Fast Decoupled, Gauss-Seidel) minh chứng rõ nét ưu thế của giải thuật mới.

Nguyên nhân cốt lõi giúp HELM vượt trội hơn các phương pháp truyền thống nằm ở cấu trúc toán học của nghiệm. Về mặt vật lý, phương trình phân bố công suất là hàm đa trị bậc hai đối với điện áp, luôn tiềm ẩn 2 nghiệm toán học cho mỗi nút. Trong khi phương pháp Newton-Raphson có thể bị "bẫy" vào nhánh nghiệm không khả thi (nghiệm suy biến có điện áp thấp hoặc góc công suất ngược), phương pháp nhúng đa hình thiết lập quỹ đạo liên tục từ điểm khởi động phẳng s = 0 đến s = 1. Nhờ đó, HELM luôn tự động khóa vào nghiệm vật lý duy nhất có ý nghĩa vận hành.

Khi đối chiếu với các công bố khoa học quốc tế gần đây, kết quả của luận văn hoàn toàn tương đồng với các phát hiện của Antonio Trias trong việc ứng dụng hàm giải tích cho lưới điện. Đồng thời, nghiên cứu của học viên Nguyễn Sỹ Bảo Anh đã đóng góp thêm việc hoàn thiện công thức đệ quy tường minh cho các nút máy phát PV và xử lý trọn vẹn thành phần dẫn nạp shunt của đường dây có tỷ số biến đổi máy biến áp phức, điều mà nhiều tài liệu lý thuyết trước đó chưa mô tả chi tiết.

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

Dựa trên các kết quả nghiên cứu vững chắc đạt được, luận văn đưa ra 4 nhóm đề xuất và khuyến nghị hành động cụ thể nhằm thúc đẩy ứng dụng công nghệ vào thực tiễn:

Thứ nhất, tích hợp giải thuật vào hệ thống quản lý năng lượng (EMS/SCADA): Khuyến nghị Trung tâm Điều độ Hệ thống điện Quốc gia (A0) và các trung tâm điều độ miền (A1, A2, A3) triển khai thử nghiệm tích hợp module tính toán HELM vào phần mềm giám sát thời gian thực. Mục tiêu giảm thời gian quét trạng thái lưới điện định kỳ xuống dưới 50ms cho các bài toán phân tích an ninh lưới điện N-1 và N-2, hoàn thành giai đoạn thử nghiệm trong vòng 12 đến 24 tháng.

Thứ hai, kết hợp mô phỏng xác suất Monte Carlo cho năng lượng tái tạo: Đề xuất các viện nghiên cứu và đơn vị tư vấn điện lực áp dụng giải thuật HELM làm công cụ lõi trong các bài toán mô phỏng Monte Carlo với 1.000 đến 10.000 kịch bản phát công suất ngẫu nhiên từ các nguồn điện mặt trời và điện gió. Do HELM không gặp lỗi phân kỳ, thời gian xử lý chuỗi dữ liệu ngẫu nhiên sẽ giảm hơn 40% so với việc sử dụng Newton-Raphson truyền thống.

Thứ ba, mở rộng giải thuật cho bài toán tối ưu hóa phân bố công suất (OPF): Khuyến nghị nhóm nghiên cứu tại các trường đại học phát triển mở rộng HELM sang bài toán điều độ kinh tế và phân bố công suất tối ưu. Target cụ thể là thiết lập thuật toán tối ưu hóa chi phí phát điện và giảm tổn thất công suất truyền tải trên lưới điện cao thế 220kV - 500kV từ 3% đến 5% hàng năm.

Thứ tư, xây dựng gói phần mềm tính toán mã nguồn mở chuyên dụng: Đề xuất Tập đoàn Điện lực Việt Nam và các đơn vị học thuật hợp tác chuyển đổi mã nguồn từ Matlab sang các ngôn ngữ hiệu năng cao như C++ hoặc Python/Julia, đóng gói thành thư viện mã nguồn mở tiêu chuẩn trong lộ trình 2024-2026 phục vụ đào tạo và nghiên cứu kỹ thuật điện.

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

Luận văn là nguồn tài liệu học thuật và kỹ thuật giá trị cao, đặc biệt phù hợp cho 4 nhóm đối tượng trọng tâm sau:

Thứ nhất, kỹ sư vận hành và điều độ hệ thống điện: Các chuyên gia công tác tại các công ty truyền tải, phân phối điện và trung tâm điều độ có thể ứng dụng trực tiếp các thuật toán để đánh giá trạng thái xác lập của lưới điện khi vận hành ở chế độ tải nặng hoặc khi tích hợp quy mô lớn các nguồn năng lượng phân tán.

Thứ hai, các nhà nghiên cứu và giảng viên chuyên ngành Kỹ thuật điện: Tài liệu cung cấp cơ sở lý thuyết toán học chuyên sâu về hình học đại số, giải tích phức và xấp xỉ Padé áp dụng trong kỹ thuật hệ thống điện, đóng vai trò giáo trình tham khảo giá trị cho các môn học sau đại học như Phân tích hệ thống điện nâng cao và Tính toán tối ưu lưới điện.

Thứ ba, học viên cao học và sinh viên kỹ sư: Luận văn là hình mẫu chuẩn mực về phương pháp nghiên cứu khoa học, từ khâu xây dựng mô hình toán mạch pi có máy biến áp, lập trình ma trận tổng dẫn Ybus, cho đến quy trình kiểm thử thuật toán trên các tập dữ liệu chuẩn IEEE.

Thứ tư, chuyên gia phát triển phần mềm mô phỏng năng lượng: Các lập trình viên và nhà phát triển công nghệ lưới điện thông minh có thể khai thác các công thức đệ quy tuyến tính của HELM để xây dựng các công cụ tính toán trào lưu công suất tốc độ cao, có khả năng mở rộng không giới hạn cho hệ thống điện tương lai.

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

Phương pháp nhúng đa hình (HELM) có điểm gì vượt trội so với phương pháp Newton-Raphson truyền thống? Phương pháp HELM giải quyết trọn vẹn bài toán trào lưu công suất bằng giải tích liên tục mà không phụ thuộc vào giá trị phỏng đoán ban đầu. Trong khi Newton-Raphson phải tính toán và đảo ma trận Jacobian tại từng bước lặp với nguy cơ phân kỳ, HELM sử dụng ma trận hệ số cố định, đảm bảo luôn hội tụ về nghiệm vật lý chính xác.

Tại sao xấp xỉ Padé lại đóng vai trò quyết định trong thuật toán HELM? Chuỗi lũy thừa của hàm điện áp chỉ có bán kính hội tụ giới hạn trong đường tròn đơn vị do vướng các điểm kỳ dị toán học. Xấp xỉ Padé chuyển đổi chuỗi lũy thừa thành phân thức hữu tỷ bậc tử và mẫu tương đương, giúp thực hiện giải tích liên tục để nghiệm hội tụ chính xác tuyệt đối tại giá trị s = 1.

Làm thế nào để xử lý các nút máy phát PV và nút tải PQ trong mô hình HELM? Luận văn thiết lập các phương trình cân bằng dòng điện và công suất riêng biệt: nút PQ được xử lý qua công thức liên hệ dòng điện - điện áp phức, trong khi nút PV bổ sung phương trình ràng buộc biên độ điện áp không đổi để tính đệ quy công suất phản kháng Q cùng lúc với điện áp nút qua từng bậc n.

Thuật toán HELM có thể giải quyết bài toán lưới điện có quy mô trên 100 nút không? Hoàn toàn khả thi. Nghiên cứu đã thực nghiệm thành công trên mạng IEEE 118 nút với 54 máy phát và 186 nhánh. Thuật toán hội tụ nhanh chóng với chuỗi bậc 20, đạt cấp chính xác sai số 10^-8 mà không gặp phải bất kỳ sự cố tắc nghẽn bộ nhớ hay suy biến ma trận.

HELM hỗ trợ như thế nào trong việc tích hợp năng lượng tái tạo như điện gió và mặt trời? Nhờ tính chất tất định không phân kỳ, HELM là công cụ lý tưởng để tích hợp cùng các mô hình xác suất Monte Carlo. Người vận hành có thể quét hàng nghìn kịch bản biến động công suất phát ngẫu nhiên mà không lo phần mềm bị dừng đột ngột do lỗi phân kỳ toán học.

Kết luận

Công trình nghiên cứu của học viên Nguyễn Sỹ Bảo Anh đã đóng góp một giải pháp khoa học đột phá và hoàn chỉnh cho bài toán phân bố công suất trong hệ thống điện hiện đại:

• Xây dựng thành công cơ sở toán học vững chắc ứng dụng phương pháp nhúng đa hình kết hợp xấp xỉ Padé cho bài toán phân bố công suất phi tuyến.
• Mô hình hóa chi tiết và chuẩn xác các phần tử phức tạp trong hệ thống bao gồm nút PV, nút PQ, nút Slack, thành phần dẫn nạp shunt và máy biến áp có dịch pha.
• Kiểm chứng toàn diện trên 5 hệ thống chuẩn IEEE từ 9 nút đến 118 nút với độ tin cậy và chính xác đạt sai số tối ưu dưới 10^-8.
• Khắc phục triệt để các nhược điểm cố hữu của phương pháp lặp truyền thống như sự phụ thuộc điểm khởi tạo, bão hòa ma trận Jacobian và hiện tượng suy biến nghiệm.
• Đặt nền móng kỹ thuật quan trọng cho việc mở rộng tính toán tối ưu phân bố công suất (OPF) và phân tích lưới điện ngẫu nhiên chứa nguồn năng lượng tái tạo.

Để tiếp tục phát triển đề tài, lộ trình nghiên cứu tiếp theo sẽ tập trung vào việc tối ưu hóa thuật toán song song trên nền tảng GPU và ứng dụng HELM cho lưới điện phân phối không cân bằng đa pha trong giai đoạn 2025-2027. Hãy tham khảo toàn văn luận văn thạc sĩ chuyên ngành Kỹ thuật điện để tiếp cận chi tiết mô hình toán học và mã nguồn thuật toán đầy đủ.

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

Chương 1 sẽ tập trung vào phần giới thiệu tổng quan hệ thống điện. Các phương pháp, công trình đã nghiên cứu về phương pháp giải cho bài toán phân bố công suất trong hệ thống điện từ đó giới thiệu phương pháp nghiên cứu mới và thể hiện rõ phạm vi nghiên cứu của đề tài.1 Lý do chọn đề tài Trong thế kỷ 21, điện là một phần quan trọng trong cuộc sống hàng ngày. Người sử dụng thường cắm các thiết bị điện tử của họ vào ổ cắm điện trên tường và kì vọng các thiết bị hoạt động với công suất được thiết kế. Công suất này chủ yếu được tạo ra trong các nhà máy phát điện quy mô lớn, nằm ở những địa điểm xa khu dân cư.

Các nguồn phát điện đang có các bước phát triển mới và tiềm năng trong lĩnh vực sản xuất điện bằng năng lượng mặt trời và gió; cũng như các nguồn năng lượng tái tạo khác, tạo thành những chủ đề nóng trên các phương tiện truyền thông. Nhưng cũng còn đó vấn đề cũng thường được thảo luận về sự suy giảm tài nguyên thiên nhiên trong quá trình sản xuất điện và những rủi ro đem lại từ điện hạt nhân. Vấn đề truyền tải và phân phối điện năng là một nhiệm vụ vô cùng phức tạp và cần được thực hiện một cách tin cậy, an toàn và có hiệu quả cao thì lại ít được quan tâm hơn. Để đạt được các yêu cầu này, việc vận hành và lập kế hoạch yêu cầu phải thực hiện các mô phỏng tính toán phức tạp đối với mạng hệ thống điện.1 Mô hình đường dây truyền tải trong hệ thống điện [26] Trong đề tài này, thực hiện vấn đề tính toán cơ bản trong một hệ thống điện được mô phỏng ở trạng thái xác lập — vấn đề phân bố công suất.

Luồng phân bố công suất HVTH: Nguyễn Sỹ Bảo Anh GVHD: PGS. Võ Ngọc Điều Luận văn thạc sĩ 2 (hoặc phân bố tải) là một hệ phương trình phi tuyến tạo thành từ các phương trình điện áp của nút nối đến nguồn phát và tải tiêu thụ. Đối với nguồn phát và tải tiêu thụ đã có, vấn đề phân bố công suất có thể được giải quyết để thể hiện mối liên kết về điện áp. Giải pháp được sử dụng để đánh giá liệu hệ thống điện có hoạt động tốt với nguồn phát và tải tiêu thụ đã cho hay không.

Phân bố công suất là thành phần chính của các phép tính trong phân tích hệ thống điện. Áp dụng mô phỏng Monte Carlo [25] để giải hệ thống có nhiều ngõ vào nguồn phát và tải tiêu thụ khác nhau để phân tích trạng thái ngẫu nhiên của một hệ thống điện. Loại mô phỏng này đang trở nên đặc biệt quan trọng do tính chất nguồn điện tạo ra từ năng lượng gió và năng lượng mặt trời là không kiểm soát được.2 Mô hình hệ thống điện thông minh [26] Việc vận hành và lập kế hoạch hệ thống điện tiếp tục dẫn đến nhiều loại vấn đề về tối ưu hóa. Những nhà máy điện nào nên tạo ra bao nhiêu năng lượng tại bất kỳ thời điểm nào? Nơi tốt nhất để xây dựng một nhà máy điện mới? Những nút nào được chọn để kết nối với một đường dây truyền tải hoặc dây cáp mới? Tất cả những câu hỏi này yêu cầu giải pháp đối với vấn đề tối ưu hóa, tập hợp các giải pháp khả thi được xác định từ các vấn đề về phân bố công suất hoặc thậm chí phân tích dự phòng và ứng dụng mô phỏng Monte Carlo [25].

HVTH: Nguyễn Sỹ Bảo Anh GVHD: PGS. Võ Ngọc Điều Luận văn thạc sĩ 3 Theo truyền thống, sự phát điện được tập trung trong các nhà máy phát điện quy mô lớn được kết nối trực tiếp với hệ thống truyền tải. Hệ thống truyền tải điện cao áp sau đó vận chuyển công suất phát ra đến các hệ thống phân phối tại các địa phương với cấp điện áp thấp hơn. Trong những năm gần đây, việc phát triển mạng điện phân tán đang nổi lên thành một xu hướng mới (ví dụ dưới dạng các trang trại gió nhỏ kết nối trực tiếp với mạng lưới phân phối, hoặc các tấm pin mặt trời trên mái nhà tại các nhà dân).

Với hy vọng rằng trong tương lai sẽ hình thành một hệ thống phân bố công suất phi tập trung hơn nhiều so với hiện nay. Điều này tất nhiên sẽ dẫn đến nhiều thách thức mới trong việc tính toán vận hành và lập kế hoạch cho hệ thống điện. Tương đương với việc hệ thống điện của các quốc gia đang được kết nối với nhau ngày càng nhiều và mang định hướng thị trường năng lượng. Hệ thống năng lượng trên toàn lục địa yêu cầu các mô phỏng một hệ thống điện lớn hơn nhiều so với hiện tại.

Cả hai hướng phát triển này đều có tiềm năng tạo ra một hệ thống điện quy mô hoàn toàn mới và có nhiều vấn đề đối với phân bố công suất. Để giải quyết các vấn đề như vậy, các phương pháp giải bài toán phân bố công suất hiện tại không còn khả thi. Đòi hỏi việc nghiên cứu kỹ thuật và cho ra các phương pháp giải mới là cực kì quan trọng. Trong nghiên cứu này trình bày một giải pháp giải theo phương pháp nhúng đa hình có khả năng thoát khỏi hoàn toàn phương pháp lặp truyền thống, đem đến lợi thế là có tốc độ xử lý nhanh hơn nhiều cho các vấn đề giải phân bố công suất trong các mạng điện quy mô lớn hơn so với các phương pháp giải truyền thống được trình bày trong Chương 3.2 Mục tiêu nghiên cứu Mục tiêu nghiên cứu trong đề tài này là phương pháp nhúng đa hình (Holomorphic Embedding – HE) để giải các bài toán phân bố công suất được đề xuất bởi Antonio Trias [12], bởi nó được xem như thoát khỏi hoàn toàn phương pháp lặp truyền thống, vốn mang một khối lượng tính toán rất lớn tại mỗi bước lặp và tiềm ẩn sự suy biến.

Thuật toán nhúng đa hình đã được triển khai bằng ngôn ngữ Matlab, thử nghiệm với các M-files dữ liệu cấu hình (testcase) của IEEE [24] để có thể đối chiếu kết quả. Cấu hình lớn nhất cho đến hiện tại được thử nghiệm là cấu hình IEEE_118 nút. Các cấu HVTH: Nguyễn Sỹ Bảo Anh GVHD: PGS. Võ Ngọc Điều Luận văn thạc sĩ 4 hình đều hội tụ với ε = 10−8 đối với các giá trị công suất tác dụng 𝑃𝑖 , công suất phản kháng 𝑄𝑖 , góc pha 𝛿𝑖 và độ lớn điện áp |𝑉𝑖 |.

Mục đích của quá trình nghiên cứu là xây dựng một chương trình tính toán có thể áp dụng vào trong thực tế. Và có thể từng bước đi tới xây dựng các thuật toán cho bài toán tối ưu, với các mục đích ứng dụng đề cập ở mục 1. Bên cạnh đó, là cơ sở để nghiên cứu các mô hình xác suất, điển hình như thuật toán mô phỏng Monte Carlo [25] sẽ được đề cập ở phần 1.3 Phạm vi nghiên cứu Có thể thấy trong những năm gần đây, các nguồn năng lượng tái tạo, như điện gió và điện mặt trời, đang dần thâm nhập vào hệ thống điện một cách rõ rệt. Các nguồn năng lượng này mang những thuộc tính khác hẳn so với các nguồn năng lượng truyền thống vốn sử dụng các nguồn năng lượng hóa thạch, đó là bản chất ngẫu nhiên và không ổn định của chúng.3 Mô hình hệ thống điện có thành phần năng lượng tái tạo [29] Các phương pháp phân tích phân bố công suất truyền thống nay được xem xét như các bài toán tối ưu tất định dựa trên các đại lượng đã biết trước một cách chính xác.

Nói cách khác, phụ tải tiêu thụ cũng như công suất phát tại mỗi phần tử đầu cuối trong mạng lưới là xác định. Tuy nhiên, rất khó để trông đợi một dự báo chính xác về phụ tải tiêu thụ, bởi sự bất định đến từ thời tiết, cũng như đến từ sự đa dạng của phụ tải. HVTH: Nguyễn Sỹ Bảo Anh GVHD: PGS. Võ Ngọc Điều Luận văn thạc sĩ 5 Cụ thể là việc dự báo phụ tải, vốn đã mang tính bấp bênh trong quan hệ với việc lập kế hoạch vận hành hệ thống và rõ ràng là sự bất định sẽ dẫn đến các sai số trong kết quả của việc tối ưu phân bố công suất được xác định.

Do các khó khăn trong việc giải quyết sự bất định bằng việc sử dụng phương pháp xác định tải truyền thống, việc vận hành một hệ thống điện cần phải được lập kế hoạch dựa trên sự ước lượng nhu cầu phụ tải một cách ngẫu nhiên, do đó cần xem xét tới việc kết hợp phương pháp xác suất [7]. Câu hỏi đặt ra là những mô hình tối ưu sẽ tiếp cận theo phương pháp xác suất nào là cần được xem xét? Và những đặc điểm cũng như yêu cầu cho lời giải bài toán phân bố công suất được đặt ra.1 Mô hình tối ưu theo phân bố xác suất Có một số mô hình tối ưu theo phương pháp xác suất được đề cập trong các tài liệu [4][8][14][15]. Trong đó, có thể thấy “giải bài toán phân bố công suất với mô hình xác xuất của công suất tác dụng và công suất phản kháng để xác định một kế hoạch điều độ tối ưu với các nguồn điện không ổn định được thực hiện bởi mô phỏng Monte- Carlo, là phương pháp phổ biến và chính xác cho các bài toán xác xuất” [8]. Cụ thể, mô phỏng Monte Carlo [25] với các lượt phân tích phân bố công suất với các dữ liệu máy phát và phụ tải đầu vào biến đổi, có thể được sử dụng để phân tích đặc tính ngẫu nhiên của một hệ thống điện.

Thuật toán mô phỏng này càng trở nên quan trọng đặc biệt bởi các đặc tính không thể điều khiển của các nguồn năng lượng tái tạo như điện gió và điện mặt trời [8][14].2 Yêu cầu cho lời giải đối với bài toán phân bố công suất Theo truyền thống, bài toán phân bố công suất được giải bằng các phương pháp Gauss-Seidel (GS), Newton-Raphson (NR) hoặc Fast Decoupled Power Flow (FDPF), một biến thể của Newton-Raphson. Phương pháp Newton-Raphson có độ hội tụ bậc hai nhưng cần một khối lượng tính toán rất lớn tại mỗi bước lặp, đặc biệt là với bài toán kích cỡ lớn. Phương pháp Fast Decoupled Power Flow có khối lượng tính toán tương đối ít hơn, nhưng chỉ là hội tụ tuyến tính. Trong thực tế, phương pháp Newton- HVTH: Nguyễn Sỹ Bảo Anh GVHD: PGS.

Võ Ngọc Điều Luận văn thạc sĩ 6 Raphson nói chung được ưa chuộng hơn do trong một số trường hợp Fast Decoupled Power Flow không thể hội tụ thì Newton-Raphson vẫn tiếp tục giải bài toán. Tuy nhiên, không phương pháp nào có thể đáp ứng với những bài toán cực lớn.