Tổng quan luận án

Vận hành hệ thống điện trong bối cảnh thị trường điện cạnh tranh đòi hỏi việc xác định chính xác công suất phát của các nhà máy điện và chi phí sản xuất điện năng theo từng chu kỳ thời gian thực tế, đặc biệt là kế hoạch trước 24 giờ của phụ tải. Sự tham gia ngày càng tăng của các nguồn năng lượng tái tạo, tiêu biểu là năng lượng gió, mang lại lợi ích giảm tiêu thụ nhiên liệu hóa thạch và cắt giảm phát thải khí nhà kính. Tuy nhiên, tính chất bất định của tốc độ gió dẫn đến sự dao động công suất phát, làm thay đổi trạng thái vận hành của lưới điện và đặt ra thách thức cho các bài toán tối ưu hóa vận hành. Các phương pháp giải tích truyền thống gặp khó khăn khi giải quyết các mô hình toán học phi tuyến, không liên tục, đa biến và có hàm chi phí không lồi của hệ thống điện tích hợp điện gió. Luận án giải quyết khoảng trống này bằng việc nghiên cứu và áp dụng các thuật toán bầy đàn cải tiến và thuật toán tìm kiếm thông minh nhân tạo vào chuỗi ba bài toán điều độ vận hành hệ thống điện có xét đến năng lượng gió theo chu kỳ phụ tải 24 giờ.

Mục tiêu nghiên cứu của luận án được xác định gồm ba nội dung chính:

  1. Ứng dụng và chọn lựa các thuật toán để giải bài toán điều độ tối ưu của hệ thống điện.
  2. Tính được công suất truyền tải trên đường dây, công suất và công suất phản kháng của các nhà máy phát nhiệt điện (NMPND) và nhà máy điện gió (NMĐG), dung lượng tụ bù, nấc điều chỉnh bộ điều áp tối ưu với chi phí phát điện cực tiểu theo phụ tải 24 giờ.
  3. Tính tối ưu tổn thất công suất, tối ưu độ lệch điện áp tại các nút tải, tối ưu ổn định điện áp của hệ thống có tham gia của NMĐG theo phụ tải 24 giờ.

Đối tượng nghiên cứu của đề tài bao gồm:

  • Cấu trúc và nguyên lý của các thuật toán tìm kiếm bầy đàn: Tối ưu hóa bầy đàn cơ bản (PSO), Tối ưu hóa bầy đàn với hệ số gia tốc biến đổi thời gian (PSO-TVAC), Tối ưu hóa bầy đàn với gradient giả (PG-PSO), Tối ưu hóa bầy đàn gradient giả kết hợp hệ số giới hạn (PG-PSOCF) và Thuật toán chim tu hú (CS).
  • Hàm công suất và mô hình chi phí của tua bin gió và nhà máy điện gió (NMĐG).
  • Mô hình toán học của ba bài toán điều độ: Điều độ kinh tế (ED), Phân bố tối ưu công suất (OPF) và Điều độ tối ưu công suất phản kháng (ORPD) khi có sự tham gia của NMĐG.
  • Thuật toán giải các bài toán ED, OPF, ORPD bằng các phương pháp bầy đàn theo yêu cầu phụ tải 24 giờ.

Phạm vi nghiên cứu tập trung vào:

  • Không gian tính toán: Lưới điện thử nghiệm chuẩn IEEE 30 nút được tích hợp nguồn điện gió.
  • Thời gian khảo sát: Chu kỳ vận hành 24 giờ trước thời điểm phát điện, sử dụng dữ liệu phụ tải 24 giờ và dữ liệu vận tốc gió đo đạc thực tế 24 giờ tại tỉnh Phú Yên vào ngày 16/04/2012.
  • Các chỉ tiêu kỹ thuật - kinh tế được tối ưu: Chi phí nhiên liệu và chi phí vận hành máy phát nhiệt điện kết hợp NMĐG, công suất truyền tải trên đường dây, điện áp các nút, dung lượng bù công suất phản kháng, nấc điều áp máy biến áp, tổn thất công suất tác dụng và độ lệch điện áp nút.

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

