Passive Islanding Detection Using Adaptive Threshold Based on Instantaneous Frequency Droop

Document Type : Research Article

Authors

1 Power Department, Faculty of Electrical Engineering, Shahrood University of Technology, Shahrood, Iran

2 Associate Professor, Power Department, Faculty of Electrical Engineering, Shahrood University of Technology, Shahrood, Iran

3 Assistant Professor, Power Department, Faculty of Electrical Engineering, Shahrood University of Technology, Shahrood, Iran

Abstract

Islanding is one of the important challenges in power networks in the presence of distributed generation which is considered an undesirable incident due to the possibility of damage to operators, network equipment and consumers. Therefore, it is necessary to quickly detect islanding and make a decision regarding the connection status of local distributed generation units in the network. In this paper, a passive islanding detection method is proposed using adaptive threshold-based instantaneous frequency droop characteristic. The proposed threshold limit is dynamically changed depending on under studied conditions in order to discriminate the island from load changes, capacitor bank switching, motor start-up and short circuit types. Two medium voltage networks of Cigre and 34-bus IEEE are used to evaluate the proposed method. The simulations are performed in Digsilent software and implementation of the proposed method carried out in MATLAB software. The simulation results indicate that the island is correctly detected with appropriate accuracy in a short period of time from other disturbances in the presence of wind, solar and diesel types of distributed generations.

Keywords

Main Subjects


  • مقدمه[1]

در دهه‌های اخیر مسائلی مانند افزایش قیمت سوخت‌های فسیلی و نگرانی‌های زیست محیطی ناشی از آن در نیروگاه‌های حرارتی، نیاز به منابع تجدیدپذیر[1] (RES) را افزایش داده است. همچنین افزایش مشارکت بخش خصوصی در تولید انرژی، مزایای زیست‌محیطی و کاهش تلفات، صنعت برق را به سمت استفاده از تولیدات پراکنده[2] (DG) به خصوص منابع انرژی تجدیدپذیر سوق داده است [1].

نیروگاه‌های تولید پراکنده از جمله نیروگاه بادی[3]، میکروتوربین[4]، خورشیدی[5] و CHP [6]، واحدهای کوچک نیروگاهی هستند که می‌توانند نزدیک به بار و بدون هزینه‌های انتقال، نصب شوند [2]. حضور تولیدات پراکنده در شبکه‌ قدرت دارای مزایای زیادی از جمله بهبود کارایی شبکه، افزایش انعطاف‌پذیری شبکه، کاهش تلفات و کاهش آلودگی محیط‌ زیست است؛ با این‌حال، افزایش تولیدات پراکنده موجب نگرانی در هماهنگی حفاظت، امنیت و پایداری شبکه، تنظیم ولتاژ و نیز جزیره‌ای شدن شبکه شده است [3].

وقوع جزیره به صورت عمدی و غیرعمدی در شبکه رخ می‌دهد. جزیره عمدی به منظور بالا بردن قابلیت اطمینان بارهای خاص در زمان قطع شبکه اصلی اتفاق می‌افتد. جزیره غیرعمدی به قطع ناگهانی شبکه اصلی، در حالی‌که انرژی قسمتی از شبکه توسط DG تأمین می‌گردد، گفته می‌شود. این نوع جزیره به دلیل داشتن خطراتی مانند تغییرات فرکانس و ولتاژ، احتمال آسیب به بهره‌برداران شبکه و از بین رفتن حالت سنکرون شبکه، یک وضعیت نامطلوب تلقی شده، لذا تشخیص جزیره غیرعمد ضروری است [4, 5].

کارایی و سودمندی یک روش تشخیص از زوایای مختلف مورد ارزیابی قرار می‌گیرد. اصول اولیه مقایسه سرعت و دقت در تشخیص است؛ با این‌حال کوچک بودن ناحیه غیرقابل تشخیص[7] (NDZ)، تشخیص در ضریب‌های مختلف کیفیت بار، عملکرد در حضور چندین واحد تولید پراکنده و عدم کاهش کیفیت توان ازجمله مواردی است که باعث ارجحیت یک روش می‌شود. به طور کلی روش‌های تشخیص جزیره‌[8] (IDM) به دو دسته محلی[9] و راه دور[10] تقسیم‌بندی می‌شوند [6, 7]. روش راه دور با استفاده از ارتباط بین پست، ژنراتور، تجهیزات و DG ها، جزیره را تشخیص می‌دهد. روش راه دور دارای قابلیت اطمینان بالا با ناحیه NDZ صفر است؛ اما به خاطر هزینه و پیچیدگی بالای این روش، IDM های محلی برای مهندسین جذاب‌تر می‌باشند [3]. روش ارائه شده در [8] از تغییرات جریان کلیدهای قدرت، فرکانس و توان اکتیو مدار در تمام شاخه‌ها برای تشخیص جزیره استفاده می‌کند. زمان تشخیص جزیره بسیار کوتاه و ناحیه NDZ صفر است؛ ولی هزینه اجرای این طرح بالاست و همچنین نیاز به اندازه‌گیری مقادیر محلی می‌باشد. روش‌های محلی از پارامترهای اندازه‌گیری‌ شده شبکه در محل DG کمک می‌گیرد. روش‌های محلی به روش‌های فعال[11]، غیرفعال[12] و مبتنی بر ریاضیات دسته‌بندی می‌شوند. IDM های فعال با تزریق یک اختلال در سیستم، جزیره را تشخیص می‌دهند که دارای قابلیت تشخیص جزیره در شرایط تعادل توان می‌باشند [9]. در روش‌های فعال با استفاده از سیگنال­های تزریقی به شبکه، در زمان وقوع جزیره پارامترهای سیستم (مانند ولتاژ و فرکانس) دچار تغییر چشمگیری می­شوند که این تغییرات از وقوع جزیره خبر می‌دهد. این روش‌ها، دارای ناحیه NDZ کوچک و حتی صفر می‌باشند و سبب بهبود سرعت تشخیص می‌گردند ولی ناخواسته موجب کاهش کیفیت توان می‌شوند. با این حال ایراد اصلی روش­های فعال، عملکرد ضعیف آن­ها در حضور منابع تولید پراکنده اینورتری است [10]. استفاده از تزریق هارمونیک به شبکه جهت تشخیص جزیره روشی است که در [9] از آن استفاده شده است. در این مقاله برخلاف سایر روش­های فعال، کیفیت توان بهبود پیدا کرده و پایداری شبکه حفظ می‌شود. هارمونیک دوم جریان با توجه به فیدبک هارمونیک دوم ولتاژ در نقطه اتصال مشترک[13] (PCC) به شبکه تزریق می‌شود. شناسایی جزیره بوسیله یک روش فعال برای نیروگاه‌های خورشیدی مبتنی بر کنترل‌کننده لغزشی در [11] مورد مطالعه قرار گرفته است. دقت این روش برای در نیروگاه‌های خورشیدی تکفاز مناسب بوده ولی در DG های غیراینورتری کاربردی ندارد. در [12] از اعوجاج توان راکتیو برای تشخیص جزیره استفاده می‌شود. این روش برای DG های اینورتری با ضریب توان غیر یک، مناسب می‌باشد.

