Security-Constraint scheduling of generation units and Demand Response

Document Type : Research Article

Authors

1 Shahid Beheshti University

2 Imam Hussein University

Abstract

In this paper, a security-constraint demand response and unit commitment problem (SCDRUC) is proposed in which demand response programs (DRPs) are modeled as commercial resources beside the generation units. Indeed, the main aim of this paper is to introduce the security-constraint unit commitment problem in which the presence of DRPs is modeled as the available and necessary resources. Therefore, by integrating the generation units with demand side resources, the SCDRUC problem is proposed. The optimization results of the proposed problem will determine the generation of units as well as the demand response participation plan in order to minimize the operational cost while security constraints are satisfied. The IEEE 24-bus test system is utilized to investigate the simulation results. All the numerical analysis approves the performance and efficiency of the proposed model.

Keywords


1- مقدمه

[1]

برنامه‌های پاسخ‌گویی بار اقدامات کوتاه‌مدتی هستند که توسط مصرف‌کننده‌ها برای تعدیل مصرف انرژی انجام می‏شوند [1]. در برنامه راهبردی آژانس بین‌المللی انرژی، برنامه‌های پاسخ‌گویی بار به علت منافعی که در سطوح بهره‌برداری و اقتصادی دارند، به عنوان نخستین انتخاب در تمامی سیاست‌های انرژی معرفی شده‌اند [2]. همان‌طور که در شکل (1) نشان داده شده است، برنامه‌های پاسخ‌گویی بار را می‌توان در قالب سه دسته کلی گروه‌بندی کرد: برنامه‌های زمان‌محور[1] (TBPs)، برنامه‌های تشویق‌محور[2] (IBPs) و برنامه‌های مبتنی بر بازار[3] (MBPs) [3]. در برنامه‌های زمان‌محور، قیمت برق در دوره‌های مختلف تغییر می‌کند‌. بنابراین، مصرف‌کننده‌ها باید میزان مصرف خود را مطابق با تعرفه‌های موجود تعدیل سازند. در برنامه‌های تشویق‌محور، مصرف‌کننده‌ها توسط بهره‌بردار مستقل سیستم و یا شرکت‌های دیگر، تشویق می‌شوند که میزان مصرف خود را کاهش دهند. در برنامه‌های مبتنی بر بازار، تمامی بازیگران بازار به دو گروه تقسیم می‌شوند: خریدارانDRو فراهم‌کنندگان DR (DRSP[4]). خریداران DR برای بهبود قابلیت اطمینان سیستم و تجارت وابسته به برق خود به DR احتیاج دارند. فراهم‌کنندگانDR ظرفیت لازم برای شرکت در این برنامه‌ها را دارند و با مشارکت در برنامه‌های DR به منافعی از نظر اقتصادی می‌رسند.

حضور برنامه‌های پاسخ‌گویی بار، می‌تواند در تولید توان و همچنین عرضه آن تأثیرگذار باشد. تأمین بار و رزرو، همواره یکی از مهم‌ترین دغدغه‌ها در سیستم قدرت بوده است. برنامه‌های پاسخ‌گویی بار می‌توانند برای تأمین مقدار بار و ذخیره مورد نیاز شبکه استفاده شوند. این برنامه‌ها ابزاری مفید برای بهره‌بردار مستقل سیستم[5] (ISO) هستند و می‌توانند در شرایط بحرانی شبکه فراخوانی شوند. مسأله مشارکت واحدها[6] (UC) و برنامه‌ریزی تولید[7] همواره به‏عنوان مسأله‌ای اساسی در سیستم قدرت مطرح بوده است. UC حالات روشن و خاموش بودن واحدها و همچنین، میزان تولید آن‌ها را به منظور تأمین بار پیش‌بینی شده در یک افق زمانی مشخص، تعیین می‌کند‌‌ [4]. برنامه بهینه تولید باید هزینه بهره‌برداری از شبکه در مدت‌زمان بهره‌برداری از آن را در حین تأمین بار، نیازهای مرتبط با ذخیره چرخان و تمامی قیود شبکه کمینه کند [5- 8]. در دهه‌های اخیر، روش‌های قطعی[8]، ابتکاری[9]، و سایر روش‌های ترکیبی[10] برای حل مسأله مشارکت واحدها ارایه شده است. در حالت کلی، در صورتی که اندازه مسأله بهینه‌سازی متوسط و یا بزرگ باشد، روش‌های قطعی قادر به یافتن جواب مسأله در مدت‌زمان لازم نخواهند بود [9 و 10]. این محدودیت‌ها به معرفی روش‌های ابتکاری منجر شد [11- 16]. الگوریتم‌های بهینه‌سازی ابتکاری برای حل مسائل بهینه‌سازی پیچیده، مزایایی نسبت به الگوریتم‌های حل قطعی دارند، اما مشکل اساسی این روش‌ها این است که لزوماً به بهترین جواب ممکن نمی‌رسند و ممکن است نقطه بهینه محلی را به جای جواب بهینه کلی ارایه دهند. برای بهبود روش‌های حل مسأله مشارکت واحدها، الگوریتم‌های ترکیبی معرفی و استفاده شده‏اند [17- 20].

اگر در حل مسأله برنامه‌ریزی مشارکت واحدها، حفظ امنیت مورد تأکید باشد، مسأله جدید را برنامه‌ریزی مقید به امنیت مشارکت واحدها[11] (SCUC) می‌نامند. در این نوع برنامه‌ریزی، تأمین بار قیدی سخت و اجباری است. حداکثرسازی امنیت نیز بیشتر توسط رزروها و خدمات جانبی مطرح می‌شود. حداقل‌سازی هزینه نیز با مشارکت دادن واحدهای ارزان‌تر و در عین حال با ارضای قیود مربوطه انجام می‏شود.

هدف از برنامه ریزی امنیت– مقید مشارکت واحدها، برنامه‌ریزی مشارکت واحدها باکم‌ترین هزینه تولید بدون تأثیر منفی بر قابلیت‌اطمینان سیستم است. قابلیت‌اطمینان یک سیستم به شکل برآورده کردن دو تابع کفایت و امنیت تفسیر می‌شود. مقدار کافی از ظرفیت منابع برای تأمین تقاضای بیشینه باید موجود باشد (کفایت) و سیستم باید بتواند تغییرات و پیشامدها را بر مبنای روزانه و ساعتی تحمل کند (امنیت).

SCUC را می‌توان به شکل یک مسأله بهینه‌سازی بیان کرد که در آن هزینه‌های بهره‌برداری بر اساس هزینه‌های افزایشی حاصل از واحدهای تولید، کمینه می‌شود. تابع هدف این مسأله همانند مسأله UC است و تنها تفاوت این مسأله با مسأله UC در تعداد قیودی است که بررسی می‏شوند.

