Bối cảnh và vấn đề nghiên cứu

Trong quá trình vận hành hệ thống điện thuộc thị trường điện cạnh tranh, nhu cầu phụ tải liên tục gia tăng trong khi khả năng mở rộng nguồn và lưới truyền tải luôn chịu nhiều giới hạn về kỹ thuật, kinh tế và môi trường. Sự mất cân bằng cục bộ giữa khả năng truyền tải và nhu cầu tiêu thụ dẫn đến hiện tượng quá tải trên một số tuyến đường dây truyền tải, gây ra nghẽn mạch hệ thống. Nghẽn mạch buộc các đơn vị điều độ phải thay đổi kế hoạch phát điện của các tổ máy, huy động các nguồn phát có chi phí cận biên cao hơn điểm cân bằng tối ưu, làm tăng chi phí sản xuất điện năng ($C_2 > C_1$) và đẩy giá thành điện năng thương phẩm lên cao.

Các giải pháp truyền thống chống quá tải như xây dựng đường dây song song, điều độ lại nguồn phát (OPF), hoặc sa thải bớt phụ tải bộc lộ nhiều hạn chế về mặt chi phí đầu tư, thời gian thực thi cũng như làm giảm độ tin cậy cung cấp điện. Việc ứng dụng các thiết bị truyền tải điện xoay chiều linh hoạt (FACTS), cụ thể là bộ bù dọc điều khiển bằng Thyristor (TCSC), cho phép kiểm soát trào lưu công suất, phân bố lại dòng công suất sang các tuyến dây chưa đầy tải. Tuy nhiên, do chi phí đầu tư thiết bị FACTS tương đối cao, không thể lắp đặt trên toàn bộ các nhánh. Vấn đề cốt lõi đặt ra là phải xác định chính xác tập hợp các nhánh xung yếu nhất (nút thắt cổ chai - bottleneck) thường xuyên bị đe dọa quá tải để chọn vị trí và dung lượng bù tối ưu.

Mục tiêu và nhiệm vụ nghiên cứu của đồ án được xác định cụ thể theo 4 nội dung:

  1. Tìm hiểu về thuật toán tìm luồng công suất cực đại và lý thuyết mặt cắt tối thiểu (Max-flow Min-cut).
  2. Nghiên cứu nguyên lý cấu tạo, hoạt động và xây dựng mô hình toán học của bộ bù dọc điều khiển bằng Thyristor (TCSC).
  3. Xác định vị trí và dung lượng bù tối ưu của TCSC trên lưới điện truyền tải nhằm phân bố lại dòng công suất chống quá tải.
  4. Xây dựng chương trình tính toán trên môi trường MATLAB (kết nối Excel), ứng dụng mô phỏng kiểm chứng trên phần mềm PowerWorld đối với các lưới điện cụ thể và so sánh hiệu quả kinh tế - kỹ thuật với các phương pháp điều độ truyền thống.

Đối tượng nghiên cứu của đồ án là lưới điện truyền tải điện xoay chiều có xét đến giới hạn truyền tải của các nhánh và thiết bị bù TCSC. Phạm vi nghiên cứu tập trung vào bài toán phân bố trào lưu công suất xác lập, thuật toán xác định mặt cắt tối thiểu trên đồ thị mạng điện chuyển đổi, và thực nghiệm tính toán trên sơ đồ lưới điện mẫu 7 nút (7-bus) và 14 nút (14-bus).

Cơ sở lý thuyết và phương pháp

Đồ án xây dựng trên các nền tảng lý thuyết tối ưu hóa kinh tế hệ thống điện, công nghệ thiết bị FACTS và lý thuyết đồ thị mạng:

Khung lý thuyết kinh tế vận hành và phân bố công suất: Chi phí phát điện của tổ máy thứ $i$ được mô hình hóa theo hàm bậc hai: $$C_i(P_{gi}) = \alpha_{0i} P_{gi}^2 + \alpha_{1i} P_{gi} + \alpha_{2i}$$ Mục tiêu vận hành tối ưu kinh tế khi chưa xét giới hạn truyền tải là tối thiểu hóa tổng chi phí phát điện: $$C_1 = \text{Minimum} \sum C_i(P_{gi})$$ Thỏa mãn cân bằng công suất $\sum P_{Gi} = \sum P_L$. Khi xảy ra nghẽn mạch đường dây, việc điều độ lại công suất phát $P'{gi}$ làm tổng chi phí tăng lên thành $C_2 = \sum C_i(P'{gi}) > C_1$.

