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

Trong bối cảnh của cuộc Cách mạng Công nghiệp lần thứ tư, việc ứng dụng các phần mềm tự động hóa nhằm hỗ trợ con người trong tính toán và thiết kế kỹ thuật trở nên cần thiết. Đối với ngành kỹ thuật điện và bộ môn Điện công nghiệp, quá trình tính toán và thiết kế hệ thống điện bao gồm nhiều bài toán phức tạp với khối lượng dữ liệu lớn, liên kết chặt chẽ. Các bài toán chuyên ngành như tính toán trở kháng ngắn mạch, lựa chọn tiết diện cáp dẫn, bù công suất phản kháng hay lựa chọn máy biến áp đòi hỏi chuỗi công thức dài; sai sót ở một thông số ban đầu có thể dẫn đến sai lệch toàn bộ kết quả phía sau. Do đó, việc ứng dụng phần mềm chuyên dụng như ETAP (Electrical Transient Analyzer Program) giúp rút ngắn thời gian tính toán, đồng thời cung cấp công cụ kiểm tra, đối chiếu và mô phỏng chính xác so với phương pháp tính tay truyền thống.

Mục tiêu tổng quát của đồ án là nghiên cứu, khám phá các tiện ích của phần mềm ETAP trong tính toán, thiết kế và phân tích hệ thống điện, kết hợp các kiến thức chuyên ngành cung cấp điện và hệ thống điện để giải quyết 5 nhiệm vụ/vấn đề kỹ thuật cụ thể:

  1. Tính toán chọn cáp điện phân phối: Xác định tiết diện dây dẫn, phương thức lắp đặt và các thông số ảnh hưởng đến việc lựa chọn cáp trong mạng điện phân phối.
  2. Ứng dụng năng lượng tái tạo trong phân tích ổn định hệ thống điện: Tích hợp nguồn phát điện mặt trời vào mạng lưới để đánh giá mức độ ảnh hưởng đến độ ổn định của hệ thống.
  3. Lựa chọn lưới nối đất: Tính toán và thiết kế hệ thống nối đất an toàn, đáp ứng các tiêu chuẩn kỹ thuật và tối ưu hóa chi phí kinh tế.
  4. Phân bố công suất: Khảo sát, phân tích cơ chế vận hành và giải thuật tính toán phân bố công suất (Load Flow Analysis).
  5. Tính toán ngắn mạch: Phân tích, mô phỏng diễn biến và các thành phần dòng điện khi xảy ra sự cố ngắn mạch trong hệ thống điện.
  • Đối tượng nghiên cứu: Các thuật toán, công cụ mô phỏng của phần mềm ETAP và các mô hình hệ thống điện (hệ thống truyền tải, phân phối, hệ thống chuẩn IEEE 9-bus, trạm biến áp hạ áp 22kV/380V/220V, hệ thống tiếp địa).
  • Phạm vi nghiên cứu: Mô phỏng tĩnh và quá độ trên nền tảng phần mềm ETAP phiên bản ứng dụng cho các bài toán phân bố công suất, ngắn mạch theo tiêu chuẩn IEC 60909, ổn định hệ thống điện có nguồn phân tán (PV), định cỡ dây dẫn và lưới nối đất theo tiêu chuẩn IEEE/IEC; dữ liệu khảo sát được thực hiện tại Trường Đại học Sư phạm Kỹ thuật TP. Hồ Chí Minh (hoàn thành tháng 7/2019).

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

Khung lý thuyết và mô hình tính toán

