Nghiên ứu áp dụng mô hình đẳng trị một máy phát trong đánh giá ổn định hệ thống điện

BỘ GIÁO DỤC VÀ ĐÀO TẠO
TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI
---------------------------------------

NGUYỄN VĂN CHIẾN

NGHIÊN CỨU ÁP DỤNG MƠ HÌNH ĐẲNG TRỊ MỘT MÁY
PHÁT TRONG ĐÁNH GIÁ ỔN ĐỊNH HỆ THỐNG ĐIỆN
Chuyên ngành mạng và hệ thống điện

LUẬN VĂN THẠC SĨ KHOA HỌC

NGƯỜI HƯỚNG DẪN KHOA HỌC :
TSKH Trần Kỳ Phúc

Hà Nội – 2011

Tai ngay!!! Ban co the xoa dong chu nay!!! 17057205300851000000


MỞ ĐẦU ........................................................................................................................................... 4

T
03

T
03

1. ĐẶC ĐIỂM CỦA HỆ THỐNG ĐIỆN HIỆN ĐẠI VÀ CÁC YÊU CẦU LIÊN QUAN ĐẾN
ĐÁNH GIÁ ỔN ĐỊNH HỆ THỐNG ĐIỆN. ................................................................................. 4

T
03

T
03

2. SỰ PHÁT TRIỂN CỦA HỆ THỐNG ĐIỆN VIÊT NAM VÀ CÁC VẤN ĐỀ TÍNH TỐN
PHÂN TÍCH ỔN ĐỊNH ĐỘNG. ................................................................................................... 4

T
03

T
03

3. LÝ DO CHỌN ĐỀ TÀI VÀ MỤC ĐÍCH NGHIÊN CỨU CỦA LUẬN ÁN ............................ 6

T
03

T
03

4. PHƯƠNG PHÁP NGHIÊN CỨU ............................................................................................. 7

T
03

T
03

5. LĨNH VỰC ÁP DỤNG .............................................................................................................. 7

T
03

T
03

CHƯƠNG 1 TỔNG QUAN VỀ ỔN ĐỊNH ĐỘNG VÀ CÁC PHƯƠNG PHÁP TÍNH TỐN ỔN
ĐỊNH ĐỘNG CỦA HỆ THỐNG ĐIỆN ........................................................................................... 8

T
03

T
03

T
03

T
03

1.1 ỔN ĐỊNH ĐỘNG TRONG HỆ THỐNG ĐIỆN ..................................................................... 8

T
03

T

03

1.2. C ÁC THÀNH PHẦN CƠ BẢN CUẢ BÀI TỐN ỔN ĐỊNH ĐỘNG ............................ 11

T
03

T
03

T
03

T
03

1.2.1.

T
03

Phương trình vi phân mô tả chuyển động quá độ roto của các máy phát điện đồng bộ. 11

T
03

T
03

T

03

1.2.2.

Hệ phương trình vi phân QTQĐ điện từ trong máy phát điện...................................... 14

1.2.3.

Hệ phương trình cân bằng cơng suất lưới................................................................... 16

T
03

T
03

T
03

T
03

T
03

T
03

T
03

T
03

1.3. CÁC PHƯƠNG PHÁP PHÂN TÍCH TÍNH TỐN ỔN ĐỊNH ĐỘNG. ............................. 18

T
03

T
03

T
03

T
03

1.3.1.Phương pháp Ơle.............................................................................................................. 18

T
03

T
03

1.3.2.Phương pháp Ơle biến thể ................................................................................................ 19

T
03

T
03

1.3.3. Phương pháp Runge - Kutta (R - K) ................................................................................ 20

T
03

T
03

4. VẤN ĐỀ NÂNG CAO ỔN ĐỊNH HTĐ .................................................................................. 21

T
03

T
03

T
03

T
03

KẾT LUẬN CHƯƠNG 1 ............................................................................................................ 22

T
03

T
03

CHƯƠNG II XÂY DỰNG CẤU HÌNH CHUNG HỆ THỐNG ĐÁNH GIÁỔN ĐỊNH ĐỘNG HỆ
THỐNG ĐIỆN ................................................................................................................................ 23

T
03

T
03

T
03

2.1. TỔNG QUAN ....................................................................................................................... 23

T
03

T
03

2.1.1. Đặt vấn đề ...................................................................................................................... 23

T
03

T

03

2.1.2. Các bước xây dựng hệ thống DSA .................................................................................... 24

T
03

T
03

1


2.1.3. Cấu trúc bài toán đánh giá ổn định HTĐ ......................................................................... 25

T
03

T
03

2.1.4. Cơng cụ tính tốn đánh giá ổn định HTĐ ......................................................................... 26

T
03

T
03

2.1.5. Hạ tầng thông tin - đo lường phục vụ đánh giá ổn định HTĐ ........................................... 27

T
03

T
03

2.2. XÂY DỰNG CẤU HÌNH CHUNG HỆ THỐNG ĐÁNH GIÁ GIÁM SÁT ỔN ĐỊNH
ĐỘNG DSA TRONG HỆ THỐNG ĐIỆN VIỆT NAM .............................................................. 28

T
03

T
03