تأمین بار و ذخیره مورد نیاز شبکه، به عنوان یک مسأله اساسی در برنامه‌ریزی تولید مطرح است. همان‌طور که اشاره شد، برنامه‌های پاسخ‌گویی بار ابزاری مفید برای ISO هستند و می‌توانند در مدت‌زمان کمابیش کوتاه فراخوانی شوند. بنابراین، استفاده از برنامه‌های پاسخ‌گویی بار می‌تواند برای رفع برخی از نگرانی‌های مربوط به تأمین بار و ذخیره در شبکه استفاده شوند. در سال‌های اخیر، مدل اقتصادی TBP و IBP ها در پژوهش‏های زیادی مورد توجه قرار گرفته‌اند [21- 28]. والرو و همکارانش [21] سیاست‌های مربوط به DR و مشتریان را در بازارهای برق بررسی کرده‌اند. در مراجع [22 و 23] مدل اقتصادی بارهای پاسخ‌گو به قیمت[12] بر اساس مقادیر ثابت الاستیسیته[13] بار ارایه شده است. برنامه‌های تسویه بازار[14] با در نظر گرفتن منافع اقتصادی آن‌ها در مرجع [24] بحث شده‌اند. در مطالعاتی که در مراجع [25- 28] انجام شده است، مدل اقتصادی بارهای پاسخ‌گو به قیمت بررسی شده‌اند. در مرجع [29]، ساختار خطی با ضرایب ثابت برای منحنی فروش DRSP ها برای شرکت در برنامه‌های مبتنی بر بازار ارایه شده است. با توجه به پتانسیل برنامه‌های DR، نیاز به طرح مسأله‌ای که طی آن، مقدار بهینه توان تولیدی واحدها و DR فراخوانی شده در کنار هم محاسبه شوند، ضروری به نظر می‌رسد؛ زیرا حل این مسأله، می‌تواند علاوه بر به تعویق انداختن نیاز به ساخت زیربنای تولیدی جدید، موجب کاهش هزینه‌های بهره‌برداری از شبکه شود.

 

 
   

 

 

 

 

 

 

 

 

 

 

 

 


شکل (1): تقسیم‌بندی برنامه‌های پاسخ‌گویی بار

 

 

 

 

 

 

 

 

 

 

 

 

 

 


شکل(2): چارچوب کلی تجارت برنامه‌های پاسخ‌گویی بار از دیدگاه بهره‌بردار مستقل سیستم

 

 

همان‌طور که بیان شد، برنامه‌های پاسخ‌گویی بار قابلیت فراخوانی و اجرا در مدت‌زمان کمابیش کوتاه را دارند [30] و از این‌رو می‌توانند به عنوان ابزاری مفید در اختیار بهره‌بردار سیستم قرار گیرند. در مرجع [25]، مدل‌های اقتصادی بارهای پاسخ‌گو براساس کشش قیمتی تقاضا[15] و تابع سود مشتریان[16] توسعه یافته‌اند. برنامه‌های قطع بار و پاسخ‌گویی بار اضطراری در این مرجع مطالعه و اثر آن‌ها در برنامه‌ریزی مشارکت واحدها بررسی شده است. در مرجع [31]، مسأله مشارکت واحدها با در نظر داشتن عدم‌قطعیت در تولید توان واحدهای بادی و نیز حضور بارهای پاسخ‌گو مدل‌سازی شده است. در این مرجع، هدف بهره‌بردار مستقل سیستم بیشینه‌سازی رفاه اجتماهی تحت بدترین حالات ممکن تولید توان بادی و مشارکت DR است. مسأله پیشنهادی در [31]، به کمک روش تجزیه بندرز حل شده است. در مرجع [32]، برنامه‌ریزی مشارکت واحدها با در نظر داشتن برنامه‌های زمان- محور پاسخ‌گویی بار و نیز مزارع بادی مدل‌سازی شده است. در مرجع [33]، شاخص‌هایی برای ارزیابی پاسخ مشتریان معرفی شده‌اند. سپس، برنامه‌ریزی تولید و پاسخ‌گویی بار در این مرجع بررسی شده است. در مرجع [34]، نحوه مشارکت فراهم‌کنندگان DR و مکانیزم قیمت‌دهی آن‌ها در شبکه مدل‌سازی شده است. در مرجع [3]، روشی پویا برای مشارکت مشتریان در پاسخ‌گویی بار پیشنهاد شده است. سپس، برنامه‌ریزی مشارکت واحدها در حضور بارهای پاسخ‌گو بحث شده است. مرجع [35]، به برنامه‌ریزی مشارکت واحدها در حضور برنامه‌های پاسخ‌گویی بار اضطراری و مقایسه آن با قراردادهای قطع بار پرداخته است. در مرجع [36]، برنامه‌ریزی تصادفی مشارکت واحدها در حضور منابع پاسخ‌گویی بار معرفی شده است. در مرجع [37] تأثیر برنامه‌های پاسخ‌گویی بار قیمت زمان پیک بحرانی در مسأله برنامه‌ریزی مشارکت واحدها بررسی شده است. در مراجع [38 و 39] نیز برنامه‌های پاسخ‌گویی بار قیمت زمان واقعی در کنار برنامه‌ریزی مشارکت واحدها لحاظ شده‌اند. مرجع [40] اثر پاسخ‌گویی بار در برنامه‌ریزی تولید در بازار روز- پیش[17] را بررسی کرده است. در مرجع [41]، اثر برنامه‌های پاسخ‌گویی بار قیمت زمان واقعی در تغییر الگوی بار مشترکین و برنامه‌ریزی مشارکت واحدها بررسی شده است. در این مرجع، استفاده از پتانسیل برنامه‌های پاسخ‌گویی بار در پوشش خطای پیش‌بینی تولید نیروگاه بادی نیز بررسی شده است. در مرجع [42]، مدل اقتصادی بارهای پاسخ‌گو با استفاده از مفهوم کشش قیمتی تقاضا استخراج و اثر منابع سمت تقاضا در برنامه‌ریزی تولید بررسی شده است. مرور و بررسی کارهای انجام شده، نشان می‏دهد که تاکنون مطالعه جامعی به شکل مدل‌سازی برنامه‌های مبتنی بر بازار پاسخ‌گویی بار[18] و تجمیع مناسب آن با مسأله برنامه‌ریزی امنیت- مقید مشارکت واحدها انجام نشده است. بنابراین، نوآوری اصلی این مقاله، مدل‌سازی برنامه‌های مبتنی بر بازار پاسخ‌گویی بار به شکل منابعی در دسترس و ترکیب آن در دل مسأله برنامه‌ریزی امنیت- مقید مشارکت واحدهاست. در این راستا، با مدل‌سازی برنامه‌های مبتنی بر بازار پاسخ‌گویی بار، منحنی عرضه پاسخ‌گویی بار حاصل شده و توابع هزینه آن‌ها به طور هم‏زمان در کنار توابع هزینه تولید واحدهای نیروگاهی قرار می‌گیرد. بنابراین، یک تابع هدف واحد برای بهره‌بردار مستقل سیستم حاصل می‌شود که با حل آن، مقدار توان تولیدی واحدها به طور هم‌زمان با مقدار فعال‌سازی DR مشخص می‌شود. بنابراین، در این مقاله، مسأله برنامه‌ریزی امنیت- مقید واحدهای تولیدی و بارهای پاسخ‌گو[19] (SCDRUC) طرح می‌شود که ضمن آن، منحنی‌های عرضه DR مدل‌سازی شده و به عنوان کالایی تجاری در کنار منحنی‌های تولید توان ژنراتورها برای حل مسأله مشارکت واحدها قرار می‌گیرند. نتایج حاصل، برنامه‌ریزی بهینه تولید واحدها را به طور همزمان با میزان بهینه اجرای برنامه‌های پاسخ‌گویی بار، در حین برآورد امنیت سیستم تعیین خواهد کرد.