Mô hình thiết bị bù dọc điều khiển bằng Thyristor (TCSC): Cấu tạo TCSC gồm tụ điện cố định $C$ mắc song song với cuộn kháng $L$ nối tiếp hai van Thyristor đấu song song ngược chiều (TCR). Trở kháng tương đương của TCSC thay đổi theo góc dẫn $\sigma = 2(\pi - \alpha)$ của Thyristor: $$X_{TCSC}(\alpha) = X_C - \frac{X_C^2}{X_C - X_L} \frac{\sigma + \sin\sigma}{\pi} + \frac{4X_C^2 \cos^2(\sigma/2)}{(X_C - X_L)(k^2 - 1)} \frac{k\tan(k\sigma/2) - \tan(\sigma/2)}{\pi}$$ Khi lắp đặt TCSC có dung kháng bù $-jx_c$ trên nhánh $i-j$, dòng công suất tác dụng truyền qua nhánh (bỏ qua điện trở đường dây $r \ll x$) được nâng cao theo biểu thức: $$P_{ij}^c = \frac{V_i V_j}{x_{ij} - x_c} \sin\theta_{ij}$$

Lý thuyết mặt cắt tối thiểu và luồng cực đại (Max-flow Min-cut): Phát triển từ nguyên lý Ford & Fulkerson (1956), khả năng truyền tải công suất cực đại từ nút nguồn ($s$) đến nút tải ($t$) bị chặn trên bởi thông lượng của lát cắt nhỏ nhất (Min-cut). Đồ án mô hình hóa mạng điện $n$ nút thành sơ đồ mạng tương đương gồm $n+2$ nút bằng cách gom toàn bộ công suất phát của các máy phát về nút nguồn giả định $s$ và gom toàn bộ phụ tải tiêu thụ về nút tải giả định $t$.

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

  • Phương pháp phân tích - tổng hợp và mô hình hóa toán học hệ thống điện.
  • Phương pháp lý thuyết đồ thị: Thành lập ma trận thông lượng $A_{n \times n}$, lập trình thuật toán tìm mặt cắt tối thiểu trên phần mềm MATLAB.
  • Phương pháp mô phỏng số trào lưu công suất bằng phần mềm PowerWorld.
  • Nguồn dữ liệu: Thông số kỹ thuật đường dây, máy phát, phụ tải của sơ đồ lưới chuẩn 7 nút và 14 nút (IEEE 14-bus), cùng dữ liệu thống kê chi phí đầu tư các thiết bị FACTS.

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

Chương 1: Tổng quan

Chương 1 trình bày bối cảnh cạnh tranh trong thị trường điện và các phương thức nâng cao khả năng truyền tải nhằm giảm chi phí sản xuất điện năng. Tác giả phân tích các giải pháp xử lý nghẽn mạch truyền thống gồm điều độ kế hoạch nguồn phát (OPF), điều độ cắt giảm tải, và mở rộng lưới truyền tải; đồng thời chỉ ra các nhược điểm về mặt kinh tế, thời gian thu hồi vốn và độ tin cậy.

Chương này tổng hợp các thiết bị FACTS phổ biến (SVC, STATCOM, SSSC, UPFC, TCSC) và đưa ra bảng so sánh chi phí đầu tư trên một đơn vị công suất bù:

Thiết bị bù Chi phí đầu tư (USD/kVAr)
Tụ bù song song 8
Tụ bù nối tiếp 20
TCSC 40
SVC 40
STATCOM 50
UPFC 50

Tác giả phân tích các nghiên cứu trước đây vốn phụ thuộc vào giải thuật Di truyền (GA) làm tăng không gian tìm kiếm, hoặc phương pháp thử sai (Trial and Error) và phân vùng chủ quan của M. Haghifam. Từ đó, đề xuất sử dụng lý thuyết lát cắt cực tiểu (Min-cut) để xác định tập nhánh nghẽn mạch. Trên mô hình 2 nút ví dụ với các lát cắt $f_1 = 70$, $f_2 = 110$, $f_3 = 85$, $f_4 = 85$, vị trí lát cắt cực tiểu chỉ ra đúng nút thắt cổ chai của hệ thống. Tác giả phân loại 5 trường hợp vị trí lát cắt:

Trường hợp Biểu thức lát cắt cực tiểu Nội dung quy hoạch / vận hành
1 $f = \sum c_{si}$ Lát cắt thuộc tập nguồn phát; cần quy hoạch nguồn
2 $f = \sum c_{it}$ Lát cắt thuộc tập đỉnh thu; cần phân phối lại phụ tải
3 $f = \sum c_{ij}$ Lát cắt thuộc các nhánh truyền tải; cần lắp FACTS hoặc mở rộng đường dây
4 $f = \sum c_{si} + \sum c_{ij}$ Lát cắt thuộc cả nguồn và đường dây; quy hoạch nguồn và đường dây
5 $f = \sum c_{it} + \sum c_{ij}$ Lát cắt thuộc cả tải và đường dây; quy hoạch tải và đường dây

Chương 2: Mô hình toán học

Chương 2 thiết lập bài toán phân bố trào lưu công suất khi có và không có thiết bị bù TCSC. Tác giả chỉ ra rằng khi lắp TCSC với giá trị điện kháng dung $-jx_c$, tổng trở nhánh giảm xuống $x_{ij} - x_c$, làm tăng trào lưu công suất qua nhánh có bù và giảm tải cho các nhánh song song khác trong cùng mạch vòng.

Chương này phân tích giải thuật xác định lát cắt cực tiểu của Mechthild Stoer và Frank Wagner trên đồ thị vô hướng, chỉ ra hai hạn chế khi áp dụng vào lưới điện:

  1. Số phép chọn cặp đỉnh $s-t$ là $A_n^2 = n(n-1)$, dẫn đến bùng nổ không gian và thời gian tính toán.
  2. Kết quả không đảm bảo phân tách hoàn toàn nguồn phát $s$ và phụ tải $t$ ở hai nửa mặt phẳng cắt.

Để giải quyết, đồ án đề xuất lưu đồ giải thuật Max-flow Min-cut cải tiến dành riêng cho hệ thống điện $n+2$ nút. Thuật toán khởi tạo ma trận thông lượng $A_{n \times n}$ (với $a_{ii}=0$, $a_{ij}=a_{ji}$ là giới hạn truyền tải nhánh $i-j$), xét tập nguồn $S={s}$ và tập tải $T={t}$, tính tổng các thông lượng nhánh $C_{si}$ và $C_{it}$ vào danh sách listcut, từ đó xác định lát cắt có tổng dung lượng nhỏ nhất đi qua các nhánh đường dây thực tế. Chương 2 cũng trình bày quy trình tính toán dung lượng bù $x_c$ của TCSC và cơ chế liên kết dữ liệu giữa MATLAB và Microsoft Excel để tự động hóa quá trình tính toán.

Chương 3: Các ví dụ mô phỏng cụ thể

Chương 3 tiến hành tính toán và mô phỏng trên hai hệ thống điện:

  • Lưới điện 7 nút (7-bus): Chương trình Mincut_Maxflow trên MATLAB nhận diện chính xác tập hợp nhánh xung yếu gây nghẽn mạch khi phụ tải tăng. Kết quả mô phỏng trên PowerWorld cho thấy khi chưa có TCSC, việc điều độ máy phát để tránh quá tải đường dây làm tăng đáng kể chi phí phát điện. Sau khi đặt TCSC có dung lượng tính toán tối ưu vào nhánh thuộc mặt cắt cực tiểu, trào lưu công suất được điều chuyển hợp lý, loại bỏ hoàn toàn quá tải và đưa chi phí vận hành về gần mức tối ưu $C_1$.
  • Lưới điện 14 nút (14-bus): Đồ án cung cấp bảng thông số công suất nguồn phát, phụ tải tiêu thụ và thông số đường dây. Khi tăng tải 10%, hệ thống xuất hiện các nhánh quá tải nghiêm trọng. So sánh giữa phương án điều độ nguồn phát chấp nhận quá tải (chi phí thấp nhưng đe dọa an ninh hệ thống), điều độ nguồn phát không chấp nhận quá tải (chi phí sản xuất điện tăng vọt), và phương án lắp đặt TCSC tại vị trí mặt cắt tối thiểu. Phương án TCSC giải quyết dứt điểm quá tải, đồng thời tác giả thiết lập đồ thị phân tích thời gian hoàn vốn đầu tư của thiết bị bù TCSC.

Chương 4: Kết luận và hướng phát triển đề tài

Chương 4 tổng kết tính khả thi của giải thuật kết hợp giữa Max-flow Min-cut và mô hình TCSC, khẳng định phương pháp giải quyết đồng thời cả hai bài toán kinh tế và kỹ thuật. Đồng thời, tác giả đề xuất các hướng mở rộng nghiên cứu trong tương lai.

Kết quả và đóng góp