2.2.1. Nhiệm vụ và phương pháp giải quyết ............................................................................... 28

T
03

T
03

2.2.1.1. Nhiệm vụ và tính năng yêu cầu hệ thống DSA........................................................... 28

T
03

T

03

2.2.1.2. Phương pháp giải quyết bài toán của hệ thống DSA .................................................. 29

T
03

T
03

2.2.2. Thiết kế định hướng hệ thống DSA trong hệ thống điện Việt Nam..................................... 30

T
03

T
03

2.2.3. Cấu hình phần cứng hệ thống đánh giá ổn định động ....................................................... 31

T
03

T
03

2.2.4. Cấu hình phần mềm ......................................................................................................... 32

T
03

T
03

2.2.4.1. Mơ tả bản chất của phần mềm DSA........................................................................... 32

T
03

T
03

2.2.4.2. Các thành phần chính của hệ thống DSA................................................................... 33

T
03

T
03

2.2.4.3. Cấu hình phần mềm hệ thống DSA ............................................................................ 34

T
03

T
03

2.2.4.4. Thu thập dữ liệu ........................................................................................................ 35

T
03

T
03

2.2.4.5. Ví dụ biểu đồ thuật toán vận hành online hệ thống DSA ............................................ 37

T
03

T
03

2.2.4.6. Vận hành off-line hệ thống DSA ................................................................................ 38

T
03

T
03

2.2.4.7 Các phần mềm chức năng ứng dụng hệ thống DSA..................................................... 39

T
03

T
03

KẾT LUẬN CHƯƠNG 2 ............................................................................................................ 41

T
03

T
03

CHƯƠNG IIIÁP DỤNG MÔ HÌNH ĐẲNG TRỊ MỘT MÁY PHÁT TRONG ĐÁNH GIÁ ỔN
ĐỊNH QUÁ ĐỘ HỆ THỐNG ĐIỆN-PHƯƠNG PHÁP SIME ...................................................... 42

T
03

3.1. PHƯƠNG PHÁP SIME ....................................................................................................... 42

T
03

T
03

3.1.1. Phân loại các tổ máy thành 2 nhóm Critical và non-Critical ............................................ 43

T
03

T
03

3.1.2. Xác định các thơng số của 2 nhóm máy phát. ................................................................... 45

T
03

T
03

3.1.3. Xác định các thông số của máy phát đẳng trị tương đương. ............................................. 47

T
03

T
03

3.1.4. Đánh giá khả năng ổn định của hệ thống dựa vào đặc tính của OMIB.............................. 49

T
03

T
03

3.1.4.1. Điều kiện không ổn định............................................................................................ 50

T
03

T

03

2


3.1.4.2 Xác định các thông số:

T
03

T
03

t u , ωu , δu ........................................................................... 52

3.1.4.3. Điều kiện ổn định ...................................................................................................... 54

T
03

T
03

3.1.4.4. Xác định các thông số: δ r , t r .................................................................................... 56

T
03

T
03

3.2. THUẬT TOÁN ĐÁNH GIÁ ỔN HỆ THỐNG ĐIỆN THEO PHƯƠNG PHÁP SIME…...57
3.3. ÁP DỤNG THUẬT TOÁN XÂY DỰNG CHƯƠNG TRÌNH ĐÁNH GIÁ PHÂN TÍCH ỔN
ĐỊNH ĐỘNG HỆ THỐNG ĐIỆN............................................................................................... 58

T
03

T
03

3.3.1. Sự cố ngừng hoạt động đường dây 100-200. .................................................................... 60

T
03

T
03

3.3.1.1. thời gian giải trừ sự cố là CT=0.5s ........................................................................... 60

T
03

T
03

3.3.1.2. thời gian giải trừ sự cố là CT=0.7s ........................................................................... 62

T

03

T
03

3.3.2. Sự cố ngắn mạch tại bus 100 đồng thời một máy phát ngừng hoạt động. .......................... 64

T
03

T
03

3.3.2.1 thời gian giải trừ sự cố là CT=0.5s ............................................................................ 64

T
03

T
03

3.3.2.2. thời gian giải trừ sự cố là CT=0.3s ........................................................................... 66

T
03

T
03

KẾT LUẬN CHƯƠNG 3 ............................................................................................................ 69

T
03

T
03

CHƯƠNG 4 ÁP DỤNG PHẦN MỀM XÂY DỰNG TỪ THUẬT TOÁN PHƯƠNG PHÁP SIME
ĐÁNH GIÁ ỔN ĐỊNH HỆ THỐNG ĐIỆN VIỆT NAM NĂM 2015. ............................................ 70

T
03

4.1. SỰ CỐ NGẮN MẠCH TRÊN ĐƯỜNG DÂY 500kV HỒ BÌNH-SƠN LA. ..................... 71

T
03

T
03

T
03

T
03

4.1.1. phương pháp kinh điển..................................................................................................... 71

T

03

T
03

T
03

T
03

4.1.2. phương pháp SIME .......................................................................................................... 71

T
03

T
03

T
03

T
03

4.2. SỰ CỐ NGẮN MẠCH TRÊN THANH CÁI 500kV CỦA THUỶ ĐIỆN SƠN LA. ............ 74

T
03

T
03

T
03