Luận án đã hệ thống hóa các công trình nghiên cứu trong và ngoài nước về phương pháp tối ưu hóa vận hành hệ thống điện và tích hợp năng lượng tái tạo:

  • Các phương pháp giải tích truyền thống: Khởi đầu từ nghiên cứu của Dommel và Tinney (1960) về bài toán phân bố tối ưu công suất (OPF), các phương pháp cổ điển như Quy hoạch tuyến tính (LP), Quy hoạch phi tuyến (NP), Phương pháp đơn hình (SA), Phương pháp nhân tử Lagrange (LMM) và Phương pháp Newton (NM) đã được ứng dụng rộng rãi. Tuy nhiên, các kỹ thuật này bộc lộ hạn chế khi giải quyết các hàm mục tiêu không lồi, có bước nhảy và hệ thống quy mô lớn.
  • Các thuật toán trí tuệ nhân tạo và tối ưu hóa bầy đàn: Từ năm 1980, các công cụ như Mạng nơ-ron (NN), Hệ mờ (FS), Tính toán tiến hóa (EC) và Thuật toán di truyền (GA) được đưa vào nghiên cứu kỹ thuật điện. Năm 1995, James Kennedy và Russell Eberhart đề xuất thuật toán Tối ưu hóa bầy đàn (PSO). Từ năm 1996 đến 2004, các biến thể PSO được phát triển nhằm cải thiện tốc độ hội tụ và tránh bẫy cực trị địa phương, bao gồm PSO-TVAC của Ratnaweera và cộng sự (2004), PG-PSO và PG-PSOCF. Đến năm 2009, Xin-She Yang và Suash Deb giới thiệu thuật toán Chim tu hú (Cuckoo Search - CS) kết hợp bước nhảy Lévy flights.
  • Các nghiên cứu giải bài toán ED: Tác giả điểm lại các thuật toán đã công bố như GA lai hệ mờ (HFGA), Thuật toán tiến hóa (EA), Thuật toán tìm kiếm hòa âm tiến hóa (EHSA), Thuật toán tiến hóa vi phân lai tìm kiếm hòa âm (HDEHS), Thuật toán tối ưu hóa đàn kiến miễn dịch tiến hóa vi phân (DEIAC), và Quy hoạch toàn phương tuần tự lai tiến hóa vi phân (HDE-SQP). Các phương pháp SGA, FGA, EA, SA, EP, WOA cũng được đối chiếu kết quả trên hệ thống chuẩn.
  • Các nghiên cứu giải bài toán OPF và ORPD: Đối với OPF, các tài liệu ghi nhận sự tham gia của thuật toán Tìm kiếm đàn kiến (ACS), Thuật toán đàn ong (BCA), Thuật toán tôi thép mô phỏng (SAM), Mạng nơ-ron kết hợp (CNN), WIPSO, MAKWO, IPSO. Đối với ORPD, các kỹ thuật như HLG, Thuật toán loài dơi (BA), EA, HPSO và QSCSA được dùng để giảm tổn thất công suất ($P_{loss}$) và độ lệch điện áp ($VD$).
  • Tích hợp năng lượng gió vào hệ thống điện: Các nghiên cứu trước xác định công suất phát của tua bin gió qua phân bố xác suất Weibull hoặc Gauss, hoặc phân bổ theo các kịch bản tỷ lệ công suất định mức (15%, 25%, 50%, 75%, 90%). Khi xây dựng hàm chi phí điện gió, các tác giả quốc tế đề xuất tính toán chi phí vận hành, bảo trì kết hợp với hàm phạt do sai lệch công suất phát thực tế so với kế hoạch đăng ký. Một số thuật toán đã áp dụng cho bài toán tích hợp gió bao gồm Thuật toán mô phỏng sinh trưởng thực vật (PGSA), GAS, PSO, GWO, KOA, BHOA, BDA, PIO và GGABC.

Khoảng trống nghiên cứu luận án lựa chọn giải quyết: Phần lớn các công trình công bố trước đây chỉ giải quyết riêng lẻ từng bài toán (chỉ xét ED, hoặc chỉ xét OPF, hoặc chỉ xét ORPD) và thường khảo sát tại một thời điểm tải tĩnh cố định. Việc mô hình hóa đồng thời cả ba bài toán ED, OPF và ORPD có tích hợp nguồn điện gió với hàm chi phí xét đến độ bất định của gió (phân bố Weibull kèm hệ số phạt thiếu/thừa công suất), giải quyết theo chu kỳ biến thiên phụ tải 24 giờ và so sánh toàn diện hiệu năng của năm thuật toán bầy đàn (PSO, PSO-TVAC, PG-PSO, PG-PSOCF, CS) chưa được thực hiện đầy đủ trong các tài liệu trước đó.

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