Đồ án xây dựng nền tảng lý thuyết dựa trên các khối kiến thức chuyên ngành hệ thống điện và cung cấp điện, bao gồm:

  • Phân loại nút trong bài toán phân bố công suất:
    • Nút cân bằng (Slack/Swing bus): Cung cấp trước độ lớn điện áp $U$ và góc lệch pha $\delta$; thường chỉ có một nút cân bằng trong hệ thống.
    • Nút máy phát ($P-U$ bus / Voltage-controlled bus): Biết trước công suất thực $P$ và điện áp $U$; cần xác định góc pha $\delta$ và công suất phản kháng $Q$.
    • Nút phụ tải ($P-Q$ bus / Load bus): Biết trước công suất $P$ và $Q$; cần tìm điện áp $U$ và góc pha $\delta$.
  • Các giải thuật giải hệ phương trình phi tuyến trong phân bố công suất:
    • Phương pháp Newton-Raphson: Sử dụng ma trận Jacobian ($J_1$ đến $J_4$) để giải lặp phương trình dòng tải; có tốc độ hội tụ bậc hai nhanh. Tiêu chuẩn hội tụ công suất thực và phản kháng thường đặt ở mức $0{,}001\text{ MW}$ và $\text{Mvar}$. Thuật toán phụ thuộc vào giá trị điện áp ban đầu, do đó ETAP thực hiện vài bước lặp Gauss-Seidel sơ bộ để khởi tạo giá trị.
    • Phương pháp Adaptive Newton-Raphson: Cải tiến từ Newton-Raphson bằng cách chia nhỏ bước tăng khi gặp điều kiện phân kỳ tiềm năng, kết hợp nội suy/ngoại suy tuyến tính.
    • Phương pháp tách nhanh (Fast-Decoupled Method): Đơn giản hóa từ Newton-Raphson dựa trên đặc tính góc pha ít ảnh hưởng đến $Q$ và biên độ điện áp ít ảnh hưởng đến $P$, tách thành hai hệ phương trình giải lặp riêng biệt.
    • Phương pháp gia tốc Gauss-Seidel (Accelerated Gauss-Seidel): Giải lặp phương trình điện áp nút dựa trên ma trận dẫn nạp $Y_{BUS}$. Tiêu chí hội tụ kiểm tra dung sai điện áp (đặt ở mức $0{,}000001\text{ pu}$), sử dụng hệ số gia tốc từ $1{,}2$ đến $1{,}7$ (thường chọn $1{,}45$).
  • Tiêu chuẩn kỹ thuật áp dụng: Tiêu chuẩn quốc tế IEC (như IEC 60909 cho tính toán dòng ngắn mạch), hệ thống tiêu chuẩn ANSI và các tiêu chuẩn IEEE trong thiết kế tiếp địa và phân tích ổn định.

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

  • Phương pháp tổng hợp lý thuyết: Phân tích các công thức giải tích trong môn học Cung cấp điện và Hệ thống điện; đối chiếu với tài liệu hướng dẫn kỹ thuật của ETAP.
  • Phương pháp mô hình hóa và mô phỏng số:
    • Xây dựng sơ đồ đơn tuyến (One-line diagram) trên môi trường ETAP.
    • Thiết lập thông số kỹ thuật cho nguồn, máy phát, đường dây, máy biến áp, tải, tấm pin mặt trời (PV array), nghịch lưu (Inverter) và hệ thống nối đất.
    • Thiết lập kịch bản tính toán (Study case), chạy mô phỏng, ghi nhận dữ liệu đồ thị và bảng báo cáo kết quả.
  • Nguồn dữ liệu: Thư viện tích hợp của phần mềm ETAP (thông số cáp, máy biến áp, rơle, cầu chì từ các nhà sản xuất thực tế), dữ liệu mô hình kiểm chuẩn IEEE 9-bus và các bài toán mẫu hệ thống điện công nghiệp.

Thiết kế và triển khai hệ thống mô phỏng kỹ thuật

Phần mềm ETAP (do Operation Technology, Inc. - OTI, Irvine, California, Mỹ phát triển) được tổ chức trên nền tảng cơ sở dữ liệu ba chiều kết nối qua chuẩn ODBC (Microsoft Access). Cấu trúc quản lý dự án bao gồm 3 thành phần trực giao:

  • Hiển thị (Presentation): Cho phép tạo không giới hạn các giao diện đồ họa độc lập của sơ đồ đơn tuyến (sơ đồ trở kháng, kết quả tính toán, đồ thị).
  • Trạng thái cấu hình (Configuration): Quản lý trạng thái đóng/cắt của các thiết bị đóng cắt và chế độ vận hành của tải.
  • Dữ liệu hiệu chỉnh (Revision Data): Theo dõi các thay đổi, điều chỉnh đặc tính kỹ thuật của phần tử mà không cần sao chép lại toàn bộ dự án.