T
03

4.2.1. phương pháp kinh điển..................................................................................................... 74

T
03

T
03

T
03

T
03

4.2.2. phương pháp SIME .......................................................................................................... 75

T
03

T
03

T
03

T
03

KẾT LUẬN CHƯƠNG 4 ............................................................................................................ 77

T
03

T
03

3


MỞ ĐẦU
1. ĐẶC ĐIỂM CỦA HỆ THỐNG ĐIỆN HIỆN ĐẠI VÀ CÁC YÊU CẦU LIÊN
QUAN ĐẾN ĐÁNH GIÁ ỔN ĐỊNH HỆ THỐNG ĐIỆN.
Cùng với sự phát triển của đời sống xã hội hệ thống điện của mỗi quốc gia có
quy mô công suất ngày càng lớn, cấu trúc lưới ngày càng phức tạp. Hầu như khơng cịn
tồn tại các hệ thống điện khu vực hoạt động độc lập với một vài nhà máy điện. Đặc
điểm này đã làm thay đổi cơ bản về chất của các hệ thống điện hiện đại: mang các đặc
trưng của hệ thống lớn phức tạp [23]. Hệ thống điện Việt Nam sau khi được hợp nhất
bằng đường dây siêu cao áp 500 KV Bắc – Trung – Nam cũng đã trở thành HTĐ hợp
nhất liên kết một số lượng lớn nguồn và phụ tải ở cả 3 miền. Nhiều bài toán thiết kế và
vận hành HTĐ trước đây nay đòi hỏi phải được xem xét lại với sự thay đổi cơ bản cả
về mơ hình cũng như nội dung các phương pháp giải. Chẳng hạn tính đa dạng chủng

loại của các nguồn (thủy điện; nhiệt điện than; dầu; nhiệt điện khí chu trình hỗn hợp…)
trong HTĐ địi hỏi giải bài tốn vận hành kinh tế hệ thống theo mơ hình tổng qt. Các
đường dây dài điện áp siêu cao bị giới hạn công suất truyển tải theo điều kiện ổn định
tĩnh cũng làm phức tạp thêm đáng kể về thuật toán giải (cần xét đến các ràng buộc về
giới hạn công suất trao đổi giữa các khu vực) … Đặc biệt, các bài toán đánh giá về độ
tin cậy, phân tích tính ổn định hệ thống sẽ cần phải có phương pháp phân tích đánh giá
thích hợp (nhanh-chính xác) trong các tài liệu trước đây [24].Trong khi đó các u cầu
đảm bảo về tính ổn định, độ tin cậy làm việc hệ thống lại đòi hỏi ngày càng cao hơn.
Những thay đổi này làm phức tạp thêm nhiều cho các bài tốn phân tích và điều khiển
hệ thống điện.
2. SỰ PHÁT TRIỂN CỦA HỆ THỐNG ĐIỆN VIÊT NAM VÀ CÁC VẤN ĐỀ
TÍNH TỐN PHÂN TÍCH ỔN ĐỊNH ĐỘNG.
Hệ thống 500kV gồm đường dây truyền tải dài gần 1500 km và các trạm biến áp
Hịa Bình, Hà Tĩnh, Đà Nẵng, Pleiku, Phú Lâm đã liên kết hệ thống điện Việt Nam
thành một hệ thống hợp nhất từ năm 1994. Hệ thống 500kV đã làm thay đổi đáng kể về

4


mặt cấu trúc hệ thống điện Việt Nam [1], [3]. Bên cạnh những ưu việt của hệ thống lớn
như tạo điều kiện hỗ trợ công suất tác dụng giữa các khu vực, nâng cao hiệu quả vận
hành kinh tế đảm bảo độ tin cậy về cung cấp điện và chất lượng điện năng đến hộ tiêu
thụ; hệ thống điện Việt Nam ngày càng trở nên hết sức phức tạp về mặt cấu trúc, đa
dạng về công nghệ phát điện và quy mô công suất nguồn, đặc biệt bao gồm những tổ
máy có cơng suất lớn. Ngồi ra tính phức tạp của sơ đồ lưới xác suất xảy ra sự cố lớn
từ các đường dây dài 500kV đã ảnh hưởng nhiều đến tính chất chung và ổn định của hệ
thống điện.
Nhận rõ vai trò quan trọng của vấn đề này ngay từ giai đoạn thiết kế đường dây
siêu cao áp 500kV nhiều cơ quan thiết kế và tư vấn trong và ngồi nước đã tiến hành
tính tốn phân tích các đặc trưng ổn định của hệ thống điện hợp nhất Việt Nam. Ngồi