در روش غیرفعال بر اساس پارامترهای اندازه­گیری شده از قبیل ولتاژ و فرکانس، اعوجاج هارمونیکی در نقطه PCC، جزیره تشخیص داده می‌شود. روش­های مبتنی بر کاهش یا افزایش ولتاژ، کاهش یا افزایش فرکانس، نرخ تغییر فرکانس، نرخ تغییر فرکانس با توان، کنترل هارمونیک در دسته روش‌های غیرفعال قرار می­گیرند. معمولاً از یک یا چند شاخص با حساسیت کافی برای تشخیص جزیره در روش‌های غیرفعال، استفاده می‌شود. سادگی و هزینه کم پیاده‌سازی و حفظ کیفیت توان از جمله مزایای روش غیرفعال است؛ با این حال ایراد اصلی این روش وجود ناحیه NDZ بزرگ‌تر نسبت به دیگر روش‌ها می‌باشد [9, 13]. در [14] به بررسی نسبت پارامترهای ولتاژ، فرکانس، توان اکتیو و توان راکتیو بر یکدیگر و استفاده از آنها در تشخیص جزیره می‌پردازد. نتایج مطالعه مقاله مذکور از 16 پارامتر بدست آمده حاکی از آن است که شاخص df/dq معیار مناسبی برای تشخیص جزیره با سرعت خوب و ناحیه NDZ کوچک می‌باشد.

علاوه بر روش‌های فوق از روش‌های ریاضی و هوش مصنوعی در تشخیص جزیره استفاده می‌شود. در [15] از شبکه عصبی احتمالاتی و تبدیل موجک برای تشخیص جزیره استفاده کرده است. با وجود دقت بالای روش در تشخیص جزیره این روش بر روی واحدهای تولید پراکنده اینورتری متمرکز است. در [16] با استفاده از یک دسته‌بندی کننده شبکه عصبی آموزش دیده و از طریق تبدیل موجک گسسته سیگنال‌های گذرای ولتاژ و جریان به تشخیص رخدادهای گذرا و جزیره‌ای شدن می‌پردازد. این روش در حالی که دقت تشخیص بالایی دارد ولی بر روی واحدهای اینورتری بررسی نشده است.

در این مقاله، یک روش غیرفعال بر اساس مشخصه افتی فرکانسی لحظه‌ای برای تشخیص جزیره از تغییرات بار، کلیدزنی خازنی، راه‌اندازی موتور و نیز وقوع انواع اتصال کوتاه‌های تکفاز، دوفاز، دوفاز به زمین و سه فاز پیشنهاد می‌گردد. در روش پیشنهادی حد آستانه تطبیق‌پذیر به صورت دینامیکی مشخیص می‌شود که باعث رفع مشکل روش‌های غیرفعال در تعیین مقدار آستانه می‌گردد. همچنین روش پیشنهادی در حضور انواع منابع تولیدات پراکنده دارای کارایی مطلوب است. نتایج به دست آمده نشان می‌دهد که جزیره با سرعت و دقت مطلوب از سایر حوادث به کمک روش پیشنهادی تشخیص داده می‌شود.

در ادامه مقاله، در بخش دوم مقاله، مشخصه افتی بیان می‌شود. در بخش سوم روش پیشنهادی به همراه فرمول‌ها و روند آن گفته می‌شود. در بخش چهارم شبیه‌سازی بر روی دو شبکه آزمون و وضعیت‌های متفاوت انجام و تحلیل و بررسی شده و درنهایت در بخش پنجم، نتیجه‌گیری شده است.

 

2- مشخصه افتی فرکانسی

هماهنگی تولیدات پراکنده برای پاسخگویی به بار ضروری است، تاکنون روش‌های مختلفی برای این کار از جمله کنترل افتی[14]، کنترل رهبر-پیرو[15] و کنترل میانگین جریان اشتراکی[16]، به منظور فعالیت تولیدات پراکنده به صورت موازی با شبکه جهت تأمین بار پیشنهاد شده است. در میان روش‌های مذکور، روش مشخصه افتی به دلیل عدم نیاز به شبکه‌های ارتباطی به‌طور گسترده مورد استفاده قرار گرفته است [17]. در حالت غیرجزیره، مقدار فرکانس توسط شبکه قدرت کنترل می‌شود تا در محدوده مجاز قرار بگیرد. طبق روش کنترل افتی فرکانسی، میزان مشارکت توان واحدهای تولید کننده در حضور شبکه قدرت به صورت شکل (1)-الف مشخص می‌شود [18]. در این جا ، ،  و  به ترتیب فرکانس نامی شبکه، توان اکتیو تولیدی منبع تولید پراکنده، توان اکتیو شبکه اصلی و توان اکتیو مصرفی بار محلی می‌باشد.

 

 

(الف)

 

(ب)

شکل (1): (الف) مشخصه افتی P-f در حالت غیرجزیره‌ای؛ (ب) مشخصه افتی P-f در حالت جزیره‌ای [18]

 

در صورت بروز خطا در حالت غیرجزیره، مشخصه افتی فرکانسی به خاطر اتصال منبع تولید پراکنده به شبکه ثابت بوده و بعد از رفع خطا، شبکه به وضعیت عادی خود برمی‌گردد. در صورت وقوع جزیره، وظیفه تنظیم فرکانس بر عهده منابع تولید پراکنده است؛ به عنوان مثال شکل (1)-ب چگونگی تنظیم فرکانس را در حضور دو منبع تولید پراکنده نشان می‌دهد. قبل از وقوع جزیره، منابع تولید پراکنده توان اکتیو  و  به شبکه تزریق می‌نمایند، در حالی که پس از وقوع جزیره، توان منابع برای حفظ فرکانس بسته به مشخصه افتی P-f به  و  تغییر می‌کند.

در یک ریزشبکه می‌توان از کنترل افتی فرکانسی برای کنترل و شارش توان به مانند یک شبکه اصلی بهره برد [19]. توان ظاهری تزریقی به شین PCC طبق شکل (2) به کمک رابطه (1) به دست می‌آید:

(1)

 

 

در این رابطه، E، ، Z، θ و V به ترتیب ولتاژ DG، زاویه ولتاژ DG، امپدانس خط، زاویه امپدانس خط و ولتاژ شین PCC را نشان می‌دهد. لازم به ذکر است که امپدانس خروجی منابع تولید پراکنده اینورتری و سنکرون به دلیل استفاده از سلف، عمدتاً القایی است [20]؛ بر این اساس، زاویه فاز θ نزدیک به 90 درجه است با فرض کوچک بودن δ، مقدار توان اکتیو با کمک رابطه‌ (2) به دست می‌آید [21]:

(2)

 

 

رابطه‌ (2) بر اساس مشخصه افتی P-f می‌تواند به صورت زیر بازنویسی شود:

(3)

 

 

که در آن  و  به ترتیب فرکانس نامی شبکه و ضریب افتی توان اکتیو می‌باشند. همچنین  را می‌توان به‌صورت زیر به دست آورد:

(4)

 

 

که در آن Δf حداکثر انحراف مجاز فرکانس و ΔP حداکثر انحراف مجاز توان اکتیو واحد تولید پراکنده‌ هستند.

 

شکل (2): مدار معادل یک واحد تولید پراکنده متصل به یک شین مشترک AC

 

3- معرفی روش پیشنهادی مبتنی بر مشخصه افتی لحظه‌ای