Luận án dựa trên cơ sở lý thuyết vận hành hệ thống điện, lý thuyết xác suất thống kê mô tả năng lượng gió và các nguyên lý tối ưu hóa metaheuristic:

  • Mô hình khí động học và công suất tua bin gió: Công suất phát điện của tua bin gió được xác định theo biểu thức: $$P_w = \frac{1}{2} \rho A C_p v_w^3$$ Trong đó $\rho$ là mật độ không khí ($\text{kg/m}^3$), $A$ là diện tích quét của cánh quạt ($\text{m}^2$), $C_p$ là hệ số công suất rotor ($C_p = 0,3 - 0,5$), $v_w$ là vận tốc gió ($\text{m/s}$). Đặc tính vận hành của tua bin được chia thành ba vùng làm việc dựa trên vận tốc gió khởi động ($v_{cut-in}$), vận tốc gió định mức ($v_{rated}$) và vận tốc gió ngắt dừng ($v_{cut-out}$).
  • Mô hình thống kê tốc độ gió và hàm chi phí gió: Vận tốc gió được mô tả bằng hàm phân phối xác suất Weibull hai thông số: hệ số hình dạng $k$ và hệ số quy mô $c$. Hàm chi phí phát điện của NMĐG gồm chi phí vận hành bảo dưỡng trực tiếp và chi phí phạt do sai lệch công suất phát so với công suất cam kết ($d_j(t)$): $$d_j(t) = C_{up,j} E(Y_{up,j}) + C_{dn,j} E(Y_{dn,j})$$ Hàm chi phí của các tổ máy nhiệt điện sử dụng hàm đa thức bậc hai truyền thống: $F_i(P_{Gi}) = a_i + b_i P_{Gi} + c_i P_{Gi}^2$.
  • Mô hình ba bài toán điều độ:
    • Bài toán ED: Cực tiểu hóa tổng chi phí nhiên liệu của các nhà máy nhiệt điện và chi phí NMĐG, thỏa mãn cân bằng công suất tác dụng $\sum P_i = P_D + P_{loss}$ (tính qua hệ số ma trận tổn thất $B$) và giới hạn công suất phát $P_{Gi}^{min} \le P_{Gi} \le P_{Gi}^{max}$.
    • Bài toán OPF: Cực tiểu hóa tổng chi phí phát điện toàn hệ thống, thỏa mãn hệ phương trình trào lưu công suất tác dụng và phản kháng phi tuyến tại các nút, giới hạn điện áp nút máy phát $V_{Gi}$, điện áp nút tải $V_L$, khả năng truyền tải công suất biểu kiến trên đường dây $S_{ij}$, và phạm vi điều chỉnh tỷ số phân áp máy biến áp $T_k$.
    • Bài toán ORPD: Tối ưu hóa phân bố công suất phản kháng nhằm cực tiểu hóa tổn thất công suất tác dụng ($P_{loss}$), cực tiểu hóa tổng độ lệch điện áp nút tải ($VD = \sum |V_k - 1,0|$), và cực đại hóa chỉ số ổn định điện áp ($L_{max}$), với các biến điều khiển là điện áp máy phát, nấc điều chỉnh đầu phân áp và dung lượng đóng cắt tụ bù $Q_{Ci}$.

Phương pháp nghiên cứu thực nghiệm mô hình toán học: Tác giả thiết lập chương trình tính toán và mô phỏng trên nền tảng phần mềm Matlab, triển khai thử nghiệm trên hệ thống điện chuẩn IEEE 30 nút. Năm thuật toán tối ưu hóa bầy đàn được xây dựng và cài đặt thông số:

  • PSO: Cập nhật vị trí và vận tốc hạt dựa trên kinh nghiệm bản thân ($pbest$) và kinh nghiệm bầy ($gbest$).
  • PSO-TVAC: Thay đổi tuyến tính hai hệ số gia tốc nhận thức ($c_1$) và xã hội ($c_2$) theo thời gian vòng lặp để chuyển trọng tâm từ thăm dò không gian sang khai thác cục bộ.
  • PG-PSO: Tích hợp thành phần hướng dẫn dốc giả (pseudo-gradient) vào công thức vận tốc nhằm định hướng hạt di chuyển nhanh về vùng nghiệm hứa hẹn.
  • PG-PSOCF: Bổ sung hệ số co hẹp (constriction factor) vào thuật toán PG-PSO nhằm kiểm soát độ phân kỳ và đảm bảo tính ổn định hội tụ.
  • CS: Mô phỏng hành vi đẻ nhờ tổ của chim tu hú, sử dụng phân phối bước nhảy Lévy flights (bước nhảy dài ngẫu nhiên) để tăng cường khả năng thoát khỏi cực trị địa phương, kết hợp cơ chế loại bỏ trứng lạ với xác suất phát hiện $P_a$.