các tính tốn của các chun gia Việt Nam, các nghiên cứu của tư vấn Nippon Koei
(Nhật Bản), PPI/SECVI (Úc) về ổn định quá độ của HTĐ Việt Nam sau khi liên kết
đường dây 500kV đã đưa ra được những kết quả quan trọng, giúp lựa chọn hợp lý các
phương án kết dây và lắp đặt hệ thống các thiết bị bù. Tuy nhiên các tính tốn đều được
thực hiện dựa trên các mơ hình có những mức độ đơn giản hóa hệ thống nhất định.
Hơn nữa, từ sau khi đường dây được đưa vào vận hành đến nay cấu trúc hệ
thống điện Việt Nam đã thay đổi nhiều cả về quy mô công suất nguồn cũng như sơ đồ
lưới điện, ảnh hưởng nhiều đến tính chất ổn định của hệ thống điện. QTQĐ sau sự cố
của hệ thống nhiều tổ máy lớn như nước ta hiện nay diễn ra hết sức phức tạp. Trong
thực tế vận hành, đã có hiện tượng của nhiều sự cố mà sau đó dẫn đến hệ thống bị tan
rã, nguyên nhân chắc chắn có liên quan nhiều đến vấn đề mất ổn định hệ thống.
Mặt khác theo tổng sơ đồ phát triển của HTĐ Việt Nam (đến năm 2015) lưới
500kV của Việt Nam sẽ tăng trưởng trên một phạm vi rất rộng: Từ nhà máy thủy điện
Sơn La đến Cà Mau với tổng chiều dài lên tới trên 5000km. Nhiều nhà máy điện lớn sẽ
đưa vào vận hành như Sơn La, Nhơn Trạch, Ô Môn…

5


Những đặc điểm vừa nêu đã thể hiện một yêu cầu cấp thiết phải có những
nghiên cứu cụ thể hơn về các vấn đề ổn định hệ thống, đặc biệt những nghiên cứu để
đưa ra phương pháp phân tích đánh giá các đặc trưng động của hệ thống điện phức tạp
tương ứng với sơ đồ HTĐ của Việt Nam phát triển trong tương lai. Có thể nói, cho đến
nay những tính tốn nghiên cứu về những nội dung này cịn rất hạn chế, thiếu những
nghiên cứu sâu về ổn định động hệ thống tương ứng với các thông số và sơ đồ đặc
trưng của HTĐ Việt Nam.
3. LÝ DO CHỌN ĐỀ TÀI VÀ MỤC ĐÍCH NGHIÊN CỨU CỦA LUẬN ÁN
Từ hiện trạng và nhiều vấn đề của hệ thống điện Việt Nam cần được quan tâm
như đã phân tích ở trên có thể thấy vấn đề nghiên cứu ổn định động và nâng cao ổn
định là nội dung có ý nghĩa hết sức quan trọng trong quá trình phát triển hệ thống điện

hiện đại (tốc độ tínhtốn- trực tuyến-nhanh- khả năng tin cậy cao)
Đề tài nghiên cứu của luận án: “Áp dụng mơ hình đẳng trị một máy phát trong
đánh giá ổn định hệ thống điện” được lựa chọn xuất phát từ mong muốn được đóng
góp một phần nhỏ bé trong hướng nghiên cứu nói trên.
Mục đích chính của luận án là:
-

Dựa trên lý thuyết hiện đại về tính tốn phân tích ổn định động của HTĐ phức
tạp, nghiên cứu đặc điểm cụ thể các trang thiết bị và sơ đồ đầy đủ của HTĐ Việt
Nam xây dựng thuật đánh giá ổn định quá độ hệ thống điện khi có sự cố theo
phương pháp SIME.

-

Áp dụng thuật toán trên với các sơ đồ cụ thể của HTĐ Việt Nam các giai đoạn,
đánh giá mức độ ổn định của hệ thống tương ứng với những tình huống sự cố
điển hình phổ biến diễn ra trong hệ thống.

6


4. PHƯƠNG PHÁP NGHIÊN CỨU
Đề tài đi theo hướng nghiên cứu là:
-

Ứng dụng các phương tiện phần mềm kết hợp máy tính hiệu năng cao để xác
định các thơng số của bài tốn phân tích ổn định động của thống điện, bước đầu
đánh giá phân tích một số đặc trưng động của hệ thống điện Việt Nam. Dự kiến
sử dụng kết hợp chương trình tính tốn phân tích hệ thống PSS/E-29.

-

Dựa trên các tài liệu tham khảo trong nước và nước ngồi về tính tốn, phân tích
ổn định động của hệ thống điện, từng bước xây dựng thuật đánh giá ổn định
theo phương pháp mơ hình đẳng trị một máy phát. Sau đó, trên cơ sở thuật tốn
xây dựng phần mềm để ứng dụng trong thực tế.

-

Áp dụng phần mềm này với mơ hình hệ thống điện đơn giản, từ đó hồn thiện
phần mềm để có thể áp dụng cho các sơ đồ cụ thể của toàn bộ hệ thống điện
Việt Nam 2015.

5. LĨNH VỰC ÁP DỤNG
Đề tài nghiên cứu xuất phát từ nhu cầu thực tế của HTĐ Việt Nam, do đó
các kết quả có ý nghĩa thực tiễn, có thể áp dụng cho các giai đoạn quy hoạch, thiết
kế và vận hành hệ thống.
Nơi ứng dụng: Các viện nghiên cứu quy hoạch phát triển HTĐ, Trung tâm
Điều độ HTĐ quốc gia, Các công ty tư vấn thiết kế xây dựng điện.

7


CHƯƠNG 1
TỔNG QUAN VỀ ỔN ĐỊNH ĐỘNG VÀ CÁC PHƯƠNG PHÁP TÍNH
TỐN ỔN ĐỊNH ĐỘNG CỦA HỆ THỐNG ĐIỆN
1.1.