شرایط عدم تطابق توان عمدتاً در ریزشبکه با حضور منابع پراکنده مشاهده می‌گردد. در زمان اتصال ریزشبکه به شبکه اصلی اختلاف توان توسط شبکه تأمین شده، ولی در زمان جزیره این توان توسط تولید پراکنده تأمین می‌شود. در صورت عدم تطابق توان زیاد امکان دارد مقادیر فرکانس و ولتاژ از حد نامی تجاوز کرده و باعث بروز عکس‌العمل توسط رله‌های ولتاژی و فرکانسی شود [22].

 

شکل (3): شبکه قدرت در حضور تولید پراکنده و با RLC موازی [14]

 

مطابق شکل (3) در شرایط بهره‌برداری عادی تأمین بار RLC مطابق روابط (5) و (6) است که R مقدار مقاومت بار، ΔR تغییرات مقاومت پس از جزیره، L مقدار اندوکتانس بار، ΔL تغییرات اندوکتانس پس از جزیره، C ظرفیت خازنی بار و ΔC تغییرات ظرفیت خازنی بار پس از جزیره است.

(5)

 

(6)

 

 

با توجه به روابط بالا عدم تطابق توان اکتیو بین تولید پراکنده و بار در شرایط جزیره‌ای منجر به تغییرات ولتاژ می‌شود به‌طور مشابه عدم تطابق توان راکتیو بین تولید پراکنده و بار باعث تغییرات فرکانس خواهد شد [23]؛ البته در زمان جزیره باوجود عدم تطابق توان اکتیو باعث تغییرات ولتاژ شده که این تغییرات خود بر روی توان راکتیو تأثیرگذار می‌باشد. برای توان مصرفی بار در زمان جزیره به دست می‌آید:

 

(7)

 

(8)

 

 

با استفاده از روابط (7) و (8) می‌توان نوشت:

(9)

 

 

از رابطه بالا دیده می‌شود تغییرات توان اکتیو با تغییرات فرکانس نسبت دارد [14] لذا:

(10)

 

 

از طرفی بیان شد که مشخصه افتی همان نسبت تغییرات فرکانس بر تغییرات توان اکتیو بوده و نیز در رابطه (10) تغییرات توان اکتیو تابعی از تغییرات فرکانس است:

(11)

 

 

بر اساس رابطه فوق، مشخصه افتی در شرایط جزیره‌ای تابعی از فرکانس می‌باشد.

(12)

 

(13)

 

 

این مطلب اساس روش پیشنهادی است که با استفاده از نسبت تغییرات فرکانس به تغییرات توان اکتیو اقدام به شناسایی جزیره در شبکه می‌نماید.

مشخصه افتی متداول در تولیدکننده‌های توان، مقدار ثابتی است که طبق آن، توان تولیدی در شرایط بهره‌برداری به شبکه تحویل داده می‌شود؛ اما مشخصه افتی لحظه‌ای تعریف شده مقداری متغیر است که با توجه به مقادیر لحظه‌ای فرکانس و توان اکتیو خروجی تولید پراکنده محاسبه می‌گردد. طبق مشخصه افتی متداول:

(14)

 

 

در این‌ رابطه  فرکانس لحظه‌ای خروجی تولید پراکنده،  فرکانس نامی شبکه،  توان اکتیو لحظه‌ای تولید پراکنده و  توان اکتیو تولید پراکنده در شرایط بهره‌برداری می‌باشد.

به منظور تعمیم رابطه بالا در شبکه‌های مختلف و شرایط گوناگون صورت و مخرج رابطه تبدیل به پریونیت می‌شود؛ لذا عبارت  بر صورت و  بر مخرج تقسیم می‌گردد. در نهایت مقدار مشخصه افتی لحظه‌ای ( ) تعریف شده طبق روابط زیر محاسبه می‌شود:

(15)

 

(16)

 

 

که   فرکانس لحظه‌ای خروجی تولید پراکنده برحسب Hz،   فرکانس نامی شبکه برحسب Hz،  توان اکتیو لحظه‌ای تولید پراکنده برحسب MW و  توان اکتیو تولید پراکنده در شرایط بهره‌برداری برحسب MW می‌باشد.

 

 

(الف)

 

(ب)

 

(ج)

شکل (4): مقدار مشخصه افتی فرکانسی لحظه‌ای (الف) پس از وقوع جزیره در صورت پایداری فرکانس ریزشبکه؛ (ب) پس از وقوع جزیره در صورت ناپایداری فرکانس ریزشبکه؛ (ج) پس از رخداد غیرجزیره‌ای

 

براساس شکل (1)، لختی شبکه سراسری قدرت از تغییرات کلی فرکانس جلوگیری می‌نماید؛ درحالی‌که در هنگام جزیره، بسته به مقدار عدم تطابق توان و ظرفیت‌های تولیدات پراکنده شبکه، این احتمال وجود دارد که فرکانس مقدار جدیدی به خود گرفته و یا تا رسیدن فرکانس به مقدار نامی شبکه مدت زمان زیادی طول بکشد. براساس شبیه‌سازی‌های صورت گرفته بر روی دو شبکه مورد مطالعه، شکل خروجی  در هنگام وقوع جزیره به دو صورت است. در صورت بازگشت شرایط نامی شبکه از جمله فرکانس و ولتاژ پس از وقوع جزیره فرم  به‌صورت شکل (4)-الف خواهد بود؛ اما اگر فرکانس و ولتاژ ریزشبکه پس از جزیره از مقدار نامی متفاوت باشد فرم  مانند شکل (4)-ب می‌باشد. در زمان بروز حوادث غیرجزیره‌ای مقدار  که در شکل (4)-ج آورده شده در زمان بروز خطا زیاد شده؛ اما به‌سرعت این مقدار به سمت صفر میل می‌نماید. حضور شبکه اصلی در زمان خطا باعث صفر شدن سریع مقدار  می‌شود که در غیاب شبکه اصلی مقدار  مدت زمان بیشتری غیرصفر است و این رفتار شبکه جهت تشخیص جزیره سودمند است.

چالش اصلی روش‌های غیرفعال انتخاب حد آستانه برای سنجش شاخص و انتخاب وضعیت پیشرو بوده است؛ بااین‌حال به دست آوردن یک مقدار به‌عنوان حد آستانه شاید در شرایط بهره‌برداری و شبکه متفاوت دچار خطا گردد. به این منظور ارائه یک حد آستانه انطباقی تا حد زیادی این مشکل را رفع می‌کند. در روش پیشنهادی پس از وقوع خطا در شبکه پنجره‌ای برای بررسی تغییرات  بنام پنجره مشاهده باز می‌گردد. این پنجره دارای زمان 2 سیکل است که موجب می‌شود در صورت بروز خطای غیرجزیره مقدار  از بیشینه مقدار خود عبور نماید. در پنجره مشاهده بیشینه مقدار  به‌عنوان حد آستانه روش پیشنهادی در نظر گرفته می‌شود. با توجه به مقادیر مختلف  در شرایط مختلف می‌توان گفت بیشینه مقدار  در هر سناریو عددی متفاوت بوده و حد آستانه روش پیشنهادی یک حد آستانه تطبیق‌پذیر خواهد بود.

(17)

 

 

