Document Type : Research Article
Authors
1 Shahid Beheshti University
2 Imam Hussein University
Abstract
Keywords
برنامههای پاسخگویی بار اقدامات کوتاهمدتی هستند که توسط مصرفکنندهها برای تعدیل مصرف انرژی انجام میشوند [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)، نشاندهنده ارتباط ایجاد شده بین منابع پاسخگویی بار[20] (DRRs) و منابع سمت تولید است، به طوری که بهرهبردار مستقل سیستم، این منابع را از سمت تولید و تقاضا جمعآوری کرده و با هدف کمینهسازی هزینه بهرهبرداری از شبکه در حین حفظ امنیت سیستم، مقدار توان تولیدی ژنراتورها و نیز مقدار DR فعالسازی شده در هر لحظه را تعیین میکند.
همانطور که در شکل (2) مشاهده میشود، طی چارچوب پیشنهادی، ابتدا رفتار فراهمکنندگان DR، به عنوان منابعی مجازی و در دسترس مدلسازی شده است. در این راستا، برنامههای مبتنی بر بازار پاسخگویی بار مد نظر قرار گرفتهاند که در آن، فراهمکنندگان DR، منحنی عرضه خود را به بهرهبردار مستقل سیستم ارایه میدهند. از طرف دیگر، بهرهبردار سیستم، منحنی عرضه توان واحدهای تولیدی را نیز دریافت میکند. طی مسأله پیشنهادی مقاله، بهرهبردار سیستم قدرت منحنی عرضه پاسخگویی بار را همانند واحدهایی تولیدی در نظر میگیرد که میتوانند در برآورد بار و قیود شبکه سهیم باشند. به این طریق، ظرفیت پاسخگویی بار را میتوان به عنوان منابعی در نظر گرفت که در کنار سایر واحدهای تولیدی در اختیار بهرهبردار سیستم قرار میگیرند.
در این بخش، بیان ریاضی مسأله پیشنهادی بیان میشود که شامل تابع هدف مسأله و قیود آن است. تابع هدف، هزینههای مربوط به تولید توان ژنراتورها و فعالسازی DR را در بر دارد. قیود مسأله نیز تمامی قیود شبکه، بهرهبرداری، DR وامنیت را شامل میشود.
در مسأله 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 و مشارکت واحدها با آنها مواجه است بیان میشود.
در این بخش، قیودی که مسأله پیشنهادی با آنها مواجه است معرفی میشوند.
شرایط اولیه
شرایط اولیه واحدهای تولیدی شامل تعداد ساعاتی که یک واحد به طور پیوسته روشن و یا خاموش بوده و همچنین، مقدار توان تولیدی واحدها یک ساعت قبل از برنامهریزی[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 ها در هر ساعت تعیین شود.
در این بخش، سیستم 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 |
در این مقاله، به منظور بررسی اثر DR بر سمت تولید در بازارهای برق، مسأله برنامهریزی امنیت- مقید واحدهای تولیدی و بارهای پاسخگو مطرح شد. در مسأله ارایه شده، منابع پاسخگویی بار به عنوان کالایی تجاری مدل شدند که میتوانند به طور همزمان در کنار ژنراتورها برای تأمین بار و ذخیره مورد نیاز شبکه به کار روند. نتایج حاصل از ایده پیشنهادی مقاله، مقدار و زمان عرضه DR را در کنار برنامه تولید سایر واحدهای تولیدی شبکه مشخص میکند. به منظور بررسی نتایج عملی از شبکه استاندارد RTS استفاده شد. مقایسه نتایج حاصل از حل مسأله SCUDRC با نتایج مسأله SCUC مرسوم، تأثیر زیاد حضور این منابع را در کاهش هزینه بهرهبرداری از شبکه نشان داد. در ادامه، نتایج عددی با بررسی مقادیر مختلف ظرفیت DR در شبکه و نیز به ازای مقادیر مختلف ظرفیت پاسخگویی بار نیز بررسی و تحلیل شد. همانطور که ملاحظه شد، مدلسازی DR به عنوان منابع مجازی و استفاده از روش پیشنهادی این مقاله، موجب کاهش قابل توجه هزینه بهرهبرداری از شبکه شده است. همچنین، تغییر مقدار هزینه بهرهبرداری بر حسب میزان نفوذ پاسخگویی بار در شبکه نیز بررسی و مشاهده شد که با افزایش مقدار پتانسیل DR، هزینه بهرهبرداری از شبکه کاهش در خور توجهی خواهد داشت. بنابراین، بهرهبردار مستقل سیستم با استفاده از طرح پیشنهادی این مقاله قادر به کاهش قابل توجه هزینه بهرهبرداری از سیستم خواهد بود.
شاخصها |
|
نشانگر باسهای شبکه |
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] |
|
هزینه روشن شدن ژنراتور به شکل گرم |
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 |
SDCi,t |
هزینه روشن شدن واحد j ام در ساعت t |
SUCi,t |
مقدار ذخیره چرخان مورد نیاز در ساعت t [MW] |
|
افق زمانی برنامهریزی [hour] |
|
بیشینه ولتاژ مجاز شین b ام [p.u] |
|
کمینه ولتاژ مجاز شین b ام [p.u] |
|
توابع |
|
تابع هزینه تولید توان ژنراتور j ام |
|
تابع هزینه عرضه DR فراهمکننده i ام |
[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