ỔN ĐỊNH ĐỘNG TRONG HỆ THỐNG ĐIỆN

Ổn định động hay ổn định quá độ (Transient Stability) của hệ thống điện (HTĐ)
là ổn định của hệ thống sau những kích động lớn. Các kích động này là các sự cố ngẫu
nhiên như sự cố ngắn mạch cắt các đường dây liên kết, cắt đột ngột các máy biến áp,
các máy phát công suất lớn… và các tác động do thao tác vận hành làm mất cân bằng
cơng suất trong hệ thống. Khi đó trạng thái cân bằng mômen quay trong máy phát bị
phá vỡ, xuất hiện gia tốc làm thay đổi mạnh góc lệch roto của máy phát. Quá trình quá
độ (QTQĐ) của hệ thống đáp ứng và điều chỉnh lại sự thay đổi trạng thái cân bằng này
(hay hấp thụ động năng hệ thống [10]) diễn ra có thể ổn định hoặc khơng ổn định tùy
thuộc mức độ của các kích động: thời gian tồn tại sự cố và các tác động giải trừ sự cố
(ngắn mạch) [23].
Điều kiện để hệ thống có ổn định động đó là:
-

Tồn tại điểm cân bằng ổn định tĩnh sau sự cố (ứng với CĐXL sau sự cố)

-

Thông số biến thiên hữu hạn trong QTQĐ và dao động tắt dần về thông số xác
lập mới
Như vậy sự tồn tại của CĐXL sau sự cố là điều kiện cần nhưng chưa đủ cho tính

ổn định động của hệ thống. Chính điều kiện 2 dẫn đến yêu cầu phải áp dụng các
phương pháp phân tích ổn định động của hệ thống điện.

8


δ

δ

δ
δ
δ
δ

δ

δ

δ

a)

b)
Hình 1.1. Đánh giá ổn định động của hệ thống điện

Đối với HTĐ đơn giản (máy phát nối với thanh cái điện áp khơng đổi) ổn định
động được giải thích bằng khả năng cân bằng được động năng do kích động sinh ra với
công hãm xuất hiện trong phản ứng hệ thống (hình 1.1.a). Khi đó góc lệch dao động
khơng vượt quá 180 [2]. Khi mất ổn định động góc dao động tăng trưởng vô hạn ngay
sau chu kỳ đầu (hình 1.1b).
Đối với HTĐ nhiều máy phát, dao động roto của các máy phát được thể hiện là
sự biến thiên góc lệch δ i (t) theo thời gian so với trục quay đồng bộ máy phát (hình
1.2a). Tuy nhiên để đánh giá ổn định động trong hệ thống nhiều máy phát, cần phải xác
định biến thiên góc lệch tương đối giữa các máy phát với nhau hoặc so với góc lệch
trung bình của hệ thống [2] ( hình 1.2b). Hệ thống ổn định khi tất cả các góc lệch tương
đối dao động trong phạm vi hữu hạn. Về lý thuyết, đối với HTĐ phức tạp nếu chỉ tính
tốn được góc lệch tương đối trong một khoảng thời gian xác định thì khơng thể kết
luận được về tính hữu hạn của phạm vi dao động góc lệch của cả QTQĐ. Tuy nhiên,
các nghiên cứu ứng dụng cho thấy HTĐ sẽ ổn định động nếu dao động hữu hạn trong

9


phạm vi 10s đầu [8]. Ngồi ra tính ổn định cịn được phán đốn theo cả xu hướng biến
thiên (dao động tắt dần hay mạnh dần theo thời gian).

δ

δ
δ31
δ1

δ2

δ3

δ12
δ23

t

t

a)

b)

Hình 1.2 . Dao động góc lệch roto (a) tuyệt đối, (b) tương đối.
Nghiên cứu tính tốn phân tích ổn định hệ thống điện, đặc biệt quan tâm đến các

kích động làm hệ thống mất ổn định. Ngoài ra, diễn biến QTQĐ của hệ thống điện
nhiều máy khi mất ổn định rất phức tạp, thể hiện qua nhiều dạng khác nhau. Nghiên
cứu đặc điểm mất ổn định của hệ thống phức tạp sẽ cho ta rút ra được tính chất chung
và khả năng duy trì đồng bộ tương đối của các máy phát, từ đó có những tác động điều
khiển hiệu quả tương thích, nhằm duy trì hệ thống tái đồng bộ.
Ổn định góc lệch các máy phát khơng những chịu ảnh hưởng bởi cấu trúc các
máy phát điện đồng bộ mà cịn ảnh hưởng bởi cấu trúc lưới điện. Thơng số các đường
dây liên kết, trị số công suất truyền tải, đặc biệt đối với các đường dây dài, ảnh hưởng
lớn đến diễn biến QTQĐ khi sự cố. Ảnh hưởng của sự cố ngắn mạch trên các đường
dây truyền tải dài đến độ ổn định nhiều lúc mạnh hơn cả ngắn mạch tại máy phát [10].
Do vậy yêu cầu về ổn định của hệ thống có cấu trúc phức tạp cao hơn và bài tốn ổn
định có vai trị quan trọng hơn.