Kết quả chính:

  • Xây dựng thành công thuật toán Max-flow Min-cut cải tiến trên môi trường MATLAB, tự động hóa nhận diện tập nhánh xung yếu nhất (nút thắt cổ chai) trong lưới điện từ ma trận thông lượng đường dây.
  • Thiết lập mô hình toán và quy trình xác định vị trí, dung lượng bù của thiết bị TCSC trên các nhánh thuộc lát cắt tối thiểu, điều khiển trào lưu công suất trên lưới điện truyền tải.
  • Kiểm chứng hiệu quả thành công trên sơ đồ 7 nút và 14 nút qua phần mềm PowerWorld: loại bỏ hoàn toàn hiện tượng quá tải đường dây khi tăng tải (bao gồm trường hợp tăng tải 10% trên lưới 14 nút), giảm chi phí vận hành từ mức $C_2$ về sát mức tối ưu $C_1$, và xác định được thời gian thu hồi vốn đầu tư TCSC.

Giải pháp và kiến nghị của tác giả:

  • Ứng dụng thiết bị bù dọc TCSC (chi phí 40 USD/kVAr) thay vì đầu tư các thiết bị chi phí cao hơn như UPFC (50 USD/kVAr) hoặc STATCOM (50 USD/kVAr) cho bài toán chống quá tải truyền tải.
  • Sử dụng mô hình mạng $n+2$ nút để khoanh vùng và cô lập phạm vi tìm kiếm nhánh nghẽn mạch, tránh việc thử nghiệm dàn trải trên toàn bộ các nhánh của hệ thống điện.

Đóng góp mới của đồ án:

  • Khắc phục nhược điểm chiếm dụng bộ nhớ lớn và thời gian xử lý chậm của giải thuật Di truyền (GA) cũng như tính mò mẫm của phương pháp thử sai (Trial and Error).
  • Khắc phục tính chủ quan của phương pháp phân vùng nguồn - tải của M. Haghifam bằng thuật toán phân chia đồ thị toán học chặt chẽ.
  • Chương trình MATLAB có tính tổng quát cao, có khả năng mở rộng để phân tích cho các mô hình lưới điện có cấu trúc bất kỳ.

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

Hạn chế của đồ án:

  • Mô hình chuyển đổi $n+2$ nút giả định nguồn phát và phụ tải quy về hai nút duy nhất $s$ và $t$, do đó thuật toán loại bỏ việc phân tích các lát cắt rơi hoàn toàn vào tập nguồn phát hoặc tập đỉnh thu (Trường hợp 1 và Trường hợp 2).
  • Các công thức tính toán dòng công suất có TCSC sử dụng giả thuyết bỏ qua điện trở đường dây ($r \ll x$) để đơn giản hóa quá trình lập trình.
  • Nghiên cứu chỉ tập trung vào trạng thái xác lập và giới hạn nhiệt tĩnh của đường dây, chưa tích hợp đầy đủ các ràng buộc động về ổn định góc pha và ổn định quá độ phức tạp.

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

  • Mở rộng thuật toán cho các hệ thống điện quy mô lớn, nhiều nút và đa vùng liên kết thực tế.
  • Tích hợp thêm các ràng buộc vận hành chi tiết về điện áp nút, góc lệch pha và giới hạn công suất phản kháng vào mô hình tối ưu.
  • Nghiên cứu phối hợp đa thiết bị FACTS (kết hợp TCSC với SVC, STATCOM) để vừa chống quá tải vừa nâng cao ổn định điện áp toàn diện.

Giá trị tham khảo

Đồ án là tài liệu tham khảo hữu ích cho các đối tượng và mục đích cụ thể sau:

  • Sinh viên và học viên ngành Kỹ thuật Điện, Kỹ thuật Hệ thống điện: Tham khảo phương pháp mô hình hóa toán học thiết bị FACTS (đặc biệt là TCSC), nguyên lý điều khiển trào lưu công suất và cách tiếp cận bài toán phân bố công suất tối ưu (OPF).
  • Kỹ sư vận hành và quy hoạch hệ thống điện: Ứng dụng quy trình chuyển đổi sơ đồ mạng $n \to n+2$ nút và giải thuật Max-flow Min-cut để sàng lọc nhanh các nhánh xung yếu có nguy cơ nghẽn mạch cục bộ khi lập phương thức vận hành hoặc nâng cấp lưới.
  • Phần nội dung đáng tham khảo nhất: Chương 2 (xây dựng ma trận thông lượng và lưu đồ giải thuật mặt cắt tối thiểu) cùng Chương 3 (phương pháp phối hợp giữa MATLAB và PowerWorld để giải bài toán kỹ thuật kết hợp phân tích kinh tế).

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