Nguồn dữ liệu phục vụ nghiên cứu:

  • Thông số kỹ thuật đường dây, máy biến áp, giới hạn máy phát của hệ thống chuẩn IEEE 30 nút.
  • Chuỗi số liệu đo đạc vận tốc gió thực tế theo 24 giờ trong ngày 16/04/2012 tại tỉnh Phú Yên.
  • Biểu đồ phụ tải hệ thống điện biến thiên qua 24 giờ.

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

Chương 1: Mở đầu

Chương 1 trình bày lý do chọn đề tài dựa trên yêu cầu vận hành tối ưu kinh tế, an toàn và an ninh của hệ thống điện khi chuyển đổi sang cơ chế thị trường điện tại Việt Nam. Nội dung nêu bật tính cần thiết của việc tích hợp năng lượng tái tạo, cụ thể là năng lượng gió, và sự cần thiết của các thuật toán tối ưu hóa bầy đàn mới có khả năng giải bài toán phi tuyến, nhiều biến và hàm chi phí phức tạp. Chương này xác định rõ mục đích nghiên cứu, đối tượng nghiên cứu (thuật toán PSO, PSO-TVAC, PG-PSO, PG-PSOCF, CS; mô hình chi phí tua bin gió; bài toán ED, OPF, ORPD có gió), phạm vi nghiên cứu (hệ thống IEEE 30 nút, chu kỳ phụ tải 24 giờ) và phương pháp nghiên cứu kết hợp giữa cơ sở lý thuyết giải tích và mô phỏng số trên Matlab.

Chương 2: Tổng quan

Chương 2 tổng hợp bức tranh tổng thể về các nghiên cứu tối ưu hóa hệ thống điện và tiềm năng năng lượng gió. Nội dung khái quát lịch sử phát triển của các thuật toán trí tuệ nhân tạo từ năm 1980 đến các thuật toán bầy đàn hiện đại. Phần đánh giá tiềm năng điện gió Việt Nam trích dẫn số liệu khảo sát: hơn 39% diện tích lãnh thổ có tốc độ gió trung bình năm trên 6 m/s ở độ cao 65 m (tương đương quy mô công suất 512 GW), trong đó hơn 8% diện tích đạt cấp tiềm năng rất tốt.

Chương này phân tích nguyên lý làm việc của tua bin gió qua 3 vùng tốc độ, cấu trúc hàm chi phí phát điện gió bao gồm chi phí vận hành và hàm phạt do sai lệch kế hoạch. Đồng thời, tác giả tổng thuật các công trình đã công bố về ba bài toán ED, OPF, ORPD (cả trường hợp lưới truyền thống và trường hợp có NMĐG), so sánh kết quả thực nghiệm từ các tài liệu tham khảo quốc tế (sử dụng các thuật toán GA, PSO, DE, ACO, GWO, BDA, BHOA...), từ đó chỉ ra các khiếm khuyết của các nghiên cứu trước để định vị hướng giải quyết của luận án.

Chương 3: Các phương pháp tìm kiếm tối ưu

Chương 3 trình bày chi tiết cơ sở toán học, các bước thuật toán và lưu đồ tính toán của năm phương pháp tối ưu hóa bầy đàn:

  • Tối ưu hóa bầy đàn cơ bản (PSO): Trình bày biểu thức cập nhật vận tốc và tọa độ hạt.
  • Tối ưu hóa bầy đàn với hệ số gia tốc biến đổi thời gian (PSO-TVAC): Trình bày công thức điều chỉnh giá trị $c_1$ giảm dần từ 2,5 về 0,5 và $c_2$ tăng dần từ 0,5 lên 2,5 theo số lần lặp.
  • PSO với gradient giả (PG-PSO): Xây dựng công thức tính độ dốc giả định hướng cho véc-tơ chuyển động của hạt.
  • PSO gradient giả có hệ số giới hạn (PG-PSOCF): Bổ sung hệ số co hẹp nhằm đảm bảo tính ổn định của quỹ đạo hạt.
  • Thuật toán chim tu hú (CS): Trình bày đặc tính bước nhảy Lévy flights mô tả bằng phân phối xác suất đuôi nặng, các bước tìm kiếm tổ mới và cơ chế phát hiện trứng lạ của chim chủ với xác suất $P_a$.