10


Để đưa ra phương pháp đánh giá ổn định động của hệ thống điện nhanh và chính
xác, thì ta cần chú ý đến các phương trình và hệ phương trình cơ bản sau [2]:
- Hệ phương trình cân bằng cơng suất lưới
-

Phương trình chuyển động quay roto máy phát

-

Hệ phương trình QTQĐ điện từ trong các cuộn dây máy phát

Khi trong hệ thống có sự cố các thơng số trong phương trình và hệ phương trình sẽ
thay đổi, căn cứ vào sự thay đổi của các thơng số đó ta có thể đánh giá được hệ thống
có ổn định hay khơng.

1.2.

CÁC THÀNH PHẦN CƠ BẢN CUẢ BÀI TỐN ỔN ĐỊNH ĐỘNG

1.2.1. Phương trình vi phân mơ tả chuyển động q độ roto của các máy phát
điện đồng bộ.
Phương trình chuyển động quay roto của máy phát điện đồng bộ được thiết lập
trên cơ sở định luật 2 Newton (viết cho vật thể quay):

d 2y
J 2 =M T - M E
dt
R

R

R

(1.1)

R

Trong đó:
MT , ME là các mơ men quay của tuabin và máy phát
R

R

R

R

J là mơmen qn tính của roto tuabin - máy phát, phụ thuộc vào cấu tạo và khối
lượng phần quay.
γ là góc quay tuyệt đối của roto máy phát, là góc giữa trục dọc gắn liền với roto và
một trục đứng yên (gắn với stato).

11


ω

δn
δ

δ12

ω0
γ

δ2

δ1

a)

b)
Hình 1.4. Chuyển động roto của máy phát

Khi nghiên cứu QTQĐ trong HTĐ phương trình thường được biến đổi về dạng

thuận tiện hơn qua một số phép biến đổi sau:
a. Tính với góc lệch tương đối.
Góc quay tuyệt đối γ thường biến thiên rất nhanh theo thời gian (ngay cả trong
chế độ xác lập) khơng tiện cho phân tích các diễn biến quá độ. Người ta thường biến
đổi phương trình viết theo góc lệch tương đối δ giữa trục dọc roto với một trục quay
quy ước, cùng chiều với roto theo tốc độ đồng bộ Ω0 .
R

R

d 2 y d 2δ
=
Theo hình 1.4a: γ(t) = Ω0. t + δ(t) nên
dt 2 dt 2
R

R

b. Sử dụng hằng số quán tính H ( hay T J ) thay cho mô men quán tính J.
R

R

Hằng số qn tính H có trị số bằng động năng tương đối định mức của roto khi
quay với tốc độ đồng bộ:

JΩ 20
H=
2.SCB

, MWS / MVA. (hoặc được viết bằng sec)
R

R

12


Trong đó: S CB là lượng cơ bản chọn cho cơng suất. Ω 0 có thể tính theo số vịng
R

R

R

R

quay định mức n (vòng /phút) của máy phát:

2πn
,1/ sec.
60

Ω0 =
R

R

Hằng số quán tính TJ = 2H, thường được tính với đơn vị sec nhưng thực chất
R

R

cũng đã tính trong hệ đơn vị tương đối của S CB ( MVA).
R

.

R

c. Tính gần đúng mômen theo công suất.
MT = Ω.P T ≈ Ω.P T ;
R

R

R

R

R

M E = Ω.P E ≈ Ω0 .P E ;

R

R

R

R

R

R

R

R

R

Thực tế, biến thiên của góc lệch tương đối dδ/dt nhỏ hơn rất nhiều so với biến
thiên góc lệch tuyệt đối, nên có thể coi Ω = Ω0 + dδ/dt ≈ Ω 0.
R

R

R

R

d. Tính theo góc độ điện.
Với máy phát có p đơi cực, khi roto quay 1 vịng các đại lượng điện đã biến
thiên p chu kỳ, nghĩa là tần số điện ω = pΩ.
Với các biến đổi trên, phương trình (1.1) có thể viết lại ở dạng sau:

Tj
d δ
SCB 2 = P T – P E

ϖ0
dt
2

(1.2a)

R

R

R

R

R

R

Trong đó: P T - công suất cơ và PE – công suất điện, tính bằng MW.
R

R

R

R

Khi tính P T và PE trong hệ tương đối cơ bản S CB (cùng với T J ) ta có:
R

R

R

R

R

R

R

Tj d 2 δ
= PT - PE
ϖ0 dt 2
R

R

R

R

(1.2b)

R

Nếu xét đến các tổn hao (ma sát chuyển động và tổn hao trong các cuộn dây
cản) phương trình cịn có thêm thành phần tỉ lệ với tốc độ biến thiên góc lệch δ:

2

Tj d δ
dδ
= P T - PE (1.3)
2 + KD .
dt
ϖ 0 dt
R

R

13

R

R

R

R


Dạng (1.3) được sử dụng phổ biến trong các chương trình tính tốn QTQĐ
điện cơ của HTĐ.
1.2.2. Hệ phương trình vi phân QTQĐ điện từ trong máy phát điện.
Hệ phương trình QTQĐ điện từ trong các cuộn dây (roto và stato) máy phát
được mô tả trong hệ tọa độ quay do Park và Gorev đề xuất như trên hình 1.5.