در ادامه این مقاله، ابتدا در بخش 2، چارچوب کلی مسأله معرفی شده بیان می‌شود. در بخش 3، مسأله پیشنهادی مقاله بیان و فرمول‌بندی می‌شود. این بخش، شامل بیان روابط مربوط به ترم‌های مختلف تابع هدف و قیود مسأله پیشنهادی است. سپس، نتایج حاصل از شبیه‌سازی و انجام مطالعات عددی در بخش 4 ارایه می‌شود. در پایان، نتایج حاصل در بخش 5 بیان می‌شود.

 

2- چارچوب کلی مسأله

در شکل (2)، تصویر کلی چارچوب ارایه شده نشان داده شده است. شکل (2)، نشان‌دهنده ارتباط ایجاد شده بین منابع پاسخ‌گویی بار[20] (DRRs) و منابع سمت تولید است، به طوری که بهره‌بردار مستقل سیستم، این منابع را از سمت تولید و تقاضا جمع‌آوری کرده و با هدف کمینه‌سازی هزینه بهره‌برداری از شبکه در حین حفظ امنیت سیستم، مقدار توان تولیدی ژنراتورها و نیز مقدار DR فعال‌سازی شده در هر لحظه را تعیین می‌کند‌.

همان‌طور که در شکل (2) مشاهده می‌شود، طی چارچوب پیشنهادی، ابتدا رفتار فراهم‌کنندگان DR، به عنوان منابعی مجازی و در دسترس مدل‌سازی شده است. در این راستا، برنامه‌های مبتنی بر بازار پاسخ‌گویی بار مد نظر قرار گرفته‌اند که در آن، فراهم‌کنندگان DR، منحنی عرضه خود را به بهره‌بردار مستقل سیستم ارایه می‌دهند. از طرف دیگر، بهره‌بردار سیستم، منحنی عرضه توان واحدهای تولیدی را نیز دریافت می‌کند. طی مسأله پیشنهادی مقاله، بهره‌بردار سیستم قدرت منحنی عرضه پاسخ‌گویی بار را همانند واحدهایی تولیدی در نظر می‌گیرد که می‌توانند در برآورد بار و قیود شبکه سهیم باشند. به این طریق، ظرفیت پاسخ‌گویی بار را می‌توان به عنوان منابعی در نظر گرفت که در کنار سایر واحدهای تولیدی در اختیار بهره‌بردار سیستم قرار می‌گیرند.

 

3- فرمول‌بندی مسأله

در این بخش، بیان ریاضی مسأله پیشنهادی بیان می‌شود که شامل تابع هدف مسأله و قیود آن است. تابع هدف، هزینه‌های مربوط به تولید توان ژنراتورها و فعال‌سازی DR را در بر دارد. قیود مسأله نیز تمامی قیود شبکه، بهره‌برداری، DR وامنیت را شامل می‌شود.

 

3-1- تابع هدف

در مسأله SCUDRC پیشنهادی، هزینه اجرای DR به عنوان بخشی از کل هزینه بهره‌برداری، در کنار هزینه مربوط به تولید توان واحدها در نظر گرفته خواهد شد. با این روش، بهترین ترکیب ممکن از مقدار توان تولیدی واحدها و مقدار DR مبادله شده در هر ساعت به دست می‏آید. بنابراین، تابع هدف، شامل هزینه سوخت واحدهای تولیدی، هزینه روشن[21] و نیز خاموش شدن[22] واحدها و همچنین، هزینه مربوط به مقدار فعال‌سازی برنامه‌های پاسخ‌گویی بار خواهد بود. هزینه روشن شدن واحدها می‌تواند تحت دو حالت مد نظر قرار گیرد: هزینه با راه‌اندازی گرم و هزینه با راه‌اندازی سرد. هزینه خارج شدن واحدها نیز ثابت در نظر گرفته می‌شود. بنابراین، تابع هدف به کمک رابطه (1) بیان می‌شود:

(1)

 

هزینه سوخت واحدهای تولیدی بیشتر به شکل تابعی درجه دوم در نظر گرفته می‌شود که در رابطه (2) نشان داده شده است [43].

(2)

 

 

،  و  ضرایب ثابت در توابع مربوط به هزینه سوخت واحد j ام هستند.

در این مقاله، منحنی خطی مطابق رابطه (3) به عنوان منحنی عرضه DR در هر کدام از باس‌های حاوی بار پاسخ‌گو در نظر گرفته می‌شود.

(3)

 

 

در رابطه (3)، شاخص DPi,t قیمت اجرای DR از سوی عرضه‌کننده i ام را بر حسب ($/MWh) در هر ساعت نشان می‌دهد و ضریب ، "ضریب تمایل" نام دارد که مشخص‌کننده نوع مصرف‌کننده است. این ضریب میزان تمایل مصرف‌کنندگان برای شرکت در برنامه‌های پاسخ‌گویی بار را نشان می‌دهد.  ضریبی بین 0 تا 1 است. افزایش مقدار ضریب ، نشان‌دهنده افزایش تمایل فراهم‌کننده i ام برای شرکت در برنامه‌های DR است؛ زیرا افزایش این ضریب، مطابق رابطه (3)، موجب کاهش هزینه فعال‌سازی DR می‌شود. همچنین، aiو bi، ضرایبی ثابت هستند [29]. با توجه به رابطه (3)، هزینه DR به شکل رابطه (4) خواهد بود [29].

(4)

 

 

هزینه در مدار آمدن واحدها نیز به شکل (5) تعریف می‌شود:

(5)

 

 

در ادامه، قیودی که مسأله عرضه DR و مشارکت واحدها با آن‌ها مواجه است بیان می‌شود.

3-2- قیود مسأله

در این بخش، قیودی که مسأله پیشنهادی با آن‌ها مواجه است معرفی می‌شوند.

 

شرایط اولیه

شرایط اولیه واحدهای تولیدی شامل تعداد ساعاتی که یک واحد به طور پیوسته روشن و یا خاموش بوده و همچنین، مقدار توان تولیدی واحدها یک ساعت قبل از برنامه‌ریزی[23] تولید است.

 

شرط تعادل توان

مجموع توان تولیدی واحدها و مقدار DR فعال‌شده در هر لحظه باید با مقدار بار شبکه در آن ساعت برابر باشد.