Chương này cũng cung cấp các lưu đồ giải thuật chi tiết (Lưu đồ 3.1 cho PSO-TVAC, Lưu đồ 3.2 cho PG-PSOCF, Lưu đồ 3.3 cho CS) và xác lập bộ thông số cài đặt tốt nhất được rút ra từ các thử nghiệm chuẩn để chuẩn bị áp dụng cho các chương tính toán chuyên đề tiếp theo.

Chương 4: Điều độ kinh tế trong hệ thống điện

Chương 4 thiết lập và giải quyết bài toán ED cho hệ thống điện có và không có sự tham gia của nguồn điện gió:

  • Xây dựng mô hình toán: Thiết lập hàm mật độ xác suất tốc độ gió Weibull, xác định kỳ vọng toán của công suất phát tua bin gió, mô hình hóa hàm chi phí của NMĐG gồm chi phí đầu tư, chi phí vận hành bảo trì và hàm phạt khi công suất phát sai lệch so với kế hoạch đăng ký. Xây dựng hàm mục tiêu tổng chi phí hệ thống kết hợp ràng buộc cân bằng công suất có xét tổn thất qua ma trận hệ số $B$.
  • Dữ liệu tính toán: Sử dụng dữ liệu tốc độ gió đo đạc thực tế 24 giờ tại tỉnh Phú Yên (ngày 16/04/2012), đồ thị phụ tải biến thiên 24 giờ và hệ thống IEEE 30 nút (gồm 6 máy phát).
  • Kết quả thực nghiệm: Lập trình chạy thử nghiệm cho 24 thời điểm tải trong ngày. Chi tiết hóa kết quả tại các thời điểm tải đặc trưng (như lúc 1 giờ và lúc 24 giờ, và tại các mức tải 289,3 MW và 189 MW). Tác giả so sánh chi phí vận hành, độ lệch chuẩn và thời gian hội tụ giữa các phương pháp PSO, PSO-TVAC, PG-PSO, PG-PSOCF và CS. Đồ thị hội tụ và độ lệch chuẩn cho thấy các thuật toán cải tiến và CS đạt chi phí phát điện thấp hơn và độ ổn định nghiệm cao hơn so với PSO cơ bản.

Chương 5: Phân bố tối ưu công suất trong hệ thống điện

Chương 5 phát triển bài toán phân bố tối ưu công suất (OPF) toàn diện trên hệ thống IEEE 30 nút có kết nối NMĐG:

  • Mô hình toán: Hàm mục tiêu là cực tiểu hóa tổng chi phí vận hành toàn hệ thống, thỏa mãn đầy đủ các ràng buộc phương trình trào lưu công suất nút (AC power flow), giới hạn công suất tác dụng và phản kháng của các máy phát nhiệt điện và NMĐG, giới hạn điện áp tại tất cả các nút, khả năng tải dòng điện/công suất biểu kiến của các nhánh đường dây và phạm vi điều chỉnh nấc phân áp của máy biến áp.
  • Kịch bản khảo sát: Khảo sát hệ thống tích hợp NMĐG với các mức công suất lắp đặt khác nhau, cụ thể là trường hợp NMĐG có công suất 26 MW và 40 MW, tính toán tương ứng với các mức phụ tải lúc 1 giờ và lúc 24 giờ.
  • Kết quả và đánh giá: So sánh các chỉ số vận hành gồm công suất phát tối ưu của từng máy phát ($P_1$ đến $P_6$), công suất NMĐG, điện áp các nút máy phát, nấc điều chỉnh đầu phân áp máy biến áp ($T_{6-9}, T_{6-10}, T_{4-12}$), tổng chi phí vận hành ($/h) và thời gian giải bài toán ($t$). Đồ thị độ hội tụ nghiệm và đồ thị biểu diễn độ lệch chuẩn khẳng định CS và PG-PSOCF cho kết quả chi phí thấp hơn và tốc độ hội tụ vượt trội hơn so với PSO, GA, DE và các thuật toán đã công bố trong các tài liệu tham khảo.

Chương 6: Điều độ tối ưu công suất phản kháng trong hệ thống điện