ud

A

id

iq q

iQ

uq

uf
iD
C

B

Hình 1.5. Biểu diễn các thơng số máy phát trong hệ tọa độ vng góc
Các cuộn dây stato được quy đổi về roto gồm hai cuộn dây vng góc nhau
quay cùng tốc độ với roto. Dòng điện trong mỗi cuộn dây tương ứng là id , iq. Do khơng
R

R

R

R

có chuyển động quay tương đối nên các hệ số tự cảm và hỗ cảm là hằng số. Hệ phương
trình Gorev - Park (G-P) mô tả QTQĐ trong các cuộn dây của máy phát đồng bộ như

sau:

14


d ψd

U
Ri
=
−
−
− ωψq
d
d

dt

U = −Ri − d ψq + ωψ
q
d
 p
dt

U = −Ri − dψ 0
0
 0
dt

U = −Ri − dψ f

f
 f
dt

dψ D
0 = R Di D +
dt


dψ
0 = R Q iQ + Q

dt

(1.4a)

3 '

ψ
=
+
+
(L
M
L)id + Mf if + MD iD
d

2

'

ψ = (L + M − 3 L)i
q
q + Mq iq

2
ψ = (L − 2M)i
0
 0

3
ψ f = Lf if + 2 Mf id + MfD iD

ψ = L i + 3 M i + M i
D D
D d
fD f
 D
2

ψ Q = LQ iQ + 3 MQ iq

2

(1.4b)

Trong đó: Ψd , Ψ q ,Ψ0 , Ψ f , ΨD , ΨQ, id , i q , i 0 , if , i D , iQ là các từ thơng móc
R

R

R

R

R

R

R

R

R

R

R

R

R

R

R

R

R

R

R

R

R

R

R

R

vịng và dịng điện trong các cuộn dây stato, cuộn dây kích từ, cuộn cản dọc trục và
ngang trục.
L, L’, M, M’ là các hệ số tự cảm, hỗ cảm của các cuộn dây stato. Các hệ số có
thể xem là hằng số nếu bỏ qua bão hòa từ.

15


Lf là hệ số tự cảm cuộn dây kích từ, LD , LQ tương ứng là các hệ số tự cảm của
R

R

R

R

R

R

cuộn cản dọc trục và ngang trục.
Mf , MD , MQ là hệ số hỗ cảm của cuộn dây kích từ, cuộn cản dọc trục và ngang
R

R

R

R

R

R

trục đối với cuộn dây stato. M fD là hỗ cảm của cuộn dây kích từ và cuộn cản dọc trục.
R

R

Rf , R D, R Q, R là điện trở các cuộn dây trên roto và dây quấn stato.
R

R

R

R

R

R

Uf: điện áp kích từ. U d , Uq : các thành phần điện áp dọc trục và ngang trục stato
R

R

R

R

R

R

máy phát. U: điện áp đầu cực máy phát
Hệ phương trình QTQĐ (G-P) trên mơ tả chi tiết trạng thái điện từ của máy
phát điện quay. Quan hệ của các phương trình QTQĐ điện từ với hệ phương trình
chuyển động tương đối roto các máy phát là biểu thức:
ME = Ψq id – Ψ d iq (hệ đơn vị tương đối)
R

R

R

R

R

R

R

R

R

R

(1.5)

Phương trình chứa các đại lượng đầu cực của các máy phát Ud , Uq , i d , iq , có
R

R

R

R

R

R

R

R

quan hệ chặt chẽ với hệ phương trình lưới.
1.2.3. Hệ phương trình cân bằng công suất lưới
Trạng thái của lưới điện được mô tả dưới dạng hệ phương trình cân bằng cơng
suất nút. Giả thiết biết rõ tổng dẫn riêng và tương hỗ các nút, hệ phương trình sẽ viết
được như sau ( cho mỗi thời điểm của QTQĐ):
n
 2
α
+
U
y
sin
U i U j sin( θi − θj − αij ) = Pi
∑
ii
 i ii
=
j 1

j≠ i

n
 2
U i yii cos αii − ∑ U iU jcos( θi − θ j − αij ) = Q i
j=1


j≠ i

i = 1,2,3...n



(1.6)

Trong đó: U i ,θ i – là mơdun và góc pha của điện áp nút i;
R

Ui =
R

R

R

R

R

U 2d + U 2q ;

θi = arctg(
R

R

Ud

)
Uq

(1.7)

16


Yii = y ii ∠ αii - 90°: tổng dẫn riêng của nút i ;
R

R

R

R

R

R

Yii = y ij ∠ αij - 90°: tổng dẫn tương hỗ giữa nút i và nút j.
R

R

R

R

R

R

Với nút đầu cực máy phát:

E Qi Ui

Pi = X sin( δi − θi ) = PEi
q


2
Q = EQi − EQi Ui cos( δ − θ ) = Q
i
i
Ei
 i Xq
Xq

Eqi = E Qi
R

R

R

(1.8)

X − X qi

Xdi − X 'di
− E 'qi di
Xqi − X 'di
Xqi − Xdi