(6)

 

 

محدودیت مقدار توان تولیدی واحدها و مشارکت DRSPها

تمامی واحدها مقدار مشخصی را به عنوان کم‌ترین و بیشترین مقدار ممکن توان تولیدی خود دارند. در حقیقت، واحدهای تولیدی قادر به تولید توان بیشتر از یک حد بالایی و همچنین، کمتر از یک حد پایین نیستند. محدودیت تولید به شکل (7) بیان می‌شود:

(7)

 

 

به طور مشابه، محدودیت در مقدار فراخوانی DR به شکل (8) بیان می‌شود:

(8)

 

 

ذخیره چرخان

مقدار ذخیره چرخان[24] مورد نیاز شبکه باید برای اجتناب از هر گونه قطع بار ناخواسته، حاصل از رخدادهای متفاوت سیستم قدرت و یا افزایش غیر منتظره تقاضای بار، کافی باشد. ذخیره چرخان مقداری از پیش تعیین شده است و بیشتر برابر با ظرفیت بزرگ‌ترین واحد و یا درصدی از بار پیش‌بینی شده است [44]. به طور ریاضی‌وار، ذخیره چرخان شبکه در هر ساعت، برابر با بیش‏ترین ظرفیت ممکن تولید توان واحدهای روشن منهای کل توان تولیدی واحدها در آن لحظه است. رابطه (9) این قید را با در نظر گرفتن حضور برنامه‌های پاسخ‌گویی بار نشان می‌دهد:

(9)

 

 

قیود مربوط به نرخ افزایش و کاهش واحدها

نرخ تغییر تولید توان ژنراتورها یکی دیگر از محدودیت‌هایی است که باید در برنامه‌ریزی تولید شبکه مد نظر قرار گیرد. طبق این قید، واحد قادر به تغییر مقدار توان خروجی با هر سرعت دلخواه نیست. مقدار تغییر در خروجی واحدها با محدودیتی روبروست که به قید مربوط به نرخ افزایش و کاهش واحدها[25] موسوم است و به ترتیب توسط معادلات (10) و (11) بیان می‌شوند:

(10)

 

(11)

 
     

 

قید مربوط به حالات روشن یا خاموش شدن واحدها

هیچ کدام از واحدها نمی‌توانند در آن واحد هم روشن و هم خاموش شوند. در نتیجه [25]:

(12)

 

 

از آن‌جا که شاخص‌های SUI و SDI متغیرهایی دودویی هستند که فقط مقادیر صفر و یا یک را دارند، رابطه (12) بیان‌کننده این موضوع است که این شاخص‌ها، نمی‌توانند به طور همزمان یک باشند؛ یعنی هر واحد نمی‌تواند به طور همزمان هم روشن و هم خاموش شود.

 

ناحیه عملکرد ممنوع

ژنراتورها با محدودیت‌هایی در مقدار تولید مواجه هستند که بر اساس آن قادر به تولید توان در مقادیر خاصی از خروجی نیستند [45]. نواحی عملکرد ممنوع[26] واحدها با رابطه (13) بیان می‌شود و در شکل (3) نشان داده شده است.

(13)

 

 

همان‌طور که در شکل (3) نشان داده شده است، ژنراتور قادر به تولید توان در بازه‌هایی که با رنگ تیره مشخص شده است، نمی‌باشد.

 

 

 

شکل (3): نواحی عملکرد ممنوع واحدها و محدودیت تولید توان ژنراتورها

 


قیود مربوط به حداقل زمان روشن و خاموش بودن واحدها (MUT/MDT[27])

حداقل زمان فعالیت یک واحد، کم‌ترین مدت‌زمانی است که واحد پس از روشن شدن، باید به طور پیوسته روشن باقی بماند. این قید توسط رابطه (14) نشان داده می‌شود.

(14)

 

 

قید حداقل زمان توقف، کم‌ترین مدت زمانی است که واحدی پس از خاموش شدن باید به طور پیوسته در وضعیت خاموش باقی بماند و با رابطه (15) بیان می‌شود.

(15)

 

 

قید حداکثر جریان مجاز خطوط

تراکم خطوط به این معناست که توان عبوری از خط نباید بیشتر از حداکثر مجاز بارپذیری خط شود.

(16)

 

 

که در آن  توان عبوری از خط l در دوره t ام و  حداکثر توان مجاز عبوری از خط l است.

 

قیود ولتاز شین‌ها

ولتاژ شین‌های شبکه باید در محدوده مجاز باشد.

(17)

 

 

پیشامد مرتبه‌ی اول (معیار N- 1)[28]

یکی از متداول‌ترین روش‌ها در بررسی امنیت شبکه، استفاده از معیار N- 1 است. این معیار نشان می‌دهد که اگر بنا به علتی، تجهیزی از شبکه خارج شود، تجهیزات دیگر هم‌چنان قادر به تأمین بار شبکه در حین حفظ قیود شبکه هستند و دچار اضافه بار نیز نمی‌شوند. در واقع، مطابق با این معیار، از دست رفتن هر تجهیزی در شبکه، موجب شکست سیستم نمی‌شود.

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

در مسأله SCUDRC پیشنهادی، تابع هدف مسأله که با رابطه (1) فرمول‌بندی شد، باید در طول دوره برنامه‌ریزی با رعایت قیود مختلف، کمینه شود. بنابراین، منابع پاسخ‌گویی بار برای کمتر کردن هزینه‌های بهره‌برداری، از سوی بهره‌بردار مستقل سیستم در کنار تابع تولید توان ژنراتورها قرار می‌گیرند تا بهترین ترکیب ممکن از مقدار توان تولیدی ژنراتورها و مقدار DR فعال شده توسط DRSP ها در هر ساعت تعیین شود.

 

4- مطالعات عددی

در این بخش، سیستم RTS (24 شینه IEEE) برای مطالعات شبیه‌سازی در طی افق برنامه‌ریزی 24 ساعته بررسی می‏شود [46]. در این مقاله، از نرم‌افزار GAMS برای انجام شبیه‌سازی‌های لازم استفاده شده است. شکل (4) تصویر کلی شبکه RTS را نشان می‌دهد که شامل24 باس، 34 خط و 33 واحد تولیدی است.

 

 

 

شکل (4): تصویر کلی شبکه RTS

 

 

 

 

جدول (1): مشخصات فراهم‌کنندگان منابع DR

بیشینه ظرفیت DR (MW)

فراهم‌کننده DR مربوطه

شین‌های حاوی بارهای پاسخ‌گو

15

DRSP1

3

5

DRSP2

5

10

DRSP3

8

30

DRSP4

15

5

DRSP5

16

20

DRSP6

18

40

DRSP7

19

10

DRSP8

20

 

جدول (2): مشخصات سناریوها

هزینه بهره‌برداری ($)

ظرفیت منابع DR

قیود امنیت

نوع مسأله

شماره حالت‌ها

47/313790

-

-

UC

1

88/296420

ü

-

UDRC[29]

2

38/327172

-

ü

SCUC