1. Tại sao đồ án lựa chọn thiết bị TCSC để chống quá tải thay vì các thiết bị FACTS khác như UPFC hay STATCOM? Theo số liệu thống kê chi phí trong đồ án (Bảng 1.1), TCSC có suất đầu tư 40 USD/kVAr, thấp hơn UPFC và STATCOM (50 USD/kVAr). Đồng thời, TCSC mắc nối tiếp trên đường dây cho phép can thiệp trực tiếp làm thay đổi điện kháng tương đương nhánh ($x_{ij} - x_c$), từ đó điều khiển hiệu quả dòng công suất tác dụng qua nhánh xung yếu mà vẫn đảm bảo tính kinh tế hơn so với các thiết bị bù phối hợp phức tạp.

2. Phương pháp chuyển đổi lưới điện $n$ nút thành sơ đồ $n+2$ nút được thực hiện theo nguyên tắc nào? Tất cả công suất phát của các máy phát tại các thanh cái được gom về một nút nguồn giả định $s$. Toàn bộ phụ tải tiêu thụ tại các thanh cái được gom về một nút tải giả định $t$. Mạng điện $n$ nút trung gian ban đầu cùng với hai nút giả định $s$ và $t$ tạo thành một mạng đồ thị khép kín $n+2$ nút, cho phép áp dụng thuật toán tìm lát cắt phân chia nguồn - tải.

3. Nhược điểm chính của giải thuật Stoer - Wagner khi áp dụng vào hệ thống điện là gì? Giải thuật Stoer - Wagner áp dụng cho đồ thị vô hướng tổng quát phải duyệt qua $A_n^2 = n(n-1)$ khả năng chọn cặp đỉnh $s-t$, gây bùng nổ thời gian tính toán. Ngoài ra, lát cắt tìm được có thể không chia tách độc lập nguồn $s$ và tải $t$, khiến luồng công suất vẫn có thể truyền qua các nhánh khác mà lát cắt không đi qua, không phản ánh đúng bản chất truyền tải điện năng.

4. Tại sao đồ án loại bỏ các lát cắt cực tiểu rơi vào tập nguồn phát ($f = \sum c_{si}$) và tập đỉnh thu ($f = \sum c_{it}$)? Vì hai nút $s$ và $t$ là các nút giả định trong mô hình tính toán, trong khi trên thực tế nguồn phát và phụ tải phân bố rải rác tại các thanh cái. Lát cắt rơi vào nguồn hoặc tải thuộc về bài toán quy hoạch nguồn và phân bổ nhu cầu phụ tải, không thuộc phạm vi xử lý nghẽn mạch trên lưới truyền tải bằng thiết bị FACTS.

5. Việc lắp đặt TCSC tại vị trí mặt cắt tối thiểu mang lại lợi ích kinh tế như thế nào so với điều độ nguồn phát (OPF)? Khi xảy ra quá tải, điều độ OPF truyền thống buộc phải giảm phát ở các tổ máy giá rẻ và tăng phát ở các tổ máy giá đắt để giảm dòng qua nhánh nghẽn, làm tăng tổng chi phí phát điện lên mức $C_2$. Lắp đặt TCSC giúp tăng khả năng truyền tải của nhánh nghẽn và san tải sang các nhánh khác, cho phép hệ thống tiếp tục vận hành ở điểm phát điện tối ưu chi phí $C_1$, khoản chênh lệch chi phí $(C_2 - C_1)$ giúp nhanh chóng hoàn vốn đầu tư TCSC.

Kết luận

Đồ án tốt nghiệp "Nghiên cứu phương pháp mặt cắt tối thiểu chống quá tải trong hệ thống điện" của tác giả Trần Xuân Bách và Nguyễn Việt Hưng đã giải quyết thành công bài toán xác định điểm nghẽn mạch và vị trí bù tối ưu trên lưới truyền tải. Việc kết hợp lý thuyết Max-flow Min-cut cải tiến trên nền tảng MATLAB với mô phỏng trào lưu công suất trên PowerWorld chứng minh tính hiệu quả vượt trội trong việc khoanh vùng nhánh xung yếu so với các phương pháp GA hay thử sai trước đây. Ứng dụng thiết bị bù TCSC tại các nhánh thuộc lát cắt tối thiểu không chỉ loại bỏ dứt điểm tình trạng quá tải đường dây khi nhu cầu phụ tải tăng cao mà còn đảm bảo tối ưu hóa chi phí sản xuất điện năng trong thị trường điện cạnh tranh.