(1.9)

R

Trong đó: E qi, E’ qi , EQi - tương ứng là các sđđ đồng bộ, sđđ quá độ và giả
R

R

R

R

R

R

tưởng của máy phát thứ i [2]. Ký hiệu θ tương ứng với góc lệch pha giữa sđđ Eq của
R

R

máy phát và điện áp đầu cực. Sức điện động E Q chỉ có ý nghĩa tính tốn, sđđ Eq tỉ lệ
R

R

R

R

với dịng kích từ cịn E’q tỉ lệ với từ thơng tổng móc vịng trong cuộn dây roto máy
R

R

phát. Q trình q độ điện từ trong cuộn dây roto xác định quan hệ giữa các sđđ này
(phương trình (1.9)).
Với các nút tải Pi và Q i xác định theo đặc tính động của thiết bị dùng điện. Một
R

R

R

R

mơ hình phụ tải động hỗn hợp là một vấn đề rất khó trong hệ thống điện lớn và khơng
có mơ hình nào phản ánh đúng đặc trưng tất cả các loại phụ tải, mức độ chi tiết cũng
như tính chính xác. Trong nghiên cứu ổn định nhất là ổn định đồng bộ, hay trong giai
đoạn đầu của QTQĐ (đến khoảng 10s sau sự cố), hay nếu điểm ngắn mạch không ở
gần nút tải người ta thường cho phép thay thế bằng đặc tính tĩnh. Theo báo cáo nghiên
cứu IEEE về nghiên cứu ổn định trong chu kỳ đầu của QTQĐ; mơ hình có thể chấp
nhận được nhất là mơ hình cơng suất tác dụng theo dịng điện khơng đổi (constant
current) và mơ hình cơng suất phản kháng theo điện dẫn không đổi (constant

admittance).
Các nút trung gian có giá trị số Pi = 0 và Q i = 0.
R

R

R

17

R


Như vậy quan hệ giữa chế độ máy phát và thông số trạng thái lưới điện được
thiết lập thông qua công suất và điện áp nút đầu cực máy phát.
1.3.

CÁC PHƯƠNG PHÁP PHÂN TÍCH TÍNH TỐN ỔN ĐỊNH ĐỘNG.
Để giải hệ phương trình vi phân của hệ thống điện, phương pháp tích phân số

thường được sử dụng phổ biến và hiệu quả [2], [4].
Hiện nay có nhiều phương pháp tích phân số tính tốn được hệ phương trình
QTQĐ hệ thống điện như phương pháp Euler, Euler - Cosi, Runge - Kutta, hay Adam Sterne. Trong đó phương pháp một bước Euler và Runge - Kutta được sử dụng phổ
biến, có thể cho ta xác định điểm tiếp theo của QTQĐ khi biết được thơng số của một
điểm trước đó. Ngồi ra cịn có những phương pháp xác định điểm tiếp theo bằng
nhiều điểm trước đó, hay cịn gọi là phương pháp nhiều bước.
Cho hệ phương trình vi phân có dạng: dx/dt=f(x, t) với x = x0 tại t = t
R

R

R

0

Trong đó: x là véc tơ trạng thái của n biến phụ thuộc, t là biến độc lập (thời gian).
1.3.1.Phương pháp Ơle
Từ

dx
= f (x, t) với x=x0 tại t=t 0
dt
R

R

R

Ta có:

∆x =

dx
dt

.∆t
x= x0

x 1 = x 0 + ∆x = x 0 +

dx
dt

.∆t
x= x0

Lấy gần đúng hai số hạng đầu của dãy Taylo tại x=x0
R

18




x1 = x 0 + ∆t(x 0 ) +

∆t 2  ∆t 3 
x0+
x 0 + ...
2!
3!

Sau khi xác định được x=x 1 tương ứng với t=t 1, ta có thể xác định tiếp theo trạng thái
R

R

R

R

x=x 2 tương ứng t=t 1 +∆t như sau:
R

R

R

R

x2 = x1 +

dx
dt

.∆t
x= x1

Vì phương pháp giải x ở đạo hàm cấp 1 nên ta gọi là phương pháp bậc 1.
Để tăng độ chính xác, ∆t phải đủ nhỏ. Áp dụng phương pháp này phải chú ý đến vệc
nâng cấp của sai số vì một sai số nhỏ bước trước có thể cho sai số lớn ở bước tiếp theo.
1.3.2.Phương pháp Ơle biến thể
Phương pháp Ơle trên thường khơng chính xác do dùng đạo hàm từ thời điểm
đầu của khoảng chia thời gian. Phuơng pháp Ơle biến thể khắc phục nhược điểm trên
bằng đạo hàm điểm trung bình khoảng chia:
-

Bước dự báo: dùng phép đạo hàm ở thời điểm đầu của khoảng chia đầu tiên, để
dự báo bước tiếp theo

x p1 = x 0 +

dx
∆t
dt x = x0

-Bước hiệu chỉnh dùng giá trị dự báo x1 p, lấy trung bình của đạo hàm điểm đầu
R

RP

P

và điểm dự báo

1  dx
x c1 = x 0 + 
2  dt

19

+
x =x 0

dx
dt


 ∆t
x =x p1 