3

24/313784

ü

ü

SCUDRC

4

 

 

در این قسمت، 8 فراهم‌کننده DR در نظر گرفته می‌شود. در جدول (1)، باس‌های حاوی منابع پاسخ‌گویی بار و مقدار بیشینه ظرفیت DR نشان داده شده است. مقدار ضرایب ai و bi برای تمامی فراهم‌کنندگان DR به ترتیب برابر با 1 و 7 فرض شده است. شایان ذکر است که در ادامه، نتایج خروجی با در نظر گرفتن مقادیر مختلف ضرایب عرضه DR نیز بررسی و تحلیل می‌شوند. زمان اجرای برنامه برای شبکه مورد نظر، با استفاده از رایانه شخصی با RAM 4 GB و CPU 2GHz 32 دقیقه بوده است. شایان ذکر است، زمان اجرا بر روی شبکه‌های بزرگ، نیاز به زمان به مراتب بیش‌تر و استفاده از سیستم‌های قوی‌تر خواهد داشت.

در جدول (2)، مقدار هزینه بهره‌برداری از سیستم تحت چهار حالت مقایسه شده‌اند. همان‌طور که در جدول (2) مشاهده می‌شود، در نظر داشتن ظرفیت منابع پاسخ‌گویی بار موجب کاهش 59/17369 و 14/13388 دلاری هزینه به ترتیب در سناریوی 2 نسبت به 1 و در سناریوی 4 نسبت به 3 می‌شود.

جدول (3)، نتایج برنامه‌ریزی عرضه پاسخ‌گویی بار حاصل از حل مسأله SCUDRC را نشان می‌دهد. همان‌طور که بیان شد، منابع DR به شکل واحدهای تولیدی مجازی با شماره‌های 34 تا 41 نشان داده شده‌اند. نتایج مربوط به نحوه فراخوانی منابع پاسخ‌گویی بار، در شکل (5) نیز به شکل گرافیکی نشان داده شده‌اند. در شکل (5)، میزان DR فعال شده و زمان فراخوانی این منابع، در طی افق برنامه‌ریزی نشان داده شده‌است. در شکل (5)، مقدار مورد نیاز مشارکت هر DRSP بر اساس مقدار بیشینه توان قابل عرضه آن نرمالیزه شده است. مقدار مبنای عرضه DRSP ها، در جدول (4) نشان داده شده‌است. همان‌طور که مشاهده می‌شود، هر کدام از فروشندگان DR با توجه به ضریب

در جدول (3) مشاهده می‌شود، با وجود این‌که منحنی عرضه تمامی DRSP ها یکسان فرض شد، عرضه‌کنندگان 2، 5 و 8، در ساعات بیش‌تری جهت شرکت در DR فراخوانی شده‌اند. این مسأله به علت مکان قرارگیری این منابع است. در واقع، با توجه به ویژگی‏های مکان هر کدام از DRSP ها، منابعی که بیش‌تر به برآورده شدن قیود سیستم کمک می‌کنند، بیش‌تر فراخوانی می‌شوند. به‌علاوه، در ساعات پرباری شبکه لازم شده که از تمامی عرضه‌کنندگان DR استفاده شود. حضور و مشارکت منابع پاسخ‌گویی بار می‌تواند موجب خاموشی واحدهای کم‌بازده انرژی و بهبود کارایی سیستم شود. در واقع، با در نظر گرفتن منابع پاسخ‌گویی بار در کنار منابع سمت تولید، واحدهای کم‌بازده خاموش می‌شوند و مقدار بار لازم توسط عرضه DR فراهم می‌شود. این در حالی است که در صورت عدم حضور منابع پاسخ‌گویی بار، این واحدها باید در ساعاتی در دوره زمانی بهره‌برداری روشن شده و با توجه به قید حداقل زمان در مدار بودن، تا ساعاتی نیز در مدار بمانند. بالا بودن هزینه تولید این واحدها، در کل موجب افزایش هزینه بهره‌برداری از شبکه می‌شود. با طرح مسأله SCUDRC و در نتیجه با در نظر گرفتن پاسخ‌گویی بار به عنوان منبعی برای تأمین بار، این قبیل شرایط نامطلوب در بهره‌برداری از سیستم حذف شده و کاهش در خور توجهی در هزینه‌ها ایجاد می‌شود.

باید دقت شود که کاهش هزینه مطابق جدول (2)، به ازای مقدار خاصی از ظرفیت DR حاصل شده است. در صورت وجود ظرفیت‌های متفاوت منابع پاسخ‌گویی بار، منافع حاصل از آن‌ها نیز تغییر می‌کند. شکل (6)، تغییر در مقدار هزینه بهره‌برداری از سیستم را بر حسب مقادیر مختلف ظرفیت DR در شبکه نشان می‌دهد. در شکل (6) مقادیر متناظر با ظرفیت DR در جدول (2) به عنوان مقادیر پایه (یک پریونیت[30]) فرض شده‌اند. در شکل (6)، تغییر در هزینه بهره‌برداری بر حسب ظرفیت‌های مختلف DR در شبکه، با در نظر گرفتن مقادیر صفر، 2 و 4 به عنوان ضرایب a فراهم‌کنندگان DR نشان داده شده است. در جدول (5)، درصد تغییرات هزینه بهره‌برداری نسبت به حالت مرسوم برنامه‌ریزی مشارکت واحدها (یعنی بدون در نظر داشتن ظرفیت منابع پاسخ‌گویی بار)، نشان داده شده است.

تغییر در مقدار هزینه بهره‌برداری از سیستم بر حسب مقادیر مختلف ضرایب عرضه DR نیز در شکل (7) رسم شده است. افزایش مقدار ضرایب ai و bi در منحنی‌های عرضه DR، موجب افزایش هزینه بهره‌برداری از سیستم خواهد شد. این مسأله به علت کاهش میزان تمایل مشتریان برای شرکت در برنامه‌های پاسخ‌گویی بار در پی افزایش مقادیر این ضرایب است که موجب افزایش هزینه اجرای DR می‌شود. در واقع، افزایش ضرایب عرضه DR نشان دهنده تمایل کمتر DRSP ها برای شرکت در برنامه‌های پاسخ‌گویی بار است. اگر ضرایب عرضه DR افزایش بیشتری داشته باشد، ممکن است DRSP ها نقشی در تأمین بار نداشته باشند. در این شرایط، هزینه بهره‌برداری با هزینه بهره‌برداری در صورت عدم حضور منابع پاسخ‌گویی بار برابر می‌شود.