Hệ thống thanh công cụ trong ETAP gồm:

  • Thanh công cụ File, Edit, View, Project, Tool: Phục vụ quản trị file (*.OTI), tùy biến góc quay phần tử ($0^\circ, 90^\circ, -90^\circ, 180^\circ$), lựa chọn hệ ký hiệu chuẩn IEC hoặc ANSI.
  • Thanh công cụ hệ thống (System Toolbar): Chuyển đổi giữa các phân hệ: Sơ đồ đơn tuyến, Hệ thống dây ngầm, Hệ thống nối đất (Ground Grid), Hệ thống cáp dây, GIS, sơ đồ điều khiển AC/DC, Study wizard và Scenario wizard.

Quy trình nhập liệu và khai báo phần tử kỹ thuật:

  • Nguồn điện (Power Grid): Khai báo điện áp định mức (Rated kV), chế độ nối đất (Grounding: TN-C, TN-S, TN-CS cho lưới $\le 1\text{ kV}$), công suất ngắn mạch 3 pha/1 pha ($MVA_{sc}$), tỷ số $X/R$ và trở kháng thứ tự nghịch/thứ tự không.
  • Máy biến áp (Transformer): Nhập điện áp sơ/thứ cấp, dòng đầy tải FLA, công suất định mức, trở kháng thứ tự thuận và thứ tự không theo phần trăm, cấp cách điện và điều kiện làm việc (độ cao, nhiệt độ môi trường).
  • Cáp dẫn (Cable): Nhập chiều dài, chọn từ thư viện thông số vật liệu (Cu/Al), cách điện, số lõi, kích thước, điện trở và điện kháng đơn vị ($\Omega/\text{km}$ hoặc $\Omega/\text{ft}$).
  • Cầu chì & Thiết bị bảo vệ: Truy xuất thư viện theo nhà sản xuất, kiểu dáng (Model), điện áp cực đại và đặc tính thời gian - dòng điện.

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

Chương 1: Giới thiệu về phần mềm ETAP

Chương này trình bày lịch sử hình thành của ETAP từ phần mềm tính toán lưới điện tĩnh (off-line) đến giải pháp thời gian thực (Real-time) từ năm 1992 do OTI phát triển. Nội dung tập trung mô tả chi tiết giao diện đồ họa, cách thức quản lý dự án, cấu trúc cơ sở dữ liệu đa chiều, các thanh công cụ chức năng (File, Edit, View, Project, Tool, System) và hướng dẫn khai báo các thuộc tính kỹ thuật cơ bản cho nguồn, máy biến áp, đường cáp và thiết bị bảo vệ.

Chương 2: Phân bố công suất

Chương 2 trình bày cơ sở toán học của bài toán dòng tải phi tuyến và so sánh 4 giải thuật: Newton-Raphson, Adaptive Newton-Raphson, Fast-Decoupled và Accelerated Gauss-Seidel. Tác giả thực hiện tính toán mô phỏng trên sơ đồ mạng điện 7 bus.

Bảng thông số đường dây và phụ tải hệ thống 7 bus mô phỏng:

Nhánh (Bus - Bus) $R$ ($\Omega$) $X$ ($\Omega$) Nút Máy phát ($MW$) Phụ tải ($MW$) Phụ tải ($M\text{var}$) Ghi chú
1 - 2 15 46 1 160 55 15 -
1 - 5 12 38 2 175 95 35 -
2 - 4 10 30 3 195 85 25 -
2 - 3 13 39 4 0 195 50 Sụt áp không ổn định
4 - 5 11 34 5 0 185 45 -
3 - 5 16 48 6 0 135 40 -
3 - 6 15 46 7 Swing 75 20 Nút cân bằng (Gen 4)
6 - 7 17 50 - - - - -
5 - 7 13 39 - - - - -

Kết quả mô phỏng cho thấy điện áp tại các bus cơ bản duy trì ổn định, riêng Bus 4 bị sụt giảm điện áp ngoài ngưỡng cho phép do phụ tải tập trung lớn ($195\text{ MW} + j50\text{ Mvar}$), đòi hỏi phải tái cấu trúc phụ tải hoặc tăng tiết diện dây dẫn.