Chương 6 tập trung giải quyết bài toán ORPD khi có NMĐG kết nối vào lưới IEEE 30 nút:

  • Xây dựng mô hình toán đa mục tiêu: Khác với ED và OPF, bài toán ORPD không đặt nặng hàm chi phí nhiên liệu mà tối ưu hóa ba chỉ tiêu kỹ thuật then chốt:
    1. Cực tiểu hóa tổn thất công suất tác dụng trên lưới điện truyền tải ($P_{loss}$).
    2. Cực tiểu hóa tổng độ lệch điện áp tại các nút tải ($VD$).
    3. Cực đại hóa chỉ số ổn định điện áp của hệ thống ($L_{max}$).
  • Ràng buộc: Ràng buộc cân bằng công suất phản kháng, giới hạn điện áp nút nguồn, giới hạn dung lượng bù tại các trạm tụ bù nút 10 đến 29, và vị trí nấc biến áp.
  • Kết quả thực nghiệm: Thực hiện tính toán tối ưu riêng biệt cho từng hàm mục tiêu ($P_{loss}$, $VD$, $L_{max}$) tại các thời điểm phụ tải 1 giờ và 24 giờ. Các bảng số liệu chi tiết và đồ thị hội tụ thể hiện rõ giá trị tổn thất công suất giảm đáng kể, chỉ số ổn định điện áp $L_{max}$ được cải thiện, và độ lệch điện áp $VD$ được thu hẹp về mức cho phép. Tác giả đối chiếu kết quả tính toán giữa các thuật toán PSO, PSO-TVAC, PG-PSO, PG-PSOCF, CS và so sánh với các tài liệu đã công bố (EA, MOPSO, ICSE, DE, QSCSA).

Chương 7: Kết luận và hướng phát triển

Chương 7 tổng kết toàn bộ các kết quả đạt được của luận án từ Chương 1 đến Chương 6. Tác giả khẳng định tính hiệu quả của các giải thuật bầy đàn cải tiến và thuật toán chim tu hú trong việc tính toán vận hành tối ưu hệ thống điện tích hợp điện gió. Chương này cũng nêu tóm lược các đóng góp mới và vạch ra các định hướng mở rộng nghiên cứu trong tương lai.

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

Bảng tổng hợp đối chiếu kết quả các bài toán qua tài liệu tổng quan và tính toán của luận án:

Bài toán Hệ thống thử nghiệm Mục tiêu tối ưu chính Các thuật toán chính được áp dụng Các biến điều khiển / Kết quả đầu ra chính
Điều độ kinh tế (ED) IEEE 30 nút + NMĐG Cực tiểu hóa tổng chi phí vận hành ($/h) PSO, PSO-TVAC, PG-PSO, PG-PSOCF, CS Công suất phát các tổ máy ($P_1 - P_6$), công suất NMĐG theo phụ tải 24 giờ.
Phân bố tối ưu công suất (OPF) IEEE 30 nút + NMĐG (26 MW, 40 MW) Cực tiểu chi phí phát điện, thỏa mãn an ninh lưới PSO, PSO-TVAC, PG-PSO, PG-PSOCF, CS Công suất $P, Q$ máy phát, điện áp nút $V$, nấc biến áp $T$, công suất truyền tải nhánh.
Điều độ tối ưu công suất phản kháng (ORPD) IEEE 30 nút + NMĐG Cực tiểu $P_{loss}$, cực tiểu $VD$, cực đại ổn định $L_{max}$ PSO, PSO-TVAC, PG-PSO, PG-PSOCF, CS Điện áp máy phát, nấc điều áp máy biến áp, dung lượng bù tụ điện ($Q_{C10}-Q_{C29}$).

Những đóng góp mới về mặt lý luận và khoa học:

  • Xây dựng mô hình toán học tích hợp hoàn chỉnh cho ba bài toán ED, OPF và ORPD trong hệ thống điện có kết nối nhà máy điện gió, trong đó mô hình chi phí điện gió kết hợp đầy đủ chi phí bảo dưỡng vận hành và hàm chi phí phạt sai lệch kế hoạch dựa trên phân bố xác suất Weibull của tốc độ gió thực tế.
  • Phát triển và cải tiến các thuật toán tối ưu hóa bầy đàn: Cải tiến PSO thông qua hệ số gia tốc biến đổi theo thời gian (PSO-TVAC), định hướng chuyển động theo gradient giả (PG-PSO) và kiểm soát hội tụ bằng hệ số giới hạn (PG-PSOCF). Luận án là một trong những công trình áp dụng thành công thuật toán Chim tu hú (Cuckoo Search - CS) kết hợp bước nhảy Lévy flights cho chuỗi ba bài toán điều độ tối ưu hệ thống điện có nguồn điện gió.