مقدار سود عرضه‌کنندگان DR بر حسب مقادیر مختلف ضرایت ai و bi نیز در شکل (8) رسم شده است. همان‌طور که مشاهده می‌شود، با افزایش مقادیر این ضرایب، ابتدا مقدار سود فراهم‌کنندگان پاسخ‌گویی بار افزایش و سپس، دوباره کاهش می‌یابد. علت این مسأله این است که، در ابتدا که ضرایب  مقدار کوچکی دارند، قیمت فروش DR کم است و در نتیجه DRSP ها سود کمی را از شرکت در برنامه‌های DR می‌برند. در ادامه با افزایش ضریب ، قیمت اجرای DR افزایش یافته و سود بیشتری نصیب فراهم‌کنندگان DR می‌شود؛ ولی در صورتی که مقدار  افزایش بیشتری داشته باشد، اجرای DR به علت هزینه کمابیش زیاد آن، از دید ISO مقرون به صرفه نبوده و مقدار کمتری از برنامه‌های DR فراخوانی می‌شوند. در واقع در این حالت، کاهش مقدار DR عرضه شده، نسبت به افزایش قیمت آن، عامل تأثیرگذارتری است که با توجه به کم شدن مقدار عرضه، سود DRSP ها نیز کاهش یافته است. مطابق شکل (8)، با افزایش زیاد مقادیر ضرایب عرضه DR، سود فراهم‌کنندگان DR به سمت صفر میل می‌کند‌.

تغییرات درآمد، هزینه و سود فراهم‌کنندگان پاسخ‌گویی بار بر حسب ظرفیت نفوذ DR در شبکه در جدول (6) نشان داده شده است.

شایان ذکر است، زمان اجرای برنامه برای شبکه مورد نظر، به کمک رایانه شخصی با RAM 4 GB و CPU 2GHz در حدود 2 ساعت بوده است. زمان اجرا بر روی شبکه‌های بزرگ، نیازمند زمان به مراتب بیش‌تر و استفاده از سیستم‌های قوی‌تر خواهد بود.

 

جدول (3): مقدار DR فعال شده در طول ساعات برنامه‌ریزی

شماره واحدها

فراهم‌کننده‌ها

ساعات (1-24)

34

DRSP 1

 

 

 

 

 

 

 

 

 

ü

ü

 

 

 

 

 

ü

ü

ü

ü

 

 

 

 

35

DRSP 2

 

 

 

 

 

 

 

 

ü

ü

ü

ü

ü

ü

 

 

ü

ü

ü

ü

 

 

 

 

36

DRSP 3

 

 

 

 

 

 

 

 

 

ü

ü

 

 

ü

 

 

ü

ü

ü

ü

 

 

 

 

37

DRSP 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ü

ü

ü

ü

 

 

 

 

38

DRSP 5

 

 

 

 

 

 

 

 

ü

ü

ü

ü

ü

ü

 

 

ü

ü

ü

ü

 

 

 

 

39

DRSP 6

 

 

 

 

 

 

 

 

 

ü

ü

 

 

 

 

 

ü

ü

ü

ü

 

 

 

 

40

DRSP 7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ü

ü

ü

 

 

 

 

 

41

DRSP 8

 

 

 

 

 

 

 

 

ü

ü

ü

ü

ü

 

 

 

ü

ü

ü

ü

 

 

 

 

 

 

شکل (5): مقدار عرضه DR در طول روز

 

 

جدول (4): مقادیر مبنا در نمودار عرضه DR

واحدها

مقادیر مبنا (MW)

واحدها

مقادیر مبنا (MW)

واحد- 34

5/7

واحد- 38

5

واحد- 35

5

واحد- 39

6/6

واحد- 36

5/7

واحد- 40

10

واحد- 37

5/7

واحد- 41

5/7

 

 

 

شکل (6): هزینه بهره‌برداری بر حسب مقادیر مختلف ظرفیت DR

 

جدول (5): تغییرات هزینه بهره‌برداری به ازای مقادیر مختلف ظرفیت DR

درصد کاهش هزینه در نتیجه حضور منابع پاسخ‌گویی بار

درصد ظرفیت DR به پیک بار

بیشینه ظرفیت DR (p.u)

4=a

2=a

صفر=a

2/3

8/3

9/4

7/4

1

8/3

4/4

5/5

5/9

2

4/4

0/5

6/5

2/14

3

 

 

شکل (7): تغییر در مقدار هزینه بهره‌برداری بر حسب ضرایب عرضه DR

 

 

شکل (8): سود عرضه‌کنندگان DR بر حسب مقادیر مختلف ضرایب عرضه DR

 

جدول (6): تغییرات درآمد، هزینه و سود فراهم‌کنندگان پاسخ‌گویی بار بر حسب ضرایب a و b

سود فراهم‌کنندگان پاسخ‌گویی بار

هزینه فراهم‌کنندگان پاسخ‌گویی بار

درآمد فراهم‌کنندگان پاسخ‌گویی بار

مقدار ضریب b

4=a

2=a

4=a

2=a

4=a

2=a

750

5/1912

750

5/1912

1500

3825

صفر=b

349

5/1087

4/548

5/1807

0/897

2895

2=b

150

375

390

1095

540

1470

4=b

 


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

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

 

6- واژه‌نامه

شاخص‌ها

 

نشانگر باس‌های شبکه

b

نشانگر تعداد فراهم‌کنندگان DR

i

نشانگر تعداد ژنراتورها

j

نشانگر تعداد نواحی عملکرد ممنوع ژنراتورها

n

نشانگر تعداد ساعات در طول دوره برنامه‌ریزی

t

متغیرهای حقیقی

 

مقدار DR فروخته شده توسط فراهم‌کننده i ام در زمان t [MW]

 

توان خروجی واحد j ام در ساعت t [MW]

 

مدت زمانی که واحد j ام به طور مداوم خاموش بوده است [hour]

 

مدت زمانی که واحد j ام به طور مداوم روشن بوده است [hour]

 

ولتاژ باس b در لحظه t[p.u]

 

متغیرهای دودویی

 

نشانگر مربوط به حالت خاموش شدن واحد

 

نشانگر مربوط به حالت روشن شدن واحد

 

نشانگر مربوط به وضعیت فراخوانی عرضه‌کننده i ام DR که در صورت فراخوانی برابر یک و در غیر این صورت برابر صفر است

 

نشانگر وضعیت روشن/خاموش بودن واحد j ام در ساعت t که در صورت روشن بودن برابر یک و در غیر این صورت برابر صفر است

 

ثابت‌ها

 

ضرایب ثابت تابع عرضه DR مربوط به فروشنده i ام ( بر حسب [$/MW2h] و  بر حسب [$/MWh])

 

هزینه روشن شدن ژنراتور به شکل سرد [$/each cold start- up]

CSCj

مدت زمان لازم برای در مدار آمدن ژنراتور j ام در صورت راه‌اندازی به شکل سرد [hour]

 

مقدار تقاضای بار در ساعت t [MW]

Dt

مقدار بیشینه ظرفیت DR فراهم‌کننده i ام در ساعت t [MW]

 

هزینه روشن شدن ژنراتور به شکل گرم
[$/each hot start- up]

HSCj

کم‌ترین مدت زمان برای خاموش ماندن ژنراتور j ام [hour]

 

کم‌ترین مدت زمان برای روشن ماندن ژنراتور j ام [hour]

 

تعداد ژنراتورها

N

تعداد فراهم کنندگان DR