Chương 3: Tính toán ngắn mạch

Nội dung chương xây dựng trên cơ sở tiêu chuẩn IEC 60909. Chương phân tích phương pháp tính toán các thành phần sự cố: dòng ngắn mạch đối xứng ban đầu ($I_k''$), dòng điện đỉnh ngắn mạch ($i_p$), thành phần không chu kỳ DC ($i_{dc}$), thành phần chu kỳ AC và dòng cắt ngắn mạch đối xứng ($I_b$). Đồ án hướng dẫn thiết lập hộp thoại Short-Circuit Study Case, chọn các chế độ sự cố 3 pha, 1 pha - đất, 2 pha chạm đất và xuất báo cáo kết quả dòng ngắn mạch để kiểm tra khả năng chịu đựng của khí cụ điện.

Chương 4: Phân tích ổn định quá độ

Chương 4 nghiên cứu các dạng ổn định trong hệ thống điện: ổn định góc rotor máy phát và ổn định điện áp khi tích hợp nguồn năng lượng tái tạo. Mô hình áp dụng là hệ thống chuẩn IEEE 9-bus gồm 3 máy phát, được tích hợp thêm các tổ hợp pin quang điện (PV array) qua bộ nghịch lưu (Inverter) và máy biến áp nâng áp tại Bus 6 và Bus 8.

Mô phỏng khảo sát các kịch bản sự cố và theo dõi:

  • Biến thiên góc rotor máy phát G2 theo thời gian.
  • Dao động điện áp tại Bus 4, Bus 5, Bus 6, Bus 7 và Bus 8 dưới các mức độ thâm nhập điện mặt trời khác nhau.
  • Tổn thất công suất thực ($MW$) và công suất phản kháng ($M\text{var}$) trên các thanh cái từ Bus 1 đến Bus 6.

Chương 5: Tính toán lựa chọn dây dẫn

Chương 5 thiết lập bài toán định cỡ cáp điện (Cable Sizing) cho mạng hạ thế và đường dây phân phối. Mô hình thử nghiệm gồm nguồn lưới qua trạm biến áp $22\text{ kV}/380\text{ V}$ và $380\text{ V}/220\text{ V}$ cấp điện cho hệ thống động cơ. Phương pháp tính toán tích hợp các bảng tra tiêu chuẩn về dòng tải liên tục (Ampacity), hệ số hiệu chỉnh nhiệt độ môi trường, cách thức lắp đặt, độ sụt áp cực đại cho phép và khả năng chịu dòng ngắn mạch để tự động đề xuất tiết diện pha (Sizing-phase) trong trang Sizing của ETAP.

Chương 6: Hệ thống nối đất

Chương 6 hướng dẫn thiết kế lưới tiếp địa bằng mô-đun Ground Grid System trong ETAP. Tác giả trình bày việc sử dụng hai công cụ biên tập chính:

  • FEM Editor: Phương pháp phần tử hữu hạn phân tích lưới nối đất phức tạp.
  • IEEE Editor: Tính toán lưới tiếp địa theo tiêu chuẩn IEEE Std 80. Quy trình thực hiện gồm: Khai báo điện trở suất của các lớp đất qua bảng Soil Editor, xác lập dòng ngắn mạch sự cố, bố trí số lượng thanh dẫn ngang (Conductors) và cọc tiếp địa thẳng đứng (Rods) qua Group Editor. Kết quả mô phỏng trả về giá trị điện trở nối đất $R_g$, điện áp tiếp xúc (Touch Voltage), điện áp bước (Step Voltage) và kiểm tra điều kiện an toàn cho người vận hành.

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