Những đóng góp mới về mặt thực tiễn:

  • Cung cấp giải pháp tính toán điều độ tối ưu theo chu kỳ phụ tải 24 giờ thực tế, cho phép các đơn vị điều độ và vận hành hệ thống điện xác định trước kế hoạch huy động công suất tối ưu, giá bán điện, dung lượng bù, nấc biến áp và thông số điện áp nhằm giảm chi phí sản xuất điện và nâng cao tính ổn định điện áp của lưới điện.
  • Ứng dụng chuỗi số liệu gió thực tế đo đạc tại tỉnh Phú Yên giúp kiểm chứng tính khả thi của mô hình toán trong điều kiện khí tượng và địa hình tại Việt Nam.

Kiến nghị và giải pháp đề xuất từ luận án:

  • Ứng dụng các thuật toán bầy đàn cải tiến, đặc biệt là PG-PSOCF và CS, vào các phần mềm tính toán điều độ tại các Trung tâm Điều độ Hệ thống điện khi thị trường điện bán buôn và cạnh tranh đi vào vận hành chính thức.
  • Cần xây dựng hệ thống cơ sở dữ liệu quan trắc tốc độ gió liên tục nhiều năm tại các khu vực tiềm năng để nâng cao độ chính xác khi xác lập các tham số phân phối Weibull và dự báo công suất phát điện gió 24 giờ trước ngày vận hành.

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

Hạn chế của nghiên cứu:

  • Toàn bộ mô phỏng tính toán mới chỉ dừng lại ở lưới điện thử nghiệm chuẩn IEEE 30 nút kết hợp với một nhà máy điện gió, chưa thử nghiệm trên mô hình lưới điện thực tế quy mô lớn của hệ thống điện lực Việt Nam với hàng trăm nút và đa điểm nút đấu nối điện gió.
  • Dữ liệu gió đo đạc phục vụ tính toán điều độ 24 giờ được thu thập trong phạm vi một ngày cụ thể (ngày 16/04/2012 tại tỉnh Phú Yên), chưa bao quát hết biến động phức tạp của thời tiết giữa các mùa mưa bão và mùa khô trong năm.
  • Nghiên cứu chỉ tập trung vào nguồn điện gió, chưa xét đến sự xuất hiện đồng thời của các nguồn năng lượng tái tạo phân tán khác như điện mặt trời hoặc hệ thống lưu trữ năng lượng pin tích năng (BESS).

Hướng phát triển tiếp theo:

  • Mở rộng ứng dụng các thuật toán PG-PSOCF và CS để tính toán điều độ tối ưu cho hệ thống điện thực tế Việt Nam có quy mô lớn hơn với nhiều nút và nhiều đường dây truyền tải phức tạp.
  • Phát triển bài toán điều độ tối ưu đa nguồn năng lượng tái tạo, kết hợp đồng thời thủy điện, nhiệt điện, điện gió, điện mặt trời và hệ thống pin lưu trữ năng lượng.
  • Nghiên cứu cải tiến thêm các thông số cài đặt của thuật toán để giảm thời gian tính toán nhằm đáp ứng yêu cầu điều độ thời gian thực (real-time dispatch).

Giá trị tham khảo

Luận án mang lại giá trị tham khảo rõ ràng cho các nhóm đối tượng:

  • Nghiên cứu sinh, học viên cao học và giảng viên chuyên ngành Mạng và Hệ thống điện: Bản luận án là tài liệu tham khảo chi tiết về phương pháp mô hình hóa toán học các bài toán điều độ ED, OPF, ORPD; cơ chế xây dựng hàm chi phí phạt sai lệch công suất gió; và kỹ thuật lập trình mô phỏng các giải thuật metaheuristic (PSO, PSO-TVAC, PG-PSO, PG-PSOCF, CS) trên phần mềm Matlab.
  • Kỹ sư điều độ và chuyên viên kỹ thuật tại Trung tâm Điều độ Hệ thống điện Quốc gia (A0), các điều độ miền và các công ty truyền tải điện: Cung cấp phương pháp luận và các biểu thức giải tích để xác định phân bố công suất kinh tế, công suất phản kháng tối ưu, nấc điều chỉnh đầu phân áp máy biến áp và dung lượng bù tụ điện nhằm giảm tổn thất lưới điện và duy trì ổn định điện áp nút.
  • Cơ quan quản lý thị trường điện và các nhà đầu tư nhà máy điện gió: Cung cấp góc nhìn về cơ chế vận hành thị trường điện ngày tới (day-ahead), phương pháp xác định chi phí vận hành và rủi ro phạt sai lệch công suất khi tham gia chào giá trên thị trường điện.

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