NDRSP

تعداد باس‌های شبکه

NT

کم‌ترین مقدار مجاز تولید توان واحد j ام [MW]

 

بیش‏ترین مقدار مجاز تولید توان واحد j ام [MW]

 

کم‌ترین مقدار مجاز تولید توان واحد j ام در زمان t [MW]

 

بیش‏ترین مقدار مجاز تولید توان واحد j ام در زمان t [MW]

 

حد پایین بازه n ام از ناحیه عملکرد ممنوع واحد j ام [MW]

 

حد بالای بازه n ام از ناحیه عملکرد ممنوع واحد j ام [MW]

 

تعداد نواحی عملکرد ممنوع واحد j ام

 

نرخ کاهش تولید توان واحد j ام [MW/hour]

RDRj

نرخ افزایش تولید توان واحد j ام [MW/hour]

RURj

هزینه خاموش شدن واحد j ام در ساعت t
[$/each switching]

SDCi,t

هزینه روشن شدن واحد j ام در ساعت t
[$/each switching]

SUCi,t

مقدار ذخیره چرخان مورد نیاز در ساعت t [MW]

 

افق زمانی برنامه‌ریزی [hour]

 

بیشینه ولتاژ مجاز شین b ام [p.u]

 

کمینه ولتاژ مجاز شین b ام [p.u]

 

توابع

 

تابع هزینه تولید توان ژنراتور j ام

 

تابع هزینه عرضه DR فراهم‌کننده i ام

 


تاریخ ارسال مقاله : 09/11/1393

تاریخ پذیرش مقاله : 13/04/1394

نام نویسنده مسئول : حبیب‌اله اعلمی

نشانی نویسنده مسئول : ایران – تهران – دانشگاه جامع امام حسین (ع) – دانشکده مهندسی برق



[1] Time-based programs

[2] Incentive-based programs

[3] Market-based programs

[4] Demand response service provider

[5] Independent system operator

[6] Unit commitment

[7] Generation scheduling

[8] Deterministic

[9] Heuristic

[10] Hybrid methods

[11] Security constraint unit commitment

[12] Price responsive loads

[13] Elasticity

[14] Market clearing programs

[15] Price elasticity of demand

[16] Customers’ benefit function

[17] Day-ahead

[18] DR Exchange

[19] Security constraint DR and unit commitment

[20] Demand response resources

[21] Start-up cost

[22] Shut-down cost

[23] Scheduling

[24] Spinning reserve

[25] Unit ramp-up/down constraint

[26] Prohibited operating zone

[27] Minimum up/down time limit

[28] N-1 contingency

[29] Unit and DR commitment