که در اینجا  بیشینه مقدار مشخصه افتی لحظه‌ای در بازه پنجره مشاهده است. مقایسه مقدار مشخصه افتی لحظه‌ای  و حد آستانه تطبیق‌پذیر  پس از بسته شدن پنجره مشاهده و شروع پنجره تشخیص صورت می‌گیرد. حداکثر زمان پنجره تشخیص برابر 8 سیکل است. درصورتی‌که مقدار مشخصه افتی لحظه‌ای  در 7 نمونه متوالی از حد آستانه تطبیق‌پذیر  بیشتر شود پدیده جزیره رخ داده است. بدیهی است در صورت حصول شرایط جزیره در کم‌تر از 8 سیکل پنجره تشخیص بسته خواهد شد.

 

 

شکل (5): روندنمای روش پیشنهادی تشخیص جزیره

 

 

4- شبیه‌سازی و تحلیل نتایج

برای ارزیابی روش پیشنهادی، از دو شبکه استاندارد شامل شبکه فشار متوسط 11 شین Cigre و شبکه فشار متوسط 34 شین IEEE با در نظر گرفتن نیروگاه‌های بادی و خورشیدی استفاده می‌شود. شبیه‌سازی‌ها با سناریوهای مختلف در نرم‌افزار دیگسایلنت انجام می‌شود و پیاده‌سازی روش پیشنهادی در نرم‌افزار MATLAB انجام می‌گیرد. سناریوهای وقوع جزیره، تغییرات ناگهانی بار، کلیدزنی بانک خازنی، خطای اتصال کوتاه، راه‌اندازی موتور و خارج شدن منابع تولید پراکنده از مدار بر روی شبکه‌های نمونه بررسی شده است. در تمامی سناریوها به جز راه‌اندازی موتور آسنکرون، وقوع خطا در زمان 4/0 ثانیه بوده و شبیه‌سازی در طول 1 ثانیه انجام می‌شود. در سناریوهای راه‌اندازی موتور آسنکرون، شروع راه‌اندازی موتور در زمان 5/1 ثانیه و مدت زمان شبیه‌سازی برابر 5 ثانیه است.

 

4-1- شبکه فشار متوسط Cigre

شبکه فشار متوسط Cigre مورد نظر یک شبکه محلی است که توانایی تأمین توان یک روستا یا یک شهر کوچک را دارد. ولتاژ این شبکه 20 کیلوولت بوده که توسط یک پست با ولتاژ نامی 20/110 کیلوولت به شبکه اصلی متصل می‌گردد [24]. منابع تولید پراکنده این شبکه شامل نیروگاه خورشیدی دارای مبدل اینورتری، نیروگاه بادی دارای ژنراتور آسنکرون و نیروگاه دیزلی با ژنراتور سنکرون هستند. نمودار شبکه شبکه شبیه‌سازی شده در شکل (6) نشان داده شده است.

 

 

شکل (6): شبکه ولتاژ متوسط Cigre اصلاح شده [24]

 

 

4-1-1- وقوع جزیره‌ای شدن

به منظور بررسی روش پیشنهادی در تشخیص جزیره، 64 سناریو با در نظر گرفتن عدم تطابق‌ توان‌های متفاوت، از صفر تا حداکثر 200 کیلووات و 200 کیلووار و نیز تغییر میزان توان خروجی منابع تولید پراکنده، مورد بررسی قرار گرفته است. جهت نمایش عملکرد روش پیشنهادی، شرایط جزیره برای عدم تطابق توان با مقدار 100 کیلووات و 100 کیلووار، در شکل (7) مشاهده می‌شود. با توجه به این شکل، جزیره در زمان 4/0 ثانیه رخ داده است و پس از وقوع جزیره، در دیزل ژنراتور مقدار متوسط نرخ تغییرات فرکانس بر زمان به 8 هرتز بر ثانیه رسیده است و بنابراین پنجره مشاهده روش پیشنهادی شروع می‌شود. پس از 2 سیکل مقدار 155/0 پریونیت را به عنوان مقدار حد آستانه تطبیق‌پذیر برای منبع تولید پراکنده دیزل ژنراتور متصل به شین 9 تعیین می‌شود. مقادیر  در پنجره تشخیص از حد آستانه بالاتر بوده و جزیره 47 میلی‌ثانیه بعد از وقوع تشخیص داده می‌شود. همچنین نتایج روش پیشنهادی در تشخیص تعدادی از سناریوهای مربوط به جزیره در جدول (1) ارائه شده است. طبق جدول (1) در سناریو اول مقدار توان اکتیو و راکتیو مبادله شده بین ریزشبکه و شبکه اصلی برابر صفر است؛ با این حال به دلیل ساختار شبکه‌ها و منابع تولید پراکنده، فرکانس و ولتاژ در شین خروجی DG ها تغییرات جزئی خواهد داشت؛ این تغییرات موجب تغییر توان خروجی تولیدات پراکنده بعد از رخداد جزیره‌ای شدن خواهد شد. با توجه به نتایج این جدول، روش پیشنهادی جزیره را برای منابع تولید پراکنده دیزلی، بادی و خورشیدی به ترتیب در زمان‌های 50، 48 و 48 میلی‌ثانیه برای عدم تطابق توان 200 کیلووات و 200 کیلووار تشخیص می‌دهد. همچنین لازم به ذکر است که روش پیشنهادی 192 حالت جزیره‌ای شدن (64 حالت برای هر کدام از منابع تولید پراکنده) را با دقت 9/97 درصد با میانگین زمان تشخیص 50 میلی‌ثانیه تشخیص داده است.

 

 

(الف)

 

(ب)

شکل (7): (الف) مشخصه افتی فرکانسی لحظه‌ای (ب) سیگنال تشخیص برای نیروگاه دیزل ژنراتور در حالت وقوع جزیره با عدم تطابق توان 100 کیلووات و 100 کیلووار

 

جدول (1): عملکرد روش پیشنهادی در تشخیص جزیره در شبکه Cigre

ΔP

(kW)

ΔQ

(kvar)

تولید پراکنده

حد آستانه

عملکرد

زمان تشخیص (ms)

0

0

دیزلی

09/0

ü

47

بادی

39/0

ü

47

خورشیدی

21/0

ü

47

50

50

دیزلی

15/0

ü

47

بادی

57/1

ü

48

خورشیدی

84/0

ü

48

200

200

دیزلی

18/0

ü

50

بادی

55/0

ü

48

خورشیدی

39/0

ü

48

 

 

4-1-2- تغییرات ناگهانی بار

به‌منظور بررسی عملکرد روش پیشنهادی در شرایط تغییرات ناگهانی بار، تمامی بارها به میزان 25 درصد در زمان 4/0 ثانیه به ‌صورت ناگهانی افزایش می‌یابد؛ به عنوان مثال مشخصه افتی فرکانسی و عملکرد روش پیشنهادی برای نیروگاه بادی در شکل (8) نشان داده شده است. با توجه به این شکل، مقدار مشخصه افتی فرکانسی لحظه‌ای در دو سیکل ابتدایی حداکثر مقدار خود که برابر 86/0 است را دارد و پس از طی شدن دو سیکل این مقدار کاهش می‌یابد. دلیل افت سریع مشخصه افتی فرکانسی لحظه‌ای، اتصال ریزشبکه به شبکه اصلی است؛ بنابراین پس از تعیین مقدار حد آستانه در پنجره مشاهده، در پنجره تشخیص مقادیر از حد آستانه مشخص شده کم‌تر خواهد بود.

 

 

