طراحی کنترلکننده میراگر بر پایه مدل شناساییشده با استفاده از دادههای حاصل از اندازهگیرهای فازوری حوزه وسیع
محورهای موضوعی : مهندسی برق و کامپیوترآذین عطاردی 1 , هیمن گلپیرا 2 * , حسن بیورانی 3
1 - دانشگاه کردستان،دانشکده مهندسی
2 - دانشگاه کردستان،دانشکده مهندسی
3 - دانشگاه کردستان،دانشکده مهندسی
کلید واژه: سیستمهای پایش و کنترل حوزه وسیع, کنترلکننده میراگر, مد بین ناحیهای, مدل مرتبه پایین,
چکیده مقاله :
تغییرات مداوم و پیچیدگی سیستمهای قدرت مدرن، باعث بروز چالشهایی در زمینه مدلسازی آنها شده است. امروزه با پیشرفت سیستمهای پایش حوزه وسیع، دادههای حاصل از این سیستمها میتوانند در شناسایی و تخمین مدل سیستمهای قدرت به کار روند. این مقاله بر تنظیم پایدارساز سیستم قدرت با استفاده از مدل شناساییشده به کمک اندازهگیریهای فازوری حوزه وسیع تمرکز دارد. مدل مرتبه پایین شناساییشده بر اساس اندازهگیریها دربرگیرنده خصوصیات دینامیکی مدهای غالب بین ناحیهای بوده و میتواند برای طراحی کنترلکننده میراگر و ارزیابی اثربخشی آن در سیستم قدرت واقعی استفاده شود. کنترلکننده از نوع پایدارساز سیستم قدرت انتخاب شده و به دو روش مقاوم و زیگلر- نیکولز طراحی شده است. نتایج عددی، اثربخشی این رویکرد را در بهبود میرایی مد بین ناحیهای سیستم دوناحیهای 4ماشینه با استفاده از دادههای فازوری و دینامیکی به دست آمده از شبیهسازی در نرمافزار MATLAB نشان میدهد.
Continuous changes besides increasing complexities of modern power systems cause emergence of new challenges in modeling of power systems. Nowadays, with development of wide-area monitoring systems, data from the overall system can be used to identify and estimate model for power systems. This paper focuses on power system stabilizer tuning using the derived measurements-based model. The derived low-order model includes dynamic characteristics of inter-area dominant modes and can be used to design the damping controller and evaluate its effectiveness in power system studies. The controller can be reinterpreted as power system stabilizer and may be designed in two different methods of i) robust and ii) Ziegler-Nichols. The numerical results show the effectiveness of this approach in improving the small signal stability behavior of two-area 4-machine system using measured data, obtained from time domain simulation in MATLAB software.
[1] م. زمانی و غ. شاهقلیان، "طراحی هماهنگ پایدارساز سیستم قدرت و ادوات امپدانس متغیر برای افزایش میرایی مدهای بین ناحیهای با استفاده از الگوریتم ژنتیک،" مهندسی برق و مهندسی کامپیوتر ایران- الف مهندسی برق، سال 17، شماره 4، صص. 278-271، پاییز 1398.
[2] H. Golpira, A. R. Messina, and H. Bevrani, Renewable Integrated Power System Stability and Control, Wiley-IEEE, USA, 2021.
[3] T. Hashiguchi, et al., "Identification of characterization factor for power system oscillation based on multiple synchronized phasor measurements," Electrical Engineering in Japan, vol. 163, no. 3, pp. 10-18, May 2008.
[4] F. Al Hasnain, A. Sahami, and S. Kamalasadan, "An online wide-area direct coordinated control architecture for power grid transient stability enhancement based on subspace identification," IEEE Trans. on Industry Applications, vol. 57, no. 3, pp. 2896-2907, May/Jun. 2021.
[5] S. N. Sarmadi and V. Venkatasubramanian, "Electromechanical mode estimation using recursive adaptive stochastic subspace identification," IEEE Trans. on Power Systems, vol. 29, no. 1, pp. 349-358, Jan. 2013.
[6] H. Liu, et al., "ARMAX-based transfer function model identification using wide-area measurement for adaptive and coordinated damping control," IEEE Trans. on Smart Grid, vol. 8, no. 3, pp. 1105-1115, May 2015.
[7] R. L. Ribler, J. S. Vetter, H. Simitci, and D. A. Reed, "Autopilot: adaptive control of distributed applications," in Proc. The Seventh Int. Symp. on High Performance Distributed Computing, pp. 172-179, Chicago, IL, USA, 31-31 Jul.1998.
[8] F. Bai, et al., "Design and implementation of a measurement-based adaptive wide-area damping controller considering time delays," Electric Power Systems Research, vol. 130, no. 2, pp. 1-9, Jan. 2016.
[9] A. Hasanovic, A. Feliachi, N. Bhatt, and A. DeGroff, "Practical robust PSS design through identification of low-order transfer functions," IEEE Trans. on Power Systems, vol. 19, no. 3, pp. 1492-1500, Aug. 2004.
[10] N. Kishor, L. Haarla, J. Turunen, M. Larsson, and S. R. Mohanty, "Controller design with model identification approach in wide area power system," IET Generation, Transmission & Distribution, vol. 8, no. 8, pp. 1430-1443, Aug. 2014.
[11] T. Prakash, V. P. Singh, and S. R. Mohanty, "A synchrophasor measurement based wide-area power system stabilizer design forinter-area oscillation damping considering variable time-delays," International J. of Electrical Power and Energy Systems, vol. 105, pp. 131-141, Feb. 2019.
[12] س. اباذری، ع. عرب دردری، م. برخورداری یزدی و م. ص. پیام، "طراحی کنترلکننده مقاوم محدوده وسیع SVC جهت میرایی نوسانات بین ناحیهای در سیستم قدرت،" مهندسی برق و مهندسی کامپیوتر ایران- الف مهندسی برق، سال 13، شماره 1، صص. 46-36، بهار 1394.
[13] A. Nayak, S. Mishra, J. Hossain, and M. Nizami, "Output feedback adaptive control for inter-area oscillation damping under power system uncertainties," in Proc. IEEE Int. Conf. on Environment and Electrical Engineering and IEEE Industrial and Commercial Power Systems Europe, EEEIC/I&CPS Europe’19, 6 pp., Genova, Italy, 11-14 Jun. 2019.
[14] A. Nayak, S. Mishra, and S. Mudaliyar, "Adaptive wide area damping control for renewable integrated system," in Proc. IEEMA Engineer Infinite Conf. (eTechNxT), 6 pp., New Delhi, India, 13-14 Mar. 2018.
[15] K. M. Sreedivya, P. Aruna Jeyanthy, and D. Devaraj, "Improved design of interval type-2 fuzzy based wide area power system stabilizer for inter-area oscillation damping," Microprocessors and Microsystems, vol. 83, Article ID: 103957, Jan. 2021.
[16] I. Abdulrahman and G. Radman, "Wide-area-based adaptive neuro-fuzzy SVC controller for damping interarea oscillations," Canadian J. of Electrical and Computer Engineering, vol. 41, no. 3, pp. 133-144, Summer. 2018.
[17] M. E. Bento, "A hybrid particle swarm optimization algorithm for the wide-area damping control design," IEEE Trans. on Industrial Informatics, vol. 18, no. 1, pp. 592-599, Jan. 2022.
[18] J. Zhang, C. Chung, C. Lu, K. Men, and L. Tu, "A novel adaptive wide area PSS based on output-only modal analysis," IEEE Trans. on Power Systems, vol. 30, no. 5, pp. 2633-2642, Sept. 2014.
[19] I. Zenelis and X. Wang, "Wide-area damping control for interarea oscillations in power grids based on PMU measurements," IEEE Control Systems Letters, vol. 2, no. 4, pp. 719-724, Oct. 2018.
[20] L. Zeng, et al., "Design and real-time implementation of data-driven adaptive wide-area damping controller for back-to-back VSC-HVDC," International J. of Electrical Power & Energy Systems, vol. 109, pp. 558-574, Jul. 2019.
[21] C. Lu, Y. Zhao, K. Men, L. Tu, and Y. Han, "Wide-area power system stabiliser based on model-free adaptive control," IET Control Theory & Applications, vol. 9, no. 13, pp. 1996-2007, 27 Aug. 2015.
[22] S. Mukherjee, A. Chakrabortty, H. Bai, A. Darvishi, and B. Fardanesh, "Scalable designs for reinforcement learning-based wide-area damping control," IEEE Trans. on Smart Grid, vol. 12, no. 3, pp. 2389-2401, May 2021.
[23] I. Zenelis and X. Wang, "A model-free sparse wide-area damping controller for inter-area oscillations," International J. of Electrical Power & Energy Systems, vol. 136, Article ID: 107609, Mar. 2022.
[24] X. Zhang, C. Lu, S. Liu, and X. Wang, "A review on wide-area damping control to restrain inter-area low frequency oscillation for large-scale power systems with increasing renewable generation," Renewable and Sustainable Energy Reviews, vol. 57, pp. 45-58, May 2016.
[25] J. C. Mantzaris, A. Metsiou, and C. D. Vournas, "Analysis of inter area oscillations including governor effects and stabilizer design in South-Eastern Europe," IEEE Trans. on Power Systems, vol. 28, no. 4, pp. 4948-4956, Nov. 2013.
[26] Y. Chompoobutrgool, Concepts for Power System Small Signal Stability Analysis and Feedback Control Design Considering Synchrophasor Measurements, Diss. KTH Royal Institute of Technology, 2012.
[27] H. Bevrani, M. Watanabe, and Y. Mitani, Power System Monitoring and Control, John Wiley & Sons, 2014.
[28] P. Kundur, N. J. Balu, and M. G. Lauby, Power System Stability and Control, McGraw-Hill New York, 1994.
[29] C. Rergis, I. Kamwa, R. Khazaka, and A. R. Messina, "A Loewner interpolation method for power system identification and order reduction," IEEE Trans. on Power Systems, vol. 34, no. 3, pp. 1834-1844, May 2018.
[30] V. V. Terzija, "Adaptive underfrequency load shedding based on the magnitude of the disturbance estimation," IEEE Trans. on Power Systems, vol. 21, no. 3, pp. 1260-1266, Aug. 2006.
[31] H. Golpira, H. Bevrani, A. R. Messina, and B. Francois, "A data-driven under frequency load shedding scheme in power systems," IEEE Trans. on Power Systems, p. 1, 10.1109/TPWRS.2022.3172279, Early Access, 2022.
نشریه مهندسی برق و مهندسی كامپیوتر ایران، الف- مهندسی برق، سال 20، شماره 3، پاییز 1401 225
مقاله پژوهشی
طراحی کنترلکننده میراگر بر پایه مدل شناساییشده با استفاده
از دادههای حاصل از اندازهگیرهای فازوری حوزه وسیع
آذین عطاردی، هیمن گلپیرا و حسن بیورانی
چكیده: تغییرات مداوم و پیچیدگی سیستمهای قدرت مدرن، باعث بروز چالشهایی در زمینه مدلسازی آنها شده است. امروزه با پیشرفت سیستمهای پایش حوزه وسیع، دادههای حاصل از این سیستمها میتوانند در شناسایی و تخمین مدل سیستمهای قدرت به کار روند. این مقاله بر تنظیم پایدارساز سیستم قدرت با استفاده از مدل شناساییشده به کمک اندازهگیریهای فازوری حوزه وسیع تمرکز دارد. مدل مرتبه پایین شناساییشده بر اساس اندازهگیریها دربرگیرنده خصوصیات دینامیکی مدهای غالب بین ناحیهای بوده و میتواند برای طراحی کنترلکننده میراگر و ارزیابی اثربخشی آن در سیستم قدرت واقعی استفاده شود. کنترلکننده از نوع پایدارساز سیستم قدرت انتخاب شده و به دو روش مقاوم و زیگلر- نیکولز طراحی شده است. نتایج عددی، اثربخشی این رویکرد را در بهبود میرایی مد بین ناحیهای سیستم دوناحیهای 4ماشینه با استفاده از دادههای فازوری و دینامیکی به دست آمده از شبیهسازی در نرمافزار MATLAB نشان میدهد.
کلیدواژه: سیستمهای پایش و کنترل حوزه وسیع، کنترلکننده میراگر، مد بین ناحیهای، مدل مرتبه پایین.
1- مقدمه
امروزه دلایلی همچون بهبود قابلیت اطمینان و بهرهوری اقتصادی، توجه بهرهبرداران را به سمت اتصال سیستمهای قدرت و تشکیل سیستمهای بههمپیوسته سوق داده است. از طرف دیگر، ارتباط و بههمپیوستگی سیستمهای قدرت، باعث به وجود آمدن چالشهایی در کنترل و پایداری سیستم مانند بروز نوسانات بین ناحیهای فرکانس پایین با میرایی ضعیف میشوند. این نوسانات با میرایی ضعیف به نوبه خود باعث به خطر افتادن پایداری سیگنال کوچک سیستم میشوند. از این رو، روشهای متنوعی برای طراحی کنترلکنندههای میراگر مانند پایدارساز سیستم قدرت (PSS)، ادوات FACTS و یا مکانیزمهایی برای عملکرد هماهنگ این دو به منظور بهبود نوسانات بین ناحیهای و محلی توسعه یافتهاند [1].
میرایی نوسانات فرکانس پایین را میتوان از طریق پایدارساز سیستم قدرت و با ایجاد گشتاور میراکننده اضافی با کنترل محلی (با استفاده از سیگنالهای بازخورد محلی مانند انحراف سرعت ژنراتور) بهبود بخشید. در حالی که این رویکرد برای میراسازی مدهای محلی مناسب است، اما بهبود نوسانات بین ناحیهای همواره از طریق به کارگیری سیگنالهای محلی امکانپذیر نیست [2]. با توسعه سیستمهای پایش و کنترل حوزه وسیع (WAMS)، دادههای نقاط مختلف سیستم در دسترس قرار گرفتهاند و پایش زمان حقیقی سیستمهای قدرت بر اساس اندازهگیریهای حوزه وسیع برای تخمین حالت و حفاظت سیستم بیش از پیش مورد توجه قرار گرفتهاند. سیستم اندازهگیری حوزه وسیع مبتنی بر سنکروفازور در شبکههای قدرت، رؤیتپذیری دینامیکهای سیستم قدرت را به صورت چشمگیری افزایش میدهد و شرایط مناسبی را برای مشاهده، شناسایی و میراسازی نوسانات بین ناحیهای فراهم میکند. مطالعات نشان میدهند که مشخصات دینامیکی مدهای نوسانی، به ویژه مدهای بین ناحیهای فرکانس پایین با میرایی ضعیف، با اندازهگیریهای فازوری قابل شناسایی است [3]. مدهای ناپایدار شناساییشده میتوانند به صورت تطبیقی با تنظیم پارامترهای کنترلی سیستم قدرت بر اساس دادههای اندازهگیری شده میرا شوند.
به دلیل تغییرات پیوسته در سیستمهای قدرت پیشرفته و پیچیدگی و توسعه آنها، مدلسازی سیستمهای قدرت به امری دشوار تبدیل شده است. سیستمهای اندازهگیری حوزه وسیع میتوانند در این زمینه خاص نیز راهگشا باشند. در سالهای اخیر، روشهای زیادی برای تخمین مدهای نوسانی سیستم با استفاده از دادههای اندازهگیری حوزه وسیع توسعه یافتهاند [4]. این روشها مدل دینامیکی سیستم را با استفاده از دادههای اندازهگیری شده و به دو فرم فضای حالت فرعی [5] و تابع تبدیل [6] شناسایی میکنند. در مقایسه با مدل فضای حالت فرعی، مدل تابع تبدیل چندورودی- چندخروجی کارایی محاسباتی بیشتری دارد.
با داشتن مدل دینامیکی سیستم و تعیین مدهای نوسانی با میرایی ضعیف، میتوان به سوی طراحی کنترلکننده میراساز مناسب گام برداشت. در این راستا در سالهای اخیر، محققان تکنیکهای کنترلی مختلفی را معرفی کردهاند. در [7] یک کنترلکننده تطبیقی بر اساس الگوریتم حداقل مربعات بازگشتی برای شبکه قدرت با نفوذ بالای منابع تجدیدپذیر پیشنهاد شده است. با وقوع اغتشاش در سیستم، یک تابع تبدیل چندورودی- چندخروجی با استفاده از دادههای حاصل از اندازهگیری شناسایی شده و سپس کنترلکننده طراحیشده با کنترلکنندههای محلی منابع تجدیدپذیر ادغام شده است. در [8] برای دستیابی به روشی سریع و برخط، مدل چندورودی- چندخروجی شناساییشده با روش مانده به یک مدل تکورودی- تکخروجی تبدیل شده است. این مدل میتواند با دقت مناسبی، دینامیکهای سیستم قدرت را به صورت تابع تبدیلی ارائه کند که دربرگیرنده رفتارهای نوسانی غالب در مد مورد نظر است.
رویکردهای بسیاری برای طراحی کنترلکننده مبتنی بر اندازهگیری از جمله تکنیک شناسایی مرتبه پایین به همراه روش مقاوم تنظیم PSS بر اساس الگوریتم ژنتیک [9] و شناسایی مدل دوورودی- دوخروجی برای ژنراتور با بیشترین مشارکت در مد و طراحی کنترلکننده مقاوم به روش H∞ [10] بررسی شدهاند. در [11] مسئله تنظیم PSS كه به عنوان يك مشكل بهینهسازی چندهدفه فرموله شده است، با استفاده از الگوريتم بهینهسازی ژنتيك حل شده است. روش پيشنهادي به شبكه چندماشينه اعمال گردیده و شبیهسازیهای غیر خطی براي تحليل پايداري سيستم قدرت انجام شده است. نتايج شبیهسازی حوزه زماني نيز نشان میدهد كه روش ارائهشده توانسته تا ميرايي قابل قبولی را فراهم آورد. در [12] با استفاده از روش حساسیت مرکب H∞ و با دو هدف حذف اثر اغتشاش و بهینهسازی اثر کنترلی، کنترلکننده حوزه وسیع مقاوم جهت اعمال سیگنال اضافی پایدارساز به SVC برای میراکردن نوسانات بین ناحیهای طراحی شده و از روش نامعادلات ماتریسی خطی (LMI) برای حل مسئله استفاده شده است.
همچنین تکنیکهای تطبیقی متعددی مانند کنترل تطبیقی مدل مرجع
[13]، کنترل خودتنظیم [14]، منطق فازی [15]، کنترل عصبی- فازی [16] و الگوریتمهای بهینهسازی مانند تجمیع ذرات [17] در طراحی کنترلکنندههای میرایی حوزه وسیع مطرح شدهاند. در روش خودتنظیم، یک مدل سادهشده خطی به صورت برخط برای سیستم قدرت شناسایی میشود و تنظیم پارامترهای کنترلکننده بر اساس مدل شناساییشده با استفاده از روشهایی مانند جایابی قطب صورت میگیرد. در [18]، فاز مانده مد بین ناحیهای مورد نظر نسبت به حلقه کنترل به صورت برخط از طریق تحلیل سیگنالهای محیطی مشخص شده و سپس کنترلکننده تطبیقی به کمک روش ماندهها تنظیم شده است.
در دهه اخیر، روشهای کنترل دادهمحور برای سیستمهای پیچیده غیر خطی پیشنهاد شدهاند تا مشکلات ناشی از تکنیکهای مدلمحور را حل کنند [19] تا [23]. در [21] یک روش کنترل دادهمحور مجزا از مدل معرفی شده که با در نظر گرفتن اغتشاشات سیستم، میرایی سیستم را بهبود میبخشد. در [22] کنترل حوزه وسیع مجزا از مدل و با استفاده از روش آموزش تقویتی (RL) طراحی شده است. همچنین یک روش بازخورد خروجی با استفاده از اندازهگیری فرکانس ژنراتورها پیشنهاد شده که در برابر تغییرات قابل توجه بار به خوبی پاسخ میدهد. در [23] یک کنترلکننده برخط بازخورد حالت بهینه مبتنی بر تنظیمکننده خطی درجه دوم (MLQR) پیشنهاد شده که میتواند چندین مد بین ناحیهای نامطلوب را به طور همزمان با حداقل تلاش کنترلی خنثی کند. کارایی رویکرد پیشنهادی علیرغم نویز اندازهگیری، تلفات PMU و محدودیتهای شبکه ارتباطی تأیید شده است.
بیشتر این روشها علیرغم این که وابستگی به مدل سیستم ندارند،
به اغتشاشات بزرگ برای تعیین پارامترهای کنترلکننده بستگی دارند و فرایند آموزش ممکن است با مشکل همگرایی مواجه شود. در زمینه تحلیل مد نوسانی، رویکردهای مبتنی بر تجزیه و تحلیل پاسخ خروجی سیستم پس از اغتشاش از جمله تبدیل موجک، تبدیل هیلبرت- هوانگ، تجزیه و تحلیل پرونی و فیلتر کالمن پرکاربرد هستند [23].
تغییرات تصادفی و پیوسته بار سیستم قدرت در شرایط عادی منبع اصلی تحریک دینامیکهای الکترومکانیکی هستند. این محرکها به نویزهای محیطی در سیگنالها و دادههای اندازهگیری شده معروف هستند. نظارت و تحلیل پایداری سیستم قدرت با استفاده از سیگنالهای محیطی ثبتشده از طریق WAMS قابل انجام است. تحلیل سیگنال محیطی به حوزه پردازش سیگنال تصادفی مرتبط بوده و شامل روشهای شناسایی زیرفضای تصادفی، روشهای سریهای زمانی (مدل ARMA)، تجزیه حوزه فرکانس، حداقل مربعات بازگشتی، حداقل مربعات میانگین و ... است [24].
از دیگر موارد مورد بحث در زمینه کنترل حوزه وسیع، انتخاب سیگنال ورودی است که عمدتاً طبق روشهای مبتنی بر شاخصهای کنترلپذیری و مشاهدهپذیری صورت میگیرد. به عنوان مثال در [25] سیگنال انحراف سرعت به عنوان سیگنال ورودی به کنترلکننده انتخاب شده است. بخشی از نتایج [26] نشان میدهند که با توجه به میرایی و بهره کنترلکننده مورد نظر، اختلاف زوایای ولتاژ در مقایسه با مقدار ولتاژ و یا سرعت ژنراتور میتوانند سیگنالهای بازخورد کاراتری باشند.
در حالی که اکثر روشهای شناسایی ذکرشده به اطلاعات سیگنالهای ورودی و خروجی سیستم نیاز دارند، مزیت روش به کار رفته در این مقاله آن است که فرایند شناسایی سیستم نیاز به اطلاعات مربوط به ورودی اغتشاش ندارد. به عبارت دیگر، روش پیشنهادی صرفاً نیازمند اطلاعات فازوری ولتاژ حالت ماندگار است. مزیت دیگر آن است که در فرایند آمادهسازی سیگنال مناسب برای شناسایی مدل از روش ساده برازش منحنی استفاده میشود. مدل به دست آمده، یک مدل مرتبه پایین بوده که برای طراحی کنترلکننده بسیار مناسب است. در مرحله بعد، پایداری مدهای نوسانی غالب مستقیماً با استفاده از مقادیر ویژه مدل شناساییشده بررسی میشود. پس از تخصیص یک مدل مرتبه پایین تقریبی به سیستم، پایدارساز سیستم قدرت به دو روش H∞ و زیگلر- نیکولز برای میراکردن مد بین ناحیهای غالب شناساییشده طراحی و تنظیم میشود.
2- مدلسازی مرتبه پایین سیستم
با استفاده از دادههای فازوری
این بخش به کمک اندازهگیریهای فازوری حوزه وسیع، یک روش شناسایی سیستم قدرت را معرفی میکند. به عبارت دیگر یک مدل مرتبه پایین- که دربرگیرنده مدهای نوسانی غالب است- که در طراحی کنترلکننده مورد استفاده قرار میگیرد، به سیستم اختصاص داده میشود.
با فرض سیستم تکماشینه متصل به شین بینهایت، نوسانات دینامیکی بین ناحیهای میتواند توسط پارامترهای زاویه فاز و سرعت زاویهای
به صورت زیر بیان شود [27]
(1)
که در آن و
اختلاف زاویه فاز و سرعت زاویهای نواحی سیستم بوده و مقادیر آنها با استفاده از فازور ولتاژ به دست آمده از دادههای PMU در دسترس است. حال چنانچه
را تفاضل سرعت زاویهای دو ناحیه انتخابی
و
و
را تفاضل زاویه فاز دو ناحیه در نظر بگیریم، آن گاه خواهیم داشت
(2)
در (2) زیروند نشانگر زاویه فاز اولیه نواحی
و
در زمان
است. با مشتقگیری از عبارات فوق خواهیم داشت
(3)
شكل 1: مدل مرتبه پایین تقریبی برای تنظیم پایدارساز سیستم قدرت جهت میراکردن مد نوسانی بین ناحیهای [27].
شكل 2: سیستم قدرت دوناحیهای.
بنابراین، یک مدل نوسان مرتبه دوم برای مد نوسانی غالب میتواند به شکل زیر تعریف شود
(4)
حال با در دست داشتن دادههای زاویه فاز و سرعت زاویهای گرفتهشده از PMUها و تشکیل بردارهای و
و مشتقات آنها، ضرایب
و
را میتوان به روش حداقل مربعات خطا مشخص نمود. ویژگیهای مد استخراجشده میتواند با مقادیر ویژه ماتریس ضرایب
ارزیابی شود. با فرض مقادیر ویژه به صورت
، مدل نوسانی مرتبه پایین تقریبی برای مد مورد نظر به صورت زیر به سیستم اختصاص داده میشود
(5)
در گام بعدی هدف، طراحی کنترلکنندهای برای میراکردن مد بین ناحیهای نوسانی موجود در مدل شناساییشده در (5) است. با در نظر گرفتن یک کنترلکننده میراگر، حلقه کنترلی با مدل نوسانی
را میتوان با استفاده از شکل 1 در نظر گرفت. کنترلکننده میراگر که در اینجا پایدارساز سیستم قدرت انتخاب شده است، شامل ساختار اصلی کنترلکننده و مدل
- که بیانگر تأثیر کنترلکننده روی مد نوسانی
است- خواهد بود. برای سادهترکردن روند طراحی کنترلکننده، مرتبه
باید تا حد ممکن کاهش داده شود. در این مقاله، مدل
به شکل زیر در نظر گرفته میشود [27]
(6)
با توجه به توضیحات فوق، فرایند تعیین به صورت زیر است:
گام 1) تعیین ثابتهای زمانی تا
به نحوی که بهره حول فرکانس مورد نظر که همان بخش موهومی مقدار ویژه تخمین زده شده است، حداکثر شود.
گام 2) تعیین ثابت زمانی در حدود صفر درجه پاسخ فرکانس
حول فرکانس مورد نظر
.
شكل 3: بلوک دیاگرام سیستم تحریک و PSS ژنراتور 1 [27].
شكل 4: بلوک دیاگرام سیستم تحریک و PSS ژنراتورهای 2 تا 4 [27].
گام 3) محاسبه بهره در راستای حداقلکردن مربعات خطا بین اختلاف فاز فیلترشده و خروجی مدل تقریبی.
3- شناسایی سیستم و طراحی کنترلکننده
به منظور درک بهتر رویکرد ارائهشده، روش پیشنهادی بر روی سیستم قدرت دوناحیهای شکل 2 اعمال شده است. اطلاعات مربوط به این سیستم شامل اطلاعات ژنراتورها، بارها، خطوط و باسها در [28] آمده است. ژنراتور 1 با AVR و PSS نشان داده شده در شکل 3 برای میراکردن نوسانات بین ناحیهای مجهز شده است. سایر ژنراتورها با AVR و PSS یکسان نشان داده شده در شکل 4 مجهز شدهاند.
3-1 شناسایی سیستم
برای نشاندادن کارایی روش پیشنهادی، در گام اول رفتار سیستم بدون در نظر گرفتن PSS برای ژنراتور 1 بررسی میشود. از آنجایی که به منظور شناسایی مدل نیازی به اعمال اغتشاش وجود ندارد، تغییرات کوچک بار در باس 4 برای شبیهسازی نوسانات فازورهای اندازهگیری شده در شرایط عادی سیستم قدرت در نظر گرفته شده است. با توجه به آنچه در بخش مقدمه آمده است، فازور ولتاژ اطلاعات مناسبی برای شناسایی سیستم به دست میدهد. با توجه به این نکته و همچنین نیاز به جمعآوری داده از دو ناحیه سیستم، دو PMU در باسهای 1 و 11 برای اندازهگیری تغییرات فازور ولتاژ در دو ناحیه در نظر گرفته شده است. شکل 5 اختلاف بین زاویه فاز ولتاژ باسهای 1 و 11 را نشان میدهد. تغییرات ناگهانی دیدهشده در این شکل به دلیل کلیدزنیهای حین شبیهسازی به منظور افزایش یا کاهش بار به وجود آمده است. در شکل 6 طیف فرکانسی موجود در شکل 5 نشان داده شده است. نتایج نشان میدهند که سیستم مورد مطالعه دارای مد نوسانی بین ناحیهای غالب با فرکانس حدود 6/0 هرتز است.
با استفاده از فیلتری با فرکانس مرکزی 6/0 هرتز با پهنای باند مناسب و سپس تبدیل معکوس فوریه، دادههای اختلاف فاز فیلترشده به دست میآیند. سپس میتوان با اعمال روش حداقل مربعات بر روی دادههای اختلاف فاز فیلترشده، یک مدل مرتبه پایین به سیستم اختصاص داد. در این مقاله به جای استفاده از فیلتر بر پایه تبدیل فوریه، از برازش منحنی دادههای اختلاف فاز استفاده شده است. به این ترتیب که با در نظر گرفتن
شکل 5: اختلاف فاز بین ولتاژ باسهای 1 و 11.
شکل 6: طیف فرکانسی حاصل از تبدیل فوریه.
یک پنجره زمانی، اختلاف فاز ولتاژهای دو شین مورد نظر به شکل سیگنالی از مجموع چند تابع سینوسی به صورت زیر تقریب زده شده است
(7)
همان طور که ملاحظه میشود، فرکانس مد مورد نظر (حدود 6/0 هرتز برابر 7/3 رادیان بر ثانیه) در جملات دوم و چهارم ظاهر شده است. در ادامه، از تابع حاصل مشتقات متوالی گرفته شده تا دادههای اختلاف سرعت و مشتق آن به دست آیند. سپس این دادهها در (4) قرار داده شده و ماتریس ضرایب به کمک روش حداقل مربعات محاسبه شده است. بر این اساس ماتریس حالت به دست آمده به قرار زیر است
(8)
مقادیر ویژه این ماتریس برابر با است. با جایگذاری مقادیر ویژه فوق در (5)، مدل نوسان به صورت زیر به دست آمده است
(9)
با بهرهگیری از گامهای 1 تا 3 در بخش 2 این مقاله مدل به صورت زیر مشخص میشود
شکل 7: دیاگرام بود مدل .
جدول 1: مقایسه مدهای نوسانی مدل استخراجشده و [29].
درصد خطا | ضریب میرایی (%) | درصد خطا | فرکانس مد غالب | روش شناسایی |
5 | 94/0 | 0 | 6/0 | روش پیشنهادی |
0 | 99/0 | 4/3 | 58/0 | [29] |
(10)
شکل 7 که دیاگرام بود (Bode) مدلدر (10) را نمایش میدهد، بیانگر این نکته است که پیک بهره در فرکانس 77/3 رادیان بر ثانیه اتفاق می افتد.. در نهایت جدول 1 به ارزیابی کارایی روش شناسایی ارائهشده از طریق مقایسه با روش ارائهشده در [29] میپردازد.
3-2 طراحی کنترلکننده
همان گونه که قبلاً اشاره شد، سیستم قدرت مدرن روزبهروز در حال تغییر است. به همین دلیل، حضور عدم قطعیتها مانند تغییرات پارامتری در فرایند مدلسازی غیر قابل اجتناب است. روشهای مقاوم و تطبیقی، دو رویکرد رایج برای مواجهه با نامعینیهای سیستم هستند. در این مقاله، با در نظر گرفتن سیگنال PSS ژنراتور شماره 1 به عنوان ورودی کنترل، طراحی و تنظیم یک کنترلکننده به روش ∞H و یک کنترلکننده PID به روش زیگلر- نیکولز برای سیستم شناساییشده جهت میراکردن مد بین ناحیهای غالب مورد بررسی قرار گرفته است.
3-2-1 کنترلکننده مقاوم
در [27] پایدارساز سیستم قدرت به روش مقاوم و با تئوری ∞H
برای مدل مرتبه پایین طراحی شده است. بدین منظور، مدل تعمیمیافته
به صورت شکل 8 در نظر گرفته شده و دو تابع وزنی به شرح زیر تعریف شدهاند
(11)
و
(12)
شکل 8: ساختار سیستم حلقه بسته جهت طراحی کنترلکننده مقاوم.
شکل 9: نمودار بلوکی روش زیگلر- نیکولز مبتنی بر فیدبک رلهای.
جدول 2: ضرایب کنترلکننده ∞H.
3161 |
| 106 × 39/1 |
|
105 × 88/4 |
| 108 × 64/1 |
|
107 × 02/2 |
| 109 × 03/4 |
|
108 × 85/1 |
| 1010 × 03/3 |
|
108 × 93/5 |
| 1010 × 94/8 |
|
108 × 10/7 |
| 1010 × 46/9 |
|
107 × 06/2 |
| 1010 × 20/1 |
|
107 × 40/1 |
| 108 × 43/3 |
|
در معادلات فوق به صورت یک فیلتر پایینگذر برای لحاظکردن خطای مدلسازی و پایداری مقاوم و
به صورت یک فیلتر بالاگذر برای لحاظکردن تأثیر سایر مدها در نظر گرفته شده است.
نرم بینهایت سیستم به صورت زیر قابل تعریف است
(13)
که در آن تابع تبدیل کلی سیستم شامل
و
بوده و
تابع تبدیل کنترلکننده است. ورودی
برای مدل تقریبی از حاصلضرب معکوس (9) و دادههای اختلاف فاز تخمین زده میشود. کنترلکننده ∞H برای مدل مرتبه پایین شناساییشده مربوط به حداقل
که قید (13) را برآورده میکند بر اساس مدل نشان داده شده در شکل 8 به دست آمده است. کنترلکننده حاصل از مرتبه 7 و به صورت زیر بوده و ضرایب آن در جدول 2 آمده است
(14)
پس از 4 درجه کاهش مرتبه، تابع تبدیل پایدارساز سیستم قدرت عبارت است از
(15)
شکل 10: ورودی و خروجی سیستم با فیدبک رلهای.
جدول 3: تنظیم پارامترهای کنترلکننده PID با استفاده
از روش زیگلر- نیکولز حلقه بسته [30].
بهره بخش مشتقگیر | بهره بخش انتگرالی | بهره بخش تناسبی | نوع کنترلکننده |
- | - |
| P |
- |
|
| PI |
|
|
| PID |
3-2-2 کنترلکننده PID
در این مقاله، جهت طراحی کنترلکننده PID از روش زیگلر- نیکولز مبتنی بر فیدبک رلهای استفاده شده است. ایده اساسی، این نگرش است که بسیاری از فرایندها تحت فیدبک رلهای دارای نوسانات سیکل حدی هستند. در این روش با استفاده از یک فیدبک رلهای و ورودی سینوسی، سیستم به نوسان درآمده و سپس دو پارامتر بهره و دوره تناوب بحرانی محاسبه میشوند [30]. شکل 9 نمودار بلوکی مربوط به این روش را نشان میدهد.
اگر دامنه خروجی رله با و دامنه خروجی فرایند با
نمایش داده شوند، بهره بحرانی
از رابطه زیر به دست میآید
(16)
با در نظر گرفتن به عنوان دوره تناوب خروجی سیستم، برای تنظیم پارامترهای کنترلکننده میتوان از جدول 3 استفاده کرد.
با اعمال روش فوق روی مدل مرتبه پایین و ورودی سینوسی با دامنه واحد و فرکانس 60 هرتز، شکل 10 خروجی
را نشان میدهد. پس از اندازهگیری دامنه و دوره تناوب خروجی سیستم، کنترلکننده PID به فرم زیر حاصل میشود
(17)
با در نظر گرفتن حلقه کنترلی شکل 1، تابع تبدیل مربوط به بخش پایدارساز سیستم قدرت از رابطه زیر قابل محاسبه است
(18)
تابع تبدیل حاصل از مرتبه 5 بوده و پس از دو درجه کاهش مرتبه عبارت است از
شكل 11: پاسخ پله مدل تقریبی مرتبه پایین.
شكل 12: توان عبوری از باس 13 به 101.
[1] این مقاله در تاریخ 6 مهر ماه 1400 دریافت و در تاریخ 18 اردیبهشت ماه 1401 بازنگری شد.
آذین عطاردی، دانشکده مهندسی، گروه مهندسی برق، دانشگاه کردستان، کردستان، ایران، (email: a.atarodi@uok.ac.ir).
هیمن گلپیرا (نویسنده مسئول)، دانشکده مهندسی، گروه مهندسی برق، دانشگاه کردستان، کردستان، ایران، (email: hemin.golpira@uok.ac.ir).
حسن بیورانی، دانشکده مهندسی، گروه مهندسی برق، دانشگاه کردستان، کردستان، ایران، (email: bevrani@uok.ac.ir).
جدول 4: مدهای نوسانی سیستم دوناحیهای.
مد محلی ناحیه 2 | مد محلی ناحیه 1 | مد بین ناحیهای | مد PSS |
|