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

Báo cáo nghiên cứu khoa học: " ĐIỀU KHIỂN LQ CHO TỐC ĐỘ TUABIN NHÀ MÁY THỦY ĐIỆN"

Chia sẻ: Nguyễn Phương Hà Linh Linh | Ngày: | Loại File: PDF | Số trang:8

Thêm vào BST
Báo xấu
134
lượt xem 19
download
 
  Download Vui lòng tải xuống để xem tài liệu đầy đủ

Mô hình nhà máy th điện (NMTĐ) là một mô hình phi tuyến, trong đó hiệu ứng đàn ủy hồi của cột nước trong ống áp lực được biểu diễn bằng hàm toán học vô tỷ. Bài viết này sẽ giới thiệu mô hình phi tuyến của một NMTĐ, xấp xỉ hàm vô tỷ trong hệ bằng phương pháp tham số tập trung, và xây dựng mô hình không gian trạng thái phi tuyến bậc 8 cho hệ. Bộ điều chỉnh tốc độ tuabin dựa trên lý thuyết điều khiển tuyến tính LQ cũng được xây dựng cho mô hình...

Chủ đề:

Bình luận(0) Đăng nhập để gửi bình luận!

Lưu

Nội dung Text: Báo cáo nghiên cứu khoa học: " ĐIỀU KHIỂN LQ CHO 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Ố 6(35).2009 ĐIỀU KHIỂN LQ CHO TỐC ĐỘ TUABIN NHÀ MÁY THỦY ĐIỆN THE CONTROL OF SPEED FOR HYDRAULIC TURBINES BASED ON THE LINEAR QUADRATIC CONTROL METHOD Đoàn Quang Vinh Lê Đức Dũng Đại học Đà Nẵng Điện lực Thừa Thiên Huế TÓM T ẮT Mô hình nhà máy th điện (NMTĐ) là một mô hình phi tuyến, trong đó hiệu ứng đàn ủy hồi của cột nước trong ống áp lực được biểu diễn bằng hàm toán học vô tỷ. Bài viết này sẽ giới thiệu mô hình phi tuyến của một NMTĐ, xấp xỉ hàm vô tỷ trong hệ bằng phương pháp tham số tập trung, và xây dựng mô hình không gian trạng thái phi tuyến bậc 8 cho hệ. Bộ điều chỉnh tốc độ tuabin dựa trên lý thuyết điều khiển tuyến tính LQ cũng được xây dựng cho mô hình hệ ở lân cận điểm làm việc xác lập. Kết quả mô phỏng so sánh đáp ứng điều chỉnh giữa các bộ điều chỉnh LQR, LQG và LQG/LTR. ABSTRACT The model of the hydro-power plant is nonlinear, and the effect of the elasticity of the water in the penstock is represented by an irrational mathematical function. This paper introduces the nonlinear model of a hydro power plant, which approximates the irrational function to reduce to a lower order using the lumped parameter approximation method, and the making of an elastic eighth order space state model. The turbine speed regulator for this model in the neighbourhood of the operating point based on Linear Control Theory is made and the comparative simulation results that meet the needs of adjustment for LQR, LQG and LQG/LTR are developed. 1. Đặt vấn đề Tính năng làm việc tuabin thủy lực chịu tác động mạnh bởi những đặc tính của hệ thống dẫn cấp nước cho tuabin: hiệu ứng quán tính, áp lực của nước và độ đàn hồi của vách đường ống áp lực. Ngoài ra nó mang đặc điểm của hệ pha không cực tiểu. Hiện nay các bộ điều chỉnh PID số đang được sử dụng rộng rãi để điều tốc cho tuabin NMTĐ. Một số nghiên cứu khác đã công b áp dụng các phương pháp đi ố ều khiển nâng cao đối với hệ điều tốc tuabin trong NMTĐ: Tối ưu trượt điểm cực, Bền vững, Tối ưu vô hướng[4], Lyapunov [3],… dựa trên nền tảng lý thuyết điều khiển tuyến tính. Hầu hết các nghiên cứu xét mô hình NMTĐ trong phạm vi đường ống áp lực - tuabin và máy phát. Ít có nghiên c u xét hệ đầy đủ, bao gồm cả đường ống dẫn, tháp ứ điều áp, và các đặc tính động học khác của hệ (tổn thất và hiệu ứng đàn hồi cột nước). 2. Mô hình động học của hệ thống thủy lực NMTĐ Hệ thống thủy lực của một NMTĐ điển hình gồm có: Hồ chứa nước, đường ống dẫn nước, tháp điều áp, đường ống áp lực, tuabin và đường ống xả nước xuống hạ lưu. 22
  2. TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ 6(35).2009 2.1. Mô hình phi tuyến 2.1.1. Các phương trình phi tuyến Mô hình phi tuyến của NMTĐ đề cập trong bài viết này được xét đầy đủ các đặc tính động học chất lỏng, chỉ bỏ qua hiệu ứng đàn hồi của cột nước trong đường ống dẫn. Các phương trình động học được tổng hợp từ các mô hình khác nhau (xét ở hệ đơn vị tương đối): [1],[2]. * Dòng chảy: (= AU / AU cb= U [pu]) Ut Uc − U s = ⇒Q (1) Q/Q . . cb * Động học của ống dẫn: dU c H Q 2 H 0 − H r − H l 2 == H l 2 = f p 2 .U c . U c (2) (3) dt Twc Twc * Động học của tháp điều áp: 1 Cs ∫ = U s dt − f 0 .U s . U s H l 0 = f 0 .U s . U s (4) (5) Hr * Động học của ống áp lực: - Trường hợp xét hiệu ứng đàn hồi cột nước: Ut = G Ht H l1 = f p1.U t 2 (6) (7) ( z p = Twp / Tep ) H Q ( s ) = z p .tanh(Tep .s ).U t ( s ) (8) H t ( s ) = H r ( s ) − H l1 ( s ) − H Q ( s ) (9) - Trường hợp không xét hiệu ứng đàn hồi cột nước, tanh(Tep .s ) ≅ Tep .s : (9) ⇔ H t ( s ) = H r ( s ) − H l1 ( s ) − z p .Tep s.U t ( s ) (10) * Công suất cơ của tuabin: Pmec = At H t (Qt − Qnl ) − Ddam .G.∆ωr = At H t (U t − U nl ) − Ddam .G.∆ωr (11) Trường hợp xét đến hiệu suất, với hàm hiệu suất là η = η (G ) : = η At H t (U t − U nl ) (12) Pmec * Chuyển động của tuabin: dω Pmec − Pload 2.H . r + D.ωr = (13) dt * Động học của động cơ servo: dG +G= Tg . uG (14) dt 23
  3. TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ 6(35).2009 Hình 1. Các tham số, biến quá trình của hệ thống thủy lực NMTĐ Trong bài viết này, tham số tính toán của hệ sẽ được tham khảo từ số liệu một tổ máy 37MW, sử dụng loại tuabin Francis của NMTĐ Susqueda (Tây Ban Nha)[2]. Bảng 1. Các biến quá trình và tham số của hệ Biến quá trình Cột nước hồ chứa; tuabin; tháp điều áp [pu]. H0 , Ht , H r Tổn thất cột nước tại tháp điều áp; ống áp lực; ống dẫn H l 0 , H l1 , H l 2 [pu]. Vận tốc, hoặc lưu lượng của dòng chảy trong ống dẫn; tháp Uc ,U s ,U p ,Ut , điều áp; ống áp lực; tuabin; không tải qua tuabin [pu]. U NL Vị trí mở van hướng [pu]. G Tín hiệu điều khiển động cơ servo [pu]. uG Công suất cơ của tuabin [pu]. Pmec ωr Tốc độ quay của tuabin (roto) [pu]. Tham số Hằng số thời gian khởi động của nước ở ống dẫn; ống áp Twc , Twp lực [s]. Twc = 9,15 ; Twp = 0,82 . Hằng số thời gian đàn hồi của ống áp lực [s]. Tep = 0, 208 . Tep Hằng số tích lũy của tháp điều áp [s]. Cs = 140 . Cs Hệ số tổn thất cột nước của tháp điều áp [pu]. f 0 = 0, 005* . f0 24
  4. TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ 6(35).2009 Hệ số tổn thất cột nước của ống áp lực [pu]. f p1 = 0, 0475 . f p1 Hệ số tổn thất cột nước của ống dẫn [pu]. f p 2 = 0, 089 . f p2 Hệ số khuyếch đại tuabin [pu]. At = 1, 67 . At Lưu lượng không tải [pu]. U NL = 0,13 . U NL Hằng số quán tính của hệ tuabin-máy phát [s]. H = 4** . H Hệ số cản của tuabin [pu]. Ddam Hệ số cản của tải (phụ thuộc tần số) [pu]. D = 0, 01*** . D Ghi chú: *,**,***: Tài liệu [2] không đề cập hoặc bỏ qua. 2.1.2. Mô hình không gian trạng thái Từ số liệu tham khảo ở [2], sử dụng phương pháp bình phương nhỏ nhất với hàm polyfit của Matlab, hàm η = η (G ) với đa thức bậc 3 được xây dựng: η (G ) = 1,7783.G 4 - 4,4784.G 3 + 3,1593.G 2 - 0,7300.G + 0,9491 Do hàm vô tỷ tanh(Tep .s ) rất khó sử dụng để nghiên cứu ổn định của hệ thống, ở đây sẽ xấp xỉ tanh(Tep .s ) theo phương pháp tham số tập trung (lumped parameter): n= n=  ∞ ∞  Tep s   2  2Tep s  −2.T s 2 1 − e ep Tep s∏   ∏ 1 +   tanh(Tep s ) =ep s =1 +    nπ   n   ( 2n − 1) π   −2.T 1+ e = 1= 1     n Mô hình tổ hợp từ các phương trình ở mục 2.1.1 có đáp ứng Pmec ở hình 3(a), cho thấy hàm tanh(Tep .s ) xấp xỉ với bậc n=1 khá chính xác để nghiên cứu động học của hệ. Mô hình không gian ạng thái phi tuyến (SS_NL) với b ậc xấp xỉ tanh(Tep .s ) | n=1 có tr dạng: Đặt biến: x1 = U t ; x2 = H r ; x3 = U c ; x4 = ωr ; x5 = G ;  dx  = f ( x, u ) (15)  dt x= x= dU c / dt ; x= x= dU t / dt ; x8 x= d 2U t / dt 2 = 7 3 1 6 7  y = g ( x, u )  (15) sẽ là hệ 8 phương trình vi phân. Đáp ứng của SS_NL| n=1 ở hình 3 (b). 2.2. Mô hình tuyến tính hóa lân cận điểm làm việc Để áp dụng lý thuyết điều khiển tuyến tính cho đối tượng phi tuyến, mô hình SS_NL| n=1 được thay thế bằng mô hình tuyến tính (SS_LN) tại điểm làm việc ( x* , uGxl ) : 25
  5. TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ 6(35).2009  dx  = Ax + Bu (16)  dt  y Cx + Du =  Các ma tr hằng A,B,C,D là các maận Jacobi của các vectơ hàm ận tr f ( x, u ) , g ( x, u ) tại điểm làm việc ( x* , uGxl ) : A ∈  8 x8 , B ∈  8 x1 , C ∈  1x8 , D ∈ Θ . 0.555 1 0.9 0.550 0.8 0.7 0.545 SS-LNn=0 Pmec/G [pu] Pmec [pu] 0.6 SS-LNn=1 0.540 0.5 goc van SS-NLn=1 0.4 Dac tinh tu nhien Xap xi tanh(Tep.s)≈Tep.s 0.535 0.3 Xap xi tanh(Tep.s)n=1 0.2 0.530 Xap xi tanh(Tep.s)n=2 0.1 Xap xi tanh(T .s) 0 ep n=3 0.525 995 1000 1005 1010 1015 1020 1025 1030 1035 0 100 200 300 400 500 600 700 800 900 1000 Thoi gian [s] Thoi gian [s] a) b) Hình 3. Đáp ứng Pmec :(a) Theo tanh(Tep .s ) và tanh(Tep .s ) |n = 0,1,2,3 ; (b) SS_NL và SS_LN Hình 3(b) cho th sai lệch đ áp ứng với G tăng 0,01pu từ trạng thái xác lập ấy G = 0,5 pu giữa mô hình SS_NL| n=1 và SS_LN| n=1 rất nhỏ. 3. Điều chỉnh LQ (Linear Quadratic) 3.1. Bộ điều chỉnh LQR (Linear Quadratic Regulator) phản hồi âm Mô hình SS_LN| n=1 đảm bảo điều kiện điều chỉnh được và quan sát được tại mọi điểm làm việc, tức là:  sI − A  Rank ( sI − A, B )= 8 ∀s ⊂  và Rank  = 8 ∀s ⊂  C Điều khiển LQR được đánh giá bởi phiếm hàm mục tiêu: [5] t =∞ 1 ∫ x (t ) E x(t ) + u (t ) Fu (t )  d t→ min = T T (17) Q( x(t ), u (t ))   2 t =0 Hệ số khuyếch đại của bộ điều chỉnh: RLQR = F −1.BT .K (18) w+ x u dx Σ = Ax + Bu Với K là ngh ệm xác đ n h b á i ị n dt - dương c a phương trình đại số Riccati ủ RLQR = F −1.BT .K (CARE). K .B.F −1.BT .K − AT .K − K . A = (19) E Phương pháp điều khiển LQR đòi hỏi tất cả các biến trạng thái của quá trình đều phải đo được, điều này vừa không thực tế, vừa không kinh tế trong mọi trường hợp. 26
  6. TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ 6(35).2009 3.2. Bộ điều chỉnh LQG (Linear Quadratic Gaussian): Mục tiêu của bài toán là ước lượng được các trạng thái của mô hình NMTĐ ở (16) xét với D ∈ Θ , loại trừ tác động của nhiễu vào hệ, đồng thời phiếm hàm mục tiêu (17) đạt giá trị nhỏ nhất.[5] nx ny +  dx  dx +  = Ax + Bu + nx (20)  = Ax + Bu + nx y w+ u + Σ Σ  dt  dt  y = Cx + ny  -  y Cx + n y =  ˆ u Trong đó nx , n y là các u  ˆ x ˆ dx RLQR = F −1.BT .K =− L.C ) x + ( B L).   ˆ (A nhiễu ồn trắng tác đ ng vào ộ  y dt RLQG quá trình và ở đầu ra, các ma Bộ quan sát Kalman RLQR trận tương quan là N x xác định bán dương, N y xác định dương. 3.2.1. Bộ quan sát trạng thái (lọc) Kalman Bộ quan sát Kalman, có nhiệm vụ ước lượng trạng thái của đối tượng tuyến tính (20), được xác định như (21)[5], với x là trạng thái ước lượng của x . ˆ u ˆ dx = Ax + Bu + L( y − Cx) = ( A − LC ) x + ( B L)   (21) ˆ ˆ ˆ  y dt − Trong đó L được tính: L = P.C T .N y 1 (22) Với P là nghiệm xác định bán dương của CARE. − P.C T .N y 1.C.P − P. AT − A.P = Nx (23) 3.2.2. Bộ điều chỉnh LQG Bộ điều chỉnh LQG là sự kết hợp giữa bộ điều chỉnh LQR và bộ quan sát trạng thái Kalman đ điều chỉnh và kháng nhiễu cho hệ tuyến tính (20), với đầu ra của bộ ể quan sát là đầu vào của bộ điều chỉnh LQR. 3.2.3. Loop Transfer Recovery (LTR) Ta có hàm truyền mạch hở của LQR: H LQR ( s ) RLQR .( sI − A) −1.B = (24) Và LQG: H LQG ( s ) RLQR .( sI − A + B.RLQR + L.C ) −1.LC.( sI − A) −1.B = (25) Với việc chọn N x = σ .B ' B hiệu chỉnh hệ số ρ để N y = ρ .N 0 ( N 0 là ma trận phổ nhiễu đo) đạt giá trị hợp lý, sao cho: lim H LQG ( s ) → H LQR ( s ) được gọi là LTR, khi đó ρ →0 đặc tính điều chỉnh của LQG sẽ tiệm cận đến đặc tính điều chỉnh của LQR, dù có sự tồn tại của nhiễu. Tuy nhiên với hệ pha không cực tiểu, thì kỹ thuật LTR ít có hiệu quả. 3.2.4. Kết quả mô phỏng 27
  7. TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ 6(35).2009 Các ma trận E, F, N x , N y có ý nghĩa rất quan trọng đến chất lượng của bộ điều chỉnh. Với hệ SISO đang xét, có biến điều chỉnh uG và biến ra ωr , đã chọn: -LQR: (i, j ) = 0 ∀i ≠ j , E 8 x8 (i, j ) = 1 (i = j = 4), E 8 x 8 (i, j ) = 10−3 (i = j ≠ 4), F = 0, 05. 8 x8 E - LQG(1): N x = I 8 x8 , N y = 10−1 ; LQG(2): N x = I 8 x8 .10−1 , N y = 10−3 . - LQG/LTR: N x = BB ' , N y = 10−5 . Hình 4 trình bàyđ áp ứng điều chỉnh tốc độ tuabin tương ứng với các bộ điều chỉnh LQR, LQG, LQG/LTR, khi có sự tăng/giảm mức tải 0,1pu ở điểm làm việc 0,5pu (10% tải định mức). Nhiễu quá trình và nhiễu đo có biên độ dao động khoảng ±0,005pu (±0,5%). Hình 5 thể hiện chi tiết đáp ứng ở hình 4 lân cận giá trị đặt ωr _ ref (1,0 pu). 1.025 1.0025 LQR 1.020 1.0020 LQG(1) LQG(2) 1.015 1.0015 LQG/LTR LQR LQG(2) 1.010 1.0010 LQR LQG(2) Toc do ωr [ pu] Toc do ωr [ pu] 1.005 1.0005 1.000 1.0000 0.995 0.9995 LQG/LTR 0.990 0.9990 LQG(1) LQG(1) LQG/LTR 0.985 0.9985 0.980 0.9980 0.975 0.9975 0 10 20 30 40 50 60 70 80 90 100 0 10 20 30 40 50 60 70 80 90 100 Thoi gian [s] Thoi gian [s] Hình 4. Đáp ứng tốc độ ωr của các bộ điều chỉnh Hình 5. Chi tiết hình 4 4. Kết luận Kết quả nghiên cứu, phân tích và đánh giá kết quả mô phỏng cho thấy: - Đáp ứng điều chỉnh của các bộ điều khiển LQ đã xét rất tốt và tốt hơn nhiều về quá điều chỉnh và thời gian xác lập so với một số phương pháp khác (PID, Lyapunov - so sánh với kế t quả nghiên cứu tại tài liệu [3]). Kết quả mô phỏng cũng cho thấy dù các bộ điều chỉnh LQG (có xét đến nhiễu Gaussian) không đạt được chất lượng đáp ứng như bộ điều chỉnh LQR (không xét nhiễu), nhưng đáp ứng ra có dao động rất nhỏ, khoảng ±0,0005pu (±0,05%), hay nói cách khác, các ộb điều chỉnh LQG đã loại đến 90% tác động của nhiễu quá trình và nhiễu đo vào hệ khi hệ xác lập. - Sử dụng kỹ thuật LTR cho bộ điều chỉnh LQG mang lại dạng đặc tính điều chỉnh gần giống đặc tính điều chỉnh LQR. Tuy nhiên, do đặc điểm của hệ pha không cực tiểu, nên khó tìm được bộ điều chỉnh LQG/LTR có thể khử hoàn toàn tác động của nhiễu để đạt được đáp ứng ”bằng phẳng” như LQR. TÀI LIỆU THAM KHẢO 28
  8. TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ 6(35).2009 [1] Hydraulic turbine and turbine control models for system dynamic studies, IEEE Trans. on Power Systems, Vol.7, No.1 pp. 167÷179, February 1992. [2] Identification of partially known models of the Susqueda hidroelectric power plant, O. Quiroga, J. Riera and C. Batlle-Latin American Applied Research v.33 n.4 Bahía Blanca oct./dic. 2003. [3] Modelling and Nonlinear Control of Voltage Frequency of Hydroelectric Power Plants, Quiroga, O. D., Thesis Dissertation, Universitat Politècnica de Catalunya, Barcelona, Spain (2000). [4] A review on hydropower plant models and control, Nand Kishor-R.P.Saini-S.P. Singh; Science Direct - Renewable and Sustainable Energy Reviews 11 (2007) [5] Lý thuyết điều khiển nâng cao - Nguyễn Do ãn Phước - NXB Khoa học và Kỹ thuật - 2007.(trang 97÷101, 172÷182). 29
ADSENSE

CÓ THỂ BẠN MUỐN DOWNLOAD

LV.11: Bộ Luận Văn Tốt Nghiệp Chuyên Ngành Tài Chính Ngân Hàng
172 tài liệu
843 lượt tải
THÔNG TIN
TRỢ GIÚP
HỖ TRỢ KHÁCH HÀNG
Theo dõi chúng tôi

Chịu trách nhiệm nội dung:

Nguyễn Công Hà - Giám đốc Công ty TNHH TÀI LIỆU TRỰC TUYẾN VI NA

LIÊN HỆ

Địa chỉ: P402, 54A Nơ Trang Long, Phường 14, Q.Bình Thạnh, TP.HCM

Hotline: 093 303 0098

Email: support@tailieu.vn

Giấy phép Mạng Xã Hội số: 670/GP-BTTTT cấp ngày 30/11/2015 Copyright © 2022-2032 TaiLieu.VN. All rights reserved.

 

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