جدول (2) عملکرد روش پیشنهادی را برای تعدادی از حالت‌های تغییرات بار نشان می‌دهد. مقادیر مثبت تغییرات بار به معنی افزایش بار و مقادیر منفی به معنی کاهش بار می‌باشد. همچنین در این مطالعه، 105 حالت تغییرات بار (35 حالت برای هر کدام از منابع تولید پراکنده) به کمک روش پیشنهادی به طور کامل تشخیص داده شده‌اند.

 

 

(الف)

 

(ب)

شکل (8): (الف) مشخصه افتی فرکانسی لحظه‌ای (ب) سیگنال تشخیص برای نیروگاه بادی در حالت افزایش 25 درصدی بار

 

 

 

جدول (2): عملکرد روش پیشنهادی در شرایط تغییرات ناگهانی بار در شبکه Cigre

درصد تغییرات بار

تولید پراکنده

حد آستانه

عملکرد

20

دیزلی

95/0

ü

بادی

12/1

ü

خورشیدی

037/0

ü

40-

دیزلی

7/0

ü

بادی

2

ü

خورشیدی

046/0

ü

60

دیزلی

52/0

ü

بادی

6/1

ü

خورشیدی

035/0

ü

 

4-1-3- کلیدزنی بانک خازنی

جهت تصحیح ضریب توان و تنظیم ولتاژ استفاده از بانک خازنی در شبکه‌های قدرت متداول است. جهت بررسی رفتار روش پیشنهادی در شرایط خازن گذاری، خازن‌ها با پله‌های 50 کیلوواری در نقاط مختلف شبکه قرار داده می‌شود؛ برای مثال با توجه به شکل (6)، یک خازن 200 کیلوواری متصل به شین B03 در زمان 4/0 ثانیه وارد مدار می‌گردد. نتایج نشان می‌دهد مقدار مشخصه افتی قرکانسی لحظه‌ای در ابتدای کلیدزنی دارای مقدار بزرگی است و سپس به سرعت کم می‌شود. حد آستانه تطبیق‌پذیر در زمان شناسایی بزرگ‌تر از  بوده و این حالت به عنوان یک حادثه غیرجزیره‌ای تشخیص داده می‌شود. طبق شکل (9) مقدار حد آستانه به‌دست آمده برابر 243/0 است که مقدار مشخصه افتی فرکانسی لحظه‌ای در پنجره تشخیص به آن نزدیک نشده است؛ لذا روش پیشنهادی به‌خوبی توانسته کلیدزنی بانک خازنی را از جزیره جدا نماید. در جدول (3) خازن‌گذاری در شبکه با سه خازن متفاوت آورده شده است. در تمامی حالت‌ها، روش پیشنهادی عملکرد صحیحی نسبت به کلیدزنی بانک خازنی داشته و خروجی سیگنال تشخیص برابر صفر می‌باشد. همچنین روش پیشنهادی، 60 حالت کلیدزنی بانک خازنی (20 حالت برای هر کدام از منابع تولید پراکنده) را در این مطالعه به طور کامل از جزیره تفکیک نموده است.

 

 

(الف)

(ب)

شکل (9): (الف) مشخصه افتی فرکانسی لحظه‌ای (ب) سیگنال تشخیص وقوع جزیره برای نیروگاه خورشیدی در حالت کلیدزنی بانک خازنی 200 کیلوواری

 

جدول (3): عملکرد روش پیشنهادی در شرایط مختلف کلیدزنی بانک خازنی در شبکه Cigre

ظرفیت خازن

شین اتصال

تولید پراکنده

حد آستانه

عملکرد

50 کیلووار

3

دیزلی

038/0

ü

بادی

123/0

ü

خورشیدی

063/0

ü

100 کیلووار

3

دیزلی

076/0

ü

بادی

193/0

ü

خورشیدی

117/0

ü

150 کیلووار

3

دیزلی

095/0

ü

بادی

211/0

ü

خورشیدی

18/0

ü

 

4-1-4- خطای اتصال کوتاه

برای بررسی عملکرد روش پیشنهادی در تشخیص خطای اتصال کوتاه از جزیره، حالت‌های مختلفی با وقوع انواع اتصال کوتاه شامل خطاهای تکفاز به زمین، دو فاز به زمین، دو فاز، سه فاز و سه فاز به زمین با مقاومت خطا صفر تا 2 اهم با گام 5/0 اهم شبیه‌سازی شده‌اند؛ برای مثال در شین B07 نزدیک به نیروگاه بادی در زمان 4/0 ثانیه، اتصال کوتاه سه فاز با مقاومت 5/0 اهم رخ می‌دهد. این خطا برای مدت 100 میلی‌ثانیه اتفاق افتاده و در زمان 5/0 ثانیه رفع می‌شود. مطابق با شکل (10)، مقدار حد آستانه در این حالت برابر با 189/0 است که مقدار مشخصه افتی فرکانسی لحظه‌ای در پنجره تشخیص پس از 8 سیکل از حد آستانه عبور نمی‌کندو خطا به درستی از جزیره تفکیک می‌شود.

 

(الف)

(ب)

شکل (10): (الف) مشخصه افتی فرکانسی لحظه‌ای (ب) سیگنال تشخیص برای نیروگاه بادی در حالت اتصال کوتاه سه فاز

جدول (4): عملکرد روش پیشنهادی در شرایط اتصال کوتاه در شبکه Cigre

نوع خطا

امپدانس خطا (Ω)

محل خطا

تولید پراکنده

حد آستانه

عملکرد

دو فاز به زمین

5/0

شین 3

دیزلی

083/0

ü

بادی

13/0

ü

خورشیدی

85/0

ü

دو فاز

5/0

شین 3

دیزلی

083/0

ü

بادی

13/0

ü

خورشیدی

85/0

ü

تکفاز به زمین

5/0

شین 3

دیزلی

088/0

ü

بادی

135/0

ü

خورشیدی

9/0

ü

 

در جدول (4) عملکرد روش پیشنهادی برای خطاهای اتصال کوتاه دو فاز، دو فاز به زمین و تکفاز به زمین و محل وقوع خطا نشان داده شده است. نتایج حاکی از آن است که روش پیشنهادی این خطاها را به درستی از جزیره تشخیص می‌دهد. همچنین در این مطالعه، 72 حالت خطای اتصال کوتاه (24 حالت خطا در حضور هر کدام از منابع تولید پراکنده) به کمک روش پیشنهادی با دقت 61/98 درصد از جزیره تشخیص داده شده‌اند.

4-1-5- راه‌اندازی موتور آسنکرون

راه‌اندازی موتور آسنکرون یک چالش برای روش‌های تشخیص جزیره بوده، لذا برای ارزیابی عملکرد روش پیشنهادی، شرایط مختلفی از راه‌اندازی موتور در شبکه‌های تحت مطالعه مورد بررسی قرار گرفته است؛ برای مثال یک موتور آسنکرون با قدرت 200 کیلووات که به وسیله یک ترانسفورماتور 250 کیلوولت آمپری به شین B03 شبکه Cigre متصل است در لحظه 5/1 ثانیه راه‌اندازی می‌شود. طبق شکل (11) مقدار حد آستانه برای نیروگاه خورشیدی برابر 945/0 پریونیت است. در پنجره تشخیص مقدار  همواره کوچک‌تر از  است؛ بنابراین، روش پیشنهادی به درستی راه‌اندازی موتور را از وقوع جزیره تشخیص می‌دهد.

 