[30] Per-unit (p.u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1- مراجع

[1]      Parvania M., Fotuhi Firuzabad M., “Demand response scheduling by stochastic SCUC”, IEEE Transactions on Smart Grid, Vol. 1, No. 1, pp. 89- 98, 2010.
[2]      IEA, ‘Strategic plan for the IEA demand- side management program 2008- 2012’, (IEA Press 2008), Available at, http://www.iea.org, accessed March 2012.
[3]      Arasteh HR., Parsa Moghaddam M., Sheikh- El- Eslami MK., Abdollahi A., “Integrating commercial demand response resources with unit commitment,” Elect. Power & Energy Systems, Vol. 51, pp. 153- 161, 2013.
[4]      Bai X., Wei H., “Semi- definite programming- based method for security- constrained unit commitment with operational and optimal power flow constraints”, IET Generation, Transmission & Distribution, Vol. 3, No. 2, pp. 182–197, 2009.
[5]      Daneshi H., Srivastava AK., “Security- constrained unit commitment with wind generation and compressed air energy storage”, IET Generation, Transmission & Distribution, Vol. 6, No. 2, pp. 167–175, 2012.
[6]     Yang L., Jian J., Zhu Y., Dong Z., “Tight Relaxation Method for Unit Commitment Problem Using Reformulation and Lift- and- Project”, IEEE Transactions on Power Systems, Vol. 30, No. 1, pp. 13 - 23, 2015.
[7]      An Y., Zeng B., “Exploring the Modeling Capacity of Two- Stage Robust Optimization: Variants of Robust Unit Commitment Model”, IEEE Transactions on Power Systems, Vol. 30, No. 1, pp. 109 - 122, 2015..
[8]      Shafie- khah M., Parsa Moghaddam M., Sheikh- El- Eslami MK., “Unified solution of a non- convex SCUC problem using combination of modified Branch- and- Bound method with Quadratic Programming”, Energy Conversion and Management, Vol. 52, No. 12, pp. 3425–3432, 2011.
[9]      Inostroza JC., Hinojosa VH., “Short- term scheduling solved with a particle swarm optimizer”, IET Generation, Transmssion & Distribution, Vol. 5, No. 11, pp. 1091–1104, 2011.
[10]      Ouyang Z., Shahidehpour SM., “An intelligent dynamic programming for unit commitment application”, IEEE Transactions on Power Systems, Vol. 6, No. 3, pp.1203–1209, 1991.
[11]      Zhang N., Hu Z., Han X., Zhang J., Zhou Y., “A fuzzy chance- constrained program for unit commitment problem considering demand response, electric vehicle and wind power”, International Journal of Electrical Power & Energy Systems, Vol. 65, pp. 201–209, 2015.
[12]      Roy PK., Sarkar R., “Solution of unit commitment problem using quasi- oppositional teaching learning based algorithm”, International Journal of Electrical Power & Energy Systems, Vol. 60, pp. 96–106, 2014.
[13]      Azarinejadian F., Mirhosseini Moghaddam S M., Marzband M., Parhizi N., “Optimal energy management system in grid connected Microgrid integrated with distributed generation by using the multi-period artificial bee colony”, ISEE, Vol. 5, No. 3, pp.1-14, 2014.
[14]      Baniassadi Z., Nezamabadi H., Maghfoori Farsangi M., “A multiobjctive Solution of Gravitational Search Algorithm for Benchmark Functions and Placement of SVC”, ISEE, Vol. 1, No. 1, pp. 59-78, 2011.
[15]      Zheng JH., Chen JJ., Wu QH., Jing ZX.,  “Reliability constrained unit commitment with combined hydro and thermal generation embedded using self- learning group search optimizer”, Energy, In Press, 2015.
[16]      Gjorgiev B., Kančev D., Čepin M., Volkanovski A., “Multi- objective unit commitment with introduction of a methodology for probabilistic assessment of generating capacities availability”, Engineering Applications of Artificial Intelligence, Vol. 37, pp. 236–249, 2015.
[17]      Samiee M, Amjady N, Sharifzadeh H, “Security constrained unit commitment of power systems by a new combinatorial solution strategy composed of enhanced harmony search algorithm and numerical optimization”, International Journal of Electrical Power & Energy Systems, Vol. 44, No. 1, pp. 471–481, 2013.
[18]      Pakkhesal A. R., Aalami H. A., Ebrahimzadeh E., “Control and Management of Electrical Propulsion and Optimal Consumption despite Multiple Electrical and Hybrid Sources”, ISEE, Vol. 5, No. 4, pp. 1-14, 2015.
[19]      Abdoli M., Hatami A., “Dynamic economic emission dispatch incorporating wind farms using modified co-evolutionary particle swarm optimization meta-heuristic algorithm”, ISEE, Vol. 5, No. 4, pp. 31-44, 2015.
[20]      Giftson Samuel G., Christober Asir Rajan C., “Hybrid: Particle Swarm Optimization–Genetic Algorithm and Particle Swarm Optimization–Shuffled Frog Leaping Algorithm for long- term generator maintenance scheduling”, International Journal of Electrical Power & Energy Systems, Vol. 65, pp. 432–442, 2015.
[21]      Valero S., Ortiz M., Senabre C., Alvarez C., Franco FJG., Gabaldon A., “Methods for customer and demand response policies selection in new electricity markets”, IET Generation, Transmission & Distribution, Vol. 1, No. 1, pp. 104- 110, 2007.
[22]      Aalami H.A., Parsa Moghaddam M., Yousefi GR., “Evaluation of nonlinear models for time- based rates demand response programs”, International Journal of Electrical Power & Energy Systems, Vol. 65, pp. 282–290, 2015.
[23]      CL S., Kirschen D., “Quantifying the effect of demand response on electricity markets”, IEEE Transactions on Power Systems, Vol. 24, No. 3, pp. 1199- 207, 2009.
[24]      Khodaei A., Shahidehpour M., Bahramirad S., “SCUC with hourly demand response considering intertemporal load characteristics”, IEEE Transactions on Smart Grid, Vol. 2, No. 3, pp. 564–571, 2011.
[25]      Abdollahi A., Parsa Moghaddam M., Rashidinejad M., Sheikh - El- Eslami MK., “Investigation of economic & environmental - driven demand response measures incorporating UC”, IEEE Transactions on Smart Grid, Vol. 3, No. 1, pp. 12- 25, 2012.
[26]      Aalami H.A., Parsa Moghaddam M., Yousefi GR., “Demand response modeling considering interruptible/curtailable loads and capacity market programs”, Applied Energy, Vol. 87, No. 1, pp. 243- 250, 2010.
[27]      Aalami H.A., Parsa Moghaddam M., Yousefi G.R., “Modeling and prioritizing demand response programs in power markets”, Electric Power Systems Research, Vol. 80, No. 4, pp. 426- 435, 2010.
[28]      Parsa Moghaddam M, Abdollahi A, Rashidinejad M, “Flexible demand response programs modeling in competitive electricity markets”, Applied Energy, Vol. 88, No. 9, pp. 3257- 3269, 2011.
[29]      Nguyen DT, Negnevitsky M, de Groot M, “Pool- based demand response exchange—Concept and modeling”, IEEE Transactions on Power Systems, Vol. 26, No. 3, pp. 1677- 1685, 2011.
[30]      Yousefi S, Yousefi G, “Retail Pricing and Day- Ahead Demand Response in Smart Distribution Networks”, Intelligence Systems in Electrical Engineering, Vol. 4, No. 4, pp. 23- 32, 2014.
[31]      Zhao C, Wang J, Watson JP, Guan Y, “Multi- Stage Robust Unit Commitment Considering Wind and Demand Response Uncertainties”, IEEE Transactions on Power Systems, Vol. 28, No. 3, pp. 2708- 2717, 2013.
[32]      kia M., setayesh nazar M., sepasian M S., “Simultaneous Implementation of TOU and UC Programs Considering Wind Power with MILP Method”, ISEE, Vol. 5, No. 3, pp. 69-80, 2014.
[33]      Kwag HG, Kim JO, “Optimal combined scheduling of generation and demand response with demand resource constraints”, Applied Energy, Vol. 96, pp. 161–170, 2012.
[34]      Kamyab F, Amini M, Sheykhha S, Hasanpour M, Jalali MM, “Demand Response Program in Smart Grid Using Supply Function Bidding Mechanism”, IEEE Transactions on Smart Grid, IEEE Early Access Articles, 2015, DOI: 10.1109/TSG.2015.2430364.
[35]      Sahebi MM, Duki EA, Kia M, Soroudi A, Ehsan M, “Simultaneous emergency demand response programming and unit commitment programming in comparison with interruptible load contracts”, IET Generation, Transmission & Distribution, Vol. 6, No. 7, pp. 605–611, 2012.
[36]      Wang Q, Wang J, Guan Y, “Stochastic Unit Commitment with Uncertain Demand Response”, IEEE Transactions on Power Systems, Vol. 28, No. 1, pp. 562- 563, 2013.
[37]      Aghaei J, Alizadeh MI, “Critical peak pricing with load control demand response program in unit commitment problem”, IET Generation, Transmission & Distribution, Vol. 7, No. 7, pp. 681–690, 2013.
[38]      Conejo AJ, Morales JM., Baringo L, “Real- time demand response model”, IEEE Transactions on Smart Grid, Vol. 1, No. 3, pp. 236–242, 2010.
[39]      Roscoe AJ, Ault G, “Supporting high penetrations of renewable generation via implementation of real time electricity pricing and demand response”, Renewable Power Generation, Vol. 4, No. 4, pp. 369–382, 2010.
[40]      Magnago FH, Alemany J, Lin J, “Impact of demand response resources on unit commitment and dispatch in a day- ahead electricity market”, International Journal of Electrical Power & Energy Systems, Vol. 68, pp. 142–149, 2015.
[41]      Jonghe CD, Hobbs BF, Belmans R, “Value of Price Responsive Load for Wind Integration in Unit Commitment”, IEEE Transactions on power systems, Vol. 29, No. 2, pp. 675–685, 2014.
[42]      Rahmani- andebili M, “Investigating effects of responsive loads models on unit commitment collaborated with demand- side resources”, IET Generation, Transmission & Distribution, Vol. 7, No. 4, pp. 420–430, 2013.
[43]      Yuan X, Ji B, Yuan Y, Ikram RM, Zhang X, Huang Y, “An efficient chaos embedded hybrid approach for hydro- thermal unit commitment problem”, Energy Conversion and Management, Vol. 91, pp. 225–237, 2015.
[44]      Shayesteh E, Yousefi A, Parsa Moghaddam M, “A probabilistic risk- based approach for spinning reserve provision using day- ahead demand response program”, Energy, Vol. 35, No. 5, pp. 1908–1915, 2010.
[45]      Sabera AY, Chakrabortyb S, Abdur Razzakb SM, Senjyub T, “Optimization of economic load dispatch of higher order general cost polynomials and its sensitivity using modified particle swarm optimization”, Electric Power Systems Research, Vol. 79, No. 1, pp. 98–106, 2009.
[46]      Saebi J, Javidi MH, Oloomi Buygi M, “Toward mitigating wind- uncertainty costs in power system operation: A demand response exchange market framework,” Electric Power Systems Research, Vol. 119, pp. 157–167, 2015.