Kết quả chính

  • Mô phỏng phân bố công suất: Xác định chính xác phân bố dòng công suất và điện áp tại từng nút trên mô hình 7 bus; phát hiện hiện tượng sụt áp bất thường tại Bus 4 để đưa ra cảnh báo kỹ thuật.
  • Phân tích ngắn mạch theo IEC 60909: Xuất biểu đồ dạng sóng và giá trị đỉnh của các thành phần dòng $AC, DC$ và dòng ngắn mạch tổng, hỗ trợ lựa chọn công suất cắt cho thiết bị đóng cắt.
  • Ổn định quá độ với năng lượng tái tạo: Đánh giá được sự ảnh hưởng của nguồn điện mặt trời phân tán tại Bus 6 và Bus 8 đối với dao động góc rotor của máy phát truyền thống và mức độ cải thiện điện áp trên các thanh cái hệ thống IEEE 9-bus.
  • Tự động hóa chọn cáp và lưới nối đất: Đưa ra phương án chọn cáp cho hệ thống $22\text{ kV}/380\text{ V}/220\text{ V}$ đáp ứng đồng thời điều kiện dòng phát nóng và sụt áp; tối ưu hóa số lượng thanh dẫn và cọc tiếp địa trong thiết kế hệ thống tiếp địa an toàn.

Bảng tổng hợp phạm vi ứng dụng các mô-đun ETAP trong đồ án:

Mô-đun trong đồ án Mô hình áp dụng Tiêu chuẩn / Phương pháp Thông số đầu ra chính
Phân bố công suất Lưới 7 bus Newton-Raphson Điện áp nút ($U, \delta$), dòng công suất nhánh ($P, Q$)
Tính toán ngắn mạch Hệ thống công nghiệp IEC 60909 Dòng ngắn mạch đối xứng $I_k''$, dòng đỉnh $i_p$, thành phần $i_{dc}$
Ổn định quá độ IEEE 9-bus tích hợp PV Mô hình miền thời gian Góc rotor $\delta(t)$, đáp ứng điện áp $V(t)$, tổn thất $MW/M\text{var}$
Chọn tiết diện cáp Lưới $22\text{ kV}/0{,}38\text{ kV}/0{,}22\text{ kV}$ Tiêu chuẩn chế tạo cáp Tiết diện pha, dòng định mức hiệu chỉnh, độ sụt áp
Thiết kế nối đất Lưới Ground Grid IEEE Std 80 / FEM Điện trở tiếp địa $R_g$, điện áp bước, điện áp tiếp xúc

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

  • Đối với các nút có nguy cơ mất ổn định điện áp (như Bus 4 trong lưới 7 bus), cần phân bố lại công suất phụ tải hoặc nâng tiết diện đường dây liên kết.
  • Khi tích hợp các nguồn năng lượng tái tạo dung lượng lớn vào lưới, cần chạy mô phỏng phân tích ổn định quá độ trước khi vận hành thực tế để đánh giá quán tính hệ thống và dao động góc rotor của các máy phát lân cận.
  • Áp dụng các thư viện phần tử tiêu chuẩn trên phần mềm để hạn chế sai số chủ quan trong khâu tra cứu thông số vật lý của cáp và máy biến áp.

Đóng góp của đề tài

  • Hệ thống hóa quy trình tính toán 5 bài toán điện công nghiệp cơ bản trên phần mềm ETAP từ cơ sở toán học, thiết lập thông số đến phân tích kết quả.
  • Cung cấp tài liệu hướng dẫn từng bước (kèm hình ảnh giao diện và ví dụ số liệu thực tế) hỗ trợ việc học tập và thực hành các môn học chuyên ngành Hệ thống điện và Cung cấp điện.

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

  • Hạn chế: Đồ án chủ yếu tập trung vào một số ứng dụng cụ thể với các mô hình thử nghiệm quy mô vừa và nhỏ (mô hình 7 bus, IEEE 9-bus, lưới hạ áp mẫu); chưa triển khai mô phỏng trên toàn bộ sơ đồ mạng điện thực tế diện rộng; chưa kết nối đồng bộ giữa phân hệ quản lý thời gian thực (Real-time) với hệ thống phần cứng điều khiển.
  • Hướng nghiên cứu tiếp:
    • Khai thác các mô-đun chuyên sâu khác trong ETAP như: Phân tích sóng hài (Harmonic Analysis), Khởi động động cơ (Motor Acceleration Analysis), Phối hợp thiết bị bảo vệ (Star - Protective Device Coordination), Phân bố công suất tối ưu (Optimal Power Flow) và Đánh giá độ tin cậy lưới điện (Reliability Assessment).
    • Nghiên cứu bài toán bù tối ưu công suất phản kháng (Optimal Capacitor Placement) khi hệ thống tích hợp tỷ lệ năng lượng tái tạo cao.