(الف)

 

(ب)

شکل (11): (الف) مشخصه افتی فرکانسی لحظه‌ای (ب) سیگنال تشخیص برای نیروگاه خورشیدی در حالت راه‌اندازی موتور 200 کیلوواتی

جدول (5) عملکرد روش را در سه حالت دیگر نشان می‌دهد؛ که روش پیشنهادی به خوبی توانسته است این حالت‌ها را از وقوع جزیره تفکیک نماید. همچنین در این مطالعه، 81 حالت راه‌اندازی موتور (27 حالت در حضور هر کدام از منابع تولید پراکنده) بررسی شد که روش پیشنهادی دارای دقت تشخیص 77/98 درصدی است.

جدول (5): عملکرد روش پیشنهادی در شرایط راه‌اندازی موتور در شبکه Cigre

توان موتور

شین اتصال

تولید پراکنده

حد آستانه

عملکرد

160 کیلووات

3

دیزلی

89/0

ü

بادی

93/0

ü

خورشیدی

79/0

ü

250 کیلووات

3

دیزلی

41/1

ü

بادی

19/1

ü

خورشیدی

29/1

ü

400 کیلووات

3

دیزلی

2/1

ü

بادی

7/1

ü

خورشیدی

99/1

ü

 

4-1-6- خارج‌شدن منبع تولید پراکنده از مدار

در شرایط بهره‌برداری، در مواقعی منابع تولید پراکنده به دلایل مختلفی از جمله وقوع خطا از مدار خارج می‌گردند. برای بررسی عملکرد روش پیشنهادی در حالت خارج شدن منبع تولید پراکنده از مدار، نیروگاه خورشیدی در زمان 4/0 ثانیه از مدار خارج می‌گردد. مطابق شکل (12) مقدار حد آستانه‌ برای نیروگاه دیزل ژنراتور برابر 478/0 پریونیت به دست می‌آید که در بازه زمانی پنجره تشخیص، این مقدار همواره بزرگ‌تر از مشخصه افتی فرکانسی لحظه‌ای است؛ بنابراین روش پیشنهادی در حالت خارج شدن منبع تولید پراکنده از مدار مانند دیگر حالت‌های غیرجزیره‌ای، عملکرد مناسبی داشته است. در جدول (6) دو حالت دیگر خروج منبع تولید پراکنده از مدار نشان داده شده است که مطابق نتایج این جدول، روش پیشنهادی دارای عملکرد مطلوب می‌باشد. همچنین در این مطالعه، 24 حالت خروج منبع تولید پراکنده از مدار (8 حالت در حضور هر کدام از منابع تولید پراکنده) به کمک روش پیشنهادی به طور کامل تشخیص داده شده‌اند.

 

 

(الف)

(ب)

شکل (12): (الف) مشخصه افتی فرکانسی لحظه‌ای (ب) سیگنال تشخیص برای دیزل ژنراتور در حالت خارج‌شدن نیروگاه خورشیدی از مدار

 

جدول (6): عملکرد روش پیشنهادی در شرایط خروج منابع تولید پراکنده از مدار در شبکه Cigre

تولید پراکنده خارج از مدار

تولید پراکنده

حد آستانه

عملکرد

دیزلی

بادی

98/0

ü

خورشیدی

225/0

ü

بادی

دیزلی

744/0

ü

خورشیدی

216/0

ü

 

4-2- شبکه 34 شین IEEE

برای نشان دادن کارایی روش پیشنهادی در شبکه بزرگتر و عدم وابستگی به شبکه مورد مطالعه، شبکه 34 شین IEEE نیز در این مطالعه مورد بررسی قرار می‌گیرد. نمودار شبکه در شکل (13) نشان داده است. این شبکه، یک شبکه توزیع نامتعادل است که به صورت شعاعی توسعه یافته است. فرکانس شبکه برابر 60 هرتز و ولتاژ آن برابر 9/24 کیلوولت است که ولتاژ در شین‌های 888 و 890 برابر 16/4 کیلوولت می‌باشد [25]. منابع تولید پراکنده این شبکه شامل نیروگاه خورشیدی با مبدل اینورتری، نیروگاه بادی با ژنراتور آسنکرون و نیروگاه دیزلی با ژنراتور سنکرون می‌باشند.

 

 

 

شکل (13): شبکه ولتاژ متوسط 34 شین IEEE [25]

 

 

 

تمامی حالت‌های وقوع جزیره، تغییرات ناگهانی بار، کلیدزنی بانک خازنی، انواع خطاهای اتصال کوتاه، راه‌اندازی موتور و خارج شدن تولید پراکنده از مدار که در بخش 4-1 بر روی شبکه Cigre شبیه‎سازی شده‌اند بر روی شبکه 34 شین IEEE نیز شبیه‌سازی می‌شوند. نتایج بعضی از این حالت‌ها در جدول‌های (7) تا (12) نشان داده شده است. با توجه به نتایج ارائه شده، روش پیشنهادی علاوه بر تشخیص مناسب وقوع جزیره‌ای شدن در شبکه 34 شین IEEE، توانسته عملکرد خوبی در تفکیک حوادث دیگر از جزیره‌ داشته باشد.

جدول (7): عملکرد روش پیشنهادی در شرایط وقوع جزیره‌ای شدن در شبکه 34 شین IEEE

ΔP

(kW)

ΔQ

(kvar)

تولید پراکنده

حد آستانه

عملکرد

زمان تشخیص (ms)

0

0

دیزلی

114/0

ü

42

بادی

114/0

ü

43

خورشیدی

085/0

ü

41

50

50

دیزلی

12/0

ü

41

بادی

48/0

ü

41

خورشیدی

41/0

ü

41

200

200

دیزلی

80/0

ü

160

بادی

37/0

ü

155

خورشیدی

85/0

ü

98

 

 

 

جدول (8): عملکرد روش پیشنهادی در شرایط مختلف تغییرات ناگهانی بار در شبکه 34 شین IEEE

درصد تغییرات بار

تولید پراکنده

حد آستانه

عملکرد

20

دیزلی

07/0

ü

بادی

15/0

ü

خورشیدی

42/0

ü

40-

دیزلی

15/0

ü

بادی

33/0

ü

خورشیدی

75/0

ü

60

دیزلی

18/0

ü

بادی

17/0

ü

خورشیدی

9/0

ü

جدول (9): عملکرد روش پیشنهادی در شرایط مختلف کلیدزنی بانک خازنی در شبکه 34 شین IEEE

ظرفیت خازن

شین اتصال

تولید پراکنده

حد آستانه

عملکرد

50 کیلووار

824

دیزلی

064/0

ü

بادی

04/0

ü

خورشیدی

015/0

ü

100 کیلووار

824

دیزلی

129/0

ü

بادی

08/0

ü

خورشیدی

03/0

ü

150 کیلووار

890

دیزلی

342/0

ü

بادی

424/0

ü

خورشیدی

255/0

ü

 

 

 

جدول (10): عملکرد روش پیشنهادی در شرایط اتصال کوتاه در شبکه 34 شین IEEE

نوع خطا

امپدانس خطا (Ω)

محل خطا

تولید پراکنده

حد آستانه

عملکرد

دو فاز به زمین

5/0

شین 816

دیزلی

045/0

ü

بادی

09/0

ü

خورشیدی

053/0

