intTypePromotion=1
zunia.vn Tuyển sinh 2024 dành cho Gen-Z zunia.vn zunia.vn
ADSENSE

Tạp chí khoa học và công nghệ: Chỉnh định mờ tham số PID cho bộ điều tốc tuabin nhà máy thủy điện

Chia sẻ: Trinh Nguyễn | Ngày: | Loại File: PDF | Số trang:8

159
lượt xem
20
download
 
  Download Vui lòng tải xuống để xem tài liệu đầy đủ

Bài báo này giới thiệu việc xây dựng mô hình một nhà máy thủy điện với sự có mặt các phần tử thủy lực, máy điện và xét đầy đủ các hiệu ứng động học của hệ thống thủy lực. Mời các bạn cùng tham khảo để nắm nội dung kiến thức cần thiết.

Chủ đề:
Lưu

Nội dung Text: Tạp chí khoa học và công nghệ: Chỉnh định mờ tham số PID cho bộ điều tốc tuabin nhà máy thủy điện

  1. TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ 4(39).2010 CHỈNH ĐỊNH MỜ THAM SỐ PID CHO BỘ ĐIỀU TỐC TUABIN NHÀ MÁY THỦY ĐIỆN FUZZY ADJUSTING PID PARAMETERS OF THE HYDROELECTRIC POWER TURBINE GOVERNOR Đoàn Quang Vinh Đặng Trung Thi Trường Đại học Bách Khoa, Đại học Đà Nẵng Công ty thủy điện Ialy TÓM TẮT Hiện nay, hầu hết các nhà máy thủy điện đều sử dụng bộ điều khiển PID để điều chỉnh tốc độ tuabin, các bộ điều chỉnh này ít có khả năng bền vững đối với sự thay đổi tham số trong quá trình vận hành. Bài báo này giới thiệu việc xây dựng mô hình một nhà máy thủy điện với sự có mặt các phần tử thủy lực, máy điện và xét đầy đủ các hiệu ứng động học của hệ thống thủy lực. Từ đó ứng dụng điều khiển mờ để điều chỉnh các tham số PID của bộ điều chỉnh tốc độ tuabin nhằm nâng cao chất lượng điều khiển và đáp ứng với sự thay đổi các tham số trong quá trình vận hành. ABSTRACT Nowadays, most of the hydroelectric power has used PID controllers to regulate turbine speed. However, these controllers are less robust with parameters changes in operation. This paper presents the modelling of a hydroelectric power with hydraulic components and electric machines in consideration of dynamic effects of the hydraulic system. Hence, the use of fuzzy control to adjust PID parameters of a turbine governor is aimed to improve the control quality and satisfy with parameters changes in operation. 1. Đặt vấn đề Trong những thập niên gần đây điều khiển mờ phát triển khá mạnh mẽ. Ưu điểm của điều khiển mờ so với các phương pháp điều khiển kinh điển là có thể tổng hợp được bộ điều khiển mà không cần biết trước cấu trúc và tham số của hệ thống một cách chính xác, giải quyết được nhiều bài toán điều khiển phức tạp mà trước đây chưa thể giải quyết trọn vẹn. Hệ điều khiển mờ sử dụng được các kinh nghiệm vận hành đối tượng và các xử lý điều khiển của chuyên gia trong thuật toán điều khiển, do vậy hệ điều khiển mờ là một bước tiến gần tới tư duy con người. Việc ứng dụng kỹ thuật mờ trong thiết kế xây dựng hệ điều khiển cho truyền động có cấu trúc và tham số biến đổi là hướng nghiên cứu mới mẻ, còn nhiều tiềm năng để khai thác, là Hình 1. Cấu trúc một nhà máy thủy điện 324
  2. TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ 4(39).2010 hướng nghiên cứu có khả năng đáp ứng được các yêu cầu chất lượng của hệ và khắc phục được các nhược điểm của các hướng nghiên cứu khác. Tuabin thủy lực cùng với các cơ cấu/bộ phận thủy lực đặc trưng động học thủy lực phi tuyến và có đặc tính pha không cực tiểu có các tham số thay đổi theo thời gian. Các bộ điều khiển PID được thiết kế cho một điểm làm việc đặc trưng, không thể tạo ra sự làm việc ổn định khi có sự thay đổi đột ngột trong điều kiện làm việc hoặc trong cấu hình của đối tượng. Sử dụng điều khiển mờ để điều chỉnh các tham số PID nhằm đạt được chất lượng điều khiển tốt trong điều kiện có sự thay đổi tham số hệ thống là mục đích của đề tài này. 2. Kết quả nghiên cứu và khảo sát 2.1. Mô hình hóa nhà máy thủy điện 2.1.1. Mô hình phi tuyến: H Q2 + H l 0 + H l 2 Hệ thống thủy lực của Cs f0 f p2 U s , Qs một NMTĐ điển hình, đầy đủ Tec U c , Qc Twc gồm có: Hồ chứa nước, đường H0 Twp f p1 Ht Hr ống dẫn nước, tháp điều áp, Tep U t , Qt đường ống áp lực, tuabin và H Q + H l1 đường ống xả nước xuống hạ lưu của nhà máy. Hình 2. Phân bố tham số và biến quá trình chính Đề tài này sẽ sử dụng các số liệu của NMTĐ Susqueda (Tây Ban Bảng 1. Trị số của các tham số NMTĐ Susqueda Nha), có công suất lắp đặt (2x37MW+1x12 Tham số Trị số Tham số Trị số MW), và xét riêng một tổ máy 37MW, sử Twc 9,15 [s] fp1 0,0475 [pu] dụng tuabin loại Francis. Twp 0,82 [s] fp2 0,089 [pu] Các tham số được tính từ số liệu của Tep 0,208 [s] At 1,67 [pu] NMTĐ theo các công thức trong tài liệu [7] Tg 0,5 [s] U nl 0,13 [pu] zp 3,95 T 225 [s] trình bày như ở bảng 1. Cs 140 [s] H 4 [s] Mô hình phi tuyến [8, tr. 77] như sau: f0 0,005 [pu] D 0,01 [pu] Trong đó η (G ) là hàm hiệu suất của tuabin xây dựng theo phương pháp bình phương nhỏ nhất với hàm polyfit trong Matlab, số liệu lấy từ [7] như bảng 2 và bậc đa thức là 4. η (G) |n=4 = 1,7783.G 4 - 4,4784.G 3 + 3,1593.G 2 - 0,7300.G + 0,9491 2.1.2. Mô hình trạng thái phi tuyến: Đặt các biến quá trình như sau: x1 = U t ; x2 = H r ; x3 = U c ; x4 = ω r ; x5 = G ; dU c ; dU t ; d 2U t . x6 = x3 = x7 = x1 = x8 = x7 = dt dt dt 2 325
  3. TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ 4(39).2010 Bảng 2. Trị số hàm hiệu suất Uc Abs Π f p2 H l2 H0 theo vị trí van hướng 1 G [ pu ] P [ MW ] η (G) [ pu ] Qc Twc .s Σ 0,130 13,10 0,8922 Abs Π f0 H l0 0,180 13,60 0,9012 Us Qs 1 Hs Hr U NL 0,250 14,20 0,9107 Σ Cs .s Σ Ut Σ 0,360 14,50 0,9111 Hr G Qt HQ Ht Pmec 1 ωr 0,411 15,30 0,9113 Π Z p . tanh(Tep .s ) Σ Π At Π Σ Σ 2 H .s Ut H l1 0,603 27,30 0,9180 Π f p1 D Pload 0,664 30,15 0,9048 Ht η (G ) 0,752 30,20 0,8410 G 1 uG 0,800 30,75 0,8174 Σ 1 s Tg 0,850 31,30 0,7874 0,896 31,80 0,7610 Hình 3. Mô hình phi tuyến NMTĐ với hàm hiệu suất η (G ) Khai triển từ các phương trình vi phân, ta được hệ 8 phương trình vi phân mô tả trạng thái phi tuyến của nhà máy điện: ⎛ x7 ⎞ ⎜ ⎟ 1 ⎜ . ( x3 − x1 ) (1 − 2 f 0 .C s . ( x6 − x7 ) ) ⎟ ⎜ Cs ⎟ ⎜ ⎟ 1 ⎜ ⎜ Twc ( H 0 − x 2 − f p 2 . x3 x 3 ) ⎟ ⎟ ⎜ ⎟ ⎜ D P η ( x ). A x 2 ⎟ − . x4 − load + 5 t .( x1 − U nl ). 12 ⎜ 2H 2H 2H x5 ⎟ ⎛ x1 ⎞ ⎜ 1 ⎟ ⎜ ⎟ ⎜ . ( u − x ) ⎟ ⎜ 2⎟ ⎜ x ⎟ G 5 Tg ⎜ x3 ⎟ ⎜ ⎟ ⎜ 1 ⎟ ⎜ ⎟ − . ⎡⎣ ( x3 − x1 ) . (1 − 2 f 0 .C s . ( x6 − x7 ) ) + 2 f p 2 .C s . x3 . x6 ⎤⎦ ⎜ x 4 ⎟ = ⎜ T .C ⎟ ⎜ x5 ⎟ ⎜ wc s ⎟ (1) ⎜ ⎟ ⎜ x8 ⎟ ⎜ x 6 ⎟ ⎜ ⎟ ⎜ x ⎟ ⎜ π 2 1 ⎛ 1 ⎞ ⎛ π . x12 2 ⎞ ⎟ . x − ⎜ 7 ⎟ ⎜ Twp .Tep2 2 Twp ⎜⎝ p1 x52 ⎟⎠ ⎜ Tep2 . f + . ⎜ + 8. ( 7 1 8 ) ⎟⎟ x 2 + x . x + ⎟ ⎜ x ⎟ ⎜ ⎝ ⎠ ⎟ ⎝ 8⎠ ⎜ ⎟ ⎜⎜ ⎛ 32 x π 2 ⎞ 8.( u − x ) x 2 ⎛ 3.( u − x ) ⎞ . 13 . ( uG − x5 ) − 2 ⎟ . x7 − G 5 . 13 . ⎜ 1 + G 5 ⎟ + ⎟ ⎜ ⎜ Twp .Tg x5 Tep ⎟⎠ Twp .Tg 2 x5 ⎝ x5 ⎠ ⎟ ⎜⎝ ⎟ ⎜ 4 ⎡(x − x ) 2 f 0 . ( x3 − x1 ) 2 ⎟ − 2 f 0 . ( x6 − x 7 ) + . (1 − 2 f 0 .C s . ( x6 − x7 ) ) ⎟ 2 ⎜ .⎢ 6 7 ⎜ Twp ⎢⎣ C s Twc .C s ⎟ ⎜ ⎟ ⎜ ⎛ 2 f p 2 . x3 . x 6 ⎞⎤ ⎟ ⎜⎜ + 2 f 0 . ( x3 − x1 ) . ⎜ + x8 ⎟ ⎥ ⎟⎟ ⎝ ⎝ T wc ⎠⎦ ⎠ 2.1.3. Mô hình trạng thái tuyến tính lân cận điểm làm việc: Hệ phi tuyến (1) được đưa về dạng phương trình trạng thái: ⎧ x = f ( x, u ) ⎨ (u = uG ) (2) ⎩ y = g ( x, u ) Trong đó: 326
  4. TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ 4(39).2010 x = f ( x, u ) là hệ các phương trình vi phân (1). y = g ( x, u ) phụ thuộc vào tín hiệu ra của hệ. Ở đây ta quan tâm đến tốc độ của tuabin, vì vậy: y = x4 (3) * Xung quanh điểm làm việc ( x , uGxl ) , với giả thiết các sai lệch so với trạng thái xác lập đủ nhỏ để có thể bỏ qua các thành phần bậc cao, khi đó có thể mô tả hệ bằng phương trình trạng thái tuyến tính [1]: ⎧ dx ⎪ = A x + B u G (4) ⎨ dt ⎪⎩ y = C x + D u G ⎧ x = x − x * ⎛ x1* ⎞ Trong đó: ⎪ ; x* = ⎜ # ⎟ và u G xl = ( u G xl ) (5) ⎨ u G = u G − u G xl ⎜ ⎟ ⎪ ⎜ x* ⎟ ⎩ y = y − g ( x , u G xl ) * ⎝ n⎠ Với các ma trận A,B,C,D là các ma trận Jacobi của các vectơ hàm f ( x, uG ) , g ( x, uG ) tại điểm làm việc đang xét (ma trận hằng): [1] B11 là ma trận nhiễu tải: ⎛ ∂f ⎞ (6) B1 1 = ⎜ ⎟ ⎝ ∂ u ∆ P L ⎠ x * ,u G xl ⎧ dx Khi đó (4) có thể mô tả lại như sau: ⎪ = Ax + BuG + B11u∆PL (7) ⎨ dt ⎪⎩ y = Cx + DuG Mô hình tuyến tính hóa lân cận điểm làm việc được giới thiệu ở hình 4 có các ma trận A, B, C, D với kích thước như sau A ∈ \ 8 x 8 , B ∈ \ 8 x1 , B11 ∈ \ 8 x1 , C ∈ \1 x 8 , D ∈ Θ . Hình 4. Cấu trúc mô hình trạng thái tuyến tính lân cận điểm làm việc 2.2. Chỉnh định mờ tham số PID bộ điều tốc: 2.2.1. Phương pháp chỉnh định mờ tham số bộ điều khiển PID: Cơ sở của phương pháp này là dựa vào việc phân tích sai lệch e(t) và đạo hàm của sai Hình 5. Phương pháp chỉnh định mờ lệch de(t)/dt, các tham số KP, TI, TD của bộ tham số bộ điều khiển PID 327
  5. TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ 4(39).2010 điều khiển PID sẽ được tự động chỉnh định theo phương pháp chỉnh định mờ của Zhao, Tomizuka và Isaka [3]. Một bộ điều khiển PID với đầu vào e(t), đầu ra u(t) có mô hình toán học như sau: ⎡ 1 t ⎤ u (t ) = K R ⎢ e (t ) + ∫ e (τ )d τ + T e (t )⎥ D (8) ⎣ TI 0 ⎦ hoặc G PID (s ) = K R + K I + K D s (9) s KR KD T1 = TD = trong đó: KI và KR 2.2.2. Thiết kế bộ điều khiển mờ chỉnh định thông số PID a) Xác định biến ngôn ngữ Bộ chỉnh định mờ có 2 đầu vào là e(t) và de(t)/dt và ba đầu ra là KR, KD, α có các giới hạn xác định theo [3]. Trong đó: α TI K R2 nên xem như gồm 3 bộ chỉnh định mờ nhỏ. = ⇒ KI = TD α .K D µ NB NM NS ZE PS PM PB Định nghĩa các tập mờ vào ra: e = {NB NM NS ZE PS PM PB} de(t ) = { NB NM NS ZE PS PM PB} e, de/dt dt -1 0 1 K R = {S B} µ µ K D = {S B} S B S MS M P α = {S MS M B} K R, KD α 0 1 1 2 3 4 5 6 Hình 6. Định nghĩa các tập mờ vào ra b) Xây dựng luật hợp thành Bảng 3. Luật chỉnh định Luật chỉnh định KR: Luật chỉnh định KD: Luật chỉnh định α: de(t)/dt de(t)/dt de(t)/dt NB NM NS ZE PS PM PB NBNMNSZEPSPMPB NBNM NS ZE PS PM PB NB B B B B B B B NB S S S S S S S NB S S S S S S S NM S B B B B B S NM B B S S S B B NM MS MS S S S MSMS e(t) NS S S B B B S S e(t) NS B B B S B B B e(t) NS M MS MS S MS MS M ZE S S S B S S S ZE B B B B B B B ZE B M MS MS MS M B PS S S B B B S S PS B B B S B B B PS M MS MS S MS MS M PM S B B B B B S PM B B S S S B B PM MS MS S S S MSMS PB B B B B B B B PB S S S S S S S PB S S S S S S S 328
  6. TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ 4(39).2010 c) Chọn luật và giải mờ: Chọn luật hợp thành max-min, phương pháp giải mờ trọng tâm. d) Mô hình và kết quả mô phỏng: Hình 7. Mô hình xây dựng trong Matlab Hình 8. Cấu trúc khối PID-FC Hình 9. Cấu trúc khối FC Hình 10. Đáp ứng tốc độ theo tải với thông Hình 11. Đáp ứng tốc độ theo tải khi số tại bảng 1 thay đổi tham số Twp Hình 12. Đáp ứng tốc độ theo tải (mô hình Hình 13. Đáp ứng tốc độ theo tải (mô hình PID-FC) khi thay đổi thông số Tep PID) khi thay đổi thông số Tep 329
  7. TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ 4(39).2010 Ghi chú: PID-FC: Mô hình điều khiển chỉnh định mờ PID PID: Mô hình điều khiển PID kinh điển Hình 14. Đáp ứng tốc độ theo tải và nhiễu 3. Bàn luận Khi thay đổi các tham số của hệ thống, cụ thể trong trường hợp này là tham số Twp, Tep; đáp ứng tốc độ của mô hình PID-FC gần như không đổi, trong khi đó đáp ứng tốc độ của mô hình PID bị thay đổi rất nhiều. Ứng với sự thay đổi này, đáp ứng tốc độ đối với mô hình PID là dao động nhiều, thời gian quá độ rất lớn (hình 11) và thậm chí mất ổn định (hình 13). Lý do là khi có sự thay đổi thông số khi mô hình hóa hoặc trong khi vận hành, làm dịch chuyển điểm làm việc xác lập mà các thông số PID không thay đổi trong trường hợp này dẫn đến dao động hoặc mất ổn định. Còn đối với mô hình PID-FC nhờ sự chỉnh định mờ thông số PID mà thích nghi được với sự thay đổi này, giúp cho bộ điều khiển bền vững trong các điều kiện trên. 4. Kết luận Bài báo đã giới thiệu mô hình hóa một nhà máy thủy điện với sự tham gia đầy đủ các phần tử thủy lực, máy điện. Trên cơ sở mô mình đó, áp dụng kỹ thuật điều khiển mờ để chỉnh định các tham số PID cho bộ điều tốc nhà máy thủy điện. Với những ưu điểm của điều khiển mờ, phương pháp điều khiển trên đã khắc phục được điểm yếu của điều khiển PID kinh điển-đó là khả năng bền vững trước sự thay đổi của các tham số của hệ thống trong quá trình mô hình hóa cũng như trong quá trình vận hành. Kết quả mô phỏng cũng cho thấy chất lượng điều khiển khá tốt trước nhiễu tải, nhiễu đo lường và nhiễu hệ thống. Đây là một hướng nghiên cứu mới triển vọng và có khả năng ứng dụng vào thực tiễn. TÀI LIỆU THAM KHẢO Tiếng Việt [1] Nguyễn Doãn Phước , Lý thuyết điều khiển phi tuyến, NXB khoa học và kỹ thuật, Hà Nội, 2007. [2] Nguyễn Doãn Phước, Lý thuyết điều khiển nâng cao, NXB khoa học và kỹ thuật, Hà Nội, 2007. [3] Phan Xuân Minh, Nguyễn Doãn Phước, Lý thuyết điều khiển mờ, NXB khoa học và kỹ thuật, Hà Nội, 2002. 330
  8. TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ 4(39).2010 Tiếng Anh [4] Glenn O.Brown, “The History of the Darcy-Weisbach Equation for Pipe Flow Resistance”, Environmental and Water Resources History, American Society of Civil Engineers, Reston, VA, pp. 34-43, 2003. [5] H. Doudna, F.P. de Mello, R.J. Koessle, J. Agee, P.M. Anderson, J.H. Fish III, P.A.L. Hamm, P. Kundur, D.C. Lee, G.J. Rogers and C. Taylor, “Hydraulic turbine and turbine control models for system dynamic studies”, IEEE Transactions on Power Systems, 1992, Vol.7 (No.1), pp.167-179. [6] Kundur, Power System Stability and Control, McGraw-Hill, New York, 1994. [7] O. Quiroga, J. Riera and C. Batlle, “Identification of partially known models of the Susqueda hidroelectric power plant”, Latin American Applied Research, 2003, Vol.33 (No. 4), pp. 387-392. [8] Quiroga, O. D, Modelling and Nonlinear Control of Voltage Frequency of Hydroelectric Power Plants, Thesis Dissertation, Universitat Politècnica de Catalunya, Barcelona, Spain, 2000. 331
ADSENSE

CÓ THỂ BẠN MUỐN DOWNLOAD

 

Đồng bộ tài khoản
2=>2