1. Luận án giải quyết những bài toán tối ưu nào trong hệ thống điện?
Luận án tập trung giải quyết ba bài toán điều độ tối ưu cốt lõi: Điều độ kinh tế (ED), Phân bố tối ưu công suất (OPF) và Điều độ tối ưu công suất phản kháng (ORPD) cho hệ thống điện có xét đến sự tham gia của nhà máy điện gió và chu kỳ phụ tải 24 giờ.

2. Hàm chi phí phát điện của nhà máy điện gió được mô hình hóa như thế nào?
Hàm chi phí của NMĐG gồm chi phí vận hành, bảo dưỡng trực tiếp và thành phần chi phí phạt $d_j(t)$. Chi phí phạt này phát sinh khi công suất phát thực tế (phụ thuộc vào vận tốc gió phân phối theo xác suất Weibull) cao hơn hoặc thấp hơn so với công suất phát đã đăng ký trong kế hoạch vận hành.

3. Năm thuật toán tìm kiếm bầy đàn được nghiên cứu và so sánh trong luận án là gì?
Năm thuật toán được nghiên cứu và lập trình mô phỏng gồm:

  • PSO: Tối ưu hóa bầy đàn cơ bản.
  • PSO-TVAC: Tối ưu hóa bầy đàn với hệ số gia tốc biến đổi thời gian.
  • PG-PSO: Tối ưu hóa bầy đàn với gradient giả.
  • PG-PSOCF: Tối ưu hóa bầy đàn gradient giả kết hợp hệ số giới hạn.
  • CS: Thuật toán chim tu hú sử dụng bước nhảy Lévy flights.

4. Dữ liệu thực nghiệm nào được tác giả sử dụng để kiểm chứng mô hình?
Nghiên cứu sử dụng hệ thống điện thử nghiệm chuẩn IEEE 30 nút, kết hợp với dữ liệu chuỗi thời gian tốc độ gió đo đạc thực tế trong 24 giờ tại tỉnh Phú Yên vào ngày 16/04/2012, các thông số kỹ thuật của tua bin gió và biểu đồ phụ tải biến thiên 24 giờ của hệ thống.

5. Mục tiêu tối ưu của bài toán điều độ công suất phản kháng (ORPD) trong luận án có điểm gì khác biệt so với bài toán ED và OPF?
Khác với bài toán ED và OPF tập trung vào mục tiêu cực tiểu hóa chi phí phát điện ($/h), bài toán ORPD trong Chương 6 không tính đến chi phí nhiên liệu mà tối ưu hóa ba chỉ tiêu kỹ thuật: cực tiểu hóa tổn thất công suất tác dụng ($P_{loss}$), cực tiểu hóa độ lệch điện áp nút tải ($VD$) và cực đại hóa chỉ số ổn định điện áp ($L_{max}$).

Kết luận

Luận án tiến sĩ của tác giả Lê Anh Dũng đã xây dựng hệ thống mô hình toán học và giải pháp thuật toán hoàn chỉnh để giải quyết ba bài toán điều độ vận hành hệ thống điện (ED, OPF, ORPD) có xét đến nguồn điện gió theo chu kỳ phụ tải 24 giờ. Việc ứng dụng thành công các thuật toán bầy đàn cải tiến (PSO-TVAC, PG-PSO, PG-PSOCF) và thuật toán chim tu hú (CS) trên hệ thống chuẩn IEEE 30 nút với số liệu gió thực tế tại Phú Yên đã chứng minh khả năng giảm thiểu chi phí phát điện, giảm tổn thất công suất và cải thiện độ ổn định điện áp. Các kết quả tính toán và phân tích trong luận án cung cấp tài liệu kỹ thuật có giá trị khoa học và thực tiễn cho công tác tính toán vận hành hệ thống điện trong giai đoạn phát triển thị trường điện cạnh tranh.