ü

دو فاز

5/0

شین 830

دیزلی

13/0

ü

بادی

109/0

ü

خورشیدی

056/0

ü

تکفاز به زمین

5/0

شین 826

دیزلی

16/0

ü

بادی

067/0

ü

خورشیدی

07/0

ü

جدول (11): عملکرد روش پیشنهادی در شرایط راه‌اندازی موتور در شبکه 34 شین IEEE

توان موتور

شین اتصال

تولید پراکنده

حد آستانه

عملکرد

160 کیلووات

858

دیزلی

33/0

ü

بادی

58/2

ü

خورشیدی

63/2

ü

250 کیلووات

816

دیزلی

305/0

ü

بادی

34/1

ü

خورشیدی

58/1

ü

400 کیلووات

840

دیزلی

18/0

ü

بادی

2/0

ü

خورشیدی

09/2

ü

 

جدول (12): عملکرد روش در شرایط خروج منبع تولید پراکنده از مدار در شبکه 34 شین IEEE

تولید پراکنده خارج از مدار

تولید پراکنده

حد آستانه

عملکرد

دیزلی

بادی

165/0

ü

خورشیدی

435/0

ü

بادی

دیزلی

66/0

ü

خورشیدی

21/0

ü

 

4-3- تحلیل نتایج

در شبیه‌سازی‌ها از حضور سه DG خورشیدی، بادی و دیزل ژنراتور استفاده شده است. این تولیدات پراکنده از فنّاوری‌های مختلف اینورتری، سنکرون و آسنکرون هستند و در هر سه مورد عملکرد روش پیشنهادی مناسب بوده است. طبق نتایج شبیه‌سازی‌ها، جمع‌بندی سناریوهای وقوع جزیره، تغییرات ناگهانی بار، کلیدزنی خازنی، خطای اتصال کوتاه، راه‌اندازی موتور و خارج شدن منبع تولید پراکنده از مدار، مطابق جدول (13) و Error! Reference source not found.) می‌باشد. روش پیشنهادی طبق سناریوهای صورت گرفته در حادثه جزیره‌ای شدن، برای دو شبکه مورد بررسی به طور میانگین دارای 4/97 درصد صحت تشخیص می‌باشد، همچنین در خصوص دیگر حوادث شبکه می‌تواند عملکرد قابل قبولی داشته باشد.

در ادامه جهت نشان دادن کارایی روش پیشنهادی مقایسه‌ای بین روش پیشنهادی با برخی روش‌‌های غیرفعال ارائه شده، صورت می‌گیرد. در [16] روش تشخیص جزیره با وجود محاسبات بیش‌تر و پیچیده‌تر برای شبکه مورد مطالعه Cigre در بهینه‌ترین حالت دارای حداکثر دقت تشخیص 8/97 درصدی بوده در حالی که روش پیشنهادی با الگوریتم ساده‌تر دارای تشخیص 9/98 درصدی برای این شبکه می‌باشد. مرجع [16] شرایط جزیره‌ای شدن، تغییرات بار، خطای اتصال کوتاه، کلیدزنی بانک خازنی و خارج شدن DG از سرویس مورد بررسی قرار گرفته که این مطالعه، علاوه بر آن‌ها شرایط راه‌اندازی موتور آسنکرون نیز بررسی شده است. همچنین تولیدات پراکنده مبتنی بر اینورتر، ژنراتور سنکرون و ژنراتور آسنکرون در این مطالعه به کار گرفته شده در صورتی که [16] تولید پراکنده اینورتری را بررسی نکرده است.

مطابق جدول (15) خطای تشخیص، زمان تشخیص و ناحیه غیرقابل تشخیص روش پیشنهادی با روش‌های متداول غیرفعال نشان داده شده است. در روش پیشنهادی زمان تشخیص حدود 50 میلی‌ثانیه بوده که این زمان نسبت به روش نرخ تغییر فرکانس بر تغییر توان اکتیو زمان کوتاه‌تری است. در خصوص خطای تشخیص این روش نسبت به دیگر روش‌های متداول غیرفعال در زمان اتفاقات شبکه از جمله اتصال کوتاه، راه‌اندازی موتور آسنکرون و کلیدزنی بانک خازنی خطای کم‌تری دارد. ناحیه غیرقابل تشخیص روش پیشنهادی نسبت به دیگر روش‌های گفته شده در جدول (15) کوچک‌تر می‌باشد که در Error! Reference source not found.) NDZ برای ولتاژ یک پریونیت، فرکانس 60 هرتز و ضریب کیفیت 5/2 نشان داده شده است. همچنین روش پیشنهادی هیچ‌گونه تأثیر منفی بر روی کیفیت توان شبکه نیز ندارد.

جدول (13): عملکرد روش پیشنهادی در حالت‌های بررسی‌شده در شبکه Cigre

نوع حادثه

تعداد حالت‌ها

درصد تشخیص

جزیره‌ای شدن

192

92/97

تغییرات ناگهانی بار

105

100

کلیدزنی بانک خازنی

60

100

اتصال کوتاه

72

61/98

راه‌اندازی موتور آسنکرون

81

77/98

خارج‌شدن تولید پراکنده

24

100

مجموع

534

88/98

 

جدول (14): عملکرد روش پیشنهادی در حالت‌های بررسی‌شده در شبکه 34 شین IEEE

نوع حادثه

تعداد حالت‌ها

درصد تشخیص

جزیره‌ای شدن

192

88/96

تغییرات ناگهانی بار

105

10/98

کلیدزنی بانک خازنی

60

100

اتصال کوتاه

72

83/95

راه‌اندازی موتور آسنکرون

81

77/98

خارج‌شدن تولید پراکنده

24

100

مجموع

534

75/97

 

جدول (15): دقت روش‌های متداول غیرفعال و روش پیشنهادی

عنوان روش

خطای تشخیص

زمان تشخیص

NDZ

اضافه/کاهش فرکانس [22]

زیاد

بین 4 میلی‌ثانیه تا 1 ثانیه

بزرگ

اضافه/کاهش ولتاژ [22]

زیاد

بین 4 میلی‌ثانیه تا 2 ثانیه

بزرگ

نرخ تغییر فرکانس [5]

کم

کمتر از 1 ثانیه

کوچک

نرخ تغییر فرکانس بر تغییر توان اکتیو [5]

کم

حدود 100 میلی‌ثانیه

کوچک

روش پیشنهادی

بسیار کم

کمتر از 50 میلی‌ثانیه

ناچیز

 

 

 

 

شکل (14): ناحیه غیرقابل تشخیص روش پیشنهادی

5- نتیجه‌گیری

در این مقاله، یک روش غیرفعال برای تشخیص جزیره از دیگر اغتشاشات مانند تغییرات بار، کلیدزنی خازنی، انواع خطاهای اتصال کوتاه، راه‌اندازی موتور و خارج شدن DG با استفاده از مشخصه افتی فرکانسی لحظه‌ای که دارای آستانه تطبیق‌پذیر دینامیکی می‌باشد، پیشنهاد شده است. روش پیشنهادی باعث حل مشکل انتخاب حد آستانه در روش‌های غیرفعال برای تشخیص جزیره می‌شود. این روش در دو شبکه Cigre و 34 شین IEEE در سناریوهای مختلف با در نظر گرفتن منابع تولید پراکنده خورشیدی، بادی و دیزل ژنراتور مورد ارزیابی قرار گرفته است. با توجه به نتایج شبیه‌سازی روش پیشنهادی دارای ناحیه NDZ کوچک است و همچنین دارای سرعت و دقت تشخیص مناسب در حضور انواع منابع تولید پراکنده می‌باشد.

 