Giá trị tham khảo

  • Đối tượng sử dụng: Sinh viên, giảng viên ngành Kỹ thuật Điện, Điện - Điện tử, Hệ thống điện và Cung cấp điện; kỹ sư thiết kế, tư vấn hệ thống điện công nghiệp và trạm biến áp.
  • Nội dung tham khảo giá trị:
    • Quy trình khai báo thông số thiết bị điện theo tiêu chuẩn quốc tế (IEC/ANSI) trong môi trường ETAP.
    • Các bước thiết lập kịch bản giải bài toán phân bố công suất bằng phương pháp Newton-Raphson.
    • Phương pháp mô hình hóa nguồn pin mặt trời (PV array + Inverter) trên lưới IEEE 9-bus để đánh giá ổn định hệ thống.
    • Trình tự thiết kế tiếp địa bằng công cụ FEM Editor và IEEE Editor trong mô-đun Ground Grid.

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

1. Trong bài toán phân bố công suất (Chương 2), tiêu chuẩn hội tụ của phương pháp Newton-Raphson được ETAP quy định ở mức nào?
Tiêu chuẩn hội tụ cho công suất thực và sai lệch công suất phản kháng trong phương pháp Newton-Raphson thường được thiết lập mặc định ở mức $0{,}001\text{ MW}$ và $0{,}001\text{ Mvar}$.

2. Tại sao ETAP lại thực hiện vài bước lặp Gauss-Seidel trước khi chạy giải thuật Newton-Raphson?
Do phương pháp Newton-Raphson rất nhạy cảm và phụ thuộc nhiều vào giá trị điện áp ban đầu của các bus, nên ETAP tạo ra vài lần lặp Gauss-Seidel sơ bộ để thiết lập tập giá trị ban đầu hợp lý cho điện áp nút, giúp quá trình hội tụ diễn ra nhanh và tránh phân kỳ.

3. Trong kết quả tính toán phân bố công suất cho mạng điện 7 bus ở Chương 2, thanh cái nào gặp sự cố bất thường?
Bus 4 ghi nhận tình trạng điện áp không ổn định (sụt áp ngoài ngưỡng cho phép) do mang phụ tải lớn ($195\text{ MW} + j50\text{ Mvar}$), đòi hỏi phải điều chỉnh lại phân bố phụ tải hoặc thay đổi tiết diện dây dẫn.

4. Khi mô phỏng hệ thống năng lượng tái tạo ở Chương 4, các tổ hợp pin mặt trời (PV array) được ghép nối tại những vị trí nào?
Các tổ hợp pin quang điện (PV array) kết hợp bộ nghịch lưu (Inverter) và máy biến áp được tích hợp tại Bus 6 và Bus 8 của mô hình chuẩn IEEE 9-bus.

5. Mô-đun Ground Grid của ETAP (Chương 6) hỗ trợ những công cụ biên tập nào để tính toán lưới tiếp địa?
Mô-đun hỗ trợ hai công cụ biên tập chính là IEEE Editor (tính toán dựa trên tiêu chuẩn IEEE Std 80) và FEM Editor (tính toán theo phương pháp phần tử hữu hạn).

Kết luận

Đồ án tốt nghiệp của tác giả Ngô Duy Hiếu đã hệ thống hóa và làm chủ quy trình ứng dụng phần mềm ETAP để giải quyết các bài toán trọng tâm trong thiết kế và phân tích hệ thống điện. Thông qua các ví dụ mô phỏng cụ thể từ phân bố công suất, ngắn mạch theo chuẩn IEC 60909, định cỡ cáp, phân tích ổn định hệ thống tích hợp năng lượng tái tạo đến thiết kế lưới nối đất, đề tài đã chứng minh tính hiệu quả và độ tin cậy của công cụ mô phỏng số. Công trình là tài liệu tham khảo kỹ thuật hữu ích, có tính ứng dụng cao cho sinh viên và kỹ sư ngành Điện - Điện tử trong công tác học tập, nghiên cứu và thiết kế thực tế.