[1] تاریخ ارسال مقاله: 29/01/1402

تاریخ پذیرش مقاله: 22/03/1402

نام نویسندۀ مسئول: مهدی بانژاد

نشانی نویسندۀ مسئول: ایران، شاهرود، دانشگاه صنعتی شاهرود، دانشکده مهندسی برق

 

[1] Renewable energy sources (RESs)

[2] Distributed generation (DG)

[3] Wind Power Plant

[4] Microturbine

[5] Solar power plant

[6] Combined heat and power (CHP)

[7] Non detection zone (NDZ)

[8] Islanding detection methods (IDMs)

[9] Local methods

[10] Remote methods

[11] Active method

[12] Passive method

[13] Point of common coupling (PCC)

[14] Droop control

[15] Master-slave control

[16] Average current-sharing control

 [1]  Raza S, Mokhlis H, Arof H, Laghari JA, Wang L, "Application of Signal Processing Techniques for Islanding Detection of Distributed Generation in Distribution Network: A Review," Energy Conversion and Management, Vol. 96, pp. 613-624, 2015.
[2]  Pouryekta A, Ramachandaramurthy VK, "A Hybrid Islanding Detection Method for Distribution Systems," Distributed Generation & Alternative Energy Journal, Vol. 33, No. 4, pp. 44-67, 2018.
[3]  Abd-Elkader AG, Saleh SM, Eiteba MM, "A Passive Islanding Detection Strategy for Multi-distributed Generations," International Journal of Electrical Power & Energy Systems, Vol. 99, pp. 146-155, 2018/07/01/ 2018.
[4]  Khamis A, Shareef H, Bizkevelci E, Khatib T, "A Review of Islanding Detection Techniques for Renewable Distributed Generation Systems," Renewable and Sustainable Energy Reviews, Vol. 28, pp. 483-493, 2013.
[5]  Vieira JC, Freitas W, Xu W, Morelato A, "Efficient Coordination of ROCOF and Frequency Relays for Distributed Generation Protection by Using the Application Region," IEEE Transactions on Power Delivery, Vol. 21, No. 4, pp. 1878-1884, 2006.
 [6] Khan MA, Haque A, Kurukuru VB, Saad M, "Islanding Detection Techniques for Grid-Connected Photovoltaic Systems-A Review," Renewable and Sustainable Energy Reviews, Vol. 154, p. 111854, 2022.
[7]  Islam MM, Nagrial M, Rizk J, Hellany A, "General Aspects, Islanding Detection, and Energy Management in Microgrids: A Review," Sustainability, Vol. 13, No. 16, p. 9301, 2021.
[8]  Bayrak G, Kabalci E, "Implementation of a New Remote Islanding Detection Method for Wind–solar Hybrid Power Plants," Renewable and Sustainable Energy Reviews, Vol. 58, pp. 1-15, 2016/05/01/ 2016.
[9]  Emadi A, Afrakhte H, Sadeh J, "Fast Active Islanding Detection Method Based on Second Harmonic Drifting for Inverter-based Distributed Generation," IET Generation, Transmission & Distribution, Vol. 10, No. 14, pp. 3470-3480, 2016.
[10]   Gupta P, Bhatia RS, Jain DK, "Active ROCOF Relay for Islanding Detection," IEEE Transactions on Power Delivery, vol. 32, No. 1, pp. 420-429, 2017.
[11]   Bakhshi R, Sadeh J, "Voltage Positive Feedback Based Active Method for Islanding Detection of Photovoltaic System with String Inverter Using Sliding Mode Controller," Solar Energy, Vol. 137, pp. 564-577, 2016/11/01/ 2016.
[12]   Chen X, Li Y. Li, "An Islanding Detection Method for Inverter-Based Distributed Generators Based on the Reactive Power Disturbance," IEEE Transactions on Power Electronics, Vol. 31, No. 5, pp. 3559-3574, 2016.
[13]   Pai FS, Huang SJ, "A Detection Algorithm for Islanding-prevention of Dispersed Consumer-owned Storage and Generating Units," IEEE Transactions on Energy Conversion, Vol. 16, No. 4, pp. 346-351, 2001.
[14]   Raza S, Mokhlis H, Arof H, Laghari JA, Mohamad H, "A Sensitivity Analysis of Different Power System Parameters on Islanding Detection," IEEE Transactions on Sustainable Energy, Vol. 7, No. 2, pp. 461-470, 2016.
[15]   Kamyab, E., Sadeh, J., Yaghobi, H, "Artificial Neural Network Approach for Islanding Detection in Inverter Based Distributed Generator with a Forced Transient in System Frequency, " Computational Intelligence in Electrical Engineering, 2015; 5(4): 103-114.
[16]   Heidari orejloo, M., Seifossadat, S., Gharibreza, E. "Design and Simulation of Intelligent based Relay for Power Islanding Detection," Computational Intelligence in Electrical Engineering, 2011; 2(3): 67-78.
[17]   Tayab UB, Roslan MA, Hwai LJ, Kashif M, "A Review of Droop Control Techniques for Microgrid," Renewable and Sustainable Energy Reviews, Vol. 76, pp. 717-727, 2017.
[18]   J Chapman S, Electric Machinery Fundamentals. McGraw-Hill, 2004.
[19]   Moradi MH, Foroutan VB, Abedini M, "Power Flow Analysis in Islanded Micro-Grids via Modeling Different Operational Modes of DGs: A Review and a New Approach," Renewable and Sustainable Energy Reviews, Vol. 69, pp. 248-262, 2017.
[20]   Mohamed YA, El-Saadany EF, "Adaptive Decentralized Droop Controller to Preserve Power Sharing Stability of Paralleled Inverters in Distributed Generation Microgrids," IEEE Transactions on Power Electronics, Vol. 23, No. 6, pp. 2806-2816, 2008.
[21]   Guerrero JM, Vasquez JC, Matas J, De Vicuña LG, Castilla M, "Hierarchical Control of Droop-controlled AC and DC Microgrids- A General Approach Toward Standardization," IEEE Transactions on Industrial Electronics, Vol. 58, No. 1, pp. 158-172, 2010.
 [22]  Freitas W, Xu W, Affonso CM, Huang Z, "Comparative Analysis Between ROCOF and Vector Surge Relays for Distributed Generation Applications," IEEE Transactions on Power Delivery, Vol. 20, No. 2, pp. 1315-1324, 2005.
[23]   Zeineldin HH, El-Saadany EF, Salama MM, "Impact of DG Interface Control on Islanding Detection and Nondetection Zones," IEEE Transactions on Power Delivery, Vol. 21, No. 3, pp. 1515-1523, 2006.
[24]   Rudion K, Orths A, Styczynski ZA, Strunz K., "Design of Benchmark of Medium Voltage Distribution Network for Investigation of DG Integration," Presented at the 2006 IEEE Power Engineering Society General Meeting, 2006.
[25]   Dugan RC, Kersting WH, "Induction Machine Test Case for the 34-bus Test Feeder-Description," Presented at the 2006 IEEE Power Engineering Society General Meeting, 2006.